Blood flow simulations in the pulmonary bifurcation in relation to adult patients with repaired tetralogy of Fallot

Advancement in computational simulation and validation of congenital heart disease: a review

The improvement in congenital heart disease (CHD) treatment and management has increased the life expectancy in infants. However, the long-term efficacy is difficult to assess and thus, computational modelling has been applied for evaluating this. Here, we provide an overview of the applications of computational modelling in CHD based on three categories; CHD involving large blood vessels only, heart chambers only, and CHD that occurs at multiple heart structures. We highlight the advancement of computational simulation of CHD that uses multiscale and multiphysics modelling to ensure a complete representation of the heart and circulation. We provide a brief future direction of computational modelling of CHD such as to include growth and remodelling, detailed conduction system, and occurrence of myocardial infarction. We also proposed validation technique using advanced three-dimensional (3D) printing and particle image velocimetry (PIV) technologies to improve the model accuracy.

Hemodynamic-based virtual surgery design of double-patch repair for pulmonary arterioplasty in tetralogy of Fallot

BACKGROUND AND OBJECTIVE

Surgical correction of pulmonary artery stenosis (PAS) is essential to the prognosis of patients with tetralogy of Fallot (TOF). The double-patch method of pulmonary arterioplasty is usually applied in case of multiple stenosis in TOF patients' pulmonary artery (PA) and when PAS cannot be relieved by the single-patch method. The surgical planning for the double-patch design remains challenging. The purpose of this study is to investigate the double-patch design with different angulations between the left pulmonary artery (LPA) and the right pulmonary artery (RPA), and to understand postoperative hemodynamic alterations by the application of computer-aided design (CAD) and computational fluid dynamics (CFD) techniques.

METHODS

The three-dimensional model of the PA was reconstructed based on preoperative computed tomography imaging data obtained from the patient with TOF. Three postoperative models with different designs of double-patch were created by "virtual surgery" using the CAD technique. Double-Patch 120 Model was created with double patches implanted in the main pulmonary artery (MPA) and the PA bifurcation and without changing the spatial position of PA. The angulation between the LPA and the RPA was defined as θ, which equaled to 120° in Pre-Operative Model and Double-Patch 120 Model. Based on Double-Patch 120 Model, Double-Patch 110 Model and Double-Patch 130 Model were generated with θ equaled to 110° and 130°, respectively. Combined with CFD, the differences of velocity streamlines, wall shear stress (WSS), flow distribution ratio (FDR), and energy loss (EL) were compared to analyze postoperative pulmonary flow characteristics.

RESULTS

The values of velocity and WSS decreased significantly after virtual surgery. Obvious vortices and swirling flows were observed downstream of the stenosis of RPA and LPA in Pre-Operative Model, while fewer vortices developed along the anterior wall of the expanded lumens of RPA, especially in Double-Patch 110 Model. With the relief of PAS, two relatively higher WSS regions were observed at the posterior walls of RPA and LPA. The maximum WSS values in these regions of Double-Patch 110 Model were lower than those in Double-Patch 120 Model and Double-Patch 130 Model. Furthermore, the FDRs were elevated and the ELs were greatly reduced. It was found that Double-Patch 110 Model with the angulation between the LPA and the RPA equaled to 110° showed relatively better properties of hemodynamics than other models.

CONCLUSIONS

The angulation between the LPA and the RPA is an important factor that should be integrated in the double-patch design for TOF repair. Virtual surgery based on patient-specific vascular model and computational hemodynamics can be used to provide assistance for individualized surgical planning of double-patch arterioplasty.

Percutaneous approach to residual pulmonary bifurcation stenosis in conotruncal diseases

