Computational fluid dynamics in cardiac surgery and perfusion: A review.

Cardiovascular diseases persist as a leading cause of mortality and morbidity, despite significant advances in diagnostic and surgical approaches. Computational Fluid Dynamics (CFD) represents a branch of fluid mechanics widely used in industrial engineering but is increasingly applied to the cardiovascular system. This review delves into the transformative potential for simulating cardiac surgery procedures and perfusion systems, providing an in-depth examination of the state-of-the-art in cardiovascular CFD modeling. The study first describes the rationale for CFD modeling and later focuses on the latest advances in heart valve surgery, transcatheter heart valve replacement, aortic aneurysms, and extracorporeal membrane oxygenation. The review underscores the role of CFD in better understanding physiopathology and its clinical relevance, as well as the profound impact of hemodynamic stimuli on patient outcomes. By integrating computational methods with advanced imaging techniques, CFD establishes a quantitative framework for understanding the intricacies of the cardiac field, providing valuable insights into disease progression and treatment strategies. As technology advances, the evolving synergy between computational simulations and clinical interventions is poised to revolutionize cardiovascular care. This collaboration sets the stage for more personalized and effective therapeutic strategies. With its potential to enhance our understanding of cardiac pathologies, CFD stands as a promising tool for improving patient outcomes in the dynamic landscape of cardiovascular medicine.

Biomechanical performance of the Bicaval Transcatheter System for the treatment of severe tricuspid regurgitation.

Introduction: Tricuspid regurgitation (TR) is a relatively common valvular disease, which can result from structural abnormalities of any anatomic part of the tricuspid valve. Severe TR is linked to congestive heart failure and hemodynamic impairment, resulting in high mortality when repaired by elective surgery. This study was undertaken to quantify the structural and hemodynamic performance of the novel Transcatheter Bicaval Valves System (TricValve) percutaneously implanted in the superior vena cava (SVC) and inferior vena cava (IVC) of two patients with severe TR and venous congestion. Methods: After developing the SVC and IVC device models, the contact pressure exerted on the vena cava wall was obtained by computational analysis. Both smoothed-particle hydrodynamics (SPH) and computational fluid dynamics were carried out to quantify caval reflux in the right atrium and the pressure field of pre- and post-TricValve scenarios, respectively. Results: Analysis of contact pressure highlighted the main anchoring area of the SVC device occurring near the SVC device belly, while the IVC device exerted pronounced forces in the device's proximal and distal parts. SPH-related flow velocities revealed the absence of caval reflux, and a decrease in time-averaged pressure was observed near the SVC and IVC after TricValve implantation. Discussion: Findings demonstrated the potential of computational tools for enhancing our understanding of the biomechanical performance of structural tricuspid valve interventions and improving the way we design next-generation transcatheter therapies to treat the tricuspid valve with heterotopic caval valve implantation.

Patient-Specific Analysis of Ascending Thoracic Aortic Aneurysm with the Living Heart Human Model.

In ascending thoracic aortic aneurysms (ATAAs), aneurysm kinematics are driven by ventricular traction occurring every heartbeat, increasing the stress level of dilated aortic wall. Aortic elongation due to heart motion and aortic length are emerging as potential indicators of adverse events in ATAAs; however, simulation of ATAA that takes into account the cardiac mechanics is technically challenging. The objective of this study was to adapt the realistic Living Heart Human Model (LHHM) to the anatomy and physiology of a patient with ATAA to assess the role of cardiac motion on aortic wall stress distribution. Patient-specific segmentation and material parameter estimation were done using preoperative computed tomography angiography (CTA) and ex vivo biaxial testing of the harvested tissue collected during surgery. The lumped-parameter model of systemic circulation implemented in the LHHM was refined using clinical and echocardiographic data. The results showed that the longitudinal stress was highest in the major curvature of the aneurysm, with specific aortic quadrants having stress levels change from tensile to compressive in a transmural direction. This study revealed the key role of heart motion that stretches the aortic root and increases ATAA wall tension. The ATAA LHHM is a realistic cardiovascular platform where patient-specific information can be easily integrated to assess the aneurysm biomechanics and potentially support the clinical management of patients with ATAAs.

Statistical Shape Analysis of Ascending Thoracic Aortic Aneurysm: Correlation between Shape and Biomechanical Descriptors.

