Source author record

Hengguang Li

Hengguang Li 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

8works
5topics
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

8 published item(s)

preprint2022arXiv

A $C^0$ finite element method for the biharmonic problem with Dirichlet boundary conditions in a polygonal domain

In this paper, we study the biharmonic equation with Dirichlet boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem into a system of two Poison equations and one Stokes equation, or a system of one Stokes equation and one Poisson equation. It is shown that the solution of each system is equivalent to that of the original fourth-order problem on both convex and non-convex polygonal domains. Two finite element algorithms are in turn proposed to solve the decoupled systems. In addition, we show the regularity of the solutions in each decoupled system in both the Sobolev space and the weighted Sobolev space, and we derive the optimal error estimates for the numerical solutions on both quasi-uniform meshes and graded meshes. Numerical test results are presented to justify the theoretical findings.

preprint2022arXiv

An adaptive finite element method for two-dimensional elliptic equations with line Dirac sources

In this paper, we propose a novel adaptive finite element method for an elliptic equation with line Dirac delta functions as a source term. We first study the well-posedness and global regularity of the solution in the whole domain. Instead of regularizing the singular source term and using the classical residual-based a posteriori error estimator, we propose a novel a posteriori estimator based on an equivalent transmission problem with zero source term and nonzero flux jumps on line fractures. The transmission problem is defined in the same domain as the original problem excluding on line fractures, and the solution is therefore shown to be more regular. The estimator relies on meshes conforming to the line fractures and its edge jump residual essentially uses the flux jumps of the transmission problem on line fractures. The error estimator is proven to be both reliable and efficient, an adaptive finite element algorithm is proposed based on the error estimator and the bisection refinement method. Numerical tests show that quasi-optimal convergence rates are achieved even for high order approximations and the adaptive meshes are only locally refined at singular points.

preprint2016arXiv

A new anisotropic finite element method on polyhedral domains: interpolation error analysis

Consider the Poisson equation with the Dirichlet boundary condition on a three-dimensional polyhedral domain. For singular solutions from the non-smoothness of the domain boundary, we propose new anisotropic tetrahedral mesh refinement algorithms to improve the convergence of finite element approximation. The proposed algorithm is simple, explicit, and requires less geometric conditions on the mesh and on the domain. Then, we develop interpolation error estimates in suitable weighted spaces for the anisotropic mesh. These estimates can be used to design optimal finite element methods approximating singular solutions. We report numerical test results to validate the method.

preprint2016arXiv

Differential operators on domains with conical points: precise uniform regularity estimates

We study families of strongly elliptic, second order differential operators with singular coefficients on domains with conical points. We obtain uniform estimates on their inverses and on the regularity of the solutions to the associated Poisson problem with mixed boundary conditions. The coefficients and the solutions belong to (suitable) weighted Sobolev spaces. The space of coefficients is a Banach space that contains, in particular, the space of smooth functions. Hence, our results extend classical well-posedness results for strongly elliptic equations in domains with conical points to problems with singular coefficients. We furthermore provide precise uniform estimates on the norms of the solution operators.

preprint2012arXiv

Analysis of Schrödinger operators with inverse square potentials {II}: FEM and approximation of eigenfunctions in the periodic case

Let $V$ be a {\em periodic} potential on $\RR^3$ that is smooth everywhere except at a discrete set $\maS$ of points, where it has singularities of the form $Z/ρ^2$, with $ρ(x) = |x - p|$ for $x$ close to $p$ and $Z$ is continuous, $Z(p) > -1/4$ for $p \in \maS$. We also assume that $ρ$ and $Z$ are smooth outside $\maS$ and $Z$ is smooth in polar coordinates around each singular point. Let us denote by $Λ$ the periodicity lattice and set $\TT := \RR^3/ Λ$. In the first paper of this series \cite{HLNU1}, we obtained regularity results in weighted Sobolev space for the eigenfunctions of the Schrödinger-type operator $H = -Δ+ V$ acting on $L^2(\TT)$, as well as for the induced $\vt k$--Hamiltonians $\Hk$ obtained by resticting the action of $H$ to Bloch waves. In this paper we present two related applications: one to the Finite Element approximation of the solution of $(L+\Hk) v = f$ and one to the numerical approximation of the eigenvalues, $λ$, and eigenfunctions, $u$, of $\Hk$. We give optimal, higher order convergence results for approximation spaces defined piecewise polynomials. Our numerical tests are in good agreement with the theoretical results.

