Influence of image segmentation on one-dimensional fluid dynamics predictions in the mouse pulmonary arteries

Dissecting contributions of pulmonary arterial remodeling to right ventricular afterload in pulmonary hypertension

Pulmonary hypertension (PH) is defined as an elevation in the right ventricle (RV) afterload, characterized by increased hemodynamic pressure in the main pulmonary artery (PA). Elevations in RV afterload increase RV wall stress, resulting in RV remodeling and potentially RV failure. From a biomechanical standpoint, the primary drivers for RV afterload elevations include increases in pulmonary vascular resistance (PVR) in the distal vasculature and decreases in vessel compliance in the proximal PA. However, the individual contributions of the various vascular remodeling events toward the progression of PA pressure elevations and altered vascular hemodynamics remain elusive. In this study, we used a subject-specific one-dimensional (1D) fluid-structure interaction (FSI) model to investigate the alteration of pulmonary hemodynamics in PH and to quantify the contributions of vascular stiffening and increased resistance towards increased main pulmonary artery (MPA) pressure. We used a combination of subject-specific hemodynamic measurements, ex-vivo mechanical testing of arterial tissue specimens, and ex-vivo X-ray micro-tomography imaging to develop the 1D-FSI model and dissect the contribution of PA remodeling events towards alterations in the MPA pressure waveform. Both the amplitude and pulsatility of the MPA pressure waveform were analyzed. Our results indicated that increased distal resistance has the greatest effect on the increase in maximum MPA pressure, while increased stiffness caused significant elevations in the characteristic impedance. The method presented in this study will serve as an essential step toward understanding the complex interplay between PA remodeling events that leads to the most severe adverse effect on RV dysfunction.

Computational framework for the generation of one-dimensional vascular models accounting for uncertainty in networks extracted from medical images

One-dimensional (1D) cardiovascular models offer a non-invasive method to answer medical questions, including predictions of wave-reflection, shear stress, functional flow reserve, vascular resistance and compliance. This model type can predict patient-specific outcomes by solving 1D fluid dynamics equations in geometric networks extracted from medical images. However, the inherent uncertainty in in vivo imaging introduces variability in network size and vessel dimensions, affecting haemodynamic predictions. Understanding the influence of variation in image-derived properties is essential to assess the fidelity of model predictions. Numerous programs exist to render three-dimensional surfaces and construct vessel centrelines. Still, there is no exact way to generate vascular trees from the centrelines while accounting for uncertainty in data. This study introduces an innovative framework employing statistical change point analysis to generate labelled trees that encode vessel dimensions and their associated uncertainty from medical images. To test this framework, we explore the impact of uncertainty in 1D haemodynamic predictions in a systemic and pulmonary arterial network. Simulations explore haemodynamic variations resulting from changes in vessel dimensions and segmentation; the latter is achieved by analysing multiple segmentations of the same images. Results demonstrate the importance of accurately defining vessel radii and lengths when generating high-fidelity patient-specific haemodynamics models. KEY POINTS: This study introduces novel algorithms for generating labelled directed trees from medical images, focusing on accurate junction node placement and radius extraction using change points to provide haemodynamic predictions with uncertainty within expected measurement error. Geometric features, such as vessel dimension (length and radius) and network size, significantly impact pressure and flow predictions in both pulmonary and aortic arterial networks. Standardizing networks to a consistent number of vessels is crucial for meaningful comparisons and decreases haemodynamic uncertainty. Change points are valuable to understanding structural transitions in vascular data, providing an automated and efficient way to detect shifts in vessel characteristics and ensure reliable extraction of representative vessel radii.

Computational framework for the generation of one-dimensional vascular models accounting for uncertainty in networks extracted from medical images

