BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Charles-Edouard Bréhier

Charles-Edouard Bréhier contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.PRNumerical Analysiscond-mat.stat-mechphysics.comp-phphysics.data-an

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 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