preprint2012arXiv

Analysis of Schrödinger operators with inverse square potentials I: regularity results in 3D

Let $V$ be a potential on $\RR^3$ that is smooth everywhere except at a discrete set $\maS$ of points, where it has singularities of the form $Z/ρ^2$, with $ρ(x) = |x - p|$ for $x$ close to $p$ and $Z$ continuous on $\RR^3$ with $Z(p) > -1/4$ for $p \in \maS$. Also assume that $ρ$ and $Z$ are smooth outside $\maS$ and $Z$ is smooth in polar coordinates around each singular point. We either assume that $V$ is periodic or that the set $\maS$ is finite and $V$ extends to a smooth function on the radial compactification of $\RR^3$ that is bounded outside a compact set containing $\maS$. In the periodic case, we let $Λ$ be the periodicity lattice and define $\TT := \RR^3/ Λ$. We obtain regularity results in weighted Sobolev space for the eigenfunctions of the Schrödinger-type operator $H = -Δ+ V$ acting on $L^2(\TT)$, as well as for the induced $\vt k$--Hamiltonians $\Hk$ obtained by restricting the action of $H$ to Bloch waves. Under some additional assumptions, we extend these regularity and solvability results to the non-periodic case. We sketch some applications to approximation of eigenfunctions and eigenvalues that will be studied in more detail in a second paper.

preprint2012arXiv

Axisymmetric Stokes equations in polygonal domains: regularity and finite element approximations

We study the regularity and finite element approximation of the axisymmetric Stokes problem on a polygonal domain $Ω$. In particular, taking into account the singular coefficients in the equation and non-smoothness of the domain, we establish the well-posedness and full regularity of the solution in new weighted Sobolev spaces $\maK^m_{μ, 1}(Ω)$. Using our a priori results, we give a specific construction of graded meshes on which the Taylor-Hood mixed method approximates singular solutions at the optimal convergence rate. Numerical tests are presented to confirm the theoretical results in the paper.

preprint2012arXiv

Uniform shift estimates for transmission problems and optimal rates of convergence for the parametric Finite Element Method

Let $Ω\subset \RR^d$, $d \geqslant 1$, be a bounded domain with piecewise smooth boundary $\partial Ω$ and let $U$ be an open subset of a Banach space $Y$. Motivated by questions in "Uncertainty Quantification," we consider a parametric family $P = (P_y)_{y \in U}$ of uniformly strongly elliptic, second order partial differential operators $P_y$ on $Ω$. We allow jump discontinuities in the coefficients. We establish a regularity result for the solution $u: Ω\times U \to \RR$ of the parametric, elliptic boundary value/transmission problem $P_y u_y = f_y$, $y \in U$, with mixed Dirichlet-Neumann boundary conditions in the case when the boundary and the interface are smooth and in the general case for $d=2$. Our regularity and well-posedness results are formulated in a scale of broken weighted Sobolev spaces $\hat\maK^{m+1}_{a+1}(Ω)$ of Babuška-Kondrat'ev type in $Ω$, possibly augmented by some locally constant functions. This implies that the parametric, elliptic PDEs $(P_y)_{y \in U}$ admit a shift theorem that is uniform in the parameter $y\in U$. In turn, this then leads to $h^m$-quasi-optimal rates of convergence (i.e. algebraic orders of convergence) for the Galerkin approximations of the solution $u$, where the approximation spaces are defined using the "polynomial chaos expansion" of $u$ with respect to a suitable family of tensorized Lagrange polynomials, following the method developed by Cohen, Devore, and Schwab (2010).