An ascending thoracic aortic aneurysm (ATAA) is a heterogeneous disease showing different patterns of aortic dilatation and valve morphologies, each with distinct clinical course. This study aimed to explore the aortic morphology and the associations between shape and function in a population of ATAA, while further assessing novel risk models of aortic surgery not based on aortic size. Shape variability of n = 106 patients with ATAA and different valve morphologies (i.e., bicuspid versus tricuspid aortic valve) was estimated by statistical shape analysis (SSA) to compute a mean aortic shape and its deformation. Once the computational atlas was built, principal component analysis (PCA) allowed to reduce the complex ATAA anatomy to a few shape modes, which were correlated to shear stress and aortic strain, as determined by computational analysis. Findings demonstrated that shape modes are associated to specific morphological features of aneurysmal aorta as the vessel tortuosity and local bulging of the ATAA. A predictive model, built with principal shape modes of the ATAA wall, achieved better performance in stratifying surgically operated ATAAs versus monitored ATAAs, with respect to a baseline model using the maximum aortic diameter. Using current imaging resources, this study demonstrated the potential of SSA to investigate the association between shape and function in ATAAs, with the goal of developing a personalized approach for the treatment of the severity of aneurysmal aorta.

Shear Stress and Aortic Strain Associations With Biomarkers of Ascending Thoracic Aortic Aneurysm.

BACKGROUND: This study aims to investigate the association of wall shear stress (WSS) and aortic strain with circulating biomarkers including matrix metalloproteinases (MMP), tissue inhibitors of metalloproteinase (TIMP), and exosomal level of microRNA (miRNA) in ascending aortic aneurysms of patients with bicuspid or tricuspid aortic valve. METHODS: A total of 76 variables from 125 patients with ascending aortic aneurysms were collected from (1) blood plasma to measure plasma levels of miRNAs and protein activity; (2) computational flow analysis to estimate peak systolic WSS and time-average WSS (TAWSS); and (3) imaging analysis of computed tomography angiography to determine aortic wall strain. Principal component analysis followed by logistic regression allowed the development of a predictive model of aortic surgery by combining biomechanical descriptors and biomarkers. RESULTS: The protein activity of MMP-1, TIMP-1, and MMP-2 was positively correlated to the systolic WSS and TAWSS observed in the proximal ascending aorta (eg, R = 0.52, P < .001, for MMP-1 with TAWSS) where local maxima of WSS were found. For bicuspid patients, aortic wall strain was associated with miR-26a (R = 0.55, P = .041) and miR-320a (R = 0.69, P < .001), which shows a significant difference between bicuspid and tricuspid patients. Receiver-operating characteristics curves revealed that the combination of WSS, MMP-1, TIMP-1, and MMP-12 is predictive of aortic surgery (area under the curve 0.898). CONCLUSIONS: Increased flow-based and structural descriptors of ascending aortic aneurysms are associated with high levels of circulating biomarkers, implicating adverse vascular remodeling in the dilated aorta by mechanotransduction. A combination of shear stress and circulating biomarkers has the potential to improve the decision-making process for ascending aortic aneurysms to a highly individualized level.

Simulation study of transcatheter heart valve implantation in patients with stenotic bicuspid aortic valve.

Bicuspid aortic valve (BAV) anatomy has routinely been considered an exclusion in the setting of transcatheter aortic valve implantation (TAVI) because of the large dimension of the aortic annulus having a more calcified, bulky, and irregular shape. The study aims to develop a patient-specific computational framework to virtually simulate TAVI in stenotic BAV patients using the Edwards SAPIEN 3 valve (S3) and its improved version SAPIEN 3 Ultra and quantify stent frame deformity as well as the severity of paravalvular leakage (PVL). Specifically, the aortic root anatomy of n.9 BAV patients who underwent TAVI was reconstructed from pre-operative CT imaging. Crimping and deployment of S3 frame were performed and then followed by fluid-solid interaction analysis to simulate valve leaflet dynamics throughout the entire cardiac cycle. Modeling revealed that the S3 stent frame expanded well on BAV anatomy with an elliptical shape at the aortic annulus. Comparison of predicted S3 deformity as assessed by eccentricity and expansion indices demonstrated a good agreement with the measurement obtained from CT imaging. Blood particle flow analysis demonstrated a backward blood jet during diastole, whereas the predicted PVL flows corresponded well with those determined by transesophageal echocardiography. This study represents a further step towards the use of personalized simulations to virtually plan TAVI, aiming at improving not only the efficacy of the implantation but also the exploration of "off-label" applications as the TAVI in the setting of BAV patients. Graphical abstract Computational frameworks of TAVI in patients with stenotic bicuspid aortic valve.

