Runge–Kutta methods for jump–diffusion differential equations

In this paper we consider Runge–Kutta methods for jump–diffusion differential equations. We present a study of their mean-square convergence properties for problems with multiplicative noise. We are concerned with two classes of Runge–Kutta methods. First, we analyse schemes where the drift is approximated by a Runge–Kutta ansatz and the diffusion and jump part by a Maruyama term and second we discuss improved methods where mixed stochastic integrals are incorporated in the approximation of the next time step as well as the stage values of the Runge–Kutta ansatz for the drift. The second class of methods are specifically developed to improve the accuracy behaviour of problems with small noise. We present results showing when the implicit stochastic equations defining the stage values of the Runge–Kutta methods are uniquely solvable. Finally, simulation results illustrate the theoretical findings.     

A-stability of Runge-Kutta methods for systems with additive noise

Numerical stability of both explicit and implicit Runge-Kutta methods for solving ordinary differential equations with an additive noise term is studied. The concept of numerical stability of deterministic schemes is extended to the stochastic case, and a stochastic analogue of Dahlquist'sA-stability is proposed. It is shown that the discretization of the drift term alone controls theA-stability of the whole scheme. The quantitative effect of implicitness uponA-stability is also investigated, and stability regions are given for a family of implicit Runge-Kutta methods with optimal order of convergence.This author was partially supported by the Italian Consiglio Nazionale delle Ricerche.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

Spatial approximation of stochastic convolutions

We study linear stochastic evolution partial differential equations driven by additive noise. We present a general and flexible framework for representing the infinite dimensional Wiener process, which drives the equation. Since the eigenfunctions and eigenvalues of the covariance operator of the process are usually not available for computations, we propose an expansion in an arbitrary frame. We show how to obtain error estimates when the truncated expansion is used in the equation. For the stochastic heat and wave equations, we combine the truncated expansion with a standard finite element method and derive a priori bounds for the mean square error. Specializing the frame to biorthogonal wavelets in one variable, we show how the hierarchical structure, support and cancelation properties of the primal and dual bases lead to near sparsity and can be used to simplify the simulation of the noise and its update when new terms are added to the expansion.     

A new approach to stochastic evolution equations with adapted drift

In this paper we develop a new approach to stochastic evolution equations with an unbounded drift A which is dependent on time and the underlying probability space in an adapted way. It is well-known that the semigroup approach to equations with random drift leads to adaptedness problems for the stochastic convolution term. In this paper we give a new representation formula for the stochastic convolution which avoids integration of non-adapted processes. Here we mainly consider the parabolic setting. We establish connections with other solution concepts such as weak solutions. The usual parabolic regularity properties are derived and we show that the new approach can be applied in the study of semilinear problems with random drift. At the end of the paper the results are illustrated with two examples of stochastic heat equations with random drift.     

Convergence and stability of implicit runge-kutta methods for systems with multiplicative noise

A class ofimplicit Runge-Kutta schemes for stochastic differential equations affected bymultiplicative Gaussian white noise is shown to be optimal with respect to global order of convergence in quadratic mean. A test equation is proposed in order to investigate the stability of discretization methods for systems of this kind. Herestability is intended in a truly probabilistic sense, as opposed to the recently introduced extension of A-stability to the stochastic context, given for systems with additive noise. Stability regions for the optimal class are also given.Partially supported by the Italian Consiglio Nazionale delle Ricerche.     

Stochastic optimization theory of backward stochastic differential equations with jumps and viscosity solutions of Hamilton–Jacobi–Bellman equations

In this paper we study stochastic optimal control problems with jumps with the help of the theory of Backward Stochastic Differential Equations (BSDEs) with jumps. We generalize the results of Peng [S. Peng, BSDE and stochastic optimizations, in: J. Yan, S. Peng, S. Fang, L. Wu, Topics in Stochastic Analysis, Science Press, Beijing, 1997 (Chapter 2) (in Chinese)] by considering cost functionals defined by controlled BSDEs with jumps. The application of BSDE methods, in particular, the use of the notion of stochastic backward semigroups introduced by Peng in the above-mentioned work allows a straightforward proof of a dynamic programming principle for value functions associated with stochastic optimal control problems with jumps. We prove that the value functions are the viscosity solutions of the associated generalized Hamilton–Jacobi–Bellman equations with integral-differential operators. For this proof, we adapt Peng’s BSDE approach, given in the above-mentioned reference, developed in the framework of stochastic control problems driven by Brownian motion to that of stochastic control problems driven by Brownian motion and Poisson random measure.     

A Lax equivalence theorem for stochastic differential equations

In this paper, a stochastic mean square version of Lax’s equivalence theorem for Hilbert space valued stochastic differential equations with additive and multiplicative noise is proved. Definitions for consistency, stability, and convergence in mean square of an approximation of a stochastic differential equation are given and it is shown that these notions imply similar results as those known for approximations of deterministic partial differential equations. Examples show that the assumptions made are met by standard approximations.     

Mean-square convergence of stochastic multi-step methods with variable step-size

We study mean-square consistency, stability in the mean-square sense and mean-square convergence of drift-implicit linear multi-step methods with variable step-size for the approximation of the solution of Itô stochastic differential equations. We obtain conditions that depend on the step-size ratios and that ensure mean-square convergence for the special case of adaptive two-step-Maruyama schemes. Further, in the case of small noise we develop a local error analysis with respect to the h$h$–ε$\epsilon$ approach and we construct some stochastic linear multi-step methods with variable step-size that have order 2 behaviour if the noise is small enough.     

Optimal expulsion and optimal confinement of a Brownian particle with a switching cost

