Adaptive outflow boundary conditions improve post-operative predictions after repair of peripheral pulmonary artery stenosis

A probabilistic neural twin for treatment planning in peripheral pulmonary artery stenosis

The substantial computational cost of high-fidelity models in numerical hemodynamics has, so far, relegated their use mainly to offline treatment planning. New breakthroughs in data-driven architectures and optimization techniques for fast surrogate modeling provide an exciting opportunity to overcome these limitations, enabling the use of such technology for time-critical decisions. We discuss an application to the repair of multiple stenosis in peripheral pulmonary artery disease through either transcatheter pulmonary artery rehabilitation or surgery, where it is of interest to achieve desired pressures and flows at specific locations in the pulmonary artery tree, while minimizing the risk for the patient. Since different degrees of success can be achieved in practice during treatment, we formulate the problem in probability, and solve it through a sample-based approach. We propose a new offline-online pipeline for probabilistic real-time treatment planning which combines offline assimilation of boundary conditions, model reduction, and training dataset generation with online estimation of marginal probabilities, possibly conditioned on the degree of augmentation observed in already repaired lesions. Moreover, we propose a new approach for the parametrization of arbitrarily shaped vascular repairs through iterative corrections of a zero-dimensional approximant. We demonstrate this pipeline for a diseased model of the pulmonary artery tree available through the Vascular Model Repository.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

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.

Computational approaches for mechanobiology in cardiovascular development and diseases

The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.

A computational growth and remodeling framework for adaptive and maladaptive pulmonary arterial hemodynamics

Hemodynamic loading is known to contribute to the development and progression of pulmonary arterial hypertension (PAH). This loading drives changes in mechanobiological stimuli that affect cellular phenotypes and lead to pulmonary vascular remodeling. Computational models have been used to simulate mechanobiological metrics of interest, such as wall shear stress, at single time points for PAH patients. However, there is a need for new approaches that simulate disease evolution to allow for prediction of long-term outcomes. In this work, we develop a framework that models the pulmonary arterial tree through adaptive and maladaptive responses to mechanical and biological perturbations. We coupled a constrained mixture theory-based growth and remodeling framework for the vessel wall with a morphometric tree representation of the pulmonary arterial vasculature. We show that non-uniform mechanical behavior is important to establish the homeostatic state of the pulmonary arterial tree, and that hemodynamic feedback is essential for simulating disease time courses. We also employed a series of maladaptive constitutive models, such as smooth muscle hyperproliferation and stiffening, to identify critical contributors to development of PAH phenotypes. Together, these simulations demonstrate an important step toward predicting changes in metrics of clinical interest for PAH patients and simulating potential treatment approaches.

A Computational Study of Blood Flow Dynamics in the Pulmonary Arteries

In this work we study the blood dynamics in the pulmonary arteries by means of a 3D-0D geometric multiscale approach, where a detailed 3D model for the pulmonary arteries is coupled with a lumped parameters (0D) model of the cardiovascular system. We propose to investigate three strategies for the numerical solution of the 3D-0D coupled problem: the Splitting-Explicit and Implicit algorithms, where information are exchanged between 3D and 0D models at each time step at the interfaces, and the One-Way algorithm, where the 0D is solved first off-line. In our numerical experiments performed in a realistic patient-specific 3D domain with a physiologically calibrated 0D model, we discuss first the issue on instabilities that may arise when not suitable connections are considered between 3D and 0D models; second we compare the performance and accuracy of the three proposed numerical strategies. Finally, we report a comparison between a healthy and a hypertensive case, providing a preliminary result highlighting how our method could be used in future for clinical purposes.

Virtual Transcatheter Interventions for Peripheral Pulmonary Artery Stenosis in Williams and Alagille Syndromes

Background

Despite favorable outcomes of surgical pulmonary artery (PA) reconstruction, isolated proximal stenting of the central PAs is common clinical practice for patients with peripheral PA stenosis in association with Williams and Alagille syndromes. Given the technical challenges of PA reconstruction and the morbidities associated with transcatheter interventions, the hemodynamic consequences of all treatment strategies must be rigorously assessed. Our study aims to model, assess, and predict hemodynamic outcomes of transcatheter interventions in these patients.

Methods and Results

