Wei Leng

Wei Leng, Lili Ju

In this paper, we propose and analyze an additive domain decomposition method (DDM) for solving the high-frequency Helmholtz equation with the Sommerfeld radiation condition. In the proposed method, the computational domain is partitioned into structured subdomains along all spatial directions, and each subdomain contains an overlapping region for source transferring. At each iteration all subdomain PML problems are solved completely in parallel, then all horizontal, vertical and corner directional residuals on each subdomain are passed to its corresponding neighbor subdomains as the source for the next iteration. This DDM method is highly scalable in nature and theoretically shown to produce the exact solution for the PML problem defined in ${\mathbb{R}}^2$ in the constant medium case. A slightly modified version of the method for bounded truncated domains is also developed for its use in practice and an error estimate is rigorously proved. Various numerical experiments in two and three dimensions are conducted on the supercomputer "Tianhe-2 Cluster" to verify the theoretical results and demonstrate excellent performance of the proposed method as an iterative solver or a preconditioner.

Thi-Thao-Phuong Hoang, Lili Ju, Wei Leng et al.

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RKDG) methods and design LTS algorithms based on strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integrations in different regions. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

Hui Liu, Tao Cui, Wei Leng et al.

In this paper a new hp-adaptive strategy for elliptic problems based on refinement history is proposed, which chooses h-, p- or hp-refinement on individual elements according to a posteriori error estimate, as well as smoothness estimate of the solution obtained by comparing the actual and expected error reduction rate. Numerical experiments show that exponential convergence can be achieved with this strategy.

Young-Ju Lee, Wei Leng, Chen-Song Zhang

In this paper, we revisit an auxiliary space preconditioning method proposed by Xu [Computing 56, 1996], in which low-order finite element spaces are employed as auxiliary spaces for solving linear algebraic systems arising from high-order finite element discretizations. We provide a new convergence rate estimate and parallel implementation of the proposed algorithm. We show that this method is user-friendly and can play an important role in a variety of Poisson-based solvers for more challenging problems such as the Navier--Stokes equation. We investigate the performance of the proposed algorithm using the Poisson equation and the Stokes equation on 3D unstructured grids. Numerical results demonstrate the advantages of the proposed algorithm in terms of efficiency, robustness, and parallel scalability.

Wei Leng, Lili Ju

In this paper, based on the overlapping domain decomposition method (DDM) proposed in \cite{Leng2015}, an one step preconditioner is proposed to solve 2D high frequency Helmholtz equation. The computation domain is decomposed in both $x$ and $y$ directions, and the local solution on each subdomain is updated simultaneously in one iteration, thus there is no sweeping along certain directions. In these ways, the overlapping DDM is similar to the popular DDM for Poisson problem. The one step preconditioner simply take the restricted source on each subdomain, solve the local problems and summarize the local solutions on all subdomains including their PML area. The complexity of solving the problem with the preconditioner is $O(N n_{\text{iter}})$, where $n_{\text{iter}}$ is the number of iteration, and it is shown numerically that $n_{\text{iter}}$ is proportional to the number of subdomains in one direction. 2D Helmholtz problem with nearly a billion unknowns are solved efficiently with the preconditioner on massively parallel machines.

Wei Leng

A fast method is proposed for solving the high frequency Helmholtz equation. The building block of the new fast method is an overlapping source transfer domain decomposition method for layered medium, which is an extension of the source transfer domain decomposition method proposed by Chen and Xiang \cite{Chen2013a,Chen2013b}. The new fast method contains a setup phase and a solving phase. In the setup phase, the computation domain is decomposed hierarchically into many subdomains of different levels, and the mapping from incident traces to field traces on all the subdomains are set up bottom-up. In the solving phase, first on the bottom level, the local problem on the subdomains with restricted source is solved, then the wave propagates on the boundaries of all the subdomains bottom-up, at last the local solutions on all the subdomains are summed up top-down. The total computation cost of the new fast method is $O(n^{\frac{3}{2}} \log n)$ for 2D problem. Numerical experiments shows that with the new fast method, Helmholtz equations with half billion unknowns could be solved efficiently on massively parallel machines.