Abstract:In this paper, we propose offline and online adaptive enrichment algorithms for the generalized multiscale approximation of a mixed finite element method with velocity elimination to solve the subsurface flow problem in high-contrast and heterogeneous porous media. In the offline adaptive method, we first derive an a-posteriori error indicator based on one weighted L2-norm of the local residual operator, where the weighted L2-norm is related to the pressure fields of the local snapshot space. Then we enrich the multiscale space by increasing the number of offline basis functions iteratively on coarse elements where the error indicator takes large values. While in the online adaptive method, we add online basis functions on selected coarse elements based on another weighted L2-norm of the local residual operator to enrich the multiscale space, here the weighted L2-norm is associated with the velocity fields of the local snapshot space. Online basis functions are constructed in the online stage depending on the solution of the previous iteration and some optimal estimates. We give theoretical analyses for the convergences of these two adaptive methods, which show that sufficient initial basis functions (belong to the offline space) lead to faster convergence rates. A series of numerical examples are provided to highlight the performances of both these two adaptive methods and also validate the theoretical analyses. Both offline and online adaptive methods are effective that can reduce the relative error substantially. In addition, the online adaptive method generally performs better than the offline adaptive method as online basis functions contain important global information such as distant effects that cannot be captured by offline basis functions. The numerical results also show that with a suitable initial multiscale space that includes all offline basis functions corresponding to relative smaller eigenvalues of each local spectral decomposition in the offline stage, the convergence rate of the online enrichment is independent of the permeability contrast.

Parallel Algorithm Combined with Mixed Element Procedure for Compressible Miscible Displacement Problem

Parallel Algorithms for Moving Boundary Problems by Local Remeshing

A parallel finite element variational multiscale method based on fully overlapping domain decomposition for incompressible flows

Hybrid mixed discontinuous Galerkin finite element method for incompressible miscible displacement problem

An upwind-mixed volume element method on changing meshes for compressible miscible displacement problem

Parallel Reservoir Simulators for Fully Implicit Complementarity Formulation of Multicomponent Compressible Flows.

Convergence analysis of a nonconforming virtual element method for compressible miscible displacement problems in porous media

Parallel Domain Decomposition method for the fully-mixed Stokes-dual-permeability fluid flow model with Beavers-Joseph interface conditions

Multi-Point Flux MFE Decoupled Method for Compressible Miscible Displacement Problem

Nonconforming virtual element method for an incompressible miscible displacement problem in porous media

Parallel finite element method for 2D free surface seepage

Parallel Computation of Compressible Flows for Implicit Schemes Without Sub-iteration and Modification of the Implicit Solver

An efficient parallel algorithm with application to computational fluid dynamics

Adaptive Space-Time Domain Decomposition for Multiphase Flow in Porous Media with Bound Constraints

Adaptive generalized multiscale approximation of a mixed finite element method with velocity elimination

Multiscale Mixed Methods with Improved Accuracy: The Role of Oversampling and Smoothing

Heterogeneous parallel computing method for 3D transient nonlinear thermomechanical problems on CPU-GPU platforms

Parallel in time partially explicit splitting scheme for high contrast multiscale problems

Parallel Computing Of Cavity Flow By Finite Element Method

Parallel Domain Decomposition Methods with Mixed Order Discretization for Fully Implicit Solution of Tracer Transport Problems on the Cubed-Sphere.