Source author record

Xuejun Xu

Xuejun Xu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2026arXiv

Neural Green's Function Accelerated Iterative Methods for Solving Indefinite Boundary Value Problems

Neural operators, which learn mappings between the function spaces, have been applied to solve boundary value problems in various ways, including learning mappings from the space of the forcing terms to the space of the solutions with the substantial requirements of data pairs. In this work, we present a data-free neural operator integrated with physics, which learns the Green kernel directly. Our method proceeds in three steps: 1. The governing equations for the Green's function are reformulated into an interface problem, where the delta Dirac function is removed; 2. The interface problem is embedded in a lifted space of higher-dimension to handle the jump in the derivative, but still solved on a two-dimensional surface without additional sampling cost; 3. Deep neural networks are employed to address the curse of dimensionality caused by this lifting operation. The approximate Green's function obtained through our approach is then used to construct preconditioners for the linear systems allowed by its mathematical properties. Furthermore, the spectral bias of it revealed through both theoretical analysis and numerical validation contrasts with the smoothing effects of traditional iterative solvers, which motivates us to propose a hybrid iterative method that combines these two solvers. Numerical experiments demonstrate the effectiveness of our approximate Green's function in accelerating iterative methods, proving fast convergence for solving indefinite problems even involving discontinuous coefficients.

preprint2021arXiv

A Two-Level Preconditioned Helmholtz-Jacobi-Davidson Method for the Maxwell Eigenvalue Problem

In this paper, based on a domain decomposition (DD) method, we shall propose an efficient two-level preconditioned Helmholtz-Jacobi-Davidson (PHJD) method for solving the algebraic eigenvalue problem resulting from the edge element approximation of the Maxwell eigenvalue problem. In order to eliminate the components in orthogonal complement space of the eigenvalue, we shall solve a parallel preconditioned system and a Helmholtz projection system together in fine space. After one coarse space correction in each iteration and minimizing the Rayleigh quotient in a small dimensional Davidson space, we finally get the error reduction of this two-level PHJD method as $γ=c(H)(1-C\frac{δ^{2}}{H^{2}})$, where $C$ is a constant independent of the mesh size $h$ and the diameter of subdomains $H$, $δ$ is the overlapping size among the subdomains, and $c(H)$ decreasing as $H\to 0$, which means the greater the number of subdomains, the better the convergence rate. Numerical results supporting our theory shall be given.

preprint2020arXiv

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

In this paper, we develop an adaptive finite element method for the nonlinear steady-state Poisson-Nernst-Planck equations, where the spatial adaptivity for geometrical singularities and boundary layer effects are mainly considered. As a key contribution, the steady-state Poisson-Nernst-Planck equations are studied systematically and rigorous analysis for a residual-based a posteriori error estimate of the nonlinear system is presented. With the help of Schauder fixed point theorem, we show the solution existence and uniqueness of the linearized system derived by taking $G-$derivatives of the nonlinear system, followed by the proof of the relationship between the error of solution and the a posteriori error estimator $η$. Numerical experiments are given to validate the efficiency of the a posteriori error estimator and demonstrate the expected rate of convergence. In the further tests, adaptive mesh refinements for geometrical singularities and boundary layer effects are successfully observed.

preprint2015arXiv

Domain Decomposition Methods for Space Fractional Partial Differential Equations

In this paper, a two-level additive Schwarz preconditioner is proposed for solving the algebraic systems resulting from the finite element approximations of space fractional partial differential equations (SFPDEs). It is shown that the condition number of the preconditioned system is bounded by C(1+H/δ), where H is the maximum diameter of subdomains and δis the overlap size among the subdomains. Numerical results are given to support our theoretical findings.

preprint2014arXiv

A FETI-DP preconditioner of discontinuous Galerkin method for multiscale problems in high constrast media

In this paper we consider second order elliptic partial differential equations with highly varying (heterogeneous) coefficients on a two-dimensional region. The problems are discretized by a composite finite element (FE) and discontinuous Galerkin (DG) Method. The fine grids are in general nonmatching across the subdomain boundaries, and the subdomain partitioning does not need to resolve the jumps in the coefficient. A FETI-DP preconditioner is proposed and analyzed to solve the resulting linear system. Numerical results are presented to support our theory.

preprint2013arXiv

A Robust Multilevel Method for Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation

A robust multilevel preconditioner based on the hybridizable discontinuous Galerkin method for the Helmholtz equation with high wave number is presented in this paper. There are two keys in our algorithm, one is how to choose a suitable intergrid transfer operator, and the other is using GMRES smoothing on coarse grids. The multilevel method is performed as a preconditioner in the outer GMRES iteration. To give a quantitative insight of our algorithm, we use local Fourier analysis to analyze the convergence property of the proposed multilevel method. Numerical results show that for fixed wave number, the convergence of the algorithm is mesh independent. Moreover, the performance of the algorithm depends relatively mildly on wave number.

preprint2013arXiv

Multilevel Preconditioner with Stable Coarse Grid Corrections for the Helmholtz Equation

In this paper we consider a class of robust multilevel precontioners for the Helmholtz equation with high wave number. The key idea in this work is to use the continuous interior penalty finite element methods (CIP-FEM) studied in \cite{Wu12,Wu12-hp} to construct the stable coarse grid correction problems. The multilevel methods, based on GMRES smoothing on coarse grids, are then served as a preconditioner in the outer GMRES iteration. In the one dimensional case, convergence property of the modified multilevel methods is analyzed by the local Fourier analysis. From our numerical results, we find that the proposed methods are efficient for a reasonable range of frequencies. The performance of the algorithms depends relatively mildly on wave number. In particular, only one GMRES smoothing step may guarantee the optimal convergence of our multilevel algorithm, which remedies the shortcoming of the multilevel algorithm in \cite{EEO01}.

preprint2012arXiv

A Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation with High Wave Number

This paper analyzes the error estimates of the hybridizable discontinuous Galerkin (HDG) method for the Helmholtz equation with high wave number in two and three dimensions. The approximation piecewise polynomial spaces we deal with are of order $p\geq 1$. Through choosing a specific parameter and using the duality argument, it is proved that the HDG method is stable without any mesh constraint for any wave number $κ$. By exploiting the stability estimates, the dependence of convergence of the HDG method on $κ,h$ and $p$ is obtained. Numerical experiments are given to verify the theoretical results.

preprint2012arXiv

Domain Decomposition Preconditioners for a Discontinuous Galerkin Formulation of a Multiscale Elliptic Problem

In this paper, we propose a domain decomposition method for multiscale second order elliptic partial differential equations with highly varying coefficients. The method is based on a discontinuous Galerkin formulation. We present both a nonoverlapping and an overlapping version of the method. We prove that the condition number bound of the preconditioned algebraic system in either case can be made independent of the coefficients under certain assumptions. Also, in our analysis, we do not need to assume that the coefficients are continuous across the coarse grid boundaries. The analysis and the condition number bounds are new, and contribute towards further extension of the theory for the discontinuous Galerkin discretization for multiscale problems.