Error Estimates In Norm Research Articles

A locally conservative staggered least squares method on polygonal meshes

<abstract><p>In this paper, we propose a novel staggered least squares method for elliptic equations on polygonal meshes. Our new method can be flexibly applied to rough grids and allows hanging nodes, which is of particular interest in practical applications. Moreover, it offers the advantage of not having to deal with inf-sup conditions and yielding positive definite discrete problems. Optimal a priori error estimates in energy norm are derived. In addition, a superconvergent estimates in energy norm are also developed by employing variational error expansion. The main difficulty involved here is to show the $ L^2 $ norm error estimates for the potential variable, where duality argument and the superconvergent estimates are the key ingredients. The single valued flux over the outer boundary of the dual partition enables us to construct a locally conservative flux. Numerical experiments confirm the theoretical findings and the performance of the adaptive mesh refinement guided by the least squares functional estimator are also displayed.</p></abstract>

Modified [formula omitted] interior penalty analysis for fourth order Dirichlet boundary control problem and a posteriori error estimates

We revisit the L2 norm error estimate for the C0-interior penalty analysis of fourth order Dirichlet boundary control problem. The L2 norm error estimate for the optimal control is derived under reduced regularity assumption and this analysis can be carried out on any convex polygonal domain. Moreover, residual based a posteriori error bounds are derived for the optimal control, state and adjoint state variables under minimal regularity assumptions. The error estimator is shown to be reliable and locally efficient. The theoretical findings are illustrated by numerical experiments.

An enriched cut finite element method for Stokes interface equations

In this paper, we consider an enriched cut finite element method (ECFEM) with interface-unfitted meshes for solving Stokes interface equations consisting of two incompressible fluids with different viscosities. By approximating the velocity with the enriched P1 element and the pressure with the P0 element, and stabilizing the Galerkin variational formulation with suitable ghost penalty terms, we propose the new ECFEM and prove that it is well-posed and has the optimal a priori error estimate in the energy norm. All derived results are independent of the interface position. Moreover, compared with other conforming finite element methods with the optimal rate in convergence, the proposed scheme here not only has the minimum degrees of freedom, but also avoids using the derivative of the pressure in the penalty term. The presented numerical examples validate the theoretical predictions.

A priori and a posteriori error analysis for semilinear problems in liquid crystals

In this paper, we develop a unified framework for the a priori and a posteriori error control of different lowest-order finite element methods for approximating the regular solutions of systems of partial differential equations under a set of hypotheses. The systems involve cubic nonlinearities in lower order terms, non-homogeneous Dirichlet boundary conditions, and the results are established under minimal regularity assumptions on the exact solution. The key contributions include (i) results for existence and local uniqueness of the discrete solutions using Newton–Kantorovich theorem, (ii) a priori error estimates in the energy norm, and (iii) a posteriori error estimates that steer the adaptive refinement process. The results are applied to conforming, Nitsche, discontinuous Galerkin, and weakly over penalized symmetric interior penalty schemes for variational models of ferronematics and nematic liquid crystals. The theoretical estimates are corroborated by substantive numerical results.

Richardson extrapolation method for solving the Riesz space fractional diffusion problem

AbstractThe Richardson extrapolation is widely utilized in numerically solving the differential equations since it enjoys both high accuracy and ease of implementation. For the Riesz space fractional diffusion equation, the fractional centered difference operator is employed to approximate the fractional derivative, then the Richardson extrapolation methods for two difference schemes are constructed with the help of the asymptotic expansions of the discrete solutions. Specifically, for the Crank–Nicolson difference scheme, the extrapolation method contains two extrapolation formulae that achieve the fourth order and the sixth order both in the temporal and spatial directions, respectively. The extrapolation method for the compact difference scheme involves one extrapolation formula by which the sixth order can be obtained when the time step size is proportional to the squares of the space step size. The maximum norm error estimates of the extrapolation solutions are proved by the discrete fractional Sobolev embedding inequalities. The extension to the high dimensional and nonlinear cases is also demonstrated. Numerical results verify the theoretical convergence orders and efficiency of our methods.

A weak Galerkin finite element method for 1D semiconductor device simulation models

In this paper, we study a weak Galerkin (WG) finite element method for semiconductor device simulations. We consider the one-dimensional drift–diffusion (DD) and high-field (HF) models, which involves not only first derivative convection terms but also second derivative diffusion terms, as well as a coupled Poisson potential equation. The main difficulties in the analysis include the treatment of the nonlinearity and coupling of the models. The weak Galerkin finite element method adopts piecewise polynomials of degree k for the approximations of electron concentration and electric potential in the interior of elements, and piecewise polynomials of degree k+1 for the discrete weak derivative space. The optimal order error estimates in a discrete H1 norm and the standard L2 norm are derived. Numerical experiments are presented to illustrate our theoretical analysis. Moreover, numerical schemes also work out for the discontinuous diffusion coefficient problems.

Finite element and isogeometric stabilized methods for the advection-diffusion-reaction equation

