Overlapping Schwarz waveform relaxation method for the solution of the forward–backward heat equation

In this article we analyzed the convergence of the Schwarz waveform relaxation method for solving the forward–backward heat equation. Numerical results are presented for a specific type of model problem.     

A priori and a posteriori error analyses in the study of viscoelastic problems

In this work, the numerical approximation of a viscoelastic problem is studied. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. Then, two numerical analyses are presented. First, a priori estimates are proved from which the linear convergence of the algorithm is derived under suitable regularity conditions. Secondly, an a posteriori error analysis is provided extending some preliminary results obtained in the study of the heat equation. Upper and lower error bounds are obtained.     

Numerical analysis for a conservative difference scheme to solve the Schrödinger-Boussinesq equation

In this paper, we present a finite difference scheme for the solution of an initial-boundary value problem of the Schrödinger-Boussinesq equation. The scheme is fully implicit and conserves two invariable quantities of the system. We investigate the existence of the solution for the scheme, give computational process for the numerical solution and prove convergence of iteration method by which a nonlinear algebra system for unknown Vⁿ⁺¹ is solved. On the basis of a priori estimates for a numerical solution, the uniqueness, convergence and stability for the difference solution is discussed. Numerical experiments verify the accuracy of our method.     

A deterministic particle method for the Vlasov–Fokker–Planck equation in one dimension

The Vlasov–Fokker–Planck equation is a model for a collisional, electrostatic plasma. The approximation of this equation in one spatial dimension is studied. The equation under consideration is linear in that the electric field is given as a known function that is not internally consistent with the phase space distribution function. The approximation method applied is the deterministic particle method described in Wollman and Ozizmir [Numerical approximation of the Vlasov–Poisson–Fokker–Planck system in one dimension, J. Comput. Phys. 202 (2005) 602–644]. For the present linear problem an analysis of the stability and convergence of the numerical method is carried out. In addition, computations are done that verify the convergence of the numerical solution. It is also shown that the long term asymptotics of the computed solution is in agreement with the steady state solution derived in Bouchut and Dolbeault [On long time asymptotics of the Vlasov–Fokker–Planck equation and of the Vlasov–Poisson–Fokker–Planck system with coulombic and Newtonian potentials, Differential Integral Equations 8(3) (1995) 487–514].     

A singular perturbation of the heat equation with memory

In this paper we consider a hyperbolic equation, with a memory term in time, which can be seen as a singular perturbation of the heat equation with memory. The qualitative properties of the solutions of the initial boundary value problems associated with both equations are studied. We propose numerical methods for the hyperbolic and parabolic models and their stability properties are analyzed. Finally, we include numerical experiments illustrating the performance of those methods.     

The solution of the two-dimensional sine-Gordon equation using the method of lines

The method of lines is used to transform the initial/boundary-value problem associated with the two-dimensional sine-Gordon equation in two space variables into a second-order initial-value problem. The finite-difference methods are developed by replacing the matrix-exponential term in a recurrence relation with rational approximants. The resulting finite-difference methods are analyzed for local truncation error, stability and convergence. To avoid solving the nonlinear system a predictor–corrector scheme using the explicit method as predictor and the implicit as corrector is applied. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given.     

Convergence analysis of the FEM approximation of the first order projection method for incompressible flows with and without the inf-sup condition

In this paper we obtain convergence results for the fully discrete projection method for the numerical approximation of the incompressible Navier–Stokes equations using a finite element approximation for the space discretization. We consider two situations. In the first one, the analysis relies on the satisfaction of the inf-sup condition for the velocity-pressure finite element spaces. After that, we study a fully discrete fractional step method using a Poisson equation for the pressure. In this case the velocity-pressure interpolations do not need to accomplish the inf-sup condition and in fact we consider the case in which equal velocity-pressure interpolation is used. Optimal convergence results in time and space have been obtained in both cases.     

Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy

Summary A fully discrete finite element method for the Cahn-Hilliard equation with a logarithmic free energy based on the backward Euler method is analysed. Existence and uniqueness of the numerical solution and its convergence to the solution of the continuous problem are proved. Two iterative schemes to solve the resulting algebraic problem are proposed and some numerical results in one space dimension are presented.     

Numerical approximation of the Cahn-Larché equation

Summary Spinodal decomposition, i.e., the separation of a homogeneous mixture into different phases, can be modeled by the Cahn-Hilliard equation - a fourth order semilinear parabolic equation. If elastic stresses due to a lattice misfit become important, the Cahn-Hilliard equation has to be coupled to an elasticity system to take this into account. Here, we present a discretization based on finite elements and an implicit Euler scheme. We first show solvability and uniqueness of solutions. Based on an energy decay property we then prove convergence of the scheme. Finally we present numerical experiments showing the impact of elasticity on the morphology of the microstructure.Research supported by DFG Priority Program Analysis, Modeling and Simulation of Multiscale Problems under AR234/5-2 and GA695/2-2     

Application of the generalized finite difference method to solve the advection-diffusion equation

The study of the advection-diffusion equation continues to be an active field of research. The subject has important applications to fluid dynamics as well as many other branches of science and engineering.This paper shows the application of the generalized finite difference method to solve the advection-diffusion equation by the explicit method. The convergence of the method has been studied and the truncation error over irregular grids is given. An example has been solved using the explicit finite difference formulae and the criterion of stability.     

Convergence properties of the Runge-Kutta-Chebyshev method

