BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Charles-Edouard Bréhier

Charles-Edouard Bréhier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.PRmath.NANumerical AnalysisComputationcond-mat.stat-mechphysics.comp-phphysics.data-an

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Analysis of a modified Euler scheme for parabolic semilinear stochastic PDEs

We propose a modification of the standard linear implicit Euler integrator for the weak approximation of parabolic semilinear stochastic PDEs driven by additive space-time white noise. The new method can easily be combined with a finite difference method for the spatial discretization. The proposed method is shown to have improved qualitative properties compared with the standard method. First, for any time-step size, the spatial regularity of the solution is preserved, at all times. Second, the proposed method preserves the Gaussian invariant distribution of the infinite dimensional Ornstein--Uhlenbeck process obtained when the nonlinearity is absent, for any time-step size. The weak order of convergence of the proposed method is shown to be equal to $1/2$ in a general setting, like for the standard Euler scheme. A stronger weak approximation result is obtained when considering the approximation of a Gibbs invariant distribution, when the nonlinearity is a gradient: one obtains an approximation in total variation distance of order $1/2$, which does not hold for the standard method. This is the first result of this type in the literature. A key point in the analysis is the interpretation of the proposed modified Euler scheme as the accelerated exponential Euler scheme applied to a modified stochastic evolution equation. Finally, it is shown that the proposed method can be applied to design an asymptotic preserving scheme for a class of slow-fast multiscale systems, and to construct a Markov Chain Monte Carlo method which is well-defined in infinite dimension. We also revisit the analysis of the standard and the accelerated exponential Euler scheme, and we prove new results with approximation in the total variation distance, which serve to illustrate the behavior of the proposed modified Euler scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Splitting schemes for FitzHugh--Nagumo stochastic partial differential equations

We design and study splitting integrators for the temporal discretization of the stochastic FitzHugh--Nagumo system. This system is a model for signal propagation in nerve cells where the voltage variable is solution of a one-dimensional parabolic PDE with a cubic nonlinearity driven by additive space-time white noise. We first show that the numerical solutions have finite moments. We then prove that the splitting schemes have, at least, the strong rate of convergence $1/4$. Finally, numerical experiments illustrating the performance of the splitting schemes are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Uniform error bounds for numerical schemes applied to multiscale SDEs in a Wong-Zakai diffusion approximation regime

We study a family of numerical schemes applied to a class of multiscale systems of stochastic differential equations. When the time scale separation parameter vanishes, a well-known homogenization or Wong--Zakai diffusion approximation result states that the slow component of the considered system converges to the solution of a stochastic differential equation driven by a real-valued Wiener process, with Stratonovich interpretation of the noise. We propose and analyse schemes for effective approximation of the slow component. Such schemes satisfy an asymptotic preserving property and generalize the methods proposed in a recent article. We fill a gap in the analysis of these schemes and prove strong error estimates, which are uniform with respect to the time scale separation parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Uniform strong and weak error estimates for numerical schemes applied to multiscale SDEs in a Smoluchowski-Kramers diffusion approximation regime

We study a family of numerical schemes applied to a class of multiscale systems of stochastic differential equations. When the time scale separation parameter vanishes, a well-known Smoluchowski--Kramers diffusion approximation result states that the slow component of the considered system converges to the solution of a standard Itô stochastic differential equation. We propose and analyse schemes for strong and weak effective approximation of the slow component. Such schemes satisfy an asymptotic preserving property and generalize the methods proposed in a recent article. We fill a gap in the analysis of these schemes and prove strong and weak error estimates, which are uniform with respect to the time scale separation parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Uniform weak error estimates for an asymptotic preserving scheme applied to a class of slow-fast parabolic semilinear SPDEs

We study an asymptotic preserving scheme for the temporal discretization of a system of parabolic semilinear SPDEs with two time scales. Owing to the averaging principle, when the time scale separation $ε$ vanishes, the slow component converges to the solution of a limiting evolution equation, which is captured when the time-step size $Δt$ vanishes by a limiting scheme. The objective of this work is to prove weak error estimates which are uniform with respect to $ε$, in terms of $Δt$: the scheme satisfies a uniform accuracy property. This is a non trivial generalization of a recent article in an infinite dimensional framework. The fast component is discretized using the modified Euler scheme for SPDEs introduced in a recent work. Proving the weak error estimates requires delicate analysis of the regularity properties of solutions of infinite dimensional Kolmogorov equations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weak error estimates of fully-discrete schemes for the stochastic Cahn-Hilliard equation

