The eigenvalues and eigenspaces of some discrete div- and curl-related operators are investigated. The discrete operators give some good discrete analogues of the continuous counterparts and play an important role in developing finite volume schemes for solving div-curl equations and electromagnetic systems. Knowledge of the eigenvalues and eigenspaces is very useful in the numerical analysis of finite volume methods for electromagetic systems in <i>nonhomogeneous</i> media.

In this paper, we consider the numerical approximation of the three-dimensional poroelastic wave equations in the spherical coordinate system. One difficulty in the design of an efficient numerical scheme is that the problem is singular in the center and the polar axes of the computational domain. Nevertheless, we develop a hybrid finite difference/control volume method for solving this problem. Our method is explicit and is second order accurate in both space and time. Numerical results are shown to confirm the convergence rate of our method and the effectiveness to simulate wave propagation in poroelastic media in the spherical coordinate system.

We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling. This poroelasticity problem suffers from rapidly oscillating material parameters, which calls for a thorough numerical treatment. In this paper, we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity. Therein, local corrector problems are constructed in line with the static equations, whereas we propose to consider the full system. This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure. We prove the optimal first-order convergence of this method and verify the result by numerical experiments.

In this paper, we develop an adaptive generalized multiscale discontinuous Galerkin method (GMsDGM) for a class of high-contrast flow problems and derive a priori and a posteriori error estimates for the method. Based on the a posteriori error estimator, we develop an adaptive enrichment algorithm for our GMsDGM and prove its convergence. The adaptive enrichment algorithm gives an automatic way to enrich the approximation space in regions where the solution requires more basis functions, which are shown to perform well compared with a uniform enrichment. We also discuss an approach that adaptively selects multiscale basis functions by correlating the residual to multiscale basis functions (cf. [S. S. Chen, D. L. Donoho, and M. A. Saunders, <i>SIAM Rev.</i>, 43 (2001), pp. 129--159]). The proposed error indicators are L_2-based and can be inexpensively computed, which makes our approach efficient

In this paper, a class of FETI-DP preconditioners is developed for a fast solution of the linear system arising from staggered discontinuous Galerkin discretization of the two-dimensional Stokes equations. The discretization has been recently developed and has the distinctive advantages that it is optimally convergent and has a good local conservation property. In order to efficiently solve the linear system, two kinds of FETI-DP preconditioners, namely, lumped and Dirichlet preconditioners, are considered and analyzed. Scalable bounds C (H/h) and C (1+ log (H/h)) 2 are proved for the lumped and Dirichlet preconditioners, respectively, with the constant C depending on the infsup constant of the discrete spaces but independent of any mesh parameters. Here H/h stands for the number of elements across each subdomain. Numerical results are presented to confirm the theoretical estimates.