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

Shangyou Zhang

Shangyou Zhang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

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

25 published item(s)

preprint2026arXiv

A Simple Weak Galerkin Finite Element Method for Convection-Diffusion-Reaction Equations on Nonconvex Polytopal Meshes

This article introduces a simple weak Galerkin (WG) finite element method for solving convection-diffusion-reaction equation. The proposed method offers significant flexibility by supporting discontinuous approximating functions on general nonconvex polytopal meshes. We establish rigorous error estimates within a suitable norm. Finally, numerical experiments are presented to validate the theoretical convergence rates and demonstrate the computational efficiency of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

A weak Galerkin least squares finite element method for linear convection equations in non-divergence form

This article develops a weak Galerkin least-squares (WG--LS) finite element method for first-order linear convection equations in non-divergence form. The method is formulated using discontinuous finite element functions and does not require any coercivity assumption on the convection vector or reaction coefficient. The resulting discrete problem leads to a symmetric and positive definite linear system and is applicable to general polygonal and polyhedral meshes. Under minimal regularity assumptions on the coefficients, optimal-order error estimates are established for the WG--LS approximation in a suitable energy norm. Numerical experiments are presented to confirm the theoretical convergence results and to demonstrate the accuracy and efficiency of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A parallel iterative procedure for weak Galerkin methods for second order elliptic problems

A parallelizable iterative procedure based on domain decomposition is presented and analyzed for weak Galerkin finite element methods for second order elliptic equations. The convergence analysis is established for the decomposition of the domain into individual elements associated to the weak Galerkin methods or into larger subdomains. A series of numerical tests are illustrated to verify the theory developed in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Integral Representation of Hydraulic Permeability

In this paper, we show that the permeability of a porous {material} and that of a bubbly fluid are limiting cases of the complexified version of the two-fluid models posed in {Lipton_Avellaneda_1990}. We assume the viscosity of the inclusion fluid is $zμ_1$ and the viscosity of the hosting fluid is $μ_1\in \mathbb{R}^+$, $z\in\mathbb{C}$. The proof is carried out by the construction of solutions for large $|z|$ and small $|z|$ with an iteration process and the analytic continuation. Moreover, we also show that for a fixed microstructure, the permeabilities of these three cases share the same integral representation formula (IRF) in Equation (99) with different values of contrast parameter $s:=1/(z-1)$, as long as $s$ is outside the interval $[-\frac{2E_2^2}{1+2E_2^2},-\frac{1}{1+2E_1^2}]$, where the positive constants $E_1$ and $E_2$ are the extension constants that depend only on the geometry of the periodic pore space of the material.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The nodal basis of $C^m$-$P_{k}^{(3)}$ and $C^m$-$P_{k}^{(4)}$ finite elements on tetrahedral and 4D simplicial grids

We construct the nodal basis of $C^m$-$P_{k}^{(3)}$ ($k \ge 2^3m+1$) and $C^m$-$P_{k}^{(4)}$ ($k \ge 2^4m+1$) finite elements on 3D tetrahedral and 4D simplicial grids, respectively. $C^m$-$P_{k}^{(n)}$ stands for the space of globally $C^m$ ($m\ge1$) and locally piecewise $n$-dimensional polynomials of degree $k$ on $n$-dimensional simplicial grids. We prove the uni-solvency and the $C^m$ continuity of the constructed $C^m$-$P_{k}^{(3)}$ and $C^m$-$P_{k}^{(4)}$ finite element spaces. A computer code is provided which generates the index set for the nodal basis of $C^m$-$P_k^{(n)}$ finite elements on $n$-dimensional simplicial grids.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A weak Galerkin-mixed finite element method for the Stokes-Darcy problem

In this paper, we propose a new numerical scheme for the coupled Stokes-Darcy model with Beavers-Joseph-Saffman interface condition. We use the weak Galerkin method to discretize the Stokes equation and the mixed finite element method to the Darcy equation. A discrete inf-sup condition is proved and optimal error estimates are also derived. Numerical experiments validate the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new weak gradient for the stabilizer free weak Galerkin method with polynomial reduction

The weak Galerkin (WG) finite element method is an effective and flexible general numerical technique for solving partial differential equations. It is a natural extension of the classic conforming finite element method for discontinuous approximations, which maintains simple finite element formulation. Stabilizer free weak Galerkin methods further simplify the WG methods and reduce computational complexity. This paper explores the possibility of optimal combination of polynomial spaces that minimize the number of unknowns in the stabilizer free WG schemes without compromising the accuracy of the numerical approximation. A new stabilizer free weak Galerkin finite element method is proposed and analyzed with polynomial degree reduction. To achieve such a goal, a new definition of weak gradient is introduced. Error estimates of optimal order are established for the corresponding WG approximations in both a discrete $H^1$ norm and the standard $L^2$ norm. The numerical examples are tested on various meshes and confirm the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stabilizer free weak Galerkin element method with supercloseness of order two

