Carathéodory approximate solutions for a class of perturbed stochastic differential equations with reflecting boundary

Abstract

In this work, we study the Carathéodory approximate solution for a class of one-dimensional perturbed stochastic differential equations with reflecting boundary (PSDERB). Based on the Carathéodory approximation procedure, we prove that PSDERB have a unique solution and show that the Carathéodory approximate solution converges to the solution of PSDERB whose both drift and diffusion coefficients are non-Lipschitz. After that, we establish an explicit rate of convergence in the case of PSDERB with Lipschitz coefficients.     

A Primal-Dual Active Set Method for Bilaterally Control Constrained Optimal Control of the Navier–Stokes Equations

Abstract

This paper is devoted to the analysis and numerical solution of distributed optimal control of the Navier–Stokes equations in presence of bilateral pointwise control constraints. The analysis of the problem involves the proof of existence of an optimal solution, as well as the presentation of necessary and sufficient conditions for optimality. For the numerical solution of the problem we apply a primal-dual active set strategy and show global and local convergence properties of the method. Finally, some numerical experiments, which illustrate the performance of the method, are discussed.     

Exact Rate of Convergence of Some Approximation Schemes Associated to SDEs Driven by a Fractional Brownian Motion

In this article, we derive the exact rate of convergence of some approximation schemes associated to scalar stochastic differential equations driven by a fractional Brownian motion with Hurst index H. We consider two cases. If H>1/2, the exact rate of convergence of the Euler scheme is determined. We show that the error of the Euler scheme converges almost surely to a random variable, which in particular depends on the Malliavin derivative of the solution. This result extends those contained in J. Complex. 22(4), 459–474, 2006 and C.R. Acad. Sci. Paris, Ser. I 340(8), 611–614, 2005. When 1/6<H<1/2, the exact rate of convergence of the Crank-Nicholson scheme is determined for a particular equation. Here we show convergence in law of the error to a random variable, which depends on the solution of the equation and an independent Gaussian random variable.     

Jacobi Spectral Methods for Volterra-Urysohn Integral Equations of Second Kind with Weakly Singular Kernels

Abstract

In this article, we discuss Jacobi spectral Galerkin and iterated Jacobi spectral Galerkin methods for Volterra-Urysohn integral equations with weakly singular kernels and obtain the convergence results in both the infinity and weighted L²-norm. We show that the order of convergence in iterated Jacobi spectral Galerkin method improves over Jacobi spectral Galerkin method. We obtain the convergence results in two cases when the exact solution is sufficiently smooth and non-smooth. For finding the improved convergence results, we also discuss Jacobi spectral multi-Galerkin and iterated Jacobi spectral multi-Galerkin method and obtain the convergence results in weighted L²-norm. In fact, we prove that the iterated Jacobi spectral multi-Galerkin method improves over iterated Jacobi spectral Galerkin method. We provide numerical results to verify the theoretical results.     

An optimal domain decomposition preconditioner for low-frequency time-harmonic Maxwell equations

The time-harmonic Maxwell equations are considered in the low-frequency case. A finite element domain decomposition approach is proposed for the numerical approximation of the exact solution. This leads to an iteration-by-subdomain procedure, which is proven to converge. The rate of convergence turns out to be independent of the mesh size, showing that the preconditioner implicitly defined by the iterative procedure is optimal. For obtaining this convergence result it has been necessary to prove a regularity theorem for Dirichlet and Neumann harmonic fields.

     

The Convergence of the Navier–Stokes–Poisson System to the Incompressible Euler Equations

ABSTRACT

The combining quasineutral and inviscid limit of the Navier–Stokes–Poisson system in the torus 𝕋 ^d , d ≥ 1 is studied. The convergence of the Navier–Stokes–Poisson system to the incompressible Euler equations is proven for the global weak solution and for the case of general initial data.     

Numerical Solution to Hybrid Stochastic Differential Systems

Abstract

In this article numerical methods for solving hybrid stochastic differential systems of Itô-type are developed by piecewise application of numerical methods for SDEs. We prove a convergence result if the corresponding method for SDEs is numerically stable with uniform convergence in the mean square sense. The Euler and Runge–Kutta methods for hybrid stochastic differential equations are specifically described and the order of the error is given for the Euler method. A numerical example is given to illustrate the theory.     

Oseen Coupling Method for the Exterior Flow Part I: Oseen Coupling Approximation