Usefulness of regional right ventricular and right atrial strain for prediction of early and late right ventricular failure following a left ventricular assist device implant: A machine learning approach.

BACKGROUND: Identifying candidates for left ventricular assist device surgery at risk of right ventricular failure remains difficult. The aim was to identify the most accurate predictors of right ventricular failure among clinical, biological, and imaging markers, assessed by agreement of different supervised machine learning algorithms. METHODS: Seventy-four patients, referred to HeartWare left ventricular assist device since 2010 in two Italian centers, were recruited. Biomarkers, right ventricular standard, and strain echocardiography, as well as cath-lab measures, were compared among patients who did not develop right ventricular failure (N = 56), those with acute-right ventricular failure (N = 8, 11%) or chronic-right ventricular failure (N = 10, 14%). Logistic regression, penalized logistic regression, linear support vector machines, and naïve Bayes algorithms with leave-one-out validation were used to evaluate the efficiency of any combination of three collected variables in an "all-subsets" approach. RESULTS: Michigan risk score combined with central venous pressure assessed invasively and apical longitudinal systolic strain of the right ventricular-free wall were the most significant predictors of acute-right ventricular failure (maximum receiver operating characteristic-area under the curve = 0.95, 95% confidence interval = 0.91-1.00, by the naïve Bayes), while the right ventricular-free wall systolic strain of the middle segment, right atrial strain (QRS-synced), and tricuspid annular plane systolic excursion were the most significant predictors of Chronic-RVF (receiver operating characteristic-area under the curve = 0.97, 95% confidence interval = 0.91-1.00, according to naïve Bayes). CONCLUSION: Apical right ventricular strain as well as right atrial strain provides complementary information, both critical to predict acute-right ventricular failure and chronic-right ventricular failure, respectively.

Patterns of ascending aortic dilatation and predictors of surgical replacement of the aorta: A comparison of bicuspid and tricuspid aortic valve patients over eight years of follow-up.

BACKGROUND: Predictors of thoracic aorta growth and early cardiac surgery in patients with bicuspid aortic valve are undefined. Our aim was to identify predictors of ascending aorta dilatation and cardiac surgery in patients with bicuspid aortic valve (BAV). METHODS: Forty-one patients with BAV were compared with 165 patients with tricuspid aortic valve (TAV). All patients had LV EF > 50%, normal LV dimensions, and similar degree of aortic root or ascending aorta dilatation at enrollment. Patients with more than mild aortic stenosis or regurgitation were excluded. A CT-scan was available on 76% of the population, and an echocardiogram was repeated every year for a median time of 4 years (range: 2 to 8 years). Patterns of aortic expansion in BAV and TAV groups were analyzed by a mixed-effects longitudinal linear model. In the time-to-event analysis, the primary end point was elective or emergent surgery for aorta replacement. RESULTS: BAV patients were younger, while the TAV group had greater LV wall thickness, arterial hypertension, and dyslipidemia than BAV patients. Growth rate was 0.46 ±â€¯0.04 mm/year, similar in BAV and TAV groups (p = 0.70). Predictors of cardiac surgery were aorta dimensions at baseline (HR 1.23, p = 0.01), severe aortic regurgitation developed during follow-up (HR 3.49, p 0.04), family history of aortic aneurysm (HR 4.16, p 1.73), and history of STEMI (HR 3.64, p < 0.001). CONCLUSIONS: Classic baseline risk factors were more commonly observed in TAV aortopathy compared with BAV aortopathy. However, it is reassuring that, though diagnosed with aneurysm on average 10 years earlier and in the absence of arterial hypertension, BAV patients had a relatively low growth rate, similar to patients with a tricuspid valve. Irrespective of aortic valve morphology, patients with a family history of aortic aneurysm, history of coronary artery disease, and those who developed severe aortic regurgitation at follow-up, had the highest chances of being referred for surgery.

Computational fluid dynamics of the ascending aorta before the onset of type A aortic dissection.

We performed a pre-dissection computational fluid analysis of an ascending aortic aneurysm associated with unicuspid aortic valve. The analysis showed an abnormal helical flow pattern inside the aneurysm and an increased wall stress on the right postero-lateral wall of the ascending aorta. These values were largely higher than the theoretical cut-off for aortic wall dissection, their topographic distribution followed the intimal tear site as subsequently diagnosed by computed tomography scan and confirmed during the operation for dissection repair.

Early distal remodeling after elephant trunk repair of thoraco-abdominal aortic aneurysms.

