Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method

TL;DR

This paper introduces a novel multiscale discontinuous Galerkin method that efficiently reduces computational costs in simulating flow and wave propagation in heterogeneous media by extracting dominant modes through energy minimization.

Contribution

It develops an energy minimization based multiscale model reduction technique with discontinuous Galerkin discretization, enhancing efficiency and accuracy in high-contrast heterogeneous media simulations.

Findings

Convergence is linear in coarse mesh size and contrast-independent.

Method effectively captures high-contrast features with minimal degrees of freedom.

Numerical results demonstrate improved performance in flow and wave simulations.

Abstract

Numerical simulation of flow problems and wave propagation in heterogeneous media has important applications in many engineering areas. However, numerical solutions on the fine grid are often prohibitively expensive, and multiscale model reduction techniques are introduced to efficiently solve for an accurate approximation on the coarse grid. In this paper, we propose an energy minimization based multiscale model reduction approach in the discontinuous Galerkin discretization setting. The main idea of the method is to extract the non-decaying component in the high conductivity regions by identifying dominant modes with small eigenvalues of local spectral problems, and define multiscale basis functions in coarse oversampled regions by constraint energy minimization problems. The multiscale basis functions are in general discontinuous on the coarse grid and coupled by interior penalty discontinuous Galerkin formulation. The minimal degree of freedom in representing high-contrast features is achieved through the design of local spectral problems, which provides the most compressed local multiscale space. We analyze the method for solving Darcy flow problem and show that the convergence is linear in coarse mesh size and independent of the contrast, provided that the oversampling size is appropriately chosen. Numerical results are presented to show the performance of the method for simulation on flow problem and wave propagation in high-contrast heterogeneous media.