Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
169 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
45 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

A new stabilized time-spectral finite element solver for fast simulation of blood flow (2311.09465v3)

Published 16 Nov 2023 in math.NA, cs.NA, and physics.flu-dyn

Abstract: The increasing application of cardiorespiratory simulations for diagnosis and surgical planning necessitates the development of computational methods significantly faster than the current technology. To achieve this objective, we leverage the time-periodic nature of these flows by discretizing equations in the frequency domain instead of the time domain. This approach markedly reduces the size of the discrete problem and, consequently, the simulation cost. With this motivation, we introduce a finite element method for simulating time-periodic flows that are physically stable. The proposed time-spectral method is formulated by augmenting the baseline Galerkin's method with a least-squares penalty term that is weighed by a positive-definite stabilization matrix. An error estimate is established for the convective-diffusive system, showing that the proposed method emulates the behavior of existing standard time methods including optimal convergence rate in diffusive regimes and stability in strong convection. This method is tested on a patient-specific Fontan model at nominal Reynolds and Womersley numbers of 500 and 10, respectively, demonstrating its ability to replicate conventional time simulation results using as few as 7 modes at 11% of the computational cost. Owing to its higher local-to-processor computation density, the proposed method also exhibits improved parallel scalability for fast simulation of time-critical applications.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (104)
  1. Charles A Taylor and CA Figueroa. Patient-specific modeling of cardiovascular mechanics. Annual review of biomedical engineering, 11:109–134, 2009.
  2. Constrained optimization of an idealized Y-shaped baffle for the Fontan surgery at rest and exercise. Computer Methods in Applied Mechanics and Engineering, 199(33-36):2135–2149, 2010.
  3. A computational framework for derivative-free optimization of cardiovascular geometries. Computer Methods in Applied Mechanics and Engineering, 197(21-24):1890–1905, 2008.
  4. Introduction of a new optimized total cavopulmonary connection. The Annals of Thoracic Surgery, 83(6):2182–2190, 2007.
  5. The surgical treatment of malformations of the heart. Journal of the American Medical Association, 128(3):189–202, 1945.
  6. Hypoplastic left heart syndrome: experience with palliative surgery. The American Journal of Cardiology, 45(1):87–91, 1980.
  7. The assisted bidirectional Glenn: a novel surgical approach for first stage single ventricle heart palliation. The Journal of Thoracic and Cardiovascular Surgery, 149(3):699–705, 2015.
  8. Optimization of shunt placement for the Norwood surgery using multi-domain modeling. Journal of Biomechanical Engineering, 134(5):051002, 2012.
  9. Optimization of the assisted bidirectional Glenn procedure for first stage single ventricle repair. World Journal for Pediatric and Congenital Heart Surgery, 9(2):157–170, 2018.
  10. Comparison of coronary computed tomography angiography, fractional flow reserve, and perfusion imaging for ischemia diagnosis. Journal of the American College of Cardiology, 73(2):161–173, 2019.
  11. Predicting the physiological effect of revascularization in serially diseased coronary arteries: clinical validation of a novel CT coronary angiography–based technique. Circulation: cardiovascular interventions, 12(2):e007577, 2019.
  12. Thomas JR Hughes. Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations. International journal for numerical methods in fluids, 7(11):1261–1275, 1987.
  13. Stabilized finite element methods: II. the incompressible Navier-Stokes equations. Computer Methods in Applied Mechanics and Engineering, 99(2-3):209–233, 1992.
  14. Thomas JR Hughes. Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods. Computer methods in applied mechanics and engineering, 127(1-4):387–401, 1995.
  15. G Hauke and TJR Hughes. A unified approach to compressible and incompressible flows. Computer Methods in Applied Mechanics and Engineering, 113(3-4):389–395, 1994.
  16. Ramon Codina. Stabilization of incompressibility and convection through orthogonal sub-scales in finite element methods. Computer methods in applied mechanics and engineering, 190(13-14):1579–1599, 2000.
  17. Stabilized methods for compressible flows. Journal of Scientific Computing, 43:343–368, 2010.
  18. Thomas JR Hughes. A multidimentional upwind scheme with no crosswind diffusion. Finite element methods for convection dominated flows, AMD 34, 1979.
  19. 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.
  20. Olga Alexandrovna Ladyzhenskaya. The mathematical theory of viscous incompressible flow, volume 2. Gordon and Breach New York, 1969.
  21. Ivo Babuska. Error-bounds for finite element method. Numerische Mathematik, 16(4):322–333, 1971.
  22. Franco Brezzi. On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. ESAIM: Mathematical Modelling and Numerical Analysis-Modélisation Mathématique et Analyse Numérique, 8(R2):129–151, 1974.
  23. A new finite element formulation for computational fluid dynamics: V. circumventing the Babuška-Brezzi condition: A stable Petrov-Galerkin formulation of the Stokes problem accommodating equal-order interpolations. Computer Methods in Applied Mechanics and Engineering, 59(1):85–99, 1986.
  24. Tayfun E Tezduyar. Stabilized finite element formulations for incompressible flow computations. In Advances in applied mechanics, volume 28, pages 1–44. Elsevier, 1991.
  25. A new finite element formulation for computational fluid dynamics: III. the generalized streamline operator for multidimensional advective-diffusive systems. Computer methods in applied mechanics and engineering, 58(3):305–328, 1986.
  26. A new finite element formulation for computational fluid dynamics: X. the compressible Euler and Navier-Stokes equations. Computer Methods in Applied Mechanics and Engineering, 89(1-3):141–219, 1991.
  27. A comparative study of different sets of variables for solving compressible and incompressible flows. Computer methods in applied mechanics and engineering, 153(1-2):1–44, 1998.
  28. A new finite element formulation for computational fluid dynamics: VIII. the Galerkin/least-squares method for advective-diffusive equations. Computer methods in applied mechanics and engineering, 73(2):173–189, 1989.
  29. Farzin Shakib. Finite element analysis of the compressible Euler and Navier-Stokes equations. Stanford University, 1989.
  30. A scalable spectral Stokes solver for simulation of time-periodic flows in complex geometries. Journal of Computational Physics, 445:110601, 2021.
  31. A stabilized formulation for the solution of the incompressible unsteady stokes equations in the frequency domain. Journal of Computational Physics, 473:111736, 2023.
  32. Stabilized finite element methods for the time-spectral convection-diffusion equation. arXiv preprint arXiv:2305.12038, 2023.
  33. Application of a non-linear frequency domain solver to the Euler and Navier-Stokes equations. In 40th AIAA aerospace sciences meeting & exhibit, page 120, 2002.
  34. Matthew Scott McMullen. The application of non-linear frequency domain methods to the Euler and Navier-Stokes equations. Citeseer, 2003.
  35. Time spectral method for periodic unsteady computations over two- and three-dimensional bodies. In 43rd AIAA aerospace sciences meeting and exhibit, page 1220, 2005.
  36. Acceleration of convergence to a periodic steady state in turbomachinery flows. In 39th aerospace sciences meeting and exhibit, page 152, 2001.
  37. Three-dimensional unsteady multi-stage turbomachinery simulations using the harmonic balance technique. In 45th AIAA Aerospace Sciences Meeting and Exhibit, page 892, 2007.
  38. Analysis of rotor-rotor and stator-stator interferences in multi-stage turbomachines. In Turbo Expo: Power for Land, Sea, and Air, volume 3610, pages 287–298, 2002.
  39. A time-domain harmonic balance method for rotor/stator interactions. J. Turbomach., 134(1):011001, 2012.
  40. Comparison of parallel time-periodic Navier-Stokes solvers. In International Conference on Parallel Processing and Applied Mathematics, pages 57–67. Springer, 2017.
  41. A parallel Navier–Stokes solver using spectral discretisation in time. International journal of computational fluid dynamics, 30(7-10):489–494, 2016.
  42. Harmonic balance methods applied to computational fluid dynamics problems. International Journal of Computational Fluid Dynamics, 27(2):52–67, 2013.
  43. Comparison of two algorithms for harmonic balance and path continuation. Mechanical Systems and Signal Processing, 136:106503, 2020.
  44. Computation of unsteady nonlinear flows in cascades using a harmonic balance technique. AIAA Journal, 40(5):879–886, 2002.
  45. Nonlinear input/output analysis: application to boundary layer transition. Journal of Fluid Mechanics, 911:A15, 2021.
  46. The clinical assessment of intraventricular flows. Annual Review of Fluid Mechanics, 47:315–342, 2015.
  47. Application of large-eddy simulation to the study of pulsatile flow in a modeled arterial stenosis. Journal of Biomechanical Engineering, 123(4):325–332, 2001.
  48. Simvascular: An open source pipeline for cardiovascular simulation. Annals of Biomedical Engineering, 45(3):525–541, 2017.
  49. Characterization of the ejector pump performance for the assisted bidirectional Glenn procedure. Fluids, 7(1):31, 2022.
  50. An efficient assisted bidirectional glenn design with lowered superior vena cava pressure for stage-one single ventricle patients. Journal of Biomechanical Engineering, 143(7):071008, 2021.
  51. Multiscale modelling as a tool to prescribe realistic boundary conditions for the study of surgical procedures. Biorheology, 39:359–364, 2002.
  52. Numerical treatment of defective boundary conditions for the Navier-Stokes equations. SIAM Journal on Numerical Analysis, 40:376–401, 2002.
  53. 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.
  54. 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.
  55. A modular numerical method for implicit 0D/3D coupling in cardiovascular finite element simulations. Journal of Computational Physics, 224:63–79, 2013.
  56. Isogeometric fluid-structure interaction: theory, algorithms, and computations. Computational Mechanics, 43:3–37, 2008.
  57. Computational fluid-structure interaction: methods and application to a total cavopulmonary connection. Computational Mechanics, 45:77–89, 2009.
  58. Fluid–structure interaction simulation of pulsatile ventricular assist devices. Computational Mechanics, 52(5):971–981, 2013.
  59. An immersogeometric variational framework for fluid–structure interaction: Application to bioprosthetic heart valves. Computer methods in applied mechanics and engineering, 284:1005–1053, 2015.
  60. Simulations reveal adverse hemodynamics in patients with multiple systemic to pulmonary shunts. Journal of Biomechanical Engineering, 137(3):031001–031001–12, 2015.
  61. A longitudinal comparison of hemodynamics and intraluminal thrombus deposition in abdominal aortic aneurysms. American Journal of Physiology-Heart and Circulatory Physiology, 307(12):H1786–H1795, 2014.
  62. A non-discrete method for computation of residence time in fluid mechanics simulations. Physics of Fluids, 25(11):110802–21, 2013.
  63. Potential fluid mechanic pathways of platelet activation. Biomechanics and modeling in mechanobiology, 12:467–474, 2013.
  64. A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows. Journal of Computational Physics, 461:111204, 2022.
  65. Coupled simulation of hemodynamics and vascular growth and remodeling in a subject-specific geometry. Annals of biomedical engineering, 43(7):1543–1554, 2015.
  66. Abhay Bangalore Ramachandra. Towards a Fluid-Structure-Growth and Remodeling Framework to Simulate Vein Graft Failure Post Coronary Artery Bypass Surgery. PhD thesis, UC San Diego, 2017.
  67. Computational simulations of hemodynamic changes within thoracic, coronary, and cerebral arteries following early wall remodeling in response to distal aortic coarctation. Biomechanics and modeling in mechanobiology, 12(1):79–93, 2013.
  68. 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.
  69. Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers? SIAM Journal on numerical analysis, 34(6):2392–2423, 1997.
  70. Design and analysis of finite element methods for the Helmholtz equation in exterior domains. Appl. Mech. Rev., 43(5):366–373, 1990.
  71. A Galerkin least-squares finite element method for the two-dimensional Helmholtz equation. International Journal for numerical methods in engineering, 38(3):371–397, 1995.
  72. Tridiagonal Toeplitz matrices: properties and novel applications. Numerical linear algebra with applications, 20(2):302–326, 2013.
  73. Finite element analysis of convection dominated reaction–diffusion problems. Applied Numerical Mathematics, 48(2):205–222, 2004.
  74. A discourse on the stability conditions for mixed finite element formulations. Computer methods in applied mechanics and engineering, 82(1-3):27–57, 1990.
  75. Error analysis of Galerkin least squares methods for the elasticity equations. SIAM Journal on Numerical Analysis, 28(6):1680–1697, 1991.
  76. Convergence analyses of Galerkin least-squares methods for symmetric advective-diffusive forms of the Stokes and incompressible Navier-Stokes equations. Computer Methods in Applied Mechanics and Engineering, 105(2):285–298, 1993.
  77. On the stability of finite-element discretizations of convection–diffusion–reaction equations. IMA Journal of numerical analysis, 31(1):147–164, 2011.
  78. A new finite element formulation for computational fluid dynamics: II. beyond SUPG. Computer Methods in Applied Mechanics and Engineering, 54(3):341–355, 1986.
  79. Finite element methods for linear hyperbolic problems. Computer methods in applied mechanics and engineering, 45:285–312, 1984.
  80. 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.
  81. John R Womersley. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. The Journal of Physiology, 127(3):553, 1955.
  82. The finite element method with Lagrange multipliers on the boundary: circumventing the Babuška-Brezzi condition. Computer Methods in Applied Mechanics and Engineering, 85(1):109–128, 1991.
  83. Variational multiscale analysis: the fine-scale Green’s function, projection, optimization, localization, and stabilized methods. SIAM Journal on Numerical Analysis, 45(2):539–557, 2007.
  84. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal on scientific and statistical computing, 7(3):856–869, 1986.
  85. ISAAC Fried. Bounds on the extremal eigenvalues of the finite element stiffness and mass matrices and their spectral condition number. Journal of Sound and Vibration, 22(4):407–418, 1972.
  86. A generalized-α𝛼\alphaitalic_α 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.
  87. A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations. Computational Mechanics, 48(3):277–291, 2011.
  88. Benchmark problems for numerical treatment of backflow at open boundaries. International Journal for Numerical Methods in Biomedical Engineering, 34(2):e2918, 2018.
  89. 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.
  90. The vascular model repository: A public resource of medical imaging data and blood flow simulation results. Journal of Medical Devices, 7(4):040923, 2013.
  91. A new multiparameter approach to computational simulation for fontan assessment and redesign. Congenital Heart Disease, 5(2):104–117, 2010.
  92. Pulmonary arteriovenous malformations after cavopulmonary anastomosis. The Annals of thoracic surgery, 76(5):1759–1766, 2003.
  93. Optimization of a Y-graft design for improved hepatic flow distribution in the fontan circulation. Journal of biomechanical engineering, 135(1):011002, 2013.
  94. Predicting hemodynamic performance of Fontan operation for Glenn physiology using computational fluid dynamics: Ten patient-specific cases. Journal of clinical images and medical case reports, 3(6), 2022.
  95. M. Esmaily. Multi-physics finite element solver (mupfes), May 2014.
  96. A bi-partitioned iterative algorithm for solving linear systems arising from incompressible flow problems. Computer Methods in Applied Mechanics and Engineering, 286(1):40–62, 2015.
  97. 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.
  98. G Rigas and PJ Schmid. Data-driven closure of the harmonic-balanced navier-stokes equations in the frequency domain. Proceedings of the Summer Program, pages 67–76, 2022.
  99. Impact of data distribution on the parallel performance of iterative linear solvers with emphasis on CFD of incompressible flows. Computational Mechanics, 55(1):93–103, 2015.
  100. Computational fluid-structure interaction: methods and applications. Wiley, 2013.
  101. Improving stability of stabilized and multiscale formulations in flow simulations at small time steps. Computer Methods in Applied Mechanics and Engineering, 199(13-16):828–840, 2010.
  102. Time dependent subscales in the stabilized finite element approximation of incompressible flow problems. Computer Methods in Applied Mechanics and Engineering, 196(21-24):2413–2430, 2007.
  103. A time-consistent stabilized finite element method for fluids with applications to hemodynamics. Scientific Reports, 13(1):19120, 2023.
  104. Residual-Based Large Eddy Simulation with Isogeometric Divergence-Conforming Discretizations, pages 91–130. Springer International Publishing, 2018.
Citations (2)
View on Semantic Scholar

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now