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

Weifeng Qiu

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.AP

Catalog footprint

What is connected

21works
3topics
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

21 published item(s)

preprint2022arXiv

A $C^{0}$ interior penalty method for $m$th-Laplace equation

In this paper, we propose a $C^{0}$ interior penalty method for $m$th-Laplace equation on bounded Lipschitz polyhedral domain in $\mathbb{R}^{d}$, where $m$ and $d$ can be any positive integers. The standard $H^{1}$-conforming piecewise $r$-th order polynomial space is used to approximate the exact solution $u$, where $r$ can be any integer greater than or equal to $m$. Unlike the interior penalty method in [T.~Gudi and M.~Neilan, {\em An interior penalty method for a sixth-order elliptic equation}, IMA J. Numer. Anal., \textbf{31(4)} (2011), pp. 1734--1753], we avoid computing $D^{m}$ of numerical solution on each element and high order normal derivatives of numerical solution along mesh interfaces. Therefore our method can be easily implemented. After proving discrete $H^{m}$-norm bounded by the natural energy semi-norm associated with our method, we manage to obtain stability and optimal convergence with respect to discrete $H^{m}$-norm. Numerical experiments validate our theoretical estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The pointwise stabilities of piecewise linear finite element method on non-obtuse tetrahedral meshes of nonconvex polyhedra

Let $Ω$ be a Lipschitz polyhedral (can be nonconvex) domain in $\mathbb{R}^{3}$, and $V_{h}$ denotes the finite element space of continuous piecewise linear polynomials. On non-obtuse quasi-uniform tetrahedral meshes, we prove that the finite element projection $R_{h}u$ of $u \in H^{1}(Ω) \cap C(\overlineΩ)$ (with $R_{h} u$ interpolating $u$ at the boundary nodes) satisfies \begin{align*} \Vert R_{h} u\Vert_{L^{\infty}(Ω)} \leq C \vert \log h \vert \Vert u\Vert_{L^{\infty}(Ω)}. \end{align*} If we further assume $u \in W^{1,\infty}(Ω)$, then \begin{align*} \Vert R_{h} u\Vert_{W^{1, \infty}(Ω)} \leq C \vert \log h \vert \Vert u\Vert_{W^{1, \infty}(Ω)}. \end{align*}

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A mixed finite element scheme for biharmonic equation with variable coefficient and von Kármán equations

In this paper, a new mixed finite element scheme using element-wise stabilization is introduced for the biharmonic equation with variable coefficient on Lipschitz polyhedral domains. The proposed scheme doesn't involve any integration along mesh interfaces. The gradient of the solution is approximated by $H({\rm div})$-conforming $BDM_{k+1}$ element or vector valued Lagrange element with order $k+1$, while the solution is approximated by Lagrange element with order $k+2$ for any $k\geq 0$.This scheme can be easily implemented and produces positive definite linear system. We provide a new discrete $H^{2}$-norm stability, which is useful not only in analysis of this scheme but also in $C^{0}$ interior penalty methods and DG methods. Optimal convergences in both discrete $H^{2}$-norm and $L^{2}$-norm are derived. This scheme with its analysis is further generalized to the von Kármán equations. Finally, numerical results verifying the theoretical estimates of the proposed algorithms are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of a semi-implicit structure-preserving finite element method for the nonstationary incompressible Magnetohydrodynamics equations

We revise the structure-preserving finite element method in [K. Hu, Y. MA and J. Xu. (2017) Stable finite element methods preserving $\nabla \cdot \mathbf{B}=0$ exactly for MHD models. Numer. Math.,135, 371-396]. The revised method is semi-implicit in time-discretization. We prove the linearized scheme preserves the divergence free property for the magnetic field exactly at each time step. Further, we showed the linearized scheme is unconditionally stable and we obtain optimal convergence in the energy norm of the revised method even for solutions with low regularity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of an interior penalty DG method for the quad-curl problem