Residual stenosis after right ventricle outflow tract surgery represents a major issue to manage in the children and adult patient with conotruncal defects. Despite a detailed multimodality imaging, the anatomy of distal pulmonary trunk and pulmonary artery bifurcation may be challenging in these patients.The aim of this study was to analyse retrospectively the outcome of the percutaneous transcatheter treatment in children with post-surgical stenosis of pulmonary artery bifurcation.We enrolled 39 patients with a median age of 6.0 years. Standard high-pressure balloon dilation was attempted in 33 patients, effective in 5 of them. Pulmonary branch stenting was performed in 10 patients, effective in 6. A kissing balloon approach was chosen in 17 patients (6 after angioplasty or stenting failure), and this technique was effective in 16 cases. Finally, a bifurcation stenting was performed in 10 patients (second step in 9 cases), effective in all the cases. None of the patients approached by kissing balloon needed a bifurcation stenting.In conclusion, standard balloon angioplasty and standard stenting might be ineffective in post-surgical stenosis involving pulmonary artery bifurcation. In this population, kissing balloon or bifurcation stenting, followed by side branch de-jailing, may be more effective in relieving the gradient.

Morphological characterisation of pediatric Turner syndrome aortae: Insights from a small cohort study

Cardiovascular disease is widespread in girls and women living with Turner syndrome (TS). Despite this prevalence, cardiovascular risk evaluation using the current guidelines has seen life-threatening aortic events occurring at dimensions classified within the normal threshold. In this study, we characterized the three-dimensional aortic geometries of Turner syndrome children and their age-matched healthy counterparts to evaluate various morphological parameters. Turner syndrome girls had overall greater values in ten out of fifteen parameters examined (p > 0.05), when compared to healthy children: the aortic arch height and width; the ascending aorta, aortic arch (2 locations), and descending aorta diameters; the ratio of the ascending to descending aorta diameter; average curvature; average torsion; and average curvature-torsion score. Additionally, significant associations were found in the TS group: body surface area and both arch height (p = 0.03) and arch height to width ratio (p = 0.05), and aortic arch diameter and both body surface area (p = 0.04) and weight (p = 0.04). The new information resulting from this small cohort study contributes to an improved understanding of the morphological parameters affecting the hemodynamic environment in TS, and the clinical assessment of the increased cardiovascular risk in this population.

Reconstruction and Validation of Arterial Geometries for Computational Fluid Dynamics Using Multiple Temporal Frames of 4D Flow-MRI Magnitude Images

PURPOSE

Segmentation and reconstruction of arterial blood vessels is a fundamental step in the translation of computational fluid dynamics (CFD) to the clinical practice. Four-dimensional flow magnetic resonance imaging (4D Flow-MRI) can provide detailed information of blood flow but processing this information to elucidate the underlying anatomical structures is challenging. In this study, we present a novel approach to create high-contrast anatomical images from retrospective 4D Flow-MRI data.

METHODS

For healthy and clinical cases, the 3D instantaneous velocities at multiple cardiac time steps were superimposed directly onto the 4D Flow-MRI magnitude images and combined into a single composite frame. This new Composite Phase-Contrast Magnetic Resonance Angiogram (CPC-MRA) resulted in enhanced and uniform contrast within the lumen. These images were subsequently segmented and reconstructed to generate 3D arterial models for CFD. Using the time-dependent, 3D incompressible Reynolds-averaged Navier-Stokes equations, the transient aortic haemodynamics was computed within a rigid wall model of patient geometries.

RESULTS

Validation of these models against the gold standard CT-based approach showed no statistically significant inter-modality difference regarding vessel radius or curvature (p > 0.05), and a similar Dice Similarity Coefficient and Hausdorff Distance. CFD-derived near-wall hemodynamics indicated a significant inter-modality difference (p > 0.05), though these absolute errors were small. When compared to the in vivo data, CFD-derived velocities were qualitatively similar.

CONCLUSION

This proof-of-concept study demonstrated that functional 4D Flow-MRI information can be utilized to retrospectively generate anatomical information for CFD models in the absence of standard imaging datasets and intravenous contrast.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

Impact of atrial fibrillation on left atrium haemodynamics: A computational fluid dynamics study