Summary The Runge-Kutta-Chebyshev method is ans-stage Runge-Kutta method designed for the explicit integration of stiff systems of ordinary differential equations originating from spatial discretization of parabolic partial differential equations (method of lines). The method possesses an extended real stability interval with a length proportional tos ². The method can be applied withs arbitrarily large, which is an attractive feature due to the proportionality of withs ². The involved stability property here is internal stability. Internal stability has to do with the propagation of errors over the stages within one single integration step. This internal stability property plays an important role in our examination of full convergence properties of a class of 1st and 2nd order schemes. Full convergence means convergence of the fully discrete solution to the solution of the partial differential equation upon simultaneous space-time grid refinement. For a model class of linear problems we prove convergence under the sole condition that the necessary time-step restriction for stability is satisfied. These error bounds are valid for anys and independent of the stiffness of the problem. Numerical examples are given to illustrate the theoretical results.Dedicated to Peter van der Houwen for his numerous contributions in the field of numerical integration of differential equations.Paper presented at the symposium Construction of Stable Numerical Methods for Differential and Integral Equations, held at CWI, March 29, 1989, in honor of Prof. Dr. P.J. van der Houwen to celebrate the twenty-fifth anniversary of his stay at CWI     

Convergence of a finite volume scheme for nonlinear degenerate parabolic equations

Summary. One approximates the entropy weak solution u of a nonlinear parabolic degenerate equation by a piecewise constant function using a discretization in space and time and a finite volume scheme. The convergence of to u is shown as the size of the space and time steps tend to zero. In a first step, estimates on are used to prove the convergence, up to a subsequence, of to a measure valued entropy solution (called here an entropy process solution). A result of uniqueness of the entropy process solution is proved, yielding the strong convergence of to{\it u}. Some on a model equation are shown. Received September 27, 2000 / Published online October 17, 2001     

Multi-dimensional option pricing using radial basis functions and the generalized Fourier transform

We show that the generalized Fourier transform can be used for reducing the computational cost and memory requirements of radial basis function methods for multi-dimensional option pricing. We derive a general algorithm, including a transformation of the Black–Scholes equation into the heat equation, that can be used in any number of dimensions. Numerical experiments in two and three dimensions show that the gain is substantial even for small problem sizes. Furthermore, the gain increases with the number of dimensions.     

A stabilized finite element method for the incompressible magnetohydrodynamic equations

Summary. We propose and analyze a stabilized finite element method for the incompressible magnetohydrodynamic equations. The numerical results that we present show a good behavior of our approximation in experiments which are relevant from an industrial viewpoint. We explain in particular in the proof of our convergence theorem why it may be interesting to stabilize the magnetic equation as soon as the hydrodynamic diffusion is small and even if the magnetic diffusion is large. This observation is confirmed by our numerical tests. Received August 31, 1998 / Revised version received June 16, 1999 / Published online June 21, 2000     

The exponential integrator scheme for stochastic partial differential equations: Pathwise error bounds

We present an error analysis for the pathwise approximation of a general semilinear stochastic evolution equation in d dimensions. We discretise in space by a Galerkin method and in time by using a stochastic exponential integrator. We show that for spatially regular (smooth) noise the number of nodes needed for the noise can be reduced and that the rate of convergence degrades as the regularity of the noise reduces (and the noise becomes rougher).     

Convergence results for an inhomogeneous system arising in various high frequency approximations

Summary. This paper is devoted to both theoretical and numerical study of a system involving an eikonal equation of Hamilton-Jacobi type and a linear conservation law as it comes out of the geometrical optics expansion of the wave equation or the semiclassical limit for the Schr?dinger equation. We first state an existence and uniqueness result in the framework of viscosity and duality solutions. Then we study the behavior of some classical numerical schemes on this problem and we give sufficient conditions to ensure convergence. As an illustration, some practical computations are provided. Received December 6, 1999 / Revised version received August 2, 2000 / Published online June 7, 2001     

IMEX Runge–Kutta schemes for reaction–diffusion equations

A fundamental research is carried out into convergence and stability properties of IMEX (implicit–explicit) Runge–Kutta schemes applied to reaction–diffusion equations. It is shown that a fully discrete scheme converges if it satisfies certain conditions using a technique of the B-convergence analysis, developed by Burrage, Hundsdorfer and Verwer in 1986. Stability of the schemes is also examined on the basis of a scalar test equation, proposed by Frank, Hundsdorfer and Verwer in 1997.     

A convergent scheme for a non local Hamilton Jacobi equation modelling dislocation dynamics

We study dislocation dynamics with a level set point of view. The model we present here looks at the zero level set of the solution of a non local Hamilton Jacobi equation, as a dislocation in a plane of a crystal. The front has a normal speed, depending on the solution itself. We prove existence and uniqueness for short time in the set of continuous viscosity solutions. We also present a first order finite difference scheme for the corresponding level set formulation of the model. The scheme is based on monotone numerical Hamiltonian, proposed by Osher and Sethian. The non local character of the problem makes it not monotone. We obtain an explicit convergence rate of the approximate solution to the viscosity solution. We finally provide numerical simulations.This work has been supported by funds from ACI JC 1041 “Mouvements d’interfaces avec termes non-locaux”, from ACI-JC 1025 “Dynamique des dislocations” and from ONERA, Office National d’Etudes et de Recherches. The second author was also supported by the ENPC-Région Ile de France.     

Error estimate for finite volume scheme

We study the convergence of a finite volume scheme for the linear advection equation with a Lipschitz divergence-free speed in R ^d . We prove a h ^1/2-error estimate in the L ^∞(0,t;L ¹)-norm for BV data. This result was expected from numerical experiments and is optimal.     

A numerical scheme for the Maxwell equations in the quasi-static regime with a non-local source

In this paper we study a non-linear evolution equation, based on quasi-static electromagnetic fields, with a non-local field-dependent source. This model occurs in transformer driven active magnetic shielding. We present a numerical scheme for both time and space discretization and prove convergence of this scheme. We also derive the corresponding error estimates.