Analysis of a Two-Scale, Locally Conservative Subgrid Upscaling for Elliptic Problems

We present a two-scale theoretical framework for approximating the solution of a second order elliptic problem. The elliptic coefficient is assumed to vary on a scale that can be resolved on a fine numerical grid, but limits on computational power require that computations be performed on a coarse grid. We consider the elliptic problem in mixed variational form over W × V ⊂ L² × H( div). We base our scale expansion on local mass conservation over the coarse grid. It is used to define a direct sum decomposition of W × V into coarse and "subgrid" subspaces Wc × Vc and δ W × δ V such that (1) ∇ · Vc=Wc and ∇ · δ V = δ W, and (2) the space δ V is locally supported over the coarse mesh. We then explicitly decompose the variational problem into coarse and subgrid scale problems. The subgrid problem gives a well-defined operator taking Wc × Vc to δ W × δ V, which is localized in space, and it is used to upscale, that is, to remove the subgrid from the coarse-scale problem. Using standard mixed finite element spaces, two-scale mixed spaces are defined. A mixed approximation is defined, which can be viewed as a type of variational multiscale method or a residual-free bubble technique. A numerical Green's function approach is used to make the approximation to the subgrid operator efficient to compute. A mixed method π-operator is defined for the two-scale approximation spaces and used to show optimal order error estimates.