A high‐order local time stepping finite volume solver for one‐dimensional blood flow simulations: application to the ADAN model

In recent years, the complexity of vessel networks for one-dimensional blood flow models has significantly increased, because of enhanced anatomical detail or automatic peripheral vasculature generation, for example. This fact, along with the application of these models in uncertainty quantification and parameter estimation poses the need for extremely efficient numerical solvers. The aim of this work is to present a finite volume solver for one-dimensional blood flow simulations in networks of elastic and viscoelastic vessels, featuring high-order space-time accuracy and local time stepping (LTS). The solver is built on (i) a high-order finite volume type numerical scheme, (ii) a high-order treatment of the numerical solution at internal vertexes of the network, often called junctions, and (iii) an accurate LTS strategy. The accuracy of the proposed methodology is verified by empirical convergence tests. Then, the resulting LTS scheme is applied to arterial networks of increasing complexity and spatial scale heterogeneity, with a number of one-dimensional segments ranging from a few tens up to several thousands and vessel lengths ranging from less than a millimeter up to tens of centimeters, in order to evaluate its computational cost efficiency. The proposed methodology can be extended to any other hyperbolic system for which network applications are relevant. Copyright © 2016 John Wiley & Sons, Ltd.

[1]  Joaquim Peiró,et al.  Physical determining factors of the arterial pulse waveform: theoretical analysis and calculation using the 1-D formulation , 2012, Journal of Engineering Mathematics.

[2]  Lucas O Müller,et al.  A global multiscale mathematical model for the human circulation with emphasis on the venous system , 2014, International journal for numerical methods in biomedical engineering.

[3]  D. F. Young,et al.  Computer simulation of arterial flow with applications to arterial and aortic stenoses. , 1992, Journal of biomechanics.

[4]  H. Boom,et al.  The use of an analog computer in a circulation model. , 1963, Progress in cardiovascular diseases.

[5]  R. Himeno,et al.  Biomechanical characterization of ventricular-arterial coupling during aging: a multi-scale model study. , 2009, Journal of biomechanics.

[6]  Thomas J. R. Hughes,et al.  On the one-dimensional theory of blood flow in the larger vessels , 1973 .

[7]  Pablo J Blanco,et al.  A black‐box decomposition approach for coupling heterogeneous components in hemodynamics simulations , 2013, International journal for numerical methods in biomedical engineering.

[8]  N. Stergiopulos,et al.  Validation of a one-dimensional model of the systemic arterial tree. , 2009, American journal of physiology. Heart and circulatory physiology.

[9]  Pablo J Blanco,et al.  Blood flow distribution in an anatomically detailed arterial network model: criteria and algorithms , 2014, Biomechanics and modeling in mechanobiology.

[10]  Pablo J Blanco,et al.  On the effect of preload and pre-stretch on hemodynamic simulations: an integrative approach , 2015, Biomechanics and Modeling in Mechanobiology.

[11]  A Noordergraaf,et al.  Analog studies of the human systemic arterial tree. , 1969, Journal of biomechanics.

[12]  Pedro A. Lemos,et al.  An Anatomically Detailed Arterial Network Model for One-Dimensional Computational Hemodynamics , 2015, IEEE Transactions on Biomedical Engineering.

[13]  F. Lazeyras,et al.  Validation of a patient-specific one-dimensional model of the systemic arterial tree. , 2011, American journal of physiology. Heart and circulatory physiology.

[14]  Giulio Antonini,et al.  Rational macromodeling of 1D blood flow in the human cardiovascular system. , 2015, International journal for numerical methods in biomedical engineering.

[15]  A. Quarteroni,et al.  Analysis of lumped parameter models for blood flow simulations and their relation with 1D models , 2004 .

[16]  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 .

[17]  Michael Dumbser,et al.  Arbitrary-Lagrangian–Eulerian ADER–WENO finite volume schemes with time-accurate local time stepping for hyperbolic conservation laws , 2014, 1402.6897.

[18]  Lucas O. Müller,et al.  A high order approximation of hyperbolic conservation laws in networks: Application to one-dimensional blood flow , 2015, J. Comput. Phys..

[19]  A. Quarteroni,et al.  On the coupling of 3D and 1D Navier-Stokes equations for flow problems in compliant vessels , 2001 .

[20]  E. Toro,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes - V. Local time stepping and p-adaptivity , 2007 .

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

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

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