We develop two new stabilized methods for the steady advection-diffusion-reaction equation, referred to as the Streamline GSC Method and the Directional GSC Method. Both are globally conservative and perform well in numerical studies utilizing linear, quadratic, cubic, and quartic Lagrange finite elements and maximally smooth B-spline elements. For the streamline GSC method we can prove coercivity, convergence, and optimal-order error estimates in a strong norm that are robust in the advective and reactive limits. The directional GSC method is designed to accurately resolve boundary layers for flows that impinge upon the boundary at an angle, a long-standing problem. The directional GSC method performs better than the streamline GSC method in the numerical studies, but it is not coercive. We conjecture it is inf-sup stable but we are unable to prove it at this time. However, calculations of the inf-sup constant support the conjecture. In the numerical studies, B-spline finite elements consistently perform better than Lagrange finite elements of the same order and number of unknowns.

A convergent finite element algorithm for mean curvature flow in arbitrary codimension

Optimal-order uniform-in-time H^1 -norm error estimates are given for semi- and full discretizations of mean curvature flow of surfaces in arbitrarily high codimension. The proposed and studied numerical method is based on a parabolic system coupling the surface flow to evolution equations for the mean curvature vector and for the orthogonal projection onto the tangent space. The algorithm uses evolving surface finite elements and linearly implicit backward difference formulas. This numerical method admits a convergence analysis in the case of finite elements of polynomial degree at least 2 and backward difference formulas of orders 2 to 5. Numerical experiments in codimension 2 illustrate and complement our theoretical results.

An upwind finite volume element method on a Shishkin mesh for singularly perturbed convection–diffusion problems

In this paper, we construct and analyze an upwind finite volume element method on a Shishkin mesh for singularly perturbed convection–diffusion problems. We prove the stability of the method under the assumption of the convection and reaction term coefficients. The error estimate in the energy norm is presented on the Shishkin mesh, and the optimal error bound O(N−1(lnN)3/2) is obtained. This error bound is uniformly valid with respect to the singular perturbation parameter. Numerical examples are provided to illustrate our theoretical results.

A modified weak Galerkin finite element method for parabolic equations on anisotropic meshes

In this paper, we study a modified weak Galerkin finite element method (MWG-FEM) on anisotropic triangular meshes for a class of parabolic equations. Different from conventional weak Galerkin methods, the MWG-FEM replaces the boundary functions by the average of interior functions, which possesses flexibility in the approximation functions and mesh generation. Moreover, MWG-FEM is compatible with anisotropic meshes and suitable for solving problems with anisotropic property. By using the anisotropic interpolation and projection operators, the optimal order error estimate in L2 norm is derived. Numerical experiments are performed to demonstrate the stability and efficiency of the method.

Higher-order error estimates for physics-informed neural networks approximating the primitive equations

Large-scale dynamics of the oceans and the atmosphere are governed by primitive equations (PEs). Due to the nonlinearity and nonlocality, the numerical study of the PEs is generally challenging. Neural networks have been shown to be a promising machine learning tool to tackle this challenge. In this work, we employ physics-informed neural networks (PINNs) to approximate the solutions to the PEs and study the error estimates. We first establish the higher-order regularity for the global solutions to the PEs with either full viscosity and diffusivity, or with only the horizontal ones. Such a result for the case with only the horizontal ones is new and required in the analysis under the PINNs framework. Then we prove the existence of two-layer tanh PINNs of which the corresponding training error can be arbitrarily small by taking the width of PINNs to be sufficiently wide, and the error between the true solution and its approximation can be arbitrarily small provided that the training error is small enough and the sample set is large enough. In particular, all the estimates are a priori, and our analysis includes higher-order (in spatial Sobolev norm) error estimates. Numerical results on prototype systems are presented to further illustrate the advantage of using the H^s norm during the training.

A posteriori error estimates for wave maps into spheres

We provide a posteriori error estimates in the energy norm for temporal semi-discretisations of wave maps into spheres that are based on the angular momentum formulation. Our analysis is based on novel weak–strong stability estimates which we combine with suitable reconstructions of the numerical solution. We present time-adaptive numerical simulations based on the a posteriori error estimators for solutions involving blow-up.

L 2 norm error estimates of BDF methods up to fifth-order for the phase field crystal model

Abstract The well-known backward difference formula (BDF) of the third, the fourth and the fifth orders are investigated for time integration of the phase field crystal model. By building up novel discrete gradient structures of the BDF-$\textrm{k}$ ($\textrm{k}=3,4,5$) formulas, we establish energy dissipation laws at the discrete levels and then obtain a priori solution estimates for the associated numerical schemes; however, we can not build any discrete energy dissipation law for the corresponding BDF-6 scheme because the BDF-6 formula itself does not have any discrete gradient structures. With the help of the discrete orthogonal convolution kernels and Young-type convolution inequalities, some concise $L^{2}$ norm error estimates (with respect to the starting data in the $L^{2}$ norm) are established via the discrete energy technique. To the best of our knowledge, this is the first time such type $L^{2}$ norm error estimates of non-A-stable BDF schemes are obtained for nonlinear parabolic equations. Numerical examples are presented to verify and support the theoretical analysis.

