1
|
Pfaller MR, Pham J, Verma A, Pegolotti L, Wilson NM, Parker DW, Yang W, Marsden AL. Automated generation of 0D and 1D reduced-order models of patient-specific blood flow. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2022; 38:e3639. [PMID: 35875875 PMCID: PMC9561079 DOI: 10.1002/cnm.3639] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Revised: 05/24/2022] [Accepted: 07/19/2022] [Indexed: 06/13/2023]
Abstract
Three-dimensional (3D) cardiovascular fluid dynamics simulations typically require hours to days of computing time on a high-performance computing cluster. One-dimensional (1D) and lumped-parameter zero-dimensional (0D) models show great promise for accurately predicting blood bulk flow and pressure waveforms with only a fraction of the cost. They can also accelerate uncertainty quantification, optimization, and design parameterization studies. Despite several prior studies generating 1D and 0D models and comparing them to 3D solutions, these were typically limited to either 1D or 0D and a singular category of vascular anatomies. This work proposes a fully automated and openly available framework to generate and simulate 1D and 0D models from 3D patient-specific geometries, automatically detecting vessel junctions and stenosis segments. Our only input is the 3D geometry; we do not use any prior knowledge from 3D simulations. All computational tools presented in this work are implemented in the open-source software platform SimVascular. We demonstrate the reduced-order approximation quality against rigid-wall 3D solutions in a comprehensive comparison with N = 72 publicly available models from various anatomies, vessel types, and disease conditions. Relative average approximation errors of flows and pressures typically ranged from 1% to 10% for both 1D and 0D models, measured at the outlets of terminal vessel branches. In general, 0D model errors were only slightly higher than 1D model errors despite requiring only a third of the 1D runtime. Automatically generated ROMs can significantly speed up model development and shift the computational load from high-performance machines to personal computers.
Collapse
Affiliation(s)
- Martin R. Pfaller
- Pediatric Cardiology, Stanford University, CA, USA
- Institute for Computational and Mathematical Engineering, Stanford University, CA, USA
- Cardiovascular Institute, Stanford University, CA, USA
| | - Jonathan Pham
- Mechanical Engineering, Stanford University, CA, USA
| | | | - Luca Pegolotti
- Pediatric Cardiology, Stanford University, CA, USA
- Institute for Computational and Mathematical Engineering, Stanford University, CA, USA
| | | | | | | | - Alison L. Marsden
- Pediatric Cardiology, Stanford University, CA, USA
- Institute for Computational and Mathematical Engineering, Stanford University, CA, USA
- Cardiovascular Institute, Stanford University, CA, USA
- Bioengineering, Stanford University, CA, USA
| |
Collapse
|
2
|
Mansilla Alvarez LA, Bulant CA, Ares GD, Feijóo RA, Blanco PJ. Feasibility of coronary blood flow simulations using mid-fidelity numeric and geometric models. Biomech Model Mechanobiol 2022; 21:317-334. [PMID: 35001231 DOI: 10.1007/s10237-021-01536-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 11/23/2021] [Indexed: 11/25/2022]
Abstract
The fractional flow reserve index (FFR) is currently used as a gold standard to quantify coronary stenosis's functional relevance. Due to its highly invasive nature, the development of noninvasive surrogates based on simulations has drawn much attention in recent years, emphasizing efficient strategies that enable translational research. The focus of this work is twofold. First, to assess the feasibility of using a mid-fidelity numerical strategy (transversally enriched pipe element method, TEPEM), placed between low- and high-fidelity models, for the estimation of flow-related quantities, such as FFR and wall shear stress (WSS). Low-fidelity models, as zero- or one-dimensional models, are computationally inexpensive but in detriment of poorer spatially detailed predictions. On the other hand, high-fidelity models, such as classical three-dimensional numerical approximations, can provide detailed predictions but their transition to clinical application is prohibitive due to high computational costs. As a second goal, we quantify the impact of the length of lateral branches in the blood flow through the interrogated vessel of interest to further reduce the computational burden. Both studies are addressed considering a cohort of 17 coronary geometries. A total of 20 locations were selected to estimate the FFR index for a wide range of Coronary Flow Reserve (CFR) scenarios. Numerical results suggest that the mid-fidelity TEPEM model is a reliable approach for the efficient estimation of the FFR index and WSS, with an error in the order of [Formula: see text] and [Formula: see text], respectively, when compared to the high-fidelity prediction. Moreover, such mid-fidelity models require much less computational resources, in compliance with infrastructure frequently available in the clinic, by achieving a speedup between 30 and 60 times compared to a conventional finite element approach. Also, we show that shortening peripheral branches does not introduce considerable perturbations either in the flow patterns, in the wall shear stress, or the pressure drop. Comparing the different geometric models, the error in the estimation of FFR index and WSS is reduced to less than [Formula: see text] and [Formula: see text], respectively.
Collapse
Affiliation(s)
- L A Mansilla Alvarez
- National Laboratory for Scientific Computing, LNCC/MCTI, Av. Getúlio Vargas, 333, Petrópolis, RJ, 25651-075, Brazil. .,National Institute of Science and Technology in Medicine Assisted by Scientific Computing, INCT-MACC, Petrópolis, Brasil.
| | - C A Bulant
- National Scientific and Technical Research Council, CONCITEC and Pladema Institute, National University of the Center of the Buenos Aires Province, Tandil, Argentina.,National Institute of Science and Technology in Medicine Assisted by Scientific Computing, INCT-MACC, Petrópolis, Brasil
| | - G D Ares
- National Scientific and Technical Research Council, CONCITEC, Universidad Nacional del Mar del Plata, UNMdP, Tandil, Argentina.,National Institute of Science and Technology in Medicine Assisted by Scientific Computing, INCT-MACC, Petrópolis, Brasil
| | - R A Feijóo
- National Laboratory for Scientific Computing, LNCC/MCTI, Av. Getúlio Vargas, 333, Petrópolis, RJ, 25651-075, Brazil.,National Institute of Science and Technology in Medicine Assisted by Scientific Computing, INCT-MACC, Petrópolis, Brasil
| | - P J Blanco
- National Laboratory for Scientific Computing, LNCC/MCTI, Av. Getúlio Vargas, 333, Petrópolis, RJ, 25651-075, Brazil.,National Institute of Science and Technology in Medicine Assisted by Scientific Computing, INCT-MACC, Petrópolis, Brasil
| |
Collapse
|
3
|
Blanco PJ, dos Santos GHV, Bulant CA, Alvarez AM, Oliveira FA, Cunha-Lima G, Lemos PA. Scaling laws and the left main coronary artery bifurcation. A combination of geometric and simulation analyses. Med Eng Phys 2022; 99:103701. [DOI: 10.1016/j.medengphy.2021.08.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 08/25/2021] [Accepted: 08/31/2021] [Indexed: 11/16/2022]
|
4
|
Pegolotti L, Pfaller MR, Marsden AL, Deparis S. Model order reduction of flow based on a modular geometrical approximation of blood vessels. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2021; 380:113762. [PMID: 34176992 PMCID: PMC8232546 DOI: 10.1016/j.cma.2021.113762] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We are interested in a reduced order method for the efficient simulation of blood flow in arteries. The blood dynamics is modeled by means of the incompressible Navier-Stokes equations. Our algorithm is based on an approximated domain-decomposition of the target geometry into a number of subdomains obtained from the parametrized deformation of geometrical building blocks (e.g., straight tubes and model bifurcations). On each of these building blocks, we build a set of spectral functions by Proper Orthogonal Decomposition of a large number of snapshots of finite element solutions (offline phase). The global solution of the Navier-Stokes equations on a target geometry is then found by coupling linear combinations of these local basis functions by means of spectral Lagrange multipliers (online phase). Being that the number of reduced degrees of freedom is considerably smaller than their finite element counterpart, this approach allows us to significantly decrease the size of the linear system to be solved in each iteration of the Newton-Raphson algorithm. We achieve large speedups with respect to the full order simulation (in our numerical experiments, the gain is at least of one order of magnitude and grows inversely with respect to the reduced basis size), whilst still retaining satisfactory accuracy for most cardiovascular simulations.
Collapse
Affiliation(s)
- Luca Pegolotti
- SCI-SB-SD, Institute of Mathematics, École Polytechnique Fédérale de Lausanne, Station 8, EPFL, CH–1015 Lausanne, Switzerland
| | - Martin R. Pfaller
- Department of Pediatrics (Cardiology), Bioengineering, Stanford University, Clark Center E1.3, 318 Campus Drive, Stanford, CA 94305, USA
| | - Alison L. Marsden
- Department of Pediatrics (Cardiology), Bioengineering, Stanford University, Clark Center E1.3, 318 Campus Drive, Stanford, CA 94305, USA
| | - Simone Deparis
- SCI-SB-SD, Institute of Mathematics, École Polytechnique Fédérale de Lausanne, Station 8, EPFL, CH–1015 Lausanne, Switzerland
| |
Collapse
|
5
|
An Automated Workflow for Hemodynamic Computations in Cerebral Aneurysms. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2020; 2020:5954617. [PMID: 32655681 PMCID: PMC7317611 DOI: 10.1155/2020/5954617] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Revised: 04/03/2020] [Accepted: 04/28/2020] [Indexed: 01/06/2023]
Abstract
In recent years, computational fluid dynamics (CFD) has become a valuable tool for investigating hemodynamics in cerebral aneurysms. CFD provides flow-related quantities, which have been shown to have a potential impact on aneurysm growth and risk of rupture. However, the adoption of CFD tools in clinical settings is currently limited by the high computational cost and the engineering expertise required for employing these tools, e.g., for mesh generation, appropriate choice of spatial and temporal resolution, and of boundary conditions. Herein, we address these challenges by introducing a practical and robust methodology, focusing on computational performance and minimizing user interaction through automated parameter selection. We propose a fully automated pipeline that covers the steps from a patient-specific anatomical model to results, based on a fast, graphics processing unit- (GPU-) accelerated CFD solver and a parameter selection methodology. We use a reduced order model to compute the initial estimates of the spatial and temporal resolutions and an iterative approach that further adjusts the resolution during the simulation without user interaction. The pipeline and the solver are validated based on previously published results, and by comparing the results obtained for 20 cerebral aneurysm cases with those generated by a state-of-the-art commercial solver (Ansys CFX, Canonsburg PA). The automatically selected spatial and temporal resolutions lead to results which closely agree with the state-of-the-art, with an average relative difference of only 2%. Due to the GPU-based parallelization, simulations are computationally efficient, with a median computation time of 40 minutes per simulation.
Collapse
|
6
|
Fleeter CM, Geraci G, Schiavazzi DE, Kahn AM, Marsden AL. Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2020; 365:113030. [PMID: 32336811 PMCID: PMC7182133 DOI: 10.1016/j.cma.2020.113030] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. We propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo (MLMF) estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models (zero- and one-dimensional) to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these two low-fidelity model alternatives and our high-fidelity three-dimensional models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. Our goal is to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to single fidelity Monte Carlo estimators and to a multilevel Monte Carlo approach based only on three-dimensional simulations, but leveraging multiple spatial resolutions. We demonstrate significant, on the order of 10 to 100 times, reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities such as outlet pressure and flow show larger reductions than local quantities, such as those relating to wall shear stress, as the latter rely more heavily on the highest fidelity model evaluations. Similarly, healthy models show larger reductions than diseased models. In all cases, our workflow coupling Dakota's MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.
Collapse
Affiliation(s)
- Casey M. Fleeter
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA, USA
| | - Gianluca Geraci
- Center for Computing Research, Sandia National Laboratories, Albuquerque, NM, USA
| | - Daniele E. Schiavazzi
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, USA
| | - Andrew M. Kahn
- Department of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Alison L. Marsden
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA, USA
- Departments of Pediatrics and Bioengineering, Stanford University, Stanford, CA, USA
| |
Collapse
|