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

Junping Wang

Junping Wang 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-phmath.MP

Catalog footprint

What is connected

32works
4topics
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

32 published item(s)

preprint2022arXiv

A Generalized Weak Galerkin method for Oseen equation

In this work, the authors introduce a generalized weak Galerkin (gWG) finite element method for the time-dependent Oseen equation. The generalized weak Galerkin method is based on a new framework for approximating the gradient operator. Both a semi-discrete and a fully-discrete numerical scheme are developed and analyzed for their convergence, stability, and error estimates. A generalized {\em{inf-sup}} condition is developed to assist the convergence analysis. The backward Euler discretization is employed in the design of the fully-discrete scheme. Error estimates of optimal order are established mathematically, and they are validated numerically with some benchmark examples.

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

An $L^p$- Primal-Dual Weak Galerkin method for div-curl Systems

This paper presents a new $L^p$-primal-dual weak Galerkin (PDWG) finite element method for the div-curl system with the normal boundary condition for $p>1$. Two crucial features for the proposed $L^p$-PDWG finite element scheme are as follows: (1) it offers an accurate and reliable numerical solution to the div-curl system under the low $W^{α, p}$-regularity ($α>0$) assumption for the exact solution; (2) it offers an effective approximation of the normal harmonic vector fields on domains with complex topology. An optimal order error estimate is established in the $L^q$-norm for the primal variable where $\frac{1}{p}+\frac{1}{q}=1$. A series of numerical experiments are presented to demonstrate the performance of the proposed $L^p$-PDWG algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A New Numerical Method for Div-Curl Systems with Low Regularity Assumptions

This paper presents a numerical method for div-curl systems with normal boundary conditions by using a finite element technique known as primal-dual weak Galerkin (PDWG). The PDWG finite element scheme for the div-curl system has two prominent features in that it offers not only an accurate and reliable numerical solution to the div-curl system under the low $H^α$-regularity ($α>0$) assumption for the true solution, but also an effective approximation of normal harmonic vector fields regardless the topology of the domain. Results of seven numerical experiments are presented to demonstrate the performance of the PDWG algorithm, including one example on the computation of discrete normal harmonic vector fields.

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

Finite element method with the total stress variable for Biot's consolidation model

In this work, semi-discrete and fully-discrete error estimates are derived for the Biot's consolidation model described using a three-field finite element formulation. The fields include displacements, total stress and pressure. The model is implemented using a backward Euler discretization in time for the fully-discrete scheme and validated for benchmark examples. Computational experiments presented verifies the convergence orders for the lowest order finite elements with discontinuous and continuous finite element appropriation for the total stress.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A primal-dual finite element method for first-order transport problems

This article devises a new numerical method for first-order transport problems by using the primal-dual weak Galerkin (PD-WG) finite element method recently developed in scientific computing. The PD-WG method is based on a variational formulation of the modeling equation for which the differential operator is applied to the test function so that low regularity for the exact solution of the original equation is sufficient for computation. The PD-WG finite element method indeed yields a symmetric system involving both the original equation for the primal variable and its dual for the dual variable (also known as Lagrangian multiplier). For the linear transport problem, it is shown that the PD-WG method offers numerical solutions that conserve mass locally on each element. Optimal order error estimates in various norms are derived for the numerical solutions arising from the PD-WG method with weak regularity assumptions on the modelling equations. A variety of numerical results are presented to demonstrate the accuracy and stability of the new method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Discrete Divergence-Free Weak Galerkin Finite Element Method for the Stokes Equations

A discrete divergence-free weak Galerkin finite element method is developed for the Stokes equations based on a weak Galerkin (WG) method introduced in the reference [15]. Discrete divergence-free bases are constructed explicitly for the lowest order weak Galerkin elements in two and three dimensional spaces. These basis functions can be derived on general meshes of arbitrary shape of polygons and polyhedrons. With the divergence-free basis derived, the discrete divergence-free WG scheme can eliminate the pressure variable from the system and reduces a saddle point problem to a symmetric and positive definite system with many fewer unknowns. Numerical results are presented to demonstrate the robustness and accuracy of this discrete divergence-free WG method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective Implementation of the Weak Galerkin Finite Element Methods for the Biharmonic Equation

