Do you want to publish a course? Click here

A shallow water with variable pressure model for blood flow simulation

151   0   0.0 ( 0 )
 Added by Olivier Delestre
 Publication date 2015
  fields
and research's language is English




Ask ChatGPT about the research

We performed numerical simulations of blood flow in arteries with a variable stiffness and cross-section at rest using a finite volume method coupled with a hydrostatic reconstruction of the variables at the interface of each mesh cell. The method was then validated on examples taken from the literature. Asymptotic solutions were computed to highlight the effect of the viscous and viscoelastic source terms. Finally, the blood flow was computed in an artery where the cross-section at rest and the stiffness were varying. In each test case, the hydrostatic reconstruction showed good results where other simpler schemes did not, generating spurious oscillations andnonphysical velocities.



rate research

Read More

174 - Olivier Delestre 2012
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.
A new two-dimensional model for blood flows in arteries with arbitrary cross sections is derived. The model consists of a system of balance laws for conservation of mass and balance of momentum in the axial and angular directions. The equations are derived by applying asymptotic analysis to the incompressible Navier-Stokes equations in narrow, large vessels and integrating in the radial direction in each cross section. The main properties of the system are discussed and a positivity-preserving well-balanced central-upwind scheme is presented. The merits of the scheme will be tested in a variety of scenarios. In particular, numerical results of simulations using an idealized aorta model are shown. We analyze the time evolution of the blood flow under different initial conditions such as perturbations to steady states consisting of a bulging in the vessels wall. We consider different situations given by distinct variations in the vessels elasticity.
In this paper, we propose a direct Eulerian generalized Riemann problem (GRP) scheme for a blood flow model in arteries. It is an extension of the Eulerian GRP scheme, which is developed by Ben-Artzi, et. al. in J. Comput. Phys., 218(2006). By using the Riemann invariants, we diagonalize the blood flow system into a weakly coupled system, which is used to resolve rarefaction wave. We also use Rankine-Hugoniot condition to resolve the local GRP formulation. We pay special attention to the acoustic case as well as the sonic case. The extension to the two dimensional case is carefully obtained by using the dimensional splitting technique. We test that the derived GRP scheme is second order accuracy.
The blood flow model maintains the steady state solutions, in which the flux gradients are non-zero but exactly balanced by the source term. In this paper, we design high order finite difference weighted non-oscillatory (WENO) schemes to this model with such well-balanced property and at the same time keeping genuine high order accuracy. Rigorous theoretical analysis as well as extensive numerical results all indicate that the resulting schemes verify high order accuracy, maintain the well-balanced property, and keep good resolution for smooth and discontinuous solutions.
192 - Emmanuel Audusse 2013
In this work we are interested in numerical simulations for bedload erosion processes. We present a relaxation solver that we apply to moving dunes test cases in one and two dimensions. In particular we retrieve the so-called anti-dune process that is well described in the experiments. In order to be able to run 2D test cases with reasonable CPU time, we also describe and apply a parallelization procedure by using domain decomposition based on the classical MPI library.
comments
Fetching comments Fetching comments
Sign in to be able to follow your search criteria
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا