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

Yunrong Zhu

Yunrong Zhu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.APmath-phmath.MPNumerical Analysis

Trust snapshot

Quick read

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

10 published item(s)

preprint2020arXiv

Stabilizer-Free Weak Galerkin Methods for Monotone Quasilinear Elliptic PDEs

In this paper, we study the stabilizer-free weak Galerkin methods on polytopal meshes for a class of second order elliptic boundary value problems of divergence form and with gradient nonlinearity in the principal coefficient. With certain assumptions on the nonlinear coefficient, we show that the discrete problem has a unique solution. This is achieved by showing that the associated operator satisfies certain continuity and monotonicity properties. With the help of these properties, we derive optimal error estimates in the energy norm. We present several numerical examples to verify the error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence of Goal-Oriented Adaptive Finite Element Methods for Semilinear Problems

In this article we develop a convergence theory for goal-oriented adaptive finite element algorithms designed for a class of second-order semilinear elliptic equations. We briefly discuss the target problem class, and introduce several related approximate dual problems that are crucial to both the analysis as well as to the development of a practical numerical method. We then review some standard facts concerning conforming finite element discretization and error-estimate-driven adaptive finite element methods (AFEM). We include a brief summary of a priori estimates for this class of semilinear problems, and then describe some goal-oriented variations of the standard approach to AFEM (GOAFEM). Following the recent approach of Mommer-Stevenson and Holst-Pollock for increasingly general linear problems, we first establish a quasi-error contraction result for the primal problem. We then develop some additional estimates that make it possible to establish contraction of the combined primal-dual quasi-error, and subsequently show convergence with respect to the quantity of interest. Finally, a sequence of numerical experiments are then carefully examined. It is observed that the behavior of the implementation follows the predictions of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multilevel Preconditioners for Reaction-Diffusion Problems with Discontinuous Coefficients

In this paper, we extend some of the multilevel convergence results obtained by Xu and Zhu in [Xu and Zhu, M3AS 2008], to the case of second order linear reaction-diffusion equations. Specifically, we consider the multilevel preconditioners for solving the linear systems arising from the linear finite element approximation of the problem, where both diffusion and reaction coefficients are piecewise-constant functions. We discuss in detail the influence of both the discontinuous reaction and diffusion coefficients to the performance of the classical BPX and multigrid V-cycle preconditioners.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local Multilevel Preconditioners for Elliptic Equations with Jump Coefficients on Bisection Grids

The goal of this paper is to design optimal multilevel solvers for the finite element approximation of second order linear elliptic problems with piecewise constant coefficients on bisection grids. Local multigrid and BPX preconditioners are constructed based on local smoothing only at the newest vertices and their immediate neighbors. The analysis of eigenvalue distributions for these local multilevel preconditioned systems shows that there are only a fixed number of eigenvalues which are deteriorated by the large jump. The remaining eigenvalues are bounded uniformly with respect to the coefficients and the meshsize. Therefore, the resulting preconditioned conjugate gradient algorithm will converge with an asymptotic rate independent of the coefficients and logarithmically with respect to the meshsize. As a result, the overall computational complexity is nearly optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Two-Grid Methods for Semilinear Interface Problems

In this article we consider two-grid finite element methods for solving semilinear interface problems in d space dimensions, for d=2 or d=3. We first describe in some detail the target problem class with discontinuous diffusion coefficients, which includes problems containing sub-critical, critical, and supercritical nonlinearities. We then establish basic quasi-optimal a priori error estimate for Galerkin approximations. In the critical and subcritical cases, we follow our recent approach to controling the nonlinearity using only pointwise control of the continuous solution and a local Lipschitz property, rather than through pointwise control of the discrete solution; this eliminates the requirement that the discrete solution satisfy a discrete form of the maximum principle, hence eliminating the need for restrictive angle conditions in the underlying mesh. The supercritical case continues to require such mesh conditions in order to control the nonlinearity. We then design a two-grid algorithm consisting of a coarse grid solver for the original nonlinear problem, and a fine grid solver for a linearized problem. We analyze the quality of approximations generated by the algorithm, and show that the coarse grid may be taken to have much larger elements than the fine grid, and yet one can still obtain approximation quality that is asymptotically as good as solving the original nonlinear problem on the fine mesh. The algorithm we describe, and its analysis in this article, combines four sets of tools: the work of Xu and Zhou on two-grid algorithms for semilinear problems; the recent results for linear interface problems due to Li, Melenk, Wohlmuth, and Zou; recent work on the Poisson-Boltzmann equation; and recent work on a priori estimates for semilinear problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of a multigrid preconditioner for Crouzeix-Raviart discretization of elliptic PDE with jump coefficient

In this paper, we present a multigrid $V$-cycle preconditioner for the linear system arising from piecewise linear nonconforming Crouzeix-Raviart discretization of second order elliptic problems with jump coefficients. The preconditioner uses standard conforming subspaces as coarse spaces. We showed that the convergence rate of the multigrid $V$-cycle algorithm will deteriorate rapidly due to large jumps in coefficient. However, the preconditioned system has only a fixed number of small eigenvalues, which are deteriorated due to the large jump in coefficient, and the effective condition number is bounded logarithmically with respect to the mesh size. As a result, the multigrid $V$-cycle preconditioned conjugate gradient algorithm converges nearly uniformly. Numerical tests show both robustness with respect to jumps in the coefficient and the mesh size.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite Element Error Estimates for Critical Growth Semilinear Problems without Angle Conditions

In this article we consider a priori error and pointwise estimates for finite element approximations of solutions to semilinear elliptic boundary value problems in d>=2 space dimensions, with nonlinearities satisfying critical growth conditions. It is well-understood how mesh geometry impacts finite element interpolant quality, and leads to the reasonable notion of shape regular simplex meshes. It is also well-known how to perform both mesh generation and simplex subdivision, in arbitrary space dimension, so as to guarantee the entire hierarchy of nested simplex meshes produced through subdivision continue to satisfy shape regularity. However, much more restrictive angle conditions are needed for basic a priori quasi-optimal error estimates, as well as for a priori pointwise estimates. These angle conditions, which are particularly difficult to satisfy in three dimensions in any type of unstructured or adaptive setting, are needed to gain pointwise control of the nonlinearity through discrete maximum principles. This represents a major gap in finite element approximation theory for nonlinear problems on unstructured meshes, and in particular for adaptive methods. In this article, we close this gap in the case of semilinear problems with critical or sub-critical nonlinear growth, by deriving a priori estimates directly, without requiring the discrete maximum principle, and hence eliminating the need for restrictive angle conditions. Our main result is a type of local Lipschitz property that relies only on the continuous maximum principle, together with the growth condition. We also show that under some additional smoothness assumptions, the a priori error estimate itself is enough to give pointwise control the discrete solution, without the need for restrictive angle conditions. Numerical experiments confirm our theoretical conclusions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients

In this paper, we present a multigrid preconditioner for solving the linear system arising from the piecewise linear nonconforming Crouzeix-Raviart discretization of second order elliptic problems with jump coefficients. The preconditioner uses the standard conforming subspaces as coarse spaces. Numerical tests show both robustness with respect to the jump in the coefficient and near-optimality with respect to the number of degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Adaptive Finite Element Modeling Techniques for the Poisson-Boltzmann Equation

We develop an efficient and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the regularization technique of Chen, Holst, and Xu; this technique made possible the first a priori pointwise estimates and the first complete solution and approximation theory for the Poisson-Boltzmann equation. It also made possible the first provably convergent discretization of the PBE, and allowed for the development of a provably convergent AFEM for the PBE. However, in practice the regularization turns out to be numerically ill-conditioned. In this article, we examine a second regularization, and establish a number of basic results to ensure that the new approach produces the same mathematical advantages of the original regularization, without the ill-conditioning property. We then design an AFEM scheme based on the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori pointwise estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. The stability advantages of the new regularization are demonstrated using an FETK-based implementation, through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local Convergence of Adaptive Methods for Nonlinear Partial Differential Equations

In this article we develop convergence theory for a general class of adaptive approximation algorithms for abstract nonlinear operator equations on Banach spaces, and use the theory to obtain convergence results for practical adaptive finite element methods (AFEM) applied to several classes of nonlinear elliptic equations. In the first part of the paper, we develop a weak-* convergence framework for nonlinear operators, whose Gateaux derivatives are locally Lipschitz and satisfy a local inf-sup condition. The framework can be viewed as extending the recent convergence results for linear problems of Morin, Siebert and Veeser to a general nonlinear setting. We formulate an abstract adaptive approximation algorithm for nonlinear operator equations in Banach spaces with local structure. The weak-* convergence framework is then applied to this class of abstract locally adaptive algorithms, giving a general convergence result. The convergence result is then applied to a standard AFEM algorithm in the case of several semilinear and quasi-linear scalar elliptic equations and elliptic systems, including: a semilinear problem with subcritical nonlinearity, the steady Navier-Stokes equations, and a quasilinear problem with nonlinear diffusion. This yields several new AFEM convergence results for these nonlinear problems. In the second part of the paper we develop a second abstract convergence framework based on strong contraction, extending the recent contraction results for linear problems of Cascon, Kreuzer, Nochetto, and Siebert and of Mekchay and Nochetto to abstract nonlinear problems. The contraction result is then applied to a standard AFEM algorithm for semilinear problems with sub- and super-critical nonlinearities and for the Hamiltonian constraint in general relativity.

0 citations0 reviewsNo code linkedOpen access