Graded mesh refinement and error estimates of higher order for DGFE solutions of elliptic boundary value problems in polygons

Error estimates for DGFE solutions are well investigated if one assumes that the exact solution is sufficiently regular. In this article, we consider a Dirichlet and a mixed boundary value problem for a linear elliptic equation in a polygon. It is well known that the first derivatives of the solutions develop singularities near reentrant corner points or points where the boundary conditions change. On the basis of the regularity results formulated in Sobolev–Slobodetskii spaces and weighted spaces of Kondratiev type, we prove error estimates of higher order for DGFE solutions using a suitable graded mesh refinement near boundary singular points. The main tools are as follows: regularity investigation for the exact solution relying on general results for elliptic boundary value problems, error analysis for the interpolation in Sobolev–Slobodetskii spaces, and error estimates for DGFE solutions on special graded refined meshes combined with estimates in weighted Sobolev spaces. Our main result is that there exist a local grading of the mesh and a piecewise interpolation by polynoms of higher degree such that we will get the same order O (h^α) of approximation as in the smooth case. © 2011 Wiley Periodicals, Inc. Numer Mehods Partial Differential Eq, 2012     

New fixed point theorems for mixed monotone operators and local existence-uniqueness of positive solutions for nonlinear boundary value problems

In this article we present a new fixed point theorem for a class of general mixed monotone operators, which extends the existing corresponding results. Moreover, we establish some pleasant properties of nonlinear eigenvalue problems for mixed monotone operators. Based on them the local existence-uniqueness of positive solutions for nonlinear boundary value problems which include Neumann boundary value problems, three-point boundary value problems and elliptic boundary value problems for Lane-Emden-Fowler equations is proved. The theorems for nonlinear boundary value problems obtained here are very general.     

Accelerating Galerkin BEM for linear elasticity using adaptive cross approximation

The adaptive cross approximation (ACA) algorithm (Numer. Math. 2000; 86 :565–589; Computing 2003; 70 (1):1–24) provides a means to compute data‐sparse approximants of discrete integral formulations of elliptic boundary value problems with almost linear complexity. ACA uses only few of the original entries for the approximation of the whole matrix and is therefore well‐suited to speed up existing computer codes. In this article we extend the convergence proof of ACA to Galerkin discretizations. Additionally, we prove that ACA can be applied to integral formulations of systems of second‐order elliptic operators without adaptation to the respective problem. The results of applying ACA to boundary integral formulations of linear elasticity are reported. Furthermore, we comment on recent implementation issues of ACA for non‐smooth boundaries. Copyright © 2006 John Wiley & Sons, Ltd.     

无界区域抛物方程自然边界元方法

本文应用自然边界元方法求解无界区域抛物型初边值问题。首先将控制方程对时间进行离散化,得到关于时间步长离散化的椭圆型问题。通过Fourier展开,导出相应问题的自然积分方程和Poisson积分公式。研究了自然积分算子的性质,并讨论了自然积分方程的数值解法,最后给出数值例子。从而解决了抛物型问题的自然边界归化和自然边界元方法。     

Error analysis of staggered finite difference finite volume schemes on unstructured meshes

This work combines the consistency in lower‐order differential operators with external approximations of functional spaces to obtain error estimates for finite difference finite volume schemes on unstructured nonuniform meshes. This combined approach is first applied to a one‐dimensional elliptic boundary value problem on nonuniform meshes, and a first‐order convergence rate is obtained, which agrees with the results previously reported. The approach is also applied to the staggered Marker‐and‐Cell scheme for the two‐dimensional incompressible Stokes problem on unstructured meshes. A first‐order convergence rate is obtained, which improves over a previously reported result in that it also holds on unstructured meshes. For both problems considered in this work, the convergence rate is one order higher on meshes satisfying special requirements. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1159–1182, 2017     

The Dual Singular Function Method for 2D Boundary Integral Equations

The accuracy of standard boundary element methods for elliptic boundary value problems deteriorates if the boundary of the domain contains corners or if the boundary conditions change along the boundary. Here we first investigate the convergence behaviour of standard spline Galerkin approximation on quasi-uniform meshes for boundary integral equations on polygonal domains. It turns out, that the order of convergence depends on some constant describing the singular behaviour of solutions near corner points of the boundary. In order to recover the full order of convergence for the Galerkin approximation we propose the dual singular function method which is often used for improving the accuracy of finite element methods. The theoretical convergence results are confirmed and illustrated by a numerical example.     

Non‐linear initial boundary value problemsof hyperbolic–parabolic type. A general investigationof admissible couplings between systems of higher order. Part 1: linear theory

