Superconvergence Properties Research Articles

Hybridized formulations of flux reconstruction schemes for advection-diffusion problems

We present the hybridization of flux reconstruction methods for advection-diffusion problems. Hybridization introduces a new variable into the problem so that it can be reduced via static condensation. This allows the solution of implicit discretizations to be done more efficiently. We derive an energy statement from a stability analysis considering a range of correction functions on hybridized and embedded flux reconstruction schemes. Then, we establish connections to standard formulations. We devise a post-processing scheme that leverages existing flux reconstruction operators to enhance accuracy for diffusion-dominated problems. Results show that the implicit convergence of these methods for advection-diffusion problems can result in performance benefits of over an order of magnitude. In addition, we observe that the superconvergence property of hybridized methods can be extended to the family of FR schemes for a range of correction functions.

Superconvergence analysis of an energy-conserving Galerkin FEM for the two-dimensional sine-Gordon equation

In this paper, the superconvergence error analysis of an energy-conserving Galerkin fully discrete scheme is proposed and investigated for the two-dimensional sine-Gordon equation. The unique solvability of the numerical scheme as well as the energy conservation are studied firstly. Then, based on the special property of the bilinear element on the rectangular mesh and the superclose estimate between interpolation and Ritz projection of the exact solution in H 1 -norm, the global superconvergence result in H 1 -norm is obtained in terms of a post-processing technique. Finally, numerical results are provided to confirm the energy conservation and superconvergence properties.

Error estimates of the direct discontinuous Galerkin method for two-dimensional nonlinear convection–diffusion equations: Superconvergence analysis

This paper is concerned with superconvergence properties of the direct discontinuous Galerkin (DDG) method for two-dimensional nonlinear convection–diffusion equations. To copy with the complexity and difficulty caused by the nonlinear term, we utilize the method of Taylor expansion and correction idea to achieve our superconvergence goal. Specifically, we first adopt Taylor expansion to decompose the error into two parts: a high-order nonlinear term and a lower-order linear term. Then by using the idea of correction function for the linear term, a high order error bound is obtained. By doing so, a unified analysis of DDG method for general nonlinear problems is finally established. We prove that, for any piecewise tensor-product polynomials of degree k≥2, the DDG solution is superconvergent at nodes and Lobatto points, with an order of O(h2k) and O(hk+2), respectively. Moreover, superconvergence properties for the derivative approximation are also studied and the superconvergence points are identified at Gauss points, with an order of O(hk+1). Numerical experiments are presented to confirm the sharpness of all the theoretical findings.

Error Approximation of the Second Order Hyperbolic Differential Equationby Using DG Finite Element Method

This article presents a simple efficient and asynchronously correcting a posteriori error approximation for discontinuous finite element solutions of the second-order hyperbolic partial differential problems on triangular meshes. This study considersthe basis functions for error spaces corresponding to some finite element spaces. The discretization error of each triangle is estimated by solving the local error problem. It also shows global super convergence for discontinuous solution on triangular lattice. In this article, the triangular elements are classify into three types: (i) elements with one inflow and two outflow edges are of type I, (ii) elements with two inflows and one outflow edges are of type II and (iii) elements with one inflow edge, one outflow edge, and one edge parallel to the characteristics are of type III. The article investigated higher-dimension discontinuous Galerkin methods for hyperbolic problems on triangular meshes and also studied the effect of finite element spaces on the superconvergence properties of DG solutions on three types of triangular elements and it showed that the DG solution is <i>O(h<sup>p+2</sup>)</i> superconvergent at Legendre points on the outflow edge on triangles having one outflow edge using three polynomial spaces. A posteriori error estimates are tested on a number of linear and nonlinear problems to show their efficiency and accuracy under lattice refinement for smooth and discontinuous solutions.

Superconvergence analysis of the nonconforming FEM for the Allen–Cahn equation with time Caputo–Hadamard derivative

