Search

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. [ABSTRACT FROM AUTHOR]

In a conforming discontinuous Galerkin (CDG) finite element method, discontinuous P k polynomials are employed. To connect discontinuous functions, the inter-element traces, { u h } and { ∇ u h } , are usually defined as some averages in discontinuous Galerkin finite element methods. But in this CDG finite element method, they are defined as projections of a lifted P k + 4 polynomial from four P k polynomials on neighboring triangles. With properly chosen weak Hessian spaces, when tested by discontinuous polynomials, the variation form can have no inter-element integral, neither any stabilizer. That is, the bilinear form is the same as that of conforming finite elements for solving the biharmonic equation. Such a conforming discontinuous Galerkin finite element method converges four orders above the optimal order, i.e., the P k solution has an O (h k + 5) convergence in L 2 -norm, and an O (h k + 3) convergence in H 2 -norm. A local post-process is defined, which lifts the P k solution to a P k + 4 quasi-optimal solution. Numerical tests are provided, confirming the theory. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We employ a radial basis function (RBF) - differential quadrature (DQ) method for the numerical solution of elliptic boundary value problems in annular domains. With an appropriate selection of collocation points, for any choice of RBF, both the coefficient and right hand side matrices in the systems appearing in this discretization possess block circulant structures. These linear systems can thus be solved efficiently using matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs). In particular, we consider problems governed by the Poisson equation, the inhomogeneous biharmonic equation and the inhomogeneous Cauchy–Navier equations of elasticity. In addition to its simplicity, the proposed method can both achieve high accuracy and solve large-scale problems. The feasibility of the proposed techniques is illustrated by several numerical examples. [ABSTRACT FROM AUTHOR]

We prove the convergence of an innovative Chebyshev–Gauss–Lobatto (CGL) pseudospectral method applied to fourth order boundary value problems. The proposed method enjoyed all the advantages of pseudospectral methods. Moreover, we can select (N − 3) interior CGL collocation points to enforce the equation, meanwhile we use the remaining four collocation points to assure the boundary conditions. [ABSTRACT FROM AUTHOR]

Abstract: The classical a posteriori error estimates are mostly oriented to the use in the finite element -methods while the contemporary higher-order -methods usually require new approaches in a posteriori error estimation. These methods hold a very important position among adaptive numerical procedures for solving ordinary as well as partial differential equations arising from various technical applications. In the paper, we are concerned with a review and comparison of error estimation procedures for the biharmonic and some more general fourth order partial differential problems with special regards to the needs of the -method. We point out some advantages and drawbacks of analytical and computational a posteriori error estimates. [Copyright &y& Elsevier]

Abstract: We propose a stabilized finite element method for the approximation of the biharmonic equation with a clamped boundary condition. The mixed formulation of the biharmonic equation is obtained by introducing the gradient of the solution and a Lagrange multiplier as new unknowns. Working with a pair of bases forming a biorthogonal system, we can easily eliminate the gradient of the solution and the Lagrange multiplier from the saddle point system leading to a positive definite formulation. Using a superconvergence property of a gradient recovery operator, we prove an optimal a priori estimate for the finite element discretization for a class of meshes. [Copyright &y& Elsevier]

Abstract: In this paper, the Trefftz method of fundamental solution (FS), called the method of fundamental solution (MFS), is used for biharmonic equations. The bounds of errors are derived for the MFS with Almansi’s fundamental solutions (denoted as the MAFS) in bounded simply connected domains. The exponential and polynomial convergence rates are obtained from highly and finitely smooth solutions, respectively. The stability analysis of the MAFS is also made for circular domains. Numerical experiments are carried out for both smooth and singularity problems. The numerical results coincide with the theoretical analysis made. When the particular solutions satisfying the biharmonic equation can be found, the method of particular solutions (MPS) is always superior to the MFS and the MAFS, based on numerical examples. However, if such singular particular solutions near the singular points do not exist, the local refinement of collocation nodes and the greedy adaptive techniques can be used for seeking better source points. Based on the computed results, the MFS using the greedy adaptive techniques may provide more accurate solutions for singularity problems. Moreover, the numerical solutions by the MAFS with Almansi’s FS are slightly better in accuracy and stability than those by the traditional MFS. Hence, the MAFS with the AFS is recommended for biharmonic equations due to its simplicity. [Copyright &y& Elsevier]

Abstract: The solution of boundary value problems (BVP) for fourth order differential equations by their reduction to BVP for second order equations, with the aim to use the achievements for the latter ones attracts attention from many researchers. In this paper, using the technique developed by ourselves in recent works, we construct iterative method for the Neumann BVP for biharmonic type equation. The convergence rate of the method is proved and some numerical experiments are performed for testing it in dependence on the choice of an iterative parameter. [Copyright &y& Elsevier]

This paper presents a mixed finite volume element scheme based on rectangular partition for solving biharmonic equations. It also gives a kind of adaptive Uzawa iteration method for the scheme. It is rigorously proved that the scheme has first-order accuracy in H1 semi-norm and L2 norm according to the characteristics of the scheme. Finally, two numerical examples illustrate the effectiveness of the method. [Copyright &y& Elsevier]

For the global superconvergence over the entire solution domain originated by Lin and his colleagues [7,8], this paper gives a framework of new estimates and new proofs for the basic estimates for bounds of ∫∫lower limit Ω (u−uI)xxvxx ds, ∫∫lower limit Ω (u−uI)xyvxy ds and ∫∫lower limit Ω (u−uI)xxvyy ds, which reveal more intrinsic characteristics and easier understanding and better readable. Suppose that the solution is smooth enough and the solution domain can be split into quasiuniform rectangular elements □ij with the maximal boundary length h. The study of [7,8] dealt with only the clamped boundary condition for biharmonic equations, to obtain the global superconvergence O(h4) in H2 norms under uniform rectangles □ij for the solution u∈H6(Ω). This paper is devoted to other kinds of important boundary conditions, such as the simple support condition, the natural boundary condition and their mixed types where the different boundary conditions are subject to different edges of ∂Ω. New error estimates are derived theoretically, and verified numerically to reach the global superconvergence O(h3.5) and O(h4) for different boundary conditions on different edges of Ω under uniform □ij. Note that the new superconvergence estimates in this paper are essential in practical applications, because different boundary conditions are needed in 3D blending surfaces [5] and in the combined methods for singularity problems [6]. [Copyright &y& Elsevier]