A modified weak Galerkin method for (curl)-elliptic problem

In this paper, we design and analyze a modified weak Galerkin (MWG) finite element method for H ( curl ) -elliptic problem. We first introduce a new discrete weak curl operator and the MWG finite element space. The modified weak Galerkin method does not require the penalty parameter by comparing with traditional DG methods. We prove an optimal error estimate in energy norm. At last, we provide the numerical results to confirm these theoretical results.

L2 estimates for weak Galerkin finite element methods for second-order wave equations with polygonal meshes

In this paper, we describe weak Galerkin finite element methods for solving hyperbolic problems on polygonal meshes. We propose both semidiscrete and fully discrete schemes to numerically solve the second-order linear wave equation. For the time discretization, we have used implicit second order Newmark scheme. For sufficiently smooth solutions, optimal order error estimate in the L2 norm is shown to hold as O(hk+1+τ2), where h is the mesh size and τ the time step. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the proposed method.

Adaptive Quadratic Finite Element Method for the Unilateral Contact Problem

In this paper, we present and analyze a posteriori error estimates in the energy norm of a quadratic finite element method for the frictionless unilateral contact problem. The reliability and the efficiency of a posteriori error estimator is discussed. The suitable decomposition of the discrete space $${\mathbf{V^h_0}}$$ and a discrete space $${\mathbf{Q^h}}$$ , where the discrete counterpart of the contact force density is defined, play crucial role in deriving a posteriori error estimates. Numerical results are presented exhibiting the reliability and the efficiency of the proposed error estimator.

A mass‐ and energy‐preserving numerical scheme for solving coupled Gross–Pitaevskii equations in high dimensions

AbstractThis article is concerned with numerical study of a coupled system of Gross–Pitaevskii equations which describes the spin‐orbit‐coupled Bose–Einstein condensates. Due to the fact that this system possesses the total mass and energy conservation property and often appears in high dimensions, it brings a significant burden in designing and analyzing a suitable numerical scheme for solving the coupled Gross–Pitaevskii equations (CGPEs). In this article, an implicit finite difference scheme is proposed to solve the CGPEs, which is proved to be uniquely solvable, mass‐ and energy‐conservative in the discrete sense. In particular, it is proved in a rigorous way that, without any grid‐ratio restriction, the scheme is stable and convergent at the rate of with time step and mesh size in the maximum norm, while previous works often require certain restriction on the grid ratio and only give the error estimates in the discrete norm or norm which could not imply the maximum error estimate. Numerical results are carried out to underline the error estimate and conservation laws, and investigate several dynamics of the CGPEs.

NIPG Finite Element Method for Convection-Dominated Diffusion Problems with Discontinuous Data

This paper presents the nonsymmetric interior penalty Galerkin (NIPG) finite element method for a class of one-dimensional convection dominated diffusion problems with discontinuous coefficients. The solution of the considered class of problem exhibits boundary and interior layers. Piecewise uniform Shishkin-type meshes are used for the spatial discretization. The error estimates in the energy norm have been derived for the proposed schemes. Theoretical results are supported by conducting numerical experiments. It is established that the errors are uniform with respect to the perturbation parameter [Formula: see text]. The uniformness of the error estimates with the perturbation parameter [Formula: see text] has also been established numerically for [Formula: see text]- norm.

A hybrid high-order scheme for the stationary, incompressible magnetohydrodynamics equations

Abstract We propose and analyse a hybrid high-order scheme for the stationary incompressible magnetohydrodynamics equations. The scheme has an arbitrary order of accuracy and is applicable on generic polyhedral meshes. For sources that are small enough, we prove error estimates in energy norm for the velocity and magnetic field, and $L^2$-norm for the pressure; these estimates are fully robust with respect to small faces, and of optimal order with respect to the mesh size. Using compactness techniques, we also prove that the scheme converges to a solution of the continuous problem, irrespective of the source being small or large. Finally, we illustrate our theoretical results through 3D numerical tests on tetrahedral and Voronoi mesh families.

Discontinuous Galerkin approximations to elliptic and parabolic problems with a Dirac line source

The analyses of interior penalty discontinuous Galerkin methods of any order k for solving elliptic and parabolic problems with Dirac line sources are presented. For the steady state case, we prove convergence of the method by deriving a priori error estimates in the L2 norm and in weighted energy norms. In addition, we prove almost optimal local error estimates in the energy norm for any approximation order. Further, almost optimal local error estimates in the L2 norm are obtained for the case of piecewise linear approximations whereas suboptimal error bounds in the L2 norm are shown for any polynomial degree. For the time-dependent case, convergence of semi-discrete and of backward Euler fully discrete scheme is established by proving error estimates in L2 in time and in space. Numerical results for the elliptic problem are added to support the theoretical results.