The quad-curl term is an essential part of the resistive magnetohydrodynamic (MHD) equation and the fourth order inverse electromagnetic scattering problem, which are both of great significance in science and engineering. It is desirable to develop efficient and practical numerical methods for the quad-curl problem. In this paper, we first present some new regularity results for the quad-curl problem on Lipschitz polyhedron domains and then propose a mixed finite element method for solving the quad-curl problem. With a {\em novel} discrete Sobolev imbedding inequality for the piecewise polynomials, we obtain stability results and derive error estimates based on a relatively low regularity assumption of the exact solution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A note on the Monge-Ampère type equations with general source terms

In this paper we consider the generalised solutions to the Monge-Amp{è}re type equations with general source terms. We firstly prove the so-called comparison principle and then give some important propositions for the border of generalised solutions. Furthermore, we design well-posed finite element methods for the generalised solutions with the classical and weak Dirichlet boundary conditions respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A superconvergent HDG method for the Maxwell equations

We present and analyze a new hybridizable discontinuous Galerkin (HDG) method for the steady state Maxwell equations. In order to make the problem well-posed, a condition of divergence is imposed on the electric field. Then a Lagrange multiplier $p$ is introduced, and the problem becomes the solution of a mixed curl-curl formulation of the Maxwell's problem. We use polynomials of degree $k+1$, $k$, $k$ to approximate $\bfu,\nabla \times \bfu$ and $p$ respectively. In contrast, we only use a non-trivial subspace of polynomials of degree $k+1$ to approximate the numerical tangential trace of the electric field and polynomials of degree $k+1$ to approximate the numerical trace of the Lagrange multiplier on the faces. On the simplicial meshes, a special choice of the stabilization parameters is applied, and the HDG system is shown to be well-posed. Moreover, we show that the convergence rates for $\boldsymbol{u}$ and $\nabla \times \boldsymbol{u}$ are independent of the Lagrange multiplier $p$. If we assume the dual operator of the Maxwell equation on the domain has adequate regularity, we show that the convergence rate for $\boldsymbol{u}$ is $O(h^{k+2})$. From the point of view of degrees of freedom of the globally coupled unknown: numerical trace, this HDG method achieves superconvergence for the electric field without postprocessing. Finally, we show that on general polyhedral elements, by a particular choice of the stabilization parameters again, the HDG system is also well-posed and the superconvergence of the HDG method is derived.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An absolutely stable $hp$-HDG method for the time-harmonic Maxwell equations with high wave number

We present and analyze a hybridizable discontinuous Galerkin (HDG) method for the time-harmonic Maxwell equations. The divergence-free condition is enforced on the electric field, then a Lagrange multiplier is introduced, and the problem becomes the solution of a mixed curl-curl formulation of the Maxwell's problem. The method is shown to be an absolutely stable HDG method for the indefinite time-harmonic Maxwell equations with high wave number. By exploiting the duality argument, the dependence of convergence of the HDG method on the wave number k, the mesh size h and the polynomial order p is obtained. Numerical results are given to verify the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An HDG method for linear elasticity with strong symmetric stresses

This paper presents a new hybridizable discontinuous Galerkin (HDG) method for linear elasticity on general polyhedral meshes, based on a strong symmetric stress formulation. The key feature of this new HDG method is the use of a special form of the numerical trace of the stresses, which makes the error analysis different from the projection-based error analyzes used for most other HDG methods. For arbitrary polyhedral elements, we approximate the stress by using polynomials of degree k>=1 and the displacement by using polynomials of degree k+1. In contrast, to approximate the numerical trace of the displacement on the faces, we use polynomials of degree k only. This allows for a very efficient implementation of the method, since the numerical trace of the displacement is the only globally-coupled unknown, but does not degrade the convergence properties of the method. Indeed, we prove optimal orders of convergence for both the stresses and displacements on the elements. In the almost incompressible case, we show the error of the stress is also optimal in the standard L2-norm. These optimal results are possible thanks to a special superconvergence property of the numerical traces of the displacement, and thanks to the use of a crucial elementwise Korn's inequality. Several numerical results are presented to support our theoretical findings in the end.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A first order system least squares method for the Helmholtz equation