This paper first introduces a high-order approximate formula for the Caputo–Hadamard fractional derivative using interpolation theory, and analyzes its convergence and accuracy. Although the formula bears resemblance to the logarithmic L2-1σ formula in form, it differs in the selection of mesh points, making it suitable for more complex models in numerical computations. Subsequently, a numerical method for the Allen–Cahn equation with time Caputo–Hadamard derivative is proposed. The temporal direction is discretized using the newly derived difference formula, while the spatial direction is approximated by the anisotropic nonconforming quasi-Wilson finite element method (FEM). It is proven that the derived scheme has optimal convergence accuracy in the L2-norm and global superconvergence property in the H1-norm. Numerical examples are provided to corroborate the theoretical claims.

Study of numerical treatment of functional first-kind Volterra integral equations

<abstract><p>First-kind Volterra integral equations have ill-posed nature in comparison to the second-kind of these equations such that a measure of ill-posedness can be described by $ \nu $-smoothing of the integral operator. A comprehensive study of the convergence and super-convergence properties of the piecewise polynomial collocation method for the second-kind Volterra integral equations (VIEs) with constant delay has been given in <sup>[<xref ref-type="bibr" rid="b1">1</xref>]</sup>. However, convergence analysis of the collocation method for first-kind delay VIEs appears to be a research problem. Here, we investigated the convergence of the collocation solution as a research problem for a first-kind VIE with constant delay. Three test problems have been fairly well-studied for the sake of verifying theoretical achievements in practice.</p></abstract>

Error estimate and superconvergence of a high-accuracy difference scheme for 2D heat equation with nonlocal boundary conditions

<p>In this work, we initially construct an implicit Euler difference scheme for a two-dimensional heat problem, incorporating both local and nonlocal boundary conditions. Subsequently, we harness the power of the discrete Fourier transform and develop an innovative transformation technique to rigorously demonstrate that our scheme attains the asymptotic optimal error estimate in the maximum norm. Furthermore, we derive a series of approximation formulas for the partial derivatives of the solution along the two spatial dimensions, meticulously proving that each of these formulations possesses superconvergence properties. Lastly, to validate our theoretical findings, we present two comprehensive numerical experiments, showcasing the efficiency and accuracy of our approach.</p>

EXtended Hybridizable Discontinuous Galerkin (X-HDG) method for linear convection-diffusion equations on unfitted domains

In this work, we propose a novel strategy for the numerical solution of linear convection diffusion equation (CDE) over unfitted domains. In the proposed numerical scheme, strategies from high order Hybridized Discontinuous Galerkin method and eXtended Finite Element method are combined with the level set definition of the boundaries. The proposed scheme and hence, is named as eXtended Hybridizable Discontinuous Galerkin (XHDG) method. In this regard, the Hybridizable Discontinuous Galerkin (HDG) method is eXtended to the unfitted domains; i.e., the computational mesh does not need to fit to the domain boundary; instead, the boundary is defined by a level set function and cuts through the background mesh arbitrarily. The original unknown structure of HDG and its hybrid nature ensuring the local conservation of fluxes is kept, while developing a modified bilinear form for the elements cut by the boundary. At every cut element, an auxiliary nodal trace variable on the boundary is introduced, which is eliminated afterwards while imposing the boundary conditions. Both stationary and time dependent CDEs are studied over a range of flow regimes from diffusion to convection dominated; using high order (p≤4) XHDG through benchmark numerical examples over arbitrary unfitted domains. Results proved that XHDG inherits optimal (p+1) and super (p+2) convergence properties of HDG while removing the fitting mesh restriction.

Error estimates with low-order polynomial dependence for the fully-discrete finite element invariant energy quadratization scheme of the Allen–Cahn equation

In this paper, for the Allen–Cahn equation, we obtain the error estimate of the temporal semi-discrete scheme, and the fully-discrete finite element numerical scheme, both of which are based on the invariant energy quadratization (IEQ) time-marching strategy. We establish the relationship between the [Formula: see text]-error bound and the [Formula: see text]-stabilities of the numerical solution. Then, by converting the numerical schemes to a form compatible with the original format of the Allen–Cahn equation, using mathematical induction, the superconvergence property of nonlinear terms, and the spectrum argument, the optimal error estimates that only depends on the low-order polynomial degree of [Formula: see text] instead of [Formula: see text] for both of the semi and fully-discrete schemes are derived. Numerical experiment also validates our theoretical convergence analysis.

