Model-Based Fluid-Structure Interaction Approach for Evaluation of Thoracic Endovascular Aortic Repair Endograft Length in Type B Aortic Dissection

Impact of Residual Intimal Flap Displacement Post-TEVAR on TBAD Haemodynamics in Compliant, Patient-specific CFD Simulations Informed by MRI

We propose a novel formulation of a moving boundary method to account for the motion of the intimal flap (IF) in a TBAD post-thoracic endovascular aortic repair using patient-specific compliant computational fluid dynamics simulations. The simulations were informed by non-invasive 4D flow MRI sequences. Predicted flow waveforms, aortic wall, and IF displacements were validated against in vivo 4D flow MRI and cine-MRI data. The patient-specific simulation showed that at peak systole, the dynamic interplay between high IF displacement and high transmural pressures promoted true lumen compression and false lumen expansion, whilst luminal patterns were reversed at the deceleration phase. High vorticity and swirling flow patterns were observed throughout the cardiac cycle at the primary entry tear, the descending aorta and proximal to the visceral aortic branches, correlating with high relative residence time, which could indicate an increased localised risk of aortic growth proximal to the IF. A rigid IF simulation revealed significant discrepancies in haemodynamic metrics, highlighting the potential mispredictions when using a rigid wall assumption to assess disease progression. Simulations assuming a more compliant IF highlighted potential increased risks of visceral branches malperfusion and localised aortic wall degeneration. The study underscores the necessity of patient-specific compliant IF simulations for accurate TBAD haemodynamic assessments. These insights can improve disease understanding and inform future treatment strategies.

Mechanisms of aortic dissection: From pathological changes to experimental and in silico models

Aortic dissection continues to be responsible for significant morbidity and mortality, although recent advances in medical data assimilation and in experimental and in silico models have improved our understanding of the initiation and progression of the accumulation of blood within the aortic wall. Hence, there remains a pressing necessity for innovative and enhanced models to more accurately characterize the associated pathological changes. Early on, experimental models were employed to uncover mechanisms in aortic dissection, such as hemodynamic changes and alterations in wall microstructure, and to assess the efficacy of medical implants. While experimental models were once the only option available, more recently they are also being used to validate in silico models. Based on an improved understanding of the deteriorated microstructure of the aortic wall, numerous multiscale material models have been proposed in recent decades to study the state of stress in dissected aortas, including the changes associated with damage and failure. Furthermore, when integrated with accessible patient-derived medical data, in silico models prove to be an invaluable tool for identifying correlations between hemodynamics, wall stresses, or thrombus formation in the deteriorated aortic wall. They are also advantageous for model-guided design of medical implants with the aim of evaluating the deployment and migration of implants in patients. Nonetheless, the utility of in silico models depends largely on patient-derived medical data, such as chosen boundary conditions or tissue properties. In this review article, our objective is to provide a thorough summary of medical data elucidating the pathological alterations associated with this disease. Concurrently, we aim to assess experimental models, as well as multiscale material and patient data-informed in silico models, that investigate various aspects of aortic dissection. In conclusion, we present a discourse on future perspectives, encompassing aspects of disease modeling, numerical challenges, and clinical applications, with a particular focus on aortic dissection. The aspiration is to inspire future studies, deepen our comprehension of the disease, and ultimately shape clinical care and treatment decisions.

Aortic stretch and recoil create wave-pumping effect: the second heart in the systemic circulation

Wave propagation in the heart tube is key to establishing an early pumping mechanism, as explained by impedance pump theory in zebrafish. Though initially proposed for embryonic blood circulation, the role of impedance-like behaviour in the mature cardiovascular system remains unclear. This study focuses on the understudied physiological mechanism of longitudinal displacement in the adult aorta caused by the long-axis motion of the heart. Using magnetic resonance imaging on 159 individuals, we compared aortic displacement profiles between a control group and those with heart failure, revealing a significant difference in aortic stretch between the two groups. Building on this clinical evidence, we conducted in vitro experiments to isolate the effects of longitudinal aortic wave pumping by eliminating the pumping action of the heart. We identified three biomechanical properties of stretch-related longitudinal wave pumping that exhibit characteristics like impedance pump: (i) a nonlinear flow-frequency relationship, (ii) bidirectional flow, and (iii) the potential for both positive and negative flow at a fixed frequency, contingent upon the aorta's wave speed dictating the wave state. Our results demonstrate for the first time that this mechanism generates a significant flow, potentially providing a supplementary pumping mechanism for the heart.

Long-Term Longitudinal Computational Study of a Marfan Syndrome Patient After Hybrid Repair of Aortic Arch Dissection With Parallel Stent-Grafts

