Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation

Efficient uncertainty quantification in a spatially multiscale model of pulmonary arterial and venous hemodynamics

Pulmonary hypertension (PH) is a debilitating disease that alters the structure and function of both the proximal and distal pulmonary vasculature. This alters pressure-flow relationships in the pulmonary arterial and venous trees, though there is a critical knowledge gap in the relationships between proximal and distal hemodynamics in disease. Multiscale computational models enable simulations in both the proximal and distal vasculature. However, model inputs and measured data are inherently uncertain, requiring a full analysis of the sensitivity and uncertainty of the model. Thus, this study quantifies model sensitivity and output uncertainty in a spatially multiscale, pulse-wave propagation model of pulmonary hemodynamics. The model includes fifteen proximal arteries and twelve proximal veins, connected by a two-sided, structured tree model of the distal vasculature. We use polynomial chaos expansions to expedite sensitivity and uncertainty quantification analyses and provide results for both the proximal and distal vasculature. We quantify uncertainty in blood pressure, blood flow rate, wave intensity, wall shear stress, and cyclic stretch. The latter two are important stimuli for endothelial cell mechanotransduction. We conclude that, while nearly all the parameters in our system have some influence on model predictions, the parameters describing the density of the microvascular beds have the largest effects on all simulated quantities in both the proximal and distal arterial and venous circulations.

Estimating pulmonary arterial remodeling via an animal-specific computational model of pulmonary artery stenosis

Pulmonary artery stenosis (PAS) often presents in children with congenital heart disease, altering blood flow and pressure during critical periods of growth and development. Variability in stenosis onset, duration, and severity result in variable growth and remodeling of the pulmonary vasculature. Computational fluid dynamics (CFD) models enable investigation into the hemodynamic impact and altered mechanics associated with PAS. In this study, a one-dimensional (1D) fluid dynamics model was used to simulate hemodynamics throughout the pulmonary arteries of individual animals. The geometry of the large pulmonary arteries was prescribed by animal-specific imaging, whereas the distal vasculature was simulated by a three-element Windkessel model at each terminal vessel outlet. Remodeling of the pulmonary vasculature, which cannot be measured in vivo, was estimated via model-fitted parameters. The large artery stiffness was significantly higher on the left side of the vasculature in the left pulmonary artery (LPA) stenosis group, but neither side differed from the sham group. The sham group exhibited a balanced distribution of total distal vascular resistance, whereas the left side was generally larger in the LPA stenosis group, with no significant differences between groups. In contrast, the peripheral compliance on the right side of the LPA stenosis group was significantly greater than the corresponding side of the sham group. Further analysis indicated the underperfused distal vasculature likely moderately decreased in radius with little change in stiffness given the increase in thickness observed with histology. Ultimately, our model enables greater understanding of pulmonary arterial adaptation due to LPA stenosis and has potential for use as a tool to noninvasively estimate remodeling of the pulmonary vasculature.

Application and reduction of a nonlinear hyperelastic wall model capturing ex vivo relationships between fluid pressure, area, and wall thickness in normal and hypertensive murine left pulmonary arteries

Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery. The model relating blood pressure and area is calibrated using ex vivo measurements of vessel diameter and wall thickness changes, under controlled pressure conditions, in left pulmonary arteries isolated from control and hypertensive mice. A two-layer, hyperelastic, and anisotropic model incorporating residual stresses is formulated using the Holzapfel-Gasser-Ogden model. Complex relations predicting vessel area and wall thickness with increasing blood pressure are derived and calibrated using the data. Sensitivity analysis, parameter estimation, subset selection, and physical plausibility arguments are used to systematically reduce the 16-parameter model to one in which a much smaller subset of identifiable parameters is estimated via solution of an inverse problem. Our final reduced one layer model includes a single set of three elastic moduli. Estimated ranges of these parameters demonstrate that nonlinear stiffening is dominated by elastin in the control animals and by collagen in the hypertensive animals. The pressure-area relation developed in this novel manner has potential impact on one-dimensional fluids network models of vessel wall remodeling in the presence of cardiovascular disease.

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.

Efficient Uncertainty Quantification in a Multiscale Model of Pulmonary Arterial and Venous Hemodynamics

Computational hemodynamics models are becoming increasingly useful in the management and prognosis of complex, multiscale pathologies, including those attributed to the development of pulmonary vascular disease. However, diseases like pulmonary hypertension are heterogeneous, and affect both the proximal arteries and veins as well as the microcirculation. Simulation tools and the data used for model calibration are also inherently uncertain, requiring a full analysis of the sensitivity and uncertainty attributed to model inputs and outputs. Thus, this study quantifies model sensitivity and output uncertainty in a multiscale, pulse-wave propagation model of pulmonary hemodynamics. Our pulmonary circuit model consists of fifteen proximal arteries and twelve proximal veins, connected by a two-sided, structured tree model of the distal vasculature. We use polynomial chaos expansions to expedite the sensitivity and uncertainty quantification analyses and provide results for both the proximal and distal vasculature. Our analyses provide uncertainty in blood pressure, flow, and wave propagation phenomenon, as well as wall shear stress and cyclic stretch, both of which are important stimuli for endothelial cell mechanotransduction. We conclude that, while nearly all the parameters in our system have some influence on model predictions, the parameters describing the density of the microvascular beds have the largest effects on all simulated quantities in both the proximal and distal circulation.

