Coupling between lumped and distributed models for blood flow problems

Abstract.In this paper we propose a method for coupling distributed and lumped models for the blood circulation. Lumped parameter models, based on an analogy between the circulatory system and an electric or a hydraulic network are widely employed in the literature to investigate different systemic responses in physiologic and pathologic situations (see e.g. [13, 24, 30, 15, 4, 27, 11, 14]). From the mathematical viewpoint these models are represented by ordinary differential equations. On the other hand, for the accurate description of local phenomena, the Navier–Stokes equations for incompressible fluids are considered. In the multiscale perspective, lumped models have been adopted (see e.g. [16]) as a numerical preprocessor to provide a quantitative estimate of the boundary conditions at the interfaces. However, the two solvers (i.e. the lumped and the distributed one) have been used separately. In the present work, we introduce a genuinely heterogeneous multiscale approach where the local model and the systemic one are coupled at a mathematical and numerical level and solved together. In this perspective, we have no longer boundary conditions on the artificial sections, but interface conditions matching the two submodels. The mathematical model and its numerical approximation are carefully addressed and several test cases are considered.

[1]  J. Womersley.  Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[2]  D. A. Mcdonald.  Blood flow in arteries , 1974 .

[3]  V. Rideout,et al.  Difference-differential equations for fluid flow in distensible tubes. , 1967, IEEE transactions on bio-medical engineering.

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

[5]  Archie Jp,et al.  Critical stenosis of the internal carotid artery. , 1981, Surgery.

[6]  Abraham Noordergraaf,et al.  Reduced Models of Arterial Systems , 1985, IEEE Transactions on Biomedical Engineering.

[7]  G. Hutchins,et al.  Shear-dependent thickening of the human arterial intima. , 1986, Atherosclerosis.

[8]  Michael M. Resch,et al.  Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model. , 1991, Journal of biomechanical engineering.

[9]  K B Campbell,et al.  Fit to diastolic arterial pressure by third-order lumped model yields unreliable estimates of arterial compliance. , 1996, Medical engineering & physics.

[10]  Alfio Quarteroni,et al.  Domain Decomposition Methods for Partial Differential Equations , 1999 .

[11]  Rolf Rannacher,et al.  ARTIFICIAL BOUNDARIES AND FLUX AND PRESSURE CONDITIONS FOR THE INCOMPRESSIBLE NAVIER–STOKES EQUATIONS , 1996 .

[12]  R Pietrabissa,et al.  A mathematical model of circulation in the presence of the bidirectional cavopulmonary anastomosis in children with a univentricular heart. , 1997, Medical engineering & physics.

[13]  Mette S Olufsen,et al.  Structured tree outflow condition for blood flow in larger systemic arteries. , 1999, American journal of physiology. Heart and circulatory physiology.

[14]  Alfio Quarteroni,et al.  Multiscale modelling of the circulatory system: a preliminary analysis , 1999 .

[15]  Alfio Quarteroni,et al.  From the global cardiovascular hemodynamics down to the local blood motion: preliminary applications of a multiscale approach , 2000 .

[16]  Alfio Quarteroni,et al.  Computational vascular fluid dynamics: problems, models and methods , 2000 .

[17]  A. Quarteroni,et al.  Factorization methods for the numerical approximation of Navier-Stokes equations , 2000 .

[18]  Stefania Ragni,et al.  Multiscale modelling of the circulatory system , 2001 .

[19]  R Pietrabissa,et al.  Multiscale modelling as a tool to prescribe realistic boundary conditions for the study of surgical procedures. , 2002, Biorheology.