The weak Galerkin (WG) finite element method is an effective and flexible general numerical techniques for solving partial differential equations. A simple weak Galerkin finite element method is introduced for second order elliptic problems. First we have proved that stabilizers are no longer needed for this WG element. Then we have proved the supercloseness of order two for the WG finite element solution. The numerical results confirm the theory

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stabilizer free weak Galerkin finite element method on polytopal mesh: Part II

A stabilizer free weak Galerkin (WG) finite element method on polytopal mesh has been introduced in Part I of this paper (J. Comput. Appl. Math, 371 (2020) 112699. arXiv:1906.06634.) Removing stabilizers from discontinuous finite element methods simplifies formulations and reduces programming complexity. The purpose of this paper is to introduce a new WG method without stabilizers on polytopal mesh that has convergence rates one order higher than optimal convergence rates. This method is the first WG method that achieves superconvergence on polytopal mesh. Numerical examples in 2D and 3D are presented verifying the theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stabilizer free weak Galerkin finite element method on polytopal mesh: Part III

A weak Galerkin (WG) finite element method without stabilizers was introduced in [J. Comput. Appl. Math., 371 (2020). arXiv:1906.06634] on polytopal mesh. Then it was improved in [arXiv:2008.13631] with order one superconvergence. The goal of this paper is to develop a new stabilizer free WG method on polytopal mesh. This method has convergence rates two orders higher than the optimal convergence rates for the corresponding WG solution in both an energy norm and the $L^2$ norm. The numerical examples are tested for low and high order elements in two and three dimensional spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stabilizer-free pressure-robust finite element method for the Stokes equations

In this paper, we introduce a new finite element method for solving the Stokes equations in the primary velocity-pressure formulation. This method employs $H(div)$ finite elements to approximate velocity, which leads to two unique advantages: exact divergence free velocity field and pressure-robustness. In addition, this method has a simple formulation without any stabilizer or penalty term. Optimal-order error estimates are established for the corresponding numerical approximation in various norms. Extensive numerical investigations are conducted to test accuracy and robustness of the method and to confirm the theory. The numerical examples cover low and high order approximations up to the degree four, and 2D and 3D cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

De Rham Complexes for Weak Galerkin Finite Element Spaces

Two de Rham complex sequences of the finite element spaces are introduced for weak finite element functions and weak derivatives developed in the weak Galerkin (WG) finite element methods on general polyhedral elements. One of the sequences uses polynomials of equal order for all the finite element spaces involved in the sequence and the other one uses polynomials of naturally decending orders. It is shown that the diagrams in both de Rham complexes commute for general polyhedral elements. The exactness of one of the complexes is established for the lowest order element.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical investigation on weak Galerkin finite elements

The weak Galerkin (WG) finite element method is an effective and flexible general numerical technique for solving partial differential equations. The novel idea of weak Galerkin finite element methods is on the use of weak functions and their weak derivatives defined as distributions. Weak functions and weak derivatives can be approximated by polynomials with various degrees. Different combination of polynomial spaces generates different weak Galerkin finite elements. The purpose of this paper is to study stability, convergence and supercloseness of different WG elements by providing many numerical experiments recorded in 31 tables. These tables serve two purposes. First it provides a detail guide of the performance of different WG elements. Second, the information in the tables opens new research territory why some WG elements outperform others.

0 citations0 reviewsNo code linkedOpen access
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
preprint2015arXiv

A family of conforming mixed finite elements for linear elasticity on triangular grids

This paper presents a family of mixed finite elements on triangular grids for solving the classical Hellinger-Reissner mixed problem of the elasticity equations. In these elements, the matrix-valued stress field is approximated by the full $C^0$-$P_k$ space enriched by $(k-1)$ $H(\d)$ edge bubble functions on each internal edge, while the displacement field by the full discontinuous $P_{k-1}$ vector-valued space, for the polynomial degree $k\ge 3$. The main challenge is to find the correct stress finite element space matching the full $C^{-1}$-$P_{k-1}$ displacement space. The discrete stability analysis for the inf-sup condition does not rely on the usual Fortin operator, which is difficult to construct. It is done by characterizing the divergence of local stress space which covers the $P_{k-1}$ space of displacement orthogonal to the local rigid-motion. The well-posedness condition and the optimal a priori error estimate are proved for this family of finite elements. Numerical tests are presented to confirm the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite element approximations of symmetric tensors on simplicial grids in Rn: the lower order case

In this paper, we construct, in a unified fashion, lower order finite element subspaces of spaces of symmetric tensors with square-integrable divergence on a domain in any dimension. These subspaces are essentially the symmetric H(div)-Pk (1=<k<=n) tensor spaces, enriched, for each n-1 dimensional simplex, by (n+1)n/2 H(div)-Pn+1 bubble functions when 1=< k<= n-1, and by (n-1)n/2 H(div)-P n+1 bubble functions when k= n. These spaces can be used to approximate the symmetric matrix field in a mixed formulation problem where the other variable is approximated by discontinuous piecewise Pk-1 polynomials. This in particular leads to first order mixed elements on simplicial grids with total degrees of freedom per element $18$ plus $3$ in 2D, 48 plus 6 in 3D. The previous record of the degrees of freedom of first order mixed elements is, 21 plus 3 in 2D, and 156 plus 6 in 3D, on simplicial grids. We also derive, in a unified way and without using any tools like differential forms, a family of auxiliary mixed finite elements in any dimension. One example in this family is the Raviart-Thomas elements in one dimension, the second example is the mixed finite elements for linear elasticity in two dimensions due to Arnold and Winther, the third example is the mixed finite elements for linear elasticity in three dimensions due to Arnold, Awanou and Winther.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superconvergence of simple conforming mixed finite elements for linear elasticity on rectangular grids in any space dimension