Hybrid repair is a valuable alternative treatment for aortic arch disease in Marfan syndrome patients after proximal aorta replacement. This study aimed to investigate the long-term durability of this technique with the use of parallel stent-grafts and evaluate strategies to prevent abdominal aortic dilation. One Marfan syndrome patient who underwent hybrid aortic repair with parallel stent-grafts for arch dissection after the Bentall procedure was admitted. Five patient-specific three-dimensional models were reconstructed based on preoperative and follow-up computed tomography angiography scans. Three hypothetical models addressing the closure of an endoleak or reentry tears were created. Hemodynamic parameters were assessed using computational fluid dynamics. Postoperatively, increased blood flow into the descending aorta and rising abdominal aortic pressure were observed. During the 5-year follow-up, no new thoracic aorta-related adverse events occurred. One early type III endoleak persisted, and three reentry tears were identified in the descending aorta. The abdominal aorta dilated from 31 to 49 mm. Simultaneously addressing both the endoleak and reentry tears was more effective in reducing false lumen pressure and flow velocity in the abdominal aorta and expanding the high-value relative residence time region. Longitudinal follow-up imaging demonstrated the long-term durability of hybrid aortic arch repair with parallel stent-grafts in a Marfan syndrome patient after ascending aorta replacement. The increased pressure resulting from blood flow redistribution was associated with downstream aortic dilation. Furthermore, computational fluid dynamics simulations can offer predictive analyses for optimizing intervention strategies in the treatment of distal aneurysmal degeneration.

The impact of blood viscosity modeling on computational fluid dynamic simulations of pediatric patients with Fontan circulation

For univentricular heart patients, the Fontan circulation presents a unique pathophysiology due to chronic non-pulsatile low-shear-rate pulmonary blood flow, where non-Newtonian effects are likely substantial. This study evaluates the influence of non-Newtonian behavior of blood on fluid dynamics and energetic efficiency in pediatric patient-specific models of the Fontan circulation. We used immersed boundary-lattice Boltzmann method simulations to compare Newtonian and non-Newtonian viscosity models. The study included models from twenty patients exhibiting a low cardiac output state (cardiac index of 2 L/min/m2). We quantified metrics of energy loss (indexed power loss and viscous dissipation), non-Newtonian importance factors, and hepatic flow distribution. We observed significant differences in flow structure between Newtonian and non-Newtonian models. Specifically, the non-Newtonian simulations demonstrated significantly higher local and average viscosity, corresponding to a higher non-Newtonian importance factor and larger energy loss. Hepatic flow distribution was also significantly different in a subset of patients. These findings suggest that non-Newtonian behavior contributes to flow structure and energetic inefficiency in the low cardiac output state of the Fontan circulation.

Longitudinal investigation of aortic dissection in mice with computational fluid dynamics

Predicting late adverse events in aortic dissections is challenging. One commonly observed risk factor is partial thrombosis of the false lumen. In this study we investigated false lumen thrombus progression over 7 days in four mice with angiotensin II-induced aortic dissection. We performed computational fluid dynamic simulations with subject-specific boundary conditions from velocity and pressure measurements. We investigated endothelial cell activation potential, mean velocity, thrombus formation potential, and other hemodynamic factors. Our findings support the hypothesis that flow stagnation is the predominant hemodynamic factor driving a large thrombus ratio in false lumina, particularly those with a single fenestration.

A preoperative planning procedure of septal myectomy for hypertrophic obstructive cardiomyopathy using image-based computational fluid dynamics simulations and shape optimization

Although septal myectomy is the preferred treatment for medication-refractory hypertrophic obstructive cardiomyopathy (HOCM), the procedure remains subjective. A preoperative planning procedure is proposed using computational fluid dynamics simulations and shape optimization to assist in the objective assessment of the adequacy of the resection. 3 patients with HOCM were chosen for the application of the proposed procedure. The geometries of the preoperative left ventricular outflow tract (LVOT) of patients in the systolic phase were reconstructed from medical images. Computaional fluid dynamics (CFD) simulations were performed to assess hemodynamics within LVOT. Sensitivity analysis was performed to determine the resection extent on the septal wall, and the depth of the resection was optimized to relieve LVOT obstruction while minimizing damage to the septum. The optimized resection was then transferred from systole to diastole to provide surgeons with instructive guidance for septal myectomy. Comparison between preoperative and postoperative hemodynamics showed an evident improvement with respect to the pressure gradient throughout the LVOT. The resected myocardium in the diastolic phase is more extended and thinner than its state in the systolic phase. The proposed preoperative planning procedure may be a viable addition to the current preoperative assessment of patients with HOCM.

