Existence of a Weak Solution to a Nonlinear Fluid–Structure Interaction Problem Modeling the Flow of an Incompressible, Viscous Fluid in a Cylinder with Deformable Walls

We study a nonlinear, unsteady, moving boundary, fluid–structure interaction (FSI) problem arising in modeling blood flow through elastic and viscoelastic arteries. The fluid flow, which is driven by the time-dependent pressure data, is governed by two-dimensional incompressible Navier–Stokes equations, while the elastodynamics of the cylindrical wall is modeled by the one-dimensional cylindrical Koiter shell model. Two cases are considered: the linearly viscoelastic and the linearly elastic Koiter shell. The fluid and structure are fully coupled (two-way coupling) via the kinematic and dynamic lateral boundary conditions describing continuity of velocity (the no-slip condition), and the balance of contact forces at the fluid–structure interface. We prove the existence of weak solutions to the two FSI problems (the viscoelastic and the elastic case) as long as the cylinder radius is greater than zero. The proof is based on a novel semi-discrete, operator splitting numerical scheme, known as the kinematically coupled scheme, introduced in Guidoboni et al. (J Comput Phys 228(18):6916–6937, 2009) to numerically solve the underlying FSI problems. The backbone of the kinematically coupled scheme is the well-known Marchuk–Yanenko scheme, also known as the Lie splitting scheme. We effectively prove convergence of that numerical scheme to a solution of the corresponding FSI problem.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

A new approach for computing the steady state fluid–structure interaction response of periodic problems

A special type of fluid–structure interaction (FSI) problems are problems with periodic boundary conditions like in turbomachinery. The steady state FSI response of these problems is usually calculated with similar techniques as used for transient FSI analyses. This means that, when the fluid and structure problem are not simultaneously solved with a monolithic approach, the problem is partitioned into a fluid and structural part and that each time step coupling iterations are performed to account for strong interactions between the two sub-domains. This paper shows that a time-partitioned FSI computation can be very inefficient to compute the steady state FSI response of periodic problems. A new approach is introduced in which coupling iterations are performed on periodic level instead of per time step. The convergence behaviour can be significantly improved by implementing existing partitioned solution methods as used for time step coupling (TSC) algorithms in the time periodic coupling (TPC) framework. The new algorithm has been evaluated by comparing the convergence behaviour to TSC algorithms. It is shown that the number of fluid–structure evaluations can be considerably reduced when a TPC algorithm is applied instead of a TSC. One of the most appealing advantages of the TPC approach is that the structural problem can be solved in the frequency domain resulting in a very efficient algorithm for computing steady state FSI responses.     

A fast method for solving fluid–structure interaction problems numerically

The paper presents a semi‐implicit algorithm for solving an unsteady fluid–structure interaction problem. The algorithm for solving numerically the fluid–structure interaction problems was obtained by combining the backward Euler scheme with a semi‐implicit treatment of the convection term for the Navier–Stokes equations and an implicit centered scheme for the structure equations. The structure is governed either by the linear elasticity or by the non‐linear St Venant–Kirchhoff elasticity models. At each time step, the position of the interface is predicted in an explicit way. Then, an optimization problem must be solved, such that the continuity of the velocity as well as the continuity of the stress hold at the interface. During the Broyden, Fletcher, Goldforb, Shano (BFGS) iterations for solving the optimization problem, the fluid mesh does not move, which reduces the computational effort. The term ‘semi‐implicit’ used for the fully algorithm means that the interface position is computed explicitly, while the displacement of the structure, velocity and the pressure of the fluid are computed implicitly. Numerical results are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

An embedded boundary framework for compressible turbulent flow and fluid–structure computations on structured and unstructured grids

