Abstract:Partial differential equations are discretized using the mortar finite element method, where the mortar space contains piecewise quadratic and -cubic functions. We prove the wellposedness of the resulting saddle point system under a theoretical framework by Braess et al. [SIAM J. Numer. Analy. (1999)]. For solving this saddle point system, we introduce two domain decomposition (DD) algorithms. Under the assumption that every local problem approximates the global problem well in the interior of its own partition, the required communication complexity is proportional to the size of the interface. With a globally conforming fine mesh, a DD algorithm by Bank et al. [SIAM Rev. (2003), SIAM J. Sci. Comput. (2004)] is used, in which both the interface solution and the interface residual of each local problem are sent to the neighboring processors. To accommodate globally nonconforming fine meshes, we propose a new algorithm where only the interface solution is communicated. Consequently, the communication complexity is reduced by 1/2. The convergence property of these DD algorithms are studied on one-dimensional two-grid meshes. If the interior residual of each local problem is also communicated, the DD algorithms are shown to converge in one iteration. Otherwise, these algorithms are shown to converge in two iterations. Numerical results are presented, on solving the Poisson equation and a convection diffusion equation with a highly adaptive mesh underlying the finite element discretization. The scalability and parallel efficiency of our DD algorithms are studied. As the number of processors increases, linear and logarithmic speed-up of the solving time is observed for the convection-diffusion equation and the Poisson equation respectively. Furthermore, several parallel algebraic multilevel algorithms are introduced. Three two-level and one three-level algorithms are proposed and implemented, together with an overlapping Schwartz method, which uses a three-level algorithm followed by a two-level one at each iteration. The scalability results indicate that these algorithms are effective in reducing the communication complexity during the solving stage. In summary, our algorithms improve the coordination among processors during the solving stage and may provide a powerful tool for scientific computing.

A finite-difference method for the variable coefficient Poisson equation on hierarchical Cartesian meshes

A Fast Direct Solver for Elliptic PDEs on a Hierarchy of Adaptively Refined Quadtrees

High-order finite-difference ghost-point methods for elliptic problems in domains with curved boundaries

High order finite-difference ghost-point methods for elliptic problems in domains with curved boundaries

An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

A numerical strategy to discretize and solve the Poisson equation on dynamically adapted multiresolution grids for time-dependent streamer discharge simulations

High Order Finite Difference Discretization for Composite Grid Hierarchy and Its Applications

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

A metric-based adaptive mesh refinement criterion under constrain for solving elliptic problems on quad/octree grids

A Fast Algorithm for Solving the Poisson Equation on a Nested Grid

A Non-Nested Multilevel Method for Meshless Solution of the Poisson Equation in Heat Transfer and Fluid Flow

Scalable parallel multilevel algorithms for solving partial differential equations

Neural-Network-Assisted Finite Difference Discretization for Numerical Solution of Partial Differential Equations

High-order Finite Element--Integral Equation Coupling on Embedded Meshes

Adaptive Finite Element Methods for Elliptic Equations over Hierarchical T-meshes.

A New Multilevel Method for Electrostatic Problems Through Hierarchical Loop Basis

A Moving Mesh Finite Element Algorithm for Singular Problems in Two and Three Space Dimensions

A sparse hierarchical $hp$-finite element method on disks and annuli

A parallel shared-memory implementation of a high-order accurate solution technique for variable coefficient Helmholtz problems

A locally second order symmetric method for discontinuous solution of Poisson's equation on uniform cartesian grids