Convergence analysis of constraint energy minimizing generalized multiscale finite element method for a linear stochastic parabolic partial differential equation driven by additive noises
Introduction
Stochastic partial differential equations (SPDEs) have recently gained much attention for their ability of fully describing intrinsic uncertainty in complex phenomena from sciences and engineering. SPDEs have been used for various applications, such as the approximative pricing of interest-rate based financial derivatives [1], nonlinear filtering [2], population biology [3], turbulence [4], and climate prediction [5]. It is rarely obtained the analytic solutions to SPDEs, so numerical methods have been developed to solve SPDEs [6], [7], [8], [9], [10], [11], [12].
Many media and materials in science and engineering have a multiple scale nature. At the same time, they are heterogeneous at a certain scale. Typical examples are human bones, composite materials, porous media, and other engineering materials. It is a great challenge to resolve all scales in the models using traditional numerical methods in a fine scale. There have been many existing approaches in the literatures to handle multiscale problems. One of these approaches is based on averaging idea and computing effective properties, such as homogenization approaches [13], numerical upscaling methods [14], [15]. Another approach is to construct multiscale basis functions on a coarse grid, such as multiscale finite element method (MsFEM) [16], [17], [18], generalized multiscale finite element method (GMsFEM) [19], [20], and constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) [21], [22] and so on. The multiscale basis functions solve some local problems, and the multiscale information is embedded in the multiscale basis functions. Thus a coarse computational model is built through a variation form and the multiscale basis functions. This can significantly reduce the computation complexity of the original fine scale model. In this paper, we will use the CEM-GMsFEM to construct multiscale basis functions. The CEM-GMsFEM shares the framework of the GMsFEM. To construct multiscale basis functions, a constraint energy minimization is imposed on auxiliary basis functions via local spectral problems. Because the basis functions have a fast decay outside the target coarse block, we can compute the basis functions in a local region.
Stochastic partial differential equations (SPDEs) with multiscale structures can be used to model accurately multiscale dynamical systems affected by noises. The purpose of this paper is to present a numerical multiscale method and solve SPDEs efficiently. We consider linear stochastic parabolic partial differential equations driven by additive noises with multiscale diffusion coefficients, and carry out the convergence analysis for the CEM-GMsFEM. The main idea of the CEM-GMsFEM is to use the multiscale basis functions which can be accurately computed to capture complex multiscale features. The challenge in the convergence analysis of multiscale finite element approximation for the SPDEs is the lack of regularity of white noises. For instance, as shown in [23], the required regularity conditions of the SPDEs solution are not satisfied for the standard error estimates of finite element methods. To overcome this difficulty, Allen, Novosel, and Zhang [23] have given a piecewise constant approximation of the white noises, but the regularity of white noises are not improved much due to the nature of standard Wiener process. In order to obtain “a more regular” noises and improve the regularity of the solution, Du and Zhang [24] represent noises term as a linear combination of Fourier modes. Following this idea to approximate white noises, the approximation solutions of SPDEs converge to the exact solutions as the white noise approximation becomes finer. We have proved that the CEM-GMsFEM achieves the first order convergence with respect to coarse grid size in -norm under some mild assumption on the SPDEs. The stability analysis and convergence analysis are also presented for fully discretizing spatial space, temporal space and white noises.
This paper is organized as follows: In Section 2, we present the white noise approximation for a stochastic parabolic equation. The error estimate between the exact solution and approximate solution is established. In Section 3, the CEM-GMsFEM is briefly introduced. In Section 4, we present a convergence analysis for the SPDEs using the CEM-GMsFEM and discretized white noises. In Section 5, the full discrete scheme is presented. In Section 6, a few numerical results are presented to confirm the error analysis. Finally, some conclusions are given.
Section snippets
Approximation of stochastic parabolic differential equations
In this paper, we will consider the following linear stochastic parabolic differential equation with form where is a given deterministic function, denotes a linear, self-adjoint, positive definite and not necessarily bounded operator, such that is compactly embedded into a Hilbert space . The solution is an -valued random process. Here denotes an infinite dimensional Brownian motion on a probability space .
Construction of CEM-GMsFEM basis functions
In this section, we follow [21] and present the construction of CEM-GMsFEM basis functions. There are two stages for it. In the first stage, we construct the auxiliary space by solving a local spectral problem on each coarse grid . Based on the auxiliary space, the second stage is to construct the multiscale basis functions by solving some local constraint energy minimization problems.
Next, we give the notions of the fine grid, coarse grid and oversampling domain for the computational domain
The semidiscrete formulation
In this section, we will present the convergence analysis of CEM-GMsFEM for a linear stochastic parabolic equation in the sense of -norm.
Fully discrete scheme
In this section, we utilize the backward Euler method to discretize the temporal variable for the stochastic parabolic problem. For convenience, for any , we define
Numerical results
In this section, we present a few representative numerical examples to show the performance of CEM-GMsFEM for solving linear stochastic parabolic equations driven by additive noises. For the time discretization, we will use the backward Euler scheme. In Section 6.1, we present some numerical examples with high-contrast diffusion field to demonstrate the convergence of CEM-GMsFEM. In Section 6.2, the effects of the stochastic parabolic equation driven by different noises are considered. For
Conclusions
In this paper, we have presented CEM-GMsFEM for solving a linear stochastic parabolic equation driven by additive noises with multiscale diffusion coefficients. The analysis showed that the convergence depends on the coarse grid size, eigenvalue decay of local spectral problems and the stochastic noise approximation, but is independent of the scale length and contrast of the media. For the additive noise, it was represented in a general basis to facilitate convergence analysis for the
Acknowledgments
The authors gratefully appreciate the valuable comments from the reviewers, which have contributed significantly to the improvement of this manuscript. This research was supported by the National Natural Science Foundation of China (11871378, 11301392), Shanghai Peak Discipline Program for Higher Education Institutions (Class I)-Civil Engineering, China, and Fundamental Research Funds for Central Universities, China (No. 22120180529).
References (27)
- et al.
Upscaling method for problems in perforated domains with non-homogeneous boundary conditions on perforations using non-local multi-continuum method (NLMC)
J. Comput. Appl. Math.
(2019) - et al.
A multiscale finite element method for elliptic problems in composite materials and porous media
J. Comput. Phys.
(1997) - et al.
Generalized multiscale finite element methods
J. Comput. Phys.
(2013) - et al.
Constraint energy minimizing generalized multiscale finite element method, oversampling strategies
Comput. Methods Appl. Mech. Engrg.
(2018) - et al.
Consistent recalibration of yield curve models
Math. Finance
(2018) On the optimal filtering of diffusion processes
Probab. Theory Related Fields
(1969)- et al.
Measure-valued processes and renormalization of branching particle systems
Functionals and the random-force method in turbulence theory
Sov. J. Exp. Theor. Phys.
(1965)- et al.
Models for stochastic climate prediction
Appl. Math.
(1996) - et al.
A numerical approximation of parabolic stochastic partial differential equations driven by a Poisson random measure
BIT
(2006)
Numerical methods for stochastic parabolic PDEs
Numer. Funct. Anal. Optim.
Numerical Approximations of Stochastic Partial Differential Equations
The numerical approximation of stochastic partial differential equations
Milan J. Math.
Cited by (3)
Stochastic heat equation: Numerical positivity and almost surely exponential stability
2022, Computers and Mathematics with ApplicationsCitation Excerpt :Utilizing the properties of truncated random variables and M-matrix, the unboundedness is overcame and the numerical positivity is obtained. The stability analysis of the stochastic heat equation is always an active research area (see [1,7–9]). Our main aim is to investigate the numerical almost surely exponential stability, which is an important conception for the stochastic system.
ONLINE MULTISCALE MODEL REDUCTION FOR NONLINEAR STOCHASTIC PDES WITH MULTIPLICATIVE NOISE
2022, Multiscale Modeling and Simulation