The search session has expired. Please query the service again.

The search session has expired. Please query the service again.

Displaying 21 – 40 of 53

Showing per page

First order second moment analysis for stochastic interface problems based on low-rank approximation

Helmut Harbrecht, Jingzhi Li (2013)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

In this paper, we propose a numerical method to solve stochastic elliptic interface problems with random interfaces. Shape calculus is first employed to derive the shape-Taylor expansion in the framework of the asymptotic perturbation approach. Given the mean field and the two-point correlation function of the random interface, we can thus quantify the mean field and the variance of the random solution in terms of certain orders of the perturbation amplitude by solving a deterministic elliptic interface...

Fully-discrete finite element approximations for a fourth-order linear stochastic parabolic equation with additive space-time white noise

Georgios T. Kossioris, Georgios E. Zouraris (2010)

ESAIM: Mathematical Modelling and Numerical Analysis

We consider an initial and Dirichlet boundary value problem for a fourth-order linear stochastic parabolic equation, in one space dimension, forced by an additive space-time white noise. Discretizing the space-time white noise a modelling error is introduced and a regularized fourth-order linear stochastic parabolic problem is obtained. Fully-discrete approximations to the solution of the regularized problem are constructed by using, for discretization in space, a Galerkin finite element method...

Implementation of optimal Galerkin and Collocation approximations of PDEs with Random Coefficients⋆⋆⋆

J. Beck, F. Nobile, L. Tamellini, R. Tempone (2011)

ESAIM: Proceedings

In this work we first focus on the Stochastic Galerkin approximation of the solution u of an elliptic stochastic PDE. We rely on sharp estimates for the decay of the coefficients of the spectral expansion of u on orthogonal polynomials to build a sequence of polynomial subspaces that features better convergence properties compared to standard polynomial subspaces such as Total Degree or Tensor Product. We consider then the Stochastic Collocation method, and use the previous estimates to introduce...

Kermack-McKendrick epidemic model revisited

Josef Štěpán, Daniel Hlubinka (2007)

Kybernetika

This paper proposes a stochastic diffusion model for the spread of a susceptible-infective-removed Kermack–McKendric epidemic (M1) in a population which size is a martingale N t that solves the Engelbert–Schmidt stochastic differential equation (). The model is given by the stochastic differential equation (M2) or equivalently by the ordinary differential equation (M3) whose coefficients depend on the size N t . Theorems on a unique strong and weak existence of the solution to (M2) are proved and computer...

Measuring the Irreversibility of Numerical Schemes for Reversible Stochastic Differential Equations

Markos Katsoulakis, Yannis Pantazis, Luc Rey-Bellet (2014)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

For a stationary Markov process the detailed balance condition is equivalent to the time-reversibility of the process. For stochastic differential equations (SDE’s), the time discretization of numerical schemes usually destroys the time-reversibility property. Despite an extensive literature on the numerical analysis for SDE’s, their stability properties, strong and/or weak error estimates, large deviations and infinite-time estimates, no quantitative results are known on the lack of reversibility...

Milstein’s type schemes for fractional SDEs

Mihai Gradinaru, Ivan Nourdin (2009)

Annales de l'I.H.P. Probabilités et statistiques

Weighted power variations of fractional brownian motion B are used to compute the exact rate of convergence of some approximating schemes associated to one-dimensional stochastic differential equations (SDEs) driven by B. The limit of the error between the exact solution and the considered scheme is computed explicitly.

Multiscale Finite Element approach for “weakly” random problems and related issues

Claude Le Bris, Frédéric Legoll, Florian Thomines (2014)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

We address multiscale elliptic problems with random coefficients that are a perturbation of multiscale deterministic problems. Our approach consists in taking benefit of the perturbative context to suitably modify the classical Finite Element basis into a deterministic multiscale Finite Element basis. The latter essentially shares the same approximation properties as a multiscale Finite Element basis directly generated on the random problem. The specific reference method that we use is the Multiscale...

Numerical algorithms for backward stochastic differential equations with 1-d brownian motion: Convergence and simulations***

Shige Peng, Mingyu Xu (2011)

ESAIM: Mathematical Modelling and Numerical Analysis

In this paper we study different algorithms for backward stochastic differential equations (BSDE in short) basing on random walk framework for 1-dimensional Brownian motion. Implicit and explicit schemes for both BSDE and reflected BSDE are introduced. Then we prove the convergence of different algorithms and present simulation results for different types of BSDEs.

Numerical algorithms for backward stochastic differential equations with 1-d brownian motion: Convergence and simulations***

Shige Peng, Mingyu Xu (2011)

ESAIM: Mathematical Modelling and Numerical Analysis

In this paper we study different algorithms for backward stochastic differential equations (BSDE in short) basing on random walk framework for 1-dimensional Brownian motion. Implicit and explicit schemes for both BSDE and reflected BSDE are introduced. Then we prove the convergence of different algorithms and present simulation results for different types of BSDEs.

Numerical analysis of parallel replica dynamics

Gideon Simpson, Mitchell Luskin (2013)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit distribution from a given well for a single process can be approximated by...

Numerical schemes for multivalued backward stochastic differential systems

Lucian Maticiuc, Eduard Rotenstein (2012)

Open Mathematics

We define approximation schemes for generalized backward stochastic differential systems, considered in the Markovian framework. More precisely, we propose a mixed approximation scheme for the following backward stochastic variational inequality: d Y t + F ( t , X t , Y t , Z t ) d t φ ( Y t ) d t + Z t d W t , where ∂φ is the subdifferential operator of a convex lower semicontinuous function φ and (X t)t∈[0;T] is the unique solution of a forward stochastic differential equation. We use an Euler type scheme for the system of decoupled forward-backward variational...

On surrogate learning for linear stability assessment of Navier-Stokes equations with stochastic viscosity

Bedřich Sousedík, Howard C. Elman, Kookjin Lee, Randy Price (2022)

Applications of Mathematics

We study linear stability of solutions to the Navier-Stokes equations with stochastic viscosity. Specifically, we assume that the viscosity is given in the form of a stochastic expansion. Stability analysis requires a solution of the steady-state Navier-Stokes equation and then leads to a generalized eigenvalue problem, from which we wish to characterize the real part of the rightmost eigenvalue. While this can be achieved by Monte Carlo simulation, due to its computational cost we study three surrogates...

On the convergence of generalized polynomial chaos expansions

Oliver G. Ernst, Antje Mugler, Hans-Jörg Starkloff, Elisabeth Ullmann (2012)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

A number of approaches for discretizing partial differential equations with random data are based on generalized polynomial chaos expansions of random variables. These constitute generalizations of the polynomial chaos expansions introduced by Norbert Wiener to expansions in polynomials orthogonal with respect to non-Gaussian probability measures. We present conditions on such measures which imply mean-square convergence of generalized polynomial chaos expansions to the correct limit and complement...

On the convergence of generalized polynomial chaos expansions

Oliver G. Ernst, Antje Mugler, Hans-Jörg Starkloff, Elisabeth Ullmann (2011)

ESAIM: Mathematical Modelling and Numerical Analysis

A number of approaches for discretizing partial differential equations with random data are based on generalized polynomial chaos expansions of random variables. These constitute generalizations of the polynomial chaos expansions introduced by Norbert Wiener to expansions in polynomials orthogonal with respect to non-Gaussian probability measures. We present conditions on such measures which imply mean-square convergence of generalized polynomial...

On the convergence of the stochastic Galerkin method for random elliptic partial differential equations

Antje Mugler, Hans-Jörg Starkloff (2013)

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique

In this article we consider elliptic partial differential equations with random coefficients and/or random forcing terms. In the current treatment of such problems by stochastic Galerkin methods it is standard to assume that the random diffusion coefficient is bounded by positive deterministic constants or modeled as lognormal random field. In contrast, we make the significantly weaker assumption that the non-negative random coefficients can be bounded strictly away from zero and infinity by random...

Currently displaying 21 – 40 of 53