A multiscale framework for defining homeostasis in distal vascular trees: applications to the pulmonary circulation

Pulmonary arteries constitute a low-pressure network of vessels, often characterized as a bifurcating tree with heterogeneous vessel mechanics. Understanding the vascular complexity and establishing homeostasis is important to study diseases such as pulmonary arterial hypertension (PAH). The onset and early progression of PAH can be traced to changes in the morphometry and structure of the distal vasculature. Coupling hemodynamics with vessel wall growth and remodeling (G&R) is crucial for understanding pathology at distal vasculature. Accordingly, the goal of this study is to provide a multiscale modeling framework that embeds the essential features of arterial wall constituents coupled with the hemodynamics within an arterial network characterized by an extension of Murray's law. This framework will be used to establish the homeostatic baseline characteristics of a pulmonary arterial tree, including important parameters such as vessel radius, wall thickness and shear stress. To define the vascular homeostasis and hemodynamics in the tree, we consider two timescales: a cardiac cycle and a longer period of vascular adaptations. An iterative homeostatic optimization, which integrates a metabolic cost function minimization, the stress equilibrium, and hemodynamics, is performed at the slow timescale. In the fast timescale, the pulsatile blood flow dynamics is described by a Womersley's deformable wall analytical solution. Illustrative examples for symmetric and asymmetric trees are presented that provide baseline characteristics for the normal pulmonary arterial vasculature. The results are compared with diverse literature data on morphometry, structure, and mechanics of pulmonary arteries. The developed framework demonstrates a potential for advanced parametric studies and future G&R and hemodynamics modeling of PAH.

Computer simulation of surgical interventions for the treatment of refractory pulmonary hypertension

This paper describes computer models of three interventions used for treating refractory pulmonary hypertension (RPH). These procedures create either an atrial septal defect, a ventricular septal defect or, in the case of a Potts shunt, a patent ductus arteriosus. The aim in all three cases is to generate a right-to-left shunt, allowing for either pressure or volume unloading of the right side of the heart in the setting of right ventricular failure, while maintaining cardiac output. These shunts are created, however, at the expense of introducing de-oxygenated blood into the systemic circulation, thereby lowering the systemic arterial oxygen saturation. The models developed in this paper are based on compartmental descriptions of human hemodynamics and oxygen transport. An important parameter included in our models is the cross-sectional area of the surgically created defect. Numerical simulations are performed to compare different interventions and various shunt sizes and to assess their impact on hemodynamic variables and oxygen saturations. We also create a model for exercise and use it to study exercise tolerance in simulated pre-intervention and post-intervention RPH patients.

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.

Cardiovascular Regulation During Acute Gravitational Changes with Exhaling on Exertion

During gravitational changes or changes in the direction of action in relation to the body, fluid displacements can be observed. In special cases different breathing maneuvers (e. g., exhaling on exertion; Ex-Ex) are used to counteract acute fluid shifts. Both factors have a significant impact on cardiovascular regulation. Eight healthy male subjects were tested on a tilt seat, long arm human centrifuge, and parabolic flight. The work aims to investigate the effect of exhaling on exertion on the cardiovascular regulation during acute gravitational changes compared to normal breathing. Possible interactions and differences between conditions (Ex-Ex, normal breathing) for the parameters V’O₂ , V’_E , HR, and SV were analysed over a 40 s period by a three-way ANOVA. Significant (p≤0.05) effects for all main factors and interactions between condition and time as well as maneuver and time were found for all variables. The exhaling on exertion maneuver had a significant influence on the cardiovascular response during acute gravitational and positional changes. For example, the significant increase of V’O2 at the end of the exhalation on exertion maneuver indicates an increased lung circulation as a result of the maneuver.

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.

Simulating Multi-Scale Pulmonary Vascular Function by Coupling Computational Fluid Dynamics With an Anatomic Network Model

The function of the pulmonary circulation is truly multi-scale, with blood transported through vessels from centimeter to micron scale. There are scale-dependent mechanisms that govern the flow in the pulmonary vascular system. However, very few computational models of pulmonary hemodynamics capture the physics of pulmonary perfusion across the spatial scales of functional importance in the lung. Here we present a multi-scale model that incorporates the 3-dimensional (3D) complexities of pulmonary blood flow in the major vessels, coupled to an anatomically-based vascular network model incorporating the multiple contributing factors to capillary perfusion, including gravity. Using the model we demonstrate how we can predict the impact of vascular remodeling and occlusion on both macro-scale functional drivers (flow distribution between lungs, and wall shear stress) and micro-scale contributors to gas exchange. The model predicts interactions between 3D and 1D models that lead to a redistribution of blood between postures, both on a macro- and a micro-scale. This allows us to estimate the effect of posture on left and right pulmonary artery wall shear stress, with predictions varying by 0.75-1.35 dyne/cm² between postures.

Closed-Loop Multiscale Computational Model of Human Blood Circulation. Applications to Ballistocardiography

