Michael Presho

Mehdi Ghommem, Michael Presho, Victor M. Calo et al.

In this paper, we combine concepts of the generalized multiscale finite element method and mode decomposition methods to construct a robust local-global approach for model reduction of flows in high-contrast porous media. This is achieved by implementing proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) techniques on a coarse grid. The resulting reduced-order approach enables a significant reduction in the flow problem size while accurately capturing the behavior of fully resolved solutions. We consider a variety of high-contrast coefficients and present the corresponding numerical results to illustrate the effectiveness of the proposed technique. This paper is a continuation of the first part where we examine the applicability of POD and DMD to derive simplified and reliable representations of flows in high-contrast porous media. In the current paper, we discuss how these global model reduction approaches can be combined with local techniques to speed-up the simulations. The speed-up is due to inexpensive, while sufficiently accurate, computations of global snapshots.

Lijian Jiang, Michael Presho

In this paper we use a splitting technique to develop new multiscale basis functions for the multiscale finite element method (MsFEM). The multiscale basis functions are iteratively generated using a Green's kernel. The Green's kernel is based on the first differential operator of the splitting. The proposed MsFEM is applied to deterministic elliptic equations and stochastic elliptic equations, and we show that the proposed MsFEM can considerably reduce the dimension of the random parameter space for stochastic problems. By combining the method with sparse grid collocation methods, the need for a prohibitive number of deterministic solves is alleviated. We rigorously analyze the convergence of the proposed method for both the deterministic and stochastic elliptic equations. Computational complexity discussions are also offered to supplement the convergence analysis. A number of numerical results are presented to confirm the theoretical findings.

Michael Presho, Juan Galvis

In this paper, we propose a method for the construction of locally conservative flux fields through a variation of the Generalized Multiscale Finite Element Method (GMsFEM). The flux values are obtained through the use of a Ritz formulation in which we augment the resulting linear system of the continuous Galerkin (CG) formulation in the higher-order GMsFEM approximation space. In particular, we impose the finite volume-based restrictions through incorporating a scalar Lagrange multiplier for each mass conservation constraint. This approach can be equivalently viewed as a constraint minimization problem where we minimize the energy functional of the equation restricted to the subspace of functions that satisfy the desired conservation properties. To test the performance of the method we consider equations with heterogeneous permeability coefficients that have high-variation and discontinuities, and couple the resulting fluxes to a two-phase flow model. The increase in accuracy associated with the computation of the GMsFEM pressure solutions is inherited by the flux fields and saturation solutions, and is closely correlated to the size of the reduced-order systems. In particular, the addition of more basis functions to the enriched multiscale space produces solutions that more accurately capture the behavior of the fine scale model. A variety of numerical examples are offered to validate the performance of the method.