One-dimensional (1D) cardiovascular models offer a non-invasive method to answer medical questions, including predictions of wave-reflection, shear stress, functional flow reserve, vascular resistance, and compliance. This model type can predict patient-specific outcomes by solving 1D fluid dynamics equations in geometric networks extracted from medical images. However, the inherent uncertainty in in-vivo imaging introduces variability in network size and vessel dimensions, affecting hemodynamic predictions. Understanding the influence of variation in image-derived properties is essential to assess the fidelity of model predictions. Numerous programs exist to render three-dimensional surfaces and construct vessel centerlines. Still, there is no exact way to generate vascular trees from the centerlines while accounting for uncertainty in data. This study introduces an innovative framework employing statistical change point analysis to generate labeled trees that encode vessel dimensions and their associated uncertainty from medical images. To test this framework, we explore the impact of uncertainty in 1D hemodynamic predictions in a systemic and pulmonary arterial network. Simulations explore hemodynamic variations resulting from changes in vessel dimensions and segmentation; the latter is achieved by analyzing multiple segmentations of the same images. Results demonstrate the importance of accurately defining vessel radii and lengths when generating high-fidelity patient-specific hemodynamics models.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstructs flow and increases pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the small pulmonary arteries due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we applied a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) approach to investigate the impact of CTEPH on pulmonary artery stiffening, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) in extralobar (main, right, and left) pulmonary arteries and intralobar (distal to the extralobar) arteries. Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar and intralobar arteries. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI, quantified by the novel parameter φ , which is related to thrombogenicity. Our analysis reveals a strong positive correlation between increases in mean pulmonary artery pressure (mPAP) and φ from baseline to CTEPH in individual subjects, which supports the suggestion that increased φ drives disease severity. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

A computational study of aortic reconstruction in single ventricle patients

Patients with hypoplastic left heart syndrome (HLHS) are born with an underdeveloped left heart. They typically receive a sequence of surgeries that result in a single ventricle physiology called the Fontan circulation. While these patients usually survive into early adulthood, they are at risk for medical complications, partially due to their lower than normal cardiac output, which leads to insufficient cerebral and gut perfusion. While clinical imaging data can provide detailed insight into cardiovascular function within the imaged region, it is difficult to use these data for assessing deficiencies in the rest of the body and for deriving blood pressure dynamics. Data from patients used in this paper include three-dimensional, magnetic resonance angiograms (MRA), time-resolved phase contrast cardiac magnetic resonance images (4D-MRI) and sphygmomanometer blood pressure measurements. The 4D-MRI images provide detailed insight into velocity and flow in vessels within the imaged region, but they cannot predict flow in the rest of the body, nor do they provide values of blood pressure. To remedy these limitations, this study combines the MRA, 4D-MRI, and pressure data with 1D fluid dynamics models to predict hemodynamics in the major systemic arteries, including the cerebral and gut vasculature. A specific focus is placed on studying the impact of aortic reconstruction occurring during the first surgery that results in abnormal vessel morphology. To study these effects, we compare simulations for an HLHS patient with simulations for a matched control patient that has double outlet right ventricle (DORV) physiology with a native aorta. Our results show that the HLHS patient has hypertensive pressures in the brain as well as reduced flow to the gut. Wave intensity analysis suggests that the HLHS patient has irregular circulatory function during light upright exercise conditions and that predicted wall shear stresses are lower than normal, suggesting the HLHS patient may have hypertension.

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.

Hemodynamic analysis of a curved artery based on microcirculation boundary

BACKGROUND: Microcirculation plays a key role in regulating blood flow but is not considered in previous research of hemodynamics. OBJECTIVE:  A curved artery model is established to study its hemodynamic characteristics based on microcirculation boundary. METHODS: The hemodynamic model of a curved artery is constructed and simulated by computational fluid dynamics. The curved artery model is simulated by fluid-structure interaction. At the same time, a porous medium is used to simulate microcirculation as the outlet boundary. RESULTS: The distribution characteristics of the blood flow velocity, the pressure and the wall shear stress in different sections at different time of the cardiac cycle are obtained. The results show that the velocities in curved arteries decrease and the pressures gradually increase. The blood flow velocity waveform and value are affected and they are sensitive to the microcirculation boundary. However, the pressure value is only affected by the microcirculation function. CONCLUSIONS:  This work is useful for researchers to deeply understand the hemodynamic characteristics of curved arteries. There is important clinical significance to analyze the pathogenesis of cardiovascular disease considering microcirculation function and its coupling effect.

Numerical predictions of shear stress and cyclic stretch in pulmonary hypertension due to left heart failure

Isolated post-capillary pulmonary hypertension (Ipc-PH) occurs due to left heart failure, which contributes to 1 out of every 9 deaths in the United States. In some patients, through unknown mechanisms, Ipc-PH transitions to combined pre-/post-capillary PH (Cpc-PH) and is associated with a dramatic increase in mortality. Altered mechanical forces and subsequent biological signaling in the pulmonary vascular bed likely contribute to the transition from Ipc-PH to Cpc-PH. However, even in a healthy pulmonary circulation, the mechanical forces in the smallest vessels (the arterioles, capillary bed, and venules) have not been quantitatively defined. This study is the first to examine this question via a computational fluid dynamics model of the human pulmonary arteries, arterioles, venules, and veins. Using this model, we predict temporal and spatial dynamics of cyclic stretch and wall shear stress with healthy and diseased hemodynamics. In the normotensive case for large vessels, numerical simulations show that large arteries have higher pressure and flow than large veins, as well as more pronounced changes in area throughout the cardiac cycle. In the microvasculature, shear stress increases and cyclic stretch decreases as vessel radius decreases. When we impose an increase in left atrial pressure to simulate Ipc-PH, shear stress decreases and cyclic stretch increases as compared to the healthy case. Overall, this model predicts pressure, flow, shear stress, and cyclic stretch that providing a way to analyze and investigate hypotheses related to disease progression in the pulmonary circulation.

Geometric Uncertainty in Patient-Specific Cardiovascular Modeling with Convolutional Dropout Networks

We propose a novel approach to generate samples from the conditional distribution of patient-specific cardiovascular models given a clinically aquired image volume. A convolutional neural network architecture with dropout layers is first trained for vessel lumen segmentation using a regression approach, to enable Bayesian estimation of vessel lumen surfaces. This network is then integrated into a path-planning patient-specific modeling pipeline to generate families of cardiovascular models. We demonstrate our approach by quantifying the effect of geometric uncertainty on the hemodynamics for three patient-specific anatomies, an aorto-iliac bifurcation, an abdominal aortic aneurysm and a sub-model of the left coronary arteries. A key innovation introduced in the proposed approach is the ability to learn geometric uncertainty directly from training data. The results show how geometric uncertainty produces coefficients of variation comparable to or larger than other sources of uncertainty for wall shear stress and velocity magnitude, but has limited impact on pressure. Specifically, this is true for anatomies characterized by small vessel sizes, and for local vessel lesions seen infrequently during network training.

A multiscale model of vascular function in chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) is caused by recurrent or unresolved pulmonary thromboemboli, leading to perfusion defects and increased arterial wave reflections. CTEPH treatment aims to reduce pulmonary arterial pressure and reestablish adequate lung perfusion, yet patients with distal lesions are inoperable by standard surgical intervention. Instead, these patients undergo balloon pulmonary angioplasty (BPA), a multisession, minimally invasive surgery that disrupts the thromboembolic material within the vessel lumen using a catheter balloon. However, there still lacks an integrative, holistic tool for identifying optimal target lesions for treatment. To address this insufficiency, we simulate CTEPH hemodynamics and BPA therapy using a multiscale fluid dynamics model. The large pulmonary arterial geometry is derived from a computed tomography (CT) image, whereas a fractal tree represents the small vessels. We model ring- and web-like lesions, common in CTEPH, and simulate normotensive conditions and four CTEPH disease scenarios; the latter includes both large artery lesions and vascular remodeling. BPA therapy is simulated by simultaneously reducing lesion severity in three locations. Our predictions mimic severe CTEPH, manifested by an increase in mean proximal pulmonary arterial pressure above 20 mmHg and prominent wave reflections. Both flow and pressure decrease in vessels distal to the lesions and increase in unobstructed vascular regions. We use the main pulmonary artery (MPA) pressure, a wave reflection index, and a measure of flow heterogeneity to select optimal target lesions for BPA. In summary, this study provides a multiscale, image-to-hemodynamics pipeline for BPA therapy planning for patients with inoperable CTEPH. NEW & NOTEWORTHY This article presents novel computational framework for predicting pulmonary hemodynamics in chronic thromboembolic pulmonary hypertension. The mathematical model is used to identify the optimal target lesions for balloon pulmonary angioplasty, combining simulated pulmonary artery pressure, wave intensity analysis, and a new quantitative metric of flow heterogeneity.

Markov chain Monte Carlo with Gaussian processes for fast parameter estimation and uncertainty quantification in a 1D fluid-dynamics model of the pulmonary circulation

The past few decades have witnessed an explosive synergy between physics and the life sciences. In particular, physical modelling in medicine and physiology is a topical research area. The present work focuses on parameter inference and uncertainty quantification in a 1D fluid-dynamics model for quantitative physiology: the pulmonary blood circulation. The practical challenge is the estimation of the patient-specific biophysical model parameters, which cannot be measured directly. In principle this can be achieved based on a comparison between measured and predicted data. However, predicting data requires solving a system of partial differential equations (PDEs), which usually have no closed-form solution, and repeated numerical integrations as part of an adaptive estimation procedure are computationally expensive. In the present article, we demonstrate how fast parameter estimation combined with sound uncertainty quantification can be achieved by a combination of statistical emulation and Markov chain Monte Carlo (MCMC) sampling. We compare a range of state-of-the-art MCMC algorithms and emulation strategies, and assess their performance in terms of their accuracy and computational efficiency. The long-term goal is to develop a method for reliable disease prognostication in real time, and our work is an important step towards an automatic clinical decision support system.

