BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zhaonan Dong

Zhaonan Dong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

7works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

7 published item(s)

preprint2022arXiv

Hybrid high-order and weak Galerkin methods for the biharmonic problem

We devise and analyze two hybrid high-order (HHO) methods for the numerical approximation of the biharmonic problem. The methods support polyhedral meshes, rely on the primal formulation of the problem, and deliver $O(h^{k+1})$ $H^2$-error estimates when using polynomials of order $k\ge0$ to approximate the normal derivative on the mesh (inter)faces. Both HHO methods hinge on a stabilization in the spirit of Lehrenfeld--Schöberl for second-order PDEs. The cell unknowns are polynomials of order $(k+2)$ that can be eliminated locally by means of static condensation. The face unknowns approximating the trace of the solution on the mesh (inter)faces are polynomials of order $(k+1)$ in the first HHO method which is valid in dimension two and uses an original stabilization involving the canonical hybrid finite element, and they are of order $(k+2)$ for the second HHO method which is valid in arbitrary dimension and uses only $L^2$-orthogonal projections in the stabilization. A comparative discussion with the weak Galerkin methods from the literature is provided, highlighting the close connections and the improvements proposed herein. Additionally, we show how the two HHO methods can be combined with a Nitsche-like boundary-penalty technique to weakly enforce the boundary conditions. An originality in the devised Nitsche's technique is to avoid any penalty parameter that must be large enough. Finally, numerical results showcase the efficiency of the proposed methods and indicate that the HHO methods can generally outperform discontinuous Galerkin methods and even be competitive with $C^0$-interior penalty methods on triangular meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local decay rates of best-approximation errors using vector-valued finite elements for fields with low regularity and integrable curl or divergence

We estimate best-approximation errors using vector-valued finite elements for fields with low regularity in the scale of fractional-order Sobolev spaces. By assuming additionally that the target field has a curl or divergence property, we establish upper bounds on these errors that can be localized to the mesh cells. These bounds are derived using the quasi-interpolation errors with or without boundary prescription derived in [A. Ern and J.-L. Guermond, ESAIM Math. Model. Numer. Anal., 51 (2017), pp.~1367--1385]. By using the face-to-cell lifting operators analyzed in [A. Ern and J.-L. Guermond, Found. Comput. Math., (2021)], and exploiting the additional assumption made on the curl or the divergence of the target field, a localized upper bound on the quasi-interpolation error is derived. As an illustration, we show how to apply these results to the error analysis of the curl-curl problem associated with Maxwell's equations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust interior penalty discontinuous Galerkin methods

Classical interior penalty discontinuous Galerkin (IPDG) methods for diffusion problems require a number of assumptions on the local variation of mesh-size, polynomial degree, and of the diffusion coefficient to determine the values of the, so-called, discontinuity-penalization parameter and/or to perform error analysis. Variants of IPDG methods involving weighted averages of the gradient of the approximate solution have been proposed in the context of high-contrast diffusion coefficients to mitigate the dependence of the contrast in the stability and in the error analysis. Here, we present a new IPDG method, involving carefully constructed weighted averages of the gradient of the approximate solution, which is shown to be robust even for the most extreme simultaneous local mesh, polynomial degree and diffusion coefficient variation scenarios, without resulting in unreasonably large penalization. The new method, henceforth termed as \emph{robust IPDG} (RIPDG), offers typically significantly better conditioning than the standard IPDG method when applied to scenarios with strong mesh/polynomial degree/diffusion local variation. On the other hand, when using uniform meshes, constant polynomial degree and for problems with constant diffusion coefficients, the RIPDG method is identical to the classical IPDG. Numerical experiments indicate the favourable performance of the new RIPDG method over the classical version in terms of conditioning and error.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

GPU-accelerated discontinuous Galerkin methods on polytopic meshes

Discontinuous Galerkin (dG) methods on meshes consisting of polygonal/polyhedral (henceforth, collectively termed as \emph{polytopic}) elements have received considerable attention in recent years. Due to the physical frame basis functions used typically and the quadrature challenges involved, the matrix-assembly step for these methods is often computationally cumbersome. To address this important practical issue, this work proposes two parallel assembly implementation algorithms on CUDA-enabled graphics cards for the interior penalty dG method on polytopic meshes for various classes of linear PDE problems. We are concerned with both single GPU parallelization, as well as with implementation on distributed GPU nodes. The results included showcase almost linear scalability of the quadrature step with respect to the number of GPU-cores used since no communication is needed for the assembly step. In turn, this can justify the claim that polytopic dG methods can be implemented extremely efficiently, as any assembly computing time overhead compared to finite elements on `standard' simplicial or box-type meshes can be effectively circumvented by the proposed algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the suboptimality of the p-version discontinuous Galerkin methods for first order hyperbolic problems

We address the issue of the suboptimality in the p-version discontinuous Galerkin (dG) methods for first order hyperbolic problems. The convergence rate is derived for the upwind dG scheme on tensor product meshes in any dimension. The standard proof in seminal work [14] leads to suboptimal convergence in terms of the polynomial degree by 3/2 order for general convection fields, with the exception of piecewise multi-linear convection fields, which rather yield optimal convergence. Such suboptimality is not observed numerically. Thus, it might be caused by a limitation of the analysis, which we partially overcome: for a special class of convection fields, we shall show that the dG method has a p-convergence rate suboptimal by 1/2 order only.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Residual-based a posteriori error estimates for $\mathbf{hp}$-discontinuous Galerkin discretisations of the biharmonic problem

We introduce a residual-based a posteriori error estimator for a novel $hp$-version interior penalty discontinuous Galerkin method for the biharmonic problem in two and three dimensions. We prove that the error estimate provides an upper bound and a local lower bound on the error, and that the lower bound is robust to the local mesh size but not the local polynomial degree. The suboptimality in terms of the polynomial degree is fully explicit and grows at most algebraically. Our analysis does not require the existence of a $\mathcal{C}^1$-conforming piecewise polynomial space and is instead based on an elliptic reconstruction of the discrete solution to the $H^2$ space and a generalised Helmholtz decomposition of the error. This is the first $hp$-version error estimator for the biharmonic problem in two and three dimensions. The practical behaviour of the estimator is investigated through numerical examples in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast multilevel sparse Gaussian kernels for high-dimensional approximation and integration

A fast multilevel algorithm based on directionally scaled tensor-product Gaussian kernels on structured sparse grids is proposed for interpolation of high-dimensional functions and for the numerical integration of high-dimensional integrals. The algorithm is based on the recent Multilevel Sparse Kernel-based Interpolation (MLSKI) method (Georgoulis, Levesley \& Subhan, \emph{SIAM J. Sci. Comput.}, 35(2), pp.~A815--A831, 2013), with particular focus on the fast implementation of Gaussian-based MLSKI for interpolation and integration problems of high-dimen-sional functions $f:[0,1]^d\to\mathbb{R}$, with $5\le d\le 10$. The MLSKI interpolation procedure is shown to be interpolatory and a fast implementation is proposed. More specifically, exploiting the tensor-product nature of anisotropic Gaussian kernels, one-dimensional cardinal basis functions on a sequence of hierarchical equidistant nodes are precomputed to machine precision, rendering the interpolation problem into a fully parallelisable ensemble of linear combinations of function evaluations. A numerical integration algorithm is also proposed, based on interpolating the (high-dimensional) integrand. A series of numerical experiments highlights the applicability of the proposed algorithm for interpolation and integration for up to 10-dimensional problems.

0 citations0 reviewsNo code linkedOpen access