Cardiac mechanical activity leads to periodic changes in the distribution of blood throughout the body, which causes micro-oscillations of the body's center of mass and can be measured by ballistocardiography (BCG). However, many of the BCG findings are based on parameters whose origins are poorly understood. Here, we generate simulated multidimensional BCG signals based on a more exhaustive and accurate computational model of blood circulation than previous attempts. This model consists in a closed loop 0D-1D multiscale representation of the human blood circulation. The 0D elements include the cardiac chambers, cardiac valves, arterioles, capillaries, venules, and veins, while the 1D elements include 55 systemic and 57 pulmonary arteries. The simulated multidimensional BCG signal is computed based on the distribution of blood in the different compartments and their anatomical position given by whole-body magnetic resonance angiography on a healthy young subject. We use this model to analyze the elements affecting the BCG signal on its different axes, allowing a better interpretation of clinical records. We also evaluate the impact of filtering and healthy aging on the BCG signal. The results offer a better view of the physiological meaning of BCG, as compared to previous models considering mainly the contribution of the aorta and focusing on longitudinal acceleration BCG. The shape of experimental BCG signals can be reproduced, and their amplitudes are in the range of experimental records. The contributions of the cardiac chambers and the pulmonary circulation are non-negligible, especially on the lateral and transversal components of the velocity BCG signal. The shapes and amplitudes of the BCG waveforms are changing with age, and we propose a scaling law to estimate the pulse wave velocity based on the time intervals between the peaks of the acceleration BCG signal. We also suggest new formulas to estimate the stroke volume and its changes based on the BCG signal expressed in terms of acceleration and kinetic energy.

A nonlinear multi-scale model for blood circulation in a realistic vascular system

In the last decade, numerical models have become an increasingly important tool in biological and medical science. Numerical simulations contribute to a deeper understanding of physiology and are a powerful tool for better diagnostics and treatment. In this paper, a nonlinear multi-scale model framework is developed for blood flow distribution in the full vascular system of an organ. We couple a quasi one-dimensional vascular graph model to represent blood flow in larger vessels and a porous media model to describe flow in smaller vessels and capillary bed. The vascular model is based on Poiseuille's Law, with pressure correction by elasticity and pressure drop estimation at vessels' junctions. The porous capillary bed is modelled as a two-compartment domain (artery and venous) using Darcy's Law. The fluid exchange between the artery and venous capillary bed compartments is defined as blood perfusion. The numerical experiments show that the proposed model for blood circulation: (i) is closely dependent on the structure and parameters of both the larger vessels and of the capillary bed, and (ii) provides a realistic blood circulation in the organ. The advantage of the proposed model is that it is complex enough to reliably capture the main underlying physiological function, yet highly flexible as it offers the possibility of incorporating various local effects. Furthermore, the numerical implementation of the model is straightforward and allows for simulations on a regular desktop computer.

High fidelity blood flow in a patient-specific arteriovenous fistula

An arteriovenous fistula, created by artificially connecting segments of a patient’s vasculature, is the preferred way to gain access to the bloodstream for kidney dialysis. The increasing power and availability of supercomputing infrastructure means that it is becoming more realistic to use simulations to help identify the best type and location of a fistula for a specific patient. We describe a 3D fistula model that uses the lattice Boltzmann method to simultaneously resolve blood flow in patient-specific arteries and veins. The simulations conducted here, comprising vasculatures of the whole forearm, demonstrate qualified validation against clinical data. Ongoing research to further encompass complex biophysics on realistic time scales will permit the use of human-scale physiological models for basic and clinical medicine.

Sensitivity analysis and uncertainty quantification of 1-D models of pulmonary hemodynamics in mice under control and hypertensive conditions

Pulmonary hypertension (PH), defined as an elevated mean blood pressure in the main pulmonary artery (MPA) at rest, is associated with vascular remodeling of both large and small arteries. PH has several sub-types that are all linked to high mortality rates. In this study, we use a one-dimensional (1-D) fluid dynamics model driven by in vivo measurements of MPA flow to understand how model parameters and network size influence MPA pressure predictions in the presence of PH. We compare model predictions with in vivo MPA pressure measurements from a control and a hypertensive mouse and analyze results in three networks of increasing complexity, extracted from micro-computed tomography (micro-CT) images. We introduce global scaling factors for boundary condition parameters and perform local and global sensitivity analysis to calculate parameter influence on model predictions of MPA pressure and correlation analysis to determine a subset of identifiable parameters. These are inferred using frequentist optimization and Bayesian inference via the Delayed Rejection Adaptive Metropolis (DRAM) algorithm. Frequentist and Bayesian uncertainty is computed for model parameters and MPA pressure predictions. Results show that MPA pressure predictions are most sensitive to distal vascular resistance and that parameter influence changes with increasing network complexity. Our outcomes suggest that PH leads to increased vascular stiffness and decreased peripheral compliance, congruent with clinical observations.

Plasticity and Enzymatic Degradation Coupled With Volumetric Growth in Pulmonary Hypertension Progression

