- 1Department of Mechanical Engineering, University of Bristol, Bristol, United Kingdom
- 2University Hospitals Bristol and Weston, NHS Foundation Trust, Bristol, United Kingdom
- 3Bristol Medical School, University of Bristol, Bristol, United Kingdom
- 4National Heart and Lung Institute, Imperial College London, London, United Kingdom
Objectives: Effective management of aortic coarctation (CoA) affects long-term cardiovascular outcomes. Full appreciation of CoA hemodynamics is important. This study aimed to analyze the relationship between aortic shape and hemodynamic parameters by means of computational simulations, purposely isolating the morphological variable.
Methods: Computational simulations were run in three aortic models. MRI-derived aortic geometries were generated using a statistical shape modeling methodology. Starting from n = 108 patients, the mean aortic configuration was derived in patients without CoA (n = 37, “no-CoA”), with surgically repaired CoA (n = 58, “r-CoA”) and with unrepaired CoA (n = 13, “CoA”). As such, the aortic models represented average configurations for each scenario. Key hemodynamic parameters (i.e., pressure drop, aortic velocity, vorticity, wall shear stress WSS, and length and number of strong flow separations in the descending aorta) were measured in the three models at three time points (peak systole, end systole, end diastole).
Results: Comparing no-CoA and CoA revealed substantial differences in all hemodynamic parameters. However, simulations revealed significant increases in vorticity at the site of CoA repair, higher WSS in the descending aorta and a 12% increase in power loss, in r-CoA compared to no-CoA, despite no clinically significant narrowing (CoA index >0.8) in the r-CoA model.
Conclusions: Small alterations in aortic morphology impact on key hemodynamic indices. This may contribute to explaining phenomena such as persistent hypertension in the absence of any clinically significant narrowing. Whilst cardiovascular events in these patients may be related to hypertension, the role of arch geometry may be a contributory factor.
Introduction
Aortic coarctation (CoA) is a discrete congenital heart defect with an incidence of 1/2,500 live births that can occur in isolation or in association with other left sided congenital lesions, such as bicuspid aortic valve (BAV) (1). The association between CoA and BAV is reported in up to 85% of CoA cases and has significant clinical implications on long-term outcomes, morbidity and mortality (1). Cardiovascular effects are mostly through new or persistent systemic hypertension; this may persist after surgery, which suggests that patients retain an abnormal vascular phenotype post-repair.
Aortic arch geometry patterns have been described as either “Gothic” (more angulated) or “Romanesque” (rounder) morphology. It has been suggested that arch geometry, together with increased arterial stiffness, contribute to the development of hypertension in CoA patients (2, 3). Abnormal aortic arch geometry has been associated with rest or exercise-induced hypertension in the long-term follow-up of patients after successful repair of CoA without residual obstruction of the arch (3–5). Hypertension has been linked to lower distensibility and elevated stiffness of the ascending aorta and greater loss of systolic wave amplitude across the aortic arch. This could partially explain the presence of late hypertension in the long-term follow-up of these patients, even in the absence of residual arch stenosis (6, 7). Furthermore, gothic arch geometry can unfavorably impact on left ventricular mass index (LVMI), independent of age and resting systolic blood pressure, suggesting a chronic increase in left ventricular afterload in the presence of gothic arch geometry (4).
The concept of “optimal surgical shape,” based on a three-dimensional (3D) analysis of aortic morphology in BAV patients with and without CoA, has been recently introduced into the literature. This showed that patients with repaired CoA were more likely to have a gothic aortic arch morphology and this has been previously associated with new or persistent systemic hypertension in the long term (8). Furthermore, the “selfish brain” hypothesis was recently described, showing that vertebral artery hypoplasia with an incomplete posterior circle of Willis in repaired CoA subjects may be important in the development of hypertension or its persistence, following CoA repair (9). This is suggestive of a more global vascular developmental phenomenon, rather than simple hemodynamics.
Undoubtedly, CoA leads to long-term ventricular and arterial implications and cardiovascular outcomes, with these patients presenting significantly lower survival than the normal population (10), even in the modern era. Therefore, a thorough understanding of CoA hemodynamics and their ramifications is important. The motivation of the study is indeed to make use of computational modeling and the insight they offer in terms of isolating a variable of interest (in this case, aortic arch shape) to contribute to the understanding of aortic pressure changes in CoA patients. The broad clinical relevance of this approach is not only related to the observation that hypertension may be an inevitable consequence of CoA, even following early effective anatomical repair, but that maladaptive processes may be at play in these patients (11). Also in the context of interventional rather than surgical repair, the importance of the aortic arch anatomy in CoA patients is increasingly recognized as a parameter to define cases at higher risk of residual hypertension, even despite optimized isthmic stent implantation (12).
This work focuses purely on fluidodynamic aspects, aiming to analyze the relationship between aortic shape and hemodynamic parameters by means of computational fluid dynamics, purposely isolating the morphological variable and comparing patients with successfully repaired CoA, with unrepaired CoA and without CoA.
Methods
Aortic Arch Models
The study is based on three aortic models derived from a previously described population of patients with BAV and CoA (8). Starting from n = 108 patients, the mean aortic arch configuration was derived in patients without CoA (n = 37, “No-CoA”), with surgically repaired CoA (n = 58, “r-CoA”) and with unrepaired CoA (n = 13, “CoA”). The aortic shapes were reconstructed from magnetic resonance imaging (MRI) data and the final three configurations were obtained by means of statistical shape modeling, thus producing three average models summarizing the geometrical features of the patients in each group. This resulted in three representative models for the three scenarios of interest (Figure 1). More details on the process of image reconstruction are available in (8) and details on the statistical shape modeling framework are available in (13). All datasets were anonymized and, in view of the retrospective study design, formal ethical approval was waived by the local Institutional Research and Innovation Department, and the rest of the study focused on computational simulations.
Computational Fluid Dynamics
Computational simulations were run in the three aortic models, purposely isolating the geo-metrical variable. The mesh generation and the numerical simulation were both carried out in STAR-CCM+ 13.04.010-r8 (Siemens). The flow was simulated using an unsteady implicit scheme, with the SIMPLE algorithm to solve the pressure and velocity in a segregated manner. Second order accuracy was chosen for both temporal and spatial discretization. The time step size and the convergence criteria were set to be t = 0.001 s and ε = 10−8, respectively. The fluid was assumed to be incompressible and Newtonian, and the density and dynamic viscosity were set to ρ = 1, 060 kg/m3 and μ = 0.004 Pa.s, respectively.
A polyhedral computational grid together with 10 prism layers was generated with an average total number of 6 M elements for the finite volume solver. In order to capture complex flow features as much as possible, the flow was modeled using Large Eddy Simulation (LES), more specifically Dynamic Smagorinsky Subgrid Scale. Literature has suggested that this model is reasonably capable of simulating aortic flows (14–16). A realistic aortic flow waveform derived from MRI (Figure 2) was prescribed as the inlet velocity profile, with average flow rate of 6 L/min, peak velocity of 53.7 cm/s and average velocity of 11.7 cm/s. The outlet boundary condition was considered to be zero pressure. Ten cardiac cycles were simulated and only the last one was used for subsequent analyzes.
In order to compare the computed hemodynamic results with respect to geometric characteristics of the aortic arch, the medial axis of each geometry was computed using Vascular Modeling Toolkit (VMTK) (17). One-hundred cross-sections along the aorta were further generated such that their normal vectors were locally tangent to the centerline. The average of flow parameters of interest (i.e., pressure, velocity, wall shear stress, vorticity) were measured on each cross-section along the aorta. The slices at positions 40, 50, and 60 were located before, across, and after the coarctation, respectively, and were selected for direct comparison between the different geometries. The slice at position 90 was also included to capture the effect of coarctation on the blood flow properties downstream of the narrowing, in the descending aorta. The location of the slices was normalized from 0 to 1, and is displayed in Figure 3.
Comparisons across the three aortic configurations were carried out at three time points, i.e., peak systole (0.13 s, 16% of the cardiac cycle), end systole (0.34 s, 43% of the cardiac cycle) and end diastole (0.80 s, 100% of the cardiac cycle). Key measurements included velocity (m/s, derived directly from the simulations), pressure drop (mmHg, measured as P(S0.9)–P(S0.4), accounting for pressure recovery) and vorticity (defined as the curl of the velocity, which is a measure of the rotation of the flow, /s). Wall shear stress (WSS, dyn/cm2) was also derived from the simulations. In order to identify areas of elevated WSS and draw comparisons between the three geometries, the average of the obtained WSS was calculated and areas with values higher than this were considered as elevated WSS. The length and number of strong separations in the descending aorta were also investigated by comparing velocity streamlines along the aortic arches.
Data Analysis
Results were compared for the three aortic configurations at the above mentioned three time instances along the length of the aorta. Whilst the aortic models summarize the geometrical features of dozens of patients, the comparison was ultimately carried out on the three statistically derived shapes. This precluded a full statistical analysis, which was beyond the scope of this study, rather taking advantage of the fact that the representative aortic geometries in themselves summarize features from multiple patients.
Results
Velocity, pressure and vorticity values for all geometries and time points are reported in Table 1.
Table 1. Average values of velocity (V), pressure (P), and vorticity (ω) on slices of 40, 50, 60, and 90.
Unsurprisingly, comparison between the no-CoA and the CoA models revealed much higher aortic pressure in the ascending aorta (four-fold increase at peak systole) and a substantial pressure drop across the coarctation site (5.20 vs. 36.50 mmHg; Figure 4) in the CoA group, as well as much higher aortic velocity past the CoA (1.14 vs. 3.61 m/s). Vorticity showed approximately a four-fold increase between the two configurations in the descending aorta (387 vs. 1,495 1/s at location 0.5 and 251 vs. 1,078 1/s at location 0.6). Vortex cores and vorticity across the aortic cross-section in the descending aorta are shown in Figure 5. Given the obtained WSS maps (Figure 6), areas of elevated WSS were identified (WSS >50 dyn/cm2) and noted in the descending aorta and no-CoA and CoA configurations were compared. There was an average 378% increase in the CoA group. Finally, a three-fold increase in power loss was observed at peak systole in the ascending aorta in the CoA group compared to the no-CoA group (0.9 vs. 3.2 W, No-CoA vs. CoA; Figure 7).
Figure 5. Result of the simulations at end systole. First row: vortical structure by iso-surface of λ2 = −30, 000 in the descending aorta. Second row: contour of vorticity (/s) along with velocity streamlines. Third row: contour of pressure along with velocity streamlines. Fourth row: velocity cross-section (m/s). The cross-sections are identified in the first row.
Figure 7. Plot of power loss along the aorta. Order of plots from top to bottom: end systole, peak systole, and end diastole.
More interestingly, differences were also observed when comparing the no-CoA and r-CoA configurations, despite the latter representing qualitatively a successful CoA repair with a coarctation index >0.8 (i.e., no significant narrowing from a clinical standpoint). The hemodynamic observations in the r-CoA were mainly:
1) increased vorticity at the site of CoA repair (no-CoA: 387 1/s vs. r-CoA: 711 1/s at location 0.5), Figure 5;
2) regions of elevated WSS in the descending aorta of the r-CoA model more comparable to the CoA configuration (Figure 6), with an increase in areas of WSS (>50 dyn/cm2) in 153% of the r-CoA model;
3) increased power loss (+22%) in the ascending aorta (no-CoA: 0.9 W vs. r-CoA: 1.1 W; Figure 7).
The velocity streamlines along the three aortic models are plotted in Figure 8. Two internal and external separations in the CoA aorta can be appreciated, leading to adverse pressure gradient in the descending aorta to such an extent that the average pressure in cross-sections S0.5 and S0.6 becomes negative at peak systole, as reported in Table 1. Internal separation was only observed in no-CoA and r-CoA aorta. The length of separations was measured and is reported in Table 2.
Figure 8. Contour of velocity and velocity streamlines along the aorta at peak systole, end systole, and end diastole. Top row: No-CoA aorta; Middle row: r-CoA; Bottom Row: CoA.
Discussion
Computational simulations have been providing an accessible and versatile tool for gathering insight into the fluid dynamics of aortic coarctation for over a decade. Computational modeling analysis suggests that the concomitant presence of BAV with CoA results in increased maximal velocity, secondary flow, pressure loss, time-averaged wall shear stress and oscillatory shear index in the descending aorta, past the CoA, compared to a tricuspid valve scenario (18). Computational fluid dynamics can highlight regions of the thoracic aorta with unfavorable hemodynamics (19) and, in the presence of CoA, reveal differences in local WSS (20). Such insight can yield clinically relevant messages. For instance, WSS changes can be related to areas of plaque formation in locations influenced by surgical CoA repair, such as resection with end-to-end anastomosis (19). The modeling literature suggests that simulations can be used to assess changes in aortic wall biomechanics after a percutaneous procedure (e.g., stenting) (21), to detect unfavorable flow patterns in patients with re-CoA (22), and to plan interventions in complex scenarios (23). In the present study, a computational approach was chosen to isolate the effect of aortic arch geometry (based on representative anatomical configurations of the aortic arch derived from a statistical shape modeling framework, rather than individual patients' anatomies) and offer new insights into the hemodynamic effect of aortic arch architecture, particularly in the presence of a repaired CoA. Simulations revealed interesting differences in the repaired CoA setting. The idealized and purposely simplified scenario adopted here in the CFD simulations (i.e., same inflow, no aortic wall deformation) allowed us to isolate the effect of the geometry of CoA repair and determine the impact of that variable. Hemodynamic differences mainly included a nearly doubled value of vorticity in the descending aorta past the repaired CoA region, increased WSS in the descending aorta, a mild increase in pressure drop compared to the isolated BAV scenario, and an increase in power loss.
These findings suggest that even in a successfully repaired CoA (i.e., no residual narrowing with a clinically negligible CoA index and a qualitatively optimal aortic arch geometry), there are important changes in the architecture of the arch. The elevated WSS in the descending aorta may relate to some degree of dilation in this area that has been previously observed in a study focusing on aortic arch morphology in BAV with/without CoA (8). In addition, power loss (which is a function of the pressure drop and the aortic flow rate) was also found to be different in BAV patients with r-CoA compared to patients with isolated BAV. Power loss was higher in the presence of CoA (both repaired and unrepaired). This concurs with existing literature reporting an association between increasing stenosis level and increasing energy dissipation (24). Furthermore, previous work in other CHD, such as the geometry of the total cavopulmonary connection in Fontan patients, has shown that power loss increases dramatically in a non-linear fashion with increasing cardiac output and is dependent on geometrical variables, suggesting the importance of studying exercise conditions as well as rest conditions (25). While our simplified model focused on isolating differences related to aortic arch morphology, the observed differences in power loss, especially in the r-CoA scenario, warrant further study and extension to the investigation of the effect of exercise physiology, as this is likely to exacerbate any differences. Another scenario that would be very relevant to include in such an analysis would be the case of transverse arch hypoplasia, particularly in light of recent work, based on principal component analysis (PCA), suggesting that arch morphology is not the major determinant of vascular load following coarctation repair and that caliber is more important than curvature (7). These observations are purely geometric and did not consider effects on the flow field, so they can be viewed as complementary alongside those presented in this article, contributing to refining our knowledge of CoA. Additionally, if a growth model of the aorta with CoA was developed following a recent methodology (26), it would be interesting to include factors such as recurrent arch obstruction, which is more prevalent in patients with a CoA index <0.7, extending the observations from this study.
The flow separation results suggest that small variations in the morphology of the aortic arch (native or resulting from surgical repair) lead to changes in flow velocity acceleration and also to changes in WSS in the descending aorta. Our observations agree with previous literature suggesting a relation between geometric irregularities and flow separation at the corners and that this increases from romanesque to crenel to gothic arches (27). It was particularly interesting to observe flow separation in the scenario of successful CoA repair with no residual significant narrowing, which reinforces previous observations that even when recreating a favorable anatomy after CoA repair, aortic hemodynamics are still not fully restored (28).
This study indeed reinforces observations that CoA repair does not restore normal hemodynamics, even if there is no residual narrowing. These patients could be assumed to have an abnormal vasculature from birth, regardless of repair with on-going cardiovascular risk, even in the absence of residual narrowing or systemic hypertension. It is crucially important to pay attention to the control of cardiovascular risk factors throughout life, such as refraining from smoking, actively managing elevated blood pressure and cholesterol and the encouragement of a healthy weight and exercise habits from an early age. Certain arch phenotypes may be most at risk and warrant intensive preventative therapy.
There are limitations to this study. The aortic models were set as rigid and the same realistic aortic inflow was used as a boundary condition for all three geometrical configurations. Whilst this was purposely done in order to isolate the geometrical variable, aortic distensibility and patient-specific flow features contribute to the hemodynamic changes in these patients. We suggest that we can further augment our knowledge of CoA hemodynamics by studying these variables in a parametric fashion in future studies. Furthermore, differences in relevant parameters such as WSS are likely to be affected by changes in the patients' physiological state, i.e., rest vs. exercise (28), which was not modeled in this study. Previous MRI-based analysis suggesting that transverse arch and isthmus hypoplasia rather than acute aortic arch angulation are involved in the pathophysiology of blood pressure response to peak exercise after CoA repair (29) and it has been observed that CoA patients without significant obstruction had higher exercise-induced systolic blood pressure changes than matched controls (30). Finally, whilst this study focused on aortic hemodynamics, it is important to note that CFD simulations in this context can also allow the comparison of relevant indices between the aorta and the brachiocephalic vessels, which is an important consideration in CoA patients, given the incidence of cerebrovascular events (28, 31). Whilst events may be related to hypertension, the role of arch geometry and anatomy of the head and neck vessels may well be a contributory factor.
Conclusion
This study adds to our appreciation of aortic hemodynamics in patients with repaired aortic coarctation, suggesting that even small alterations in the aortic morphology (such as those resulting following surgical repair) impact on key hemodynamic indices. This can in part contribute to explaining phenomena such as persistent hypertension even in the absence of any clinically significant narrowing.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
VGA, MO, AG, and GB conceived the study. VGA and HG ran the computational simulations. FS provided the anatomical models. VGA and GB drafted the manuscript. All authors contributed to data interpretation and critically revising.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
VGA acknowledges support from the University of Bristol to carry out his Ph.D. Also, the authors gratefully acknowledge the support of the British Heart Foundation (Prof. Caputo's Chair and the Bristol BHF Accelerator Award) and the Bristol NIHR Biomedical Research Center (BRC).
References
1. Kaemmerer H. Aortic coarctation and interrupted aortic arch. In: Gatzoulis MA, Webb GD, Daubeney PEF, editors. Diagnosis and Management of Adult Congenital Heart Disease. 2nd ed. Churchill Livingstone (2011). p. 261–70. doi: 10.1016/B978-0-7020-3426-8.00036-8
2. Ou P, Celermajer DS, Raisky O, Jolivet O, Buyens F, Herment A, et al. Angular (gothic) aortic arch leads to enhanced systolic wave reflection, central aortic stiffness, and increased left ventricular mass late after aortic coarctation repair: evaluation with magnetic resonance flow mapping. J Thorac Cardiovasc Surg. (2008) 135:62–8. doi: 10.1016/j.jtcvs.2007.03.059
3. Ou P, Celermajer DS, Mousseaux E, Giron A, Aggoun Y, Szezepanski I, et al. Vascular remodeling after “successful” repair of coarctation: impact of aortic arch geometry. J Am Coll Cardiol. (2007) 49:883–90. doi: 10.1016/j.jacc.2006.10.057
4. Ou P, Bonnet D, Auriacombe L, Pedroni E, Balleux F, Sidi D, et al. Late systemic hypertension and aortic arch geometry after successful repair of coarctation of the aorta. Eur Heart J. (2004) 25:1853–59. doi: 10.1016/j.ehj.2004.07.021
5. de Divitiis M, Rubba P, Calabrò R. Arterial hypertension and cardiovascular prognosis after successful repair of aortic coarctation: a clinical model for the study of vascular function. Nutr Metab Cardiovasc Dis. (2005) 15:382–94. doi: 10.1016/j.numecd.2005.08.002
6. Donazzan L, Crepaz R, Stuefer J, Stellin G. Abnormalities of aortic arch shape, central aortic flow dynamics, and distensibility pre-dispose to hypertension after successful repair of aortic coarctation. World J Pediatr Congenit Heart Surg. (2014) 5:546–53. doi: 10.1177/2150135114551028
7. Quail MA, Segers P, Steeden JA, Muthurangu V. The aorta after coarctation repair–effects of calibre and curvature on arterial haemodynamics. J Cardiovasc Magn Res. (2019) 21:1–9. doi: 10.1186/s12968-019-0534-7
8. Sophocleous F, Biffi B, Milano EG, Bruse J, Caputo M, Rajakaruna C, et al. Aortic morphological variability in patients with bicuspid aortic valve and aortic coarctation. Eur J Cardiothorac Surg. (2019) 55:704–13. doi: 10.1093/ejcts/ezy339
9. Rodrigues JC, Jaring MF, Werndle MC, Mitrousi K, Lyen SM, Nightingale AK, et al. Repaired coarctation of the aorta, persistent arterial hypertension and the selfish brain. J Cardiovasc Magn Reson. (2019) 21:1–10. doi: 10.1186/s12968-019-0578-8
10. Lee MG, Babu-Narayan SV, Kempny A, Uebing A, Montanaro C, Shore DF, et al. Long-term mortality and cardiovascular burden for adult survivors of coarctation of the aorta. Heart. (2019) 105:1190–96. doi: 10.1136/heartjnl-2018-314257
11. Kenny D, Polson JW, Martin RP, Paton JFR, Wolf AR. Hypertension and coarctation of the aorta: an inevitable consequence of developmental pathophysiology. Hypertens Res. (2011) 34:543–7. doi: 10.1038/hr.2011.22
12. Iriart X, Laik J, Cremer A, Martin C, Pillois X, Jalal Z, et al. Predictive factors for residual hypertension following aortic coarctation stenting. J Clin Hypertens. (2019) 21:291–8. doi: 10.1111/jch.13452
13. Bruse JL, McLeod K, Biglino G, Ntsinjana HN, Capelli C, Hsia T-Y, et al. A statistical shape modelling framework to extract 3d shape biomarkers from medical imaging data: assessing arch morphology of repaired coarctation of the aorta. BMC Med Imaging. (2016) 16:1–19. doi: 10.1186/s12880-016-0142-z
14. Tan F, Wood N, Tabor G, Xu X. Comparison of LES of steady transitional flow in an idealized stenosed axisymmetric artery model with a RANS transitional model. J Biomech Eng. (2011) 133:051001. doi: 10.1115/1.4003782
15. Varghese SS, Frankel SH, Fischer PF. Modeling transition to turbulence in eccentric stenotic flows. J Biomech Eng. (2008) 130:014503. doi: 10.1115/1.2800832
16. Mittal R, Simmons S, Udaykumar H. Application of large-eddy simulation to the study of pulsatile flow in a modeled arterial stenosis. J Biomech Eng. (2001) 123:325–32. doi: 10.1115/1.1385840
17. Antiga L, Piccinelli M, Botti L, Ene-Iordache B, Remuzzi A, Steinman DA. An image-based modeling framework for patient-specific computational hemodynamics. Med Biol Eng Comput. (2008) 46:1097–112. doi: 10.1007/s11517-008-0420-1
18. Keshavarz-Motamed Z, Garcia J, Kadem L. Fluid dynamics of coarctation of the aorta and effect of bicuspid aortic valve. PLoS ONE. (2013) 8:e72394. doi: 10.1371/journal.pone.0072394
19. Wendell DC, Samyn MM, Cava JR, Ellwein LM, Krolikowski MM, Gandy KL, et al. Including aortic valve morphology in computational fluid dynamics simulations: initial findings and application to aortic coarctation. Med Eng Phys. (2013) 35:723–35. doi: 10.1016/j.medengphy.2012.07.015
20. LaDisa JF Jr, Dholakia RJ, Figueroa CA, Vignon-Clementel IE, Chan FP, Samyn MM, et al. Computational simulations demonstrate altered wall shear stress in aortic coarctation patients treated by resection with end-to-end anastomosis. Congenit Heart Dis. (2011) 6:432–43. doi: 10.1111/j.1747-0803.2011.00553.x
21. Caimi A, Pasquali M, Sturla F, Pluchinotta FR, Giugno L, Carminati M, et al. Prediction of post-stenting biomechanics in coarcted aortas: a pilot finite element study. J Biomech. (2020) 105:109796. doi: 10.1016/j.jbiomech.2020.109796
22. Guillot M, Ascuitto R, Ross-Ascuitto N, Mallula K, Siwik E. Computational fluid dynamics simulations as a complementary study for transcatheter endovascular stent implantation for re-coarctation of the aorta associated with minimal pressure drop: an aneurysmal ductal ampulla with aortic isthmus narrowing. Cardiol Young. (2019) 29:768–76. doi: 10.1017/S1047951119000751
23. Tossas-Betancourt C, van Bakel TM, Arthurs CJ, Coleman DM, Eliason JL, Figueroa CA, et al. Computational analysis of renal artery flow characteristics by modeling aortoplasty and aortic bypass interventions for abdominal aortic coarctation. J Vasc Surg. (2020) 71:505–16. e4. doi: 10.1016/j.jvs.2019.02.063
24. Yap C-H, Dasi LP, Yoganathan AP. Dynamic hemodynamic energy loss in normal and stenosed aortic valves. J Biomech Eng. (2010) 132:021005. doi: 10.1115/1.4000874
25. Whitehead KK, Pekkan K, Kitajima HD, Paridon SM, Yoganathan AP, Fogel MA. Non-linear power loss during exercise in single-ventricle patients after the fontan: insights from computational fluid dynamics. Circulation. (2007) 116:I165–71. doi: 10.1161/CIRCULATIONAHA.106.680827
26. Sophocleous F, Bone A, Shearn AIU, Velasco Forte MN, Bruse JL, Caputo M, et al. Feasibility of a longitudinal statistical atlas model to study aortic growth in congenital heart disease. Comput Biol Med. (2022) 144:105326. doi: 10.1016/j.compbiomed.2022.105326
27. Olivieri LJ, de Zélicourt DA, Haggerty CM, Ratnayaka K, Cross RR, Yoganathan AP. Hemodynamic modeling of surgically repaired coarctation of the aorta. Cardiovasc Eng Technol. (2011) 2:288–95. doi: 10.1007/s13239-011-0059-1
28. LaDisa JF, Alberto Figueroa C, Vignon-Clementel IE, Jin Kim H, Xiao N, Ellwein LM, et al. Computational simulations for aortic coarctation: representative results from a sampling of patients. J Biomech Eng. (2011) 133:091008. doi: 10.1115/1.4004996
29. Ntsinjana HN, Biglino G, Capelli C, Tann O, Giardini A, Derrick G, et al. Aortic arch shape is not associated with hypertensive response to exercise in patients with repaired congenital heart diseases. J Cardiovasc Magn Res. (2013) 15:1–7. doi: 10.1186/1532-429X-15-101
30. Egbe AC, Allison TG, Ammash NM. Mild coarctation of aorta is an independent risk factor for exercise-induced hypertension. Hypertension. (2019) 74:1484–9. doi: 10.1161/HYPERTENSIONAHA.119.13726
Keywords: aortic coarctation, aortic hemodynamics, computational modeling, computational fluid dynamics, power loss, wall shear stress
Citation: Goodarzi Ardakani V, Goordoyal H, Ordonez MV, Sophocleous F, Curtis S, Bedair R, Caputo M, Gambaruto A and Biglino G (2022) Isolating the Effect of Arch Architecture on Aortic Hemodynamics Late After Coarctation Repair: A Computational Study. Front. Cardiovasc. Med. 9:855118. doi: 10.3389/fcvm.2022.855118
Received: 14 January 2022; Accepted: 30 May 2022;
Published: 24 June 2022.
Edited by:
Ruth Heying, University Hospital Leuven, BelgiumReviewed by:
Mario Carminati, IRCCS San Donato Polyclinic, ItalyDaniel De Wolf, Ghent University Hospital, Belgium
Copyright © 2022 Goodarzi Ardakani, Goordoyal, Ordonez, Sophocleous, Curtis, Bedair, Caputo, Gambaruto and Biglino. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giovanni Biglino, Zy5iaWdsaW5vJiN4MDAwNDA7YnJpc3RvbC5hYy51aw==