Source author record

Gang Bao

Gang Bao 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

18works
7topics
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

18 published item(s)

preprint2026arXiv

Harnessing AI for Inverse Partial Differential Equation Problems: Past, Present, and Prospects

Solving inverse partial differential equation (PDE) problems is a fundamental topic in scientific research due to its broad significance across a wide range of real-world applications. Inverse PDE problems arise across medical imaging, geophysics, materials science, and aerodynamics, where the goal is to infer hidden causes, design structures, or control physical states. In this paper, we provide a comprehensive review of recent advances in solving inverse PDE problems using artificial intelligence (AI). We first introduce the basic formulation, key challenges, and traditional numerical foundations of inverse PDE problems, and then organize it into three major categories: inverse problems, inverse design, and control problems. For each category, we further present a methodological paradigms, and review representative state-of-the-art approaches from recent years. We then summarize representative applications across scientific and industrial domains, including mechanical systems, aerodynamic problems, thermal systems, full-waveform inversion, system identification, and medical imaging. Finally, we discuss open challenges and future prospects, such as physics-informed architectures, limited real-world data, uncertainty quantification, and inverse foundation models. This survey aims to provide the first unified and systematic perspective on AI for inverse PDE problems, demonstrating how modern learning-based methods are reshaping inverse problems, inverse design, and control problems in PDE-governed systems.

preprint2022arXiv

An Adaptive Finite Element DtN Method for Maxwell's Equations

This paper is concerned with a numerical solution to the scattering of a time-harmonic electromagnetic wave by a bounded and impenetrable obstacle in three dimensions. The electromagnetic wave propagation is modeled by a boundary value problem of Maxwell's equations in the exterior domain of the obstacle. Based on the Dirichlet-to-Neumann (DtN) operator, which is defined by an infinite series, an exact transparent boundary condition is introduced and the scattering problem is reduced equivalently into a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is developed to solve the discrete variational problem, where the DtN operator is truncated into a sum of finitely many terms. The a posteriori error estimate takes into account both the finite element approximation error and the truncation error of the DtN operator. The latter is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented to illustrate the effectiveness of the proposed method.

preprint2021arXiv

An adaptive finite element DtN method for the elastic wave scattering problem in three dimensions

Consider the elastic scattering of an incident wave by a rigid obstacle in three dimensions, which is formulated as an exterior problem for the Navier equation. By constructing a Dirichlet-to-Neumann (DtN) operator and introducing a transparent boundary condition, the scattering problem is reduced equivalently to a boundary value problem in a bounded domain. The discrete problem with the truncated DtN operator is solved by using the a posteriori error estimate based adaptive finite element method. The estimate takes account of both the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to converge exponentially with respect to the truncation parameter. Moreover, the generalized Woodbury matrix identity is utilized to solve the resulting linear system efficiently. Numerical experiments are presented to demonstrate the superior performance of the proposed method.

preprint2020arXiv

An adaptive finite element DtN method for the elastic wave scattering by biperiodic structures

Consider the scattering of a time-harmonic elastic plane wave by a bi-periodic rigid surface. The displacement of elastic wave motion is modeled by the three-dimensional Navier equation in an open domain above the surface. Based on the Dirichlet-to-Neumann (DtN) operator, which is given as an infinite series, an exact transparent boundary condition is introduced and the scattering problem is formulated equivalently into a boundary value problem in a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is proposed to solve the discrete variational problem where the DtN operator is truncated into a finite number of terms. The a posteriori error estimate takes account of the finite element approximation error and the truncation error of the DtN operator which is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented to illustrate the effectiveness of the proposed method.

preprint2020arXiv

An adaptive finite element DtN method for the open cavity scattering problems

Consider the scattering of a time-harmonic electromagnetic plane wave by an open cavity which is embedded in a perfectly electrically conducting infinite ground plane. This paper concerns the numerical solutions of the open cavity scattering problems in both transverse magnetic and transverse electric polarizations. Based on the Dirichlet-to-Neumann (DtN) map for each polarization, a transparent boundary condition is imposed to reduce the scattering problem equivalently into a boundary value problem in a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is proposed. The estimate consists of the finite element approximation error and the truncation error of the DtN operator, which is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented for both polarizations to illustrate the competitive behavior of the adaptive method.

preprint2020arXiv

An adaptive finite element DtN method for the three-dimensional acoustic scattering problem

This paper is concerned with a numerical solution of the acoustic scattering by a bounded impenetrable obstacle in three dimensions. The obstacle scattering problem is formulated as a boundary value problem in a bounded domain by using a Dirichlet-to-Neumann (DtN) operator. An a posteriori error estimate is derived for the finite element method with the truncated DtN operator. The a posteriori error estimate consists of the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to decay exponentially with respect to the truncation parameter. Based on the a posteriori error estimate, an adaptive finite element method is developed for the obstacle scattering problem. The truncation parameter is determined by the truncation error of the DtN operator and the mesh elements for local refinement are marked through the finite element approximation error. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.

preprint2020arXiv

