A stable partitioned FSI algorithm for rigid bodies and incompressible flow in three dimensions

This paper describes a novel partitioned algorithm for fluid–structure interaction (FSI) problems that couples the motion of rigid bodies and incompressible flow. This is the first partitioned algorithm that remains stable and second-order accurate, without sub-time-step iterations, for very light,...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Journal of computational physics Ročník 373; číslo C; s. 455 - 492
Hlavní autori: Banks, J.W., Henshaw, W.D., Schwendeman, D.W., Tang, Qi
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Cambridge Elsevier Inc 15.11.2018
Elsevier Science Ltd
Elsevier
Predmet:
Added-mass
Algorithms
Benchmarks
Computational fluid dynamics
Computational physics
Computer simulation
Damping
Dimensional stability
Fluid flow
Fluid mechanics
Fluid pressure
Fluid-structure interaction
Heart valves
Incompressible flow
Incompressible fluids
Incompressible Navier–Stokes
Interface stability
Mathematical models
Moving overlapping grids
Navier-Stokes equations
Partitioned schemes
Rigid bodies
Rigid structures
Stability analysis
Tensors
Three dimensional analysis
Three dimensional flow
Three dimensional models
Rigid bodies
Moving overlapping grids
Fluid–structure interaction
Incompressible Navier–Stokes
Added-mass
Partitioned schemes
ISSN:0021-9991, 1090-2716
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí:This paper describes a novel partitioned algorithm for fluid–structure interaction (FSI) problems that couples the motion of rigid bodies and incompressible flow. This is the first partitioned algorithm that remains stable and second-order accurate, without sub-time-step iterations, for very light, and even zero-mass, bodies in three dimensions. This new added-mass partitioned (AMP) algorithm extends the previous developments in [1,2] by generalizing the added-damping tensors to account for arbitrary three-dimensional rotations, and by employing a general quadrature for the surface integral over a rigid body to derive the discrete AMP interface condition for the fluid pressure. Stability analyses for two three-dimensional model problems show that the algorithm remains stable for bodies of any mass when applied to the relevant model problems. The resulting AMP algorithm is implemented in parallel using a moving composite grid framework to treat one or more rigid bodies in complex three-dimensional configurations. The new three-dimensional algorithm is verified and validated though several benchmark problems, including the motion of a sphere in a viscous incompressible fluid and the interaction of a bi-leaflet mechanical heart valve and a pulsating fluid. Numerical simulations confirm the predictions of the stability analysis even for complex problems, and show that the AMP algorithm remains stable, without sub-iterations, for light and even zero-mass three-dimensional rigid bodies of general shape. These benchmark problems are further used to examine the parallel performance of the algorithm and to investigate the conditioning of the linear system for the pressure including the newly derived AMP interface conditions. •A partitioned FSI scheme for incompressible flow and rigid bodies is extended to 3D.•Stability for simple geometries is proved using 3D model problems.•Numerical results confirm stability and accuracy for light bodies of general shapes.•A benchmark for a bi-leaflet mechanical heart valve is designed and conducted.•Refinement studies for the heart valve are provided.
Bibliografia:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2018.06.072