A "well-balanced" finite volume scheme for blood flow simulation

Abstract

We are interested in simulating blood flow in arteries with a one dimensional model. Thanks to recent developments in the analysis of hyperbolic system of conservation laws (in the Saint-Venant/ shallow water equations context) we will perform a simple finite volume scheme. We focus on conservation properties of this scheme which were not previously considered. To emphasize the necessity of this scheme, we present how a too simple numerical scheme may induce spurious flows when the basic static shape of the radius changes. On contrary, the proposed scheme is "well-balanced": it preserves equilibria of Q = 0. Then examples of analytical or linearized solutions with and without viscous damping are presented to validate the calculations. The influence of abrupt change of basic radius is emphasized in the case of an aneurism.

36 pages

References in corpus (1)

• Why many theories of shock waves are necessary. Convergence error in formally path-consistent schemes

Cited by in corpus (7)

• SWASHES: a compilation of Shallow Water Analytic Solutions for Hydraulic and Environmental Studies
• Numerical method of characteristics for one-dimensional blood flow
• Low-Shapiro hydrostatic reconstruction technique for blood flow simulation in large arteries with varying geometrical and mechanical properties
• Well-balanced finite difference WENO schemes for the blood flow model
• A well-balanced finite volume scheme for 1D hemodynamic simulations
• A fully well-balanced hydrodynamic reconstruction
• Well-balanced high-order method for non-conservative hyperbolic PDEs with source terms: application to one-dimensional blood flow equations with gravity