Low-Shapiro hydrostatic reconstruction technique for blood flow simulation in large arteries with varying geometrical and mechanical properties

JOURNAL OF COMPUTATIONAL PHYSICS

卷:331

Article

Ghigo, A. R.^1,2  Delestre, O.^3,4  Fullana, J. -M.^1,2  Lagree, P. -Y.^1,2 

[1] Sorbonne Univ, CNRS, 4 Pl Jussieu,Boite 162, F-75005 Paris, France

[2] UPMC Univ Paris 06, UMR 7190, Inst Jean Le Rond dAlembert, 4 Pl Jussieu,Boite 162, F-75005 Paris, France

[3] Univ Nice Sophia Antipolis UNS, UMR CNRS 7351, Lab JA Dieudonne, Parc Valrose, F-06108 Nice 02, France

[4] Univ Nice Sophia Antipolis UNS, Polytech Nice Sophia, Parc Valrose, F-06108 Nice 02, France

关键词: Blood flow;    1D model;    Finite volume;    Well-balanced;    Hydrostatic reconstruction;   

DOI  :  10.1016/j.jcp.2016.11.032

来源: Elsevier

【 摘 要 】

The purpose of this work is to construct a simple, efficient and accurate well-balanced numerical scheme for one-dimensional (1D) blood flow in large arteries with varying geometrical and mechanical properties. As the steady states at rest are not relevant for blood flow, we construct two well-balanced hydrostatic reconstruction techniques designed to, preserve low-Shapiro number steady states that may occur in large network simulations. The Shapiro number S-h = mu/C is the equivalent of the Froude number for shallow water equations and the Mach number for compressible Euler equations. The first is the low Shapiro hydrostatic reconstruction (HR-LS), which is a simple and efficient method, inspired from the hydrostatic reconstruction technique (HR). The second is the subsonic hydrostatic reconstruction (HR-S), adapted here to blood flow and designed to exactly preserve all subcritical steady states. We systematically compare HR, HR-LS and HR-S in a series of single artery and arterial network numerical tests designed to evaluate their well-balanced and wave-capturing properties. The results indicate that HR is not adapted to compute blood flow in large arteries as it is unable to capture wave reflections and transmissions when large variations of the arteries' geometrical and mechanical properties are considered. On the contrary, HR-S is exactly well-balanced and is the most accurate hydrostatic reconstruction technique. However, HR-LS is able to compute low-Shapiro number steady states as well as wave reflections and transmissions with satisfying accuracy and is simpler and computationally less expensive than HR-S. We therefore recommend using HR-LS for 1D blood flow simulations in large arterial network simulations. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free