Analysis of the MAC scheme for the three dimensional Stokes problem

In this paper, we propose a marker and cell method (MAC) with a proper quadrature scheme over non-uniform cuboid grids by using a staggered finite volume element method (FVEM) for 3D Stokes equations. The stability of the proposed MAC scheme is proven under Petrov-Galerkin method. By establishing a connection between MAC and FVEM, we rigorously prove the superconvergence property and the optimal order L2 error estimate. Finally, numerical results verify our theoretical findings.

Convergence of numerical schemes for convection–diffusion–reaction equations on generic meshes

Using the gradient discretisation method (GDM), we provide a generic convergence of numerical approximations for convection diffusion model with non linear reaction term. The convergence of numerical schemes is proved based on a limited number of natural assumptions and properties. The validity of the proposed gradient scheme is examined using the finite volume method performed on different general meshes. We note that the numerical results indicate that a super-convergence property can be expected for the finite volume methods.

Nonconforming modified Quasi-Wilson finite element method for convection–diffusion-reaction equation

The nonconforming modified Quasi-Wilson finite element method (FEM) is employed to investigate the unconditional superconvergence behavior of the convection–diffusion-reaction equation for the linearized energy stable Backward-Euler (BE) and the second order Backward differentiation formula (BDF2) fully discrete schemes on arbitrary quadrilateral meshes. Then, the energy stabilities of the discrete solutions for the proposed schemes are demonstrated through the mathematical induction, which is crucial to proving the unique solvabilities with Brouwer fixed point theorem. Based on the above achievements and the two typical characteristics of the element: one is that its consistency error in the broken H1-norm can be estimated as order O(h2) when the exact solution belongs to H3(Ω), which is one order higher than its interpolation error, and the other is that the nonconforming part of the function in the finite element space is orthogonal to the bilinear polynomials on each element of the subdivision in the integration sense, the superclose estimates without any restriction between spatial partition parameter h and the time step τ are obtained. Furthermore, the superconvergence properties for the aforementioned schemes are derived by the interpolation post-processing technique. Finally, the theoretical analysis is justified by the provided numerical experiments.

Mixed finite elements based on superconvergent patch recovery for strain gradient theory

Strain gradient theory is one of the most powerful candidates to simulate micro/nano-sized structures considering size effect. However, too many degrees of freedom (DOFs) have been placed to satisfy the C1-continuity condition for the development of finite elements based on the strain gradient theory, greatly reducing the practicality of these elements. In this work, four-node quadrilateral and eight-node hexahedral elements are developed to simulate size effect based on strain gradient theory. Since only displacement DOFs are placed on each node without using displacement gradient DOFs, there is dramatic reduction in computational cost. The displacements are approximated with shape functions of standard finite elements. The displacement gradient field is approximated with linear polynomials and the coefficients are determined using the superconvergent patch recovery. It was found that the superconvergent properties of Barlow points are still valid even though the gradient recovery is performed during the calculation of the element stiffness matrix. The preconditioned conjugate gradient method is utilized to solve the system equations. Numerical examples demonstrate the accuracy and efficiency of the proposed elements. Size effect observed in experiments is captured by using the developed finite elements.

Superconvergence and postprocessing of the continuous Galerkin method for nonlinear Volterra integro-differential equations

We propose a novel postprocessing technique for improving the global accuracy of the continuous Galerkin (CG) method for nonlinear Volterra integro-differential equations. The key idea behind the postprocessing technique is to add a higher order Lobatto polynomial of degree k + 1 to the CG approximation of degree k. We first show that the CG method superconverges at the nodal points of the time partition. We further prove that the postprocessed CG approximation converges one order faster than the unprocessed CG approximation in the L2-, H1- and L∞-norms. As a by-product of the postprocessed superconvergence results, we construct several a posteriori error estimators and prove that they are asymptotically exact. Numerical examples are presented to highlight the superconvergence properties of the postprocessed CG approximations and the robustness of the a posteriori error estimators.

