NAApr 24, 2018
A Unified Study of Continuous and Discontinuous Galerkin MethodsQingguo Hong, Fei Wang, Shuonan Wu et al.
A unified study is presented in this paper for the design and analysis of different finite element methods (FEMs), including conforming and nonconforming FEMs, mixed FEMs, hybrid FEMs,discontinuous Galerkin (DG) methods, hybrid discontinuous Galerkin (HDG) methods and weak Galerkin (WG) methods. Both HDG and WG are shown to admit inf-sup conditions that hold uniformly with respect to both mesh and penalization parameters. In addition, by taking the limit of the stabilization parameters, a WG method is shown to converge to a mixed method whereas an HDG method is shown to converge to a primal method. Furthermore, a special class of DG methods, known as the mixed DG methods, is presented to fill a gap revealed in the unified framework.
LGAug 9, 2022
On the Activation Function Dependence of the Spectral Bias of Neural NetworksQingguo Hong, Jonathan W. Siegel, Qinyang Tan et al.
Neural networks are universal function approximators which are known to generalize well despite being dramatically overparameterized. We study this phenomenon from the point of view of the spectral bias of neural networks. Our contributions are two-fold. First, we provide a theoretical explanation for the spectral bias of ReLU neural networks by leveraging connections with the theory of finite element methods. Second, based upon this theory we predict that switching the activation function to a piecewise linear B-spline, namely the Hat function, will remove this spectral bias, which we verify empirically in a variety of settings. Our empirical studies also show that neural networks with the Hat activation function are trained significantly faster using stochastic gradient descent and ADAM. Combined with previous work showing that the Hat activation function also improves generalization accuracy on image classification tasks, this indicates that using the Hat activation provides significant advantages over the ReLU on certain problems.
NADec 31, 2018
Two-grid economical algorithms for parabolic integro-differential equations with nonlinear memoryWansheng Wang, Qingguo Hong
In this paper, several two-grid finite element algorithms for solving parabolic integro-differential equations (PIDEs) with nonlinear memory are presented. Analysis of these algorithms is given assuming a fully implicit time discretization. It is shown that these algorithms are as stable as the standard fully discrete finite element algorithm, and can achieve the same accuracy as the standard algorithm if the coarse grid size $H$ and the fine grid size $h$ satisfy $H=O(h^{\frac{r-1}{r}})$. Especially for PIDEs with nonlinear memory defined by a lower order nonlinear operator, our two-grid algorithm can save significant storage and computing time. Numerical experiments are given to confirm the theoretical results.
55.8NAApr 27
A correction adaptive two-grid finite element method for nonselfadjoint or indefinite elliptic problemsFei Li, Qingguo Hong, Ming Tang et al.
We propose, analyze, and numerically validate a correction adaptive two-grid finite element method (CAT-GFEM) for nonselfadjoint or indefinite elliptic problems. In contrast to the adaptive two-grid finite element method (ATGFEM) of Li and Zhang [SIAM J. Sci. Comput., 43 (2021), pp. A908-A928], which is restricted to symmetric positive-definite problems, the proposed method introduces an additional correction step that solves a small-scale discrete residual problem on the coarse mesh. This step entails negligible additional computational cost and allows us to show that the L2-norm error of the corrected discrete solution is a higher-order of the energy-norm error of the discrete solution. Using this result, we prove a contraction property for a suitable sum of quasi-errors on two successive adaptive meshes and establish convergence of the method. Numerical experiments illustrate the improved effectiveness and robustness of our method in comparison with ATGFEM.
NAMar 29, 2019
Parameter-robust convergence analysis of fixed-stress split iterative method for multiple-permeability poroelasticity systemsQingguo Hong, Johannes Kraus, Maria Lymbery et al.
We consider flux-based multiple-porosity/multiple-permeability poroelasticity systems describing multiple-network flow and deformation in a poro-elastic medium, sometimes also referred to as MPET models. The focus of the paper is on the convergence analysis of the fixed-stress split iteration, a commonly used coupling technique for the flow and mechanics equations in poromechanics. We formulate the fixed-stress split method in the present context and prove its linear convergence. The contraction rate of this fixed-point iteration does not depend on any of the physical parameters appearing in the model. The theory is confirmed by numerical results which further demonstrate the advantage of the fixed-stress split scheme over a fully implicit method relying on norm-equivalent preconditioning.
NAOct 7, 2018
Uniform Stability and Error Analysis for Some Discontinuous Galerkin MethodsQingguo Hong, Jinchao Xu
In this paper, we provide a number of new estimates on the stability and convergence of both hybrid discontinuous Galerkin (HDG) and weak Galerkin (WG) methods. By using the standard Brezzi theory on mixed methods, we carefully define appropriate norms for the various discretization variables and then establish that the stability and error estimates hold uniformly with respect to stabilization and discretization parameters. As a result, by taking appropriate limit of the stabilization parameters, we show that the HDG method converges to a primal conforming method and the WG method converge to a mixed conforming method.
NAJun 13, 2017
Parameter-robust stability of classical three-field formulation of Biot's consolidation modelQingguo Hong, Johannes Kraus
This paper is devoted to the stability analysis of a classical three-field formulation of Biot's consolidation model where the unknown variables are the displacements, fluid flux (Darcy velocity), and pore pressure. Specific parameter-dependent norms provide the key in establishing the full parameter-robust inf-sup stability of the continuous problem. Therefore, stability results presented here are uniform not only with respect to the Lamé parameter $λ$, but also with respect to all the other model parameters. This allows for the construction of a uniform block diagonal preconditioner within the framework of operator preconditioning. Stable discretizations that meet the required conditions for full robustness and guarantee mass conservation, both locally and globally, are discussed and corresponding optimal error estimates proven.