DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Abstract

Abstract 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 flowsmore » and the wall shear stress in the spatial and temporal dimensions.« less

Authors:
 [1];  [2];  [3]; ORCiD logo [4]
  1. Idaho National Lab. (INL), Idaho Falls, ID (United States). Modeling and Simulation
  2. Univ. of Colorado, Denver, CO (United States). School of Medicine
  3. Univ. of Texas, San Antonio, TX (United States). Dept. of Mechanical Engineering
  4. Univ. of Colorado, Boulder, CO (United States). Dept. of Computer Science
Publication Date:
Research Org.:
Idaho National Laboratory (INL), Idaho Falls, ID (United States)
Sponsoring Org.:
USDOE; National Science Foundation (NSF)
OSTI Identifier:
1484705
Alternate Identifier(s):
OSTI ID: 1419356
Report Number(s):
INL/JOU-17-42412-Rev001
Journal ID: ISSN 2040-7939
Grant/Contract Number:  
AC07-05ID14517; DMS-1720366
Resource Type:
Accepted Manuscript
Journal Name:
International Journal for Numerical Methods in Biomedical Engineering
Additional Journal Information:
Journal Volume: 34; Journal Issue: 4; Journal ID: ISSN 2040-7939
Publisher:
Wiley
Country of Publication:
United States
Language:
English
Subject:
60 APPLIED LIFE SCIENCES; unsteady blood flows; patient-specific pulmonary artery; finite element; domain decomposition; parallel processing

Citation Formats