The finite volume method with exact two‐phase Riemann problems (FIVER) is a two‐faceted computational method for compressible multi‐material (fluid–fluid, fluid–structure, and multi‐fluid–structure) problems characterized by large density jumps, and/or highly nonlinear structural motions and deformations. For compressible multi‐phase flow problems, FIVER is a Godunov‐type discretization scheme characterized by the construction and solution at the material interfaces of local, exact, two‐phase Riemann problems. For compressible fluid–structure interaction (FSI) problems, it is an embedded boundary method for computational fluid dynamics (CFD) capable of handling large structural deformations and topological changes. Originally developed for inviscid multi‐material computations on nonbody‐fitted structured and unstructured grids, FIVER is extended in this paper to laminar and turbulent viscous flow and FSI problems. To this effect, it is equipped with carefully designed extrapolation schemes for populating the ghost fluid values needed for the construction, in the vicinity of the fluid–structure interface, of second‐order spatial approximations of the viscous fluxes and source terms associated with Reynolds averaged Navier–Stokes (RANS)‐based turbulence models and large eddy simulation (LES). Two support algorithms, which pertain to the application of any embedded boundary method for CFD to the robust, accurate, and fast solution of FSI problems, are also presented in this paper. The first one focuses on the fast computation of the time‐dependent distance to the wall because it is required by many RANS‐based turbulence models. The second algorithm addresses the robust and accurate computation of the flow‐induced forces and moments on embedded discrete surfaces, and their finite element representations when these surfaces are flexible. Equipped with these two auxiliary algorithms, the extension of FIVER to viscous flow and FSI problems is first verified with the LES of a turbulent flow past an immobile prolate spheroid, and the computation of a series of unsteady laminar flows past two counter‐rotating cylinders. Then, its potential for the solution of complex, turbulent, and flexible FSI problems is also demonstrated with the simulation, using the Spalart–Allmaras turbulence model, of the vertical tail buffeting of an F/A‐18 aircraft configuration and the comparison of the obtained numerical results with flight test data. Copyright © 2014 John Wiley & Sons, Ltd.     

An immersed smoothed point interpolation method (IS‐PIM) for fluid‐structure interaction problems

An immersed smoothed point interpolation method using 3‐node triangular background cells is proposed to solve 2D fluid‐structure interaction problems for solids with large deformation/displacement placed in incompressible viscous fluid. In the framework of immersed‐type method, the governing equations can be decomposed into 3 parts on the basis of the fictitious fluid assumption. The incompressible Navier‐Stokes equations are solved using the semi‐implicit characteristic‐based split scheme, and solids are simulated using the newly developed edge‐based smoothed point interpolation method. The fictitious fluid domain can be used to calculate the coupling force. The numerical results show that immersed smoothed point interpolation method can avoid remeshing for moving solid based on immersed operation and simulate the contact phenomenon without an additional treatment between the solid and the fluid boundary. The influence from information transfer between solid domain and fluid domain on fluid‐structure interaction problems has been investigated. The numerical results show that the proposed interpolation schemes will generally improve the accuracy for simulating both fluid flows and solid structures.     

Numerical stabilities of loosely coupled methods for robust modeling of lightweight and flexible structures in incompressible and viscous flows

The growing interest to examine the hydroelastic dynamics and stabilities of lightweight and flexible materials requires robust and accurate fluid–structure interaction(FSI)models. Classically, partitioned fluid and structure solvers are easier to implement compared to monolithic methods;however, partitioned FSI models are vulnerable to numerical("virtual added mass") instabilities for cases when the solid to fluid density ratio is low and if the flow is incompressible.As a partitioned method, the loosely hybrid coupled(LHC)method, which was introduced and validated in Young et al.(Acta Mech. Sin. 28:1030–1041, 2012), has been successfully used to efficiently and stably model lightweight and flexible structures. The LHC method achieves its numerical stability by, in addition to the viscous fluid forces, embedding potential flow approximations of the fluid induced forces to transform the partitioned FSI model into a semi-implicit scheme. The objective of this work is to derive and validate the numerical stability boundary of the LHC. The results show that the stability boundary of the LHC is much wider than traditional loosely coupled methods for a variety of numerical integration schemes. The results also show that inclusion of an estimate of the fluid inertial forces is the most critical to ensure the numerical stability when solving for fluid–structure interaction problems involving cases with a solid to fluid-added mass ratio less than one.     

Coupled solution of the species conservation equations using unstructured finite‐volume method

