×

A stable partitioned FSI algorithm for incompressible flow and deforming beams. (English) Zbl 1351.74157

Summary: An added-mass partitioned (AMP) algorithm is described for solving fluid-structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier-Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler-Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. The key AMP interface condition is a generalized Robin (mixed) condition on the fluid pressure. This condition, which is derived at a continuous level, has no adjustable parameters and is applied at the discrete level to couple the partitioned domain solvers. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet-Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.

MSC:

74S20 Finite difference methods applied to problems in solid mechanics
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
74F10 Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.)
74K10 Rods (beams, columns, shafts, arches, rings, etc.)
76D05 Navier-Stokes equations for incompressible viscous fluids

Software:

Ogen; CMPGRD
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Causin, P.; Gerbeau, J. F.; Nobile, F., Added-mass effect in the design of partitioned algorithms for fluid-structure problems, Comput. Methods Appl. Mech. Eng., 194, 4506-4527 (2005) · Zbl 1101.74027
[2] van Brummelen, E. H., Added mass effects of compressible and incompressible flows in fluid-structure interaction, J. Appl. Mech., 76, 2, 021206-021212 (2009)
[3] Yu, Y.; Baek, H.; Karniadakis, G. E., Generalized fictitious methods for fluid-structure interactions: analysis and simulations, J. Comput. Phys., 245, 317-346 (2013) · Zbl 1349.76577
[4] Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., An analysis of a new stable partitioned algorithm for FSI problems. Part I: incompressible flow and elastic solids, J. Comput. Phys., 269, 108-137 (2014) · Zbl 1349.74373
[5] Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., An analysis of a new stable partitioned algorithm for FSI problems. Part II: incompressible flow and structural shells, J. Comput. Phys., 268, 399-416 (2014) · Zbl 1349.74372
[6] Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., Deforming composite grids for solving fluid structure problems, J. Comput. Phys., 231, 9, 3518-3547 (2012) · Zbl 1402.74030
[7] Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., An added-mass partitioned algorithm for fluid-structure interactions of compressible fluids and nonlinear solids, J. Comput. Phys., 305, 1037-1064 (2016) · Zbl 1349.76428
[8] Fernández, M. A., Coupling schemes for incompressible fluid-structure interaction: implicit, semi-implicit and explicit, SeMa J., 55, 1, 59-108 (2011) · Zbl 1242.76201
[9] Degroote, J., Partitioned simulation of fluid-structure interaction, Arch. Comput. Methods Eng., 20, 3, 185-238 (2013) · Zbl 1354.74066
[10] Bukač, M.; Čanič, S.; Muha, B., A partitioned scheme for fluid-composite structure interaction problems, J. Comput. Phys., 281, 493-517 (2015) · Zbl 1351.76321
[11] Nobile, F.; Pozzoli, M.; Vergara, C., Time accurate partitioned algorithms for the solution of fluid-structure interaction problems in haemodynamics, Comput. Fluids, 86, 470-482 (2013) · Zbl 1290.76166
[12] Liu, J.; Jaiman, R. K.; Gurugubelli, P. S., A stable second-order scheme for fluid-structure interaction with strong added-mass effects, J. Comput. Phys., 270, 687-710 (2014) · Zbl 1349.76236
[13] Guidoboni, G.; Glowinski, R.; Cavallini, N.; Čanić, S., Stable loosely-coupled-type algorithm for fluid-structure interaction in blood flow, J. Comput. Phys., 228, 18, 6916-6937 (2009) · Zbl 1261.76056
[14] Čanić, S.; Muha, B.; Bukač, M., Stability of the kinematically coupled \(β\)-scheme for fluid-structure interaction problems in hemodynamics · Zbl 1332.35281
[15] Bukač, M.; Čanič, S.; Glowinski, R.; Tambača, J.; Quaini, A., Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement, J. Comput. Phys., 235, 515-541 (2013)
[16] Lukáčová-Medvid’ová, M.; Rusnáková, G.; Hundertmark-Zaušková, A., Kinematic splitting algorithm for fluid-structure interaction in hemodynamics, Comput. Methods Appl. Mech. Eng., 265, 83-106 (2013) · Zbl 1286.76176
[17] Fernández, M. A., Incremental displacement-correction schemes for the explicit coupling of a thin structure with an incompressible fluid, C. R. Math., 349, 7-8, 473-477 (2011) · Zbl 1344.76026
[18] Fernández, M. A., Incremental displacement-correction schemes for incompressible fluid-structure interaction: stability and convergence analysis, Numer. Math., 123, 1, 21-65 (2013) · Zbl 1331.76068
[19] Fernández, M. A.; Mullaert, J.; Vidrascu, M., Explicit Robin-Neumann schemes for the coupling of incompressible fluids with thin-walled structures, Comput. Methods Appl. Mech. Eng., 267, 566-593 (2013) · Zbl 1286.76041
[20] Fernández, M. A.; Landajuela, M., A fully decoupled scheme for the interaction of a thin-walled structure with an incompressible fluid, C. R. Math., 351, 3-4, 161-164 (2013) · Zbl 1307.74031
[21] Fernández, M. A.; Landajuela, M., Fully decoupled time-marching schemes for incompressible fluid/thin-walled structure interaction (Jan. 2014), INRIA, Rapport de Recherche RR-8425
[22] Fernández, M. A.; Landajuela, M.; Vidrascu, M., Fully decoupled time-marching schemes for incompressible fluid/thin-walled structure interaction (Dec. 2014), Inria, Research Report RR-8425
[23] Henshaw, W. D., A fourth-order accurate method for the incompressible Navier-Stokes equations on overlapping grids, J. Comput. Phys., 113, 1, 13-25 (1994) · Zbl 0808.76059
[24] Banks, J. W.; Sjögreen, B., A normal mode stability analysis of numerical interface conditions for fluid/structure interaction, Commun. Comput. Phys., 10, 2, 279-304 (2011), publications/BanksSjogreenFSI2011.pdf · Zbl 1364.74096
[25] Henshaw, W. D.; Petersson, N. A., A split-step scheme for the incompressible Navier-Stokes equations, (Hafez, M. M., Numerical Simulation of Incompressible Flows (2003), World Scientific), 108-125 · Zbl 1074.76036
[26] Miller, J. J.H., On the location of zeros of certain classes of polynomials with applications to numerical analysis, IMA J. Appl. Math., 8, 3, 397-406 (1971) · Zbl 0232.65070
[27] Strikwerda, J. C., Finite Difference Schemes and Partial Differential Equations (1989), Wadsworth and Brooks/Cole · Zbl 0681.65064
[28] Henshaw, W. D., The Overture hyperbolic grid generator, user guide, version 1.0 (1999), Lawrence Livermore National Laboratory, Research Report UCRL-MA-134240
[29] Henshaw, W. D., Ogen: an overlapping grid generator for Overture (1998), Lawrence Livermore National Laboratory, Research Report UCRL-MA-132237
[30] Henshaw, W. D.; Schwendeman, D. W., Moving overlapping grids with adaptive mesh refinement for high-speed reactive and non-reactive flow, J. Comput. Phys., 216, 2, 744-779 (2006) · Zbl 1220.76052
[31] Appelö, D.; Banks, J. W.; Henshaw, W. D.; Schwendeman, D. W., Numerical methods for solid mechanics on overlapping grids: linear elasticity, J. Comput. Phys., 231, 18, 6012-6050 (2012) · Zbl 1277.74005
[32] Chesshire, G. S.; Henshaw, W. D., Composite overlapping meshes for the solution of partial differential equations, J. Comput. Phys., 90, 1, 1-64 (1990) · Zbl 0709.65090
[33] Roache, P. J., Verification and Validation in Computational Science and Engineering (1998), Hermosa Publishers: Hermosa Publishers New Mexico
[34] Formaggia, L.; Gerbeau, J.; Nobile, F.; Quarteroni, A., On the coupling of 3D and 1D Navier-Stokes equations for flow problems in compliant vessels, Comput. Methods Appl. Mech. Eng., 191, 6-7, 561-582 (2001) · Zbl 1007.74035
[35] Henshaw, W. D.; Schwendeman, D. W., Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement, J. Comput. Phys., 227, 16, 7469-7502 (2008) · Zbl 1213.76138
[36] Henshaw, W. D., Mappings for Overture, a description of the Mapping class and documentation for many useful Mappings (1998), Lawrence Livermore National Laboratory, Research Report UCRL-MA-132239
[37] Graff, K. F., Wave Motion in Elastic Solids (1991), Dover Publications
[38] Banks, J. W.; Henshaw, W. D.; Sjögreen, B., A stable FSI algorithm for light rigid bodies in compressible flow, J. Comput. Phys., 245, 399-430 (2013) · Zbl 1349.76429
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.