Researcher profile

Hailiang Liu

Hailiang Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2022arXiv

A global convergence theory for deep ReLU implicit networks via over-parameterization

Implicit deep learning has received increasing attention recently due to the fact that it generalizes the recursive prediction rules of many commonly used neural network architectures. Its prediction rule is provided implicitly based on the solution of an equilibrium equation. Although a line of recent empirical studies has demonstrated its superior performances, the theoretical understanding of implicit neural networks is limited. In general, the equilibrium equation may not be well-posed during the training. As a result, there is no guarantee that a vanilla (stochastic) gradient descent (SGD) training nonlinear implicit neural networks can converge. This paper fills the gap by analyzing the gradient flow of Rectified Linear Unit (ReLU) activated implicit neural networks. For an $m$-width implicit neural network with ReLU activation and $n$ training samples, we show that a randomly initialized gradient descent converges to a global minimum at a linear rate for the square loss function if the implicit neural network is \textit{over-parameterized}. It is worth noting that, unlike existing works on the convergence of (S)GD on finite-layer over-parameterized neural networks, our convergence results hold for implicit neural networks, where the number of layers is \textit{infinite}.

preprint2022arXiv

Global Dynamics and Photon Loss in the Kompaneets Equation

The Kompaneets equation governs dynamics of the photon energy spectrum in certain high temperature (or low density) plasmas. We prove several results concerning the long-time convergence of solutions to Bose--Einstein equilibria and the failure of photon conservation. In particular, we show the total photon number can decrease with time via an outflux of photons at the zero-energy boundary. The ensuing accumulation of photons at zero energy is analogous to Bose--Einstein condensation. We provide two conditions that guarantee that photon loss occurs, and show that once loss is initiated then it persists forever. We prove that as $t\to \infty$, solutions necessarily converge to equilibrium and we characterize the limit in terms of the total photon loss. Additionally, we provide a few results concerning the behavior of the solution near the zero-energy boundary, an Oleinik inequality, a comparison principle, and show that the solution operator is a contraction in $L^1$. None of these results impose a boundary condition at the zero-energy boundary.

preprint2022arXiv

SGEM: stochastic gradient with energy and momentum

In this paper, we propose SGEM, Stochastic Gradient with Energy and Momentum, to solve a large class of general non-convex stochastic optimization problems, based on the AEGD method that originated in the work [AEGD: Adaptive Gradient Descent with Energy. arXiv: 2010.05109]. SGEM incorporates both energy and momentum at the same time so as to inherit their dual advantages. We show that SGEM features an unconditional energy stability property, and derive energy-dependent convergence rates in the general nonconvex stochastic setting, as well as a regret bound in the online convex setting. A lower threshold for the energy variable is also provided. Our experimental results show that SGEM converges faster than AEGD and generalizes better or at least as well as SGDM in training some deep neural networks.

preprint2021arXiv

Energy stable Runge-Kutta discontinuous Galerkin schemes for fourth order gradient flows

We present unconditionally energy stable Runge-Kutta (RK) discontinuous Galerkin (DG) schemes for solving a class of fourth order gradient flows. Our algorithm is geared toward arbitrarily high order approximations in both space and time, while energy dissipation remains preserved without imposing any restriction on time steps and meshes. We achieve this in two steps. First, taking advantage of the penalty free DG method introduced by Liu and Yin [J Sci. Comput. 77:467--501, 2018] for spatial discretization, we reformulate an extended linearized ODE system by the energy quadratization (EQ) approach. Second, we apply an s-stage algebraically stable RK method for temporal discretization. The resulting fully discrete DG schemes are linear and unconditionally energy stable. In addition, we introduce a prediction-correction procedure to improve both the accuracy and stability of the scheme. We illustrate the effectiveness of the proposed schemes by numerical tests with benchmark problems.

preprint2021arXiv

Necessary conditions for blow-up solutions to the restricted Euler--Poisson equations

