The application of the method of differential relations to the equations of continuum mechanics
The correct use of the Lax–Friedrichs method
We are concerned with the structure of the operator corresponding to the Lax–Friedrichs method. At first, the phenomenae which may arise by the naive use of the Lax–Friedrichs scheme are analyzed. In particular, it turns out that the correct definition of the method has to include the details of the discretization of the initial condition and the computational domain. Based on the results of the discussion, we give a recipe that ensures that the number of extrema within the discretized version of...
The correct use of the Lax–Friedrichs method
We are concerned with the structure of the operator corresponding to the Lax–Friedrichs method. At first, the phenomenae which may arise by the naive use of the Lax–Friedrichs scheme are analyzed. In particular, it turns out that the correct definition of the method has to include the details of the discretization of the initial condition and the computational domain. Based on the results of the discussion, we give a recipe that ensures that the number of extrema within the discretized version...
The finite difference scheme for 3D mathematical modeling of a wood drying process.
The Method of Fractional Steps for Conservation Laws.
The Neumann stability of high-order symmetric schemes for convection-diffusion problems.
The numerical convergence of the Landau-Lifshitz equations and its simulation.
The numerical interface coupling of nonlinear hyperbolic systems of conservation laws : II. The case of systems
We study the theoretical and numerical coupling of two hyperbolic systems of conservation laws at a fixed interface. As already proven in the scalar case, the coupling preserves in a weak sense the continuity of the solution at the interface without imposing the overall conservativity of the coupled model. We develop a detailed analysis of the coupling in the linear case. In the nonlinear case, we either use a linearized approach or a coupling method based on the solution of a Riemann problem. We...
The numerical interface coupling of nonlinear hyperbolic systems of conservation laws: II. The case of systems
We study the theoretical and numerical coupling of two hyperbolic systems of conservation laws at a fixed interface. As already proven in the scalar case, the coupling preserves in a weak sense the continuity of the solution at the interface without imposing the overall conservativity of the coupled model. We develop a detailed analysis of the coupling in the linear case. In the nonlinear case, we either use a linearized approach or a coupling method based on the solution of a Riemann problem....
The Rothe method for solving the first initial-boundary value problem for the equation of filtration in a cracked porous medium.
The splitting in potential Crank–Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip
We consider an initial-boundary value problem for a generalized 2D time-dependent Schrödinger equation (with variable coefficients) on a semi-infinite strip. For the Crank–Nicolson-type finite-difference scheme with approximate or discrete transparent boundary conditions (TBCs), the Strang-type splitting with respect to the potential is applied. For the resulting method, the unconditional uniform in time L2-stability is proved. Due to the splitting, an effective direct algorithm using FFT is developed...
Theoretical analysis of the upwind finite volume scheme on the counter-example of Peterson
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 ...
Theoretical and numerical study of a quasi-linear Zakharov system describing Landau damping
In this paper, we study a Zakharov system coupled to an electron diffusion equation in order to describe laser-plasma interactions. Starting from the Vlasov-Maxwell system, we derive a nonlinear Schrödinger like system which takes into account the energy exchanged between the plasma waves and the electrons via Landau damping. Two existence theorems are established in a subsonic regime. Using a time-splitting, spectral discretizations for the Zakharov system and a finite difference scheme for the...
Theoretical and numerical study of a quasi-linear Zakharov system describing Landau damping
In this paper, we study a Zakharov system coupled to an electron diffusion equation in order to describe laser-plasma interactions. Starting from the Vlasov-Maxwell system, we derive a nonlinear Schrödinger like system which takes into account the energy exchanged between the plasma waves and the electrons via Landau damping. Two existence theorems are established in a subsonic regime. Using a time-splitting, spectral discretizations for the Zakharov system and a finite difference scheme for...
Theory and numerical approximations for a nonlinear 1 + 1 Dirac system
We consider a nonlinear Dirac system in one space dimension with periodic boundary conditions. First, we discuss questions on the existence and uniqueness of the solution. Then, we propose an implicit-explicit finite difference method for its approximation, proving optimal order a priori error estimates in various discrete norms and showing results from numerical experiments.
Theory and numerical approximations for a nonlinear 1 + 1 Dirac system
We consider a nonlinear Dirac system in one space dimension with periodic boundary conditions. First, we discuss questions on the existence and uniqueness of the solution. Then, we propose an implicit-explicit finite difference method for its approximation, proving optimal order a priori error estimates in various discrete norms and showing results from numerical experiments.
Three-level difference schemes on non-uniform in time grids.
Time splitting for wave equations in random media
Numerical simulation of high frequency waves in highly heterogeneous media is a challenging problem. Resolving the fine structure of the wave field typically requires extremely small time steps and spatial meshes. We show that capturing macroscopic quantities of the wave field, such as the wave energy density, is achievable with much coarser discretizations. We obtain such a result using a time splitting algorithm that solves separately and successively propagation and scattering in the simplified...
Time splitting for wave equations in random media
Numerical simulation of high frequency waves in highly heterogeneous media is a challenging problem. Resolving the fine structure of the wave field typically requires extremely small time steps and spatial meshes. We show that capturing macroscopic quantities of the wave field, such as the wave energy density, is achievable with much coarser discretizations. We obtain such a result using a time splitting algorithm that solves separately and successively propagation and scattering in the...
Transport of pollutant in shallow water : a two time steps kinetic method
The aim of this paper is to present a finite volume kinetic method to compute the transport of a passive pollutant by a flow modeled by the shallow water equations using a new time discretization that allows large time steps for the pollutant computation. For the hydrodynamic part the kinetic solver ensures – even in the case of a non flat bottom – the preservation of the steady state of a lake at rest, the non-negativity of the water height and the existence of an entropy inequality. On an other...