The Morley element for fourth order elliptic equations in any dimensions

In this paper, the well-known nonconforming Morley element for biharmonic equations in two spatial dimensions is extended to any higher dimensions in a canonical fashion. The general n-dimensional Morley element consists of all quadratic polynomials defined on each n-simplex with degrees of freedom given by the integral average of the normal derivative on each (n-1)-subsimplex and the integral average of the function value on each (n-2)-subsimplex. Explicit expressions of nodal basis functions are also obtained for this element on general n-simplicial grids. Convergence analysis is given for this element when it is applied as a nonconforming finite element discretization for the biharmonic equation. The work was supported by the National Natural Science Foundation of China (10571006). This work was supported in part by NSF DMS-0209497 and NSF DMS-0215392 and the Changjiang Professorship through Peking University.     

A second order characteristic finite element scheme for convection-diffusion problems

Summary. A new characteristic finite element scheme is presented for It is of second order accuracy in time increment, symmetric, and unconditionally stable. Optimal error estimates are proved in the framework of -theory. Numerical results are presented for two examples, which show the advantage of the scheme. Received November 22, 2000 / Revised version received July 11, 2001 / Published online October 17, 2001     

Convergence of solutions of a second order Cauchy problem

Let j_νk, y_νk and c_νk denote the kth positive zeros of the Bessel functions J_ν(x), Y_ν(x) and of the general cylinder function C_ν(x), respectively. We show, among other things, that, for k = 2, 3,… and 0 < ν < ∞, c_νk is a concave function of ν, c_νk > ν + c_0k and $\text{c}{}_{\mathrm{\nu k}}\text{[v + (}\text{2}\text{π}\text{)c}{}_{\mathrm{0k}}\text{]}$ decreases as ν increases. In the cases of j_νk and y_νk, these results hold also for k = 1.     

Convergence and complexity of arbitrary order adaptive mixed element methods for the Poisson equation

This paper discusses convergence and complexity of arbitrary,but fixed,order adaptive mixed element methods for the Poisson equation in two and three dimensions.The two main ingredients in the analysis,namely the quasi-orthogonality and the discrete reliability,are achieved by use of a discrete Helmholtz decomposition and a discrete inf-sup condition.The adaptive algorithms are shown to be contractive for the sum of the error of flux in L2-norm and the scaled error estimator after each step of mesh refinement and to be quasi-optimal with respect to the number of elements of underlying partitions.The methods do not require a separate treatment for the data oscillation.     

Convergence and optimality of the adaptive Morley element method

This paper is devoted to the convergence and optimality analysis of the adaptive Morley element method for the fourth order elliptic problem. A new technique is developed to establish a quasi-orthogonality which is crucial for the convergence analysis of the adaptive nonconforming method. By introducing a new parameter-dependent error estimator and further establishing a discrete reliability property, sharp convergence and optimality estimates are then fully proved for the fourth order elliptic problem.     

A second order characteristics finite element scheme for natural convection problems

In this paper a second order characteristics finite element scheme is applied to the numerical solution of natural convection problems. Firstly, after recalling the mathematical model, a second order time discretization of the material time derivative is introduced. Next, fully discretized schemes are proposed by using finite element methods. Numerical results for the two-dimensional problem of buoyancy-driven flow in a square cavity with differentially heated side walls are given and compared with a reference solution.     

A second order difference scheme with nonuniform rectangular meshes for nonlinear parabolic system

In this paper, a difference scheme with nonuniform meshes is proposed for the initial-boundary problem of the nonlinear parabolic system. It is proved that the difference scheme is second order convergent in both space and time.     

Convergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions

In this paper we present a finite element discretization of the Joule-heating problem. We prove existence of solution to the discrete formulation and strong convergence of the finite element solution to the weak solution, up to a sub-sequence. We also present numerical examples in three spatial dimensions. The first example demonstrates the convergence of the method in the second example we consider an engineering application.     

Convergence analysis of a covolume scheme for Maxwell's equations in three dimensions

This paper contains error estimates for covolume discretizations of Maxwell's equations in three space dimensions. Several estimates are proved. First, an estimate for a semi-discrete scheme is given. Second, the estimate is extended to cover the classical interlaced time marching technique. Third, some of our unstructured mesh results are specialized to rectangular meshes, both uniform and nonuniform. By means of some additional analysis it is shown that the spatial convergence rate is one order higher than for the unstructured case.

     

Streamline-diffusion method of a lowest order nonconforming rectangular finite element for convection-diffusion problem

The streamline-diffusion method of the lowest order nonconforming rectangular finite element is proposed for convection-diffusion problem. By making full use of the element’s special property, the same convergence order as the previous literature is obtained. In which, the jump terms on the boundary are added to bilinear form with simple user-chosen parameter δKwhich has nothing to do with perturbation parameter εappeared in the problem under considered, the subdivision mesh size hKand the inverse estimate coefficient μin finite element space.     

Convergence analysis of some first order and second order time accurate gradient schemes for semilinear second order hyperbolic equations

