skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: 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, e.g., 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.  more » « less
Award ID(s):
1720366
PAR ID:
10057874
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
International journal for numerical methods in biomedical engineering
ISSN:
2040-7939
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Pulmonary arterial hypertension (PAH) is a vasculopathy characterized by sustained elevated pulmonary arterial pressures in which the pulmonary vasculature undergoes significant structural and functional remodeling. To better understand disease mechanisms, in this review article we highlight recent insights into the regulation of pulmonary arterial cells by mechanical cues associated with PAH. Specifically, the mechanobiology of pulmonary arterial endothelial cells (PAECs), smooth muscle cells (PASMCs) and adventitial fibroblasts (PAAFs) has been investigated in vivo, in vitro, and in silico. Increased pulmonary arterial pressure increases vessel wall stress and strain and endothelial fluid shear stress. These mechanical cues promote vasoconstriction and fibrosis that contribute further to hypertension and alter the mechanical milieu and regulation of pulmonary arterial cells. 
    more » « less
  2. AbstractOne‐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 inin vivoimaging 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.image Key pointsThis 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. 
    more » « less
  3. 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. 
    more » « less
  4. IntroductionPrimary pulmonary vein stenosis (PVS) is a rare congenital heart disease that proves to be a clinical challenge due to the rapidly progressive disease course and high rates of treatment complications. PVS intervention is frequently faced with in-stent restenosis and persistent disease progression despite initial venous recanalization with balloon angioplasty or stenting. Alterations in wall shear stress (WSS) have been previously associated with neointimal hyperplasia and venous stenosis underlying PVS progression. Thus, the development of patient-specific three-dimensional (3D)in vitromodels is needed to further investigate the biomechanical outcomes of endovascular and surgical interventions. MethodsIn this study, deidentified computed tomography images from three patients were segmented to generate perfusable phantom models of pulmonary veins before and after catheterization. These 3D reconstructions were 3D printed using a clear resin ink and used in a benchtop experimental setup. Computational fluid dynamic (CFD) analysis was performed on modelsin silicoutilizing Doppler echocardiography data to represent thein vivoflow conditions at the inlets. Particle image velocimetry was conducted using the benchtop perfusion setup to analyze WSS and velocity profiles and the results were compared with those predicted by the CFD model. ResultsOur findings indicated areas of undesirable alterations in WSS before and after catheterization, in comparison with the published baseline levels in the healthyin vivotissues that may lead to regional disease progression. DiscussionThe established patient-specific 3Din vitromodels and the developedin vitro–in silicoplatform demonstrate great promise to refine interventional approaches and mitigate complications in treating patients with primary PVS. 
    more » « less
  5. Carotid artery diseases, such as atherosclerosis, are a major cause of death in the United States. Wall shear stresses are known to prompt plaque formation, but there is limited understanding of the complex flow structures underlying these stresses and how they differ in a pre-disposed high-risk patient cohort. A ‘healthy’ and a novel ‘pre-disposed’ carotid artery bifurcation model was determined based on patient-averaged clinical data, where the ‘pre-disposed’ model represents a pathological anatomy. Computational fluid dynamic simulations were performed using a physiological flow based on healthy human subjects. A main hairpin vortical structure in the internal carotid artery sinus was observed, which locally increased instantaneous wall shear stress. In the pre-disposed geometry, this vortical structure starts at an earlier instance in the cardiac flow cycle and persists over a much shorter period, where the second half of the cardiac cycle is dominated by perturbed secondary flow structures and vortices. This coincides with weaker favorable axial pressure gradient peaks over the sinus for the ‘pre-disposed’ geometry. The findings reveal a strong correlation between vortical structures and wall shear stress and imply that an intact internal carotid artery sinus hairpin vortical structure has a physiologically beneficial role by increasing local wall shear stresses. The deterioration of this beneficial vortical structure is expected to play a significant role in atherosclerotic plaque formation. 
    more » « less