scholarly journals A partitioned numerical scheme for fluid-structure interaction with slip

Author(s):  
Martina Bukac ◽  
Suncica Canic

We present a loosely coupled, partitioned scheme for fluid-structure interaction problems with the Navier slip boundary condition. The fluid flow is modeled by the Navier-Stokes equations for an incompressible, viscous fluid, interacting with a thin elastic structure modeled by the membrane or Koiter shell type equations. The fluid and structure are coupled via two sets of coupling conditions: a dynamic coupling condition describing balance of forces, and a kinematic coupling condition describing fluid slipping tangentially to the moving fluid-structure interface, with no penetration in the normal direction. We propose a novel, efficient partitioned scheme where the fluid sub-problem is solved separately from the structure sub-problem, and there is no need for sub-iterations to achieve stability, convergence, and its optimal, first-order accuracy. We derive energy estimates, which prove that the proposed scheme is unconditionally stable, and present convergence analysis which shows that the method is first-order accurate in time and optimally convergent in space. The theoretical rates of convergence in time are confirmed numerically on an example with an explicit solution using the method of manufactured solutions. The effects of the slip rate and fluid viscosity on the FSI solution are numerically investigated in two additional examples.

2016 ◽  
Vol 08 (08) ◽  
pp. 1650095 ◽  
Author(s):  
H. Devaraj ◽  
Kean C. Aw ◽  
E. Haemmerle ◽  
R. Sharma

3D printed hair-like micro-structures have been previously demonstrated in a novel micro-fluidic flow sensor aimed at sensing air flows down to rates of a few milliliters per second. However, there is a lack of in-depth understanding of the structural response of these ‘micro-hairs' under a fluid flow field. This paper demonstrates the use of lattice Boltzmann methods (LBM) to understand this structural response towards a better optimization of the micro-hair flow sensors designed to suit the end applications' needs. The LBM approach was chosen as an efficient alternative to simulate Navier–Stokes equations for modeling fluid flow around complex geometries primarily for improved accuracy and simplicity with lesser computational costs. As the spatial dimensions of the sensor's flow channel are much larger in comparison to the actual micro-hairs (the sensing element), a multidimensional approach of combining two-dimensional (D2Q9) and three-dimensional (D3Q19) lattice configurations were implemented for improved computational speeds and efficiency. The drag force on the micro-hairs was estimated using the momentum-exchange method in the D3Q19 configuration and this drag force is transferred to the structural analysis model which determines the micro-hair deformation using Euler–Bernoulli beam theory. The entirety of the LBM Fluid–Structure Interaction (FSI) model was implemented within MATLAB and the obtained results are compared against the numerical model implemented on a commercially available software package.


Author(s):  
M. Benaouicha ◽  
S. Guillou ◽  
A. Santa Cruz ◽  
H. Trigui

The study deals with a 3D Fluid-Structure Interaction (FSI) numerical model of a rectangular cantilevered flexible hydrofoil subjected to a turbulent fluid flow regime. The structural response and dynamic deformations are studied by analyzing the oscillations frequencies and amplitudes, under a hydrodynamics loads. The obtained numerical results are confronted with experimental ones, for validation. The numerical model is performed in the same geometric, physical and material conditions as the experimental set-up carried out in a hydrodynamic tunnel. A polyacetal (POM) flexible hydrofoil NACA0015 with an angle of attack of 8° is considered to be immersed in a fluid flow at a Reynold number of 3 × 105. The structure is initially at rest and then moved by the action of the fluid flow. The numerical model is based on a strong coupling procedure for solving the Fluid-Structure Interaction problem. The Arbitrary Lagrangian-Eulerian (ALE) formulation of the Navier-Stokes equations is used and an anisotropic diffusion equation is solved to compute the fluid mesh velocity and position at each time step. The finite volume method is used for the numerical resolution of the fluid dynamics equations. The structure deformations are described by the linear elasticity equation which is solved by the finite elements method. The Fluid-Structure coupled problem is solved by using the partitioned FSI implicit algorithm. A good agreement between numerical and experimental results for the hydrodynamics coefficients and hydrofoil deformations, maximum deflection and frequencies is obtained. The added mass and damping are analyzed and then the FSI effect on the dynamic deformations of the structure is highlighted.


