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.

Parallel Fully Coupled Methods for Bound-Preserving Solution of Subsurface Flow and Transport in Porous Media

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

Scalable Semismooth Newton Methods with Multilevel Domain Decomposition for Subsurface Flow and Reactive Transport in Porous Media

Parallel Algorithms for Moving Boundary Problems by Local Remeshing

Nonlinear Parallel-in-time Simulations of Multiphase Flow in Porous Media.

Parallel multilevel restricted Schwarz preconditioners for implicit simulation of subsurface flows with Peng-Robinson equation of state

Adaptive Space-Time Domain Decomposition for Multiphase Flow in Porous Media with Bound Constraints

Parallel Reservoir Simulators for Fully Implicit Complementarity Formulation of Multicomponent Compressible Flows.

A Scalable Fully Implicit Framework for Reservoir Simulation on Parallel Computers

Adaptively Coupled Domain Decomposition Method for Multiphase and Multicomponent Porous Media Flows

Active-Set Reduced-Space Methods with Nonlinear Elimination for Two-Phase Flow Problems in Porous Media.

A Schur Complement Preconditioner for Scalable Parallel Fluid Simulation

A parallel solver framework for fully implicit monolithic fluid-structure interaction

An efficient parallel algorithm with application to computational fluid dynamics

Parallel Multilevel Methods for Implicit Solution of Shallow Water Equations with Nonsmooth Topography on the Cubed-Sphere

Multilevel Field-Split Preconditioners with Domain Decomposition for Steady and Unsteady Flow Problems

Scalable parallel multilevel algorithms for solving partial differential equations

Parallel Primal-Dual Active-Set Algorithm with Nonlinear and Linear Preconditioners

Parallel Domain Decomposition method for the fully-mixed Stokes-dual-permeability fluid flow model with Beavers-Joseph interface conditions

A Parallel Unstructured Finite-Volume Method for All-Speed Flows

On Robust and Efficient Parallel Reservoir Simulation on Tianhe 2