A fast algorithm for the simulation of arterial pulse waves

Tao Du, Dan Hu, David Cai

Research output: Contribution to journalArticle

Abstract

One-dimensional models have been widely used in studies of the propagation of blood pulse waves in large arterial trees. Under a periodic driving of the heartbeat, traditional numerical methods, such as the Lax-Wendroff method, are employed to obtain asymptotic periodic solutions at large times. However, these methods are severely constrained by the CFL condition due to large pulse wave speed. In this work, we develop a new numerical algorithm to overcome this constraint. First, we reformulate the model system of pulse wave propagation using a set of Riemann variables and derive a new form of boundary conditions at the inlet, the outlets, and the bifurcation points of the arterial tree. The new form of the boundary conditions enables us to design a convergent iterative method to enforce the boundary conditions. Then, after exchanging the spatial and temporal coordinates of the model system, we apply the Lax-Wendroff method in the exchanged coordinate system, which turns the large pulse wave speed from a liability to a benefit, to solve the wave equation in each artery of the model arterial system. Our numerical studies show that our new algorithm is stable and can perform ~15 times faster than the traditional implementation of the Lax-Wendroff method under the requirement that the relative numerical error of blood pressure be smaller than one percent, which is much smaller than the modeling error.

Original languageEnglish (US)
Pages (from-to)450-464
Number of pages15
JournalJournal of Computational Physics
Volume314
DOIs
StatePublished - Jun 1 2016

Fingerprint

Boundary conditions
boundary conditions
pulses
simulation
liabilities
blood pressure
exchanging
Blood pressure
Wave equations
outlets
arteries
Iterative methods
Wave propagation
wave equations
blood
wave propagation
Numerical methods
Blood
requirements
propagation

Keywords

  • Blood pulse wave
  • Fast algorithm
  • Large wave speed
  • Lax-Wendroff method
  • Riemann variables

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

A fast algorithm for the simulation of arterial pulse waves. / Du, Tao; Hu, Dan; Cai, David.

In: Journal of Computational Physics, Vol. 314, 01.06.2016, p. 450-464.

Research output: Contribution to journalArticle

Du, Tao ; Hu, Dan ; Cai, David. / A fast algorithm for the simulation of arterial pulse waves. In: Journal of Computational Physics. 2016 ; Vol. 314. pp. 450-464.
@article{abdc85c08ba840e4b2a86633d53e35a5,
title = "A fast algorithm for the simulation of arterial pulse waves",
abstract = "One-dimensional models have been widely used in studies of the propagation of blood pulse waves in large arterial trees. Under a periodic driving of the heartbeat, traditional numerical methods, such as the Lax-Wendroff method, are employed to obtain asymptotic periodic solutions at large times. However, these methods are severely constrained by the CFL condition due to large pulse wave speed. In this work, we develop a new numerical algorithm to overcome this constraint. First, we reformulate the model system of pulse wave propagation using a set of Riemann variables and derive a new form of boundary conditions at the inlet, the outlets, and the bifurcation points of the arterial tree. The new form of the boundary conditions enables us to design a convergent iterative method to enforce the boundary conditions. Then, after exchanging the spatial and temporal coordinates of the model system, we apply the Lax-Wendroff method in the exchanged coordinate system, which turns the large pulse wave speed from a liability to a benefit, to solve the wave equation in each artery of the model arterial system. Our numerical studies show that our new algorithm is stable and can perform ~15 times faster than the traditional implementation of the Lax-Wendroff method under the requirement that the relative numerical error of blood pressure be smaller than one percent, which is much smaller than the modeling error.",
keywords = "Blood pulse wave, Fast algorithm, Large wave speed, Lax-Wendroff method, Riemann variables",
author = "Tao Du and Dan Hu and David Cai",
year = "2016",
month = "6",
day = "1",
doi = "10.1016/j.jcp.2016.03.036",
language = "English (US)",
volume = "314",
pages = "450--464",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - A fast algorithm for the simulation of arterial pulse waves

AU - Du, Tao

AU - Hu, Dan

AU - Cai, David

PY - 2016/6/1

Y1 - 2016/6/1

N2 - One-dimensional models have been widely used in studies of the propagation of blood pulse waves in large arterial trees. Under a periodic driving of the heartbeat, traditional numerical methods, such as the Lax-Wendroff method, are employed to obtain asymptotic periodic solutions at large times. However, these methods are severely constrained by the CFL condition due to large pulse wave speed. In this work, we develop a new numerical algorithm to overcome this constraint. First, we reformulate the model system of pulse wave propagation using a set of Riemann variables and derive a new form of boundary conditions at the inlet, the outlets, and the bifurcation points of the arterial tree. The new form of the boundary conditions enables us to design a convergent iterative method to enforce the boundary conditions. Then, after exchanging the spatial and temporal coordinates of the model system, we apply the Lax-Wendroff method in the exchanged coordinate system, which turns the large pulse wave speed from a liability to a benefit, to solve the wave equation in each artery of the model arterial system. Our numerical studies show that our new algorithm is stable and can perform ~15 times faster than the traditional implementation of the Lax-Wendroff method under the requirement that the relative numerical error of blood pressure be smaller than one percent, which is much smaller than the modeling error.

AB - One-dimensional models have been widely used in studies of the propagation of blood pulse waves in large arterial trees. Under a periodic driving of the heartbeat, traditional numerical methods, such as the Lax-Wendroff method, are employed to obtain asymptotic periodic solutions at large times. However, these methods are severely constrained by the CFL condition due to large pulse wave speed. In this work, we develop a new numerical algorithm to overcome this constraint. First, we reformulate the model system of pulse wave propagation using a set of Riemann variables and derive a new form of boundary conditions at the inlet, the outlets, and the bifurcation points of the arterial tree. The new form of the boundary conditions enables us to design a convergent iterative method to enforce the boundary conditions. Then, after exchanging the spatial and temporal coordinates of the model system, we apply the Lax-Wendroff method in the exchanged coordinate system, which turns the large pulse wave speed from a liability to a benefit, to solve the wave equation in each artery of the model arterial system. Our numerical studies show that our new algorithm is stable and can perform ~15 times faster than the traditional implementation of the Lax-Wendroff method under the requirement that the relative numerical error of blood pressure be smaller than one percent, which is much smaller than the modeling error.

KW - Blood pulse wave

KW - Fast algorithm

KW - Large wave speed

KW - Lax-Wendroff method

KW - Riemann variables

UR - http://www.scopus.com/inward/record.url?scp=84961785324&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84961785324&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2016.03.036

DO - 10.1016/j.jcp.2016.03.036

M3 - Article

VL - 314

SP - 450

EP - 464

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -