39.5NAMay 20
Quasi-optimal polytopal finite element methods for biharmonic equationNgoc Tien Tran
This paper establishes quasi-optimal and lower-order error estimates for weak Galerkin, discontinuous Galerkin, and hybrid-high order finite element methods for the biharmonic equation under minimal regularity assumptions on general polytopal meshes. Furthermore, it is shown that the stabilization is an efficient contribution in a~posteriori error estimators.
NADec 12, 2024
Locking-free hybrid high-order method for linear elasticityCarsten Carstensen, Ngoc Tien Tran
The hybrid-high order (HHO) scheme has many successful applications including linear elasticity as the first step towards computational solid mechanics. The striking advantage is the simplicity among other higher-order nonconforming schemes and its geometric flexibility as a polytopal method on the expanse of a parameter-free refined stabilization. This paper utilizes just one reconstruction operator for the linear Green strain and therefore does not rely on a split in deviatoric and spherical behaviour as in the classical HHO discretization. The a priori error analysis provides quasi-best approximation with $λ$-independent equivalence constants. The reliable and (up to data oscillations) efficient a posteriori error estimates are stabilization-free and $λ$-robust. The error analysis is carried out on simplicial meshes to allow conforming piecewise polynomials finite elements in the kernel of the stabilization terms. Numerical benchmarks provide empirical evidence for optimal convergence rates of the a posteriori error estimator in some associated adaptive mesh-refining algorithm also in the incompressible limit, where this paper provides corresponding assertions for the Stokes problem.
NAApr 5, 2024
Discrete weak duality of hybrid high-order methods for convex minimization problemsNgoc Tien Tran
This paper derives a discrete dual problem for a prototypical hybrid high-order method for convex minimization problems. The discrete primal and dual problem satisfy a weak convex duality that leads to a priori error estimates with convergence rates under additional smoothness assumptions. This duality holds for general polyhedral meshes and arbitrary polynomial degrees of the discretization. A novel postprocessing is proposed and allows for a~posteriori error estimates on regular triangulations into simplices using primal-dual techniques. This motivates an adaptive mesh-refining algorithm, which performs superiorly compared to uniform mesh refinements.
NAFeb 2, 2024
Finite element approximation for uniformly elliptic linear PDE of second order in nondivergence formNgoc Tien Tran
This paper proposes a novel technique for the approximation of strong solutions $u \in C(\overlineΩ) \cap W^{2,n}_\mathrm{loc}(Ω)$ to uniformly elliptic linear PDE of second order in nondivergence form with continuous leading coefficient in nonsmooth domains by finite element methods. These solutions satisfy the Alexandrov-Bakelman-Pucci (ABP) maximum principle, which provides an a~posteriori error control for $C^1$ conforming approximations. By minimizing this residual, we obtain an approximation to the solution $u$ in the $L^\infty$ norm. Although discontinuous functions do not satisfy the ABP maximum principle, this approach extends to nonconforming FEM as well thanks to well-established enrichment operators. Convergence of the proposed FEM is established for uniform mesh-refinements. The built-in a~posteriori error control (even for inexact solve) can be utilized in adaptive computations for the approximation of singular solutions, which performs superiorly in the numerical benchmarks in comparison to the uniform mesh-refining algorithm.
52.1NAApr 3
A hybrid high-order method for the biharmonic problemYizhou Liang, Ngoc Tien Tran
This paper proposes a new hybrid high-order discretization for the biharmonic problem and the corresponding eigenvalue problem. The discrete ansatz space includes degrees of freedom in $n-2$ dimensional submanifolds (e.g., nodal values in 2D and edge values in 3D), in addition to the typical degrees of freedom in the mesh and on the hyperfaces in the HHO literature. This approach enables the characteristic commuting property of the hybrid high-order methodology in any space dimension. The main results are guaranteed lower eigenvalue bounds of higher order. Furthermore, we derive quasi-best approximation estimates as well as reliable and efficient a~posteriori error estimators under minimal regularity assumptions on the exact solution. The latter motivates an adaptive mesh-refining algorithm that empirically recovers optimal convergence rates for singular solutions.
NADec 1, 2025
Lower eigenvalue bounds with hybrid high-order methodsNgoc Tien Tran
This paper proposes hybrid high-order eigensolvers for the computation of guaranteed lower eigenvalue bounds. These bounds display higher order convergence rates and are accessible to adaptive mesh-refining algorithms. The involved constants arise from local embeddings and are available for all polynomial degrees. Applications include the linear elasticity and Steklov eigenvalue problem.
7.4NAApr 9
Local discontinuous Galerkin FEM for convex minimizationCarsten Carstensen, Ngoc Tien Tran
The heart of the a priori and a posteriori error control in convex minimization problems is the sharp control of the approximation of the respective discrete and exact minimal energies. Conforming finite element discretizations for p-Laplace type minimization problems provide upper bounds of the energy difference with optimal convergence rates. Proven convergence rates for higher-order non-conforming finite element discretizations for the same problem class, however, are exclusively suboptimal. Thus the popular a posteriori error control within the two-energy principle, that generalize hyper-circle identities, appears unbalanced. The innovative point of departure in a refined analysis of two discontinuous Galerkin (dG) schemes exploits duality relations between a discrete primal and a semi-discrete dual problem. The infinite-dimensional dual problem leads to a tiny duality gap that even vanishes for polynomial low-order terms. For a class of degenerated convex minimization problems with two-sided $p$ growth, the novel duality provides improved a priori convergence rates for the error in the minimal energies. The motivating two-energy principle and some post-processing for a Raviart-Thomas dual variable provides an a posteriori error control, that also may drive adaptive mesh-refining. Computational benchmarks provide striking numerical evidence for improved convergence rates of the adaptive beyond uniform mesh-refining.