1
|
McCarthy RP, Mason PJ, Marks DS, LaDisa JF. Influence of boundary conditions and blood rheology on indices of wall shear stress from IVUS-based patient-specific stented coronary artery simulations. Sci Rep 2025; 15:15868. [PMID: 40335561 PMCID: PMC12059177 DOI: 10.1038/s41598-025-99066-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2024] [Accepted: 04/15/2025] [Indexed: 05/09/2025] Open
Abstract
The long-term clinical efficacy of coronary stents is limited by restenosis. Coronary stenting results in altered arterial geometry, local blood flow patterns, and wall shear stress (WSS), all of which can influence restenosis. Computational fluid dynamics (CFD) simulations employ assumptions about blood properties and boundary conditions, which also influence WSS alterations from stenting. Our objective was to evaluate three common assumptions applied with stented coronary artery CFD simulations (inlet velocity profile, outlet boundary conditions, and viscosity) to provide insight for future studies. A patient-specific right coronary artery was reconstructed from intravascular ultrasound and computed tomography imaging. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) were compared for CFD simulations using parabolic and Womersley (α = 2.5) velocity profiles and Newtonian versus non-Newtonian (Carreau-Yasuda) viscosity. TAWSS and OSI differences were quantified for 5-element lumped parameter network (LPN) outlet boundary conditions as compared to a 3-element Windkessel approach previously applied by neglecting ventricular contraction. Differences in velocity profiles were negligible beyond two diameters from the inlet. Differences between Womersley inlet 3-element and 5-element LPNs were ~ 5% for TAWSS and ~ 200% for OSI, and most pronounced near stent struts. OSI differences were due to pressure differences between outlet boundary conditions inducing different near-wall velocity gradients. TAWSS differences between viscosity models were greatest near struts (~ 100%). Collectively these results suggest certain assumptions commonly applied for simulations of stented coronary arteries have a greater impact on TAWSS and OSI than others, with outlet boundary conditions being paramount, followed by viscosity model and inlet velocity profile.
Collapse
Affiliation(s)
- R Patrick McCarthy
- Department of Biomedical Engineering, Medical College of Wisconsin and Marquette University, Milwaukee, WI, 53226, USA.
| | - Peter J Mason
- Department of Medicine; Division of Cardiovascular Medicine, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
- Department of Medicine; Cardiovascular Medicine, University of Wisconsin, Madison, WI, 53792, USA
| | - David S Marks
- Department of Medicine; Division of Cardiovascular Medicine, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
| | - John F LaDisa
- Department of Biomedical Engineering, Medical College of Wisconsin and Marquette University, Milwaukee, WI, 53226, USA
- Department of Medicine; Division of Cardiovascular Medicine, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
- Department of Pediatrics; Section of Cardiology, Herma Heart Institute, Children's Wisconsin and the Medical College of Wisconsin, Milwaukee, WI, 53226, USA
- Department of Physiology, Medical College of Wisconsin, Milwaukee, WI, 53226, USA
| |
Collapse
|
2
|
Blanch-Granada A, LaDisa JF, Samyn MM, Cava JR, Handler SS, Gerardin JF, Goot B, Maadooliat M, Hraška V. Hemodynamic Evaluation of Norwood Aortic Arch Geometry Compared to Native Arch Controls. J Biomech Eng 2025; 147:051003. [PMID: 39992380 DOI: 10.1115/1.4068009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2024] [Accepted: 02/19/2025] [Indexed: 02/25/2025]
Abstract
The Norwood procedure creates a reconstructed neo-aorta to provide unobstructed systemic cardiac output for hypoplastic left heart syndrome patients. We used patient-specific computational fluid dynamics (CFD) simulations incorporating physiologic boundary conditions to quantify hemodynamics for reconstructed aortic arch geometries versus native aortic arches from a control group of single ventricle patients. We hypothesized that reconstructed arches from Norwood patients (n = 5) would experience significant differences in time-averaged wall shear stress normalized to body surface area (TAWSSnBSA), oscillatory shear index (OSI), energy efficiency (Eeff), and energy loss (EL) versus controls (n = 3). CFD simulations were conducted using 3 T cardiac magnetic resonance imaging, blood flow, and pressure data. Simulations incorporated downstream vascular resistance and compliance to replicate patient physiology. TAWSSnBSA and OSI were quantified axially and circumferentially. Global differences in Eeff and EL were compared. Significance was assessed by Mann-Whitney U test. Norwood patients had higher TAWSSnBSA distal to the transverse arch at locations of residual narrowing presenting following coarctation correction, as well as higher OSI within ascending aorta and transverse arch regions (p < 0.05). EL correlated with patient features including cardiac output (r = 0.9) and BT-shunt resistance (r = -0.63) but did not correlate with arch measurements or morphology. These results indicate reconstructed arches from Norwood patients are exposed to altered wall shear stress and energy indices linked to cellular proliferation and inefficiency in prior studies. These results may help clinicians further understand what constitutes an optimally reconstructed arch after confirmation in larger studies.
Collapse
Affiliation(s)
- Aloma Blanch-Granada
- Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI 53233
| | - John F LaDisa
- Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI 53233; Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226; Department of Physiology, Medical College of Wisconsin, Milwaukee, WI 53226
| | - Margaret M Samyn
- Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI 53233; Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226
| | - Joseph R Cava
- Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI 53233; Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226
| | - Stephanie S Handler
- Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226
- Medical College of Wisconsin
| | - Jennifer F Gerardin
- Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226
- Medical College of Wisconsin
| | - Benjamin Goot
- Section of Pediatric Cardiology, Department of Pediatrics, Medical College of Wisconsin, Milwaukee, WI 53226; Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226
- Medical College of Wisconsin
| | - Mehdi Maadooliat
- Department of Mathematical and Statistical Sciences, Marquette University, Milwaukee, WI 53233
- Marquette University
| | - Viktor Hraška
- Herma Heart Institute, Children's Wisconsin, Milwaukee, WI 53226; Department of Surgery - Division of Pediatric Cardiothoracic Surgery, Children's Wisconsin and the Medical College of Wisconsin, Milwaukee, WI 53226
- Medical College of Wisconsin
| |
Collapse
|
3
|
Puelz C, Rusin CG, Lior D, Sachdeva S, Doan TT, Eilers LF, Reaves-O'Neal D, Akula M, Molossi S. Fluid-structure interaction simulations for the prediction of fractional flow reserve in pediatric patients with anomalous aortic origin of a coronary artery. Med Eng Phys 2025; 136:104293. [PMID: 39979010 DOI: 10.1016/j.medengphy.2025.104293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2024] [Revised: 01/02/2025] [Accepted: 01/08/2025] [Indexed: 02/22/2025]
Abstract
Computer simulations of blood flow in patients with anomalous aortic origin of a coronary artery (AAOCA) have the promise to provide insight into this complex disease. They provide an in silico experimental platform to explore possible mechanisms of myocardial ischemia, a potentially deadly complication for patients with this defect. This paper focuses on the question of model calibration for fluid-structure interaction models of pediatric AAOCA patients. Imaging and cardiac catheterization data provide partial information for model construction and calibration. However, parameters for downstream boundary conditions needed for these models are difficult to estimate. Further, important model predictions, like fractional flow reserve (FFR), are sensitive to these parameters. We describe an approach to calibrate downstream boundary condition parameters to clinical measurements of resting FFR. The calibrated models are then used to predict FFR at stress, an invasively measured quantity that can be used in the clinical evaluation of these patients. We find reasonable agreement between the model predicted and clinically measured FFR at stress, indicating the credibility of this modeling framework for predicting hemodynamics of pediatric AAOCA patients. This approach could lead to important clinical applications since it may serve as a tool for risk stratifying children with AAOCA.
Collapse
Affiliation(s)
- Charles Puelz
- Department of Mathematics, University of Houston, United States of America; Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America.
| | - Craig G Rusin
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Dan Lior
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Shagun Sachdeva
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Tam T Doan
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Lindsay F Eilers
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Dana Reaves-O'Neal
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Monisha Akula
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| | - Silvana Molossi
- Department of Pediatrics, Division of Cardiology, Baylor College of Medicine and Texas Children's Hospital, United States of America
| |
Collapse
|
4
|
Doan TT, Puelz C, Rusin C, Molossi S. Anomalous Aortic Origin of a Coronary Artery in Pediatric Patients. CURRENT PEDIATRICS REPORTS 2024; 12:69-80. [PMID: 39816672 PMCID: PMC11729077 DOI: 10.1007/s40124-024-00317-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 05/15/2024] [Indexed: 01/18/2025]
Abstract
Purpose of Review We present a contemporary approach to risk assessment and management of patients with anomalous aortic origin of a coronary artery (AAOCA). Recent Findings Anomalous left coronary artery from the right aortic sinus (L-AAOCA) with interarterial course carries a high risk of sudden cardiac death (SCD); therefore, current guidelines recommend exercise restriction and surgical intervention. Recent data in intraseptal and juxtacommissural L-AAOCA showed inducible perfusion abnormalities, leading to consideration of surgical intervention. Anomalous right coronary artery from the left aortic sinus (R-AAOCA) carries a much lower risk and stress perfusion imaging is helpful in identifying patients with inducible ischemia. Perfusion abnormalities resolve following successful surgical intervention of AAOCA. Computational modeling techniques identifying risk features shows promise in the evaluation of AAOCA. Summary Stress perfusion imaging is helpful in assessing AAOCA upon presentation and following surgical intervention. Computational modeling has potential in bridging knowledge gaps in AAOCA. Supplementary Information The online version contains supplementary material available at 10.1007/s40124-024-00317-7.
Collapse
Affiliation(s)
- Tam T. Doan
- Coronary Artery Anomalies Program, Division of Cardiology, Texas Children’s Hospital, 6651 Main Street MC-E1920, Houston, TX 77030 USA
- Department of Pediatrics, Baylor College of Medicine, Houston, TX 77030 USA
| | - Charles Puelz
- Department of Pediatrics, Baylor College of Medicine, Houston, TX 77030 USA
| | - Craig Rusin
- Department of Pediatrics, Baylor College of Medicine, Houston, TX 77030 USA
| | - Silvana Molossi
- Coronary Artery Anomalies Program, Division of Cardiology, Texas Children’s Hospital, 6651 Main Street MC-E1920, Houston, TX 77030 USA
- Department of Pediatrics, Baylor College of Medicine, Houston, TX 77030 USA
| |
Collapse
|
5
|
Milovanovic A, Saveljic I, Filipovic N. Numerical vs analytical comparison with experimental fractional flow reserve values of right coronary artery stenosis. Technol Health Care 2022; 31:977-990. [PMID: 36442165 DOI: 10.3233/thc-220435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
BACKGROUND: The fractional flow reserve (FFR) index has been widely accepted as a standard diagnostic method for identifying functional relevance of coronary stenosis. Since the invasive techniques used for its determination are associated with a certain risk of vascular injury, as well as with an increased cost, several non-invasive procedures have been developed. OBJECTIVE: The aim of this study was to compare FFR values for the coronary artery obtained by computational fluid dynamics (CFD) and coronary computed tomography angiography (CCTA). METHODS: Computation of FFR has been performed using both numerical and the analytical method. The numerical method employs CFD to solve the governing equations which relate to mass and momentum conservation (the continuity equation and the Navier-Stokes equations) as well as CCTA to generate the three-dimensional computational domain. After imposing the appropriate boundary conditions, the values of the pressure change are calculated and the FFR index is determined. Based on Bernoulli’s law, the analytical method calculates the overall pressure drop across the stenosis in the coronary artery, enabling FFR determination. RESULTS: The clinical data for twenty patients who underwent invasive coronary angiography are used to validate the results obtained by using CFD (together with CCTA) simulation and analytical solution. The medically measured FFR compared to the analytical one differs by about 4%, while, the difference is about 2.6% when compared to the numerical FFR. For FFR values below 0.8 (which are considered to be associated with myocardial ischemia) the standard error has a value of 0.01201, while the standard deviation is 0.02081. For FFR values above 0.80, these values are slightly higher. Bland-Altman analysis showed that medical measurement and numerical FFR were in good agreement (SD = 0.0292, p< 0.0001). CONCLUSIONS: The analytically calculated FFR has a slightly lower coefficient of determination than the numerically computed FFR when compared with experimental one. However, it can still give a reliable answer to the question of whether patients need a stent, bypass surgery or only drug treatment and it requires a significantly lower computation time.
Collapse
Affiliation(s)
| | - Igor Saveljic
- Institute for Information Technologies, University of Kragujevac, Kragujevac, Serbia
- Bioengineering Research and Development Center, Kragujevac, Serbia
| | - Nenad Filipovic
- Faculty of Engineering, University of Kragujevac, Kragujevac, Serbia
- Bioengineering Research and Development Center, Kragujevac, Serbia
| |
Collapse
|
6
|
Chidyagwai SG, Vardhan M, Kaplan M, Chamberlain R, Barker P, Randles A. Characterization of hemodynamics in anomalous aortic origin of coronary arteries using patient-specific modeling. J Biomech 2022; 132:110919. [PMID: 35063831 PMCID: PMC10712838 DOI: 10.1016/j.jbiomech.2021.110919] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2021] [Revised: 12/14/2021] [Accepted: 12/15/2021] [Indexed: 12/01/2022]
Abstract
The anomalous aortic origin of coronary arteries (AAOCA) is a congenital disease that can lead to sudden cardiac death (SCD) during strenuous physical activity. Despite AAOCA being the second leading cause of SCD among young athletes, the mechanism behind sudden cardiac death remains mostly unknown. Computational fluid dynamics provides a powerful tool for studying how pathologic anatomy can affect different hemodynamic states. The present study investigates the effect of AAOCA on patient hemodynamics. We performed patient-specific hemodynamic simulations of interarterial AAOCA at baseline and in the exercise state using our massively parallel flow solver. Additionally, we investigate how surgical correction via coronary unroofing impacts patient blood flow. Results show that patient-specific AAOCA models exhibited higher interarterial time-averaged wall shear stress (TAWSS) values compared to the control patients. The oscillatory shear index had no impact on AAOCA. Finally, the coronary unroofing procedure normalized the elevated TAWSS by decreasing TAWSS in the postoperative patient. The present study provides a proof of concept for the potential hemodynamic factors underlying coronary ischemia in AAOCA during exercise state.
Collapse
Affiliation(s)
- Simbarashe G Chidyagwai
- Department of Biomedical Engineering, Duke University, Durham, NC, United States of America.
| | - Madhurima Vardhan
- Department of Biomedical Engineering, Duke University, Durham, NC, United States of America.
| | - Michael Kaplan
- Duke University School of Medicine, Duke University, Durham, NC, United States of America.
| | - Reid Chamberlain
- Department of Medicine, Duke University, Durham, NC, United States of America; Department of Pediatrics, Duke University, Durham, NC, United States of America.
| | - Piers Barker
- Department of Medicine, Duke University, Durham, NC, United States of America; Department of Pediatrics, Duke University, Durham, NC, United States of America.
| | - Amanda Randles
- Department of Biomedical Engineering, Duke University, Durham, NC, United States of America.
| |
Collapse
|
7
|
Patient-Specific Fluid-Structure Simulations of Anomalous Aortic Origin of Right Coronary Arteries. JTCVS Tech 2022; 13:144-162. [PMID: 35711199 PMCID: PMC9196314 DOI: 10.1016/j.xjtc.2022.02.022] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 02/16/2022] [Indexed: 01/10/2023] Open
Abstract
Objectives Anomalous aortic origin of the right coronary artery (AAORCA) may cause ischemia and sudden death. However, the specific anatomic indications for surgery are unclear, so dobutamine-stress instantaneous wave-free ratio (iFR) is increasingly used. Meanwhile, advances in fluid–structure interaction (FSI) modeling can simulate the pulsatile hemodynamics and tissue deformation. We sought to evaluate the feasibility of simulating the resting and dobutamine-stress iFR in AAORCA using patient-specific FSI models and to visualize the mechanism of ischemia within the intramural geometry and associated lumen narrowing. Methods We developed 6 patient-specific FSI models of AAORCA using SimVascular software. Three-dimensional geometries were segmented from coronary computed tomography angiography. Vascular outlets were coupled to lumped-parameter networks that included dynamic compression of the coronary microvasculature and were tuned to each patient's vitals and cardiac output. Results All cases were interarterial, and 5 of 6 had an intramural course. Measured iFRs ranged from 0.95 to 0.98 at rest and 0.80 to 0.95 under dobutamine stress. After we tuned the distal coronary resistances to achieve a stress flow rate triple that at rest, the simulations adequately matched the measured iFRs (r = 0.85, root-mean-square error = 0.04). The intramural lumen remained narrowed with simulated stress and resulted in lower iFRs without needing external compression from the pulmonary root. Conclusions Patient-specific FSI modeling of AAORCA is a promising, noninvasive method to assess the iFR reduction caused by intramural geometries and inform surgical intervention. However, the models’ sensitivity to distal coronary resistance suggests that quantitative stress-perfusion imaging may augment virtual and invasive iFR studies.
Collapse
|