An inverse spectral problem for a damped wave operator

This paper proposes a new and efficient numerical algorithm for recovering the damping coefficient from the spectrum of a damped wave operator, which is a classical Borg-Levinson inverse spectral problem. The algorithm is based on inverting a sequence of trace formulas, which are deduced by a recursive formula, bridging geometrical and spectrum information explicitly in terms of Fredholm integral equations. Numerical examples are presented to illustrate the efficiency of the proposed algorithm.

preprint2020arXiv

Inverse scattering by a random periodic structure

This paper develops an efficient numerical method for the inverse scattering problem of a time-harmonic plane wave incident on a perfectly reflecting random periodic structure. The method is based on a novel combination of the Monte Carlo technique for sampling the probability space, a continuation method with respect to the wavenumber, and the Karhunen-Lo$\grave{e}$ve expansion of the random structure, which reconstructs key statistical properties of the profile for the unknown random periodic structure from boundary measurements of the scattered fields away from the structure. Numerical results are presented to demonstrate the reliability and efficiency of the proposed method.

preprint2020arXiv

Recovering simultaneously a potential and a point source from Cauchy data

This paper is devoted to the inverse problem of recovering simultaneously a potential and a point source in a Shrödinger equation from the associated nonlinear Dirichlet to Neumann map. The uniqueness of the inversion is proved and logarithmic stability estimates are derived. It is well known that the inverse problem of determining only the potential while knowing the source, is ill-posed. In contrast the problem of identifying a point source when the potential is given is well posed. The obtained results show that the nonlinear Dirichlet to Neumann map contains enough information to determine simultaneously the potential and the point source. However recovering a point source imbedded in an unknown background medium becomes an ill-posed inversion.

preprint2020arXiv

Stability for the multifrequency inverse medium problem

The solution of a multi-frequency 1d inverse medium problem consists of recovering the refractive index of a medium from measurements of the scattered waves for multiple frequencies. In this paper, rigorous stability estimates are derived when the frequency takes value in a bounded interval.It is showed that the ill-posedness of the inverse medium problem decreases as the width of the frequency interval becomes larger. More precisely, under certain regularity assumptions on the refractive index, the estimates indicate that the power in Hölder stability is an increasing function of the largest value in the frequency band. Finally, a Lipschitz stability estimate is obtained for the observable part of the medium function defined through a truncated trace formula.

preprint2020arXiv

Weak Adversarial Networks for High-dimensional Partial Differential Equations

Solving general high-dimensional partial differential equations (PDE) is a long-standing challenge in numerical mathematics. In this paper, we propose a novel approach to solve high-dimensional linear and nonlinear PDEs defined on arbitrary domains by leveraging their weak formulations. We convert the problem of finding the weak solution of PDEs into an operator norm minimization problem induced from the weak formulation. The weak solution and the test function in the weak formulation are then parameterized as the primal and adversarial networks respectively, which are alternately updated to approximate the optimal network parameter setting. Our approach, termed as the weak adversarial network (WAN), is fast, stable, and completely mesh-free, which is particularly suitable for high-dimensional PDEs defined on irregular domains where the classical numerical methods based on finite differences and finite elements suffer the issues of slow computation, instability and the curse of dimensionality. We apply our method to a variety of test problems with high-dimensional PDEs to demonstrate its promising performance.

preprint2016arXiv

Direct and Inverse Elastic Scattering From Anisotropic Media

Assume a time-harmonic elastic wave is incident onto a penetrable anisotropic body embedded into a homogeneous isotropic background medium. We propose an equivalent variational formulation in a truncated bounded domain and show the uniqueness and existence of weak solutions by applying the Fredholm alternative and using properties of the Dirichlet-to-Neumann map in both two and three dimensions. The Fréchet derivative of the near-field solution operator with respect to the boundary of the scatterer is derived. As an application, we design a descent algorithm for recovering the interface from the near-field data of one or several incident directions and frequencies. Numerical examples in 2D are presented to show the validity and accuracy of our methods.

preprint2016arXiv

Inverse Random Source Scattering for Elastic Waves

This paper is concerned with the direct and inverse random source scattering problems for elastic waves where the source is assumed to be driven by an additive white noise. Given the source, the direct problem is to determine the displacement of the random wave field. The inverse problem is to reconstruct the mean and variance of the random source from the boundary measurement of the wave field at multiple frequencies. The direct problem is shown to have a unique mild solution by using a constructive proof. Based on the explicit mild solution, Fredholm integral equations of the first kind are deduced for the inverse problem. The regularized Kaczmarz method is presented to solve the ill-posed integral equations. Numerical experiments are included to demonstrate the effectiveness of the proposed method.

preprint2014arXiv

Sensitivity Analysis of an Inverse Problem for the Wave Equation with Caustics