In this work, we study the behavior of blow-up solutions to the multidimensional restricted Euler--Poisson equations which are the localized version of the full Euler--Poisson system. We provide necessary conditions for the existence of finite-time blow-up solutions in terms of the initial data, and describe the asymptotic behavior of the solutions near blow up times. We also identify a rich set of the initial data which yields global bounded solutions.

preprint2021arXiv

Positivity-preserving third order DG schemes for Poisson--Nernst--Planck equations

In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson--Nernst--Planck (PNP) equations, which has found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.

preprint2021arXiv

Unconditionally energy stable discontinuous Galerkin schemes for the Cahn-Hilliard equation

In this paper, we introduce novel discontinuous Galerkin (DG) schemes for the Cahn-Hilliard equation, which arises in many applications. The method is designed by integrating the mixed DG method for the spatial discretization with the \emph{Invariant Energy Quadratization} (IEQ) approach for the time discretization. Coupled with a spatial projection, the resulting IEQ-DG schemes are shown to be unconditionally energy dissipative, and can be efficiently solved without resorting to any iteration method. Both one and two dimensional numerical examples are provided to verify the theoretical results, and demonstrate the good performance of IEQ-DG in terms of efficiency, accuracy, and preservation of the desired solution properties.

preprint2020arXiv

A positivity-preserving and energy stable scheme for a quantum diffusion equation

We propose a new fully-discretized finite difference scheme for a quantum diffusion equation, in both one and two dimensions. This is the first fully-discretized scheme with proven positivity-preserving and energy stable properties using only standard finite difference discretization. The difficulty in proving the positivity-preserving property lies in the lack of a maximum principle for fourth order PDEs. To overcome this difficulty, we reformulate the scheme as an optimization problem based on variational structure and use the singular nature of the energy functional near the boundary values to exclude the possibility of non-positive solutions. The scheme is also shown to be mass conservative and consistent.

preprint2020arXiv

An energy stable and positivity-preserving scheme for the Maxwell-Stefan diffusion system

We develop a new finite difference scheme for the Maxwell-Stefan diffusion system. The scheme is conservative, energy stable and positivity-preserving. These nice properties stem from a variational structure and are proved by reformulating the finite difference scheme into an equivalent optimization problem. The solution to the scheme emerges as the minimizer of the optimization problem, and as a consequence energy stability and positivity-preserving properties are obtained.

preprint2020arXiv

Dynamics of many species through competition for resources

This paper is concerned with a mathematical model of competition for resource where species consume noninteracting resources. This system of differential equations is formally obtained by renormalizing the MacArthur's competition model at equilibrium, and agrees with the trait-continuous model studied by Mirrahimi S, Perthame B, Wakano JY [J. Math. Biol. 64(7): 1189-1223, 2012]. As a dynamical system, self-organized generation of distinct species occurs. The necessary conditions for survival are given. We prove the existence of the evolutionary stable distribution (ESD) through an optimization problem and present an independent algorithm to compute the ESD directly. Under certain structural conditions, solutions of the system are shown to approach the discrete ESD as time evolves. The time discretization of the system is proven to satisfy two desired properties: positivity and energy dissipation. Numerical examples are given to illustrate certain interesting biological phenomena.

preprint2020arXiv

Efficient, positive, and energy stable schemes for multi-D Poisson-Nernst-Planck systems

In this paper, we design, analyze, and numerically validate positive and energy-dissipating schemes for solving the time-dependent multi-dimensional system of Poisson-Nernst-Planck (PNP) equations, which has found much use in the modeling of biological membrane channels and semiconductor devices. The semi-implicit time discretization based on a reformulation of the system gives a well-posed elliptic system, which is shown to preserve solution positivity for arbitrary time steps. The first order (in time) fully-discrete scheme is shown to preserve solution positivity and mass conservation unconditionally, and energy dissipation with only a mild $O(1)$ time step restriction. The scheme is also shown to preserve the steady-state. For the fully second order (in both time and space) scheme with large time steps, solution positivity is restored by a local scaling limiter, which is shown to maintain the spatial accuracy. These schemes are easy to implement. Several three-dimensional numerical examples verify our theoretical findings and demonstrate the accuracy, efficiency, and robustness of the proposed schemes, as well as the fast approach to steady states.

