A well‐balanced positivity‐preserving central‐upwind scheme for one‐dimensional blood flow models
- 29 June 2020
- journal article
- research article
- Published by Wiley in International Journal for Numerical Methods in Fluids
- Vol. 93 (2) , 369-395
Abstract
Summary: In this work, we consider a hyperbolic one‐dimensional (1D) model for blood flow through compliant axisymmetric tilted vessels. The pressure is a function of the cross‐sectional area and other model parameters. Important features of the model are inherited at the discrete level by the numerical scheme. For instance, the existence of steady states may provide important information about the flow properties at low computational cost. Here, we characterize a large class of smooth equilibrium solutions by means of quantities that remain invariant. At the discrete level, the well‐balanced property in the numerical scheme leads to accurate results when steady states are perturbed. On the other hand, the model is equipped with an entropy function and an entropy inequality that can help us choose the physically relevant weak solutions. A large class of semidiscrete entropy‐satisfying numerical schemes is described. In addition, preservation of positivity for the cross‐sectional area is achieved. Numerical results show the scheme is robust, stable, and accurate. The ultimate goal of this article is the numerical application to cases that are more relevant from the medical viewpoint. In particular, a numerical simulation of cardiac cycles with appropriate parameters shows that increasing the rigidity of the artery walls delays the formation of shock waves. Gravity effects are also measured in tilted vessels, and a simulation using an idealized aorta model was conducted.Keywords
Funding Information
- Secretaría de Ciencia, Humanidades, Tecnología e Innovación (A1‐S‐17634)
This publication has 49 references indexed in Scilit:
- Patient-Specific Modeling of Cardiovascular MechanicsAnnual Review of Biomedical Engineering, 2009
- A central scheme for shallow water flows along channels with irregular geometryESAIM: Mathematical Modelling and Numerical Analysis, 2008
- An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatmentESAIM: Mathematical Modelling and Numerical Analysis, 2008
- A 1D arterial blood flow model incorporating ventricular pressure, aortic valve and regional coronary flow using the locally conservative Galerkin (LCG) methodCommunications in Numerical Methods in Engineering, 2008
- Blood Flow in Compliant Arteries: An Effective Viscoelastic Reduced Model, Numerics, and Experimental ValidationAnnals of Biomedical Engineering, 2006
- Mathematical analysis of the quasilinear effects in a hyperbolic model blood flow through compliant axi‐symmetric vesselsMathematical Methods in the Applied Sciences, 2003
- Solution of two‐dimensional Riemann problems for gas dynamics without Riemann problem solversNumerical Methods for Partial Differential Equations, 2002
- Central-Upwind Schemes for the Saint-Venant SystemESAIM: Mathematical Modelling and Numerical Analysis, 2002
- Blood flow through compliant vessels after endovascular repair: wall deformations induced by the discontinuous wall propertiesComputing and Visualization in Science, 2002
- A steady-state capturing method for hyperbolic systems with geometrical source termsESAIM: Mathematical Modelling and Numerical Analysis, 2001