A BDDC (Balancing Domain Decomposition by Constraints) algorithm is developed and analyzed for a staggered discontinuous Galerkin (DG) finite element approximation of second order scalar elliptic problems. On a quite irregular subdomain partition, an optimal condition number bound is proved for two-dimensional problems. In addition, a sub-optimal but scalable condition number bound is obtained for three-dimensional problems. These bounds are shown to be independent of coefficient jumps in the subdomain partition. Numerical results are also included to show the performance of the algorithm.

In this paper, we develop an adaptive Generalized Multiscale Discontinuous Galerkin Method (GMs-DGM) 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.[4]). The proposed error indicators are L2-based and can be inexpensively computed which makes our approach efficient. Numerical results are presented that demonstrate the robustness of the proposed error indicators.

Convergence theories and a deluxe dual and primal finite element tearing and interconnecting algorithm are developed for a hybrid staggered DG finite element approximation of H(curl) elliptic problems in two dimensions. In addition to the advantages of staggered DG methods, the basis functions of the new hybrid staggered DG method are all locally supported in the triangular elements, and a Lagrange multiplier approach is applied to enforce the global connections of these basis functions. The interface problem on the Lagrange multipliers is further reduced to the resulting problem on the subdomain interfaces, and a dual and primal finite element tearing and interconnecting algorithm with an enriched weight factor is then applied to the resulting problem. Our algorithm is shown to give a condition number bound of <i>C</i>(1+log(<i>H</i><i>h</i>))², independent of the two parameters, where <i>H</i><i>h</i> is the number of

In this paper, we present an efficient computational methodology for diffusion and convection-diffusion problems in highly heterogeneous media as well as convection-dominated diffusion problem. It is well known that the numerical computation for these problems requires a significant amount of computer memory and time. Nevertheless, the solutions to these problems typically contain a coarse component, which is usually the quantity of interest and can be represented with a small number of degrees of freedom. There are many methods that aim at the computation of the coarse component without resolving the full details of the solution. Our proposed method falls into the framework of interior penalty discontinuous Galerkin method, which is proved to be an effective and accurate class of methods for numerical solutions of partial differential equations. A distinctive feature of our method is that the solution space

Discontinuous Galerkin (DG) methods are a class of efficient tools for solving fluid flow problems. There are in the literature many greatly successful DG methods. In this paper, a new staggered DG method for the Stokes system is developed and analyzed. The key feature of our method is that the discrete system preserves the structures of the continuous problem, which results from the use of our new staggered DG spaces. This also provides local and global conservation properties, which are desirable for fluid flow applications. The method is based on the first order mixed formulation involving pressure, velocity, and velocity gradient. The velocity and velocity gradient are approximated by polynomials of the same degree while the choice of polynomial degree for pressure is flexible, namely, the approximation degree for pressure can be chosen as either that of velocity or one degree lower than that of velocity. In any