In this article (which is divided in three parts) we investigate the non‐linear initial boundary value problems (1.2) and (1.3). In both cases we consider coupled systems where each system is of higher order and of hyperbolic or parabolic type. Our goal is to characterize systematically all admissible couplings between systems of higher order and different type. By an admissible coupling we mean a condition that guarantees the existence, uniqueness and regularity of solutions to the respective initial boundary value problem. In part 1 at hand, we develop the underlying theory of linear hyperbolic and parabolic initial boundary value problems. Testing the PDEs with suitable functions we obtain a priori estimates for the respective solutions. In particular, we make use of the regularity theory for linear elliptic boundary value problems that was previously developed by the author. In part 2, we prove the local in time existence, uniqueness and regularity of solutions to the quasilinear initial boundary value problem (1.2) using the so‐called energy method. In the above sense, the regularity assumptions about the coefficients and right‐hand sides define the admissible couplings. In part 3, we extend the results of part 2 to the non‐linear initial boundary value problem (1.3). In particular, the assumptions about the respective parameters correspond to the previous regularity assumptions and hence define the admissible couplings now. Moreover, we exploit the assumptions about the respective parameters for the case of two coupled systems. Copyright © 2002 John Wiley & Sons, Ltd.     

A fractional splitting algorithm for nonoverlapping domain decomposition for parabolic problem

In this article we study the convergence of the nonoverlapping domain decomposition for solving large linear system arising from semi‐discretization of two‐dimensional initial value problem with homogeneous boundary conditions and solved by implicit time stepping using first and two alternatives of second‐order FS‐methods. The interface values along the artificial boundary condition line are found using explicit forward Euler's method for the first‐order FS‐method, and for the second‐order FS‐method to use extrapolation procedure for each spatial variable individually. The solution by the nonoverlapping domain decomposition with FS‐method is applicable to problems that requires the solution on nonuniform meshes for each spatial variable, which will enable us to use different time‐stepping over different subdomains and with the possibility of extension to three‐dimensional problem. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 609–624, 2002     

A hybrid high-order formulation for a Neumann problem on polytopal meshes

In this work, we study a hybrid high-order (HHO) method for an elliptic diffusion problem with Neumann boundary condition. The proposed method has several features, such as: (a) the support of arbitrary approximation order polynomial at mesh elements and faces on polytopal meshes, (b) the design of a local (element-wise) potential reconstruction operator and a local stabilization term, that weakly enforces the matching between local element- and face-based on degrees of freedom, and (c) cheap computational cost, thanks to static condensation and compact stencil. We prove the well-posedness of our HHO formulation, and obtain the optimal error estimates, according to previous study. Implementation aspects are thoroughly discussed. Finally, some numerical examples are provided, which are in agreement with our theoretical results.     

Convergence and superconvergence analysis of a new quadratic Hermite-type triangular element on anisotropic meshes

A new quadratic Hermite-type triangular finite element is conceived to solve a class of two-dimensional second-order elliptic boundary value problems. Its error estimates on anisotropic meshes are developed. Furthermore, we verify that some conditions set to the meshes contribute to the proof of its superconvergence properties, which can improve the approximation results. Numerical examples are given to confirm our theoretical analysis.     

Error estimation of a quadratic finite volume method on right quadrangular prism grids

In this paper, we develop a finite volume element method with affine quadratic bases on right quadrangular prism meshes for three-dimensional elliptic boundary value problems. The optimal H¹$H^1$-norm error estimate of second order accuracy is proved under certain assumptions about the meshes. Numerical results are presented to illustrate the theoretical analysis.     

An iterative method based on equation decomposition for the fourth‐order singular perturbation problem

In this article, we propose an iterative method based on the equation decomposition technique ( 1 ) for the numerical solution of a singular perturbation problem of fourth‐order elliptic equation. At each step of the given method, we only need to solve a boundary value problem of second‐order elliptic equation and a second‐order singular perturbation problem. We prove that our approximate solution converges to the exact solution when the domain is a disc. Our numerical examples show the efficiency and accuracy of our method. Our iterative method works very well for singular perturbation problems, that is, the case of 0 < ε ? 1, and the convergence rate is very fast. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

一个Crouzeix-Raviart型有限元方法关于二阶椭圆问题超逼近的一个新证明

In this paper,a new proof of superclose of a Crouzeix-Raviart type finite element is given for second order elliptic boundary value problem by Bramble-Hilbert lemma on anisotropic meshes.     

A domain decomposition preconditioner for hermite collocation problems

We propose a preconditioning method for linear systems of equations arising from piecewise Hermite bicubic collocation applied to two‐dimensional elliptic PDEs with mixed boundary conditions. We construct an efficient, parallel preconditioner for the GMRES method. The main contribution of the article is a novel interface preconditioner derived in the framework of substructuring and employing a local Hermite collocation discretization for the interface subproblems based on a hybrid fine‐coarse mesh. Interface equations based on this mesh depend only weakly on unknowns associated with subdomains. The effectiveness of the proposed method is highlighted by numerical experiments that cover a variety of problems. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 135–151, 2003     