We analyse the haemodynamics of the left atrium, highlighting differences between healthy individuals and patients affected by atrial fibrillation. The computational study is based on patient-specific geometries of the left atria to simulate blood flow dynamics. We design a novel procedure to compute the boundary data for the 3D haemodynamic simulations, which are particularly useful in absence of data from clinical measurements. With this aim, we introduce a parametric definition of atrial displacement, and we use a closed-loop lumped parameter model of the whole cardiovascular circulation conveniently tuned on the basis of the patient's characteristics. We evaluate several fluid dynamics indicators for atrial haemodynamics, validating our numerical results in terms of clinical measurements; we investigate the impact of geometric and clinical characteristics on the risk of thrombosis. To highlight the correlation of thrombus formation with atrial fibrillation, according to medical evidence, we propose a novel indicator: age stasis. It arises from the combination of Eulerian and Lagrangian quantities. This indicator identifies regions where slow flow cannot properly rinse the chamber, accumulating stale blood particles, and creating optimal conditions for clots formation.

Numerical analysis of hemodynamics in pulmonary artery stenosis

BACKGROUND: Pulmonary artery stenosis is a serious threat to people’s life and health. OBJECTIVE: The hydrodynamic mechanism of pulmonary artery stenosis is investigated. METHODS: Numerical analysis of hemodynamics in pulmonary artery stenosis using computational fluid dynamics techniques is a well-established method. An idealized model of pulmonary stenosis is established, and the model is divided into main pulmonary artery, right and left pulmonary arteries, and their branches. The sections at different positions are intercepted to study the distribution trend of maximum velocity, pressure and wall shear stress. RESULTS: The numerical simulation results show that the pressure drop at both ends of the narrow is large. High velocity and wall shear stress exist in the center of stenosis, and the wall shear stress at the distal end of stenosis gradually decreases, resulting in endothelial dysfunction. CONCLUSIONS:  To some extent, this study helps clinicians make diagnosis and treatment plans in advance and improve prognosis. This method could be used in the numerical simulation of practical models.

Elastin Insufficiency Confers Proximal and Distal Pulmonary Vasculopathy in Mice, Partially Remedied by the KATP Channel Opener Minoxidil: Considerations and Cautions for the Treatment of People With Williams-Beuren Syndrome

Background

Williams Beuren syndrome (WBS) is a recurrent microdeletion disorder that removes one copy of elastin (ELN), resulting in large artery vasculopathy. Early stenosis of the pulmonary vascular tree is common, but few data are available on longer-term implications of the condition.

Methods

Computed tomography (CT) angiogram (n = 11) and echocardiogram (n = 20) were performed in children with WBS aged 3.4–17.8 years. Controls (n = 11, aged 4.4–16.8 years) also underwent echocardiogram. Eln^+/− mice were analyzed by invasive catheter, echocardiogram, micro-CT (μCT), histology, and pressure myography. We subsequently tested whether minoxidil resulted in improved pulmonary vascular endpoints.

Results

WBS participants with a history of main or branch pulmonary artery (PA) stenosis requiring intervention continued to exhibit increased right ventricular systolic pressure (RVSP, echocardiogram) relative to their peers without intervention (p < 0.01), with no clear difference in PA size. Untreated Eln^+/− mice also show elevated RVSP by invasive catheterization (p < 0.0001), increased normalized right heart mass (p < 0.01) and reduced caliber branch PAs by pressure myography (p < 0.0001). Eln^+/− main PA medias are thickened histologically relative to Eln^+/+ (p < 0.0001). Most Eln^+/− phenotypes are shared by both sexes, but PA medial thickness is substantially greater in Eln^+/− males (p < 0.001). Eln^+/− mice showed more acute proximal branching angles (p < 0.0001) and longer vascular segment lengths (p < 0.0001) (μCT), with genotype differences emerging by P7. Diminished PA acceleration time (p < 0.001) and systolic notching (p < 0.0001) were also observed in Eln^+/− echocardiography. Vascular casting plus μCT revealed longer generation-specific PA arcade length (p < 0.0001), with increased PA branching detectable by P90 (p < 0.0001). Post-weaning minoxidil decreased RVSP (p < 0.01) and normalized PA caliber (p < 0.0001) but not early-onset proximal branching angle or segment length, nor later-developing peripheral branch number.

Conclusions

