The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

For stochastic differential equations (SDEs) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient, the classical explicit Euler scheme fails to converge strongly to the exact solution. Recently, an explicit strongly convergent numerical scheme, called the tamed Euler method, has been proposed in [8] for such SDEs. Motivated by their work, we here introduce a tamed version of the Milstein scheme for SDEs with commutative noise. The proposed method is also explicit and easily implementable, but achieves higher strong convergence order than the tamed Euler method does. In recovering the strong convergence order one of the new method, new difficulties arise and kind of a bootstrap argument is developed to overcome them. Finally, an illustrative example confirms the computational efficiency of the tamed Milstein method compared to the tamed Euler method.     

随机微分方程欧拉格式算法分析

首先给出了线性随机微分方程的欧拉格式算法,然后给出了非线性随机微分方程变步长的欧拉格式算法,接着讨论了其对初值的连续依赖性和收敛性.     

An Energy Stable Filtered Backward Euler Scheme for the MBE Equation with Slope Selection

As a promising strategy to adjust the order in the variable-order BDF algorithm, a time filtered backward Euler scheme is investigated for the molecular beam epitaxial equation with slope selection. The temporal second-order convergence in the $L^2$ norm is established under a convergence-solvability-stability (CSS)-consistent time-step constraint. The CSS-consistent condition means that the maximum step-size limit required for convergence is of the same order to that for solvability and stability (in certain norms) as the small interface parameter $ε → 0^+.$ Similar to the backward Euler scheme, the time filtered backward Euler scheme preserves some physical properties of the original problem at the discrete levels, including the volume conservation, the energy dissipation law and $L^2$ norm boundedness. Numerical tests are included to support the theoretical results.     

Weak Euler Approximation for Itô Diffusion and Jump Processes

This article studies the rate of convergence of the weak Euler approximation for Itô diffusion and jump processes with Hölder-continuous generators. It covers a number of stochastic processes including the nondegenerate diffusion processes and a class of stochastic differential equations driven by stable processes. To estimate the rate of convergence, the existence of a unique solution to the corresponding backward Kolmogorov equation in Hölder space is first proved. It then shows that the Euler scheme yields positive weak order of convergence.     

On the rate of convergence of simple and jump-adapted weak Euler schemes for Lévy driven SDEs

The paper studies the rate of convergence of a weak Euler approximation for solutions to possibly completely degenerate SDEs driven by Lévy processes, with Hölder-continuous coefficients. It investigates the dependence of the rate on the regularity of coefficients and driving processes and its robustness to the approximation of the increments of the driving process. A convergence rate is derived for some approximate jump-adapted Euler scheme as well.     

A difference equation approach to parameter estimation for differential-delay equations