Patient-specific compliant simulation framework informed by 4DMRI-extracted pulse wave Velocity: Application post-TEVAR

We introduce a new computational framework that utilises Pulse Wave Velocity (PWV) extracted directly from 4D flow MRI (4DMRI) to inform patient-specific compliant computational fluid dynamics (CFD) simulations of a Type-B aortic dissection (TBAD), post-thoracic endovascular aortic repair (TEVAR). The thoracic aortic geometry, a 3D inlet velocity profile (IVP) and dynamic outlet boundary conditions are derived from 4DMRI and brachial pressure patient data. A moving boundary method (MBM) is applied to simulate aortic wall displacement. The aortic wall stiffness is estimated through two methods: one relying on area-based distensibility and the other utilising regional pulse wave velocity (RPWV) distensibility, further fine-tuned to align with in vivo values. Predicted pressures and outlet flow rates were within 2.3 % of target values. RPWV-based simulations were more accurate in replicating in vivo hemodynamics than the area-based ones. RPWVs were closely predicted in most regions, except the endograft. Systolic flow reversal ratios (SFRR) were accurately captured, while differences above 60 % in in-plane rotational flow (IRF) between the simulations were observed. Significant disparities in predicted wall shear stress (WSS)-based indices were observed between the two approaches, especially the endothelial cell activation potential (ECAP). At the isthmus, the RPWV-driven simulation indicated a mean ECAP>1.4 Pa-1 (critical threshold), indicating areas potentially prone to thrombosis, not captured by the area-based simulation. RPWV-driven simulation results agree well with 4DMRI measurements, validating the proposed pipeline and facilitating a comprehensive assessment of surgical decision-making scenarios and potential complications, such as thrombosis and aortic growth.

An efficient procedure for the blood flow computer simulation of patient-specific aortic dissections

In this work we present a novel methodology for the numerical simulation of patient-specific aortic dissections. Our proposal, which targets the seamless virtual prototyping of customized scenarios, combines an innovative two-step segmentation procedure with a CutFEM technique capable of dealing with thin-walled bodies such as the intimal flap. First, we generate the fluid mesh from the outer aortic wall disregarding the intimal flap, similarly to what would be done in a healthy aorta. Second, we create a surface mesh from the approximate midline of the intimal flap. This approach allows us to decouple the segmentation of the fluid volume from that of the intimal flap, thereby bypassing the need to create a volumetric mesh around a thin-walled body, an operation widely known to be complex and error-prone. Once the two meshes are obtained, the original configuration of the dissection into true and false lumen is recovered by embedding the surface mesh into the volumetric one and calculating a level set function that implicitly represents the intimal flap in terms of the volumetric mesh entities. We then leverage the capabilities of unfitted mesh methods, specifically relying on a CutFEM technique tailored for thin-walled bodies, to impose the wall boundary conditions over the embedded intimal flap. We tested the method by simulating the flow in four patient-specific aortic dissections, all involving intricate geometrical patterns. In all cases, the preprocess is greatly simplified with no impact on the computational times. Additionally, the obtained results are consistent with clinical evidence and previous research.

A new methodology for determining the central pressure waveform from peripheral measurement using Fourier-based machine learning

Radial applanation tonometry is a well-established technique for hemodynamic monitoring and is becoming popular in affordable non-invasive wearable healthcare electronics. To assess the central aortic pressure using radial-based measurements, there is an essential need to develop mathematical approaches to estimate the central pressure waveform. In this study, we propose a new Fourier-based machine learning (F-ML) methodology to transfer non-invasive radial pressure measurements to the central pressure waveform. To test the method, collection of tonometry recordings of the radial and carotid pressure measurements are used from the Framingham Heart Study (2640 individuals, 55 % women) with mean (range) age of 66 (40-91) years. Method-derived estimates are significantly correlated with the measured ones for three major features of the pressure waveform (systolic blood pressure, r=0.97, p < 0.001; diastolic blood pressure, r=0.99, p < 0.001; and mean blood pressure, r=0.99, p < 0.001). In all cases, the Bland-Altman analysis shows negligible bias in the estimations and error is bounded to 5.4 mmHg. Findings also suggest that the F-ML approach reconstructs the shape of the central pressure waveform accurately with the average normalized root mean square error of 5.5 % in the testing population which is blinded to all stages of machine learning development. The results show that the F-ML transfer function outperforms the conventional generalized transfer function, particularly in terms of reconstructing the shape of the central pressure waveform morphology. The proposed F-ML transfer function can provide accurate estimates for the central pressure waveform, and ultimately expand the usage of non-invasive devices for central hemodynamic assessment.

