A 2D nonlinear multiring model for blood flow in large elastic arteries

In this paper, we propose a two-dimensional nonlinear "multiring" model for blood flow in axisymmetric elastic arteries. It is designed to overcome the numerical difficulties of three-dimensional fluid-structure interaction simulations of blood flow without using the oversimplifications necessary to obtain one-dimensional models of blood flow. This multiring model is derived by integrating over concentric rings of fluid the simplified long-wave Navier-Stokes equations coupled to an elastic model of the arterial wall. The resulting system of balance laws provides a unified framework in which both the motion of the fluid and the displacement of the wall are dealt with simultaneously. The mathematical structure of the multiring model allows us to use a finite volume method that guarantees the conservation of mass and the positivity of the numerical solution and can deal with nonlinear flows and large deformations of the arterial wall. We show that the finite volume numerical solution of the multiring model provides at a reasonable computational cost an asymptotically valid description of blood flow velocity profiles and other averaged quantities (wall shear stress, flow rate, ...) in large elastic and quasi-rigid arteries. In particular, we validate the multiring model against well-known solutions such as the Womersley or the Poiseuille solutions as well as against steady boundary layer solutions in quasi-rigid constricted and expanded tubes.

[1]  D. Korteweg Ueber die Fortpflanzungsgeschwindigkeit des Schalles in elastischen Rhren , 1878 .

[2]  Suncica Canic,et al.  Self-Consistent Effective Equations Modeling Blood Flow in Medium-to-Large Compliant Arteries , 2005, Multiscale Model. Simul..

[3]  F. Smith FLOW THROUGH CONSTRICTED OR DILATED PIPES AND CHANNELS: PART 1 , 1976 .

[4]  F. Bouchut Construction of BGK Models with a Family of Kinetic Entropies for a Given System of Conservation Laws , 1999 .

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

[6]  Nan Xiao,et al.  Multi-scale computational model of three-dimensional hemodynamics within a deformable full-body arterial network , 2013, J. Comput. Phys..

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

[8]  Ender A Finol,et al.  Flow-induced Wall Shear Stress in Abdominal Aortic Aneurysms: Part I - Steady Flow Hemodynamics , 2002, Computer methods in biomechanics and biomedical engineering.

[9]  S. Ling,et al.  A nonlinear analysis of pulsatile flow in arteries , 1972, Journal of Fluid Mechanics.

[10]  Emmanuel Audusse,et al.  A well-balanced positivity preserving second-order scheme for shallow water flows on unstructured meshes , 2005 .

[11]  J. Womersley XXIV. Oscillatory motion of a viscous liquid in a thin-walled elastic tube—I: The linear approximation for long waves , 1955 .

[12]  Raimund Bürger,et al.  A multilayer shallow water system for polydisperse sedimentation , 2013, J. Comput. Phys..

[13]  S. Wille Pulsatile pressure and flow in an arterial aneurysm simulated in a mathematical model. , 1981, Journal of biomedical engineering.

[14]  M. Toffolon,et al.  Steady analysis of transcritical flows in collapsible tubes with discontinuous mechanical properties: implications for arteries and veins , 2013, Journal of Fluid Mechanics.

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

[16]  Charles A. Taylor,et al.  Patient-Specific Modeling of Blood Flow and Pressure in Human Coronary Arteries , 2010, Annals of Biomedical Engineering.

[17]  Olivier Delestre,et al.  A shallow water with variable pressure model for blood flow simulation , 2015, Networks Heterog. Media.

[18]  Pierre-Yves Lagrée,et al.  The dicrotic notch analyzed by a numerical model , 2016, Comput. Biol. Medicine.

[19]  G. Barrenechea,et al.  A finite element method for the resolution of the Reduced Navier‐Stokes/Prandtl equations , 2009 .

[20]  F. N. Vosse,et al.  A wave propagation model of blood flow in large vessels using an approximate velocity profile function , 2007, Journal of Fluid Mechanics.

[21]  J. D. Humphrey,et al.  A haemodynamic predictor of intraluminal thrombus formation in abdominal aortic aneurysms , 2014, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  H. Fernholz Boundary Layer Theory , 2001 .

[23]  M. Dumbser,et al.  An efficient semi-implicit finite volume method for axially symmetric compressible flows in compliant tubes , 2015 .

[24]  Lei Wang,et al.  A lattice Boltzmann study of the asymmetry effect on the hemodynamics in stented fusiform aneurysms , 2016, Comput. Math. Appl..

[25]  Subin Solomen,et al.  Biomechanics of Circulation , 2018 .

[26]  S. Sherwin,et al.  Analysing the pattern of pulse waves in arterial networks: a time-domain study , 2009 .

[27]  Fayssal Benkhaldoun,et al.  A fast finite volume solver for multi-layered shallow water flows with mass exchange , 2014, J. Comput. Phys..

[28]  P. Lagrée,et al.  Comparison of computations of asymptotic flow models in a constricted channel , 2012 .

[29]  Pierre-Yves Lagrée,et al.  One-dimensional model for propagation of a pressure wave in a model of the human arterial network: comparison of theoretical and experimental results. , 2011, Journal of biomechanical engineering.

[30]  D F Young,et al.  An integral method for the analysis of flow in arterial stenoses. , 1974, Bulletin of mathematical biology.

