TY - JOUR
T1 - Study of micro-macro acceleration schemes for linear slow-fast stochastic differential equations with additive noise
AU - Debrabant, Kristian
AU - Samaey, Giovanni
AU - Zieliński, Przemysław
PY - 2020
Y1 - 2020
N2 - Computational multi-scale methods capitalize on a large time-scale separation to efficiently simulate slow dynamics over long time intervals. For stochastic systems, one often aims at resolving the statistics of the slowest dynamics. This paper looks at the efficiency of a micro-macro acceleration method that couples short bursts of stochastic path simulation with extrapolation of spatial averages forward in time. To have explicit derivations, we elicit an amenable linear test equation containing multiple time scales. We make derivations and perform numerical experiments in the Gaussian setting, where only the evolution of mean and variance matters. The analysis shows that, for this test model, the stability threshold on the extrapolation step is largely independent of the time-scale separation. In consequence, the micro-macro acceleration method increases the admissible time steps far beyond those for which a direct time discretization becomes unstable.
AB - Computational multi-scale methods capitalize on a large time-scale separation to efficiently simulate slow dynamics over long time intervals. For stochastic systems, one often aims at resolving the statistics of the slowest dynamics. This paper looks at the efficiency of a micro-macro acceleration method that couples short bursts of stochastic path simulation with extrapolation of spatial averages forward in time. To have explicit derivations, we elicit an amenable linear test equation containing multiple time scales. We make derivations and perform numerical experiments in the Gaussian setting, where only the evolution of mean and variance matters. The analysis shows that, for this test model, the stability threshold on the extrapolation step is largely independent of the time-scale separation. In consequence, the micro-macro acceleration method increases the admissible time steps far beyond those for which a direct time discretization becomes unstable.
KW - math.NA
KW - 65C30, 60H35, 65L20
KW - Kullback–Leibler divergence
KW - Stability
KW - Micro–macro simulations
KW - Stiff stochastic differential equations
KW - Entropy optimisation
U2 - 10.1007/s10543-020-00804-5
DO - 10.1007/s10543-020-00804-5
M3 - Journal article
SN - 0006-3835
VL - 60
SP - 959
EP - 998
JO - BIT Numerical Mathematics
JF - BIT Numerical Mathematics
ER -