Abstract:In this paper, we consider the numerical approximation of a time-fractional stochastic Cahn--Hilliard equation driven by an additive fractionally integrated Gaussian noise. The model involves a Caputo fractional derivative in time of order $\alpha\in(0,1)$ and a fractional time-integral noise of order $\gamma\in[0,1]$. The numerical scheme approximates the model by a piecewise linear finite element method in space and a convolution quadrature in time (for both time-fractional operators), along with the $L^2$-projection for the noise. We carefully investigate the spatially semidiscrete and fully discrete schemes, and obtain strong convergence rates by using clever energy arguments. The temporal Hölder continuity property of the solution played a key role in the error analysis. Unlike the stochastic Allen--Cahn equation, the presence of the unbounded elliptic operator in front of the cubic nonlinearity in the underlying model adds complexity and challenges to the error analysis. To overcome these difficulties, several new techniques and error estimates are developed. The study concludes with numerical examples that validate the theoretical findings.
What problem does this paper attempt to address?
The paper aims to address the numerical approximation problem of the time-fractional stochastic Cahn-Hilliard equation. Specifically:
1. **Research Background**: The Cahn-Hilliard equation is typically used to describe phase separation and coarsening phenomena in alloys, as well as spinodal decomposition in binary mixtures. The introduction of noise is to better simulate random disturbances or unknown parameters in actual physical processes.
2. **Research Object**: The paper considers a time-fractional stochastic Cahn-Hilliard equation driven by a fractional time derivative (Caputo derivative), and the equation includes a fractional integral noise.
3. **Numerical Method**: The paper proposes a spatial discretization scheme based on the finite element method, combined with the convolution quadrature method to handle the time-fractional derivative. Additionally, an $L^2$ projection is used to handle the noise term.
4. **Main Contribution**: Through a clever energy estimation method, the authors prove the strong convergence of the numerical solution and obtain the convergence rate. In particular, the paper overcomes the complexity brought by the presence of nonlinear terms and the elliptic operator.
5. **Theoretical Results**: Through detailed error analysis, the time Hölder continuity property is obtained, and based on this, the convergence rate of the numerical solution is derived.
6. **Experimental Verification**: Finally, numerical experiments are conducted to verify the effectiveness of the theoretical results.
In summary, the core issue of this paper is to develop an effective numerical method to approximate the time-fractional stochastic Cahn-Hilliard equation and to prove its theoretical convergence.