Isolated proximal and “extensive” interventions (stenting and/or balloon angioplasty of proximal and lobar vessels) were performed in silico on 6 patient‐specific PA models. Autoregulatory adaptation of the cardiac output and downstream arterial resistance was modeled in response to intervention‐induced hemodynamic perturbations. Postintervention computational fluid dynamics predictions were validated in 2 stented patients and quantitatively assessed in 4 surgical patients. Our computational methods accurately predicted postinterventional PA pressures, the primary indicators of success for treatment of peripheral PA stenosis. Proximal and extensive treatment achieved median reductions of 14% and 40% in main PA systolic pressure, 27% and 56% in pulmonary vascular resistance, and 10% and 45% in right ventricular stroke work, respectively.

Conclusions

In patients with Williams and Alagille syndromes, extensive transcatheter intervention is required to sufficiently reduce PA pressures and right ventricular stroke work. Transcatheter therapy was shown to be ineffective for long‐segment stenosis and pales hemodynamically in comparison with published outcomes of surgical reconstruction. Regardless of the chosen strategy, a virtual treatment planning platform could identify lesions most critical for optimizing right ventricular afterload.

Multiscale modelling of Potts shunt as a potential palliative treatment for suprasystemic idiopathic pulmonary artery hypertension: a paediatric case study

Potts shunt (PS) was suggested as palliation for patients with suprasystemic pulmonary arterial hypertension (PAH) and right ventricular (RV) failure. PS, however, can result in poorly understood mortality. Here, a patient-specific geometrical multiscale model of PAH physiology and PS is developed for a paediatric PAH patient with stent-based PS. In the model, 7.6mm-diameter PS produces near-equalisation of the aortic and PA pressures and [Formula: see text] (oxygenated vs deoxygenated blood flow) ratio of 0.72 associated with a 16% decrease of left ventricular (LV) output and 18% increase of RV output. The flow from LV to aortic arch branches increases by 16%, while LV contribution to the lower body flow decreases by 29%. Total flow in the descending aorta (DAo) increases by 18% due to RV contribution through the PS with flow into the distal PA branches decreasing. PS induces 18% increase of RV work due to its larger stroke volume pumped against lower afterload. Nonetheless, larger RV work does not lead to increased RV end-diastolic volume. Three-dimensional flow assessment demonstrates the PS jet impinging with a high velocity and wall shear stress on the opposite DAo wall with the most of the shunt flow being diverted to the DAo. Increasing the PS diameter from 5mm up to 10mm results in a nearly linear increase in post-operative shunt flow and a nearly linear decrease in shunt pressure-drop. In conclusion, this model reasonably represents patient-specific haemodynamics pre- and post-creation of the PS, providing insights into physiology of this complex condition, and presents a predictive tool that could be useful for clinical decision-making regarding suitability for PS in PAH patients with drug-resistant suprasystemic PAH.

A Personalized Pulmonary Circulation Model to Non-Invasively Calculate Fractional Flow Reserve for Artery Stenosis Detection

OBJECTIVE

Fractional Flow Reserve (FFR) is regarded as a fundamental index to assess pulmonary artery stenosis. The application of FFR can increase the accuracy of detection of pulmonary artery stenosis. However, the invasive examination may carry a number of physiological risks for patients. Therefore, we propose a personalized pulmonary circulation model to non- invasively calculate FFR of pulmonary artery stenosis. Method- ology: We employed a personalized pulmonary circulation model to non-invasively calculate FFR using only computed tomography angiogram (CTA) data. This model combined boundary conditions estimation and 3D pulmonary artery morphology reconstruction for CFD simulation. First, we obtained patient-specific boundary conditions by adapting the right ventricle stroke volume and main pulmonary artery pressure feature points (systolic, diastolic, and mean pressure). Secondly, the 3D pulmonary artery morphology was reconstructed by threshold segmentation. The CFD simulation was then performed to obtain pressure distribution in the entire pulmonary artery. Finally, the FFR in pulmonary artery stenoses was calculated as the ratio of distal pressure and proximal pres- sure.

RESULTS

To validate our model, we compared the calculated FFR with measured FFR by pressure guide wires examination of 8 patients. The FFR calculated by our model showed a good agreement with measured FFR by pressure guide wires exami- nation. The average accuracy rate was 91.41%.

CONCLUSION

The proposed personalized pulmonary model is capable of reasonably non-invasively calculating FFR with sufficient accuracy.

SIGNIFICANCE

FFR calculated in our model may contribute to non-invasive detection of pulmonary artery stenosis and to the assessment of invasive interventions.