We present a first order system least squares (FOSLS) method for the Helmholtz equation at high wave number k, which always deduces Hermitian positive definite algebraic system. By utilizing a non-trivial solution decomposition to the dual FOSLS problem which is quite different from that of standard finite element method, we give error analysis to the hp-version of the FOSLS method where the dependence on the mesh size h, the approximation order p, and the wave number k is given explicitly. In particular, under some assumption of the boundary of the domain, the L2 norm error estimate of the scalar solution from the FOSLS method is shown to be quasi optimal under the condition that kh/p is sufficiently small and the polynomial degree p is at least O(\log k). Numerical experiments are given to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A high order HDG method for curved-interface problems via approximations from straight triangulations

We generalize the technique of [Solving Dirichlet boundary-value problems on curved domains by extensions from subdomains, SIAM J. Sci. Comput. 34, pp. A497--A519 (2012)] to elliptic problems with mixed boundary conditions and elliptic interface problems involving a non-polygonal interface. We study first the treatment of the Neumann boundary data since it is crucial to understand the applicability of the technique to curved interfaces. We provide numerical results showing that, in order to obtain optimal high order convergence, it is desirable to construct the computational domain by interpolating the boundary/interface using piecewise linear segments. In this case the distance of the computational domain to the exact boundary is only $O(h^2)$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A superconvergent HDG method for the Incompressible Navier-Stokes Equations on general polyhedral meshes

We present a superconvergent hybridizable discontinuous Galerkin (HDG) method for the steady-state incompressible Navier-Stokes equations on general polyhedral meshes. For arbitrary conforming polyhedral mesh, we use polynomials of degree k+1, k, k to approximate the velocity, velocity gradient and pressure, respectively. In contrast, we only use polynomials of degree k to approximate the numerical trace of the velocity on the interfaces. Since the numerical trace of the velocity field is the only globally coupled unknown, this scheme allows a very efficient implementation of the method. For the stationary case, and under the usual smallness condition for the source term, we prove that the method is well defined and that the global L2-norm of the error in each of the above-mentioned variables and the discrete H1-norm of the error in the velocity converge with the order of k+1 for k>=0. We also show that for k>=1, the global L2-norm of the error in velocity converges with the order of k+2. From the point of view of degrees of freedom of the globally coupled unknown: numerical trace, this method achieves optimal convergence for all the above-mentioned variables in L2-norm for k>=0, superconvergence for the velocity in the discrete H1-norm without postprocessing for k>=0, and superconvergence for the velocity in L2-norm without postprocessing for k>=1.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Direct computation of stresses in linear elasticity

We present a new finite element method based on the formulation introduced by Philippe G.~Ciarlet and Patrick Ciarlet, Jr. in [{\em Math. Models Methods Appl. Sci., 15 (2005), pp. 259--571}], which approximates strain tensor directly. We also show the convergence rate of strain tensor is optimal. This work is a non-trivial generalization of its two dimensional analogue in [{\em Math. Models Methods Appl. Sci., 19 (2009), pp. 1043--1064}]

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An analysis of HDG methods for convection dominated diffusion problems

In this paper, we establish an a priori error analysis of HDG methods with two types of stabilization parameter applied to convection dominated diffusion problem. We show that, using polynomials of degree no greater than k, L2 error of the scalar variable of HDG methods converges in order k + 1/2 if the meshes are not aligned with the flux. This result implies that the HDG approximation converges as fast as conventional DG methods. In the case that the meshes are aligned with the flux, we show that HDG methods with the second type of stabilization parameter achieves optimal rate of convergence. We also discuss preferred form for implementation of HDG methods. In appendix, we show that the spectral condition number of stiffness matrix of HDG methods is independent of the diffusion coefficient. Numerical experiments are presented to verify our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First order least squares method with weakly imposed boundary condition for convection dominated diffusion problems