[24]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[25]  Lucas O. Müller,et al.  A benchmark study of numerical schemes for one‐dimensional arterial blood flow modelling , 2015, International journal for numerical methods in biomedical engineering.

[26]  Paolo Crosetto,et al.  Implicit Coupling of One-Dimensional and Three-Dimensional Blood Flow Models with Compliant Vessels , 2013, Multiscale Model. Simul..

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

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

[29]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

[30]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements. , 2007, Journal of biomechanics.

[31]  Eleuterio F. Toro,et al.  Derivative Riemann solvers for systems of conservation laws and ADER methods , 2006, J. Comput. Phys..

[32]  A. P. Avó.  Multi-branched model of the human arterial system , 2006 .

[33]  Lucas O. Müller,et al.  Consistent treatment of viscoelastic effects at junctions in one-dimensional blood flow models , 2016, J. Comput. Phys..

[34]  Lucas O. Müller,et al.  Enhanced global mathematical model for studying cerebral venous blood flow. , 2014, Journal of biomechanics.

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

[36]  Paris Perdikaris,et al.  An Effective Fractal-Tree Closure Model for Simulating Blood Flow in Large Arterial Networks , 2014, Annals of Biomedical Engineering.

[37]  Michael Dumbser,et al.  ADER-WENO finite volume schemes with space-time adaptive mesh refinement , 2012, J. Comput. Phys..

[38]  Ryutaro Himeno,et al.  Multi-scale modeling of the human cardiovascular system with applications to aortic valvular and arterial stenoses , 2009, Medical & Biological Engineering & Computing.

[39]  M. J. Castro,et al.  FORCE schemes on unstructured meshes II: Non-conservative hyperbolic systems , 2010 .

[40]  Pablo J. Blanco,et al.  Mathematical Model of Blood Flow in an Anatomically Detailed Arterial Network of the Arm , 2013 .

[41]  L. Formaggia,et al.  Numerical modeling of 1D arterial networks coupled with a lumped parameters description of the heart , 2006, Computer methods in biomechanics and biomedical engineering.

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

[43]  Armin Iske,et al.  ADER schemes on adaptive triangular meshes for scalar conservation laws , 2005 .

[44]  E. Toro,et al.  Solution of the generalized Riemann problem for advection–reaction equations , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[45]  Lucas O. Müller,et al.  Hyperbolic reformulation of a 1D viscoelastic blood flow model and ADER finite volume schemes , 2014, J. Comput. Phys..

[46]  Eleuterio F. Toro,et al.  Space–time adaptive numerical methods for geophysical applications , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[47]  L. Dagum,et al.  OpenMP: an industry standard API for shared-memory programming , 1998 .

[48]  Eleuterio F. Toro,et al.  Reformulations for general advection-diffusion-reaction equations and locally implicit ADER schemes , 2014, J. Comput. Phys..

[49]  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..

[50]  Arne Taube,et al.  A high-order discontinuous Galerkin method with time-accurate local time stepping for the Maxwell equations , 2009 .

[51]  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.

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

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

[54]  M. Dumbser,et al.  High order space–time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems , 2013, 1304.5408.

[55]  Lucas O. Müller,et al.  Well‐balanced high‐order solver for blood flow in networks of vessels with variable properties , 2013, International journal for numerical methods in biomedical engineering.

[56]  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.

[57]  J. Mynard,et al.  One-Dimensional Haemodynamic Modeling and Wave Dynamics in the Entire Adult Circulation , 2015, Annals of Biomedical Engineering.

[58]  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..

[59]  Gianluigi Rozza,et al.  Simulation‐based uncertainty quantification of human arterial network hemodynamics , 2013, International journal for numerical methods in biomedical engineering.

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

[61]  Pablo J. Blanco,et al.  A two-level time step technique for the partitioned solution of one-dimensional arterial networks , 2012 .

[62]  Wolfgang Dauber,et al.  Pocket atlas of human anatomy , 1985 .

[63]  P J Blanco,et al.  On the integration of the baroreflex control mechanism in a heterogeneous model of the cardiovascular system , 2012, International journal for numerical methods in biomedical engineering.

[64]  S. Sherwin,et al.  Reduced modelling of blood flow in the cerebral circulation: Coupling 1‐D, 0‐D and cerebral auto‐regulation models , 2008 .

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