BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Zhaonan Dong

Zhaonan Dong contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
2topics
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

6 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