Author(s):  
Fande Kong ◽  
Xiao-Chuan Cai

Fluid-structure interaction (FSI) problems are computationally very challenging. In this paper we consider the monolithic approach for solving the fully coupled FSI problem. Most existing techniques, such as multigrid methods, do not work well for the coupled system since the system consists of elliptic, parabolic and hyperbolic components all together. Other approaches based on direct solvers do not scale to large numbers of processors. In this paper, we introduce a multilevel unstructured mesh Schwarz preconditioned Newton–Krylov method for the implicitly discretized, fully coupled system of partial differential equations consisting of incompressible Navier–Stokes equations for the fluid flows and the linear elasticity equation for the structure. Several meshes are required to make the solution algorithm scalable. This includes a fine mesh to guarantee the solution accuracy, and a few isogeometric coarse meshes to speed up the convergence. Special attention is paid when constructing and partitioning the preconditioning meshes so that the communication cost is minimized when the number of processor cores is large. We show numerically that the proposed algorithm is highly scalable in terms of the number of iterations and the total compute time on a supercomputer with more than 10,000 processor cores for monolithically coupled three-dimensional FSI problems with hundreds of millions of unknowns.


Author(s):  
Yogesh Ramesh Patel

This paper provides a brief overview of the research in the field of Fluid-structure interaction in Wind Turbines. Fluid-Structure Interaction (FSI) is the interplay of some movable or deformable structure with an internal or surrounding fluid flow. Flow brought about vibrations of two airfoils used in wind turbine blades are investigated by using a strong coupled fluid shape interplay approach. The approach is based totally on a regularly occurring Computational Fluid Dynamics (CFD) code that solves the Navier-Stokes equations defined in Arbitrary Lagrangian-Eulerian (ALE) coordinates by way of a finite extent method. The need for the FSI in the wind Turbine system is studied and comprehensively presented.


Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 94 ◽  
Author(s):  
Cornel Marius Murea

A monolithic semi-implicit method is presented for three-dimensional simulation of fluid–structure interaction problems. The updated Lagrangian framework is used for the structure modeled by linear elasticity equation and, for the fluid governed by the Navier–Stokes equations, we employ the Arbitrary Lagrangian Eulerian method. We use a global mesh for the fluid–structure domain where the fluid–structure interface is an interior boundary. The continuity of velocity at the interface is automatically satisfied by using globally continuous finite element for the velocity in the fluid–structure mesh. The method is fast because we solve only a linear system at each time step. Three-dimensional numerical tests are presented.


2019 ◽  
Vol 2019 ◽  
pp. 1-14
Author(s):  
Ricardo Faria ◽  
Sérgio Oliveira ◽  
Ana L. Silvestre

We consider a dam-water system modeled as a fluid-structure interaction, specifically, a coupled hyperbolic second-order problem, formulated in terms of the displacement of the structure and the fluid pressure. Firstly, we investigate the well posedness of the corresponding variational formulation using Galerkin approximations, energy estimates, and mollification. Then, we apply the finite element method along with the state-space representation of the discrete problem in order to perform a 3D numerical simulation of Cabril arch dam (Zêzere river, Portugal). The numerical model is validated by comparison with available experimental data from a monitoring vibration system installed in Cabril dam.


Author(s):  
Lucia Sargentini ◽  
Benjamin Cariteau ◽  
Morena Angelucci

