skip to main content

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.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
International journal for numerical methods in biomedical engineering
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Numerical simulations for computational hemodynamics in clinical settings require a combination of many ingredients, mathematical models, solvers and patient-specific data. The sensitivity of the solutions to these factors may be critical, particularly when we have a partial or noisy knowledge of data. Uncertainty quantification is crucial to assess the reliability of the results. We present here an extensive sensitivity analysis in aortic flow simulations, to quantify the dependence of clinically relevant quantities to the patient-specific geometry and the inflow boundary conditions. Geometry and inflow conditions are generally believed to have a major impact on numerical simulations. We resort tomore »a global sensitivity analysis, (i.e., not restricted to a linearization around a working point), based on polynomial chaos expansion (PCE) and the associated Sobol' indices. We regard the geometry and the inflow conditions as the realization of a parametric stochastic process. To construct a physically consistent stochastic process for the geometry, we use a set of longitudinal-in-time images of a patient with an abdominal aortic aneurysm (AAA) to parametrize geometrical variations. Aortic flow is highly disturbed during systole. This leads to high computational costs, even amplified in a sensitivity analysis -when many simulations are needed. To mitigate this, we consider here a large Eddy simulation (LES) model. Our model depends in particular on a user-defined parameter called filter radius. We borrowed the tools of the global sensitivity analysis to assess the sensitivity of the solution to this parameter too. The targeted quantities of interest (QoI) include: the total kinetic energy (TKE), the time-average wall shear stress (TAWSS), and the oscillatory shear index (OSI). The results show that these indexes are mostly sensitive to the geometry. Also, we find that the sensitivity may be different during different instants of the heartbeat and in different regions of the domain of interest. This analysis helps to assess the reliability of in silico tools for clinical applications.« less
  2. 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 inmore »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.« less
  3. 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 andmore »fibrosis that contribute further to hypertension and alter the mechanical milieu and regulation of pulmonary arterial cells.« less
  4. Advances in mechanobiology have evolved through insights from multiple disciplines including structural engineering, biomechanics, vascular biology, and orthopaedics. In this paper, we reviewed the impact of key reports related to the study of applied loads on tissues and cells and the resulting signal transduction pathways. We addressed how technology has helped advance the burgeoning field of mechanobiology (over 33,600 publications from 1970 to 2016). We analyzed the impact of critical ideas and then determined how these concepts influenced the mechanobiology field by looking at the citation frequency of these reports as well as tracking how the overall number of citationsmore »within the field changed over time. These data allowed us to understand how a key publication, idea, or technology guided or enabled the field. Initial observations of how forces acted on bone and soft tissues stimulated the development of computational solutions defining how forces affect tissue modeling and remodeling. Enabling technologies, such as cell and tissue stretching, compression, and shear stress devices, allowed more researchers to explore how deformation and fluid flow affect cells. Observation of the cell as a tensegrity structure and advanced methods to study genetic regulation in cells further advanced knowledge of specific mechanisms of mechanotransduction. The future of the field will involve developing gene and drug therapies to simulate or augment beneficial load regimens in patients and in mechanically conditioning organs for implantation. Here, we addressed a history of the field, but we limited our discussions to advances in musculoskeletal mechanobiology, primarily in bone, tendon, and ligament tissues.« less
  5. Cerebral aneurysm clips are biomedical implants applied by neurosurgeons to re-approximate arterial vessel walls and prevent catastrophic aneurysmal hemorrhages in patients. Current methods of aneurysm clip production are labor intensive and time-consuming, leading to high costs per implant and limited variability in clip morphology. Metal additive manufacturing is investigated as an alternative to traditional manufacturing methods that may enable production of patient-specific aneurysm clips to account for variations in individual vascular anatomy and possibly reduce surgical complication risks. Relevant challenges to metal additive manufacturing are investigated for biomedical implants, including material choice, design limitations, postprocessing, printed material properties, and combinedmore »production methods. Initial experiments with additive manufacturing of 316 L stainless steel aneurysm clips are carried out on a selective laser melting (SLM) system. The dimensions of the printed clips were found to be within 0.5% of the dimensions of the designed clips. Hardness and density of the printed clips (213 ± 7 HV1 and 7.9 g/cc, respectively) were very close to reported values for 316 L stainless steel, as expected. No ferrite and minimal porosity is observed in a cross section of a printed clip, with some anisotropy in the grain orientation. A clamping force of approximately 1 N is measured with a clip separation of 1.5 mm. Metal additive manufacturing shows promise for use in the creation of custom aneurysm clips, but some of the challenges discussed will need to be addressed before clinical use is possible.« less