Error Analysis of Compact Finite Difference Schemes for the Biharmonic Equation

Lian, Xiaopeng,Cheng, Xiaoliang

DOI: https://doi.org/10.1109/icise.2010.5688642

2010-01-01

Abstract:In this paper, we borrow some techniques in finite element method to prove some error estimates of some compact finite difference schemes for biharmonic equation with Dirichlet boundary condition.

What problem does this paper attempt to address?

Engineering Analysis Error Estimation when Removing Finite-Sized Features in Nonlinear Elliptic Problems.

Ming Li,Shuming Gao,Ralph R. Martin

DOI: https://doi.org/10.1016/j.cad.2012.10.019

IF: 3.652

2013-01-01

Computer-Aided Design

Abstract:The paper provides novel approaches for a posteriori estimation of goal-oriented engineering analysis error caused by removing finite-sized negative features from a complex model, in the case of analysis of nonlinear elliptic physical phenomena. The features may lie within the model's interior or along its boundary, and may be constrained with either Neumann or Dirichlet boundary conditions. The main use is for deciding whether detail design features can be removed from a model, to simplify meshing and engineering analysis, without unduly affecting analysis results. Error estimates are found using adjoint theory. Using a rigorous mathematical derivation, the error is first reformulated as a local quantity defined over the boundary of the feature to be suppressed, via linearization and Green's theorem. This intermediate result still involves unknown terms, which we overcome in three ways. In one, an approximate upper bound of the error is obtained rigorously utilizing classical theories of differential operators; the others are heuristic practical approaches. The performance and the effectivity of these three different approaches are examined on 2D and 3D internal and boundary features, with Neumann and Dirichlet boundary conditions.
Some Mixed Finite Element Methods For Biharmonic Equation

Xl Cheng,Wm Han,Hc Huang

DOI: https://doi.org/10.1016/S0377-0427(99)00342-8

IF: 2.872

2000-01-01

Journal of Computational and Applied Mathematics

Abstract:Some perturbed mixed finite element methods related to the reduced integration technique are considered for solving the biharmonic equation problem. On a rectangular mesh, a similar scheme was proposed in Malkus and Hughes (Comput. Methods Appl. Mech. Eng. 15 (1978) 63-81) and its convergence was analyzed in Johnson and Pitkaranta (Math. Comp. 38 (1982) 375-400). Here we modify the scheme proposed in Malkus and Hughes (1978) and prove the optimal order error estimate without the extra smoothness assumption on the solution made in Johnson and Pitkaranta (1982). On a triangular mesh, an analogous scheme is studied, and an order error estimate is proved. Some numerical results are given to show the convergence behavior of the numerical solutions. (C) 2000 Elsevier Science B.V. All rights reserved. MSG: 65N30.
Fourth-order compact finite difference schemes for solving biharmonic equations with Dirichlet boundary conditions

Kejia Pan,Jin Li,Zhilin Li,Kang Fu

2024-09-02

Abstract:In this study, we propose a genuine fourth-order compact finite difference scheme for solving biharmonic equations with Dirichlet boundary conditions in both two and three dimensions. In the 2D case, we build upon the high-order compact (HOC) schemes for flux-type boundary conditions originally developed by Zhilin Li and Kejia Pan [SIAM J. Sci. Comput., 45 (2023), pp. A646-A674] to construct a high order compact discretization for coupled boundary conditions. When considering the 3D case, we modify carefully designed undetermined coefficient methods of Li and Pan to derive the finite difference approximations of coupled boundary conditions. The resultant FD discretization maintains the global fourth order convergence and compactness. Unlike the very popular Stephenson method, the number of unknows do not increase with dimensions. Besides, it is noteworthy that the condition number of the coefficient matrix increases at a rate of $O(h^{-2})$ in both 2D and 3D. We also validate the performance of the proposed genuine HOC methods through nontrivial examples.

Numerical Analysis
Accuracy Analysis of the Adini Element for Biharmonic Equation

Ping Luo,Qun Lin

DOI: https://doi.org/10.1007/s10114-004-0317-3

2004-01-01

Abstract:In this paper, we consider the solution of the biharmonic equation using Adini nonconforming finite element, and report new results of the asymptotic error expansions of the interpolation error functionals and nonconforming remainder. These expansions are used to develop two extrapolation formulas and a series of sharp error estimates. Finally, the numerical results have verified the extrapolation theory.
Multiscale Finite Element Discretizations Based on Local Defect Correction for the Biharmonic Eigenvalue Problem of Plate Buckling

Shijie Wang,Yidu Yang,Hai Bi

DOI: https://doi.org/10.1002/mma.5409

2018-01-01

Mathematical Methods in the Applied Sciences

