We present one- and two-dimensional central-upwind schemes for approximating solutions of the Saint-Venant system with source terms due to bottom topography. The Saint-Venant system has steady-state solutions in which nonzero flux gradients are exactly balanced by the source terms. It is a challenging problem to preserve this delicate balance with numerical schemes. Small perturbations of these states are also very difficult to compute. Our approach is based on extending semi-discrete central...

This paper concerns numerical methods for two-phase flows. The governing equations are the compressible 2-velocity, 2-pressure flow model. Pressure and velocity relaxation are included as source terms. Results obtained by a Godunov-type central scheme and a Roe-type upwind scheme are presented. Issues of preservation of pressure equilibrium, and positivity of the partial densities are addressed.

This paper concerns numerical methods for two-phase flows. The governing equations are the compressible 2-velocity, 2-pressure flow model. Pressure and velocity relaxation are included as source terms. Results obtained by a Godunov-type central scheme and a Roe-type upwind scheme are presented. Issues of preservation of pressure equilibrium, and positivity of the partial densities are addressed.

Numerical approximation schemes are discussed for the solution of contaminant transport with adsorption in dual-well flow. The method is based on time stepping and operator splitting for the transport with adsorption and diffusion. The nonlinear transport is solved by Godunov’s method. The nonlinear diffusion is solved by a finite volume method and by Newton’s type of linearization. The efficiency of the method is discussed.

We present and analyse in this paper a novel cell-centered collocated finite volume scheme for incompressible flows. Its definition involves a partition of the set of control volumes; each element of this partition is called a cluster and consists in a few neighbouring control volumes. Under a simple geometrical assumption for the clusters, we obtain that the pair of discrete spaces associating the classical cell-centered approximation for the velocities and cluster-wide constant pressures is inf-sup...

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...

Modeling the kinetics of a precipitation dissolution reaction occurring in a porous medium where diffusion also takes place leads to a system of two parabolic equations and one ordinary differential equation coupled with a stiff reaction term. This system is discretized by a finite volume scheme which is suitable for the approximation of the discontinuous reaction term of unknown sign. Discrete solutions are shown to exist and converge towards a weak solution of the continuous problem. Uniqueness...

This contribution deals with the numerical simulation of dislocation dynamics. Dislocations are described by means of the evolution of a family of closed or open smooth curves $\Gamma \left(t\right):S\to {ℝ}^2$, $t\geqq 0$. The curves are driven by the normal velocity $v$ which is the function of curvature $\kappa$ and the position. The evolution law reads as: $v=-\kappa +F$. The motion law is treated using direct approach numerically solved by two schemes, i. e., backward Euler semi-implicit and semi-discrete method of lines. Numerical stability is improved...

We consider the lowest-order Raviart–Thomas mixed finite element method for second-order elliptic problems on simplicial meshes in two and three space dimensions. This method produces saddle-point problems for scalar and flux unknowns. We show how to easily and locally eliminate the flux unknowns, which implies the equivalence between this method and a particular multi-point finite volume scheme, without any approximate numerical integration. The matrix of the final linear system is sparse, positive...

The subject of the paper is the derivation of error estimates for the combined finite volume-finite element method used for the numerical solution of nonstationary nonlinear convection-diffusion problems. Here we analyze the combination of barycentric finite volumes associated with sides of triangulation with the piecewise linear nonconforming Crouzeix-Raviart finite elements. Under some assumptions on the regularity of the exact solution, the $L^2\left(L^2\right)$ and $L^2\left(H^1\right)$ error estimates are established. At the end...

We propose a bi-dimensional finite volume extension of a continuous ALE method on unstructured cells whose edges are parameterized by rational quadratic Bezier curves. For each edge, the control point possess a weight that permits to represent any conic (see for example [LIGACH]) and thanks to [WAGUSEDE,WAGU], we are able to compute the exact area of our cells. We then give an extension of scheme for remapping step based on volume fluxing [MARSHA]...

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...

In this paper we present recent results for the bicharacteristic based finite volume schemes, the so-called finite volume evolution Galerkin (FVEG) schemes. These methods were proposed to solve multi-dimensional hyperbolic conservation laws. They combine the usually conflicting design objectives of using the conservation form and following the characteristics, or bicharacteristics. This is realized by combining the finite volume formulation with approximate evolution operators, which use bicharacteristics...

We present a fully adaptive multiresolution scheme for spatially one-dimensional quasilinear strongly degenerate parabolic equations with zero-flux and periodic boundary conditions. The numerical scheme is based on a finite volume discretization using the Engquist-Osher numerical flux and explicit time stepping. An adaptive multiresolution scheme based on cell averages is then used to speed up the CPU time and the memory requirements of the underlying finite volume scheme, whose first-order...

We study a family of non linear schemes for the numerical solution of linear advection on arbitrary grids in several space dimension. A proof of weak convergence of the family of schemes is given, based on a new Longitudinal Variation Diminishing (LVD) estimate. This estimate is a multidimensional equivalent to the well-known TVD estimate in one dimension. The proof uses a corollary of the Perron-Frobenius theorem applied to a generalized Harten formalism.

We study a family of non linear schemes for the numerical solution of linear advection on arbitrary grids in several space dimension. A proof of weak convergence of the family of schemes is given, based on a new Longitudinal Variation Diminishing (LVD) estimate. This estimate is a multidimensional equivalent to the well-known TVD estimate in one dimension. The proof uses a corollary of the Perron-Frobenius theorem applied to a generalized Harten formalism.

The aim of this article is a qualitative analysis of two modern finite volume (FVM) schemes. First one is the so called Modified Causon’s scheme, which is based on the classical MacCormack FVM scheme in total variation diminishing (TVD) form, but is simplified in such a way that the demands on computational power are much smaller without loss of accuracy. Second one is implicit WLSQR (Weighted Least Square Reconstruction) scheme combined with various types of numerical fluxes (AUSMPW+ and HLLC)....

In this work we describe an efficient model for the simulation of a two-phase flow made of a gas and a granular solid. The starting point is the two-velocity two-pressure model of Baer and Nunziato [Int. J. Multiph. Flow16 (1986) 861–889]. The model is supplemented by a relaxation source term in order to take into account the pressure equilibrium between the two phases and the granular stress in the solid phase. We show that the relaxation process can be made thermodynamically coherent with an...