| Mesh | $h_{\mathrm{NOR}}$ | $h_{\mathrm{TAN}}$ | Growth rate | Boundary layer height ($H$) | Internal refinement | Total number of elements |
| IM1 | 0.02 | 0.4 | 1.1 | 0.5 | 0.2 | 378,867 |
| IM2 | 0.01 | 0.2 | 1.1 | 0.5 | 0.1 | 1,074,581 |
| IM3 | 0.005 | 0.1 | 1.1 | 0.5 | 0.05 | 4,180,767 |
Cardiovascular disease (CVD) remains one of the leading causes of mortality worldwide. Computational medicine and digital twins hold promise in mitigating the impact and prevalence of CVD. Recent advances in image-based computational methods have enabled the quantification of functional and biologically important metrics that would otherwise be difficult to obtain from the standard of care. However, significant challenges remain due to the manual/semi-automated nature of the processes and the domain expertise required to perform them. This paper addresses these challenges by proposing a novel framework that builds on our recently developed direct point cloud-to-CFD approach using immersogeometric analysis. The proposed method leverages advanced auto-segmentation techniques to extract medically relevant geometries as point clouds, which are then directly used for CFD simulations. The framework is validated using benchmark flow problems with analytical and computational solutions and is subsequently applied to patient-specific images to demonstrate its capabilities. The results highlight the method's ability to facilitate rapid CFD simulations directly on point clouds derived from patient scans, underscoring its potential to accelerate the image-to-simulation pipeline and enable the tractability of cardiovascular digital twins.
| Citation: |
Figure 2. Diagram of the schematic representation of the proposed pipeline for aorta segmentation and 3D mesh reconstruction. The image segmentation architecture utilizes 3D U-Net for segmenting the aorta and its branches. The encoder component of 3D U-Net is responsible for extracting high-dimensional features from the small field of view. These features are then combined with the image gradient obtained from the MRI image within the shape stream
Figure 4. (A) An example of flow within an object, illustrated using a 2D cross-section of a 3D channel represented by a point cloud. The object with boundary $ \Gamma $ is immersed into the domain $ \Omega $. The immersed boundary $ \Gamma $ divides the domain $ \Omega $ into a physical part $ \Omega_{\text{phys}} $, where the flow occurs, and a fictitious part $ \Omega_{\text{fict}} $. (B) Background mesh discretizing the domain $ \Omega $. (C) The faces ($ \mathcal{F}^f $) where the ghost penalty is applied are highlighted in red. These faces are shared by two neighboring elements that lie fully or partially within the physical domain, with at least one of them intersected by the boundary
Figure 5. Background mesh domain containing an aorta model represented by a point cloud. (A) A domain whose boundaries do not coincide with the flow inlet and outlet locations. (B) A clipped domain in which both the domain and the point cloud are trimmed to align with the desired inlet and outlet faces
Figure 12. (A) Velocity magnitude and pressure contours along the center plane. (B) Pressure drop versus distance along the centerline of the pipe before the branches split, with the maximum distance being the bifurcation apex at 13.95 mm. The present work is compared with those from previous studies by Chung [21] and Weddell et al. [96]
Figure 13. Resampled point clouds colored by normalized minimum and maximum depth values for enhanced visualization. From left to right: Case 1 – Adult patient with single ventricle disease who had undergone the Fontan procedure. Case 2 – Pediatric patient with a diseased ventricle and healthy aorta. Case 3 – Healthy lamb aorta
Figure 16. Dimensions of the simulation domain boundary with the point cloud aorta immersed inside. The smaller boxes visualized are Boolean subtractions of the original rectangular domain using an estimated centerline normal to align the background mesh surface with the desired inlet/outlet surfaces. No-slip condition on the arterial wall is enforced weekly using the point cloud
Table 1. Anisotropic mesh refinement metrics for IM1, IM2, and IM3
| Mesh | $h_{\mathrm{NOR}}$ | $h_{\mathrm{TAN}}$ | Growth rate | Boundary layer height ($H$) | Internal refinement | Total number of elements |
| IM1 | 0.02 | 0.4 | 1.1 | 0.5 | 0.2 | 378,867 |
| IM2 | 0.01 | 0.2 | 1.1 | 0.5 | 0.1 | 1,074,581 |
| IM3 | 0.005 | 0.1 | 1.1 | 0.5 | 0.05 | 4,180,767 |
Table 2. Resampling, CT resolution, and percent deviation statistic
| Aorta | Original points | Mean resolution (cm/voxel) | Resampled points | Percent deviation |
| Case 1 | 6,736 | 0.125 | 177,129 | 0.003952 |
| Case 2 | 19,150 | 0.076 | 405,688 | 0.002712 |
| Case 3 | 25,237 | 0.025 | 376,365 | 0.000797 |
| [1] |
F. Alauzet and A. Loseille, A decade of progress on anisotropic mesh adaptation for computational fluid dynamics, Computer-Aided Design, 72 (2016), 13-39.
doi: 10.1016/j.cad.2015.09.005.
|
| [2] |
D. An, P. Du, J.-X. Wang and C. Wang, Hierarchical LoG Bayesian neural network for enhanced aorta segmentation, in IEEE International Symposium on Biomedical Imaging 2025 (ISBI 2025), IEEE, (2025), 1571076368.
|
| [3] |
C. J. Arthurs, R. Khlebnikov, A. Melville, M. Marčan, A. Gomez, D. Dillon-Murphy, F. Cuomo, M. Silva Vieira, J. Schollenberger, S. R. Lynch et al., CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation, PLoS Computational Biology, 17 (2021), e1008881.
doi: 10.1371/journal.pcbi.1008881.
|
| [4] |
W. Bai, H. Suzuki, C. Qin, G. Tarroni, O. Oktay, P. M. Matthews and D. Rueckert, Recurrent neural networks for aortic image sequence segmentation with sparse annotations, in Medical Image Computing and Computer Assisted Intervention-MICCAI 2018: 21st International Conference, Granada, Spain, September 16-20, 2018, Proceedings, Part IV 11, Springer, (2018), 586-594.
doi: 10.1007/978-3-030-00937-3_67.
|
| [5] |
A. Balu, M. R. Rajanna, J. Khristy, F. Xu, A. Krishnamurthy and M.-C. Hsu, Direct immersogeometric fluid flow and heat transfer analysis of objects represented by point clouds, Computer Methods in Applied Mechanics and Engineering, 404 (2023), 115742.
doi: 10.1016/j.cma.2022.115742.
|
| [6] |
C. F. Baumgartner, L. M. Koch, M. Pollefeys and E. Konukoglu, An exploration of 2D and 3D deep learning techniques for cardiac MR image segmentation, in Statistical Atlases and Computational Models of the Heart. ACDC and MMWHS Challenges: 8th International Workshop, STACOM 2017, Held in Conjunction with MICCAI 2017, Quebec City, Canada, September 10-14, 2017, Revised Selected Papers 8, Springer, (2018), 111-119.
doi: 10.1007/978-3-319-75541-0_12.
|
| [7] |
G. Barill, N. G. Dickson, R. Schmidt, D. I. W. Levin and A. Jacobson, Fast winding numbers for soups and clouds, ACM Transactions on Graphics, 37 (2018), 43.
doi: 10.1145/3197517.3201337.
|
| [8] |
Y. Bazilevs and I. Akkerman, Large eddy simulation of turbulent Taylor-Couette flow using isogeometric analysis and the residual-based variational multiscale method, Journal of Computational Physics, 229 (2010), 3402-3414.
doi: 10.1016/j.jcp.2010.01.008.
|
| [9] |
Y. Bazilevs, V. M. Calo, J. A. Cottrell, T. J. R. Hughes, A. Reali and G. Scovazzi, Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering, 197 (2007), 173-201.
doi: 10.1016/j.cma.2007.07.016.
|
| [10] |
Y. Bazilevs, J. R. Gohean, T. J. R. Hughes, R. D. Moser and Y. Zhang, 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 (2009), 3534-3550.
doi: 10.1016/j.cma.2009.04.015.
|
| [11] |
Y. Bazilevs, M.-C. Hsu, D. J. Benson, S. Sankaran and A. L. Marsden, Computational fluid–structure interaction: Methods and application to a total cavopulmonary connection, Computational Mechanics, 45 (2009), 77-89.
doi: 10.1007/s00466-009-0419-y.
|
| [12] |
Y. Bazilevs and T. J. R. Hughes, Weak imposition of Dirichlet boundary conditions in fluid mechanics, Computers and Fluids, 36 (2007), 12-26.
doi: 10.1016/j.compfluid.2005.07.012.
|
| [13] |
Y. Bazilevs, C. Michler, V. M. Calo and T. J. R. Hughes, Weak Dirichlet boundary conditions for wall-bounded turbulent flows, Computer Methods in Applied Mechanics and Engineering, 196 (2007), 4853-4862.
doi: 10.1016/j.cma.2007.06.026.
|
| [14] |
E. Burman, Ghost penalty, Comptes Rendus Mathematique, 348 (2010), 1217-1220.
doi: 10.1016/j.crma.2010.10.006.
|
| [15] |
E. Burman, S. Claus, P. Hansbo, M. G. Larson and A. Massing, CutFEM: Discretizing geometry and partial differential equations, International Journal for Numerical Methods in Engineering, 104 (2015), 472-501.
doi: 10.1002/nme.4823.
|
| [16] |
E. Burman and P. Hansbo, Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method, Applied Numerical Mathematics, 62 (2012), 328-341.
doi: 10.1016/j.apnum.2011.01.008.
|
| [17] |
F. Cazals and M. Pouget, Estimating differential quantities using polynomial fitting of osculating jets, Computer Aided Geometric Design, 22 (2005), 121-146.
doi: 10.1016/j.cagd.2004.09.004.
|
| [18] |
F. Cazals and M. Pouget, Algorithm 889: Jet_fitting_3: A generic C++ package for estimating the differential properties on sampled surfaces via polynomial fitting, ACM Transactions on Mathematical Software, 35 (2008), 24.
doi: 10.1145/1391989.1404582.
|
| [19] |
L.-C. Chen, G. Papandreou, I. Kokkinos, K. Murphy and A. L. Yuille, DeepLab: Semantic image segmentation with deep convolutional nets, atrous convolution, and fully connected CRFs, IEEE Transactions on Pattern Analysis and Machine Intelligence, 40 (2018), 834-848.
doi: 10.1109/TPAMI.2017.2699184.
|
| [20] |
C. Chen, C. Qin, H. Qiu, G. Tarroni, J. Duan, W. Bai and D. Rueckert, Deep learning for cardiac image segmentation: A review, Frontiers in Cardiovascular Medicine, 7 (2020), 25.
doi: 10.3389/fcvm.2020.00025.
|
| [21] |
B. J. Chung, Studies of Blood Flow in Arterial Bifurcations: From Influence of Hemodynamics on Endothelial Cell Response to Vessel Wall Mechanics, PhD thesis, University of Pittsburgh, 2004.
|
| [22] |
Ö. Çiçek, A. Abdulkadir, S. S. Lienkamp, T. Brox and O. Ronneberger, 3D U-Net: Learning dense volumetric segmentation from sparse annotation, in Medical Image Computing and Computer-Assisted Intervention, (2016), 424-432.
doi: 10.1007/978-3-319-46723-8_49.
|
| [23] |
G. Coorey, G. A. Figtree, D. F. Fletcher, V. J. Snelson, S. T. Vernon, D. Winlaw, S. M. Grieve, A. McEwan, J. Y. H. Yang, P. Qian et al., The health digital twin to tackle cardiovascular disease—a review of an emerging interdisciplinary field, NPJ Digital Medicine, 5 (2022), 126.
doi: 10.1038/s41746-022-00640-7.
|
| [24] |
C. Dapogny, C. Dobrzynski and P. Frey, Three-dimensional adaptive domain remeshing, implicit domain meshing, and applications to free and moving boundary problems, Journal of Computational Physics, 262 (2014), 358-378.
doi: 10.1016/j.jcp.2014.01.005.
|
| [25] |
M. Davey, C. Puelz, S. Rossi, M. A. Smith, D. R. Wells, G. M. Sturgeon, W. P. Segars, J. P. Vavalle, C. S. Peskin and B. E. Griffith, Simulating cardiac fluid dynamics in the human heart, PNAS Nexus, 3 (2024), pgae392.
doi: 10.1093/pnasnexus/pgae392.
|
| [26] |
W. G. Dettmer, C. Kadapa and D. Perić, A stabilised immersed boundary method on hierarchical b-spline grids, Computer Methods in Applied Mechanics and Engineering, 311 (2016), 415-437.
doi: 10.1016/j.cma.2016.08.027.
|
| [27] |
M. M. Dicorato, P. Basile, G. Muscogiuri, M. C. Carella, M. L. Naccarati, I. Dentamaro, M. Guglielmo, A. Baggiano, S. Mushtaq, L. Fusini et al., Novel insights into non-invasive diagnostic techniques for cardiac amyloidosis: A critical review, Diagnostics, 14 (2024), 2249.
doi: 10.3390/diagnostics14192249.
|
| [28] |
S. C. Divi, P. H. van Zuijlen, T. Hoang, F. de Prenter, F. Auricchio, A. Reali, E. H. van Brummelen and C. V. Verhoosel, Residual-based error estimation and adaptivity for stabilized immersed isogeometric analysis using truncated hierarchical b-splines, Journal of Mechanics, 38 (2022), 204-237.
doi: 10.1093/jom/ufac015.
|
| [29] |
P. Du, D. An, C. Wang and J.-X. Wang, AI-powered automated model construction for patient-specific CFD simulations of aortic flows, 2025, https://arXiv.org/abs/2503.12515.
|
| [30] |
A. Düster, J. Parvizian, Z. Yang and E. Rank, The finite cell method for three-dimensional problems of solid mechanics, Computer Methods in Applied Mechanics and Engineering, 197 (2008), 3768-3782.
doi: 10.1016/j.cma.2008.02.036.
|
| [31] |
A. Embar, J. Dolbow and I. Harari, Imposing Dirichlet boundary conditions with Nitsche's method and spline-based finite elements, International Journal for Numerical Methods in Engineering, 83 (2010), 877-898.
doi: 10.1002/nme.2863.
|
| [32] |
M. Esmaily Moghadam, Y. Bazilevs, T.-Y. Hsia, I. E. Vignon-Clementel and A. L. Marsden, A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations, Computational Mechanics, 48 (2011), 277-291.
doi: 10.1007/s00466-011-0599-0.
|
| [33] |
J. E. Fromm, N. Wunsch, K. Maute, J. A. Evans and J.-S. Chen, Interpolation-based immersogeometric analysis methods for multi-material and multi-physics problems, Computational Mechanics, 75 (2025), 301-325.
doi: 10.1007/s00466-024-02506-z.
|
| [34] |
C. Geuzaine and J.-F. Remacle, Gmsh: A 3D finite element mesh generator with built-in pre-and post-processing facilities, International Journal for Numerical Methods in Engineering, 79 (2009), 1309-1331.
doi: 10.1002/nme.2579.
|
| [35] |
I. Harari and E. Grosu, A unified approach for embedded boundary conditions for fourth-order elliptic problems, International Journal for Numerical Methods in Engineering, 104 (2015), 655-675.
doi: 10.1002/nme.4813.
|
| [36] |
T. Hoang, C. V. Verhoosel, C.-Z. Qin, F. Auricchio, A. Reali and E. H. van Brummelen, Skeleton-stabilized immersogeometric analysis for incompressible viscous flow problems, Computer Methods in Applied Mechanics and Engineering, 344 (2019), 421-450.
doi: 10.1016/j.cma.2018.10.015.
|
| [37] |
D. L. Hoyert, Centers for Disease Control and Prevention, National Center for Health Statistics, NCHS Data Brief, No. 88, 2012, 75 Years of Mortality in the United States, 1935–2010 https://www.cdc.gov/nchs/data/databriefs/db88.pdf.
|
| [38] |
M.-C. Hsu, Y. Bazilevs, V. M. Calo, T. E. Tezduyar and T. J. R. Hughes, Improving stability of stabilized and multiscale formulations in flow simulations at small time steps, Computer Methods in Applied Mechanics and Engineering, 199 (2010), 828-840.
doi: 10.1016/j.cma.2009.06.019.
|
| [39] |
M.-C. Hsu, C. Wang, F. Xu, A. J. Herrema and A. Krishnamurthy, Direct immersogeometric fluid flow analysis using B-rep CAD models, Computer Aided Geometric Design, 43 (2016), 143-158.
doi: 10.1016/j.cagd.2016.02.007.
|
| [40] |
W. Huang, Metric tensors for anisotropic mesh generation, Journal of Computational Physics, 204 (2005), 633-665.
doi: 10.1016/j.jcp.2004.10.024.
|
| [41] |
A. Jacobson, L. Kavan and O. Sorkine-Hornung, Robust inside-outside segmentation using generalized winding numbers, ACM Transactions on Graphics, 32 (2013), 33.
doi: 10.1145/2461912.2461916.
|
| [42] |
M. Jaiswal, A. M. Corpuz and M.-C. Hsu, Mesh-driven resampling and regularization for robust point cloud-based flow analysis directly on scanned objects, Computer Methods in Applied Mechanics and Engineering, 432 (2024), 117426.
doi: 10.1016/j.cma.2024.117426.
|
| [43] |
D. Jia and M. Esmaily, A time-consistent stabilized finite element method for fluids with applications to hemodynamics, Scientific Reports, 13 (2023), 19120.
doi: 10.1038/s41598-023-46316-4.
|
| [44] |
D. Jia and X. Zhuang, Learning-based algorithms for vessel tracking: A review, Computerized Medical Imaging and Graphics, 89 (2021), 101840.
doi: 10.1016/j.compmedimag.2020.101840.
|
| [45] |
W. Jiang, C. Annavarapu, J. Dolbow and I. Harari, A robust Nitsche's formulation for interface problems with spline-based finite elements, International Journal for Numerical Methods in Engineering, 104 (2015), 676-696.
doi: 10.1002/nme.4766.
|
| [46] |
D. Kamensky, M.-C. Hsu, D. Schillinger, J. A. Evans, A. Aggarwal, Y. Bazilevs, M. S. Sacks and T. J. R. Hughes, An immersogeometric variational framework for fluid–structure interaction: Application to bioprosthetic heart valves, Computer Methods in Applied Mechanics and Engineering, 284 (2015), 1005-1053.
doi: 10.1016/j.cma.2014.10.040.
|
| [47] |
Y. Kato, B. Ambale-Venkatesh, Y. Kassai, L. Kasuboski, J. Schuijf, K. Kapoor, S. Caruthers and J. A. Lima, Non-contrast coronary magnetic resonance angiography: Current frontiers and future horizons, Magnetic Resonance Materials in Physics, Biology and Medicine, 33 (2020), 591-612.
doi: 10.1007/s10334-020-00834-8.
|
| [48] |
H. J. Kim, I. E. Vignon-Clementel, J. S. Coogan, C. A. Figueroa, K. E. Jansen and C. A. Taylor, Patient-specific modeling of blood flow and pressure in human coronary arteries, Annals of Biomedical Engineering, 38 (2010), 3195-3209.
doi: 10.1007/s10439-010-0083-6.
|
| [49] |
I. S. Lan, W. Yang, J. A. Feinstein, J. Kreutzer, R. T. Collins, M. Ma, G. T. Adamson and A. L. Marsden, Virtual transcatheter interventions for peripheral pulmonary artery stenosis in Williams and Alagille syndromes, Journal of the American Heart Association, 11 (2022), e023532.
doi: 10.1161/JAHA.121.023532.
|
| [50] |
A. S. Les, S. C. Shadden, C. A. Figueroa, J. M. Park, M. M. Tedesco, R. J. Herfkens, R. L. Dalman and C. A. Taylor, Quantification of hemodynamics in abdominal aortic aneurysms during rest and exercise using magnetic resonance imaging and computational fluid dynamics, Annals of Biomedical Engineering, 38 (2010), 1288-1313.
doi: 10.1007/s10439-010-9949-x.
|
| [51] |
J. Liu and A. L. Marsden, A unified continuum and variational multiscale formulation for fluids, solids, and fluid–structure interaction, Computer Methods in Applied Mechanics and Engineering, 337 (2018), 549-597.
doi: 10.1016/j.cma.2018.03.045.
|
| [52] |
G. Maher, D. Parker, N. Wilson and A. Marsden, Neural network vessel lumen regression for automated lumen cross-section segmentation in cardiovascular image-based modeling, Cardiovascular Engineering and Technology, 11 (2020), 621-635.
doi: 10.1007/s13239-020-00497-5.
|
| [53] |
G. Maher, N. Wilson and A. Marsden, Accelerating cardiovascular model building with convolutional neural networks, Medical and Biological Engineering and Computing, 57 (2019), 2319-2335.
doi: 10.1007/s11517-019-02029-3.
|
| [54] |
A. L. Marsden, I. E. Vignon-Clementel, F. Chan, J. A. Feinstein and C. A. Taylor, Effects of exercise and respiration on hemodynamic efficiency in CFD simulations of the total cavopulmonary connection, Annals of Biomedical Engineering, 35 (2007), 250-263.
doi: 10.1007/s10439-006-9224-3.
|
| [55] |
S. S. Martin, A. W. Aday, Z. I. Almarzooq, C. A. Anderson, P. Arora, C. L. Avery, C. M. Baker-Smith, B. B. Gibbs, A. Z. Beaton, A. K. Boehme et al., 2024 heart disease and stroke statistics: A report of US and global data from the American Heart Association, Circulation, 149 (2024), e347–e913.
doi: 10.1161/CIR.0000000000001247.
|
| [56] |
H. Meng, Z. Wang, Y. Hoi, L. Gao, E. Metaxa, D. D. Swartz and J. Kolega, Complex hemodynamics at the apex of an arterial bifurcation induces vascular remodeling resembling cerebral aneurysm initiation, Stroke, 38 (2007), 1924-1931.
doi: 10.1161/STROKEAHA.106.481234.
|
| [57] |
M. S. Olufsen, C. S. Peskin, W. Y. Kim, E. M. Pedersen, A. Nadim and J. Larsen, Numerical simulation and experimental validation of blood flow in arteries with structured-tree outflow conditions, Annals of Biomedical Engineering, 28 (2000), 1281-1299.
doi: 10.1114/1.1326031.
|
| [58] |
M. Oshima, R. Torii, T. Kobayashi, N. Taniguchi and K. Takagi, Finite element simulation of blood flow in the cerebral artery, Computer Methods in Applied Mechanics and Engineering, 191 (2001), 661-671.
doi: 10.1016/S0045-7825(01)00307-3.
|
| [59] |
K. Psiuk-Maksymowicz, D. Borys, B. Melka, M. Gracka, W. P. Adamczyk, M. Rojczyk, J. Wasilewski, J. Głowacki, M. Kruk, M. Nowak et al., Methodology of generation of CFD meshes and 4D shape reconstruction of coronary arteries from patient-specific dynamic CT, Scientific Reports, 14 (2024), 2201.
doi: 10.1038/s41598-024-52398-5.
|
| [60] |
A. B. Ramachandra, A. M. Kahn and A. L. Marsden, Patient-specific simulations reveal significant differences in mechanical stimuli in venous and arterial coronary grafts, Journal of Cardiovascular Translational Research, 9 (2016), 279-290.
doi: 10.1007/s12265-016-9706-0.
|
| [61] |
M. Restrepo Tique, O. Araque and L. Sanchez-Echeverri, Technological advances in the diagnosis of cardiovascular disease: A public health strategy, International Journal of Environmental Research and Public Health, 21 (2024), 1083.
doi: 10.3390/ijerph21081083.
|
| [62] |
T. Rüberg and F. Cirak, Subdivision-stabilised immersed b-spline finite elements for moving boundary flows, Computer Methods in Applied Mechanics and Engineering, 209-212 (2012), 266-283.
doi: 10.1016/j.cma.2011.10.007.
|
| [63] |
O. Sahni, J. Müller, K. E. Jansen, M. S. Shephard and C. A. Taylor, Efficient anisotropic adaptive discretization of the cardiovascular system, Computer Methods in Applied Mechanics and Engineering, 195 (2006), 5634-5655.
doi: 10.1016/j.cma.2005.10.018.
|
| [64] |
S. Sankaran, M. E. Moghadam, A. M. Kahn, E. E. Tseng, J. M. Guccione and A. L. Marsden, Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery, Annals of Biomedical Engineering, 40 (2012), 2228-2242.
doi: 10.1007/s10439-012-0579-3.
|
| [65] |
T. Sasaki, K. Takizawa and T. E. Tezduyar, Medical-image-based aorta modeling with zero-stress-state estimation, Computational Mechanics, 64 (2019), 249-271.
doi: 10.1007/s00466-019-01669-4.
|
| [66] |
K. Saurabh, B. Gao, M. Fernando, S. Xu, M. A. Khanwale, B. Khara, M.-C. Hsu, A. Krishnamurthy, H. Sundar and B. Ganapathysubramanian, Industrial scale Large Eddy Simulations with adaptive octree meshes using immersogeometric analysis, Computers and Mathematics with Applications, 97 (2021), 28-44.
|
| [67] |
D. Schillinger and M. Ruess, The Finite Cell Method: A review in the context of higher-order structural analysis of CAD and image-based geometric models, Archives of Computational Methods in Engineering, 22 (2015), 391-455.
doi: 10.1007/s11831-014-9115-y.
|
| [68] |
B. Schott, U. Rasthofer, V. Gravemeier and W. Wall, A face-oriented stabilized Nitsche-type extended variational multiscale method for incompressible two-phase flow, International Journal for Numerical Methods in Engineering, 104 (2015), 721–748.
doi: 10.1002/nme.4789.
|
| [69] |
B. Schott, S. Shahmiri, R. Kruse and W. Wall, A stabilized Nitsche-type extended embedding mesh approach for 3D low- and high-Reynolds-number flows, International Journal for Numerical Methods in Fluids, 82 (2016), 289-315.
doi: 10.1002/fld.4218.
|
| [70] |
K. Sel, D. Osman, F. Zare, S. Masoumi Shahrbabak, L. Brattain, J.-O. Hahn, O. T. Inan, R. Mukkamala, J. Palmer, D. Paydarfar et al., Building digital twins for cardiovascular health: From principles to clinical impact, Journal of the American Heart Association, 13 (2024), e031981.
doi: 10.1161/JAHA.123.031981.
|
| [71] |
J. Seo, A. B. Ramachandra, J. Boyd, A. L. Marsden and A. M. Kahn, Computational evaluation of venous graft geometries in coronary artery bypass surgery, Seminars in Thoracic and Cardiovascular Surgery, 34 (2022), 521-532.
doi: 10.1053/j.semtcvs.2021.03.007.
|
| [72] |
M. Shojima, M. Oshima, K. Takagi, R. Torii, M. Hayakawa, K. Katada, A. Morita and T. Kirino, Magnitude and role of wall shear stress on cerebral aneurysm, Stroke, 35 (2004), 2500-2505.
doi: 10.1161/01.STR.0000144648.89172.0f.
|
| [73] |
A. A. Soares, S. Gonzaga, C. Oliveira, A. Simões and A. I. Rouboa, Computational fluid dynamics in abdominal aorta bifurcation: Non-Newtonian versus Newtonian blood flow in a real case study, Computer Methods in Biomechanics and Biomedical Engineering, 20 (2017), 822-831.
doi: 10.1080/10255842.2017.1302433.
|
| [74] |
T. Sorgente, S. Biasotti, G. Manzini and M. Spagnuolo, A survey of indicators for mesh quality assessment, Computer Graphics Forum, 42 (2023), 461-483.
doi: 10.1111/cgf.14779.
|
| [75] |
N. Sveinsson Cepero and S. C. Shadden, SeqSeg: Learning local segments for automatic vascular model construction, Annals of Biomedical Engineering, 53 (2025), 158-179.
doi: 10.1007/s10439-024-03611-z.
|
| [76] |
J. M. Szafron, E. E. Heng, J. Boyd, J. D. Humphrey and A. L. Marsden, Hemodynamics and wall mechanics of vascular graft failure, Arteriosclerosis, Thrombosis, and Vascular Biology, 44 (2024), 1065-1085.
doi: 10.1161/ATVBAHA.123.318239.
|
| [77] |
T. Takikawa, D. Acuna, V. Jampani and S. Fidler, Gated-SCNN: Gated shape CNNs for semantic segmentation, in Proceedings of the IEEE/CVF International Conference on Computer Vision, (2019), 5229-5238.
|
| [78] |
K. Takizawa, Y. Bazilevs and T. E. Tezduyar, Space-time and ALE-VMS techniques for patient-specific cardiovascular fluid–structure interaction modeling, Archives of Computational Methods in Engineering, 19 (2012), 171-225.
doi: 10.1007/s11831-012-9071-3.
|
| [79] |
K. Takizawa, Y. Bazilevs, T. E. Tezduyar, C. C. Long, A. L. Marsden and K. Schjodt, ST and ALE-VMS methods for patient-specific cardiovascular fluid mechanics modeling, Mathematical Models and Methods in Applied Sciences, 24 (2014), 2437-2486.
doi: 10.1142/S0218202514500250.
|
| [80] |
K. Takizawa, C. Moorman, S. Wright, J. Christopher and T. E. Tezduyar, Wall shear stress calculations in space-time finite element computation of arterial fluid–structure interactions, Computational Mechanics, 46 (2010), 31-41.
doi: 10.1007/s00466-009-0425-0.
|
| [81] |
K. Takizawa, Y. Otoguro and T. E. Tezduyar, Variational multiscale method stabilization parameter calculated from the strain-rate tensor, Mathematical Models and Methods in Applied Sciences, 33 (2023), 1661-1691.
doi: 10.1142/S0218202523500380.
|
| [82] |
K. Takizawa, T. E. Tezduyar and Y. Otoguro, Stabilization and discontinuity-capturing parameters for space-time flow computations with finite element and isogeometric discretizations, Computational Mechanics, 62 (2018), 1169-1186.
doi: 10.1007/s00466-018-1557-x.
|
| [83] |
K. Takizawa, T. E. Tezduyar, H. Uchikawa, T. Terahara, T. Sasaki, K. Shiozaki, A. Yoshida, K. Komiya and G. Inoue, Aorta flow analysis and heart valve flow and structure analysis, in Frontiers in Computational fluid–structure Interaction and Flow Simulation (ed. T. E. Tezduyar), Springer International Publishing, Cham, (2018), 29-89.
doi: 10.1007/978-3-319-96469-0_2.
|
| [84] |
K. Takizawa, T. E. Tezduyar, H. Uchikawa, T. Terahara, T. Sasaki and A. Yoshida, Mesh refinement influence and cardiac-cycle flow periodicity in aorta flow analysis with isogeometric discretization, Computers and Fluids, 179 (2019), 790-798.
doi: 10.1016/j.compfluid.2018.05.025.
|
| [85] |
C. A. Taylor, K. Petersen, N. Xiao, M. Sinclair, Y. Bai, S. R. Lynch, A. UpdePac and M. Schaap, Patient-specific modeling of blood flow in the coronary arteries, Computer Methods in Applied Mechanics and Engineering, 417 (2023), 116414.
doi: 10.1016/j.cma.2023.116414.
|
| [86] |
T. Terahara, T. Kuraishi, K. Takizawa and T. E. Tezduyar, Computational flow analysis with boundary layer and contact representation: II. Heart valve flow with leaflet contact, Journal of Mechanics, 38 (2022), 185-194.
doi: 10.1093/jom/ufac013.
|
| [87] |
T. E. Tezduyar and Y. Osawa, Finite element stabilization parameters computed from element matrices and vectors, Computer Methods in Applied Mechanics and Engineering, 190 (2000), 411-430.
doi: 10.1016/S0045-7825(00)00211-5.
|
| [88] |
G. Trimarchi, F. Pizzino, U. Paradossi, I. A. Gueli, M. Palazzini, P. Gentile, F. Di Spigno, E. Ammirati, A. Garascia, A. Tedeschi et al., Charting the unseen: How non-invasive imaging could redefine cardiovascular prevention, Journal of Cardiovascular Development and Disease, 11 (2024), 245.
doi: 10.3390/jcdd11080245.
|
| [89] |
A. Updegrove, N. M. Wilson, J. Merkow, H. Lan, A. L. Marsden and S. C. Shadden, Simvascular: An open source pipeline for cardiovascular simulation, Annals of Biomedical Engineering, 45 (2017), 525-541.
doi: 10.1007/s10439-016-1762-8.
|
| [90] |
B. Urick, T. M. Sanders, S. S. Hossain, Y. J. Zhang and T. J. R. Hughes, Review of patient-specific vascular modeling: Template-based isogeometric framework and the case for CAD, Archives of Computational Methods in Engineering, 26 (2019), 381-404.
doi: 10.1007/s11831-017-9246-z.
|
| [91] |
T. M. van Bakel, C. J. Arthurs, F. J. Nauta, K. A. Eagle, J. A. van Herwaarden, F. L. Moll, S. Trimarchi, H. J. Patel and C. A. Figueroa, Cardiac remodelling following thoracic endovascular aortic repair for descending aortic aneurysms, European Journal of Cardio-Thoracic Surgery, 55 (2019), 1061-1070.
doi: 10.1093/ejcts/ezy399.
|
| [92] |
V. Varduhn, M.-C. Hsu, M. Ruess and D. Schillinger, The tetrahedral finite cell method: Higher-order immersogeometric analysis on adaptive non-boundary-fitted meshes, International Journal for Numerical Methods in Engineering, 107 (2016), 1054-1079.
doi: 10.1002/nme.5207.
|
| [93] |
D. M. Vigneault, W. Xie, C. Y. Ho, D. A. Bluemke and J. A. Noble, Ω-Net (Omega-Net): Fully automatic, multi-view cardiac MR detection, orientation, and segmentation with deep neural networks, Medical Image Analysis, 48 (2018), 95-106.
doi: 10.1016/j.media.2018.05.008.
|
| [94] |
C. Wang, F. Xu, M.-C. Hsu and A. Krishnamurthy, Rapid B-rep model preprocessing for immersogeometric analysis using analytic surfaces, Computer Aided Geometric Design, 52/53 (2017), 190-204.
doi: 10.1016/j.cagd.2017.03.002.
|
| [95] |
X. Wang, M. Jaiswal, A. M. Corpuz, S. Paudel, A. Balu, A. Krishnamurthy, J. Yan and M.-C. Hsu, Photogrammetry-based computational fluid dynamics, Computer Methods in Applied Mechanics and Engineering, 417 (2023), 116311.
doi: 10.1016/j.cma.2023.116311.
|
| [96] |
J. C. Weddell, J. Kwack, P. Imoukhuede and A. Masud, Hemodynamic analysis in an idealized artery tree: Differences in wall shear stress between Newtonian and non-Newtonian blood models, PLoS ONE, 10 (2015), e0124575.
doi: 10.1371/journal.pone.0124575.
|
| [97] |
N. M. Wilson, A. K. Ortiz and A. B. Johnson, The vascular model repository: A public resource of medical imaging data and blood flow simulation results, Journal of medical devices, 7 (2013), 040923.
doi: 10.1115/1.4025983.
|
| [98] |
J. M. Wolterink, T. Leiner, M. A. Viergever and I. Išgum, Automatic segmentation and disease classification using cardiac cine MR images, in Statistical Atlases and Computational Models of the Heart. ACDC and MMWHS Challenges: 8th International Workshop, STACOM 2017, Held in Conjunction with MICCAI 2017, Quebec City, Canada, September 10-14, 2017, Revised Selected Papers 8, Springer, (2018), 101-110.
doi: 10.1007/978-3-030-12029-0_23.
|
| [99] |
J. M. Wolterink, T. Leiner, M. A. Viergever and I. Išgum, Dilated convolutional neural networks for cardiovascular MR segmentation in congenital heart disease, in International Workshop on Reconstruction and Analysis of Moving Body Organs, Springer, (2016), 95-102.
doi: 10.1007/978-3-319-52280-7_9.
|
| [100] |
Q. Xia, Y. Yao, Z. Hu and A. Hao, Automatic 3D atrial segmentation from GE-MRIs using volumetric fully convolutional networks, in Statistical Atlases and Computational Models of the Heart. Atrial Segmentation and LV Quantification Challenges: 9th International Workshop, STACOM 2018, Held in Conjunction with MICCAI 2018, Granada, Spain, September 16, 2018, Revised Selected Papers 9, Springer, (2019), 211-220.
doi: 10.1007/978-3-030-12029-0_23.
|
| [101] |
F. Xu, D. Schillinger, D. Kamensky, V. Varduhn, C. Wang and M.-C. Hsu, The tetrahedral finite cell method for fluids: Immersogeometric analysis of turbulent flow around complex geometries, Computers and Fluids, 141 (2016), 135-154.
doi: 10.1016/j.compfluid.2015.08.027.
|
| [102] |
S. Xu, F. Xu, A. Kommajosula, M.-C. Hsu and B. Ganapathysubramanian, Immersogeometric analysis of moving objects in incompressible flows, Computers and Fluids, 189 (2019), 24-33.
doi: 10.1016/j.compfluid.2019.05.018.
|
| [103] |
Y. Zhang, Y. Bazilevs, S. Goswami, C. L. Bajaj and T. J. R. Hughes, Patient-specific vascular NURBS modeling for isogeometric analysis of blood flow, Computer Methods in Applied Mechanics and Engineering, 196 (2007), 2943-2959.
doi: 10.1016/j.cma.2007.02.009.
|
| [104] |
G. Zieff, L. C. Bates-Fraser, P. P. Lassalle, C. Paterson, K. Heffernan and L. Stoner, Editorial: Non-invasive physiological measurements: From discovery to implementation, Frontiers in Physiology, 14 (2023), 1356664.
doi: 10.3389/fphys.2023.1356664.
|
| [105] |
Centers for Disease Control and Prevention, National Center for Health Statistics, Leading Causes of Death, 1900–1998, 2001, https://stacks.cdc.gov/view/cdc/53236.
|
| [106] |
HeartFlow, Inc., Heartflow: Revolutionary non-invasive coronary artery disease test, 2024, https://www.heartflow.com/.
|
| [107] |
National Academies of Sciences, Engineering, and Medicine, Foundational Research Gaps and Future Directions for Digital Twins, The National Academies Press, Washington, DC, 2024.
doi: 10.17226/26894.
|
The pipeline of the proposed framework from initial CT scan to simulation results
Diagram of the schematic representation of the proposed pipeline for aorta segmentation and 3D mesh reconstruction. The image segmentation architecture utilizes 3D U-Net for segmenting the aorta and its branches. The encoder component of 3D U-Net is responsible for extracting high-dimensional features from the small field of view. These features are then combined with the image gradient obtained from the MRI image within the shape stream
Illustration of (A) parameters used to define the anisotropic mesh, and a comparison showcasing the difference between (B) the isotropic mesh and (C) the anisotropic mesh
(A) An example of flow within an object, illustrated using a 2D cross-section of a 3D channel represented by a point cloud. The object with boundary
Background mesh domain containing an aorta model represented by a point cloud. (A) A domain whose boundaries do not coincide with the flow inlet and outlet locations. (B) A clipped domain in which both the domain and the point cloud are trimmed to align with the desired inlet and outlet faces
Dimensions of the simulated pipe using 500,000 points. The overall domain is a cuboid of size
Visualization of the pipe flow IM3 mesh with zoomed-in figures of the side and front of the mesh
(A) Velocity magnitude and (B) pressure at the central cross-section along the streamwise direction of the pipe for
(A) Pressure drop along the centerline of the pipe in the streamwise direction and (B) streamwise velocity profile across the outlet of the pipe for all three background meshes for
Pressure drop along the centerline of the pipe in the streamwise direction with
The left shows the point cloud, dimensions, flow properties, and boundary conditions for the bifurcation flow example at
(A) Velocity magnitude and pressure contours along the center plane. (B) Pressure drop versus distance along the centerline of the pipe before the branches split, with the maximum distance being the bifurcation apex at 13.95 mm. The present work is compared with those from previous studies by Chung [21] and Weddell et al. [96]
Resampled point clouds colored by normalized minimum and maximum depth values for enhanced visualization. From left to right: Case 1 – Adult patient with single ventricle disease who had undergone the Fontan procedure. Case 2 – Pediatric patient with a diseased ventricle and healthy aorta. Case 3 – Healthy lamb aorta
From left to right: Voxel representation, NIMBUS point cloud representation, and both representations overlaid with semi-translucent voxels for Case 2
(A) The domain of the CT scan accompanied by a slice across the midsection to reveal the aorta. (B) Point output from performing auto-segmentation on the CT scan
Dimensions of the simulation domain boundary with the point cloud aorta immersed inside. The smaller boxes visualized are Boolean subtractions of the original rectangular domain using an estimated centerline normal to align the background mesh surface with the desired inlet/outlet surfaces. No-slip condition on the arterial wall is enforced weekly using the point cloud
From left to right: Volume rendering of the blood flow speed. Velocity streamlines representation with point cloud overlaid and a zoomed-in portion of disturbed flow in the aorta
(A) Volume rendering of pressure for the full aorta and zoomed-in sections of branching arteries. (B) Velocity magnitude visualization at user-selected slices throughout the simulation domain
Illustration of feedback loops in a cardiovascular digital twin and the role of high-fidelity simulation