In this paper, we discuss the application of generalized multiscale finite element method (GMsFEM) to elasticity equation in heterogeneous media. We consider steady state elasticity equations though some of our applications are motivated by elastic wave propagation in subsurface where the subsurface properties can be highly heterogeneous and have high contrast. We present the construction of main ingredients for GMsFEM such as the snapshot space and offline spaces. The latter is constructed using local spectral decomposition in the snapshot space. The spectral decomposition is based on the analysis which is provided in the paper. We consider both continuous Galerkin and discontinuous Galerkin coupling of basis functions. Both approaches have their cons and pros. Continuous Galerkin methods allow avoiding penalty parameters though they involve partition of unity functions which can alter

In this paper, a new type of staggered discontinuous Galerkin methods for the three dimensional Maxwells equations is developed and analyzed. The spatial discretization is based on staggered Cartesian grids so that many good properties are obtained. First of all, our method has the advantages that the numerical solution preserves the electromagnetic energy and automatically fulfills a discrete version of the Gauss law. Moreover, the mass matrices are diagonal, thus time marching is explicit and is very efficient. Our method is high order accurate and the optimal order of convergence is rigorously proved. It is also very easy to implement due to its Cartesian structure and can be regarded as a generalization of the classical Yees scheme as well as the quadrilateral edge finite elements. Furthermore, a superconvergence result, that is the convergence rate is one order higher at interpolation nodes, is proved

Numerical modeling of wave propagation in heterogeneous media is important in many applications. Due to their complex nature, direct numerical simulations on the fine grid are prohibitively expensive. It is therefore important to develop efficient and accurate methods that allow the use of coarse grids. In this paper, we present a multiscale finite element method for wave propagation on a coarse grid. The proposed method is based on the generalized multiscale finite element method (GMsFEM) (see [Y. Efendiev, J. Galvis, and T. Hou, <i>J. Comput. Phys.</i>, 251 (2012), pp. 116--135]). To construct multiscale basis functions, we start with two snapshot spaces in each coarse-grid block, where one represents the degrees of freedom on the boundary and the other represents the degrees of freedom in the interior. We use local spectral problems to identify important modes in each snapshot space. These local spectral

In this paper, we derive an a-posteriori error indicator for the Generalized Multiscale Finite Element Method (GMsFEM) framework. This error indicator is further used to develop an adaptive enrichment algorithm for the linear elliptic equation with multiscale high-contrast coefficients. The GMsFEM, which has recently been introduced in [13], allows solving multiscale parameter-dependent problems at a reduced computational cost by constructing a reduced-order representation of the solution on a coarse grid. The main idea of the method consists of (1) the construction of snapshot space,(2) the construction of the offline space, and (3) the construction of the online space (the latter for parameter-dependent problems). In [13], it was shown that the GMsFEM provides a flexible tool to solve multiscale problems with a complex input space by generating appropriate snapshot, offline, and online spaces. In this paper, we

We introduce an irrational factor of order <i>k</i> defined by I k ( n ) = i = 1 l p i i , where I k ( n ) = i = 1 l p i i is the factorization of <i>n</i> and $${\beta_{i} = \left\{\begin{array}{ll}\alpha_i, \quad \quad {\rm if} \quad \alpha_i &lt; k \\ \frac{1}{\alpha_i},\quad \quad {\rm if} \quad \alpha_i \geqq k \end{array}\right.}$$. It turns out that the function I k ( n ) = i = 1 l p i i well approximates the characteristic function of <i>k</i>-free integers. We also derive asymptotic formulas for I k ( n ) = i = 1 l p i i and I k ( n ) = i = 1 l p i i .