We solve two stochastic control problems in which a player tries to minimize or maximize the exit time from an interval of a Brownian particle, by controlling its drift. The player can change from one drift to another but is subject to a switching cost. In each problem, the value function is written as the solution of a free boundary problem involving second order ordinary differential equations, in which the unknown boundaries are found by applying the principle of smooth fit. For both problems, we compute the value function, we exhibit the optimal strategy and we prove its generic uniqueness.     

Improved linear multi-step methods for stochastic ordinary differential equations

We consider linear multi-step methods for stochastic ordinary differential equations and study their convergence properties for problems with small noise or additive noise. We present schemes where the drift part is approximated by well-known methods for deterministic ordinary differential equations. In previous work, we considered Maruyama-type schemes, where only the increments of the driving Wiener process are used to discretize the diffusion part. Here, we suggest the improvement of the discretization of the diffusion part by also taking into account mixed classical-stochastic integrals. We show that the relation of the applied step sizes to the smallness of the noise is essential in deciding whether the new methods are worthwhile. Simulation results illustrate the theoretical findings.     

Three-stage stochastic Runge–Kutta methods for stochastic differential equations

In this paper we discuss three-stage stochastic Runge–Kutta (SRK) methods with strong order 1.0 for a strong solution of Stratonovich stochastic differential equations (SDEs). Higher deterministic order is considered. Two methods, a three-stage explicit (E3) method and a three-stage semi-implicit (SI3) method, are constructed in this paper. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of several standard test problems.     

Exponential stability of numerical solutions for a class of stochastic age-dependent capital system with Poisson jumps

Recently, numerical solutions of stochastic differential equations have received a great deal of attention. Numerical approximation schemes are invaluable tools for exploring their properties. In this paper, we introduce a class of stochastic age-dependent (vintage) capital system with Poisson jumps. We also give the discrete approximate solution with an implicit Euler scheme in time discretization. Using Gronwall’s lemma and Barkholder-Davis-Gundy’s inequality, some criteria are obtained for the exponential stability of numerical solutions to the stochastic age-dependent capital system with Poisson jumps. It is proved that the numerical approximation solutions converge to the analytic solutions of the equations under the given conditions, where information on the order of approximation is provided. These error bounds imply strong convergence as the timestep tends to zero. A numerical example is used to illustrate the theoretical results.     

The role of coefficients of a general SPDE on the stability and convergence of a finite difference method

In this paper for the approximate solution of stochastic partial differential equations (SPDEs) of Itô-type, the stability and application of a class of finite difference method with regard to the coefficients in the equations is analyzed. The finite difference methods discussed here will be either explicit or implicit and a comparison between them will be reported. We prove the consistency and stability of these methods and investigate the influence of the multiplier (particularly multiplier of the random noise) in mean square stability. From stochastic version of Lax-Richtmyer the convergence of these methods under some conditions are established. Numerical experiments are included to show the efficiency of the methods.     

Strong convergence of split-step backward Euler method for stochastic differential equations with non-smooth drift

In this paper, we are concerned with the numerical approximation of stochastic differential equations with discontinuous/nondifferentiable drifts. We show that under one-sided Lipschitz and general growth conditions on the drift and global Lipschitz condition on the diffusion, a variant of the implicit Euler method known as the split-step backward Euler (SSBE) method converges with strong order of one half to the true solution. Our analysis relies on the framework developed in [D. J. Higham, X. Mao and A. M. Stuart, Strong convergence of Euler-type methods for nonlinear stochastic differential equations, SIAM Journal on Numerical Analysis, 40 (2002) 1041-1063] and exploits the relationship which exists between explicit and implicit Euler methods to establish the convergence rate results.     

Existence of global mild and strong solutions to stochastic hyperbolic evolution equations driven by a spatially homogeneous Wiener process

Semilinear second order stochastic hyperbolic equations driven by a spatially homogeneous Wiener process are studied. Sufficient conditions in terms of Lyapunov functions for the equation to have global mild or strong solutions are found. In particular, the results apply to equations with polynomial drift and diffusion coefficients.     

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

Stochastic generalized porous media and fast diffusion equations

We present a generalization of Krylov-Rozovskii's result on the existence and uniqueness of solutions to monotone stochastic differential equations. As an application, the stochastic generalized porous media and fast diffusion equations are studied for σ-finite reference measures, where the drift term is given by a negative definite operator acting on a time-dependent function, which belongs to a large class of functions comparable with the so-called N-functions in the theory of Orlicz spaces.     

Stochastic Cahn–Hilliard equation with singular nonlinearity and reflection

We consider a stochastic partial differential equation with logarithmic (or negative power) nonlinearity, with one reflection at 0 and with a constraint of conservation of the space average. The equation, driven by the derivative in space of a space–time white noise, contains a bi-Laplacian in the drift. The lack of the maximum principle for the bi-Laplacian generates difficulties for the classical penalization method, which uses a crucial monotonicity property. Being inspired by the works of Debussche and Zambotti, we use a method based on infinite dimensional equations, approximation by regular equations and convergence of the approximated semigroup. We obtain existence and uniqueness of a solution for nonnegative initial conditions, results on the invariant measures, and on the reflection measures.     

A new adaptive Runge–Kutta method for stochastic differential equations

In this paper, we will present a new adaptive time stepping algorithm for strong approximation of stochastic ordinary differential equations. We will employ two different error estimation criteria for drift and diffusion terms of the equation, both of them based on forward and backward moves along the same time step. We will use step size selection mechanisms suitable for each of the two main regimes in the solution behavior, which correspond to domination of the drift-based local error estimator or diffusion-based one. Numerical experiments will show the effectiveness of this approach in the pathwise approximation of several standard test problems.