Modeling and Reducing the Effect of Geometric Uncertainties in Intracranial Aneurysms with Polynomial Chaos Expansion, Data Decomposition, and 4D-Flow MRI

PURPOSE

Variations in the vessel radius of segmented surfaces of intracranial aneurysms significantly influence the fluid velocities given by computer simulations. It is important to generate models that capture the effect of these variations in order to have a better interpretation of the numerically predicted hemodynamics. Also, it is highly relevant to develop methods that combine experimental observations with uncertainty modeling to get a closer approximation to the blood flow behavior.

METHODS

This work applies polynomial chaos expansion to model the effect of geometric uncertainties on the simulated fluid velocities of intracranial aneurysms. The radius of the vessel is defined as the uncertainty variable. Proper orthogonal decomposition is applied to characterize the solution space of fluid velocities. Next, a process of projecting the 4D-Flow MRI velocities on the basis vectors followed by coefficient mapping using generalized dynamic mode decomposition enables the merging of 4D-Flow MRI with the uncertainty propagated fluid velocities.

RESULTS

Polynomial chaos expansion propagates the fluid velocities with an error of 2% in velocity magnitude relative to computer simulations. Also, the bifurcation region (or impingement location) shows a standard deviation of 0.17 m/s (since an available reported variance in the vessel radius is adopted to model the uncertainty, the expected standard deviation may be different). Numerical phantom experiments indicate that the proposed approach reconstructs the fluid velocities with 0.3% relative error in presence of geometric uncertainties.

CONCLUSION

Polynomial chaos expansion is an effective approach to propagate the effect of the uncertainty variable in the blood flow velocities of intracranial aneurysms. Merging 4D-Flow MRI and uncertainty propagated fluid velocities leads to more realistic flow trends relative to ignoring the uncertainty in the vessel radius.

Assessing model mismatch and model selection in a Bayesian uncertainty quantification analysis of a fluid-dynamics model of pulmonary blood circulation

This study uses Bayesian inference to quantify the uncertainty of model parameters and haemodynamic predictions in a one-dimensional pulmonary circulation model based on an integration of mouse haemodynamic and micro-computed tomography imaging data. We emphasize an often neglected, though important source of uncertainty: in the mathematical model form due to the discrepancy between the model and the reality, and in the measurements due to the wrong noise model (jointly called 'model mismatch'). We demonstrate that minimizing the mean squared error between the measured and the predicted data (the conventional method) in the presence of model mismatch leads to biased and overly confident parameter estimates and haemodynamic predictions. We show that our proposed method allowing for model mismatch, which we represent with Gaussian processes, corrects the bias. Additionally, we compare a linear and a nonlinear wall model, as well as models with different vessel stiffness relations. We use formal model selection analysis based on the Watanabe Akaike information criterion to select the model that best predicts the pulmonary haemodynamics. Results show that the nonlinear pressure-area relationship with stiffness dependent on the unstressed radius predicts best the data measured in a control mouse.

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.

Structural and hemodynamic properties of murine pulmonary arterial networks under hypoxia-induced pulmonary hypertension

Detection and monitoring of patients with pulmonary hypertension, defined as a mean blood pressure in the main pulmonary artery above 25 mmHg, requires a combination of imaging and hemodynamic measurements. This study demonstrates how to combine imaging data from microcomputed tomography images with hemodynamic pressure and flow waveforms from control and hypertensive mice. Specific attention is devoted to developing a tool that processes computed tomography images, generating subject-specific arterial networks in which one-dimensional fluid dynamics modeling is used to predict blood pressure and flow. Each arterial network is modeled as a directed graph representing vessels along the principal pathway to ensure perfusion of all lobes. The one-dimensional model couples these networks with structured tree boundary conditions representing the small arteries and arterioles. Fluid dynamics equations are solved in this network and compared to measurements of pressure in the main pulmonary artery. Analysis of microcomputed tomography images reveals that the branching ratio is the same in the control and hypertensive animals, but that the vessel length-to-radius ratio is significantly lower in the hypertensive animals. Fluid dynamics predictions show that in addition to changed network geometry, vessel stiffness is higher in the hypertensive animal models than in the control models.