In this article, we study the numerical approximation of stochastic differential equations driven by a multidimensional fractional Brownian motion (fBm) with Hurst parameter greater than 1/3. We introduce an implementable scheme for these equations, which is based on a second-order Taylor expansion, where the usual Lévy area terms are replaced by products of increments of the driving fBm. The convergence of our scheme is shown by means of a combination of rough paths techniques and error bounds...

In this paper, we propose a method for the approximation of the solution of high-dimensional weakly coercive problems formulated in tensor spaces using low-rank approximation formats. The method can be seen as a perturbation of a minimal residual method with a measure of the residual corresponding to the error in a specified solution norm. The residual norm can be designed such that the resulting low-rank approximations are optimal with respect to particular norms of interest, thus allowing to take...

This special volume of the ESAIM Journal, Mathematical Modelling and Numerical Analysis, contains a collection of articles on probabilistic interpretations of some classes of nonlinear integro-differential equations. The selected contributions deal with a wide range of topics in applied probability theory and stochastic analysis, with applications in a variety of scientific disciplines, including physics, biology, fluid mechanics, molecular chemistry, financial mathematics and bayesian statistics....

We are interested in the approximation and simulation of solutions for the backward stochastic differential equations. We suggest two approximation schemes, and we study the ${𝕃}^2$ induced error.

This work is concerned with the asymptotic analysis of a time-splitting scheme for the Schrödinger equation with a random potential having weak amplitude, fast oscillations in time and space, and long-range correlations. Such a problem arises for instance in the simulation of waves propagating in random media in the paraxial approximation. The high-frequency limit of the Schrödinger equation leads to different regimes depending on the distance of propagation, the oscillation pattern of the initial...

This paper describes the extension of a recently developed numerical solver for the Landau-Lifshitz Navier-Stokes (LLNS) equations to binary mixtures in three dimensions. The LLNS equations incorporate thermal fluctuations into macroscopic hydrodynamics by using white-noise fluxes. These stochastic PDEs are more complicated in three dimensions due to the tensorial form of the correlations for the stochastic fluxes and in mixtures due to couplings of energy and concentration fluxes (e.g., Soret...

Under the key assumption of finite $\rho$-variation, $\rho \in [1,2)$, of the covariance of the underlying Gaussian process, sharp a.s. convergence rates for approximations of Gaussian rough paths are established. When applied to Brownian resp. fractional Brownian motion (fBM), $\rho =1$ resp. $\rho =1/\left(2H\right)$, we recover and extend the respective results of (Trans. Amer. Math. Soc.361 (2009) 2689–2718) and (Ann. Inst. Henri Poincasé Probab. Stat.48(2012) 518–550). In particular, we establish an a.s. rate $k^{-(1/\rho -1/2-\epsilon )}$, any $\epsilon gt;0$, for Wong–Zakai and Milstein-type...

We study $L^p$ convergence for the Euler scheme for stochastic differential equations reflecting on the boundary of a general convex domain D ⊆ ℝd. We assume that the equation has the pathwise uniqueness property and its coefficients are measurable and continuous almost everywhere with respect to the Lebesgue measure. In the case D=[0,∞) new sufficient conditions ensuring pathwise uniqueness for equations with possibly discontinuous coefficients are given.

Much of uncertainty quantification to date has focused on determining the effect of variables modeled probabilistically, and with a known distribution, on some physical or engineering system. We develop methods to obtain information on the system when the distributions of some variables are known exactly, others are known only approximately, and perhaps others are not modeled as random variables at all.The main tool used is the duality between risk-sensitive integrals and relative entropy, and we...

The goal of this contribution is to introduce some approaches to uncertainty modeling in a way accessible to non-specialists. Elements of the Monte Carlo method, polynomial chaos method, Dempster-Shafer approach, fuzzy set theory, and the worst (case) scenario method are presented.

This paper discusses analytical and numerical issues related to elliptic equations with random coefficients which are generally nonlinear functions of white noise. Singularity issues are avoided by using the Itô-Skorohod calculus to interpret the interactions between the coefficients and the solution. The solution is constructed by means of the Wiener Chaos (Cameron-Martin) expansions. The existence and uniqueness of the solutions are established under rather weak assumptions, the main of which...

We study ergodic properties of stochastic geometric wave equations on a particular model with the target being the 2D sphere while considering only solutions which are independent of the space variable. This simplification leads to a degenerate stochastic equation in the tangent bundle of the 2D sphere. Studying this equation, we prove existence and non-uniqueness of invariant probability measures for the original problem and obtain also results on attractivity towards an invariant measure. We also...

We consider one-dimensional stochastic differential equations in the particular case of diffusion coefficient functions of the form ${\left|x\right|}^{\alpha }$, $\alpha \in [1/2,1)$. In that case, we study the rate of convergence of a symmetrized version of the Euler scheme. This symmetrized version is easy to simulate on a computer. We prove its strong convergence and obtain the same rate of convergence as when the coefficients are Lipschitz.

We consider one-dimensional stochastic differential equations in the particular case of diffusion coefficient functions of the form |x|α, α ∈ [1/2,1). In that case, we study the rate of convergence of a symmetrized version of the Euler scheme. This symmetrized version is easy to simulate on a computer. We prove its strong convergence and obtain the same rate of convergence as when the coefficients are Lipschitz.

In this note we propose an exact simulation algorithm for the solution of (1)$$\mathrm{d}{X}_t=\mathrm{d}{W}_t+\overline{b}\left(X_t\right)\mathrm{d}t,X_0=x,$$d X t = d W t + b̅ ( X t ) d t,   X 0 = x, where $\overline{b}$b̅is a smooth real function except at point 0 where $\overline{b}(0+)\ne \overline{b}(0-)$b̅(0 + ) ≠ b̅(0 −) . The main idea is to sample an exact skeleton of Xusing an algorithm deduced from the convergence of the solutions of the skew perturbed equation (2)$$\mathrm{d}{X}_t^{\beta }=\mathrm{d}{W}_t+\overline{b}\left(X_t^{\beta }\right)\mathrm{d}t+\beta \mathrm{d}{L}_t^0\left(X^{\beta }\right),X_0=x$$d X t β = d W t + b̅ ( X t β ) d t + β d L t 0 ( X β ) ,   X 0 = x towardsX solution of (1) as β ≠ 0 tends to 0. In this note, we show that this convergence...

We study a free energy computation procedure, introduced in [Darve and Pohorille, J. Chem. Phys.115 (2001) 9169–9183; Hénin and Chipot, J. Chem. Phys.121 (2004) 2904–2914], which relies on the long-time behavior of a nonlinear stochastic differential equation. This nonlinearity comes from a conditional expectation computed with respect to one coordinate of the solution. The long-time convergence of the solutions to this equation has been proved in [Lelièvre et al., Nonlinearity21 (2008) 1155–1181],...