Adaptive reduced basis finite element heterogeneous multiscale method

An adaptive reduced basis finite element heterogeneous multiscale method (RB-FE-HMM) is proposed for elliptic problems with multiple scales. The multiscale method is based on the RB-FE-HMM introduced in [A. Abdulle, Y. Bai, Reduced basis finite element heterogeneous multiscale method for high-order discretizations of elliptic homogenization problems, J. Comput. Phys. 231 (21) (2012) 7014-7036]. It couples a macroscopic solver with effective data recovered from the solution of micro problems solved on sampling domains. Unlike classical numerical homogenization methods, the micro problems are computed in a finite dimensional space spanned by a small number of accurately computed representative micro solutions (the reduced basis) obtained by a greedy algorithm in an offline stage. In this paper we present a residual-based a posteriori error analysis in the energy norm as well as an a posteriori error analysis in quantities of interest. For both type of adaptive strategies, rigorous a posteriori error estimates are derived and corresponding error estimators are proposed. In contrast to the adaptive finite element heterogeneous multiscale method (FE-HMM), there is no need to adapt the micro mesh simultaneously to the macroscopic mesh refinement. Up to an offline preliminary stage, the RB-FE-HMM has the same computational complexity as a standard adaptive FEM for the effective problem. Two and three dimensional numerical experiments confirm the efficiency of the RB-FE-HMM and illustrate the improvements compared to the adaptive FE-HMM. (C) 2013 Elsevier B.V. All rights reserved

Reduced basis finite element heterogeneous multiscale method for high-order discretizations of elliptic homogenization problems

A new finite element method for the efficient discretization of elliptic homogenization problems is proposed. These problems, characterized by data varying over a wide range of scales cannot be easily solved by classical numerical methods that need mesh resolution down to the finest scales and multiscale methods capable of capturing the large scale components of the solution on macroscopic meshes are needed. Recently, the finite element heterogeneous multiscale method (FE-HMM) has been proposed for such problems, based on a macroscopic solver with effective data recovered from the solution of micro problems on sampling domains at quadrature points of a macroscopic mesh. Departing from the approach used in the FE-HMM, we show that interpolation techniques based on the reduced basis methodology (an offline-online strategy) allow one to design an efficient numerical method relying only on a small number of accurately computed micro solutions. This new method, called the reduced basis finite element heterogeneous multiscale method (RB-FE-HMM) is significantly more efficient than the FE-HMM for high order macroscopic discretizations and for three-dimensional problems, when the repeated computation of micro problems over the whole computational domain is expensive. A priori error estimates of the RB-FE-HMM are derived. Numerical computations for two and three dimensional problems illustrate the applicability and efficiency of the numerical method. (C) 2012 Elsevier Inc. All rights reserved

Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods

In this paper, we discuss a general multiscale model reduction framework based on multiscale finite element methods. We give a brief overview of related multiscale methods. Due to page limitations, the overview focuses on a few related methods and is not intended to be comprehensive. We present a general adaptive multiscale model reduction framework, the Generalized Multiscale Finite Element Method. Besides the method's basic outline, we discuss some important ingredients needed for the method's success. We also discuss several applications. The proposed method allows performing local model reduction in the presence of high contrast and no scale separation

An isogeometric analysis for elliptic homogenization problems

A novel and efficient approach which is based on the framework of isogeometric analysis for elliptic homogenization problems is proposed. These problems possess highly oscillating coefficients leading to extremely high computational expenses while using traditional finite element methods. The isogeometric analysis heterogeneous multiscale method (IGA-HMM) investigated in this paper is regarded as an alternative approach to the standard Finite Element Heterogeneous Multiscale Method (FE-HMM) which is currently an effective framework to solve these problems. The method utilizes non-uniform rational B-splines (NURBS) in both macro and micro levels instead of standard Lagrange basis. Beside the ability to describe exactly the geometry, it tremendously facilitates high-order macroscopic/microscopic discretizations thanks to the flexibility of refinement and degree elevation with an arbitrary continuity level provided by NURBS basis functions. A priori error estimates of the discretization error coming from macro and micro meshes and optimal micro refinement strategies for macro/micro NURBS basis functions of arbitrary orders are derived. Numerical results show the excellent performance of the proposed method

Variational Multiscale Stabilization and the Exponential Decay of Fine-scale Correctors

This paper addresses the variational multiscale stabilization of standard finite element methods for linear partial differential equations that exhibit multiscale features. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to local cell problems is rigorously justified. The stabilization eliminates scale-dependent pre-asymptotic effects as they appear for standard finite element discretizations of highly oscillatory problems, e.g., the poor $L^2$ approximation in homogenization problems or the pollution effect in high-frequency acoustic scattering

Nonlinear nonlocal multicontinua upscaling framework and its applications

In this paper, we discuss multiscale methods for nonlinear problems. The main idea of these approaches is to use local constraints and solve problems in oversampled regions for constructing macroscopic equations. These techniques are intended for problems without scale separation and high contrast, which often occur in applications. For linear problems, the local solutions with constraints are used as basis functions. This technique is called Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM). GMsFEM identifies macroscopic quantities based on rigorous analysis. In corresponding upscaling methods, the multiscale basis functions are selected such that the degrees of freedom have physical meanings, such as averages of the solution on each continuum. This paper extends the linear concepts to nonlinear problems, where the local problems are nonlinear. The main concept consists of: (1) identifying macroscopic quantities; (2) constructing appropriate oversampled local problems with coarse-grid constraints; (3) formulating macroscopic equations. We consider two types of approaches. In the first approach, the solutions of local problems are used as basis functions (in a linear fashion) to solve nonlinear problems. This approach is simple to implement; however, it lacks the nonlinear interpolation, which we present in our second approach. In this approach, the local solutions are used as a nonlinear forward map from local averages (constraints) of the solution in oversampling region. This local fine-grid solution is further used to formulate the coarse-grid problem. Both approaches are discussed on several examples and applied to single-phase and two-phase flow problems, which are challenging because of convection-dominated nature of the concentration equation