Pulmonary hypertension (PH) is one of the least understood and highly elusive cardiovascular conditions associated with elevated pulmonary arterial pressure. Although the disease mechanisms are not completely understood, evidence has accumulated from human and animal studies that irreversible processes of pulmonary arterial wall damage, compensated by stress-mediated growth, play critical roles in eliciting the mechanisms of disease progression. The aim of this study is to develop a thermodynamic modeling structure of the pulmonary artery to consider coupled plastic-degradation-growth irreversible processes to investigate the mechanical roles of the dissipative phenomena in the disease progression. The proposed model performs a model parameter study of plastic deformation and degradation processes coupled with dissipative growth subjected to elevated pulmonary arterial pressure and computationally generates in silico simulations of PH progression using the clinical features of PH, found in human morphological and mechanical data. The results show that considering plastic deformation can provide a much better fitting of the ex vivo inflation tests than a widely used pure hyperelastic model in higher pressure conditions. In addition, the parameter sensitivity study illustrates that arterial damage and growth cause the increased stiffness, and the full simulation (combining elastic-plastic-degradation-growth models) reveals a key postpathological recovery process of compensating vessel damage by vascular adaptation by reducing the rate of vessel dilation and mediating vascular wall stress. Finally, the simulation results of luminal enlargement, arterial thickening, and arterial stiffness for an anisotropic growth are found to be close to the values from the literature.

Analysis of flow resistance in the pulmonary arterial circulation: implications for hypoxic pulmonary vasoconstriction

Hypoxic pulmonary vasoconstriction (HPV) plays an essential role in distributing blood in the lung to enhance ventilation-perfusion matching and blood oxygenation. In this study, a theoretical model of the pulmonary vasculature is used to predict the effects of vasoconstriction over specified ranges of vessel diameters on pulmonary vascular resistance (PVR). The model is used to evaluate the ability of hypothesized mechanisms of HPV to account for observed levels of PVR elevation during hypoxia. The vascular structure from pulmonary arteries to capillaries is represented using scaling laws. Vessel segments are modeled as resistive elements and blood flow rates are computed from physical principles. Direct vascular responses to intravascular oxygen levels have been proposed as a mechanism of HPV. In the lung, significant changes in oxygen level occur only in vessels less than 60 μm in diameter. The model shows that observed levels of hypoxic vasoconstriction in these vessels alone cannot account for the elevation of PVR associated with HPV. However, the elevation in PVR associated with HPV can be accounted for if larger upstream vessels also constrict. These results imply that upstream signaling by conducted responses to engage constriction of arterioles plays an essential role in the elevation of PVR during HPV.NEW & NOTEWORTHY A theoretical model of the pulmonary vasculature is used to predict the effects of vasoconstriction over specified ranges of vessel diameters on pulmonary vascular resistance (PVR). The model shows that observed levels of hypoxic vasoconstriction in terminal vessels cannot account for the elevation of PVR associated with hypoxic pulmonary vasoconstriction (HPV). Upstream signaling by conducted responses to engage constriction of arterioles, therefore, plays an essential role in the elevation of PVR during HPV.

A computational model of contributors to pulmonary hypertensive disease: impacts of whole lung and focal disease distributions

Pulmonary hypertension has multiple etiologies and so can be difficult to diagnose, prognose, and treat. Diagnosis is typically made via invasive hemodynamic measurements in the main pulmonary artery and is based on observed elevation of mean pulmonary artery pressure. This static mean pressure enables diagnosis, but does not easily allow assessment of the severity of pulmonary hypertension, nor the etiology of the disease, which may impact treatment. Assessment of the dynamic properties of pressure and flow data obtained from catheterization potentially allows more meaningful assessment of the strain on the right heart and may help to distinguish between disease phenotypes. However, mechanistic understanding of how the distribution of disease in the lung leading to pulmonary hypertension impacts the dynamics of blood flow in the main pulmonary artery and/or the pulmonary capillaries is lacking. We present a computational model of the pulmonary vasculature, parameterized to characteristic features of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension to help understand how the two conditions differ in terms of pulmonary vascular response to disease. Our model incorporates key features known to contribute to pulmonary vascular function in health and disease, including anatomical structure and multiple contributions from gravity. The model suggests that dynamic measurements obtained from catheterization potentially distinguish between distal and proximal vasculopathy typical of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension. However, the model suggests a non-linear relationship between these data and vascular structural changes typical of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension which may impede analysis of these metrics to distinguish between cohorts.

Computational simulation-derived hemodynamic and biomechanical properties of the pulmonary arterial tree early in the course of ventricular septal defects

Untreated ventricular septal defects (VSDs) can lead to pulmonary arterial hypertension (PAH) characterized by elevated pulmonary artery (PA) pressure and vascular remodeling, known as PAH associated with congenital heart disease (PAH-CHD). Though previous studies have investigated hemodynamic effects on vascular mechanobiology in late-stage PAH, hemodynamics leading to PAH-CHD initiation have not been fully quantified. We hypothesize that abnormal hemodynamics from left-to-right shunting in early stage VSDs affects PA biomechanical properties leading to PAH initiation. To model PA hemodynamics in healthy, small, moderate, and large VSD conditions prior to the onset of vascular remodeling, computational fluid dynamics simulations were performed using a 3D finite element model of a healthy 1-year-old's proximal PAs and a body-surface-area-scaled 0D distal PA tree. VSD conditions were modeled with increased pulmonary blood flow to represent degrees of left-to-right shunting. In the proximal PAs, pressure, flow, strain, and wall shear stress (WSS) increased with increasing VSD size; oscillatory shear index decreased with increasing VSD size in the larger PA vessels. WSS was higher in smaller diameter vessels and increased with VSD size, with the large VSD condition exhibiting WSS >100 dyn/cm[Formula: see text], well above values typically used to study dysfunctional mechanotransduction pathways in PAH. This study is the first to estimate hemodynamic and biomechanical metrics in the entire pediatric PA tree with VSD severity at the stage leading to PAH initiation and has implications for future studies assessing effects of abnormal mechanical stimuli on endothelial cells and vascular wall mechanics that occur during PAH-CHD initiation and progression.

A Probabilistic Approach for Contact Stability and Contact Safety Analysis of Robotic Intracardiac Catheter

The disturbances caused by the blood flow and tissue surface motions are major concerns during the motion planning of an intracardiac robotic catheter. Maintaining a stable and safe contact on the desired ablation point is essential for achieving effective lesions during the ablation procedure. In this paper, a probabilistic formulation of the contact stability and the contact safety for intravascular cardiac catheters under the blood flow and surface motion disturbances is presented. Probabilistic contact stability and contact safety metrics, employing a sample-based representation of the blood flow velocity distribution and the heart motion trajectory, are introduced. Finally, the contact stability and safety for an magnetic resonance imaging-actuated robotic catheter under main pulmonary artery blood flow disturbances and left ventricle surface motion disturbances are analyzed in simulation as example scenarios.

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.

Fluid-structure interaction in a fully coupled three-dimensional mitral-atrium-pulmonary model

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

Pulmonary hypertension in fibrosing idiopathic interstitial pneumonia: Uncertainties, challenges and opportunities

Pulmonary hypertension is a serious complication of chronic fibrosing idiopathic interstitial pneumonia (PH-fIIP) leading to greater morbidity and mortality. The pathophysiologic basis for PH in fIIP is not completely understood, but microvascular rarefaction may play a key role. Severe hypoxemia and reduced diffusion capacity are characteristic. Doppler echocardiography has limited diagnostic utility and right heart catheterization is required to confirm the diagnosis. Lung volumes can be minimally affected, and radiographic findings can be subtle, making the distinction from pulmonary arterial hypertension challenging. Several randomized controlled trials of pulmonary arterial hypertension targeted therapies have recently been completed. Endothelin-receptor antagonists have shown either no benefit or harm. Sildenafil may have some favorable short-term effects but does not appear to impact long-term outcomes. Riociguat treatment increased hospitalizations and mortality. A recent trial of inhaled treprostinil demonstrated improved exercise capacity, but the impact on long-term morbidity and mortality are unknown. Currently, the only viable option for improved survival is lung transplantation. Early referral is imperative to optimize post-transplant outcomes.

Integrative Computational Models of Lung Structure-Function Interactions

Anatomically based integrative models of the lung and their interaction with other key components of the respiratory system provide unique capabilities for investigating both normal and abnormal lung function. There is substantial regional variability in both structure and function within the normal lung, yet it remains capable of relatively efficient gas exchange by providing close matching of air delivery (ventilation) and blood delivery (perfusion) to regions of gas exchange tissue from the scale of the whole organ to the smallest continuous gas exchange units. This is despite remarkably different mechanisms of air and blood delivery, different fluid properties, and unique scale-dependent anatomical structures through which the blood and air are transported. This inherent heterogeneity can be exacerbated in the presence of disease or when the body is under stress. Current computational power and data availability allow for the construction of sophisticated data-driven integrative models that can mimic respiratory system structure, function, and response to intervention. Computational models do not have the same technical and ethical issues that can limit experimental studies and biomedical imaging, and if they are solidly grounded in physiology and physics they facilitate investigation of the underlying interaction between mechanisms that determine respiratory function and dysfunction, and to estimate otherwise difficult-to-access measures. © 2021 American Physiological Society. Compr Physiol 11:1501-1530, 2021.

Towards blood flow in the virtual human: efficient self-coupling of HemeLB

Many scientific and medical researchers are working towards the creation of a virtual human—a personalized digital copy of an individual—that will assist in a patient’s diagnosis, treatment and recovery. The complex nature of living systems means that the development of this remains a major challenge. We describe progress in enabling the HemeLB lattice Boltzmann code to simulate 3D macroscopic blood flow on a full human scale. Significant developments in memory management and load balancing allow near linear scaling performance of the code on hundreds of thousands of computer cores. Integral to the construction of a virtual human, we also outline the implementation of a self-coupling strategy for HemeLB. This allows simultaneous simulation of arterial and venous vascular trees based on human-specific geometries.

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.

Analysis of Contact Stability and Contact Safety of a Robotic Intravascular Cardiac Catheter under Blood Flow Disturbances

This paper studies the contact stability and contact safety of a robotic intravascular cardiac catheter under blood flow disturbances while in contact with tissue surface. A probabilistic blood flow disturbance model, where the blood flow drag forces on the catheter body are approximated using a quasi-static model, is introduced. Using this blood flow disturbance model, probabilistic contact stability and contact safety metrics, employing a sample based representation of the blood flow velocity distribution, are proposed. Finally, the contact stability and contact safety of a MRI-actuated robotic catheter are analyzed using these models in a specific example scenario under left pulmonary inferior vein (LIV) blood flow disturbances.

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.