preprint2020arXiv

On the SAV-DG method for a class of fourth order gradient flows

For a class of fourth order gradient flow problems, integration of the scalar auxiliary variable (SAV) time discretization with the penalty-free discontinuous Galerkin (DG) spatial discretization leads to SAV-DG schemes. These schemes are linear and shown unconditionally energy stable. But the reduced linear systems are rather expensive to solve due to the dense coefficient matrices. In this paper, we provide a procedure to pre-evaluate the auxiliary variable in the piecewise polynomial space. As a result, the computational complexity of $O(\mathcal{N}^2)$ reduces to $O(\mathcal{N})$ when exploiting the conjugate gradient (CG) solver. This hybrid SAV-DG method is more efficient and able to deliver satisfactory results of high accuracy. This was also compared with solving the full augmented system of the SAV-DG schemes.

preprint2020arXiv

Radially symmetric solutions of the ultra-relativistic Euler equations

The ultra-relativistic Euler equations for an ideal gas are described in terms of the pressure $p$, the spatial part $\underline{u} \in \R^3$ of the dimensionless four-velocity and the particle density $n$. Radially symmetric solutions of these equations are studied. Analytical solutions are presented for the linearized system. For the original nonlinear equations we design and analyze a numerical scheme for simulating radially symmetric solutions in three space dimensions. The good performance of the scheme is demonstrated by numerical examples. In particular, it was observed that the method has the capability to capture accurately the pressure singularity formation caused by shock wave reflections at the origin.

preprint2020arXiv

Selection dynamics for deep neural networks

This paper presents a partial differential equation framework for deep residual neural networks and for the associated learning problem. This is done by carrying out the continuum limits of neural networks with respect to width and depth. We study the wellposedness, the large time solution behavior, and the characterization of the steady states of the forward problem. Several useful time-uniform estimates and stability/instability conditions are presented. We state and prove optimality conditions for the inverse deep learning problem, using standard variational calculus, the Hamilton-Jacobi-Bellmann equation and the Pontryagin maximum principle. This serves to establish a mathematical foundation for investigating the algorithmic and theoretical connections between neural networks, PDE theory, variational analysis, optimal control, and deep learning.

preprint2019arXiv

Positive and free energy satisfying schemes for diffusion with interaction potentials

In this paper, we design and analyze second order positive and free energy satisfying schemes for solving diffusion equations with interaction potentials. The semi-discrete scheme is shown to conserve mass, preserve solution positivity, and satisfy a discrete free energy dissipation law for nonuniform meshes. These properties for the fully-discrete scheme (first order in time) remain preserved without a strict restriction on time steps. For the fully second order (in both time and space) scheme, we use a local scaling limiter to restore solution positivity when necessary. It is proved that such limiter does not destroy the second order accuracy. In addition, these schemes are easy to implement, and efficient in simulations over long time. Both one and two dimensional numerical examples are presented to demonstrate the performance of these schemes.

preprint2010arXiv

Global Well-Posedness for the Microscopic FENE Model with a Sharp Boundary Condition

We prove global well-posedness for the microscopic FENE model under a sharp boundary requirement. The well-posedness of the FENE model that consists of the incompressible Navier-Stokes equation and the Fokker-Planck equation has been studied intensively, mostly with the zero flux boundary condition. Recently it was illustrated by C. Liu and H. Liu [2008, SIAM J. Appl. Math., 68(5):1304--1315] that any preassigned boundary value of a weighted distribution will become redundant once the non-dimensional parameter $b>2$. In this article, we show that for the well-posedness of the microscopic FENE model ($b>2$) the least boundary requirement is that the distribution near boundary needs to approach zero faster than the distance function. Under this condition, it is shown that there exists a unique weak solution in a weighted Sobolev space. Moreover, such a condition still ensures that the distribution is a probability density. The sharpness of this boundary requirement is shown by a construction of infinitely many solutions when the distribution approaches zero as fast as the distance function.