A modular numerical method for implicit 0D/3D coupling in cardiovascular finite element simulations. (English) Zbl 1377.76041
Summary: Implementation of boundary conditions in cardiovascular simulations poses numerical challenges due to the complex dynamic behavior of the circulatory system. The use of elaborate closed-loop lumped parameter network (LPN) models of the heart and the circulatory system as boundary conditions for computational fluid dynamics (CFD) simulations can provide valuable global dynamic information, particularly for patient specific simulations. In this paper, the necessary formulation for coupling an arbitrary LPN to a finite element Navier-Stokes solver is presented. A circuit analogy closed-loop LPN is solved numerically, and pressure and flow information is iteratively passed between the 0D and 3D domains at interface boundaries, resulting in a time-implicit scheme. For Neumann boundaries, an implicit method, regardless of the LPN, is presented to achieve the desired stability and convergence properties. Numerical procedures for passing flow and pressure information between the 0D and 3D domains are described, and implicit, semi-implicit, and explicit quasi-Newton formulations are compared. The issue of divergence in the presence of backflow is addressed via a stabilized boundary formulation. The requirements for coupling Dirichlet boundary conditions are also discussed and this approach is compared in detail to that of the Neumann coupled boundaries. Having the option to select between Dirichlet and Neumann coupled boundary conditions increases the flexibility of current framework by allowing a wide range of components to be used at the 3D-0D interface.
MSC:
| 76Z05 | Physiological flows |
| 65M60 | Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs |
| 76D05 | Navier-Stokes equations for incompressible viscous fluids |
| 92C35 | Physiological flow |
Keywords:
lumped parameter network; Neumann; Dirichlet; boundary condition; backflow stabilization; patient-specific blood flow; Navier-Stokes FEM solver; multi-domain method; multi-scale modelingReferences:
| [1] | Taylor, C.; Figueroa, C., Patient-specific modeling of cardiovascular mechanics, Annual Review of Biomedical Engineering, 11, 109-134 (2009) |
| [2] | Taylor, C.; Steinman, D., Image-based modeling of blood flow and vessel wall dynamics: applications, methods and future directions, Annals of Biomedical Engineering, 38, 1188-1203 (2010) |
| [3] | Pennati, G.; Migliavacca, F.; Dubini, G.; Bove, E., Modeling of systemic-to-pulmonary shunts in newborns with a univentricular circulation: state of the art and future directions, Progress in Pediatric Cardiology, 30, 1-2, 23-29 (2010) |
| [4] | Vignon-Clementel, I.; Marsden, A.; Feinstein, J., A primer on computational simulation in congenital heart disease for the clinician, Progress in Pediatric Cardiology, 30, 3-13 (2010) |
| [5] | Westerhof, N.; Bosman, F.; Vries, C. D.; Noordergraaf, A., Analog studies of the human systemic arterial tree, Journal of Biomechanics, 2, 2, 121-143 (1969) |
| [6] | Formaggia, L.; Gerbeau, J.; Nobile, F.; Quarteroni, A., Numerical treatment of defective boundary conditions for the Navier-Stokes equations, SIAM Journal on Numerical Analysis, 40, 376-401 (2002) · Zbl 1020.35070 |
| [7] | Quarteroni, A.; Veneziani, A., Analysis of a geometrical multiscale model based on the coupling of ODEs and PDEs for blood flow simulations, Multiscale Model Simulations, 1, 2, 173-195 (2003) · Zbl 1060.35003 |
| [8] | Vignon-Clementel, I.; Figueroa, C.; Jansen, K.; Taylor, C., Outflow boundary conditions for three-dimensional finite element modeling of blood flow and pressure in arteries, Computer Methods in Applied Mechanics and Engineering, 195, 29-32, 3776-3796 (2006) · Zbl 1175.76098 |
| [9] | Vignon-Clementel, I.; Figueroa, C.; Jansen, K.; Taylor, C., Outflow boundary conditions for three-dimensional simulations of non-periodic blood flow and pressure fields in deformable arteries, Computer Methods in Biomechanics and Biomedical Engineering, 13, 5, 625-640 (2010) |
| [10] | Kim, H.; Vignon-Clementel, I.; Figueroa, C.; Ladisa, J.; Jansen, K.; Feinstein, J.; Taylor, C., On coupling a lumped parameter heart model and a three-dimensional finite element aorta model, Annual Review of Biomedical Engineering, 37, 2153-2169 (2009) |
| [11] | Kim, H.; Vignon-Clementel, I.; Coogan, J.; Figueroa, C.; Jansen, K.; Taylor, C., Patient-specific modeling of blood flow and pressure in human coronary arteries, Annals of Biomedical Engineering, 38, 3195-3209 (2010) |
| [13] | Urquiza, S.; Blanco, P.; Venere, M.; Feijoo, R., Multidimensional modelling for the carotid artery blood flow, Computer Methods in Applied Mechanics and Engineering, 195, 33-36, 4002-4017 (2006) · Zbl 1178.76395 |
| [14] | Blanco, P.; Feijoo, R.; Urquiza, S., A unified variational approach for coupling 3D-1D models and its blood flow applications, Computer Methods in Applied Mechanics and Engineering, 196, 41-44, 4391-4410 (2007) · Zbl 1173.76430 |
| [15] | Quarteroni, A.; Ragni, S.; Veneziani, A., Coupling between lumped and distributed models for blood flow problems, Computing and Visualization in Science, 4, 111-124 (2001) · Zbl 1097.76615 |
| [16] | Leiva, J.; Blanco, P.; Buscaglia, G., Iterative strong coupling of dimensionally-heterogeneous models, International Journal of Numerical Methods in Engineering, 81, 12, 1558-1580 (2010) · Zbl 1183.76838 |
| [17] | Malossi, A.; Blanco, P.; Deparis, S.; Quarteroni, A., Algorithms for the partitioned solution of weakly coupled fluid models for cardiovascular flows, International Journal of Numerical Methods in Biomedical Engineering, 27, 12, 2035-2057 (2011) · Zbl 1452.76287 |
| [18] | Johnson, D.; Naik, P.; Beris, A., Efficient implementation of the proper outlet flow conditions in blood flow simulations through asymmetric arterial bifurcations, International Journal for Numerical Methods in Fluids, 66, 11, 1383-1408 (2011) · Zbl 1218.92031 |
| [19] | Moghadam, M. E.; Bazilevs, Y.; Hsia, T.; Vignon-Clementel, I.; Marsden, A., A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations, Computational Mechanics, 48, 3, 277-291 (2011) · Zbl 1398.76102 |
| [20] | Whiting, C.; Jansen, K., A stabilized finite element method for the incompressible Navier-Stokes equations using a hierarchical basis, International Journal for Numerical Methods in Fluids, 35, 1, 93-116 (2001) · Zbl 0990.76048 |
| [21] | Franca, L.; Frey, S., Stabilized finite element methods: II. The incompressible Navier-Stokes equations, Computer Methods in Applied Mechanics and Engineering, 99, 2-3, 209-233 (1992) · Zbl 0765.76048 |
| [22] | Brooks, A.; Hughes, T., Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations, Computer Methods in Applied Mechanics and Engineering, 32, 1-3, 199-259 (1982) · Zbl 0497.76041 |
| [23] | Bazilevs, Y.; Calo, V.; Cottrell, J.; Hughes, T.; Reali, A.; Scovazzi, G., Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering, 197, 1-4, 173-201 (2007) · Zbl 1169.76352 |
| [24] | Jansen, K.; Whiting, C.; Hulbert, G., A generalized-[alpha] method for integrating the filtered Navier-Stokes equations with a stabilized finite element method, Computer Methods in Applied Mechanics and Engineering, 190, 3-4, 305-319 (2000) · Zbl 0973.76048 |
| [25] | Shakib, F.; Hughes, T.; Johan, Z., A multi-element group preconditioned gmres algorithm for nonsymmetric systems arising in finite element analysis, Computer Methods in Applied Mechanics and Engineering, 75, 1-3, 415-456 (1989) · Zbl 0687.76065 |
| [27] | Kung, E.; Les, A.; Figueroa, C.; Medina, F.; Arcaute, K.; Wicker, R.; McConnell, M.; Taylor, C., In vitro validation of finite element analysis of blood flow in deformable models, Annals of Biomedical Engineering, 39, 1947-1960 (2011) |
| [28] | Kundu, P.; Cohen, I., Fluid Mechanics (2008), Academic press |
| [29] | Kim, H.; Figueroa, C.; Hughes, T.; Jansen, K.; Taylor, C., Augmented lagrangian method for constraining the shape of velocity profiles at outlet boundaries for three-dimensional finite element simulations of blood flow, Computer Methods in Applied Mechanics and Engineering, 198, 45-46, 3551-3566 (2009) · Zbl 1229.76118 |
| [30] | de Zelicourt, D.; Ge, L.; Wang, C.; Sotiropoulos, F.; Gilmanov, A.; Yoganathan, A., Flow simulations in arbitrarily complex cardiovascular anatomies-an unstructured cartesian grid approach, Computers and Fluids, 38, 9, 1749-1762 (2009) · Zbl 1177.76481 |
| [31] | Bazilevs, Y.; Gohean, J.; Hughes, T.; Moser, R.; Zhang, Y., Patient-specific isogeometric fluid-structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device, Computer Methods in Applied Mechanics and Engineering, 198, 45-46, 3534-3550 (2009) · Zbl 1229.74096 |
| [32] | Migliavacca, F.; Pennati, G.; Dubini, G.; Fumero, R.; Pietrabissa, R.; Urcelay, G.; Bove, E.; Hsia, T.; de Leval, M., Modeling of the Norwood circulation: effects of shunt size, vascular resistances, and heart rate, American Journal of Physiology-Heart and Circulatory Physiology, 280, H2076-H2086 (2001) |
| [33] | Sankaran, S.; Moghadam, M. E.; Kahn, A.; Tseng, E.; Guccione, J.; Marsden, A., Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery, Annals of Biomedical Engineering (2012) |
| [34] | Lagana, K.; Balossino, R.; Migliavacca, F.; Pennati, G.; Bove, E.; de Leval, M.; Dubini, G., Multiscale modeling of the cardiovascular system: application to the study of pulmonary and coronary perfusions in the univentricular circulation, Journal of Biomechanics, 38, 5, 1129-1141 (2005) |
| [35] | Migliavacca, F.; Balossino, R.; Pennati, G.; Dubini, G.; Hsia, T.; de Leval, M.; Bove, E., Multiscale modelling in biofluiddynamics: application to reconstructive paediatric cardiac surgery, Journal of Biomechanics, 39, 6, 1010-1020 (2006) |
| [36] | Baretta, A.; Corsini, C.; Yang, W.; Vignon-Clementel, I.; Marsden, A.; Feinstein, J.; Hsia, T.; Dubini, G.; Migliavacca, F.; Mocha, Virtual surgeries in patients with congenital heart disease: a multi-scale modelling test case, Philosophical Transactions. A. Mathematical, Physical, and Engineering Sciences, 369, 1954, 4316-4330 (2011) |
| [37] | Moghadam, M. E.; Migliavacca, F.; Vignon-Clementel, I.; Hsia, T.; Marsden, A., Optimization of shunt placement for the norwood surgery using multi-domain modeling, Journal of Biomechanical Engineering, 134, 5, 051002 (2012) |
| [38] | Bazilevs, Y.; Calo, V.; Zhang, Y.; Hughes, T., Isogeometric fluid structure interaction analysis with applications to arterial blood flow, Computational Mechanics, 38, 310-322 (2006) · Zbl 1161.74020 |
| [39] | Figueroa, C.; Vignon-Clementel, I.; Jansen, K.; Hughes, T.; Taylor, C., A coupled momentum method for modeling blood flow in three-dimensional deformable arteries, Computer Methods in Applied Mechanics and Engineering, 195, 41-43, 5685-5706 (2006) · Zbl 1126.76029 |
| [40] | Olufsen, M.; Tran, H.; Ottesen, J.; Program, R.; Lipsitz, L.; Novak, V., Modeling baroreflex regulation of heart rate during orthostatic stress, American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, 291, 5, R1355-R1368 (2006) |
| [41] | Pope, S.; Ellwein, L.; Zapata, C.; Novak, V.; Kelley, M.; Olufsen, M., Estimation and identification of parameters in a lumped cerebrovascular model, Mathematical Biosciences Engineering, 6, 1, 93-115 (2009) · Zbl 1159.92025 |
| [42] | Troianowski, G.; Taylor, C.; Feinstein, J.; Vignon-Clementel, I., Three-dimensional simulations in glenn patients: clinically based boundary conditions, hemodynamic results and sensitivity to input data, Journal of Biomechanical Engineering, 133, 11, 111006 (2011) |
| [43] | Pennati, G.; Corsini, C.; Cosentino, D.; Hsia, T.; Luisi, V.; Dubini, G.; Migliavacca, F., Boundary conditions of patient-specific fluid dynamics modelling of cavopulmonary connections: possible adaptation of pulmonary resistances results in a critical issue for a virtual surgical planning, Interface Focus, 1, 297-307 (2011) |
| [44] | Bertoglio, C.; Moireau, P.; Gerbeau, J., Sequential parameter estimation for fluid structure problems: Application to hemodynamics, International Journal of Numerical Methods in Biomedical Engineering, 28, 4, 434-455 (2012) |
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.