[31]  Sylvie Lorthois,et al.  The RNS/Prandtl equations and their link with other asymptotic descriptions: Application to the wall shear stress scaling in a constricted pipe , 2005 .

[32]  T. Tezduyar Computation of moving boundaries and interfaces and stabilization parameters , 2003 .

[33]  N R Cholvin,et al.  Pressure Drop across Artificially Induced Stenoses in the Femoral Arteries of Dogs , 1975, Circulation research.

[34]  L. Walløe,et al.  Pulsatile pressure and flow in arterial stenoses simulated in a mathematical model. , 1981, Journal of biomedical engineering.

[35]  K. Pohlhausen,et al.  Zur näherungsweisen Integration der Differentialgleichung der Iaminaren Grenzschicht , 1921 .

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

[37]  J. Lambert,et al.  On the nonlinearities of fluid flow in nonrigid tubes , 1958 .

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

[39]  A. Shapiro Steady Flow in Collapsible Tubes , 1977 .

[40]  Pierre-Yves Lagrée,et al.  An inverse technique to deduce the elasticity of a large artery , 2000 .

[41]  R. Budwig,et al.  Steady flow in abdominal aortic aneurysm models. , 1993, Journal of biomechanical engineering.

[42]  Giovanna Guidoboni,et al.  Blood Flow in Compliant Arteries: An Effective Viscoelastic Reduced Model, Numerics, and Experimental Validation , 2006, Annals of Biomedical Engineering.

[43]  Laurent Gosse,et al.  Un schéma-équilibre adapté aux lois de conservation scalaires non-homogènes , 1996 .

[44]  J. Alastruey,et al.  A Novel Analytical Approach to Pulsatile Blood Flow in the Arterial Network , 2016, Annals of Biomedical Engineering.

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

[46]  Alfio Quarteroni,et al.  Geometric multiscale modeling of the cardiovascular system, between theory and practice , 2016 .

[47]  K. Perktold On the paths of fluid particles in an axisymmetrical aneurysm. , 1987, Journal of biomechanics.

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

[49]  L Prandtl,et al.  Motion of fluids with very little viscosity , 1928 .

[50]  Charbel Farhat,et al.  The discrete geometric conservation law and the nonlinear stability of ALE schemes for the solution of flow problems on moving grids , 2001 .

[51]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[52]  R. Ogden,et al.  A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models , 2000 .

[53]  Jacques Sainte-Marie,et al.  A multilayer Saint-Venant system with mass exchanges for shallow water flows. Derivation and numerical validation* , 2009, 0901.3887.

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

[55]  Nasir Ali,et al.  Numerical simulation of unsteady micropolar hemodynamics in a tapered catheterized artery with a combination of stenosis and aneurysm , 2015, Medical & Biological Engineering & Computing.

[56]  Tayfun E. Tezduyar,et al.  Modelling of fluid–structure interactions with the space–time finite elements: Arterial fluid mechanics , 2007 .

[57]  Thomas J. R. Hughes,et al.  Finite element modeling of blood flow in arteries , 1998 .

[58]  Pierre-Yves Lagrée,et al.  Verification and comparison of four numerical schemes for a 1D viscoelastic blood flow model , 2013, Computer methods in biomechanics and biomedical engineering.

[59]  A. Marsden,et al.  A Primer on Computational Simulation in Congenital Heart Disease for the Clinician , 2010, 1101.3726.

[60]  J. M. Siegel,et al.  A scaling law for wall shear rate through an arterial stenosis. , 1994, Journal of biomechanical engineering.

[61]  R. Courant,et al.  On the Partial Difference Equations, of Mathematical Physics , 2015 .

[62]  Matthias Mayr,et al.  A Temporal Consistent Monolithic Approach to Fluid-Structure Interaction Enabling Single Field Predictors , 2015, SIAM J. Sci. Comput..

[63]  Alison L. Marsden,et al.  Patient-Specific Multiscale Modeling of Blood Flow for Coronary Artery Bypass Graft Surgery , 2012, Annals of Biomedical Engineering.

[64]  Christopher J. Arthurs,et al.  A mathematical model of coronary blood flow control: simulation of patient-specific three-dimensional hemodynamics during exercise , 2016, American journal of physiology. Heart and circulatory physiology.

[65]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[66]  Alison L Marsden,et al.  Patient-Specific Simulations Reveal Significant Differences in Mechanical Stimuli in Venous and Arterial Coronary Grafts , 2016, Journal of Cardiovascular Translational Research.

[67]  Marie-Odile Bristeau,et al.  Boundary Conditions for the Shallow Water Equations solved by Kinetic Schemes , 2001 .

[68]  Charles A. Taylor,et al.  A coupled momentum method for modeling blood flow in three-dimensional deformable arteries , 2006 .

[69]  Jean-Frédéric Gerbeau,et al.  Kinetic scheme for arterial and venous blood flow, and application to partial hepatectomy modeling. , 2017 .

[70]  Lawrence C. Cerny,et al.  Leonhardi Euleri's “Principia pro motu sanguinis per arterias determinando” , 1974 .

[71]  Olivier Delestre,et al.  Low-Shapiro hydrostatic reconstruction technique for blood flow simulation in large arteries with varying geometrical and mechanical properties , 2017, J. Comput. Phys..