Non-invasive MRI Derived Hemodynamic Simulation to Predict Successful vs. Unsuccessful Catheter Interventions for Branch Pulmonary Artery Stenosis: Proof-of-Concept and Experimental Validation in Swine

OBJECTIVE

This study assessed the ability of hemodynamic simulations to predict the success of catheter interventions in a swine model of branch pulmonary artery stenosis (bPAS).

BACKGROUND

bPAS commonly occurs in congenital heart disease and is often managed with catheter based interventions. However, despite technical success, bPAS interventions do not lead to improved distal pulmonary blood flow (PBF) distribution in approximately 1/3rd of patients. New tools are needed to better identify which patients with bPAS would most benefit from catheter interventions.

METHODS

For 13 catheter intervention cases in swine with surgically created left PAS (LPAS), PA pressures from right heart catheterization (RHC) and PBF distributions from MRI were measured before and after catheter interventions. Hemodynamic simulations with a reduced order computational fluid dynamics (CFD) model were performed using non-invasive PBF measurements derived from MRI, and then correlated with changes in invasive measures of hemodynamics and PBF distributions before and after catheter intervention to relieve bPAS.

RESULTS

Compared to experimentally measured changes in left PBF distribution, simulations had a small bias (3.4 ± 11.1%), moderate agreement (ICC = 0.69 [0.24-0.90], 0.71 [0.23-0.91]), and good diagnostic capability to predict successful interventions (> 20% PBF increase) (AUC 0.83 [0.59-1.0]). Simulations had poorer prediction of changes in stenotic pressure gradient (ICC = 0.28 [- 0.33 to 0.73], r = 0.57 [- 0.04 to 0.87]) and MPA systolic pressure (ICC = 0.00 [- 0.52 to 0.53], r = 0.29 [- 0.32 to 0.72]).

CONCLUSION

While there was only weak to moderate agreement between predicted and measured changes in PA pressures and pulmonary blood flow distributions, hemodynamic simulations did show good diagnostic value for predicting successful versus unsuccessful catheter based interventions to relieve bPAS. The results of this proof of concept study are promising and should encourage future development for using hemodynamic models in planning interventions for patients with bPAS.

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.

Ex Vivo Models to Decipher the Molecular Mechanisms of Genetic Notch Cardiovascular Disorders

Notch is an evolutionary, conserved, cell-cell signaling pathway that is central to several biological processes, from tissue morphogenesis to homeostasis. It is therefore not surprising that several genetic mutations of Notch components cause inherited human diseases, especially cardiovascular disorders. Despite numerous efforts, current in vivo models are still insufficient to unravel the underlying mechanisms of these pathologies, hindering the development of utmost needed medical therapies. In this perspective review, we discuss the limitations of current murine models and outline how the combination of microphysiological systems (MPSs) and targeted computational models can lead to breakthroughs in this field. In particular, while MPSs enable the experimentation on human cells in controlled and physiological environments, in silico models can provide a versatile tool to translate the in vitro findings to the more complex in vivo setting. As a showcase example, we focus on Notch-related cardiovascular diseases, such as Alagille syndrome, Adams-Oliver syndrome, and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). Impact statement In this review, a comprehensive overview of the limitations of current in vivo models of genetic Notch cardiovascular diseases is provided, followed by a discussion over the potential of microphysiological systems and computational models in overcoming these limitations and in potentiating drug testing and modeling of these pathologies.

Longitudinal Evolution of Pulmonary Artery Wall Shear Stress in a Swine Model of Pulmonary Artery Stenosis and Stent Interventions

Branch pulmonary artery stenosis (PAS) commonly occurs in congenital heart disease and it has previously been hypothesized that in branch PAS the pulmonary arteries (PAs) remodel their lumen diameter to maintain constant wall shear stress (WSS). We quantified the longitudinal progression of PA WSS in a swine model of unilateral PAS and two different intervention time courses to test this hypothesis. To quantify WSS in the entire pulmonary tree we used 4D Flow MRI for the large-proximal PAs and a structured tree model for the small-distal PAs. Our results only partially supported the hypothesis that in branch PAS the PAs remodel their lumen diameter to maintain WSS homeostasis. Proximal PA WSS was similar between groups at the final study time-point but WSS of mid-sized (5 mm to 500 μm) PA segments was found to be different between the sham and LPAS groups. This suggests that WSS homeostasis may only be achieved for the large-proximal PAs. Additionally, our results do not show WSS homeostasis being achieved over shorter periods of time suggesting that any potential WSS dependent changes in PA lumen diameter were a long-term remodeling response rather than a short-term vasodilation response. Future studies should confirm if these findings hold true in humans and investigate the impacts of WSS at different levels of the pulmonary tree on growth.