This paper is related to fluid-structure interaction analysis of sodium cooled fast reactors core (Na-FBR). Sudden liquid evacuation between assemblies could lead to overall core movements (flowering and compaction) causing variations of core reactivity. The comprehension of the structure behavior during the evacuation could improve the knowledge about some SCRAMs for negative reactivity occurred in PHÉNIX reactor and could contribute on the study of the dynamic behavior of a FBR core. An experimental facility (PISE-2c) is designed composed by a Poly-methyl methacrylate hexagonal rods (2D-plan similitude with PHÉNIX assembly) with a very thin gap between assemblies. Another experimental device (PISE-1a) is designed and composed by a single hexagonal rod for testing the dynamic characteristics. Different experiments are envisaged: free vibrations and oscillations during water injection. A phenomenological analysis is reported showing the flow behavior in the gap and the structure response. Also computational simulations are presented in this paper. An efficient numerical method is used to solve Navier-Stokes equations coupled with structure dynamic equation. The numerical method is verified by the comparison of analytic models and experiments.


2008 ◽  
Vol 130 (3) ◽  
Author(s):  
C. G. Giannopapa ◽  
G. Papadakis

In the conventional approach for fluid-structure-interaction problems, the fluid and solid components are treated separately and information is exchanged across their interface. According to the conventional terminology, the current numerical methods can be grouped in two major categories: partitioned methods and monolithic methods. Both methods use separate sets of equations for fluid and solid that have different unknown variables. A unified solution method has been presented in the previous work of Giannopapa and Papadakis (2004, “A New Formulation for Solids Suitable for a Unified Solution Method for Fluid-Structure Interaction Problems,” ASME PVP 2004, San Diego, CA, July, PVP Vol. 491–1, pp. 111–117), which is different from these methods. The new approach treats both fluid and solid as a single continuum; thus, the whole computational domain is treated as one entity discretized on a single grid. Its behavior is described by a single set of equations, which are solved fully implicitly. In this paper, the elastodynamic equations are reformulated so that they contain the same unknowns as the Navier–Stokes equations, namely, velocities and pressure. Two time marching and one spatial discretization scheme, widely used for fluid equations, are applied for the solution of the reformulated equations for solids. Using linear stability analysis, the accuracy and dissipation characteristics of the resulting difference equations are examined. The aforementioned schemes are applied to a transient structural problem (beam bending) and the results compare favorably with available analytic solutions and are consistent with the conclusions of the stability analysis. A parametric investigation using different meshes, time steps, and beam dimensions is also presented. For all cases examined, the numerical solution was stable and robust and therefore is suitable for the next stage of application to full fluid-structure-interaction problems.


2015 ◽  
Vol 07 (02) ◽  
pp. 1550021 ◽  
Author(s):  
Tao He

This paper presents a partitioned strong coupling algorithm for fluid–structure interaction in the arbitrary Lagrangian–Eulerian finite element framework. The incompressible Navier–Stokes equations are solved by the semi-implicit characteristic-based split (CBS) scheme while the structural equations are temporally advanced by the Bathe method. The celled-based smoothed finite element method is adopted for the solution of a geometrically nonlinear solid. To update the dynamic mesh, the moving submesh approach is performed in conjunction with the ortho-semi-torsional spring analogy method. A mass source term is implanted into the pressure Poisson equation to respect the geometric conservation law for the fractional-step-type CBS fluid solver. The iterative solution is achieved by fixed-point method with Aitken's Δ2 accelerator. The proposed methodology is validated against flow-induced oscillations of a bluff body and a flexible body. The overall numerical results agree well with the available data. Some important flow phenomena have been disclosed successfully.


1998 ◽  
Vol 08 (04) ◽  
pp. 543-572 ◽  
Author(s):  
F. FLORI ◽  
P. ORENGA

We present in this paper an existence result of weak solutions as well as some regularity results for a fluid-structure interaction problem when a fluid velocity-structure displacement formulation is used. This formulation induces a difficulty connected with the coupling condition. Indeed, unlike the pressure-displacement formulation, we do not have a Neumann condition and consequently it is more difficult to obtain weak solutions.


Sign in / Sign up

Export Citation Format

Share Document