The weak Galerkin (WG) methods have been introduced in the references [11, 16] for solving the biharmonic equation. The purpose of this paper is to develop an algorithm to implement the WG methods effectively. This can be achieved by eliminating local unknowns to obtain a global system with significant reduction of size. In fact, this reduced global system is equivalent to the Schur complements of the WG methods. The unknowns of the Schur complement of the WG method are those defined on the element boundaries. The equivalence of the WG method and its Schur complement is established. The numerical results demonstrate the effectiveness of this new implementation technique.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hybridized Formulation for the Weak Galerkin Mixed Finite Element Method

This paper presents a hybridized formulation for the weak Galerkin mixed finite element method (WG-MFEM) which was introduced and analyzed for second order elliptic equations. The WG-MFEM method was designed by using discontinuous piecewise polynomials on finite element partitions consisting of polygonal or polyhedral elements of arbitrary shape. The key to WG-MFEM is the use of a discrete weak divergence operator which is defined and computed by solving inexpensive problems locally on each element. The hybridized formulation of this paper leads to a significantly reduced system of linear equations involving only the unknowns arising from the Lagrange multiplier in hybridization. Optimal-order error estimates are derived for the hybridized WG-MFEM approximations. Some numerical results are reported to confirm the theory and a superconvergence for the Lagrange multiplier.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Locking-Free Weak Galerkin Finite Element Method for Elasticity Problems in the Primal Formulation

This paper presents an arbitrary order locking-free numerical scheme for linear elasticity on general polygonal/polyhedral partitions by using weak Galerkin (WG) finite element methods. Like other WG methods, the key idea for the linear elasticity is to introduce discrete weak strain and stress tensors which are defined and computed by solving inexpensive local problems on each element. Such local problems are derived from weak formulations of the corresponding differential operators through integration by parts. Locking-free error estimates of optimal order are derived in a discrete $H^1$-norm and the usual $L^2$-norm for the approximate displacement when the exact solution is smooth. Numerical results are presented to demonstrate the efficiency, accuracy, and the locking-free property of the weak Galerkin finite element method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form

This article proposes a new numerical algorithm for second order elliptic equations in non-divergence form. The new method is based on a discrete weak Hessian operator locally constructed by following the weak Galerkin strategy. The numerical solution is characterized as a minimization of a non-negative quadratic functional with constraints that mimic the second order elliptic equation by using the discrete weak Hessian. The resulting Euler-Lagrange equation offers a symmetric finite element scheme involving both the primal and a dual variable known as the Lagrange multiplier, and thus the name of primal-dual weak Galerkin finite element method. Optimal order error estimates are derived for the finite element approximations in a discrete $H^2$-norm, as well as the usual $H^1$- and $L^2$-norms. Some numerical results are presented for smooth and non-smooth coefficients on convex and non-convex domains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discretization of div-curl Systems by Weak Galerkin Finite Element Methods on Polyhedral Partitions

In this paper, the authors devise a new discretization scheme for div-curl systems defined in connected domains with heterogeneous media by using the weak Galerkin finite element method. Two types of boundary value problems are considered in the algorithm development: (1) normal boundary condition, and (2) tangential boundary condition. A new variational formulation is developed for the normal boundary value problem by using the Helmholtz decomposition which avoids the computation of functions in the harmonic fields. Both boundary value problems are reduced to a general saddle-point problem involving the curl and divergence operators, for which the weak Galerkin finite element method is devised and analyzed. The novelty of the technique lies in the discretization of the divergence operator applied to vector fields with heterogeneous media. Error estimates of optimal order are established for the corresponding finite element approximations in various discrete Sobolev norms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Hybridized Weak Galerkin Finite Element Method for the Biharmonic Equation

This paper presents a hybridized formulation for the weak Galerkin finite element method for the biharmonic equation. The hybridized weak Galerkin scheme is based on the use of a Lagrange multiplier defined on the element boundaries. The Lagrange multiplier is verified to provide a numerical approximation for certain derivatives of the exact solution. An optimal order error estimate is established for the numerical approximations arising from the hybridized weak Galerkin finite element method. The paper also derives a computational algorithm (Schur complement) by eliminating all the unknown variables on each element, yielding a significantly reduced system of linear equations for unknowns on the boundary of each element.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An auxiliary space multigrid preconditioner for the weak Galerkin method