Abstract:In this paper, we study the multiscale finite element discretizations about the biharmonic eigenvalue problem of plate buckling. On the basis of the work of Dai and Zhou (SIAM J. Numer. Anal. 46[1] [2008] 295‐324), we establish a three‐scale scheme, a multiscale discretization scheme, and the associated parallel version based on local defect correction. We first prove a local priori error estimate of finite element approximations, then give the error estimates of multiscale discretization schemes. Theoretical analysis and numerical experiments indicate that our schemes are suitable and efficient for eigenfunctions with local low smoothness.
Finite Difference Method for a Semilinear Parabolic Equations with Variable Coefficients

WANG Hai-ming

DOI: https://doi.org/10.3969/j.issn.1009-4822.2007.03.001

2007-01-01

Abstract:A finite difference method is proposed for a semilinear parabolic equation with variable coefficients.The scheme resulting form this method is linear in the sense that it requires only to solve a single linear algebraic system at each time step.The existence-uniqueness and convergence of difference solution are proved.The stability of the difference scheme is also proved.It is shown that the order of accuracy of the difference solution in L∞ and L2 norms is O(h2+τ2).
Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

Lavanya V Salian,Samala Rathan,Rakesh Kumar

2024-11-27

Abstract:This study aims to construct a stable, high-order compact finite difference method for solving Sobolev-type equations with Dirichlet boundary conditions in one-space dimension. Approximation of higher-order mixed derivatives in some specific Sobolev-type equations requires a bigger stencil information. One can approximate such derivatives on compact stencils, which are higher-order accurate and take less stencil information but are implicit and sparse. Spatial derivatives in this work are approximated using the sixth-order compact finite difference method (Compact6), while temporal derivatives are handled with the explicit forward Euler difference scheme. We examine the accuracy and convergence behavior of the proposed scheme. Using the von Neumann stability analysis, we establish $L_2-$stability theory for the linear case. We derive conditions under which fully discrete schemes are stable. Also, the amplification factor $\mathcal{C}(\theta)$ is analyzed to ensure the decay property over time. Real parts of $\mathcal{C}(\theta)$ lying on the negative real axis confirm the exponential decay of the solution. A series of numerical experiments were performed to verify the effectiveness of the proposed scheme. These tests include advection-free flow, and applications to the equal width equation, such as single solitary wave propagation, interactions of two and three solitary waves, undular bore formation, and the Benjamin-Bona-Mahony-Burgers equation.

Numerical Analysis
The Adaptive Ciarlet–Raviart Mixed Method for Biharmonic Problems with Simply Supported Boundary Condition

Yidu Yang,Hai Bi,Yu Zhang

DOI: https://doi.org/10.1016/j.amc.2018.07.014

IF: 4.397

2018-01-01

Applied Mathematics and Computation

Abstract:In this paper, we study the adaptive fashion of the Ciarlet–Raviart mixed method for biharmonic equation/eigenvalue problem with simply supported boundary condition in Rd. We propose an a posteriori error indicator of the Ciarlet–Raviart approximate solution for the biharmonic equation and an a posteriori error indicator of the Ciarlet–Raviart approximate eigenfuction, and prove the reliability and efficiency of the indicators. We also give an a posteriori error indicator for the approximate eigenvalue and prove its reliability. We design an adaptive Ciarlet–Raviart mixed method with piecewise polynomials of degree less than or equal to m, and numerical experiments show that numerical eigenvalues obtained by the method can achieve the optimal convergence order O(dof−2md)(d=2,m=2,3;d=3,m=3).
Continuous Linear Finite Element Method for Biharmonic Problems on Surfaces

Ying Cai,Hailong Guo,Zhimin Zhang

2024-04-28

Abstract:This paper presents an innovative continuous linear finite element approach to effectively solve biharmonic problems on surfaces. The key idea behind this method lies in the strategic utilization of a surface gradient recovery operator to compute the second-order surface derivative of a piecewise continuous linear function defined on the approximate surface, as conventional notions of second-order derivatives are not directly applicable in this context. By incorporating appropriate stabilizations, we rigorously establish the stability of the proposed formulation. Despite the presence of geometric error, we provide optimal error estimates in both the energy norm and $L^2$ norm. Theoretical results are supported by numerical experiments.

Numerical Analysis
A Hybridized Weak Galerkin Finite Element Method for the Biharmonic Equation

Chunmei Wang,Junping Wang

DOI: https://doi.org/10.48550/arXiv.1402.1157

2014-02-06

Abstract: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.

Numerical Analysis
A novel difference equation approach for the stability and robustness of compact schemes for variable coefficient PDEs

Anindya Goswami,Kuldip Singh Patel,Pradeep Kumar Sahu

2024-01-28

Abstract:Fourth-order accurate compact schemes for variable coefficient convection diffusion equations are considered. A sufficient condition for the stability of the fully discrete problem is derived using a difference equation based approach. The constant coefficient problems are considered as a special case, and the unconditional stability of compact schemes for such case is proved theoretically. The condition number of the amplification matrix is also analyzed, and an estimate for the same is derived. The examples are provided to support the assumption taken to assure stability.