A finite difference scheme for estimating parameters in linear differential-delay equations is investigated. Convergence results and rates of convergence are obtained for a simple explicit (Euler's) method. Numerical examples are given to illustrate the convergence for the Euler method. Numerical results for a “higher order” scheme are also discussed.     

The Euler Scheme for Feller Processes

We consider the Euler scheme for stochastic differential equations with jumps, whose intensity might be infinite and the jump structure may depend on the position. This general type of SDE is explicitly given for Feller processes and a general convergence condition is presented.

In particular, the characteristic functions of the increments of the Euler scheme are calculated in terms of the symbol of the Feller process in a closed form. These increments are increments of Lévy processes and, thus, the Euler scheme can be used for simulation by applying standard techniques from Lévy processes.     

Stochastic Theta Method for a Reflected Stochastic Differential Equation

In this article, a stochastic theta method for a reflected stochastic differential equation is proposed. When the parameter θ = 0, this method coincides with the projection Euler scheme; while when the parameter θ = 1, it is called an implicit projection Euler scheme which is first proposed in this article. Under some conditions, the strong convergence and the A-stability of this numerical scheme are proved.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

A Finite Volume Implicit Euler Scheme for the Linearized Shallow Water Equations: Stability and Convergence

The linearized shallow water equations are discretized in space by a finite volume method and in time by an implicit Euler scheme. Stability and convergence of the scheme are proved.     

A note on the convergence of discretized dynamic iteration

In this note we prove convergence results, including error estimates, for the dynamic iteration scheme where the forward Euler and backward Euler method are used to compute the iterates. The proofs are interesting in that they are exact analogues of the proof for the continuous case, using discrete versions of Gronwall's inequality.     

Error estimates for the finite volume discretization for the porous medium equation

We analyze the convergence of a numerical scheme for a class of degenerate parabolic problems modelling reactions in porous media, and involving a nonlinear, possibly vanishing diffusion. The scheme involves the Kirchhoff transformation of the regularized nonlinearity, as well as an Euler implicit time stepping and triangle based finite volumes. We prove the convergence of the approach by giving error estimates in terms of the discretization and regularization parameter.     

A Lagrangian Scheme à la Brenier for the Incompressible Euler Equations

We approximate the regular solutions of the incompressible Euler equations by the solution of ODEs on finite-dimensional spaces. Our approach combines Arnold’s interpretation of the solution of the Euler equations for incompressible and inviscid fluids as geodesics in the space of measure-preserving diffeomorphisms, and an extrinsic approximation of the equations of geodesics due to Brenier. Using recently developed semi-discrete optimal transport solvers, this approach yields a numerical scheme which is able to handle problems of realistic size in 2D. Our purpose in this article is to establish the convergence of this scheme towards regular solutions of the incompressible Euler equations, and to provide numerical experiments on a few simple test cases in 2D.     

Rate of convergence in the Euler scheme for stochastic differential equations with non-Lipschitz diffusion and Poisson measure

We study the rate of convergence and some other properties of the Euler scheme for stochastic differential equations with non-Lipschitz diffusion and Poisson measure.     

On the convergence of a regularizing Levenberg–Marquardt scheme for nonlinear ill-posed problems

In this note, we study the convergence of the Levenberg–Marquardt regularization scheme for nonlinear ill-posed problems. We consider the case that the initial error satisfies a source condition. Our main result shows that if the regularization parameter does not grow too fast (not faster than a geometric sequence), then the scheme converges with optimal convergence rates. Our analysis is based on our recent work on the convergence of the exponential Euler regularization scheme (Hochbruck et al. in Inverse Probl 25(7):075009, 2009).     

Strong Convergence of a Fully Discrete Scheme for Multiplicative Noise Driving SPDEs with Non-Globally Lipschitz Continuous Coefficients

This work investigates strong convergence of numerical schemes for nonlinear multiplicative noise driving stochastic partial differential equations under some weaker conditions imposed on the coefficients avoiding the commonly used global Lipschitz assumption in the literature. Space-time fully discrete scheme is proposed, which is performed by the finite element method in space and the implicit Euler method in time. Based on some technical lemmas including regularity properties for the exact solution of the considered problem, strong convergence analysis with sharp convergence rates for the proposed fully discrete scheme is rigorously established.     

A random Euler scheme for Carathéodory differential equations

We study a random Euler scheme for the approximation of Carathéodory differential equations and give a precise error analysis. In particular, we show that under weak assumptions, this approximation scheme obtains the same rate of convergence as the classical Monte–Carlo method for integration problems.     

On the uniform convergence of a finite difference scheme for time dependent singularly perturbed reaction-diffusion problems

In this work we are interested in the numerical approximation of 1D parabolic singularly perturbed problems of reaction-diffusion type. To approximate the multiscale solution of this problem we use a numerical scheme combining the classical backward Euler method and central differencing. The scheme is defined on some special meshes which are the tensor product of a uniform mesh in time and a special mesh in space, condensing the mesh points in the boundary layer regions. In this paper three different meshes of Shishkin, Bahkvalov and Vulanovic type are used, proving the uniform convergence with respect to the diffusion parameter. The analysis of the uniform convergence is based on a new study of the asymptotic behavior of the solution of the semidiscrete problems, which are obtained after the time discretization by the Euler method. Some numerical results are showed corroborating in practice the theoretical results on the uniform convergence and the order of the method.     

Product integration methods for solving a system of nonlinear Volterra integral equations

In this paper the technique of subtracting out singularities is used to derive explicit and implicit product Euler schemes with order one convergence and a product trapezoidal scheme with order two convergence for a system of Volterra integral equations with a weakly singular kernel. The convergence proofs of the numerical schemes are presented; these are nonstandard since the nonlinear function involved in the integral equation system does not satisfy a global Lipschitz condition.     

LONG-TIME CONVERGENCE OF GENERALIZED DIFFERENCE METHOD FOR NAVIER-STOKES EQUATIONS

1　IntroductionIn this paper,we firstprovide a generalized difference method for the two-dimension-al Navier-Stokes equations by combining the ideas of staggered scheme[6] and generalizedupwind scheme [4 ] in space,and by backward Euler time-stepping.Then we apply theabstractframework of[7] to prove its long-time convergence.The outline of this paper is as follows:In§ 2 we state the generalized differencemethod.In§ 3 we provide some lemmas.In§ 4 we study the one-sided Lipschitz condi-tio…