We present and analyze a first order least squares method for convection dominated diffusion problems, which provides robust L2 a priori error estimate for the scalar variable even if the given data f in L2 space. The novel theoretical approach is to rewrite the method in the framework of discontinuous Petrov - Galerkin (DPG) method, and then show numerical stability by using a key equation discovered by J. Gopalakrishnan and W. Qiu [Math. Comp. 83(2014), pp. 537-552]. This new approach gives an alternative way to do numerical analysis for least squares methods for a large class of differential equations. We also show that the condition number of the global matrix is independent of the diffusion coefficient. A key feature of the method is that there is no stabilization parameter chosen empirically. In addition, Dirichlet boundary condition is weakly imposed. Numerical experiments verify our theoretical results and, in particular, show our way of weakly imposing Dirichlet boundary condition is essential to the design of least squares methods - numerical solutions on subdomains away from interior layers or boundary layers have remarkable accuracy even on coarse meshes, which are unstructured quasi-uniform.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mixed finite elements for elasticity on quadrilateral meshes

We present stable mixed finite elements for planar linear elasticity on general quadrilateral meshes. The symmetry of the stress tensor is imposed weakly and so there are three primary variables, the stress tensor, the displacement vector field, and the scalar rotation. We develop and analyze a stable family of methods, indexed by an integer $r \geq 2$ and with rate of convergence in the $L^2$ norm of order $r$ for all the variables. The methods use Raviart-Thomas elements for the stress, piecewise tensor product polynomials for the displacement, and piecewise polynomials for the rotation. We also present a simple first order element, not belonging to this family. It uses the lowest order BDM elements for the stress, and piecewise constants for the displacement and rotation, and achieves first order convergence for all three variables.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Robust a Posteriori Error Estimates for HDG method for Convection-Diffusion Equations

We propose a robust a posteriori error estimator for the hybridizable discontinuous Galerkin (HDG) method for convection-diffusion equations with dominant convection. The reliability and efficiency of the estimator are established for the error measured in an energy norm. The energy norm is uniformly bounded even when the diffusion coefficient tends to zero. The estimators are robust in the sense that the upper and lower bounds of error are uniformly bounded with respect to the diffusion coefficient. A weighted test function technique and the Oswald interpolation are key ingredients in the analysis. Numerical results verify the robustness of the proposed a posteriori error estimator. In numerical experiments, optimal convergence is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A locking free hp DPG method for linear elasticity with symmetric stresses

We present two new methods for linear elasticity with simultaneously yield stress and displacement approximations of optimal accuracy in both the mesh size h and polynomial degree p. This is achieved within the recently developed discontinuous Petrov-Galerkin (DPG) framework. In this framework, both the stress and the displacement approximations are discontinuous across element interfaces. We study locking-free convergence properties and the interrelationships between the two DPG methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An analysis of the practical DPG method

In this work we give a complete error analysis of the Discontinuous Petrov Galerkin (DPG) method, accounting for all the approximations made in its practical implementation. Specifically, we consider the DPG method that uses a trial space consisting of polynomials of degree $p$ on each mesh element. Earlier works showed that there is a "trial-to-test" operator $T$, which when applied to the trial space, defines a test space that guarantees stability. In DPG formulations, this operator $T$ is local: it can be applied element-by-element. However, an infinite dimensional problem on each mesh element needed to be solved to apply $T$. In practical computations, $T$ is approximated using polynomials of some degree $r > p$ on each mesh element. We show that this approximation maintains optimal convergence rates, provided that $r\ge p+N$, where $N$ is the space dimension (two or more), for the Laplace equation. We also prove a similar result for the DPG method for linear elasticity. Remarks on the conditioning of the stiffness matrix in DPG methods are also included.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Partial expansion of a Lipschitz domain and some applications

We show that a Lipschitz domain can be expanded solely near a part of its boundary, assuming that the part is enclosed by a piecewise C1 curve. The expanded domain as well as the extended part are both Lipschitz. We apply this result to prove a regular decomposition of standard vector Sobolev spaces with vanishing traces only on part of the boundary. Another application in the construction of low-regularity projectors into finite element spaces with partial boundary conditions is also indicated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Variable Order Mixed H-Finite Element Method for Linear Elasticity with Weakly Imposed Symmetry. Ii. Affine and Curvilinear Elements in 2D

We continue our study on variable order Arnold-Falk-Winther elements for 2D elasticity in context of both affine and parametric curvilinear elements. We present an $h$-stability result for affine elements, and an asymptotic stability result for curvilinear elements. Both theoretical results are confirmed with numerical experiments.

0 citations0 reviewsNo code linkedOpen access