Search

This paper presents a novel approach to the construction of the lowest order $H(\mathrm{curl})$ and $H(\mathrm{div})$ exponentially-fitted finite element spaces ${\mathcal{S}_{1^-}^{k}}~(k=1,2)$ on 3D simplicial mesh for corresponding convection-diffusion problems. It is noteworthy that this method not only facilitates the construction of the functions themselves but also provides corresponding discrete fluxes simultaneously. Utilizing this approach, we successfully establish a discrete convection-diffusion complex and employ a specialized weighted interpolation to establish a bridge between the continuous complex and the discrete complex, resulting in a coherent framework. Furthermore, we demonstrate the commutativity of the framework when the convection field is locally constant, along with the exactness of the discrete convection-diffusion complex. Consequently, these types of spaces can be directly employed to devise the corresponding discrete scheme through a Petrov-Galerkin method.

This paper introduces an extension of the Morley element for approximating solutions to biharmonic equations. Traditionally limited to piecewise quadratic polynomials on triangular elements, the extension leverages weak Galerkin finite element methods to accommodate higher degrees of polynomials and the flexibility of general polytopal elements. By utilizing the Schur complement of the weak Galerkin method, the extension allows for fewest local degrees of freedom while maintaining sufficient accuracy and stability for the numerical solutions. The numerical scheme incorporates locally constructed weak tangential derivatives and weak second order partial derivatives, resulting in an accurate approximation of the biharmonic equation. Optimal order error estimates in both a discrete $H^2$ norm and the usual $L^2$ norm are established to assess the accuracy of the numerical approximation. Additionally, numerical results are presented to validate the developed theory and demonstrate the effectiveness of the proposed extension., Comment: 22 pages, 4 tables, 4 figures

The generalized weak Galerkin (gWG) finite element method is proposed and analyzed for the biharmonic equation. A new generalized discrete weak second order partial derivative is introduced in the gWG scheme to allow arbitrary combinations of piecewise polynomial functions defined in the interior and on the boundary of general polygonal or polyhedral elements. The error estimates are established for the numerical approximation in a discrete H^2 norm and a L^2 norm. The numerical results are reported to demonstrate the accuracy and flexibility of our proposed gWG method for the biharmonic equation., Comment: 27 pages, 7 tables

This paper introduces an extension of the well-known Morley element for the biharmonic equation, extending its application from triangular elements to general polytopal elements using the weak Galerkin finite element methods. By leveraging the Schur complement of the weak Galerkin method, this extension not only preserves the same degrees of freedom as the Morley element on triangular elements but also expands its applicability to general polytopal elements. The numerical scheme is devised by locally constructing weak tangential derivatives and weak second-order partial derivatives. Error estimates for the numerical approximation are established in both the energy norm and the $L^2$ norm. A series of numerical experiments are conducted to validate the theoretical developments., Comment: 21 pages, 4 tables, 1 figure

A simplified primal-dual weak Galerkin (S-PDWG) finite element method is designed for the Fokker-Planck type equation with non-smooth diffusion tensor and drift vector. The discrete system resulting from S-PDWG method has significantly fewer degrees of freedom compared with the one resulting from the PDWG method proposed by Wang-Wang \cite{WW-fp-2018}. Furthermore, the condition number of the S-PDWG method is smaller than the PDWG method \cite{WW-fp-2018} due to the introduction of a new stabilizer, which provides a potential for designing fast algorithms. Optimal order error estimates for the S-PDWG approximation are established in the $L^2$ norm. A series of numerical results are demonstrated to validate the effectiveness of the S-PDWG method., Comment: 23 pages, 17 tables

Two types of finite element spaces on triangles are constructed for div-div conforming symmetric tensors. Besides the normal-normal continuity, the stress tensor is continuous at vertices and another trace involving combination of derivatives of stress is identified. Polynomial complex, finite element complex, and Hilbert complex are presented and a commuting diagram between them is given. The constructed div-div conforming elements are exploited to discretize the mixed formulation of the biharmonic equation. Optimal order and superconvergence error analysis is provided. By rotation, finite elements for rot-rot conforming symmetric strain are also obtained., Comment: 21 pages

We consider a surface Stokes problem in stream function formulation on a simply connected oriented surface $\Gamma \subset \mathbb{R}^3$ without boundary. This formulation leads to a coupled system of two second order scalar surface partial differential equations (for the stream function and an auxiliary variable). To this coupled system a trace finite element discretization method is applied. The main topic of the paper is an error analysis of this discretization method, resulting in optimal order discretization error bounds. The analysis applies to the surface finite element method of Dziuk-Elliott, too. We also investigate methods for reconstructing velocity and pressure from the stream function approximation. Results of numerical experiments are included.

In this paper we analyze a class of trace finite element methods (TraceFEM) for the discretization of vector-Laplace equations. A key issue in the finite element discretization of such problems is the treatment of the constraint that the unknown vector field must be tangential to the surface (``tangent condition''). We study three different natural techniques for treating the tangent condition, namely a consistent penalty method, a simpler inconsistent penalty method and a Lagrange multiplier method. A main goal of the paper is to present an analysis that reveals important properties of these three different techniques for treating the tangent constraint. A detailed error analysis is presented that takes the approximation of both the geometry of the surface and the solution of the partial differential equation into account. Error bounds in the energy norm are derived that show how the discretization error depends on relevant parameters such as the degree of the polynomials used for the approximation of the solution, the degree of the polynomials used for the approximation of the level set function that characterizes the surface, the penalty parameter and the degree of the polynomials used for the approximation of Lagrange multiplier.