Image-based scaling laws for somatic growth and pulmonary artery morphometry from infancy to adulthood

Pulmonary artery (PA) morphometry has been extensively explored in adults, with particular focus on intra-acinar arteries. However, scaling law relationships for length and diameter of extensive preacinar PAs by age have not been previously reported for in vivo human data. To understand preacinar PA growth spanning children to adults, we performed morphometric analyses of all PAs visible in the computed tomography (CT) and magnetic resonance (MR) images from a healthy subject cohort [n = 16; age: 1-51 yr; body surface area (BSA): 0.49-2.01 m²]. Subject-specific anatomic PA models were constructed from CT and MR images, and morphometric information-diameter, length, tortuosity, bifurcation angle, and connectivity-was extracted and sorted into diameter-defined Strahler orders. Validation of Murray's law, describing optimal scaling exponents of radii for branching vessels, was performed to determine how closely PAs conform to this classical relationship. Using regression analyses of vessel diameters and lengths against orders and patient metrics (BSA, age, height), we found that diameters increased exponentially with order and allometrically with patient metrics. Length increased allometrically with patient metrics, albeit weakly. The average tortuosity index of all vessels was 0.026 ± 0.024, average bifurcation angle was 28.2 ± 15.1°, and average Murray's law exponent was 2.92 ± 1.07. We report a set of scaling laws for vessel diameter and length, along with other morphometric information. These provide an initial understanding of healthy structural preacinar PA development with age, which can be used for computational modeling studies and comparison with diseased PA anatomy.NEW & NOTEWORTHY Pulmonary artery (PA) morphometry studies to date have focused primarily on large arteries and intra-acinar arteries in either adults or children, neglecting preacinar arteries in both populations. Our study is the first to quantify in vivo preacinar PA morphometry changes spanning infants to adults. For preacinar arteries > 1 mm in diameter, we identify scaling laws for vessel diameters and lengths with patient metrics of growth and establish a healthy PA morphometry baseline for most preacinar PAs.

Computed Tomography Pulmonary Vascular Volume Ratio Can Be Used to Evaluate the Effectiveness of Pulmonary Angioplasty in Peripheral Pulmonary Artery Stenosis

OBJECTIVE

To explore whether computed tomography (CT) pulmonary vascular volume ratio can be used to evaluate the effectiveness of pulmonary artery angioplasty in patients with peripheral pulmonary artery stenosis.

MATERIALS AND METHODS

Changes in CT pulmonary vascular volume ratio between serial cardiothoracic CT examinations were calculated in 38 patients. Fifteen patients underwent interim pulmonary artery angioplasty (group 1), while 23 did not (group 2). According to the effectiveness of pulmonary artery angioplasty, patients in group 1 were further divided into group 1A (improved or aggravated) and group 1B (ineffective). Changes in the pulmonary vascular volume percentages among the three groups (group 1A, group 1B, and group 2) on serial CT examinations were compared.

RESULTS

Pulmonary artery angioplasty on serial CT examinations was successful in seven patients, ineffective in seven patients, and aggravated in one patient. As a result, eight patients were included in group 1A and seven were included in group 1B. Changes in the CT pulmonary vascular volume percentages in group 1A were statistically significantly greater than those in group 1B (11.6 ± 5.6% vs. 2.7 ± 1.6%, p < 0.003) and group 2 (11.6 ± 5.6% vs. 1.9 ± 1.4%, p < 0.002), while no statistically significant difference was found between group 1B and group 2 (2.7 ± 1.6% vs. 1.9 ± 1.4%, p > 0.1).

CONCLUSION

CT pulmonary vascular volume ratio can be used to evaluate the effectiveness of pulmonary artery angioplasty in patients with peripheral pulmonary artery stenosis.

4D Flow MRI Estimation of Boundary Conditions for Patient Specific Cardiovascular Simulation