Patient-Specific Haemodynamic Analysis of Virtual Grafting Strategies in Type-B Aortic Dissection: Impact of Compliance Mismatch

INTRODUCTION

Compliance mismatch between the aortic wall and Dacron Grafts is a clinical problem concerning aortic haemodynamics and morphological degeneration. The aortic stiffness introduced by grafts can lead to an increased left ventricular (LV) afterload. This study quantifies the impact of compliance mismatch by virtually testing different Type-B aortic dissection (TBAD) surgical grafting strategies in patient-specific, compliant computational fluid dynamics (CFD) simulations.

MATERIALS AND METHODS

A post-operative case of TBAD was segmented from computed tomography angiography data. Three virtual surgeries were generated using different grafts; two additional cases with compliant grafts were assessed. Compliant CFD simulations were performed using a patient-specific inlet flow rate and three-element Windkessel outlet boundary conditions informed by 2D-Flow MRI data. The wall compliance was calibrated using Cine-MRI images. Pressure, wall shear stress (WSS) indices and energy loss (EL) were computed.

RESULTS

Increased aortic stiffness and longer grafts increased aortic pressure and EL. Implementing a compliant graft matching the aortic compliance of the patient reduced the pulse pressure by 11% and EL by 4%. The endothelial cell activation potential (ECAP) differed the most within the aneurysm, where the maximum percentage difference between the reference case and the mid (MDA) and complete (CDA) descending aorta replacements increased by 16% and 20%, respectively.

CONCLUSION

This study suggests that by minimising graft length and matching its compliance to the native aorta whilst aligning with surgical requirements, the risk of LV hypertrophy may be reduced. This provides evidence that compliance-matching grafts may enhance patient outcomes.

Assessing pressure wave components for aortic stiffness monitoring through spectral regression learning

Aims

The ageing process notably induces structural changes in the arterial system, primarily manifesting as increased aortic stiffness, a precursor to cardiovascular events. While wave separation analysis is a robust tool for decomposing the components of blood pressure waveform, its relationship with cardiovascular events, such as aortic stiffening, is incompletely understood. Furthermore, its applicability has been limited due to the need for concurrent measurements of pressure and flow. Our aim in this study addresses this gap by introducing a spectral regression learning method for pressure-only wave separation analysis.

Methods and results

Leveraging data from the Framingham Heart Study (2640 individuals, 55% women), we evaluate the accuracy of pressure-only estimates, their interchangeability with a reference method based on ultrasound-derived flow waves, and their association with carotid-femoral pulse wave velocity (PWV). Method-derived estimates are strongly correlated with the reference ones for forward wave amplitude ( R 2 = 0.91 ), backward wave amplitude ( R 2 = 0.88 ), and reflection index ( R 2 = 0.87 ) and moderately correlated with a time delay between forward and backward waves ( R 2 = 0.38 ). The proposed pressure-only method shows interchangeability with the reference method through covariate analysis. Adjusting for age, sex, body size, mean blood pressure, and heart rate, the results suggest that both pressure-only and pressure-flow evaluations of wave separation parameters yield similar model performances for predicting carotid-femoral PWV, with forward wave amplitude being the only significant factor (P < 0.001; 95% confidence interval, 0.056-0.097).

Conclusion

We propose an interchangeable pressure-only wave separation analysis method and demonstrate its clinical applicability in capturing aortic stiffening. The proposed method provides a valuable non-invasive tool for assessing cardiovascular health.

Stress Load and Ascending Aortic Aneurysms: An Observational, Longitudinal, Single-Center Study Using Computational Fluid Dynamics

Ascending aortic aneurysm (AAoA) is a silent disease with high mortality; however, the factors associated with a worse prognosis are not completely understood. The objective of this observational, longitudinal, single-center study was to identify the hemodynamic patterns and their influence on AAoA growth using computational fluid dynamics (CFD), focusing on the effects of geometrical variations on aortic hemodynamics. Personalized anatomic models were obtained from angiotomography scans of 30 patients in two different years (with intervals of one to three years between them), of which 16 (53%) showed aneurysm growth (defined as an increase in the ascending aorta volume by 5% or more). Numerically determined velocity and pressure fields were compared with the outcome of aneurysm growth. Through a statistical analysis, hemodynamic characteristics were found to be associated with aneurysm growth: average and maximum high pressure (superior to 100 Pa); average and maximum high wall shear stress (superior to 7 Pa) combined with high pressure (>100 Pa); and stress load over time (maximum pressure multiplied by the time interval between the exams). This study provides insights into a worse prognosis of this serious disease and may collaborate for the expansion of knowledge about mechanobiology in the progression of AAoA.