This work is devoted to the convergence analysis of finite volume schemes for a model of semilinear second order hyperbolic equations. The model includes for instance the so‐called Sine‐Gordon equation which appears for instance in Solid Physics (cf. Fang and Li, Adv Math (China) 42 (2013), 441–457; Liu et al., Numer Methods Partial Differ Equ 31 (2015), 670–690). We are motivated by two works. The first one is Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043) where a recent class of nonconforming finite volume meshes is introduced. The second one is Eymard et al. (Numer Math 82 (1999), 91–116) where a convergence of a finite volume scheme for semilinear elliptic equations is provided. The mesh considered in Eymard et al. (Numer Math 82 (1999), 91–116) is admissible in the sense of Eymard et al. (Elsevier, Amsterdam, 2000, 723–1020) and a convergence of a family of approximate solutions toward an exact solution when the mesh size tends to zero is proved. This article is also a continuation of our previous two works (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1278–1321; Bradji, Numer Methods Partial Differ Equ 29 (2013), 1–39) which dealt with the convergence analysis of implicit finite volume schemes for the wave equation. We use as discretization in space the generic spatial mesh introduced in Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043), whereas the discretization in time is performed using a uniform mesh. Two finite volume schemes are derived using the discrete gradient of Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043). The unknowns of these two schemes are the values at the center of the control volumes, at some internal interfaces, and at the mesh points of the time discretization. The first scheme is inspired from the previous work (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1–39), whereas the second one (in which the discretization in time is performed using a Newmark method) is inspired from the work (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1278–1321). Under the assumption that the mesh size of the time discretization is small, we prove the existence and uniqueness of the discrete solutions. If we assume in addition to this that the exact solution is smooth, we derive and prove three error estimates for each scheme. The first error estimate is concerning an estimate for the error between a discrete gradient of the approximate solution and the gradient of the exact solution whereas the second and the third ones are concerning the estimate for the error between the exact solution and the discrete solution in the discrete seminorm of and in the norm of . The convergence rate is proved to be for the first scheme and for the second scheme, where (resp. k) is the mesh size of the spatial (resp. time) discretization. The existence, uniqueness, and convergence results stated above do not require any relation between k and . The analysis presented in this work is also applicable in the gradient schemes framework. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 5–33, 2017     

二阶椭圆问题新混合元模型的超收敛分析及外推

对二阶椭圆问题通过＂增补＂办法导出一个新的混合模型.在各向异性网格下,利用积分恒等式技巧得到了真解与ECHL元近似解的超逼近性质.同时基于插值后处理技术导出了整体超收敛.进一步,通过渐进误差展开和分裂外推,得到了比通常的误差估计更高一阶的收敛速度.     

Convergence of the nonconforming Wilson element for arbitrary quadrilateral meshes

Summary A modified variational formulation, recently introduced by Taylor, Beresford and Wilson for solving second order problems, using the nonconforming Wilson element is here analysed. It is shown that the Patch Test is satisfied and that stresses and displacements are respectively first and second order accurate for arbitrary quadrilateral meshes.     

A quadrature finite element Galerkin scheme for a biharmonic problem on a rectangular polygon

A quadrature Galerkin scheme with the Bogner–Fox–Schmit element for a biharmonic problem on a rectangular polygon is analyzed for existence, uniqueness, and convergence of the discrete solution. It is known that a product Gaussian quadrature with at least three‐points is required to guarantee optimal order convergence in Sobolev norms. In this article, optimal order error estimates are proved for a scheme based on the product two‐point Gaussian quadrature by establishing a relation with an underdetermined orthogonal spline collocation scheme. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Lax-Wendroff-type schemes of arbitrary order in several space dimensions

^{The second-order accurate Lax–Wendroff scheme is basedon the first three terms of a Taylor expansion in time in whichthe time derivatives are replaced by space derivatives usingthe governing evolution equations. The space derivatives arethen approximated by central differencing. In this paper, weextend this idea and truncate the Taylor expansion at an arbitraryorder. One main building block is the so-called Cauchy–Kovalevskayaprocedure to replace all the time derivatives by space derivativeswhich can be formulated for a general system of linear equationswith arbitrary order and in two- or three-space dimensions.The linear case is the main focus of this paper because theproposed high-order schemes are good candidates for the approximationof linear wave motion over long distances and times with importantapplications in aeroacoustics and electromagnetics. The stabilityand the efficiency of Lax–Wendroff-type schemes are examined.The numerical results are compared with a standard scheme foraeroacoustical applications with respect to their quality andthe computational effort. The extensions of the schemes to generalgrids, nonconstant and nonlinear cases are alsoaddressed.}     

Lower bounds of eigenvalues of the biharmonic operators by the rectangular Morley element methods

In this article, we analyze the lower bound property of the discrete eigenvalues by the rectangular Morley elements of the biharmonic operators in both two dimension (2D) and three dimension (3D). The analysis relies on an identity for the errors of eigenvalues. We explore a refined property of the canonical interpolation operators and use it to analyze the key term in this identity. In particular, we show that such a term is of higher order for 2D, and is negative and of second order for 3D, which causes a main difficulty. To overcome it, we propose a novel decomposition of the first term in the aforementioned identity. Finally, we establish a saturation condition to show that the discrete eigenvalues are smaller than the exact ones. We present some numerical results to demonstrate the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1623–1644, 2015     

Interior superconvergence error estimates for mixed finite element methods for second order elliptic problem

1.IntroductionLetfibeaplanedomainwithsmoothboundaryonandWm,p(fl)betheusualSobolevspaceonnwithnormWhenp=2,pisusuallyomitted.WeshalldenotetheusualinnerproductinL'(fl)orLa(O)'by','),andinL'(ofl)by't').Weshallusethesamenotationstoindicatethedualltiesbetw...     

An almost second order uniformly convergent scheme for a singularly perturbed initial value problem

In this paper a nonlinear singularly perturbed initial problem is considered. The behavior of the exact solution and its derivatives is analyzed, and this leads to the construction of a Shishkin-type mesh. On this mesh a hybrid difference scheme is proposed, which is a combination of the second order difference schemes on the fine mesh and the midpoint upwind scheme on the coarse mesh. It is proved that the scheme is almost second-order convergent, in the discrete maximum norm, independently of singular perturbation parameter. Numerical experiment supports these theoretical results.