In this paper, the authors constructed an auxiliary space multigrid preconditioner for the weak Galerkin finite element method for second-order diffusion equations, discretized on simplicial 2D or 3D meshes. The idea of the auxiliary space multigrid preconditioner is to use an auxiliary space as a "coarse" space in the multigrid algorithm, where the discrete problem in the auxiliary space can be easily solved by an existing solver. In this construction, the authors conveniently use the $H^1$ conforming piecewise linear finite element space as an auxiliary space. The main technical difficulty is to build the connection between the weak Galerkin discrete space and the $H^1$ conforming piecewise linear finite element space. The authors successfully constructed such an auxiliary space multigrid preconditioner for the weak Galerkin method, as well as a reduced system of the weak Galerkin method involving only the degrees of freedom on edges/faces. The preconditioned systems are proved to have condition numbers independent of the mesh size. Numerical experiments are conducted to support the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation

A weak Galerkin (WG) method is introduced and numerically tested for the Helmholtz equation. This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property. At the same time, the WG finite element formulation is symmetric and parameter free. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Challenging problems with high wave numbers are also examined. Our numerical experiments indicate that the weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Stable Numerical Algorithm for the Brinkman Equations by Weak Galerkin Finite Element Methods

This paper presents a stable numerical algorithm for the Brinkman equations by using weak Galerkin (WG) finite element methods. The Brinkman equations can be viewed mathematically as a combination of the Stokes and Darcy equations which model fluid flow in a multi-physics environment, such as flow in complex porous media with a permeability coefficient highly varying in the simulation domain. In such applications, the flow is dominated by Darcy in some regions and by Stokes in others. It is well known that the usual Stokes stable elements do not work well for Darcy flow and vise versa. The challenge of this study is on the design of numerical schemes which are stable for both the Stokes and the Darcy equations. This paper shows that the WG finite element method is capable of meeting this challenge by providing a numerical scheme that is stable and accurate for both Darcy and the Stokes dominated flows. Error estimates of optimal order are established for the corresponding WG finite element solutions. The paper also presents some numerical experiments that demonstrate the robustness, reliability, flexibility and accuracy of the WG method for the Brinkman equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Weak Galerkin Finite Element Method for the Stokes Equations

This paper introduces a weak Galerkin (WG) finite element method for the Stokes equations in the primary velocity-pressure formulation. This WG method is equipped with stable finite elements consisting of usual polynomials of degree $k\ge 1$ for the velocity and polynomials of degree $k-1$ for the pressure, both are discontinuous. The velocity element is enhanced by polynomials of degree $k-1$ on the interface of the finite element partition. All the finite element functions are discontinuous for which the usual gradient and divergence operators are implemented as distributions in properly-defined spaces. Optimal-order error estimates are established for the corresponding numerical approximation in various norms. It must be emphasized that the WG finite element method is designed on finite element partitions consisting of arbitrary shape of polygons or polyhedra which are shape regular.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Weak Galerkin Finite Element Method with Polynomial Reduction

The novel idea of weak Galerkin (WG) 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 leads to different weak Galerkin finite element methods, which makes WG methods highly flexible and efficient in practical computation. This paper explores the possibility of optimal combination of polynomial spaces that minimize the number of unknowns in the numerical scheme, yet without compromising the accuracy of the numerical approximation. For illustrative purpose, the authors use second order elliptic problems to demonstrate the basic idea of polynomial reduction. A new weak Galerkin finite element method is proposed and analyzed. This new finite element scheme features piecewise polynomials of degree $k\ge 1$ on each element plus piecewise polynomials of degree $k-1\ge 0$ on the edge or face of each element. 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. In addition, the paper presents a great deal of numerical experiments to demonstrate the power of the WG method in dealing with finite element partitions consisting of arbitrary polygons in two dimensional spaces or polyhedra in three dimensional spaces. The numerical examples include various finite element partitions such as triangular mesh, quadrilateral mesh, honey comb mesh in 2d and mesh with deformed cubes in 3d. The numerical results show a great promise of the robustness, reliability, flexibility and accuracy of the WG method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Weak Galerkin Mixed Finite Element Method for Second-Order Elliptic Problems

A new weak Galerkin (WG) method is introduced and analyzed for the second order elliptic equation formulated as a system of two first order linear equations. This method, called WG-MFEM, is designed by using discontinuous piecewise polynomials on finite element partitions with arbitrary shape of polygons/polyhedra. The WG-MFEM is capable of providing very accurate numerical approximations for both the primary and flux variables. Allowing the use of discontinuous approximating functions on arbitrary shape of polygons/polyhedra makes the method highly flexible in practical computation. Optimal order error estimates in both discrete $H^1$ and $L^2$ norms are established for the corresponding weak Galerkin mixed finite element solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Efficient Numerical Scheme for the Biharmonic Equation by Weak Galerkin Finite Element Methods on Polygonal or Polyhedral Meshes

