A convergent stochastic scalar auxiliary variable method
arXiv:2308.07060 · doi:10.1093/imanum/drae065
Abstract
We discuss an extension of the scalar auxiliary variable approach, which was originally introduced by Shen et al. ([Shen, Xu, Yang, J. Comput. Phys., 2018]) for the discretization of deterministic gradient flows. By introducing an additional scalar auxiliary variable, this approach allows to derive a linear scheme, while still maintaining unconditional stability. Our extension augments the approximation of the evolution of this scalar auxiliary variable with higher order terms, which enables its application to stochastic partial differential equations. Using the stochastic Allen--Cahn equation as a prototype for nonlinear stochastic partial differential equations with multiplicative noise, we propose an unconditionally energy stable, linear, fully discrete finite element scheme based on our augmented scalar auxiliary variable method. Recovering a discrete version of the energy estimate and establishing Nikolskii estimates with respect to time, we are able to prove convergence of discrete solutions towards pathwise unique martingale solutions by applying Jakubowski's generalization of Skorokhod's theorem. A generalization of the Gyöngy--Krylov characterization of convergence in probability to quasi-Polish spaces finally provides convergence of fully discrete solutions towards strong solutions of the stochastic Allen--Cahn equation. Finally, we present numerical simulations underlining the practicality of the scheme and the importance of the introduced augmentation terms.
References in corpus (7)
- A generalized SAV approach with relaxation for dissipative systems
- Random numbers for large scale distributed Monte Carlo simulations
- Rate of Convergence of Space Time Approximations for stochastic evolution equations
- Optimal Strong Rates of Convergence for a Space-Time Discretization of the Stochastic Allen-Cahn Equation with multiplicative noise
- A convergent SAV scheme for Cahn--Hilliard equations with dynamic boundary conditions
- Semi-implicit energy-preserving numerical schemes for stochastic wave equation via SAV approach
- Stability and convergence of relaxed scalar auxiliary variable schemes for Cahn-Hilliard systems with bounded mass source