The numerical simulations for the blood flow in arteries by high order accurate schemes have a wide range of applications in medical engineering. The blood flow model admits the steady state solutions, in which the flux gradient is non-zero and is exactly balanced by the source term. In this paper, we present a high order finite volume weighted essentially non-oscillatory (WENO) scheme, which preserves the steady state solutions and maintains genuine high order accuracy for general solutions. The well-balanced property is obtained by a novel source term reformulation and discretisation, combined with well-balanced numerical fluxes. Extensive numerical experiments are carried out to verify well-balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.

We construct a nonlinear monotone finite volume scheme for three-dimensional diffusion equation on tetrahedral meshes. Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme, we present a new efficient eliminating method. The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously. The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes, and also show that our scheme appears to be approximate second-order accuracy for solution.

The extension of diamond scheme for diffusion equation to three dimensions is presented. The discrete normal flux is constructed by a linear combination of the directional flux along the line connecting cell-centers and the tangent flux along the cell-faces. In addition, it treats material discontinuities by a new iterative method. The stability and first-order convergence of the method is proved on distorted meshes. The numerical results illustrate that the method appears to be approximate second-order accuracy for solution.

In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws. High order accuracy in space is obtained with a standard WENO reconstruction algorithm and high order in time is obtained using the local space-time discontinuous Galerkin method recently proposed in [20]. In the Lagrangian framework considered here, the local space-time DG predictor is based on a weak formulation of the governing PDE on a moving space-time element. For the space-time basis and test functions we use Lagrange interpolation polynomials defined by tensor-product Gauss-Legendre quadrature points. The moving space-time elements are mapped to a reference element using an isoparametric approach, i.e. the space-time mapping is defined by the same basis functions as the weak solution of the PDE. We show some computational examples in one space-dimension for non-stiff and for stiff balance laws, in particular for the Euler equations of compressible gas dynamics, for the resistive relativistic MHD equations, and for the relativistic radiation hydrodynamics equations. Numerical convergence results are presented for the stiff case up to sixth order of accuracy in space and time and for the non-stiff case up to eighth order of accuracy in space and time.

We propose a numerical solution to incorporate in the simulation of a system of conservation laws boundary conditions that come from a microscopic modeling in the small mean free path regime. The typical example we discuss is the derivation of the Euler system from the BGK equation. The boundary condition relies on the analysis of boundary layers formation that accounts from the fact that the incoming kinetic flux might be far from the thermodynamic equilibrium.

We study a one-dimensional model for two-phase flows in heterogeneous media, in which the capillary pressure functions can be discontinuous with respect to space. We first give a model, leading to a system of degenerated nonlinear parabolic equations spatially coupled by nonlinear transmission conditions. We approximate the solution of our problem thanks to a monotonous finite volume scheme. The convergence of the underlying discrete solution to a weak solution when the discretization step tends to 0 is then proven. We also show, under assumptions on the initial data, a uniform estimate on the flux, which is then used during the uniqueness proof. A density argument allows us to relax the assumptions on the initial data and to extend the existence-uniqueness frame to a family of solution obtained as limit of approximations. A numerical example is then given to illustrate the behavior of the model.

This paper deals with the problem of numerical approximation in the Cauchy-Dirichlet problem for a scalar conservation law with a flux function having finitely many discontinuities. The well-posedness of this problem was proved by Carrillo [J. Evol. Eq. 3 (2003) 687–705]. Classical numerical methods do not allow us to compute a numerical solution (due to the lack of regularity of the flux). Therefore, we propose an implicit Finite Volume method based on an equivalent formulation of the initial problem. We show the well-posedness of the scheme and the convergence of the numerical solution to the entropy solution of the continuous problem. Numerical simulations are presented in the framework of Riemann problems related to discontinuous transport equation, discontinuous Burgers equation, discontinuous LWR equation and discontinuous non-autonomous Buckley-Leverett equation (lubrication theory).

We introduce a finite volume scheme for multi-dimensional drift-diffusion equations. Such equations arise from the theory of semiconductors and are composed of two continuity equations coupled with a Poisson equation. In the case that the continuity equations are non degenerate, we prove the convergence of the scheme and then the existence of solutions to the problem. The key point of the proof relies on the construction of an approximate gradient of the electric potential which allows us to deal with coupled terms in the continuity equations. Finally, a numerical example is given to show the efficiency of the scheme.

Our aim here is to study the thermal diffusion phenomenon in a forced convective flow. A system of nonlinear parabolic equations governs the evolution of the mass fractions in multicomponent mixtures. Some existence and uniqueness results are given under suitable conditions onstate functions. Then, we present a numerical scheme based on a "mixed finite element"method adapted to a finite volume scheme, of which we give numerical analysis. In a last part, we apply an homogenization technique to the studied equations in order to obtain an efficient modelling of Soret effect and adsorption in a porous medium at a macroscopic scale.