We study a class of fully-discrete schemes for the numerical approximation of solutions of stochastic Cahn--Hilliard equations with cubic nonlinearity and driven by additive noise. The spatial (resp. temporal) discretization is performed with a spectral Galerkin method (resp. a tamed exponential Euler method). We consider two situations: space-time white noise in dimension $d=1$ and trace-class noise in dimensions $d=1,2,3$. In both situations, we prove weak error estimates, where the weak order of convergence is twice the strong order of convergence with respect to the spatial and temporal discretization parameters. To prove these results, we show appropriate regularity estimates for solutions of the Kolmogorov equation associated with the stochastic Cahn--Hilliard equation, which have not been established previously and may be of interest in other contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computing return times or return periods with rare event algorithms

The average time between two occurrences of the same event, referred to as its return time (or return period), is a useful statistical concept for practical applications. For instance insurances or public agency may be interested by the return time of a 10m flood of the Seine river in Paris. However, due to their scarcity, reliably estimating return times for rare events is very difficult using either observational data or direct numerical simulations. For rare events, an estimator for return times can be built from the extrema of the observable on trajectory blocks. Here, we show that this estimator can be improved to remain accurate for return times of the order of the block size. More importantly, we show that this approach can be generalised to estimate return times from numerical algorithms specifically designed to sample rare events. So far those algorithms often compute probabilities, rather than return times. The approach we propose provides a computationally extremely efficient way to estimate numerically the return times of rare events for a dynamical system, gaining several orders of magnitude of computational costs. We illustrate the method on two kinds of observables, instantaneous and time-averaged, using two different rare event algorithms, for a simple stochastic process, the Ornstein-Uhlenbeck process. As an example of realistic applications to complex systems, we finally discuss extreme values of the drag on an object in a turbulent flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence of Adaptive Biasing Potential methods for diffusions

We prove the consistency of an adaptive importance sampling strategy based on biasing the potential energy function $V$ of a diffusion process $dX_t^0=-\nabla V(X_t^0)dt+dW_t$; for the sake of simplicity, periodic boundary conditions are assumed, so that $X_t^0$ lives on the flat $d$-dimensional torus. The goal is to sample its invariant distribution $μ=Z^{-1}\exp\bigl(-V(x)\bigr)\,dx$. The bias $V_t-V$, where $V_t$ is the new (random and time-dependent) potential function, acts only on some coordinates of the system, and is designed to flatten the corresponding empirical occupation measure of the diffusion $X$ in the large time regime. The diffusion process writes $dX_t=-\nabla V_t(X_t)dt+dW_t$, where the bias $V_t-V$ is function of the key quantity $\overlineμ_t$: a probability occupation measure which depends on the past of the process, {\it i.e.} on $(X_s)_{s\in [0,t]}$. We are thus dealing with a self-interacting diffusion. In this note, we prove that when $t$ goes to infinity, $\overlineμ_t$ almost surely converges to $μ$. Moreover, the approach is justified by the convergence of the bias to a limit which has an intepretation in terms of a free energy. The main argument is a change of variables, which formally validates the consistency of the approach. The convergence is then rigorously proven adapting the ODE method from stochastic approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-order integrator for sampling the invariant distribution of a class of parabolic SPDEs with additive space-time noise

We introduce a time-integrator to sample with high order of accuracy the invariant distribution for a class of semilinear SPDEs driven by an additive space-time noise. Combined with a postprocessor, the new method is a modification with negligible overhead of the standard linearized implicit Euler-Maruyama method. We first provide an analysis of the integrator when applied for SDEs (finite dimension), where we prove that the method has order $2$ for the approximation of the invariant distribution, instead of $1$. We then perform a stability analysis of the integrator in the semilinear SPDE context, and we prove in a linear case that a higher order of convergence is achieved. Numerical experiments, including the semilinear heat equation driven by space-time white noise, confirm the theoretical findings and illustrate the efficiency of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Weak error estimates for trajectories of SPDEs for Spectral Galerkin discretization

We consider stochastic semi-linear evolution equations which are driven by additive, spatially correlated, Wiener noise, and in particular consider problems of heat equation (analytic semigroup) and damped-driven wave equations (bounded semigroup) type. We discretize these equations by means of a spectral Galerkin projection, and we study the approximation of the probability distribution of the trajectories: test functions are regular, but depend on the values of the process on the interval $[0,T]$. We introduce a new approach in the context of quantative weak error analysis for discretization of SPDEs. The weak error is formulated using a deterministic function (Itô map) of the stochastic convolution found when the nonlinear term is dropped. The regularity properties of the Itô map are exploited, and in particular second-order Taylor expansions employed, to transfer the error from spectral approximation of the stochastic convolution into the weak error of interest. We prove that the weak rate of convergence is twice the strong rate of convergence in two situations. First, we assume that the covariance operator commutes with the generator of the semigroup: the first order term in the weak error expansion cancels out thanks to an independence property. Second, we remove the commuting assumption, and extend the previous result, thanks to the analysis of a new error term depending on a commutator.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large deviations principle for the Adaptive Multilevel Splitting Algorithm in an idealized setting

The Adaptive Multilevel Splitting (AMS) algorithm is a powerful and versatile method for the simulation of rare events. It is based on an interacting (via a mutation-selection procedure) system of replicas, and depends on two integer parameters: n $\in$ N * the size of the system and the number k $\in$ {1, . . . , n -- 1} of the replicas that are eliminated and resampled at each iteration. In an idealized setting, we analyze the performance of this algorithm in terms of a Large Deviations Principle when n goes to infinity, for the estimation of the (small) probability P(X \textgreater{} a) where a is a given threshold and X is real-valued random variable. The proof uses the technique introduced in [BLR15]: in order to study the log-Laplace transform, we rely on an auxiliary functional equation. Such Large Deviations Principle results are potentially useful to study the algorithm beyond the idealized setting, in particular to compute rare transitions probabilities for complex high-dimensional stochastic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of Adaptive Multilevel Splitting algorithms in an idealized case

The Adaptive Multilevel Splitting algorithm is a very powerful and versatile method to estimate rare events probabilities. It is an iterative procedure on an interacting particle system, where at each step, the $k$ less well-adapted particles among $n$ are killed while $k$ new better adapted particles are resampled according to a conditional law. We analyze the algorithm in the idealized setting of an exact resampling and prove that the estimator of the rare event probability is unbiased whatever $k$. We also obtain a precise asymptotic expansion for the variance of the estimator and the cost of the algorithm in the large $n$ limit, for a fixed $k$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of the Monte-Carlo error in a hybrid semi-lagrangian scheme

We consider Monte-Carlo discretizations of partial differential equations based on a combination of semi-lagrangian schemes and probabilistic representations of the solutions. We study the Monte-Carlo error in a simple case, and show that under an anti-CFL condition on the time-step $δt$ and on the mesh size $δx$ and for $N$ - the number of realizations - reasonably large, we control this error by a term of order $\mathcal{O}(\sqrt{δt /N})$. We also provide some numerical experiments to confirm the error estimate, and to expose some examples of equations which can be treated by the numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Approximation of the invariant law of SPDEs: error analysis using a Poisson equation for a full-discretization scheme

We study the long-time behavior of fully discretized semilinear SPDEs with additive space-time white noise, which admit a unique invariant probability measure $μ$. We show that the average of regular enough test functions with respect to the (possibly non unique) invariant laws of the approximations are close to the corresponding quantity for $μ$. More precisely, we analyze the rate of the convergence with respect to the different discretization parameters. Here we focus on the discretization in time thanks to a scheme of Euler type, and on a Finite Element discretization in space. The results rely on the use of a Poisson equation; we obtain that the rates of convergence for the invariant laws are given by the weak order of the discretization on finite time intervals: order $1/2$ with respect to the time-step and order $1$ with respect to the mesh-size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of a HMM time-discretization scheme for a system of Stochastic PDE's

We consider the discretization in time of a system of parabolic stochastic partial differential equations with slow and fast components; the fast equation is driven by an additive space-time white noise. The numerical method is inspired by the Averaging Principle satisfied by this system, and fits to the framework of Heterogeneous Multiscale Methods.The slow and the fast components are approximated with two coupled numerical semi-implicit Euler schemes depending on two different timestep sizes. We derive bounds of the approximation error on the slow component in the strong sense - approximation of trajectories - and in the weak sense - approximation of the laws. The estimates generalize the results of \cite{E-L-V} in the case of infinite dimensional processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approximation of the invariant measure with an Euler scheme for Stochastic PDE's driven by Space-Time White Noise

In this article, we consider a stochastic PDE of parabolic type, driven by a space-time white-noise, and its numerical discretization in time with a semi-implicit Euler scheme. When the nonlinearity is assumed to be bounded, then a dissipativity assumption is satisfied, which ensures that the SDPE admits a unique invariant probability measure, which is ergodic and strongly mixing - with exponential convergence to equilibrium. Considering test functions of class $\mathcal{C}^2$, bounded and with bounded derivatives, we prove that we can approximate this invariant measure using the numerical scheme, with order 1/2 with respect to the time step.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong and weak order in averaging for SPDEs

We show an averaging result for a system of stochastic evolution equations of parabolic type with slow and fast time scales. We derive explicit bounds for the approximation error with respect to the small parameter defining the fast time scale. We prove that the slow component of the solution of the system converges towards the solution of the averaged equation with an order of convergence is 1/2 in a strong sense - approximation of trajectories - and 1 in a weak sense - approximation of laws. These orders turn out to be the same as for the SDE case.

0 citations0 reviewsNo code linkedOpen access