This paper presents a new and efficient numerical algorithm for the biharmonic equation by using weak Galerkin (WG) finite element methods. The WG finite element scheme is based on a variational form of the biharmonic equation that is equivalent to the usual $H^2$-semi norm. Weak partial derivatives and their approximations, called discrete weak partial derivatives, are introduced for a class of discontinuous functions defined on a finite element partition of the domain consisting of general polygons or polyhedra. The discrete weak partial derivatives serve as building blocks for the WG finite element method. The resulting matrix from the WG method is symmetric, positive definite, and parameter free. An error estimate of optimal order is derived in an $H^2$-equivalent norm for the WG finite element solutions. Error estimates in the usual $L^2$ norm are established, yielding optimal order of convergence for all the WG finite element algorithms except the one corresponding to the lowest order (i.e., piecewise quadratic elements). Some numerical experiments are presented to illustrate the efficiency and accuracy of the numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Weak Galerkin Finite Element Methods for Parabolic Equations

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal order error estimates in both H^1 and L^2 norms are established. Numerical tests are performed and reported.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Weak Galerkin Finite Element Methods for the Biharmonic Equation on Polytopal Meshes

A new weak Galerkin (WG) finite element method is introduced and analyzed in this paper for the biharmonic equation in its primary form. This method is highly robust and flexible in the element construction by using discontinuous piecewise polynomials on general finite element partitions consisting of polygons or polyhedra of arbitrary shape. The resulting WG finite element formulation is symmetric, positive definite, and parameter-free. Optimal order error estimates in a discrete $H^2$ norm is established for the corresponding WG finite element solutions. Error estimates in the usual $L^2$ norm are also derived, yielding a sub-optimal order of convergence for the lowest order element and an optimal order of convergence for all high order of elements. Numerical results are presented to confirm the theory of convergence under suitable regularity assumptions.

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

A Weak Galerkin Mixed Finite Element Method for Biharmonic Equations

This article introduces and analyzes a weak Galerkin mixed finite element method for solving the biharmonic equation. The weak Galerkin method, first introduced by two of the authors (J. Wang and X. Ye) in an earlier publication for second order elliptic problems, is based on the concept of discrete weak gradients. The method allows the use of completely discrete finite element functions on partitions of arbitrary polygon or polyhedron. In this article, the idea of weak Galerkin method is applied to discretize the Ciarlet-Raviart mixed formulation for the biharmonic equation. In particular, an a priori error estimation is given for the corresponding finite element approximations. The error analysis essentially follows the framework of Babuska, Osborn, and Pitkaranta and uses specially designed mesh-dependent norms. The proof is technically tedious due to the discontinuous nature of the weak Galerkin finite element functions. Some computational results are presented to demonstrate the efficiency of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interior penalty discontinuous Galerkin method on very general polygonal and polyhedral meshes

This paper focuses on interior penalty discontinuous Galerkin methods for second order elliptic equations on very general polygonal or polyhedral meshes. The mesh can be composed of any polygons or polyhedra which satisfies certain shape regularity conditions characterized in a recent paper by two of the authors in [17]. Such general meshes have important application in computational sciences. The usual $H^1$ conforming finite element methods on such meshes are either very complicated or impossible to implement in practical computation. However, the interior penalty discontinuous Galerkin method provides a simple and effective alternative approach which is efficient and robust. This article provides a mathematical foundation for the use of interior penalty discontinuous Galerkin methods in general meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak Galerkin Finite Element Methods on Polytopal Meshes

This paper introduces a new weak Galerkin (WG) finite element method for second order elliptic equations on polytopal meshes. This method, called WG-FEM, is designed by using a discrete weak gradient operator applied to discontinuous piecewise polynomials on finite element partitions of arbitrary polytopes with certain shape regularity. The paper explains how the numerical schemes are designed and why they provide reliable numerical approximations for the underlying partial differential equations. In particular, optimal order error estimates are established for the corresponding WG-FEM approximations in both a discrete $H^1$ norm and the standard $L^2$ norm. Numerical results are presented to demonstrate the robustness, reliability, and accuracy of the WG-FEM. All the results are derived for finite element partitions with polytopes. Allowing the use of discontinuous approximating functions on arbitrary polytopal elements is a highly demanded feature for numerical algorithms in scientific computing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak Galerkin Methods for Second Order Elliptic Interface Problems