Abstract In this paper, we consider the bidimensional exterior unsteady Navier-Stokes equations with nonhomogeneous boundary conditions and present an Oseen coupling problem which approximates the Navier-Stokes problem, obtained by coupling the Navier-Stokes equations in the inner region and the Oseen equations in the outer region. Moreover, we prove the existence, uniqueness and the approximate accuracy of the weak solution of the Oseen coupling equations. Project supported by NSF of China & State Major Key Project of Basic Research     

A COMPLETE ANALYSIS FOR SOME A POSTERIORI ERROR ESTIMATES WITH THE FINITE ELEMENT METHOD OF LINES FOR A NONLINEAR PARABOLIC EQUATION

ABSTRACT

A posteriori error estimates for semidiscrete finite element methods for a nonlinear parabolic initial-boundary value problem are considered. The error estimates are obtained by solving local parabolic or elliptic equations for corrections to the solution on each element. The convergence results improve previous results where unnecessary assumptions are imposed on the approximate solution and the elliptic projection of the exact solution.     

Preconditioners Based on Fundamental Solutions

We consider a new preconditioning technique for the iterative solution of linear systems of equations that arise when discretizing partial differential equations. The method is applied to finite difference discretizations, but the ideas apply to other discretizations too. If E is a fundamental solution of a differential operator P, we have E*(Pu) = u. Inspired by this, we choose the preconditioner to be a discretization of an approximate inverse K, given by a convolution-like operator with E as a kernel. We present analysis showing that if P is a first order differential operator, KP is bounded, and numerical results show grid independent convergence for first order partial differential equations, using fixed point iterations. For the second order convection-diffusion equation convergence is no longer grid independent when using fixed point iterations, a result that is consistent with our theory. However, if the grid is chosen to give a fixed number of grid points within boundary layers, the number of iterations is independent of the physical viscosity parameter. AMS subject classification (2000) 65F10, 65N22     

Sur l'existence de solutions d'équations différentielles stochastiques progréssives rétrogrades couplées

Nonlinear BSDEs were first introduced by Pardoux and Peng, 1990, Adapted solutions of backward stochastic differential equations, Systems and Control Letters, 14, 51–61, who proved the existence and uniqueness of a solution under suitable assumptions on the coefficient. Fully coupled forward–backward stochastic differential equations and their connection with PDE have been studied intensively by Pardoux and Tang, 1999, Forward–backward stochastic differential equations and quasilinear parabolic PDE's, Probability Theory and Related Fields, 114, 123–150; Antonelli and Hamadène, 2006, Existence of the solutions of backward–forward SDE's with continuous monotone coefficients, Statistics and Probability Letters, 76, 1559–1569; Hamadème, 1998, Backward–forward SDE's and stochastic differential games, Stochastic Processes and their Applications, 77, 1–15; Delarue, 2002, On the existence and uniqueness of solutions to FBSDEs in a non-degenerate case, Stochastic Processes and Their Applications, 99, 209–286, amongst others.

Unfortunately, most existence or uniqueness results on solutions of forward–backward stochastic differential equations need regularity assumptions. The coefficients are required to be at least continuous which is somehow too strong in some applications. To the best of our knowledge, our work is the first to prove existence of a solution of a forward–backward stochastic differential equation with discontinuous coefficients and degenerate diffusion coefficient where, moreover, the terminal condition is not necessary bounded.

The aim of this work is to find a solution of a certain class of forward–backward stochastic differential equations on an arbitrary finite time interval. To do so, we assume some appropriate monotonicity condition on the generator and drift coefficients of the equation.

The present paper is motivated by the attempt to remove the classical condition on continuity of coefficients, without any assumption as to the non-degeneracy of the diffusion coefficient in the forward equation.

The main idea behind this work is the approximating lemma for increasing coefficients and the comparison theorem. Our approach is inspired by recent work of Boufoussi and Ouknine, 2003, On a SDE driven by a fractional brownian motion and with monotone drift, Electronic Communications in Probability, 8, 122–134; combined with that of Antonelli and Hamadène, 2006, Existence of the solutions of backward–forward SDE's with continuous monotone coefficients, Statistics and Probability Letters, 76, 1559–1569. Pursuing this idea, we adopt a one-dimensional framework for the forward and backward equations and we assume a monotonicity property both for the drift and for the generator coefficient.

At the end of the paper we give some extensions of our result.     

Fourier spectral projection method and nonlinear convergence analysis for Navier-Stokes equations

In this paper, we investigate the convergence rate of the Fourier spectral projection methods for the periodic problem of n-dimensional Navier-Stokes equations. Based on some alternative formulations of the Navier-Stokes equations and the related projection methods, the error estimates are carried out by a global nonlinear error analysis. It simplifies the analysis, relaxes the restriction on the time step size, weakens the regularity requirements on the genuine solution, and leads to some improved convergence results. A new correction technique is proposed for improving the accuracy of the numerical pressure.     