Kong, Fande, Kheyfets, Vitaly, Finol, Ender, and Cai, Xiao-Chuan. An efficient parallel simulation of unsteady blood flows in patient-specific pulmonary artery. United States: N. p., 2017. Web. doi:10.1002/cnm.2952.
Kong, Fande, Kheyfets, Vitaly, Finol, Ender, & Cai, Xiao-Chuan. An efficient parallel simulation of unsteady blood flows in patient-specific pulmonary artery. United States. https://doi.org/10.1002/cnm.2952
Kong, Fande, Kheyfets, Vitaly, Finol, Ender, and Cai, Xiao-Chuan. Fri . "An efficient parallel simulation of unsteady blood flows in patient-specific pulmonary artery". United States. https://doi.org/10.1002/cnm.2952. https://www.osti.gov/servlets/purl/1484705.
@article{osti_1484705,
title = {An efficient parallel simulation of unsteady blood flows in patient-specific pulmonary artery},
author = {Kong, Fande and Kheyfets, Vitaly and Finol, Ender and Cai, Xiao-Chuan},
abstractNote = {Abstract 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.},
doi = {10.1002/cnm.2952},
journal = {International Journal for Numerical Methods in Biomedical Engineering},
number = 4,
volume = 34,
place = {United States},
year = {Fri Dec 15 00:00:00 EST 2017},
month = {Fri Dec 15 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 16 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Computational modelling of 1D blood flow with variable mechanical properties and its application to the simulation of wave propagation in the human arterial system: COMPUTATIONAL MODELLING OF 1D BLOOD FLOW
journal, October 2003

  • Sherwin, S. J.; Formaggia, L.; Peiró, J.
  • International Journal for Numerical Methods in Fluids, Vol. 43, Issue 6-7
  • DOI: 10.1002/fld.543

Splitting Methods Based on Algebraic Factorization for Fluid-Structure Interaction
journal, January 2008

  • Badia, Santiago; Quaini, Annalisa; Quarteroni, Alfio
  • SIAM Journal on Scientific Computing, Vol. 30, Issue 4
  • DOI: 10.1137/070680497

Semi-implicit numerical modeling of axially symmetric flows in compliant arterial systems: AXIALLY SYMMETRIC FLOWS IN COMPLIANT ARTERIES
journal, September 2011

  • Casulli, Vincenzo; Dumbser, Michael; Toro, Eleuterio F.
  • International Journal for Numerical Methods in Biomedical Engineering, Vol. 28, Issue 2
  • DOI: 10.1002/cnm.1464

Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation
journal, March 2014

  • Qureshi, M. Umar; Vaughan, Gareth D. A.; Sainsbury, Christopher
  • Biomechanics and Modeling in Mechanobiology, Vol. 13, Issue 5
  • DOI: 10.1007/s10237-014-0563-y

A stabilized finite element method for the incompressible Navier-Stokes equations using a hierarchical basis
journal, January 2001


A Restricted Additive Schwarz Preconditioner for General Sparse Linear Systems
journal, January 1999


A Fast Implementation of the Minimum Degree Algorithm Using Quotient Graphs
journal, September 1980

  • George, Alan; Liu, Joseph W. H.
  • ACM Transactions on Mathematical Software, Vol. 6, Issue 3
  • DOI: 10.1145/355900.355906

An Automatic One-Way Dissection Algorithm for Irregular Finite Element Problems
journal, December 1980

  • George, Alan
  • SIAM Journal on Numerical Analysis, Vol. 17, Issue 6
  • DOI: 10.1137/0717062

A Highly Scalable Multilevel Schwarz Method with Boundary Geometry Preserving Coarse Spaces for 3D Elasticity Problems on Domains with Complex Geometry
journal, January 2016

  • Kong, Fande; Cai, Xiao-Chuan
  • SIAM Journal on Scientific Computing, Vol. 38, Issue 2
  • DOI: 10.1137/15M1010567

Use of computational fluid dynamics in the design of surgical procedures: Application to the study of competitive flows in cavopulmonary connections
journal, March 1996

  • de Leval, M. R.; Dubini, G.; Migliavacca c§, F.
  • The Journal of Thoracic and Cardiovascular Surgery, Vol. 111, Issue 3
  • DOI: 10.1016/S0022-5223(96)70302-1

Well-balanced high-order numerical schemes for one-dimensional blood flow in vessels with varying mechanical properties
journal, June 2013

  • Müller, Lucas O.; Parés, Carlos; Toro, Eleuterio F.
  • Journal of Computational Physics, Vol. 242
  • DOI: 10.1016/j.jcp.2013.01.050

Nested Dissection of a Regular Finite Element Mesh
journal, April 1973

  • George, Alan
  • SIAM Journal on Numerical Analysis, Vol. 10, Issue 2
  • DOI: 10.1137/0710032

A Scalable Newton–Krylov–Schwarz Method for the Bidomain Reaction-Diffusion System
journal, January 2009

  • Munteanu, M.; Pavarino, L. F.; Scacchi, S.
  • SIAM Journal on Scientific Computing, Vol. 31, Issue 5
  • DOI: 10.1137/08074355X

Patient-specific computational modeling of blood flow in the pulmonary arterial circulation
journal, July 2015

  • Kheyfets, Vitaly O.; Rios, Lourdes; Smith, Triston
  • Computer Methods and Programs in Biomedicine, Vol. 120, Issue 2
  • DOI: 10.1016/j.cmpb.2015.04.005

Impact of pulmonary vascular stiffness and vasodilator treatment in pediatric pulmonary hypertension: 21 patient-specific fluid–structure interaction studies
journal, November 2012

  • Su, Zhenbi; Hunter, Kendall S.; Shandas, Robin
  • Computer Methods and Programs in Biomedicine, Vol. 108, Issue 2
  • DOI: 10.1016/j.cmpb.2011.09.002

Fluid-structure interaction in a pulmonary arterial bifurcation
journal, January 2007


Globally Convergent Inexact Newton Methods
journal, May 1994

  • Eisenstat, Stanley C.; Walker, Homer F.
  • SIAM Journal on Optimization, Vol. 4, Issue 2
  • DOI: 10.1137/0804022

Simulations of Congenital Septal Defect Closure and Reactivity Testing in Patient-Specific Models of the Pediatric Pulmonary Vasculature: A 3D Numerical Study With Fluid-Structure Interaction
journal, January 2006

  • Hunter, Kendall S.; Lanning, Craig J.; Chen, Shiuh-Yung J.
  • Journal of Biomechanical Engineering, Vol. 128, Issue 4
  • DOI: 10.1115/1.2206202

Patient-Specific Modeling of Cardiovascular Mechanics
journal, August 2009


Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries
journal, February 2007

  • Spilker, Ryan L.; Feinstein, Jeffrey A.; Parker, David W.
  • Annals of Biomedical Engineering, Vol. 35, Issue 4
  • DOI: 10.1007/s10439-006-9240-3

Three-Dimensional Hemodynamics in the Human Pulmonary Arteries Under Resting and Exercise Conditions
journal, July 2010

  • Tang, Beverly T.; Fonte, Tim A.; Chan, Frandics P.
  • Annals of Biomedical Engineering, Vol. 39, Issue 1
  • DOI: 10.1007/s10439-010-0124-1

Finite element modeling of blood flow in arteries
journal, May 1998

  • Taylor, Charles A.; Hughes, Thomas J. R.; Zarins, Christopher K.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 158, Issue 1-2
  • DOI: 10.1016/S0045-7825(98)80008-X

A two-dimensional effective model describing fluid–structure interaction in blood flow: analysis, simulation and experimental validation
journal, December 2005

  • Čanić, Sunčica; Mikelić, Andro; Tambača, Josip
  • Comptes Rendus Mécanique, Vol. 333, Issue 12
  • DOI: 10.1016/j.crme.2005.10.005

Parallel Newton--Krylov--Schwarz Algorithms for the Transonic Full Potential Equation
journal, January 1998

  • Cai, Xiao-Chuan; Gropp, William D.; Keyes, David E.
  • SIAM Journal on Scientific Computing, Vol. 19, Issue 1
  • DOI: 10.1137/S1064827596304046

Works referencing / citing this record:

Simulation of unsteady blood flows in a patient‐specific compliant pulmonary artery with a highly parallel monolithically coupled fluid‐structure interaction algorithm
journal, May 2019

  • Kong, Fande; Kheyfets, Vitaly; Finol, Ender
  • International Journal for Numerical Methods in Biomedical Engineering, Vol. 35, Issue 7
  • DOI: 10.1002/cnm.3208