Hemodynamic alterations occur when the elephant trunk (ET) technique is adopted to treat extensive aortic aneurysms. In planning the 2nd stage operation to complete ET repair, surgeons must weigh an adequate recovery time after initial surgery against the risk of postoperative ET-related complications. The purpose of this study was to understand the mechanistic link between the flow alteration caused by the ET graft and the development of premature aortic rupture before the 2nd stage operation. Specifically, fluid-structure interaction (FSI) analysis was performed using the CT imaging data of aorta at different stages of ET repair, and then computational variables were compared to those observed in patients who underwent a prophylactic 2nd stage operation to complete aortic repair. Results show that intramural stress exerted near the distal ET anastomosis (IMS=37.5kPa) at the time of urgent intervention was comparable to that of the extensive aortic aneurysm (IMS=47.4kPa) at initial in-hospital admission, but was considerably higher than that occurring after the 1st stage procedure (IMS=3.5kPa). Pressure index suggested higher peri-graft pressurization than aortic lumen pressure during diastole, imparting an apparent risk of aortic dilatation. These critical hemodynamic and structural parameters are related to the impending rupture of descending aorta observed clinically and can thus guide prophylactic intervention and optimal timing for the 2nd stage operation of a ET technique.

Constitutive modeling of ascending thoracic aortic aneurysms using microstructural parameters.

Ascending thoracic aortic aneurysm (ATAA) has been associated with diminished biomechanical strength and disruption in the collagen fiber microarchitecture. Additionally, the congenital bicuspid aortic valve (BAV) leads to a distinct extracellular matrix structure that may be related to ATAA development at an earlier age than degenerative aneurysms arising in patients with the morphological normal tricuspid aortic valve (TAV). The purpose of this study was to model the fiber-reinforced mechanical response of ATAA specimens from patients with either BAV or TAV. This was achieved by combining image-analysis derived parameters of collagen fiber dispersion and alignment with tensile testing data. Then, numerical simulations were performed to assess the role of anisotropic constitutive formulation on the wall stress distribution of aneurysmal aorta. Results indicate that both BAV ATAA and TAV ATAA have altered collagen fiber architecture in the medial plane of experimentally-dissected aortic tissues when compared to normal ascending aortic specimens. The study findings highlight that differences in the collagen fiber distribution mostly influences the resulting wall stress distribution rather than the peak stress. We conclude that fiber-reinforced constitutive modeling that takes into account the collagen fiber defect inherent to the aneurysmal ascending aorta is paramount for accurate finite element predictions and ultimately for biomechanical-based indicators to reliably distinguish the more from the less 'malignant' ATAAs.

Role of computational modeling in thoracic aortic pathology: a review.

Thoracic aortic diseases are life-threatening conditions causing significant mortality and morbidity despite advances in diagnostic and surgical treatments. Computational methods combined with imaging techniques provide quantitative information of disease progression, which may improve clinical treatments and therapeutic strategies for clinical practice. Since hemodynamic and wall mechanics play important roles in the natural history and progression of aortic diseases, we reviewed the potential application of computational modeling of the thoracic aorta. We placed emphasis on the clinical relevance of these techniques for the assessment of aortic dissection, thoracic aortic aneurysm, and aortic coarctation. Current clinical guidelines and treatment are also described.

Haemodynamic predictors of a penetrating atherosclerotic ulcer rupture using fluid-structure interaction analysis.

We present preliminary data on the flow-induced haemodynamic and structural loads exerted on a penetrating atherosclerotic aortic ulcer (PAU). Specifically, one-way fluid-structure interaction analysis was performed on the aortic model reconstructed from a 66-year-old male patient with a PAU that evolved into an intramural haematoma and rupture of the thoracic aorta. The results show that elevated blood pressure (117 mmHg) and low flow velocity at the aortic wall (0.15 m/s(2)) occurred in the region of the PAU. We also found a low value of time-averaged wall shear stress (1.24 N/m(2)) and a high value of the temporal oscillation in the wall shear stress (oscillatory shear index = 0.13) in the region of the PAU. After endovascular treatment, these haemodynamic parameters were distributed uniformly on the luminal surface of the stent graft. These findings suggest that wall shear stress could be considered one of the major haemodynamic factors indicating the structural fragility of the PAU wall, which ultimately lead to PAU growth and rupture.

Computer modeling for the prediction of thoracic aortic stent graft collapse.