This paper is devoted to the construction and analysis of the finite element approximations for the $H(D)$ convection-diffusion problems, where $D$ can be chosen as ${\rm grad}$, ${\rm curl}$ or ${\rm div}$ in 3D case. An essential feature of these constructions is to properly average the PDE coefficients on the sub-simplexes. The schemes are of the class of exponential fitting methods that result in special upwind schemes when the diffusion coefficient approaches to zero. Their well-posedness are established for sufficiently small mesh size assuming that the convection-diffusion problems are uniquely solvable. Convergence of first order is derived under minimal smoothness of the solution. Some numerical examples are given to demonstrate the robustness and effectiveness for general convection-diffusion problems.

This short note reports a new derivation of the optimal order of the a priori error estimates for conforming virtual element methods (VEM) on 3D polyhedral meshes based on an error equation. The geometric assumptions, which are necessary for the optimal order of the conforming VEM error estimate in the $H^1$-seminorm, are relaxed for that in a bilinear form-induced energy norm.

Compared to conforming P1 finite elements, nonconforming P1 finite element discretizations are thought to be less sensitive to the appearance of distorted triangulations. E.g., optimal-order discrete $H^1$ norm best approximation error estimates for $H^2$ functions hold for arbitrary triangulations. However, similar estimates for the error of the Galerkin projection for second-order elliptic problems show a dependence on the maximum angle of all triangles in the triangulation. We demonstrate on the example of a special family of distorted triangulations that this dependence is essential, and due to the deterioration of the consistency error. We also provide examples of sequences of triangulations such that the nonconforming P1 Galerkin projections for a Poisson problem with polynomial solution do not converge or converge at arbitrarily slow speed. The results complement analogous findings for conforming P1 elements., Comment: 23 pages, 10 figures

This paper generalizes the non-conforming FEM of Crouzeix and Raviart and its fundamental projection property by a novel mixed formulation for the Poisson problem based on the Helmholtz decomposition. The new formulation allows for ansatz spaces of arbitrary polynomial degree and its discretization coincides with the mentioned non-conforming FEM for the lowest polynomial degree. The discretization directly approximates the gradient of the solution instead of the solution itself. Besides the a priori and medius analysis, this paper proves optimal convergence rates for an adaptive algorithm for the new discretization. These are also demonstrated in numerical experiments. Furthermore, this paper focuses on extensions of this new scheme to quadrilateral meshes, mixed FEMs, and three space dimensions.

We apply the local discontinuous Galerkin (LDG for short) method to solve a mixed boundary value problems for the Helmholtz equation in bounded polygonal domain in 2D. Under some assumptions on regularity of the solution of an adjoint problem, we prove that: (a) the corresponding indefinite discrete scheme is well posed; (b) there is convergence with the expected convergence rates as long as the meshsize h is small enough. We give precise information on how small h has to be in terms of the size of the wavenumber and its distance to the set of eigenvalues for the same boundary value problem for the Laplacian. We also present an a posteriori error estimator showing both the reliability and efficiency of the estimator complemented with detailed information on the dependence of the constants on the wavenumber. We finish presenting extensive numerical experiments which illustrate the theoretical results proven in this paper and suggest that stability and convergence may occur under less restrictive assumptions than those taken in the present work.

We present a simple discretization by radial basis functions for the Poisson equation with Dirichlet boundary condition. A Lagrangian multiplier using piecewise polynomials is used to accommodate the boundary condition. This simplifies previous attempts to use radial basis functions in the interior domain to approximate the solution and on the boundary to approximate the multiplier, which technically requires that the mesh norm in the interior domain is significantly smaller than that on the boundary. Numerical experiments confirm theoretical results.

In this note we show that the non-symmetric version of the classical Nitsche's method for the weak imposition of boundary conditions is stable without penalty term. We prove optimal $H^1$-error estimates and $L^2$-estimates that are suboptimal with half an order in $h$. Both the pure diffusion and the convection--diffusion problems are discussed., Comment: 20 pages, 9 figures

In this paper we develop two conforming finite element methods for a fourth order bi-wave equation arising as a simplified Ginzburg-Landau-type model for d-wave superconductors in absence of applied magnetic field. Unlike the biharmonic operator $\Delta^2$, the bi-wave operator $\Box^2$ is not an elliptic operator, so the energy space for the bi-wave equation is much larger than the energy space for the biharmonic equation. This then makes it possible to construct low order conforming finite elements for the bi-wave equation. However, the existence and construction of such finite elements strongly depends on the mesh. In the paper, we first characterize mesh conditions which allow and not allow construction of low order conforming finite elements for approximating the bi-wave equation. We then construct a cubic and a quartic conforming finite element. It is proved that both elements have the desired approximation properties, and give optimal order error estimates in the energy norm, suboptimal (and optimal in some cases) order error estimates in the $H^1$ and $L^2$ norm. Finally, numerical experiments are presented to guage the efficiency of the proposed finite element methods and to validate the theoretical error bounds., Comment: 21 pages, 19 figures and 1 table