Existence,uniqueness, and numerical approximations for stochastic Burgers equations

Abstract

In this article, we propose an all-in-one statement which includes existence, uniqueness, regularity, and numerical approximations of mild solutions for a class of stochastic partial differential equations (SPDEs) with non-globally monotone nonlinearities. The proof of this result exploits the properties of an existing fully explicit space-time discrete approximation scheme, in particular the fact that it satisfies suitable a priori estimates. We also obtain almost sure and strong convergence of the approximation scheme to the mild solutions of the considered SPDEs. We conclude by applying the main result of the article to the stochastic Burgers equations with additive space-time white noise.     

Stochastic Evolution Equations Driven by a Fractional White Noise

Abstract

In this paper we study stochastic evolution equations driven by a fractional white noise with arbitrary Hurst parameter in infinite dimension. We establish the existence and uniqueness of a mild solution for a nonlinear equation with multiplicative noise under Lipschitz condition by using a fixed point argument in an appropriate inductive limit space. In the linear case with additive noise, a strong solution is obtained. Those results are applied to stochastic parabolic partial differential equations perturbed by a fractional white noise.     

Numerical Solutions of Stochastic Differential Delay Equations with Jumps

Abstract

In this article, we investigate the strong convergence of the Euler–Maruyama method and stochastic theta method for stochastic differential delay equations with jumps. Under a global Lipschitz condition, we not only prove the strong convergence, but also obtain the rate of convergence. We show strong convergence under a local Lipschitz condition and a linear growth condition. Moreover, it is the first time that we obtain the rate of the strong convergence under a local Lipschitz condition and a linear growth condition, i.e., if the local Lipschitz constants for balls of radius R are supposed to grow not faster than log R.     

Splitting methods for second‐order initial value problems

We consider implicit integration methods for the solution of stiff initial value problems for second-order differential equations of the special form y' = f(y). In implicit methods, we are faced with the problem of solving systems of implicit relations. This paper focuses on the construction and analysis of iterative solution methods which are effective in cases where the Jacobian of the right‐hand side of the differential equation can be split into a sum of matrices with a simple structure. These iterative methods consist of the modified Newton method and an iterative linear solver to deal with the linear Newton systems. The linear solver is based on the approximate factorization of the system matrix associated with the linear Newton systems. A number of convergence results are derived for the linear solver in the case where the Jacobian matrix can be split into commuting matrices. Such problems often arise in the spatial discretization of time‐dependent partial differential equations. Furthermore, the stability matrix and the order of accuracy of the integration process are derived in the case of a finite number of iterations. This revised version was published online in August 2006 with corrections to the Cover Date.     

A constrained optimization approach to solving certain systems of convex equations

This research presents a new constrained optimization approach for solving systems of nonlinear equations. Particular advantages are realized when all of the equations are convex. For example, a global algorithm for finding the zero of a convex real-valued function of one variable is developed. If the algorithm terminates finitely, then either the algorithm has computed a zero or determined that none exists; if an infinite sequence is generated, either that sequence converges to a zero or again no zero exists. For solving n-dimensional convex equations, the constrained optimization algorithm has the capability of determining that the system of equations has no solution. Global convergence of the algorithm is established under weaker conditions than previously known and, in this case, the algorithm reduces to Newton’s method together with a constrained line search at each iteration. It is also shown how this approach has led to a new algorithm for solving the linear complementarity problem.     

Quasi-Neutral Limit of the Euler–Poisson and Euler–Monge–Ampère Systems

ABSTRACT

This paper studies the pressureless Euler–Poisson system and its fully nonlinear counterpart, the Euler–Monge–Ampère system, where the fully nonlinear Monge–Ampère equation substitutes for the linear Poisson equation. While the first is a model of plasma physics, the second is derived as a geometric approximation to the Euler incompressible equations. Using energy estimates, convergence of both systems to the Euler incompressible equations is proved.     

Nonlinear Stochastic Difference Equations Driven by Martingales

Abstract

We provide in this paper a systematic development of nonlinear stochastic difference equations driven by martingales (that depend on a spatial parameter); three such equations are considered. We begin with the existence and uniqueness of solutions and continue with the study of stochastic properties, such as the martingale and Markov properties, along with ? irreducibility and recurrence. We discuss in the final section the discrete-time flow and asymptotic flow properties of the solution process.