Well‐balanced high‐order solver for blood flow in networks of vessels with variable properties

We present a well-balanced, high-order non-linear numerical scheme for solving a hyperbolic system that models one-dimensional flow in blood vessels with variable mechanical and geometrical properties along their length. Using a suitable set of test problems with exact solution, we rigorously assess the performance of the scheme. In particular, we assess the well-balanced property and the effective order of accuracy through an empirical convergence rate study. Schemes of up to fifth order of accuracy in both space and time are implemented and assessed. The numerical methodology is then extended to realistic networks of elastic vessels and is validated against published state-of-the-art numerical solutions and experimental measurements. It is envisaged that the present scheme will constitute the building block for a closed, global model for the human circulation system involving arteries, veins, capillaries and cerebrospinal fluid.

[1]  Eleuterio F. Toro,et al.  Exact solution of the Riemann problem for the shallow water equations with discontinuous bottom geometry , 2008, J. Comput. Phys..

[2]  M. J. Castro,et al.  ADER schemes on unstructured meshes for nonconservative hyperbolic systems: Applications to geophysical flows , 2009 .

[3]  A. Quarteroni,et al.  One-dimensional models for blood flow in arteries , 2003 .

[4]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[5]  Laurent Dumas,et al.  An optimal reconstruction of the human arterial tree from doppler echotracking measurements , 2012, GECCO '12.

[6]  Michael Dumbser,et al.  Well-Balanced High-Order Centred Schemes for Non-Conservative Hyperbolic Systems. Applications to Shallow Water Equations with Fixed and Mobile Bed , 2009 .

[7]  P. Blanco,et al.  On the potentialities of 3D-1D coupled models in hemodynamics simulations. , 2009, Journal of biomechanics.

[8]  E. Attinger,et al.  Wall properties of veins. , 1969, IEEE transactions on bio-medical engineering.

[9]  Sunčica Čanić,et al.  Blood flow through compliant vessels after endovascular repair: wall deformations induced by the discontinuous wall properties , 2002 .

[10]  Michael Dumbser,et al.  ADER Schemes for Nonlinear Systems of Stiff Advection–Diffusion–Reaction Equations , 2011, J. Sci. Comput..

[11]  T J Pedley,et al.  Blood pressure and flow rate in the giraffe jugular vein. , 1996, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[12]  Michael Dumbser,et al.  A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems , 2011, J. Sci. Comput..

[13]  Timothy J. Pedley,et al.  Numerical solutions for unsteady gravity-driven flows in collapsible tubes: evolution and roll-wave instability of a steady state , 1999, Journal of Fluid Mechanics.

[14]  Eleuterio F. Toro,et al.  Flow in Collapsible Tubes with Discontinuous Mechanical Properties: Mathematical Model and Exact Solutions , 2013 .

[15]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[16]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[17]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[18]  Olivier Delestre,et al.  A ‘well‐balanced’ finite volume scheme for blood flow simulation , 2012, ArXiv.

[19]  C. Parés.  Numerical methods for nonconservative hyperbolic systems: a theoretical framework. , 2006 .

[20]  Michael Dumbser,et al.  Semi-implicit numerical modeling of axially symmetric flows in compliant arterial systems. , 2012, International journal for numerical methods in biomedical engineering.

[21]  S. Sherwin,et al.  Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows. , 2007, Journal of biomechanics.

[22]  M. Vázquez-Cendón.  Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .

[23]  F. Liang,et al.  A computational model study of the influence of the anatomy of the circle of willis on cerebral hyperperfusion following carotid artery surgery , 2011, Biomedical engineering online.

[24]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[25]  S. Sherwin,et al.  One-dimensional modelling of a vascular network in space-time variables , 2003 .

[26]  S. Sherwin,et al.  Lumped parameter outflow models for 1-D blood flow simulations: Effect on pulse waves and parameter estimation , 2008 .

[27]  Manuel Jesús Castro Díaz,et al.  Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes , 2008, J. Comput. Phys..

[28]  Philip L. Roe,et al.  Upwind differencing schemes for hyperbolic conservation laws with source terms , 1987 .

[29]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[30]  Michael Dumbser,et al.  Comparison of solvers for the generalized Riemann problem for hyperbolic systems with source terms , 2012, J. Comput. Phys..

[31]  Carlos Parés,et al.  On the well-balance property of Roe?s method for nonconservative hyperbolic systems , 2004 .

[32]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[33]  J. Karemaker,et al.  Human cerebral venous outflow pathway depends on posture and central venous pressure , 2004, The Journal of physiology.

[34]  Michael Dumbser,et al.  On Universal Osher-Type Schemes for General Nonlinear Hyperbolic Conservation Laws , 2011 .

[35]  Spencer J. Sherwin,et al.  Computational modelling of 1D blood flow with variable mechanical properties and its application to the simulation of wave propagation in the human arterial system , 2003 .

[36]  Lucas O. Müller,et al.  Well-balanced high-order numerical schemes for one-dimensional blood flow in vessels with varying mechanical properties , 2013, J. Comput. Phys..

[37]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[38]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements , 2011, Journal of biomechanics.

[39]  Eleuterio F. Toro,et al.  ADER scheme on unstructured meshes for shallow water: simulation of tsunami waves , 2012 .

[40]  Carlos Parés Madroñal,et al.  On the Convergence and Well-Balanced Property of Path-Conservative Numerical Schemes for Systems of Balance Laws , 2011, J. Sci. Comput..