OBJECTIVE: To assess the biomechanical implications of excessive stent protrusion into the aortic arch in relation to thoracic aortic stent graft (TASG) collapse by simulating the structural load and quantifying the fluid dynamics on the TASG wall protrusion extended into a model arch. METHODS: One-way coupled fluid-solid interaction analyses were performed to investigate the flow-induced hemodynamic and structural loads exerted on the proximal protrusion of the TASG and aortic wall reconstructed from a patient who underwent traumatic thoracic aortic injury repair. Mechanical properties of a Gore TAG thoracic endoprosthesis (W. L. Gore and Assoc, Flagstaff, Ariz) were assessed via experimental radial compression testing and incorporated into the computational modeling. The TASG wall protrusion geometry was characterized by the protrusion extension (PE) and by the angle (θ) between the TASG and the lesser curvature of the aorta. The effect of θ was explored with the following four models with PE fixed at 1.1 cm: θ = 10 degrees, 20 degrees, 30 degrees, and 40 degrees. The effect of PE was evaluated with the following four models with θ fixed at 10 degrees: PE = 1.1 cm, 1.4 cm, 1.7 cm and 2.0 cm. RESULTS: The presence of TASG wall protrusion into the aortic arch resulted in the formation of swirling, complex flow regions in the proximal luminal surface of the endograft. High PE values (PE = 2.0 cm) led to a markedly reduced left subclavian flow rate (0.27 L/min), low systolic perfusion pressure (98 mm Hg), and peak systolic TASG diameter reduction (2 mm). The transmural pressure load across the TASG was maximum for the model with the highest PE and θ, 15.2 mm Hg for the model with PE = 2.0 cm and θ = 10 degrees, and 11.6 mm Hg for PE = 1.1 cm and θ = 40 degrees. CONCLUSIONS: The findings of this study suggest that increased PE imparts an apparent risk of distal end-organ malperfusion and proximal hypertension and that both increased PE and θ lead to a markedly increased transmural pressure across the TASG wall, a load that would portend TASG collapse. Patient-specific computational modeling may allow for identification of patients with high risk of TASG collapse and guide preventive intervention. CLINICAL RELEVANCE: A potentially devastating complication that may occur after endovascular repair of traumatic thoracicaortic injuries is stent graft collapse. Although usually asymptomatic, stent graft collapse may be accompanied by adverse hemodynamic consequences. Numerous anatomic and device-related factors contribute to the development of collapse, but predictive factors have not yet been clearly defined. In the present study, we assessed the relevant hemodynamics and solid mechanics underlying stent graft collapse using a computational fluid-structure interaction framework of stent graft malapposition. Our findings suggest that both increased stent graft angle and extension into the aortic arch lead to a markedly increased transmural pressure across the stent graft wall, portending collapse. Patient-specific computational modeling may allow for identification of patients at high risk for collapse and aid in planning for an additional, prophylactic intervention to avert its occurrence.

Effect of aneurysm on the mechanical dissection properties of the human ascending thoracic aorta.

OBJECTIVES: The acute dissection of an ascending thoracic aortic aneurysm (ATAA) represents a devastating separation of elastic layers occurring when the hemodynamic loads on the diseased wall exceed the adhesive strength between layers. At present, the mechanics underlying aortic dissection are largely unclear, and the biomechanical delamination properties of the aneurysmal aorta are not defined. Individuals with bicuspid aortic valve (BAV) are particularly predisposed to ascending aortic aneurysm formation, with a marked risk of aortic dissection. The purpose of this study was to evaluate and compare the dissection properties of nonaneurysmal and aneurysmal human ascending thoracic aorta from patients with BAV morphology or normal tricuspid aortic valve (TAV) morphology using biomechanical delamination testing. METHODS: The influence on the delamination strength (S(d)) of the aorta associated with BAV was compared with that in patients with TAV. After complete delamination of ATAA tissue samples, tensile tests were performed on each delaminated half for comparison of their tensile strengths. RESULTS: The results showed that the aneurysmal aortas with BAV and TAV have lower S(d) than nonaneurysmal aortas and that ATAA with BAV has a lower S(d) than that with TAV. We have found a significant difference in S(d) between longitudinal and circumferential directions of the nondiseased aorta, suggesting anisotropic dissection properties. CONCLUSIONS: The tensile testing results suggest that the weaker intimal half of the aortic wall might fail before the outer adventitial half. Scanning electron microscope analyses suggest different failure modalities of dissection between the two morphologies, and the lower S(d) in ATAAs appears to be associated with a disorganized microstructure. BAV ATAAs have a lower S(d) than TAV ATAAs, suggesting a greater propensity for aortic dissection.