Accurate image based cardiovascular simulations require patient specific boundary conditions (BCs) for inlets, outlets and vessel wall mechanical properties. While inlet BCs are typically determined non-invasively, invasive pressure catheterization is often used to determine patient specific outlet BCs and vessel wall mechanical properties. A method using 4D Flow MRI to non-invasively determine both patient specific outlet BCs and vessel wall mechanical properties is presented and results for both in vitro validation with a latex tube and an in vivo pulmonary artery stenosis (PAS) stent intervention are presented. For in vitro validation, acceptable agreement is found between simulation using BCs from 4D Flow MRI and benchtop measurements. For the PAS virtual intervention, simulation correctly predicts flow distribution with 9% error compared to MRI. Using 4D Flow MRI to noninvasively determine patient specific BCs increases the ability to use image based simulations as pressure catheterization is not always performed.

Computational Study of Pulmonary Flow Patterns After Repair of Transposition of Great Arteries

Patients that undergo the arterial switch operation (ASO) to repair transposition of great arteries (TGA) can develop abnormal pulmonary trunk morphology with significant long-term complications. In this study, cardiovascular magnetic resonance was combined with computational fluid dynamics to investigate the impact of the postoperative layout on the pulmonary flow patterns. Three ASO patients were analyzed and compared to a volunteer control. Results showed the presence of anomalous shear layer instabilities, vortical and helical structures, and turbulent-like states in all patients, particularly as a consequence of the unnatural curvature of the pulmonary bifurcation. Streamlined, mostly laminar flow was instead found in the healthy subject. These findings shed light on the correlation between the post-ASO anatomy and the presence of altered flow features, and may be useful to improve surgical planning as well as the long-term care of TGA patients.

Computational simulation of postoperative pulmonary flow distribution in Alagille patients with peripheral pulmonary artery stenosis

BACKGROUND

Up to 90% of individuals with Alagille syndrome have congenital heart diseases. Peripheral pulmonary artery stenosis (PPS), resulting in right ventricular hypertension and pulmonary flow disparity, is one of the most common abnormalities, yet the hemodynamic effects are ill-defined, and optimal patient management and treatment strategies are not well established. The purpose of this pilot study is to use recently refined computational simulation in the setting of multiple surgical strategies, to examine the influence of pulmonary artery reconstruction on hemodynamics in this population.

MATERIALS AND METHODS

Based on computed tomography angiography and cardiac catheterization data, preoperative pulmonary artery models were constructed for 4 patients with Alagille syndrome with PPS (all male, age range: 0.6-2.9 years), and flow simulations with deformable walls were performed. Surgeon directed virtual surgery, mimicking the surgical procedure, was then performed to derive postoperative models. Postoperative simulation-derived hemodynamics and blood flow distribution were then compared with the clinical results.

RESULTS

Simulations confirmed substantial resistance, resulting from preoperative severe ostial stenoses, and the use of newly developed adaptive outflow boundary conditions led to excellent agreement with in vivo measurements. Relief of PPS decreased pulmonary artery pressures and improved pulmonary flow distribution both in vivo and in silico with good correlation.

CONCLUSIONS

Using adaptive outflow boundary conditions, computational simulations can estimate postoperative overall pulmonary flow distribution in patients with Alagille syndrome after pulmonary artery reconstruction. Obstruction relief along with pulmonary artery vasodilation determines postoperative pulmonary flow distribution and newer methods can incorporate these physiologic changes. Evolving blood flow simulations may be useful in surgical or transcatheter planning and in understanding the complex interplay among various obstructions in patients with peripheral pulmonary stenosis.

Cardiovascular mechanics in the early stages of pulmonary hypertension: a computational study

We formulate and study a new mathematical model of pulmonary hypertension. Based on principles of fluid and elastic dynamics, we introduce a model that quantifies the stiffening of pulmonary vasculature (arteries and arterioles) to reproduce the hemodynamics of the pulmonary system, including physiologically consistent dependence between compliance and resistance. This pulmonary model is embedded in a closed-loop network of the major vessels in the body, approximated as one-dimensional elastic tubes, and zero-dimensional models for the heart and other organs. Increasingly severe pulmonary hypertension is modeled in the context of two extreme scenarios: (1) no cardiac compensation and (2) compensation to achieve constant cardiac output. Simulations from the computational model are used to estimate cardiac workload, as well as pressure and flow traces at several locations. We also quantify the sensitivity of several diagnostic indicators to the progression of pulmonary arterial stiffening. Simulation results indicate that pulmonary pulse pressure, pulmonary vascular compliance, pulmonary RC time, luminal distensibility of the pulmonary artery, and pulmonary vascular impedance are much better suited to detect the early stages of pulmonary hypertension than mean pulmonary arterial pressure and pulmonary vascular resistance, which are conventionally employed as diagnostic indicators for this disease.