Vascular deficiencies beyond arterial caliber persist in individuals with WBS who have undergone PA stenosis intervention. Evaluation of Eln^+/− mice reveals complex vascular changes that affect the proximal and distal vasculatures. Minoxidil, given post-weaning, decreases RVSP and improves lumen diameter, but does not alter other earlier-onset vascular patterns. Our data suggest additional therapies including minoxidil could be a useful adjunct to surgical therapy, and future trials should be considered.

Characterization of Flow Dynamics in the Pulmonary Bifurcation of Patients With Repaired Tetralogy of Fallot: A Computational Approach

The hemodynamic environment of the pulmonary bifurcation is of great importance for adult patients with repaired tetralogy of Fallot (rTOF) due to possible complications in the pulmonary valve and narrowing of the left pulmonary artery (LPA). The aim of this study was to computationally investigate the effect of geometrical variability and flow split on blood flow characteristics in the pulmonary trunk of patient-specific models. Data from a cohort of seven patients was used retrospectively and the pulmonary hemodynamics was investigated using averaged and MRI-derived patient-specific boundary conditions on the individualized models, as well as a statistical mean geometry. Geometrical analysis showed that curvature and tortuosity are higher in the LPA branch, compared to the right pulmonary artery (RPA), resulting in complex flow patterns in the LPA. The computational analysis also demonstrated high time-averaged wall shear stress (TAWSS) at the outer wall of the LPA and the wall of the RPA proximal to the junction. Similar TAWSS patterns were observed for averaged boundary conditions, except for a significantly modified flow split assigned at the outlets. Overall, this study enhances our understanding about the flow development in the pulmonary bifurcation of rTOF patients and associates some morphological characteristics with hemodynamic parameters, highlighting the importance of patient-specificity in the models. To confirm these findings, further studies are required with a bigger cohort of patients.

Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls

Congenital abnormalities in girls and women with Turner syndrome (TS), alongside an underlying predisposition to obesity and hypertension, contribute to an increased risk of cardiovascular disease and ultimately reduced life expectancy. We observe that children with TS present a greater variance in aortic arch morphology than their healthy counterparts, and hypothesize that their hemodynamics is also different. In this study, computational fluid dynamic (CFD) simulations were performed for four TS girls, and three age-matched healthy girls, using patient-specific inlet boundary conditions, obtained from phase-contrast MRI data. The visualization of multidirectional blood flow revealed an increase in vortical flow in the arch, supra-aortic vessels, and descending aorta, and a correlation between the presence of aortic abnormalities and disturbed flow. Compared to the relatively homogeneous pattern of time-averaged wall shear stress (TAWSS) on the healthy aortae, a highly heterogeneous distribution with elevated TAWSS values was observed in the TS geometries. Visualization of further shear stress parameters, such as oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS), revealed dissimilar heterogeneity in the oscillatory and multidirectional nature of the aortic flow. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% were obese or overweight), which is believed to accelerate the initiation and progression of endothelial dysfunction, these findings may be an indication of atherosclerotic disease manifesting earlier in life in TS patients. Age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk in TS children.

Computational Analysis of the Pulmonary Arteries in Congenital Heart Disease: A Review of the Methods and Results

With the help of computational fluid dynamics (CFD), hemodynamics of the pulmonary arteries (PA's) can be studied in detail and varying physiological circumstances and treatment options can be simulated. This offers the opportunity to improve the diagnostics and treatment of PA stenosis in biventricular congenital heart disease (CHD). The aim of this review was to evaluate the methods of computational studies for PA's in biventricular CHD and the level of validation of the numerical outcomes. A total of 34 original research papers were selected. The literature showed a great variety in the used methods for (re) construction of the geometry as well as definition of the boundary conditions and numerical setup. There were 10 different methods identified to define inlet boundary conditions and 17 for outlet boundary conditions. A total of nine papers verified their CFD outcomes by comparing results to clinical data or by an experimental mock loop. The diversity in used methods and the low level of validation of the outcomes result in uncertainties regarding the reliability of numerical studies. This limits the current clinical utility of CFD for the study of PA flow in CHD. Standardization and validation of the methods are therefore recommended.