Numerical Analysis,Computational Finance
The a posteriori error estimates of the Ciarlet-Raviart mixed finite element method for the biharmonic eigenvalue problem

Jinhua Feng,Shixi Wang,Hai Bi,Yidu Yang

DOI: https://doi.org/10.3934/math.2024163

2024-01-01

AIMS Mathematics

Abstract:The biharmonic equation/eigenvalue problem is one of the fundamental model problems in mathematics and physics and has wide applications. In this paper, for the biharmonic eigenvalue problem, based on the work of Gudi [ Numer. Methods Partial Differ. Equ. , 27 (2011), 315-328], we study the a posteriori error estimates of the approximate eigenpairs obtained by the Ciarlet-Raviart mixed finite element method. We prove the reliability and efficiency of the error estimator of the approximate eigenfunction and analyze the reliability of the error estimator of the approximate eigenvalues. We also implement the adaptive calculation and exhibit the numerical experiments which show that our method is efficient and can get an approximate solution with high accuracy.

mathematics, applied
Error analysis of mixed finite element methods for nonlinear parabolic equations

Huadong Gao,Weifeng Qiu

DOI: https://doi.org/10.48550/arXiv.1707.02677

2017-07-10

Abstract:In this paper, we prove a discrete embedding inequality for the Raviart--Thomas mixed finite element methods for second order elliptic equations, which is analogous to the Sobolev embedding inequality in the continuous setting. Then, by using the proved discrete embedding inequality, we provide an optimal error estimate for linearized mixed finite element methods for nonlinear parabolic equations. Several numerical examples are provided to confirm the theoretical analysis.

Numerical Analysis
Error Inhibiting Schemes for Initial Boundary Value Heat Equation

Adi Ditkowski,Paz Fink Shustin

DOI: https://doi.org/10.48550/arXiv.2010.00476

2021-02-16

Abstract:Finite Difference (FD) schemes are widely used in science and engineering for approximating solutions of partial differential equations (PDEs). Error analysis of FD schemes relies on estimating the truncation error at each time step. This approach usually leads to a global error whose order is of the same order of the truncation error. For classical FD schemes the global error is indeed of the same order as the truncation error. A particular class of FD schemes is the Block Finite Difference (BFD) schemes, in which the grid is divided into blocks. The structure of such schemes is similar to the structure of the Discontinuous Galerkin (DG) method, and allows inhabitation of the truncation errors. Recently, much effort was devoted to design BFD schemes whose global error converges faster than the truncation error. In this paper, we elaborate the approach presented in <a class="link-https" data-arxiv-id="1711.07926" href="https://arxiv.org/abs/1711.07926">arXiv:1711.07926</a> for the heat equation with periodic boundary conditions. We generalize this methodology to design BFD schemes for the heat equation with Dirichlet or Neumann boundary conditions, whose global error converges faster than the truncation error. Such schemes are henceforth called Error Inhibiting Schemes. We provide an explicit error analysis, including proofs of stability and convergence of the proposed schemes. We illustrate our approach using several numerical examples, which demonstrate the efficiency of our method in comparison to standard FD schemes.

Numerical Analysis
Lowest-degree robust finite element schemes for inhomogeneous bi-Laplace problems

Bin Dai,Huilan Zeng,Chensong Zhang,Shuo Zhang

2024-04-21

Abstract:In this paper, we study the numerical method for the bi-Laplace problems with inhomogeneous coefficients; particularly, we propose finite element schemes on rectangular grids respectively for an inhomogeneous fourth-order elliptic singular perturbation problem and for the Helmholtz transmission eigenvalue problem. The new methods use the reduced rectangle Morley (RRM for short) element space with piecewise quadratic polynomials, which are of the lowest degree possible. For the finite element space, a discrete analogue of an equality by Grisvard is proved for the stability issue and a locally-averaged interpolation operator is constructed for the approximation issue. Optimal convergence rates of the schemes are proved, and numerical experiments are given to verify the theoretical analysis.

Numerical Analysis
Error analysis of a fully discrete finite element method for variable density incompressible flows in two dimensions

Wentao Cai,Buyang Li,Ying Li

DOI: https://doi.org/10.1051/m2an/2020029

2021-01-01

ESAIM: Mathematical Modelling and Numerical Analysis

Abstract:An error estimate is presented for a fully discrete, linearized and stabilized finite element method for solving the coupled system of nonlinear hyperbolic and parabolic equations describing incompressible flow with variable density in a two-dimensional convex polygon. In particular, the error of the numerical solution is split into the temporal and spatial components, separately. The temporal error is estimated by applying discrete maximal L p -regularity of time-dependent Stokes equations, and the spatial error is estimated by using energy techniques based on the uniform regularity of the solutions given by semi-discretization in time.
The stabilizer free weak Galerkin mixed finite elements method for the biharmonic equation

