A Time-consistent Stabilized Finite Element Method for Fluids with Applications to Hemodynamics

Overview

Journal Sci Rep

Specialty Science

Date 2023 Nov 5

PMID 37926732

Authors

Dongjie Jia

Mahdi Esmaily

Affiliations

Soon will be listed here.

Abstract

Several finite element methods for simulating incompressible flows rely on the streamline upwind Petrov-Galerkin stabilization (SUPG) term, which is weighted by [Formula: see text]. The conventional formulation of [Formula: see text] includes a constant that depends on the time step size, producing an overall method that becomes exceedingly less accurate as the time step size approaches zero. In practice, such method inconsistency introduces significant error in the solution, especially in cardiovascular simulations, where small time step sizes may be required to resolve multiple scales of the blood flow. To overcome this issue, we propose a consistent method that is based on a new definition of [Formula: see text]. This method, which can be easily implemented on top of an existing streamline upwind Petrov-Galerkin and pressure stabilizing Petrov-Galerkin method, involves the replacement of the time step size in [Formula: see text] with a physical time scale. This time scale is calculated in a simple operation once every time step for the entire computational domain from the ratio of the L-norm of the acceleration and the velocity. The proposed method is compared against the conventional method using four cases: a steady pipe flow, a blood flow through vascular anatomy, an external flow over a square obstacle, and a fluid-structure interaction case involving an oscillatory flexible beam. These numerical experiments, which are performed using linear interpolation functions, show that the proposed formulation eliminates the inconsistency issue associated with the conventional formulation in all cases. While the proposed method is slightly more costly than the conventional method, it significantly reduces the error, particularly at small time step sizes. For the pipe flow where an exact solution is available, we show the conventional method can over-predict the pressure drop by a factor of three. This large error is almost completely eliminated by the proposed formulation, dropping to approximately 1% for all time step sizes and Reynolds numbers considered.

References

Taylor C, Hughes T, Zarins C . Finite element modeling of three-dimensional pulsatile flow in the abdominal aorta: relevance to atherosclerosis. Ann Biomed Eng. 1998; 26(6):975-87. DOI: 10.1114/1.140. View

Verma A, Esmaily M, Shang J, Figliola R, Feinstein J, Hsia T . Optimization of the Assisted Bidirectional Glenn Procedure for First Stage Single Ventricle Repair. World J Pediatr Congenit Heart Surg. 2018; 9(2):157-170. PMC: 6668725. DOI: 10.1177/2150135117745026. View

Esmaily Moghadam M, Migliavacca F, Vignon-Clementel I, Hsia T, Marsden A . Optimization of shunt placement for the Norwood surgery using multi-domain modeling. J Biomech Eng. 2012; 134(5):051002. DOI: 10.1115/1.4006814. View

Jia D, Peroni M, Khalapyan T, Esmaily M . An Efficient Assisted Bidirectional Glenn Design With Lowered Superior Vena Cava Pressure for Stage-One Single Ventricle Patients. J Biomech Eng. 2021; 143(7). DOI: 10.1115/1.4050170. View

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. Ann Biomed Eng. 2010; 38(10):3195-209. DOI: 10.1007/s10439-010-0083-6. View

Shang J, Esmaily M, Verma A, Reinhartz O, Figliola R, Hsia T . Patient-Specific Multiscale Modeling of the Assisted Bidirectional Glenn. Ann Thorac Surg. 2018; 107(4):1232-1239. PMC: 6516076. DOI: 10.1016/j.athoracsur.2018.10.024. View

Arthurs C, Khlebnikov R, Melville A, Marcan M, Gomez A, Dillon-Murphy D . CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation. PLoS Comput Biol. 2021; 17(5):e1008881. PMC: 8148362. DOI: 10.1371/journal.pcbi.1008881. View

Esmaily-Moghadam M, Hsia T, Marsden A . The assisted bidirectional Glenn: a novel surgical approach for first-stage single-ventricle heart palliation. J Thorac Cardiovasc Surg. 2014; 149(3):699-705. PMC: 6089632. DOI: 10.1016/j.jtcvs.2014.10.035. View

Dyverfeldt P, Bissell M, Barker A, Bolger A, Carlhall C, Ebbers T . 4D flow cardiovascular magnetic resonance consensus statement. J Cardiovasc Magn Reson. 2015; 17:72. PMC: 4530492. DOI: 10.1186/s12968-015-0174-5. View

10.

Steinman D, Hoi Y, Fahy P, Morris L, Walsh M, Aristokleous N . Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge. J Biomech Eng. 2013; 135(2):021016. DOI: 10.1115/1.4023382. View

11.

Arbia G, Corsini C, Esmaily Moghadam M, Marsden A, Migliavacca F, Pennati G . Numerical blood flow simulation in surgical corrections: what do we need for an accurate analysis?. J Surg Res. 2013; 186(1):44-55. DOI: 10.1016/j.jss.2013.07.037. View

12.

Sankaran S, Esmaily Moghadam M, Kahn A, Tseng E, Guccione J, Marsden A . Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery. Ann Biomed Eng. 2012; 40(10):2228-42. PMC: 3570226. DOI: 10.1007/s10439-012-0579-3. View

13.

Schulz-Menger J, Bluemke D, Bremerich J, Flamm S, Fogel M, Friedrich M . Standardized image interpretation and post-processing in cardiovascular magnetic resonance - 2020 update : Society for Cardiovascular Magnetic Resonance (SCMR): Board of Trustees Task Force on Standardized Post-Processing. J Cardiovasc Magn Reson. 2020; 22(1):19. PMC: 7066763. DOI: 10.1186/s12968-020-00610-6. View

14.

Updegrove A, Wilson N, Merkow J, Lan H, Marsden A, Shadden S . SimVascular: An Open Source Pipeline for Cardiovascular Simulation. Ann Biomed Eng. 2016; 45(3):525-541. PMC: 6546171. DOI: 10.1007/s10439-016-1762-8. View

15.

Taylor C, Draney M, Ku J, Parker D, Steele B, Wang K . Predictive medicine: computational techniques in therapeutic decision-making. Comput Aided Surg. 1999; 4(5):231-47. DOI: 10.1002/(SICI)1097-0150(1999)4:5<231::AID-IGS1>3.0.CO;2-Z. View

16.

Markl M, Frydrychowicz A, Kozerke S, Hope M, Wieben O . 4D flow MRI. J Magn Reson Imaging. 2012; 36(5):1015-36. DOI: 10.1002/jmri.23632. View

17.

Markl M, Chan F, Alley M, Wedding K, Draney M, Elkins C . Time-resolved three-dimensional phase-contrast MRI. J Magn Reson Imaging. 2003; 17(4):499-506. DOI: 10.1002/jmri.10272. View

18.

Taylor C, Figueroa C . Patient-specific modeling of cardiovascular mechanics. Annu Rev Biomed Eng. 2009; 11:109-34. PMC: 4581431. DOI: 10.1146/annurev.bioeng.10.061807.160521. View

19.

Kamensky D, Hsu M, Schillinger D, Evans J, Aggarwal A, Bazilevs Y . An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves. Comput Methods Appl Mech Eng. 2014; 284:1005-1053. PMC: 4274080. DOI: 10.1016/j.cma.2014.10.040. View

20.

Taylor C, Fonte T, Min J . Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis. J Am Coll Cardiol. 2013; 61(22):2233-41. DOI: 10.1016/j.jacc.2012.11.083. View