A symmetric integrated radial basis function method for solving differential equations

In this article, integrated radial basis functions (IRBFs) are used for Hermite interpolation in the solution of differential equations, resulting in a new meshless symmetric RBF method. Both global and local approximation‐based schemes are derived. For the latter, the focus is on the construction of compact approximation stencils, where a sparse system matrix and a high‐order accuracy can be achieved together. Cartesian‐grid‐based stencils are possible for problems defined on nonrectangular domains. Furthermore, the effects of the RBF width on the solution accuracy for a given grid size are fully explored with a reasonable computational cost. The proposed schemes are numerically verified in some elliptic boundary‐value problems governed by the Poisson and convection‐diffusion equations. High levels of the solution accuracy are obtained using relatively coarse discretisations.     

Numerical methods and error analysis for one dimensional elliptic problems containing singularities

The method of auxiliary mapping (MAM), introduced by Babu?ka and Oh, was proven to be very successful in dealing with monotone singularities arising in two‐dimensional problems. In this article, in the framework of the p‐version of FEM, MAM is presented for one‐dimensional elliptic boundary value problems containing singularities. Moreover, in order to show the effectiveness of MAM, a detailed proof of an error estimate is also presented, which gives a sharp error bound of MAM. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 399–420, 2003.     

Non‐linear initial boundary value problems of hyperbolic–parabolic type. A general investigation of admissible couplings between systems of higher order. Part 2: abstract quasilinear theory

This is the second part of an article that is devoted to the theory of non‐linear initial boundary value problems. We consider coupled systems where each system is of higher order and of hyperbolic or parabolic type. Our goal is to characterize systematically all admissible couplings between systems of higher order and different type. By an admissible coupling we mean a condition that guarantees the existence, uniqueness and regularity of solutions to the respective initial boundary value problem. In part 1, we develop the underlying theory of linear hyperbolic and parabolic initial boundary value problems. Testing the PDEs with suitable functions we obtain a priori estimates for the respective solutions. In particular, we make use of the regularity theory for linear elliptic boundary value problems that was previously developed by the author. In part 2 at hand, we prove the local in time existence, uniqueness and regularity of solutions to the quasilinear initial boundary value problem (3.4) using the so‐called energy method. In the above sense the regularity assumptions (A6) and (A7) about the coefficients and right‐hand sides define the admissible couplings. In part 3, we extend the results of part 2 to non‐linear initial boundary value problems. In particular, the assumptions about the respective parameters correspond to the previous regularity assumptions and hence define the admissible couplings now. Moreover, we exploit the assumptions about the respective parameters for the case of two coupled systems. Copyright © 2002 John Wiley & Sons, Ltd.     

A new superconvergence property of Wilson nonconforming finite element

Summary. In this paper the Wilson nonconforming finite element method is considered to solve a class of two-dimensional second-order elliptic boundary value problems. A new superconvergence property at the vertices and the midpoints of four edges of rectangular meshes is obtained. Received May 5, 1995 / Revised version received November 11, 1996     

各向异性网格下Wilson元的超收敛性分析

在各向异性网格下研究了二阶椭圆边值问题的Wilson有限元方法,利用单元构造的特殊性和一些新的技巧得到相应的超逼近和超收敛结果.数值算例的结果与理论分析是相吻合的.     

Non‐linear initial boundary value problems of hyperbolic–parabolic type. A general investigation of admissible couplings between systems of higher order. Part 3: applications

This is the third part of an article that is devoted to the theory of non‐linear initial boundary value problems. We consider coupled systems where each system is of higher order and of hyperbolic or parabolic type. Our goal is to characterize systematically all admissible couplings between systems of higher order and different type. By an admissible coupling we mean a condition that guarantees the existence, uniqueness and regularity of solutions to the respective initial boundary value problem. In part 1, we develop the underlying theory of linear hyperbolic and parabolic initial boundary value problems. Testing the PDEs with suitable functions we obtain a priori estimates for the respective solutions. In particular, we make use of the regularity theory for linear elliptic boundary value problems that was previously developed by the author. In part 2, we prove the local in time existence, uniqueness and regularity of solutions to quasilinear initial boundary value problems using the so‐called energy method. In the above sense the regularity assumptions about the coefficients and right‐hand sides define the admissible couplings. In part 3 at hand, we extend the results of part 2 to the nonlinear initial boundary value problem (4.2). In particular, assumptions (B8) and (B9) about the respective parameters correspond to the previous regularity assumptions and hence define the admissible couplings now. Moreover, we exploit assumptions (B8) and (B9) for the case of two coupled systems. Copyright © 2002 John Wiley & Sons, Ltd.