Superconvergence error analysis of discontinuous Galerkin method with interior penalties for 2D elliptic convection – diffusion – reaction problems

This article focuses on constructing and analyzing a non-symmetric interior penalty Galerkin method (NIPG) on Shishkin mesh for solving singularly perturbed 2D elliptic boundary-value problems (BVPs). We use piecewise Lagrange interpolation at Gaussian points to improve the order of convergence of the interpolation error. We then study the superconvergence properties of the NIPG method and prove order of convergence in the discrete energy–norm. Various numerical experiments are provided to validate the theoretical results.

Superconvergence analysis of a conservative mixed finite element method for the nonlinear Klein–Gordon–Schrödinger equations

AbstractIn this paper, a linearized mass and energy conservative mixed finite element method (MFEM) is proposed for solving the nonlinear Klein–Gordon–Schrödinger equations. Optimal error estimates without grid‐ratio condition are derived by some rigorous analysis and an error splitting technique, that is, one is the temporal error which is only τ‐dependent and the other is the spatial error which is only h‐dependent. Furthermore, the superclose results are obtained by using the idea of combination of interpolation and projection. Besides, a so‐called “lifting” approach also play an important role to obtain the superclose results. With the above achievements, the global superconvergent properties are deduced through the interpolated post processing operators. Finally, three numerical examples are given to validate the convergence order, unconditional stability, mass, and energy conservation.

Superconvergence of Legendre-Gauss-Lobatto interpolation and spectral collocation

In this paper, we mainly study the superconvergence properties in approximation of the second- and third-order BVPs by spectral collocation methods. The theoretical analyses identify the superconvergence points of the interpolation function. Ample numerical experiments are carried out which perfectly match the theoretical results. In addition, numerical results show that the superconvergence properties also hold in solving PDEs.

A Review of Collocation Approximations to Solutions of Differential Equations

This review considers piecewise polynomial functions, that have long been known to be a useful and versatile tool in numerical analysis, for solving problems which have solutions with irregular features, such as steep gradients and oscillatory behaviour. Examples of piecewise polynomial functions used include splines, in particular B-splines, and Hermite functions. Spline functions are useful for obtaining global approximations whilst Hermite functions are useful for approximation over finite elements. Our aim in this review is to study quintic Hermite functions and develop a numerical collocation scheme for solving ODEs and PDEs. This choice of basis functions is further motivated by the fact that we are interested in solving problems having solutions with steep gradients and oscillatory properties, for which this approximation basis seems to be a suitable choice. We derive the quintic Hermite basis and use it to formulate the orthogonal collocation on finite element (OCFE) method. We present the error analysis for third order ODEs and derive both global and nodal error bounds to illustrate the super-convergence property at the nodes. Numerical simulations using the Julia programming language are performed for both ODEs and PDEs and enhance the theoretical results.

A Hessian recovery-based finite difference method for biharmonic problems

In this paper, we study a new finite difference method by combining Hessian recovery techniques and the ghost points method for biharmonic equations. Numerical results validate that our proposed method has optimal convergence orders in the L2 norm and H1 seminorm. Moreover, superconvergence properties of the recovered gradient and Hessian have been observed.

On mixed finite element approximations of shape gradients in shape optimization with the Navier–Stokes equation

AbstractFor shape optimization of fluid flows governed by the Navier–Stokes equation, we investigate effectiveness of shape gradient algorithms by analyzing convergence and accuracy of mixed finite element approximations to both the distributed and boundary types of shape gradients. We present convergence analysis with a priori error estimates for the two approximate shape gradients. The theoretical analysis shows that the distributed formulation has superconvergence property. Numerical results with comparisons are presented to verify theory and show that the shape gradient algorithm based on the distributed formulation is highly effective and robust for shape optimization.