Abstract:We study the effect of adaptive mesh refinement on a parallel domain decomposition solver of a linear system of algebraic equations. These concepts need to be combined within a parallel adaptive finite element software. A prototype implementation is presented for this purpose. It uses adaptive mesh refinement with one level of hanging nodes. Two and three-level versions of the Balancing Domain Decomposition based on Constraints (BDDC) method are used to solve the arising system of algebraic equations. The basic concepts are recalled and components necessary for the combination are studied in detail. Of particular interest is the effect of disconnected subdomains, a typical output of the employed mesh partitioning based on space-filling curves, on the convergence and solution time of the BDDC method. It is demonstrated using a large set of experiments that while both refined meshes and disconnected subdomains have a negative effect on the convergence of BDDC, the number of iterations remains acceptable. In addition, scalability of the three-level BDDC solver remains good on up to a few thousands of processor cores. The largest presented problem using adaptive mesh refinement has over 10^9 unknowns and is solved on 2048 cores.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to combine Adaptive Mesh Refinement (AMR) with non - overlapping domain decomposition solvers in the challenges of large - scale parallel computing. Specifically, the authors focus on how to effectively integrate these two techniques into an efficient parallel finite - element software and evaluate the impact of AMR on the Balancing Domain Decomposition based on Constraints (BDDC) solver. ### Core of the Problem 1. **Adaptive Mesh Refinement (AMR)**: - AMR is a technique for locally improving computational accuracy, especially near complex regions or boundary layers. - In large - scale parallel computing, the implementation of AMR is very challenging because it needs to handle load balancing, hanging node management, and the selection of appropriate algebraic solvers. 2. **Non - overlapping domain decomposition solver (BDDC)**: - The BDDC method is an effective solver for linear equations, especially suitable for large - scale parallel computing. - The authors studied the influence of AMR on the convergence and solution time of the BDDC method, especially when the sub - domains are not connected. ### Main Objectives - **Develop and implement large - scale parallel AMR code**: This code is coupled with a non - overlapping domain decomposition solver to evaluate the impact of AMR on the solver's performance. - **Study the handling of hanging nodes**: Prove that the handling of hanging nodes can be relatively simple and straightforward, rather than the common belief that non - conforming meshes should be avoided. - **Evaluate the impact of disconnected sub - domains**: Explore the influence of disconnected sub - domains caused by space - filling curve partitioning on the convergence and solution time of the BDDC method. ### Key Technologies and Methods - **Space - filling curve**: Used to partition the computational grid to ensure that each sub - domain contains approximately an equal number of elements. - **Multilevel BDDC method**: Extends the BDDC method to enable it to better handle complex problems in large - scale parallel computing. - **Hanging node handling**: Only allow first - level hanging nodes and use equal - order elements to simplify the algorithm implementation. ### Experimental Results Through a large number of experiments, the authors showed that although the refined grid and disconnected sub - domains have a negative impact on the convergence of the BDDC method, the number of iterations still remains within an acceptable range. In addition, the three - level BDDC solver shows good scalability on thousands of processor cores, the largest test problem contains more than \(10^9\) unknowns, and was successfully solved on 2048 cores. ### Summary This paper provides an effective large - scale parallel computing framework by combining the AMR and BDDC methods, which can handle complex grid structures and large - scale data sets while ensuring computational accuracy.

Coupling parallel adaptive mesh refinement with a nonoverlapping domain decomposition solver

Domain decomposition approach for parallel improvement of tetrahedral meshes.

Parallel Algorithms for Moving Boundary Problems by Local Remeshing

Fine-Grained Parallel Algorithm for Unstructured Surface Mesh Generation.

Scalable Parallel Quadrilateral Mesh Generation Coupled with Mesh Partitioning

Fine-grained parallel algorithm for unstructured surface mesh generation

A Parallel Local Reconnection Approach for Tetrahedral Mesh Improvement

Scalable generation of large-scale unstructured meshes by a novel domain decomposition approach.

Implementation of balancing domain decomposition method for parallel finite element analysis involving inactive elements

Multispace and multilevel BDDC

Balancing Domain Decomposition Method on Additive Schwartz Framework for Multi-Level Implementation

Parallel domain discretization algorithm for RBF-FD and other meshless numerical methods for solving PDEs

Generalized bridging domain method for coupling finite elements with discrete elements

Newton-Krylov-BDDC deluxe solvers for non-symmetric fully implicit time discretizations of the Bidomain model

Parallel CFD-Algorithm on Unstructured Adaptive Meshes

A Highly Solid Model Boundary Preserving Method for Large-Scale Parallel 3D Delaunay Meshing on Parallel Computers

Domain Decomposition Based High Performance Parallel Computing

An explicit substructuring method for overlapping domain decomposition based on stochastic calculus

Iteration-merged Solver for Linear Banded Matrices Based on Domain Decomposition and MPI

Parallel Scheme for Multi-Layer Refinement Non-Uniform Grid Lattice Boltzmann Method Based on Load Balancing

A parallel dual-grid multiscale approach to CFD-DEM couplings