In the present work we investigate the numerical simulation of liquid-vapor phase change in compressible flows. Each phase is modeled as a compressible fluid equipped with its own equation of state (EOS). We suppose that inter-phase equilibrium processes in the medium operate at a short time-scale compared to the other physical phenomena such as convection or thermal diffusion. This assumption provides an implicit definition of an equilibrium EOS...

Inspired by the growing use of non linear discretization techniques for the linear diffusion equation in industrial codes, we construct and analyze various explicit non linear finite volume schemes for the heat equation in dimension one. These schemes are inspired by the Le Potier’s trick [C. R. Acad. Sci. Paris, Ser. I 348 (2010) 691–695]. They preserve the maximum principle and admit a finite volume formulation. We provide a original functional setting for the analysis of convergence of such methods....

This paper presents and summarize our results concerning the nonlinear tensor diffusion which enhances image structure coherence. The core of the paper comes from [3, 2, 4, 5]. First we briefly describe the diffusion model and provide its basic properties. Further we build a semi-implicit finite volume scheme for the above mentioned model with the help of a co-volume mesh. This strategy is well-known as diamond-cell method owing to the choice of co-volume as a diamondshaped polygon, see [1]. We...

Stability and convergence of the linear semi-implicit discrete duality finite volume (DDFV) numerical scheme in 2D for the solution of the regularized curvature driven level set equation is proved. Numerical experiments concerning comparison with exact solution and image filtering problem using proposed scheme are included.

We discuss the numerical modeling of heat exchange between the infiltrated water and porous media matrix. An unsaturated-saturated flow is considered with boundary conditions reflecting the external driven forces. The developed numerical method is efficient and can be used for solving the inverse problems concerning determination of transmission coefficients for heat energy exchange inside and also on the boundary of porous media. Numerical experiments support our method.

We propose a new numerical scheme based on the finite volumes to simulate the urethra flow based on hyperbolic balance law. Our approach is based on the Riemann solver designed for the augmented quasilinear homogeneous formulation. The scheme has general semidiscrete wave–propagation form and can be extended to arbitrary high order accuracy. The first goal is to construct the scheme, which is well balanced, i.e. maintains not only some special steady states but all steady states which can occur....

Analyzing the validity and success of a data assimilation algorithmwhen some state variable observations are not available is an important problem in meteorology and engineering. We present an improved data assimilation algorithm for recovering the exact full reference solution (i.e. the velocity and temperature) of the 3D Planetary Geostrophic model, at an exponential rate in time, by employing coarse spatial mesh observations of the temperature alone. This provides, in the case of this paradigm,...

We study a two-phase pipe flow model with relaxation terms in the momentum and energy equations, driving the model towards dynamic and thermal equilibrium. These equilibrium states are characterized by the velocities and temperatures being equal in each phase. For each of these relaxation processes, we consider the limits of zero and infinite relaxation times. By expanding on previously established results, we derive a formulation of the mixture sound velocity for the thermally relaxed model. This...

We study a two-phase pipe flow model with relaxation terms in the momentum and energy equations, driving the model towards dynamic and thermal equilibrium. These equilibrium states are characterized by the velocities and temperatures being equal in each phase. For each of these relaxation processes, we consider the limits of zero and infinite relaxation times. By expanding on previously established results, we derive a formulation of the mixture sound velocity for the thermally relaxed model. This...

In this paper, we study a postprocessing procedure for improving accuracy of the finite volume element approximations of semilinear parabolic problems. The procedure amounts to solve a source problem on a coarser grid and then solve a linear elliptic problem on a finer grid after the time evolution is finished. We derive error estimates in the L2 and H1 norms for the standard finite volume element scheme and an improved error estimate in the H1 norm. Numerical results demonstrate the accuracy...

A general class of nonconforming meshes has been recently studied for stationary anisotropic heterogeneous diffusion problems, see Eymard et al. (IMA J. Numer. Anal. 30 (2010), 1009–1043). Thanks to the basic ideas developed in the stated reference for stationary problems, we derive a new discretization scheme in order to approximate the nonstationary heat problem. The unknowns of this scheme are the values at the centre of the control volumes, at some internal interfaces, and at the mesh points...

There are many methods and approaches to solving convection--diffusion problems. For those who want to solve such problems the situation is very confusing and it is very difficult to choose the right method. The aim of this short overview is to provide basic guidelines and to mention the common features of different methods. We place particular emphasis on the concept of linear and non-linear stabilization and its implementation within different approaches.

We show stability and consistency of the linear semi-implicit complementary volume numerical scheme for solving the regularized, in the sense of Evans and Spruck, mean curvature flow equation in the level set formulation. The numerical method is based on the finite volume methodology using the so-called complementary volumes to a finite element triangulation. The scheme gives the solution in an efficient and unconditionally stable way.

We study the use of a GPU for the numerical approximation of the curvature dependent flows of graphs - the mean-curvature flow and the Willmore flow. Both problems are often applied in image processing where fast solvers are required. We approximate these problems using the complementary finite volume method combined with the method of lines. We obtain a system of ordinary differential equations which we solve by the Runge-Kutta-Merson solver. It is a robust solver with an automatic choice of the...

We demonstrate some a priori estimates of a scheme using stabilization and hybrid interfaces applying to partial differential equations describing miscible displacement in porous media. This system is made of two coupled equations: an anisotropic diffusion equation on the pressure and a convection-diffusion-dispersion equation on the concentration of invading fluid. The anisotropic diffusion operators in both equations require special care while discretizing by a finite volume method SUSHI. Later,...

When applied to the linear advection problem in dimension two, the upwind finite volume method is a non consistent scheme in the finite differences sense but a convergent scheme. According to our previous paper [Bouche et al., SIAM J. Numer. Anal.43 (2005) 578–603], a sufficient condition in order to complete the mathematical analysis of the finite volume scheme consists in obtaining an estimation of order p, less or equal to one, of a quantity that depends only on the mesh and on the advection ...

This work aims at evaluating in practical situations the capability of the mesh refinement technique based on the multiresolution adaptive method coupled with high resolution spatial and temporal approximations, to recover elementary physical mechanisms by achieving gains in both CPU time and memory use compared to single grid computations. We first present a summary of the multiresolution procedure. We then describe MR algorithms. Finally, the evaluation of the method is presented on several well...

This paper concerns the discretization of multiphase Darcy flows, in the case of heterogeneous anisotropic porous media and general 3D meshes used in practice to represent reservoir and basin geometries. An unconditionally coercive and symmetric vertex centred approach is introduced in this paper. This scheme extends the Vertex Approximate Gradient scheme (VAG), already introduced for single phase diffusive problems in [9], to multiphase Darcy flows....