Coronary arteries hemodynamics: effect of arterial geometry on hemodynamic parameters causing atherosclerosis

Coronary arteries have high curvatures, and hence, flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing realistic curvatures and stenosis and (ii) patient-specific coronary arteries with stenoses. Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Secondly, we have similarly determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), which were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis; our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves. This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient's coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is in our showing how in actual coronary arteries (based on their CT imaging) curvilinearity and narrowing complications affect the computed WSS and WPG, associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery. We want to go one step further and provide clinical application of our research work. For that, we are offering to cardiologists worldwide to carry out hemodynamic analysis of the medically imaged coronary arteries of their patients and compute the values of the hemodynamic parameters of WSS and WPG, so as to provide them an assessment of the risk of atherosclerosis for their patients. Graphical abstract Theme and aims: Coronary arteries have high curvatures, and hence flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing curvatures and stenosis, and (ii) patient-specific coronary arteries with stenoses. Methods and results: Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Then, we have determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), that were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis, as illustrated in Figure 1, in which the locations of the stenoses are highlighted by arrows. Figure 1: Three-dimensional CT visualization of arteries in patients with suspected coronary disease. The arteries can be seen as a combination of various curved segments with stenoses at unspecific locations highlighted by arrows. Our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves, as depicted in Figure 2. Therein, the encapsulations show (i) high WSS, and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Figure 2: WSS and WPG surface plot of realistic arteries (a), (b), (c) and (d), wherein the small squared parts are enlarged to show the detailed localized contour plots at the stenotic regions. Therein, the circular encapsulations show (i) high WSS and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Conclusion and novelty: This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is our showing how in actual coronary arteries (based on their CT imaging), curvilinearity and narrowing complications affect the computed WSS and WPG associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery.

Paediatric haemodynamic modelling: development and experimental validation using quantitative flow MRI

Background

Congenital vascular disease is one of the leading causes of death in paediatric age. Despite the importance of paediatric haemodynamics, large investigations have been devoted to the evaluation of circulation in adults. The novelty of this study consists in the development of a well calibrated mathematical model of cardiovascular circulation in paediatric subjects. To reach the purpose, a model for adult circulation was modified and recalibrated with experimental data and literature from children to be able to calculate the flow rates and pressures in the brain and neck.

Methods

The haemodynamic model simulates the 76 main arteries, together with the main veins in brain and neck. A proper magnetic resonance imaging (MRI) dataset of 29 volunteers aged 12 ± 5 years (mean ± standard deviation) was used to extract age-dependent physiological and clinical parameters such as heart rate, flow rate, vessel cross section area, and blood pressure. The computational model was calibrated using such experimental data. The paediatric and adult model results were compared.

Results

Increase of the vessels stiffness due to aging contributes to a flow rate decrease while blood pressure increases. In accordance, our simulation results show about 16% decrease in mean pressure of internal jugular vein in paediatric rather than adult subjects. The model outcomes indicated about 88% correlation with MRI data.

Conclusions

The mathematical model simulates the paediatric head and neck blood circulation. The model provides detailed information of human haemodynamics including arterial and venous network to study both paediatric and adult blood circulation.

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

Computational fluid dynamics (CFD) models are emerging tools for assisting in diagnostic assessment of cardiovascular disease. Recent advances in image segmentation have made subject-specific modelling of the cardiovascular system a feasible task, which is particularly important in the case of pulmonary hypertension, requiring a combination of invasive and non-invasive procedures for diagnosis. Uncertainty in image segmentation propagates to CFD model predictions, making the quantification of segmentation-induced uncertainty crucial for subject-specific models. This study quantifies the variability of one-dimensional CFD predictions by propagating the uncertainty of network geometry and connectivity to blood pressure and flow predictions. We analyse multiple segmentations of a single, excised mouse lung using different pre-segmentation parameters. A custom algorithm extracts vessel length, vessel radii and network connectivity for each segmented pulmonary network. Probability density functions are computed for vessel radius and length and then sampled to propagate uncertainties to haemodynamic predictions in a fixed network. In addition, we compute the uncertainty of model predictions to changes in network size and connectivity. Results show that variation in network connectivity is a larger contributor to haemodynamic uncertainty than vessel radius and length.

Anatomically based simulation of hepatic perfusion in the human liver

Liver structures of a healthy subject are digitised and segmented from computed tomography (CT) images, and hepatic perfusion is modelled in the hepatic artery and portal vein of the healthy subject with structured tree-based outflow boundary conditions. This self-similar structured tree is widely used in the literature, eg, blood flow simulation in larger systemic arteries and cerebral circulation, and is used in this study to model the effect of the smaller hepatic arteries and arterioles, as well as the smaller hepatic portal veins and portal venules. Physiologically reasonable results are obtained. Since the structured tree terminates at the size of the microvasculature system in liver lobules, the structured tree boundary condition will enable the proposed organ-level model of hepatic arterial flow to be easily connected to tissue-level models of liver lobules. Blood flow in the hepatic vein is also modelled in this subject with three-element Windkessel model as outflow boundary conditions. The benefit of integrating the perfusion in all hepatic vascular vessels is that it helps us analyse some complicated clinical phenomenon more efficiently, eg, one possible application is to obtain the portal pressure gradient (PPG) to help examine the reliability of hepatic venous pressure gradient (HVPG) as an indirect measure of portal pressure. Moreover, since four to six generations of hepatic vessels, which are sufficient for liver classification analysis, were employed in the model, this study is setting the computational foundation of a potentially handy surgical tool.

Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid-only to coupled fluid-structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

Pulmonary Vascular Dynamics

The pulmonary circulation carries deoxygenated blood from the systemic veins through the pulmonary arteries to be oxygenated in the capillaries that line the walls of the pulmonary alveoli. The pulmonary circulation carries the cardiac output with a relatively low driving pressure, and so differs considerably in structure and function from the systemic circulation to maintain a low-resistance vascular system. The pulmonary circulation is often considered to be a quasi-static system in both experimental and computational studies of pulmonary perfusion and its matching to ventilation (air flow) for exchange. However, the system is highly dynamic, with cardiac output and regional perfusion changing with posture, exercise, and over time. Here we review this dynamic system, with a focus on understanding the physiology of pulmonary vascular dynamics across spatial and temporal scales, and the changes to these dynamics that are reflective of disease. © 2019 American Physiological Society. Compr Physiol 9:1081-1100, 2019.

EXPRESS: Statement on imaging and pulmonary hypertension from the Pulmonary Vascular Research Institute (PVRI)

Pulmonary hypertension (PH) is highly heterogeneous and despite treatment advances it remains a life-shortening condition. There have been significant advances in imaging technologies, but despite evidence of their potential clinical utility, practice remains variable, dependent in part on imaging availability and expertise. This statement summarizes current and emerging imaging modalities and their potential role in the diagnosis and assessment of suspected PH. It also includes a review of commonly encountered clinical and radiological scenarios, and imaging and modeling-based biomarkers. An expert panel was formed including clinicians, radiologists, imaging scientists, and computational modelers. Section editors generated a series of summary statements based on a review of the literature and professional experience and, following consensus review, a diagnostic algorithm and 55 statements were agreed. The diagnostic algorithm and summary statements emphasize the key role and added value of imaging in the diagnosis and assessment of PH and highlight areas requiring further research.

Zero-dimensional lumped approach to incorporate the dynamic part of the pressure at vessel junctions in a 1D wave propagation model

A benchmark study by Boileau et al tested 6 commonly used numerical schemes for 1D wave propagation, for their ability to capture the main features of pressure, flow, and area waveforms in large arteries. While all numerical schemes showed good agreement in pressure and flow waveforms for smaller arterial networks, the simplified trapezium rule method proposed by Kroon et al showed an overestimation for the systolic pressure of 1% in proximal regions and an underestimation of 3% in distal regions in comparison with the 5 other schemes when using a larger arterial network, published as the ADAN56 model. The authors attributed this difference to the neglection of the dynamic part of the pressure at vessel junctions. Carson et al resolved these differences by proposing 2 methods to implement the dynamic part of the pressure in the simplified trapezium rule method scheme. In the present study, an alternative method is introduced extending the work by Kroon et al. This alternative method consists of a new 0D element, which is placed at vessel junctions. The strength of this new element is the ease of implementation and its flexible coupling with other elements, without introducing additional degrees of freedom or the need of a penalty function. This new approach is compared with 5 other numerical schemes, which already have the dynamic part of the pressure incorporated. The new method shows excellent agreement with these schemes for the ADAN56 model.

An efficient parallel simulation of unsteady blood flows in patient-specific pulmonary artery

Simulation of blood flows in the pulmonary artery provides some insight into certain diseases by examining the relationship between some continuum metrics, eg, the wall shear stress acting on the vascular endothelium, which responds to flow-induced mechanical forces by releasing vasodilators/constrictors. V. Kheyfets, in his previous work, studies numerically a patient-specific pulmonary circulation to show that decreasing wall shear stress is correlated with increasing pulmonary vascular impedance. In this paper, we develop a scalable parallel algorithm based on domain decomposition methods to investigate an unsteady model with patient-specific pulsatile waveforms as the inlet boundary condition. The unsteady model offers tremendously more information about the dynamic behavior of the flow field, but computationally speaking, the simulation is a lot more expensive since a problem which is similar to the steady-state problem has to be solved many times, and therefore, the traditional sequential approach is not suitable anymore. We show computationally that simulations using the proposed parallel approach with up to 10 000 processor cores can be obtained with much reduced compute time. This makes the technology potentially usable for the routine study of the dynamic behavior of blood flows in the pulmonary artery, in particular, the changes of the blood flows and the wall shear stress in the spatial and temporal dimensions.

A patient-specific lumped-parameter model of coronary circulation