Mechanistic insights on age-related changes in heart-aorta-brain hemodynamic coupling using a pulse wave model of the entire circulatory system

Age-related changes in aortic biomechanics can impact the brain by reducing blood flow and increasing pulsatile energy transmission. Clinical studies have shown that impaired cardiac function in patients with heart failure is associated with cognitive impairment. Although previous studies have attempted to elucidate the complex relationship between age-associated aortic stiffening and pulsatility transmission to the cerebral network, they have not adequately addressed the effect of interactions between aortic stiffness and left ventricle (LV) contractility (neither on energy transmission nor on brain perfusion). In this study, we use a well-established and validated one-dimensional blood flow and pulse wave computational model of the circulatory system to address how age-related changes in cardiac function and vasculature affect the underlying mechanisms involved in the LV-aorta-brain hemodynamic coupling. Our results reveal how LV contractility affects pulsatile energy transmission to the brain, even with preserved cardiac output. Our model demonstrates the existence of an optimal heart rate (near the normal human heart rate) that minimizes pulsatile energy transmission to the brain at different contractility levels. Our findings further suggest that the reduction in cerebral blood flow at low levels of LV contractility is more prominent in the setting of age-related aortic stiffening. Maintaining optimal blood flow to the brain requires either an increase in contractility or an increase in heart rate. The former consistently leads to higher pulsatile power transmission, and the latter can either increase or decrease subsequent pulsatile power transmission to the brain.NEW & NOTEWORTHY We investigated the impact of major aging mechanisms of the arterial system and cardiac function on brain hemodynamics. Our findings suggest that aging has a significant impact on heart-aorta-brain coupling through changes in both arterial stiffening and left ventricle (LV) contractility. Understanding the underlying physical mechanisms involved here can potentially be a key step for developing more effective therapeutic strategies that can mitigate the contributions of abnormal LV-arterial coupling toward neurodegenerative diseases and dementia.

On the Longitudinal Wave Pumping in Fluid-filled Compliant Tubes

This study investigates the physics of the longitudinal stretching-based wave pumping mechanism, a novel extension of the traditional impedance pump. In its simplest form, an impedance pump consists of a fluid-filled elastic tube connected to rigid tubes with a wave generator. These valveless pumps operate based on the principles of wave propagation in a fluid-filled compliant tube. Cardiovascular magnetic resonance imaging of the human circulatory system has shown substantial stretching of the aorta (the largest compliant artery of the body carrying blood) during the heart contraction and recoil of the aorta during the relaxation. Inspired by this dynamic mechanism, a comprehensive analysis of a longitudinal impedance pump is conducted in this study where waves are generated by stretching of the elastic wall and its recoil. We developed a fully coupled fluid-structure interaction computational model consisting of a straight fluid-filled elastic tube with longitudinal stretch at one end and fixed reflection site at the other end. The pump's behavior is quantified as a function of stretching frequency and tube wall characteristics. Our results indicate that stretch-related wave propagation and reflection can induce frequency-dependent pumping. Findings suggest a non-linear pattern for the mean flow-frequency relationship. Based on the analysis of the propagated waveforms, the underlying physical mechanism in the longitudinal impedance pump is discussed. It is shown that both the direction and magnitude of the net flow strongly depend on the wave characteristics. These findings provide a fundamental understanding of stretch-related wave pumping and can inform the future design of such pumps.

Direct 0D-3D coupling of a lattice Boltzmann methodology for fluid-structure aortic flow simulations

This work introduces a numerical approach and implementation for the direct coupling of arbitrary complex ordinary differential equation- (ODE-)governed zero-dimensional (0D) boundary conditions to three-dimensional (3D) lattice Boltzmann-based fluid-structure systems for hemodynamics studies. In particular, a most complex configuration is treated by considering a dynamic left ventricle- (LV-)elastance heart model which is governed by (and applied as) a nonlinear, non-stationary hybrid ODE-Dirichlet system. Other ODE-based boundary conditions, such as lumped parameter Windkessel models for truncated vasculature, are also considered. Performance studies of the complete 0D-3D solver, including its treatment of the lattice Boltzmann fluid equations and elastodynamics equations as well as their interactions, is conducted through a variety of benchmark and convergence studies that demonstrate the ability of the coupled 0D-3D methodology in generating physiological pressure and flow waveforms-ultimately enabling the exploration of various physical and physiological parameters for hemodynamics studies of the coupled LV-arterial system. The methods proposed in this paper can be easily applied to other ODE-based boundary conditions as well as to other fluid problems that are modeled by 3D lattice Boltzmann equations and that require direct coupling of dynamic 0D boundary conditions.