A coupled solver was developed to solve the species conservation equations on an unstructured mesh with implicit spatial as well as species‐to‐species coupling. First, the computational domain was decomposed into sub‐domains comprised of geometrically contiguous cells—a process similar to additive Schwarz decomposition. This was done using the binary spatial partitioning algorithm. Following this step, for each sub‐domain, the discretized equations were developed using the finite‐volume method, and solved using an iterative solver based on Krylov sub‐space iterations, that is, the pre‐conditioned generalized minimum residual solver. Overall (outer) iterations were then performed to treat explicitness at sub‐domain interfaces and nonlinearities in the governing equations. The solver is demonstrated for both two‐dimensional and three‐dimensional geometries for laminar methane–air flame calculations with 6 species and 2 reaction steps, and for catalytic methane–air combustion with 19 species and 24 reaction steps. It was found that the best performance is manifested for sub‐domain size of 2000 cells or more, the exact number depending on the problem at hand. The overall gain in computational efficiency was found to be a factor of 2–5 over the block (coupled) Gauss–Seidel procedure. All calculations were performed on a single processor machine. The largest calculations were performed for about 355 000 cells (4.6 million unknowns) and required 900 MB of peak runtime memory and 19 h of CPU on a single processor. Copyright © 2009 John Wiley & Sons, Ltd.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

A fixed‐grid b‐spline finite element technique for fluid–structure interaction

We present a fixed‐grid finite element technique for fluid–structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b‐spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision‐stabilisation technique is used to ensure inf–sup stability. The beam equations are discretised with b‐splines and the shell equations with subdivision basis functions, both leading to a rotation‐free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet–Robin partitioning scheme, and the fluid equations are solved with a pressure–correction method. Auxiliary techniques employed for improving numerical robustness include the level‐set based implicit representation of the structure interface on the fluid grid, a cut‐cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. Copyright © 2013 John Wiley & Sons, Ltd.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A hybrid immersed boundary and material point method for simulating 3D fluid–structure interaction problems

A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd.     

Time‐related element‐free Taylor–Galerkin method with non‐splitting decoupling process for incompressible steady flow

The time‐related element‐free Taylor–Galerkin method with non‐splitting decoupling process (EFTG‐NSD) is proposed for the simulation of steady flows. The goal of the present paper is twofold. One is to raise the efficiency of the time‐related methods for solving steady flow problems, and the other is to obtain a good stability. The EFTG‐NSD method, which uses the time‐related Navier–Stokes equations to describe steady flows, does not care about the intermediate process and obtains solution of steady flows through time marching. Different from the classical time‐related fractional step methods, the EFTG‐NSD method decouples the Navier–Stokes equations without any operator‐splitting and correction. Because the elimination of correction at each iteration step reduces the computation cost, the EFTG‐NSD method possesses higher computation efficiency. In addition, the EFTG‐NSD method has a good stability due to the use of the Taylor–Galerkin formula in time and space discretization. Furthermore, the method combining element‐free Galerkin method with Taylor–Galerkin method is an important supplement of the element‐free Galerkin method for solving flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Strongly coupled simulation of fluid–structure interaction in a Francis hydroturbine

This work simulates a complex fluid flow in fluid–structure interaction (FSI). The flow under consideration is governed by Navier–Stokes equations for incompressible viscous fluids and modeled with the finite volume method. Large eddy simulation is used to simulate the unsteady turbulent flow. The structure is represented by a finite element formulation. The present work introduces a strongly coupled partitioned approach that is applied to complex flow in fluid machinery. In this approach, the fluid and structure equations are solved separately using different solvers, but are implicitly coupled into one single module based on sensitivity analysis of the important displacement and stress modes. The applied modes and their responses are used to build up a reduced‐order model. The proposed model is used to predict the unsteady flow fields of a 3D complete passage, involving in stay, guide vanes, and runner blades, for a Francis hydro turbine and FSI is considered. The computational results show that a fairly good convergence solution is achieved by using the reduced‐order model that is based on only a few displacement and stress modes, which largely reduces the computational cost, compared with traditional approaches. At the same time, a comparison of the numerical results of the model with available experimental data validates the methodology and assesses its accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

Improving the CBS‐based partitioned semi‐implicit coupling algorithm for fluid‐structure interaction

We detail in this work 2 simple but effective alternatives to improve the characteristic‐based split–based partitioned semi‐implicit coupling algorithm for fluid‐structure interaction. The basic idea lies in introducing the end‐of‐step velocity into the implicit stages of the 2 algorithms integrating different splits. The algorithm built upon the second‐order pressure split is further stabilized via the pressure gradient projection with particular emphasis on the extremely low mass ratio. The smoothed finite element method is exploited for spatial discretization of fluid and solid equations. Even without any accelerators, both the semi‐implicit solvers incorporating fixed‐point iterations engender visible improvements versus the previously published data for several benchmarks.