This paper is to prove superconvergence of a family of simple conforming mixed finite elements of first orderfor the linear elasticity problem with the Hellinger--Reissner variational formulation. The analysis is based on three main ingredients: a new interpolation operator, a new expansion method, and a new iterative argument for superconvergence analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonconforming finite element methods on quadrilateral meshes

It is well-known that it is comparatively difficult to design nonconforming finite elements on quadrilateral meshes by using Gauss-Legendre points on each edge of triangulations. One reason lies in that these degrees of freedom associated to these Gauss-Legendre points are not all linearly independent for usual expected polynomial spaces, which explains why only several lower order nonconforming quadrilateral finite elements can be found in literature. The present paper proposes two families of nonconforming finite elements of any odd order and one family of nonconforming finite elements of any even order on quadrilateral meshes. Degrees of freedom are given for these elements, which are proved to be well-defined for their corresponding shape function spaces in a unifying way. These elements generalize three lower order nonconforming finite elements on quadrilaterals to any order. In addition, these nonconforming finite element spaces are shown to be full spaces which is somehow not discussed for nonconforming finite elements in literature before.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The simplest mixed finite element method for linear elasticity in the symmetric formulation on $n$-rectangular grids

A family of mixed finite elements is proposed for solving the first order system of linear elasticity equations in any space dimension, where the stress field is approximated by symmetric finite element tensors. This family of elements has a perfect matching between the stress components and the displacement. The discrete spaces for the normal stress $σ_{ii}$, the shear stress $σ_{ij}$ and the displacement $u_i$ are $\operatorname{span}\{1,x_i\}$, $\operatorname{span}\{1,x_i,x_j\}$ and $\operatorname{span}\{1\}$, respectively, on rectangular grids. In particular, the definition remains the same for all space dimensions. As a result of these choices, the theoretical analysis is independent of the spatial dimension as well. In 1D, this element is nothing else but the 1D Raviart-Thomas element, which is the only conforming element in this family. In 2D and higher dimensions, they are new elements but of the minimal degrees of freedom. The total degrees of freedom per element is 2 plus 1 in 1D, 7 plus 2 in 2D, and 15 plus 3 in 3D. The previous record of the least degrees of freedom is, 13 plus 4 in 2D, and 54 plus 12 in 3D, on the rectangular grid. These elements are the simplest element for any space dimension. The well-posedness condition and the optimal a priori error estimate of the family of finite elements are proved for both pure displacement and traction problems. Numerical tests in 2D and 3D are presented to show a superiority of the new element over others, as a superconvergence is surprisingly exhibited.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A C^0-Weak Galerkin Finite Element Method for the Biharmonic Equation

A C^0-weak Galerkin (WG) method is introduced and analyzed for solving the biharmonic equation in 2D and 3D. A weak Laplacian is defined for C^0 functions in the new weak formulation. This WG finite element formulation is symmetric, positive definite and parameter free. Optimal order error estimates are established in both a discrete H^2 norm and the L^2 norm, for the weak Galerkin finite element solution. Numerical results are presented to confirm the theory. As a technical tool, a refined Scott-Zhang interpolation operator is constructed to assist the corresponding error estimate. This refined interpolation preserves the volume mass of order (k+1-d) and the surface mass of order (k+2-d) for the P_{k+2} finite element functions in d-dimensional space.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flux-conserving finite element methods

We analyze the flux conservation property of the finite element method. It is shown that the finite element solution does approximate the flux locally in the optimal order, i.e., the same order as that of the nodal interpolation operator. We propose two methods, post-processing the finite element solutions locally. The new solutions, remaining as optimal-order solutions, are flux-conserving elementwise. In one of our methods, the processed solution also satisfies the original finite element equations. While the high-order finite volume schemes are still under construction, our methods produce finite-volume-like finite element solution of any order. In particular, our methods avoid solving non-symmetric finite volume equations. Numerical tests in 2D and 3D verify our findings.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Superconvergence of the $Q_{k+1,k}$-$Q_{k,k+1}$ divergence-free finite element

By the standard theory, the stable $Q_{k+1,k}$-$Q_{k,k+1}/Q_{k}^{dc&#39;}$ divergence-free element converges with the optimal order of approximation for the Stokes equations, but only order $k$ for the velocity in $H^1$-norm and the pressure in $L^2$-norm. This is due to one polynomial degree less in $y$ direction for the first component of velocity, a $Q_{k+1,k}$ polynomial. In this manuscript, we will show a superconvergence of the divergence free element that the order of convergence is truly $k+1$, for both velocity and pressure. Numerical tests are provided confirming the sharpness of the theory.

0 citations0 reviewsNo code linkedOpen access