Shanshan Gu,Fuchang Huo,Shicheng Liu

2024-03-19

Abstract:In this article, the stabilizer free weak Galerkin (SFWG) finite element method is applied to the Ciarlet-Raviart mixed form of the Biharmonic equation. We utilize the SFWG solutions of the second elliptic problems to define projection operators, build error equations, and further derive the error estimates. Finally, numerical examples support the results reached by the theory.

Numerical Analysis
Anisotropic Error Analysis of Weak Galerkin finite element method for Singularly Perturbed Biharmonic Problems

Aayushman Raina,Srinivasan Natesan,Şuayip Toprakseven

2024-09-11

Abstract:We consider the Weak Galerkin finite element approximation of the Singularly Perturbed Biharmonic elliptic problem on a unit square domain with clamped boundary conditions. Shishkin mesh is used for domain discretization as the solution exhibits boundary layers near the domain boundary. Error estimates in the equivalent $H^{2}-$ norm have been established and the uniform convergence of the proposed method has been proved. Numerical examples are presented corroborating our theoretical findings.

Numerical Analysis
An optimal error estimate for a mixed finite element method for a biharmonic problem with clamped boundary conditions

Bishnu P Lamichhane

2023-05-17

Abstract:We consider a mixed finite element method for a biharmonic equation with clamped boundary conditions based on biorthogonal systems with weakly imposed Dirichlet boundary condition. We show that the weak imposition of the boundary condition arising from a natural minimisation formulation allows to get an optimal a priori error estimate for the finite element scheme improving the existing error estimate for such a formulation without weakly imposed Dirichlet boundary condition. We also briefly outline the algebraic formulation arising from the finite element method.

Numerical Analysis
Fast compact finite difference schemes on graded meshes for fourth-order multi-term fractional sub-diffusion equations with the first Dirichlet boundary conditions

Zhibo Wang,Caixia Ou,Dakang Cen

DOI: https://doi.org/10.1080/00207160.2022.2119080

IF: 1.75

2022-09-20

International Journal of Computer Mathematics

Abstract:In this paper, fast compact finite difference schemes are derived for fourth-order multi-term fractional sub-diffusion equations with initial singularity under the first Dirichlet boundary conditions. In contrast to the direct scheme, the fast algorithm adopted to approximate the Caputo derivative reduces the computation costs effectively. Sharp error estimate of the proposed scheme for linear models is rigorously presented by the energy method. Furthermore, to handle more intricate nonlinear models, a crucial Grönwall inequality is then deduced to analyse the stability and convergence of the linearized scheme. It is worth pointing out that the Grönwall inequality is also helpful in numerical analysis of multi-step schemes for other problems, such as, integro-differential equations with multiple fractional derivatives. Ultimately, numerical examples are provided to verify the efficiency of the established difference schemes.

Error Analysis of Compact Finite Difference Schemes for the Biharmonic Equation

Engineering Analysis Error Estimation when Removing Finite-Sized Features in Nonlinear Elliptic Problems.

Some Mixed Finite Element Methods For Biharmonic Equation

Fourth-order compact finite difference schemes for solving biharmonic equations with Dirichlet boundary conditions

Accuracy Analysis of the Adini Element for Biharmonic Equation

Multiscale Finite Element Discretizations Based on Local Defect Correction for the Biharmonic Eigenvalue Problem of Plate Buckling

Finite Difference Method for a Semilinear Parabolic Equations with Variable Coefficients

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

The Adaptive Ciarlet–Raviart Mixed Method for Biharmonic Problems with Simply Supported Boundary Condition

Continuous Linear Finite Element Method for Biharmonic Problems on Surfaces

A Hybridized Weak Galerkin Finite Element Method for the Biharmonic Equation

A novel difference equation approach for the stability and robustness of compact schemes for variable coefficient PDEs

The a posteriori error estimates of the Ciarlet-Raviart mixed finite element method for the biharmonic eigenvalue problem

Error analysis of mixed finite element methods for nonlinear parabolic equations

Error Inhibiting Schemes for Initial Boundary Value Heat Equation

Lowest-degree robust finite element schemes for inhomogeneous bi-Laplace problems

Error analysis of a fully discrete finite element method for variable density incompressible flows in two dimensions

The stabilizer free weak Galerkin mixed finite elements method for the biharmonic equation

Anisotropic Error Analysis of Weak Galerkin finite element method for Singularly Perturbed Biharmonic Problems

An optimal error estimate for a mixed finite element method for a biharmonic problem with clamped boundary conditions

Fast compact finite difference schemes on graded meshes for fourth-order multi-term fractional sub-diffusion equations with the first Dirichlet boundary conditions