Volume 31, Issue 1

Backward stochastic differential equations (BSDEs) belong nowadays to the most frequently studied equations in stochastic analysis and computational stochastics. BSDEs in applications are often nonlinear and high-dimensional. In nearly all cases such nonlinear high-dimensional BSDEs cannot be solved explicitly and it has been and still is a very active topic of research to design and analyze numerical approximation methods to approximatively solve nonlinear high-dimensional BSDEs. Although there are a large number of research articles in the scientific literature which analyze numerical approximation methods for nonlinear BSDEs, until today there has been no numerical approximation method in the scientific literature which has been proven to overcome the curse of dimensionality in the numerical approximation of nonlinear BSDEs in the sense that the number of computational operations of the numerical approximation method to approximatively compute one sample path of the BSDE solution grows at most polynomially in both the reciprocal 1/ε of the prescribed approximation accuracy ε ∈ (0, ∞) and the dimension d ∈ N = {1, 2, 3, . . .} of the BSDE. It is the key contribution of this article to overcome this obstacle by introducing a new Monte Carlo-type numerical approximation method for high-dimensional BSDEs and by proving that this Monte Carlo-type numerical approximation method does indeed overcome the curse of dimensionality in the approximative computation of solution paths of BSDEs.

In this note we consider a parabolic evolution equation in a polygonal space-time cylinder. We show, that the elliptic part is given by a m -accretive mapping from L q ( Ω ) → L q ( Ω ). Therefore we can apply the theory of nonlinear semigroups in Banach spaces in order to get regularity results in time and space. The second part of the paper deals with the numerical solution of the problem. It is dedicated to the application of the space-time discontinuous Galerkin method (STDGM). It means that both in space as well as in time discontinuous piecewise polynomial approximations of the solution are used. We concentrate to the theoretical analysis of the error estimation.

The convergence analysis with rates for adaptive least-squares finite element methods (ALSFEMs) combines arguments from the a posteriori analysis of conforming and mixed finite element schemes. This paper provides an overview of the key arguments for the verification of the axioms of adaptivity for an ALSFEM for the solution of a linear model problem. The formulation at hand allows for the simultaneous analysis of first-order systems of the Poisson model problem, the Stokes equations, and the linear elasticity equations. Following [Carstensen and Park, SIAM J. Numer. Anal. 53 (1), 2015], the adaptive algorithm is driven by an alternative residual-based error estimator with exact solve and includes a separate marking strategy for quasi-optimal data resolution of the right-hand side. This presentation covers conforming discretisations for an arbitrary polynomial degree and mixed homogeneous boundary conditions.

We present a loosely coupled, non-iterative time-splitting scheme based on Robin–Robin coupling conditions. We apply a novel unified analysis for this scheme applied to both a parabolic/parabolic coupled system and a parabolic/hyperbolic coupled system. We show for both systems that the scheme is stable, and the error converges as O(ΔtT+log(1Δt)), $\mathcal{O}\big({\Delta t} \sqrt{T +\log(\frac{1}{{\Delta t}})}\big),$where Δt is the time step.