Abstract:We analyze the random fluctuations of several multiscale algorithms, such as the multiscale finite element method (MsFEM) and the finite element heterogeneous multiscale method (HMM), that have been developed to solve partial differential equations with highly heterogeneous coefficients. Such multiscale algorithms are often shown to correctly capture the homogenization limit when the highly oscillatory random medium is stationary and ergodic. This paper is concerned with the random fluctuations of the solution about the deterministic homogenization limit. We consider the simplified setting of the one-dimensional elliptic equation, where the theory of random fluctuations is well understood. We develop a fluctuation theory for the multiscale algorithms in the presence of random environments with short-range and long-range correlations. For a given mesh size h, we show that the fluctuations converge in distribution in the space of continuous paths to Gaussian processes as the correlation length epsilon -> 0. We next derive the limit of such Gaussian processes as h -> 0 and compare this limit with the distribution of the random fluctuations of the continuous model. When such limits agree, we conclude that the multiscale algorithm captures the random fluctuations accurately and passes the corrector test. This property serves as an interesting benchmark to assess the behavior of the multiscale algorithm in practical situations where the assumptions necessary for the theory of homogenization are not met. What we find is that the computationally more expensive methods MsFEM, and HMM with a choice of parameter delta = h, correctly capture the random fluctuations both for short-range and long-range oscillations in the medium. The less expensive method HMM with delta < h correctly captures the fluctuations for long-range oscillations and strongly amplifies their size in media with short-range oscillations. We present a modified scheme with an intermediate computational cost that captures the random fluctuations in all cases.

Corrector theory for elliptic equations with long-range correlated random potential.

Corrector Theory for Elliptic Equations in Random Media with Singular Green’s Function. Application to Random Boundaries

Corrector Analysis of a Heterogeneous Multi-scale Scheme for Elliptic Equations with Random Potential

Homogenization and Corrector Theory for Linear Transport in Random Media

Limiting distribution of elliptic homogenization error with periodic diffusion and random potential

Scaling Limits for One-Dimensional Long-Range Percolation: Using the Corrector Method

Artificial boundary conditions for random ellitpic systems with correlated coefficient field

Corrector Theory for MsFEM and HMM in Random Media

Optimal artificial boundary conditions based on second-order correctors for three dimensional random elliptic media

Fluctuations in the Homogenization of Semilinear Equations with Random Potentials

Localization of Random Surfaces with Monotone Potentials and an FKG-Gaussian Correlation Inequality

Corrections to scaling in random resistor networks and diluted continuous spin models near the percolation threshold

Weak Mean Random Attractor of Reversible Selkov Lattice Systems Driven By Locally Lipschitz Lévy Noises

Discrete random stabilities of attractors for nonlocal lattice equations with implicit or colored noise

Approximation diffusion for the Nonlinear Schrödinger equation with a random potential

Stability estimates of inverse random source problems for the wave equations by using correlation-based data

Calderon-Zygmund estimates for stochastic elliptic systems on bounded Lipschitz domains

Limit Theorem for a Rough Differential Equation with a Negative Long-Range Random Coefficient

Random numerical stability of attractors for nonlinear Schrodinger equations on infinite lattices

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods