312 research outputs found

    A Unified Analysis of Balancing Domain Decomposition by Constraints for Discontinuous Galerkin Discretizations

    Get PDF
    The BDDC algorithm is extended to a large class of discontinuous Galerkin (DG) discretizations of second order elliptic problems. An estimate of C(1 + log(H/h))2 is obtained for the condition number of the preconditioned system where C is a constant independent of h or H or large jumps in the coefficient of the problem. Numerical simulations are presented which confirm the theoretical results. A key component for the development and analysis of the BDDC algorithm is a novel perspective presenting the DG discretization as the sum of element-wise “local” bilinear forms. The element-wise perspective allows for a simple unified analysis of a variety of DG methods and leads naturally to the appropriate choice for the subdomain-wise local bilinear forms. Additionally, this new perspective enables a connection to be drawn between the DG discretization and a related continuous finite element discretization to simplify the analysis of the BDDC algorithm.Boeing CompanyMassachusetts Institute of Technology (Zakhartchenko Fellowship

    Balancing domain decomposition by constraints algorithms for incompressible Stokes equations with nonconforming finite element discretizations

    Get PDF
    Hybridizable Discontinuous Galerkin (HDG) is an important family of methods, which combine the advantages of both Discontinuous Galerkin in terms of flexibility and standard finite elements in terms of accuracy and efficiency. The impact of this method is partly evidenced by the prolificacy of research work in this area. Weak Galerkin (WG) is a relatively newly proposed method by introducing weak functions and generalizing the differential operator for them. This method has also drawn remarkable interests from both numerical practitioners and analysts recently. HDG and WG are different but closely related. BDDC algorithms are developed for numerical solution of elliptic problems with both methods. We prove that the optimal condition number estimate for BDDC operators with standard finite element methods can be extended to the counterparts arising from the HDG and WG methods, which are nonconforming finite element methods. Numerical experiments are conducted to verify the theoretical analysis. Further, we propose BDDC algorithms for the saddle point system arising from the Stokes equations using both HDG and WG methods. By design of the preconditioner, the iterations are restricted to a benign subspace, which makes the BDDC operator effectively positive definite thus solvable by the conjugate gradient method. We prove that the algorithm is scalable in the number of subdomains with convergence rate only dependent on subdomain problem size. The condition number bound for the BDDC preconditioned Stokes system is the same as the optimal bound for the elliptic case. Numerical results confirm the theoretical analysis
    • …
    corecore