A new lumped-parameter model for coronary hemodynamics is developed. This model is developed for the whole coronary network based on CT scans of a patient-specific geometry including the right coronary tree, which is absent in many previous mathematical models. The model adopts the structured tree model boundary conditions similar to the work of Olufsen et al., thus avoiding the necessity of invasive perfusion measurements. In addition, we also incorporated the effects of the head loss at the two inlets of the large coronary arteries for the first time. The head loss could explain the phenomenon of a sudden increase of the resistance at the inlet of coronary vessel. The estimated blood pressure and flow rate results from the model agree well with the clinical measurements. The computed impedances also match the experimental perfusion measurement. The effects of coronary arterial stenosis are considered and the fractional flow reserve and relative flow in the coronary vessels for a stenotic vessel computed in this model show good agreement with published experimental data. It is believed that the approach could be readily translated to clinical practice to facilitate real time clinical diagnosis.

Determinant Factors for Arterial Hemodynamics in Hypertension: Theoretical Insights From a Computational Model-Based Study

Hypertension is a well-documented predictive factor for cardiovascular events. Clinical studies have extensively demonstrated the differential hemodynamic consequences of various antihypertensive drugs, but failed to clearly elucidate the underlying mechanisms due to the difficulty in performing a quantitative deterministic analysis based on clinical data that carry confounding information stemming from interpatient differences and the nonlinearity of cardiovascular hemodynamics. In the present study, a multiscale model of the cardiovascular system was developed to quantitatively investigate the relationships between hemodynamic variables and cardiovascular properties under hypertensive conditions, aiming to establish a theoretical basis for assisting in the interpretation of clinical observations or optimization of therapy. Results demonstrated that heart period, central arterial stiffness, and arteriolar radius were the major determinant factors for blood pressures and flow pulsatility indices both in large arteries and in the microcirculation. These factors differed in the degree and the way in which they affect hemodynamic variables due to their differential effects on wave reflections in the vascular system. In particular, it was found that the hemodynamic effects of varying arteriolar radius were considerably influenced by the state of central arterial stiffness, and vice versa, which implied the potential of optimizing antihypertensive treatment by selecting proper drugs based on patient-specific cardiovascular conditions. When analyzed in relation to clinical observations, the simulated results provided mechanistic explanations for the beneficial pressure-lowering effects of vasodilators as compared to β-blockers, and highlighted the significance of monitoring and normalizing arterial stiffness in the treatment of hypertension.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

Arterial pulse attenuation prediction using the decaying rate of a pressure wave in a viscoelastic material model

The present study examines the possibility of attenuating blood pulses by means of introducing prosthetic viscoelastic materials able to absorb energy and damp such pulses. Vascular prostheses made of polymeric materials modify the mechanical properties of blood vessels. The effect of these materials on the blood pulse propagation remains to be fully understood. Several materials for medical applications, such as medical polydimethylsiloxane or polytetrafluoroethylene, show viscoelastic behavior, modifying the original vessel stiffness and affecting the propagation of blood pulses. This study focuses on the propagation of pressure waves along a pipe with viscoelastic materials using the Maxwell and the Zener models. An expression of exponential decay has been obtained for the Maxwell material model and also for low viscous coefficient values in the Zener model. For relatively high values of the viscous term in the Zener model, the steepest part of the pulse can be damped quickly, leaving a smooth, slowly decaying wave. These mathematical models are critical to tailor those materials used in cardiovascular implants to the mechanical environment they are confronted with to repair or improve blood vessel function.

SIMULATION ANALYSIS OF BLOOD FLOW IN ARTERIES OF THE HUMAN ARM

Arteries in the upper limb play important roles in the circulation system of the human body. In particular, the radial artery has received considerable attention in traditional Chinese medicine for thousands of years. Here, a 3D model for the arm arteries has been created uncomplicated, in a Chinese adult's left hand, from the magnetic resonance imaging data, using professional modeling software to restore the basic structure of the arm artery in human body, before being imported to Ansys software for simulation. Blood model has been only simulated, and using the blood density of constant parameter and viscosity using the Carreau fluid model, and using viscous-laminar model of Fluent to obtain the velocity profile, static pressure and shear stress in the brachial, interosseous, ulnar, radial and palmar arch arteries. In particular, the brachial and bifurcations have the high pressure and velocity profiles. The simulation results obtained here are also validated by those published in the literature and proved the ulnar artery prevails over the radial artery as a blood supplier to the vessels in the wrist and hand.

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.

Flow optimization in vascular networks

The development of mathematical models for studying phenomena observed in vascular networks is very useful for its potential applications in medicine and physiology. Detailed 3D studies of flow in the arterial system based on the Navier-Stokes equations require high computational power, hence reduced models are often used, both for the constitutive laws and the spatial domain. In order to capture the major features of the phenomena under study, such as variations in arterial pressure and flow velocity, the resulting PDE models on networks require appropriate junction and boundary conditions. Instead of considering an entire network, we simulate portions of the latter and use inflow and outflow conditions which realistically mimic the behavior of the network that has not been included in the spatial domain. The resulting PDEs are solved numerically using a discontinuous Galerkin scheme for the spatial and Adam-Bashforth method for the temporal discretization. The aim is to study the effect of truncation to the flow in the root edge of a fractal network, the effect of adding or subtracting an edge to a given network, and optimal control strategies on a network in the event of a blockage or unblockage of an edge or of an entire subtree.