Weak Galerkin methods refer to general finite element methods for PDEs in which differential operators are approximated by their weak forms as distributions. Such weak forms give rise to desirable flexibilities in enforcing boundary and interface conditions. A weak Galerkin finite element method (WG-FEM) is developed in this paper for solving elliptic partial differential equations (PDEs) with discontinuous coefficients and interfaces. The paper also presents many numerical tests for validating the WG-FEM for solving second order elliptic interface problems. For such interface problems, the solution possesses a certain singularity due to the nonsmoothness of the interface. A challenge in research is to design high order numerical methods that work well for problems with low regularity in the solution. The best known numerical scheme in the literature is of order one for the solution itself in $L_\infty$ norm. It is demonstrated that the WG-FEM of lowest order is capable of delivering numerical approximations that are of order 1.75 in the usual $L_\infty$ norm for $C^1$ or Lipschitz continuous interfaces associated with a $C^1$ or $H^2$ continuous solutions. Theoretically, it is proved that high order of numerical schemes can be designed by using the WG-FEM with polynomials of high order on each element.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Computational Study of the Weak Galerkin Method for Second-Order Elliptic Equations

The weak Galerkin finite element method is a novel numerical method that was first proposed and analyzed by Wang and Ye for general second order elliptic problems on triangular meshes. The goal of this paper is to conduct a computational investigation for the weak Galerkin method for various model problems with more general finite element partitions. The numerical results confirm the theory established by Wang and Ye. The results also indicate that the weak Galerkin method is efficient, robust, and reliable in scientific computing.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation with Large Wave Numbers

Weak Galerkin (WG) refers to general finite element methods for partial differential equations in which differential operators are approximated by weak forms through the usual integration by parts. In particular, WG methods allow the use of discontinuous finite element functions in the algorithm design. One of such examples was recently introduced by Wang and Ye for solving second order elliptic problems. The goal of this paper is to apply the WG method of Wang and Ye to the Helmholtz equation with high wave numbers. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Our numerical experiments indicate that weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Weak Galerkin Finite Element Method for Second-Order Elliptic Problems

In this paper, authors shall introduce a finite element method by using a weakly defined gradient operator over discontinuous functions with heterogeneous properties. The use of weak gradients and their approximations results in a new concept called {\em discrete weak gradients} which is expected to play important roles in numerical methods for partial differential equations. This article intends to provide a general framework for operating differential operators on functions with heterogeneous properties. As a demonstrative example, the discrete weak gradient operator is employed as a building block to approximate the solution of a model second order elliptic problem, in which the classical gradient operator is replaced by the discrete weak gradient. The resulting numerical approximation is called a weak Galerkin (WG) finite element solution. It can be seen that the weak Galerkin method allows the use of totally discontinuous functions in the finite element procedure. For the second order elliptic problem, an optimal order error estimate in both a discrete $H^1$ and $L^2$ norms are established for the corresponding weak Galerkin finite element solutions. A superconvergence is also observed for the weak Galerkin approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Maximum Principles for P1-Conforming Finite Element Approximations of Quasi-Linear Second Order Elliptic Equations

This paper derives some discrete maximum principles for $P1$-conforming finite element approximations for quasi-linear second order elliptic equations. The results are extensions of the classical maximum principles in the theory of partial differential equations to finite element methods. The mathematical tools are based on the variational approach that was commonly used in the classical PDE theory. The discrete maximum principles are established by assuming a property on the discrete variational form that is of global nature. In particular, the assumption on the variational form is verified when the finite element partition satisfies some angle conditions. For the general quasi-linear elliptic equation, these angle conditions indicate that each triangle or tetrahedron needs to be $\mathcal{O}(h^α)$-acute in the sense that each angle $α_{ij}$ (for triangle) or interior dihedral angle $α_{ij}$ (for tetrahedron) must satisfy $α_{ij}\le π/2-γh^α$ for some $α\ge 0$ and $γ>0$. For the Poisson problem where the differential operator is given by Laplacian, the angle requirement is the same as the existing ones: either all the triangles are non-obtuse or each interior edge is non-negative. It should be pointed out that the analytical tools used in this paper are based on the powerful De Giorgi's iterative method that has played important roles in the theory of partial differential equations. The mathematical analysis itself is of independent interest in the finite element analysis.

0 citations0 reviewsNo code linkedOpen access