The paper investigates the sensitivity of the inverse problem of recovering the velocity field in a bounded domain from the boundary dynamic Dirichlet-to-Neumann map (DDtN) for the wave equation. Three main results are obtained: (1) assuming that two velocity fields are non-trapping and are equal to a constant near the boundary, it is shown that the two induced scattering relations must be identical if their corresponding DDtN maps are sufficiently close; (2) a geodesic X-ray transform operator with matrix-valued weight is introduced by linearizing the operator which associates each velocity field with its induced Hamiltonian flow. A selected set of geodesics whose conormal bundle can cover the cotangent space at an interior point is used to recover the singularity of the X-ray transformed function at the point; a local stability estimate is established for this case. Although fold caustics are allowed along these geodesics, it is required that these caustics contribute to a smoother term in the transform than the point itself. The existence of such a set of geodesics is guaranteed under some natural assumptions in dimension greater than or equal to three by the classification result on caustics and regularity theory of Fourier Integral Operators. The interior point with the above required set of geodesics is called "fold-regular"; (3) assuming that a background velocity field with every interior point fold-regular is fixed and another velocity field is sufficiently close to it and satisfies a certain orthogonality condition, it is shown that if the two corresponding DDtN maps are sufficiently close then they must be equal.

preprint2013arXiv

Nearly Cloaking the Full Maxwell Equations

The approximate cloaking is investigated for time-harmonic Maxwell's equations via the approach of transformation optics. The problem is reduced to certain boundary effect estimates due to an inhomogeneous electromagnetic inclusion with an asymptotically small support but an arbitrary content enclosed by a thin high-conducting layer. Sharp estimates are established in terms of the asymptotic parameter, which are independent of the material tensors of the small electromagnetic inclusion. The result implies that the `blow-up-a-small-region' construction via the transformation optics approach yields a near-cloak for the electromagnetic waves. A novelty lies in the fact that the geometry of the cloaking construction of this work can be very general. Moreover, by incorporating the conducting layer developed in the present paper right between the cloaked region and the cloaking region, arbitrary electromagnetic contents can be nearly cloaked. Our mathematical technique extends the general one developed in [30] for nearly cloaking scalar optics. In order to investigate the approximate electromagnetic cloaking for general geometries with arbitrary cloaked contents, new techniques and analysis tools must be developed for this more challenging vector optics case.

preprint2013arXiv

Nearly Cloaking the Full Maxwell Equations: Cloaking Active Contents with General Conducting Layers

The regularized near-cloak via the transformation optics approach in the time-harmonic electromagnetic scattering is considered. This work extends the existing studies mainly in two aspects. First, it presents a near-cloak construction by incorporating a much more general conducting layer between the cloaked and cloaking regions. This might be of significant practical importance when production fluctuations occur. Second, it allows the cloaked contents to be both passive and active with an applied current inside. In assessing the near-cloaking performance, comprehensive and sharp estimates are derived for the scattering amplitude in terms of the asymptotic regularization parameter and the material tensors of the conducting layer. The scattering estimates are independent of the passive/active contents being cloaked, which implies that one could nearly cloak arbitrary contents by using the proposed near-cloak construction.

preprint2007arXiv

Improved Simulation of the Mass Charging for ASTROD I

The electrostatic charging of the test mass in ASTROD I (Astrodynamical Space Test of Relativity using Optical Devices I) mission can affect the quality of the science data as a result of spurious Coulomb and Lorentz forces. To estimate the size of the resultant disturbances, credible predictions of charging rates and the charging noise are required. Using the GEANT4 software toolkit, we present a detailed Monte Carlo simulation of the ASTROD I test mass charging due to exposure of the spacecraft to galactic cosmic-ray (GCR) protons and alpha particles (3He, 4He) in the space environment. A positive charging rate of 33.3 e+/s at solar minimum is obtained. This figure reduces by 50% at solar maximum. Based on this charging rate and factoring in the contribution of minor cosmic-ray components, we calculate the acceleration noise and stiffness associated with charging. We conclude that the acceleration noise arising from Coulomb and Lorentz effects are well below the ASTROD I acceleration noise limit at 0.1 mHz both at solar minimum and maximum. The coherent Fourier components due to charging are investigated, it needs to be studied carefully in order to ensure that these do not compromise the quality of science data in the ASTROD I mission.

preprint2000arXiv

Models for Accretion Disk Fluctuations through Self-Organized Criticality Including Relativistic Effects

The possibility that some of the observed X-ray and optical variability in active galactic nuclei and galactic black hole candidates is produced in accretion disks through the development of a self-organized critical state is reconsidered. New simulations, including more complete calculations of relativistic effects, do show that this model can produce light-curves and power-spectra for the variability which agree with the range observed in optical and X-ray studies of AGN and X-ray binaries. However, the universality of complete self-organized criticality is not quite achieved. This is mainly because the character of the variations depend quite substantially on the extent of the unstable disk region. If it extends close to the innermost stable orbit then a physical scale is introduced and the scale-free character of self-organized criticality is vitiated. Significant dependence of the power spectrum density slope on the type of diffusion within the disk and a weaker dependence on the amount of differential rotation are noted. When general relativistic effects are incorporated in the models, additional substantial differences are produced if the disk is viewed from directions far from the accretion disk axis.