Y. Efendiev,W.T. Leung,S. W. Cheung,N. Guha,V. H. Hoang,B. Mallick
Abstract:In this paper, we develop a Bayesian multiscale approach based on a multiscale finite element method. Because of scale disparity in many multiscale applications, computational models can not resolve all scales. Various subgrid models are proposed to represent un-resolved scales. Here, we consider a probabilistic approach for modeling un-resolved scales using the Multiscale Finite Element Method (cf., [1, 2]). By representing dominant modes using the Generalized Multiscale Finite Element, we propose a Bayesian framework, which provides multiple inexpensive (computable) solutions for a deterministic problem. These approximate probabilistic solutions may not be very close to the exact solutions and, thus, many realizations are needed. In this way, we obtain a rigorous probabilistic description of approximate solutions. In the paper, we consider parabolic and wave equations in heterogeneous media. In each time interval, the domain is divided into subregions. Using residual information, we design appropriate prior and posterior distributions. The likelihood consists of the residual minimization. To sample from the resulting posterior distribution, we consider several sampling strategies. The sampling involves identifying important regions and important degrees of freedom beyond permanent basis functions, which are used in residual computation. Numerical results are presented. We consider two sampling algorithms. The first algorithm uses sequential sampling and is inexpensive. In the second algorithm, we perform full sampling using the Gibbs sampling algorithm, which is more accurate compared to the sequential sampling. The main novel ingredients of our approach consist of: defining appropriate permanent basis functions and the corresponding residual; setting up a proper posterior distribution; and sampling the posteriors.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is that in multi - scale applications, due to scale differences, computational models are unable to resolve all scales. Specifically, mesoscale information in many multi - scale problems is difficult to fully capture and represent, which leads to inaccuracies in simulation results. To solve this problem, the author proposes a Bayesian multi - scale method based on the multi - scale finite element method to model unresolved sub - grid information in a probabilistic manner.
### Specific problem description
1. **Scale differences in multi - scale problems**: In many practical applications (such as porous media, materials science, etc.), problems involve multiple spatio - temporal scales. A detailed description of the finest scale often contains uncertainties, and the information at these mesoscales is usually difficult to obtain completely.
2. **Modeling of sub - grid information**: Due to the limitations of grid resolution, computational models are unable to resolve all mesoscales. To弥补 this deficiency, traditional methods use techniques such as the representative volume element (RVE) to calculate macroscopic quantities, but these methods have uncertainties in RVE size and boundary conditions, resulting in macroscopic parameters being unable to be modeled deterministically.
3. **Quantification of uncertainty**: For multi - scale problems, especially in the case of missing mesoscale information, the computational solution and its associated uncertainty are very important. Traditional deterministic methods are difficult to handle these problems, so a probability framework needs to be introduced for uncertainty quantification.
### Solutions proposed in the paper
The paper proposes a Bayesian framework based on the multi - scale finite element method for modeling missing sub - grid information. The main innovations include:
- **Defining appropriate permanent basis functions and corresponding residuals**: By selecting some dominant modes as permanent basis functions and using the generalized multi - scale finite element method (GMsFEM) to represent these basis functions.
- **Setting appropriate posterior distributions**: Construct prior and posterior distributions based on residual information, where the likelihood function is composed of residual minimization.
- **Sampling methods**: Design several sampling strategies, including sequential sampling and a full - sampling algorithm based on Gibbs sampling, to explore the posterior distribution and obtain a probabilistic description of the approximate solution.
### Mathematical expression
Consider the time - dependent equation:
\[
\frac{\partial u}{\partial t}=L(\kappa(x, t), u, \nabla u),
\]
where \(\kappa(x, t)\) is a heterogeneous spatio - temporal function and \(L\) is a differential operator. By constructing multi - scale basis functions and using a small number of permanent basis functions, combined with the Bayesian framework, probabilistic modeling of unresolved sub - grid information is carried out.
### Summary
The main objective of the paper is to develop a new Bayesian multi - scale method for probabilistically modeling missing sub - grid information in multi - scale problems and quantifying the associated uncertainty. This method not only improves computational efficiency but also provides a rigorous probabilistic description framework, which is suitable for practical problems with complex multi - scale characteristics.