Asymptotic expansions of finite element solutions to Robin problems in H 3 and their application in extrapolation cascadic multigrid method

For the Poisson equation with Robin boundary conditions,by using a few techniques such as orthogonal expansion(M-type),separation of the main part and the finite element projection,we prove for the first time that the asymptotic error expansions of bilinear finite element have the accuracy of O(h3)for u∈H3.Based on the obtained asymptotic error expansions for linear finite elements,extrapolation cascadic multigrid method(EXCMG)can be used to solve Robin problems effectively.Furthermore,by virtue of Richardson not only the accuracy of the approximation is improved,but also a posteriori error estimation is obtained.Finally,some numerical experiments that confirm the theoretical analysis are presented.     

EQ 1 rot nonconforming finite element method for nonlinear dual phase lagging heat conduction equations

EQ rot 1 nonconforming finite element approximation to a class of nonlinear dual phase lagging heat conduction equations is discussed for semi-discrete and fully-discrete schemes. By use of a special property, that is, the consistency error of this element is of order O(h2 ) one order higher than its interpolation error O(h), the superclose results of order O(h2 ) in broken H1 -norm are obtained. At the same time, the global superconvergence in broken H1 -norm is deduced by interpolation postprocessing technique. Moreover, the extrapolation result with order O(h4 ) is derived by constructing a new interpolation postprocessing operator and extrapolation scheme based on the known asymptotic expansion formulas of EQ rot 1 element. Finally, optimal error estimate is gained for a proposed fully-discrete scheme by different approaches from the previous literature.     

High accuracy analysis of the lowest order <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-Galerkin mixed finite element method for nonlinear sine-Gordon equations

The lowest order H¹-Galerkin mixed finite element method (for short MFEM) is proposed for a class of nonlinear sine-Gordon equations with the simplest bilinear rectangular element and zero order Raviart-Thomas element. Base on the interpolation operator instead of the traditional Ritz projection operator which is an indispensable tool in the traditional FEM analysis, together with mean-value technique and high accuracy analysis, the superclose properties of order O(h²)/O(h² + τ²) in H¹-norm and H(div;Ω)-norm are deduced for the semi-discrete and the fully-discrete schemes, where h, τ denote the mesh size and the time step, respectively, which improve the results in the previous literature.     

Supraconvergence and Supercloseness of a Scheme for Elliptic Equations on Nonuniform Grids

In this paper, we study the convergence of a finite difference scheme on nonuniform grids for the solution of second-order elliptic equations with mixed derivatives and variable coefficients in polygonal domains subjected to Dirichlet boundary conditions. We show that the scheme is equivalent to a fully discrete linear finite element approximation with quadrature. It exhibits the phenomenon of supraconvergence, more precisely, for s ∈ [1,2] order O(h ^s )-convergence of the finite difference solution, and its gradient is shown if the exact solution is in the Sobolev space H ^1+s (Ω). In the case of an equation with mixed derivatives in a domain containing oblique boundary sections, the convergence order is reduced to O(h ^3/2?ε) with ε > 0 if u ∈ H ³(Ω). The second-order accuracy of the finite difference gradient is in the finite element context nothing else than the supercloseness of the gradient. For s ∈ {1,2}, the given error estimates are strictly local.     

Convergence of H 1-Galerkin mixed finite element method for parabolic problems with reduced regularity on initial data

We study the convergence of H ¹-Galerkin mixed finite element method for parabolic problems in one space dimension. Both semi-discrete and fully discrete schemes are analyzed assuming less regularity on initial data. More precisely, for the spatially discrete scheme, error estimates of order \(\mathcal{O}\) (h ² t ^?1/2) for positive time are established assuming the initial function p ₀ ∈ H ²(Ω) ∩ H ₀ ¹ (Ω). Further, we use energy technique together with parabolic duality argument to derive error estimates of order \(\mathcal{O}\) (h ² t ^?1) when p ₀ is only in H ₀ ¹ (Ω). A discrete-in-time backward Euler method is analyzed and almost optimal order error bounds are established.     

A NONCONFORMING QUADRILATERAL FINITE ELEMENT APPROXIMATION TO NONLINEAR SCHRDINGER EQUATION

In this article, a nonconforming quadrilateral element(named modified quasiWilson element) is applied to solve the nonlinear schr¨odinger equation(NLSE). On the basis of a special character of this element, that is, its consistency error is of order O(h~3) for broken H1-norm on arbitrary quadrilateral meshes, which is two order higher than its interpolation error, the optimal order error estimate and superclose property are obtained. Moreover,the global superconvergence result is deduced with the help of interpolation postprocessing technique. Finally, some numerical results are provided to verify the theoretical analysis.     

The hull and geodetic numbers of orientations of graphs

For an oriented graph D, let I_D[u,v] denote the set of all vertices lying on a u-v geodesic or a v-u geodesic. For S⊆V(D), let I_D[S] denote the union of all I_D[u,v] for all u,v∈S. Let [S]_D denote the smallest convex set containing S. The geodetic number g(D) of an oriented graph D is the minimum cardinality of a set S with I_D[S]=V(D) and the hull number h(D) of an oriented graph D is the minimum cardinality of a set S with [S]_D=V(D). For a connected graph G, let O(G) be the set of all orientations of G, define g⁻(G)=min{g(D):D∈O(G)}, g⁺(G)=max{g(D):D∈O(G)}, h⁻(G)=min{h(D):D∈O(G)}, and h⁺(G)=max{h(D):D∈O(G)}. By the above definitions, h⁻(G)≤g⁻(G) and h⁺(G)≤g⁺(G). In the paper, we prove that g⁻(G)<h⁺(G) for a connected graph G of order at least 3, and for any nonnegative integers a and b, there exists a connected graph G such that g⁻(G)−h⁻(G)=a and g⁺(G)−h⁺(G)=b. These results answer a problem of Farrugia in [A. Farrugia, Orientable convexity, geodetic and hull numbers in graphs, Discrete Appl. Math. 148 (2005) 256-262].     

Unconditional superconvergent analysis of a new mixed finite element method for Ginzburg–Landau equation

In this article, unconditional superconvergent analysis of a linearized fully discrete mixed finite element method is presented for a class of Ginzburg–Landau equation based on the bilinear element and zero‐order Nédélec's element pair (Q₁₁/Q₀₁ × Q₁₀). First, a time‐discrete system is introduced to split the error into temporal error and spatial error, and the corresponding error estimates are deduced rigorously. Second, the unconditional superclose and optimal estimate of order O(h² + τ) for u in H¹‐norm and p = ?u in L²‐norm are derived respectively without the restrictions on the ratio between h and τ, where h is the subdivision parameter and τ, the time step. Third, the global superconvergent results are obtained by interpolated postprocessing technique. Finally, some numerical results are carried out to confirm the theoretical analysis.     

Extrapolation of the Nédélec element for the Maxwell equations by the mixed finite element method

In this paper, we use the integral-identity argument to obtain asymptotic error expansions for the mixed finite element approximation of the Maxwell equations on a rectangular mesh. The extrapolation method is applied to improve the accuracy of the approximation via an interpolation postprocessing technique. With the extrapolation, the approximation accuracy can be improved from O(h) to O(h ⁴) in the L ²-norm. Illustrative numerical results are given to demonstrate the higher order accuracy of the extrapolation method. This research was supported by the National Natural Science Foundation of China (No.10471103), Social Science Foundation of the Ministry of Education of China (06JA630047), Tianjin Natural Science Foundation (07JCYBJC14300).     

Extrapolation and superconvergence of the Steklov eigenvalue problem

On the basis of a transform lemma, an asymptotic expansion of the bilinear finite element is derived over graded meshes for the Steklov eigenvalue problem, such that the Richardson extrapolation can be applied to increase the accuracy of the approximation, from which the approximation of O(h ^3.5) is obtained. In addition, by means of the Rayleigh quotient acceleration technique and an interpolation postprocessing method, the superconvergence of the bilinear finite element is presented over graded meshes for the Steklov eigenvalue problem, and the approximation of O(h ³) is gained. Finally, numerical experiments are provided to demonstrate the theoretical results.     

Regularity of the solutions of the steady‐state Boussinesq equations with thermocapillarity effects on the surface of the liquid

In this paper we show that every variational solution of the steady‐state Boussinesq equations ( u , p, θ) with thermocapillarity effect on the surface of the liquid has the following regularity: u ∈ H²(Ω)², p ∈ H¹(Ω), θ ∈ H²(Ω) under appropriate hypotheses on the angles of the ‘2‐D’ container (a cross‐section of the 3‐D container in fact) and of the horizontal surface of the liquid with the inner surface of the container. The difficulty comes from the boundary condition on the surface of the liquid (e.g. water) which modelizes the thermocapillarity effect on the surface of the liquid (equation (68.10) of Levich [7]). More precisely we will show that u ∈ P²₂(Ω)² and that θ ∈ P²₂(Ω), where P²₂(Ω) denotes the usual Kondratiev space. This result will be used in a forthcoming paper to prove convergence results for finite element methods intended to compute approximations of a non‐singular solution [1] of this problem. Copyright © 1999 John Wiley & Sons, Ltd.     

A second order nonconforming rectangular finite element method for approximating Maxwell’s equations

The main objective of this paper is to present a new rectangular nonconforming finite element scheme with the second order convergence behavior for approximation of Maxwell’s equations.Then the corresponding optimal error estimates are derived.The difficulty in construction of this finite element scheme is how to choose a compatible pair of degrees of freedom and shape function space so as to make the consistency error due to the nonconformity of the element being of order O(h 3 ) ,properly one order higher than that of its interpolation error O(h 2 ) in the broken energy norm,where h is the subdivision parameter tending to zero.     

A low order anisotropic nonconforming characteristic finite element method for a convection-dominated transport problem

A low order anisotropic nonconforming rectangular finite element method for the convection-diffusion problem with a modified characteristic finite element scheme is studied in this paper. The O(h²) order error estimate in L²-norm with respect to the space, one order higher than the expanded characteristic-mixed finite element scheme with order O(h), and the same as the conforming case for a modified characteristic finite element scheme under regular meshes, is obtained by use of some distinct properties of the interpolation operator and the mean value technique, instead of the so-called elliptic projection, which is an indispensable tool in the convergence analysis of the previous literature. Lastly, some numerical results of the element are provided to verify our theoretical analysis.     

On the error bounds of nonconforming finite elements

We prove that the error estimates of a large class of nonconforming finite elements are dominated by their approximation errors, which means that the well-known Cea's lemma is still valid for these nonconforming finite element methods. Furthermore, we derive the error estimates in both energy and L2 norms under the regularity assumption u ∈ H1+s(Ω) with any s 0. The extensions to other related problems are possible.     

Blow-up of critical norms for the 3-D Navier-Stokes equations

Let u =(uh, u3) be a smooth solution of the 3-D Navier-Stokes equations in R3× [0, T). It was proved that if u3 ∈ L∞(0, T;˙B-1+3/p p,q(R3)) for 3 p, q ∞ and uh∈ L∞(0, T; BMO-1(R3)) with uh(T) ∈ VMO-1(R3), then u can be extended beyond T. This result generalizes the recent result proved by Gallagher et al.(2016), which requires u ∈ L∞(0, T;˙B-1+3/pp,q(R3)). Our proof is based on a new interior regularity criterion in terms of one velocity component, which is independent of interest.     

On the computation of the hull number of a graph

Let G be a graph. If u,v∈V(G), a u-vshortest path of G is a path linking u and v with minimum number of edges. The closed interval I[u,v] consists of all vertices lying in some u-v shortest path of G. For S⊆V(G), the set I[S] is the union of all sets I[u,v] for u,v∈S. We say that S is a convex set if I[S]=S. The convex hull of S, denoted I_h[S], is the smallest convex set containing S. A set S is a hull set of G if I_h[S]=V(G). The cardinality of a minimum hull set of G is the hull number of G, denoted by hn(G). In this work we prove that deciding whether hn(G)≤k is NP-complete.We also present polynomial-time algorithms for computing hn(G) when G is a unit interval graph, a cograph or a split graph.     

Averaging of gradient in the space of linear triangular and bilinear rectangular finite elements

An averaging method for the second-order approximation of the values of the gradient of an arbitrary smooth function u = u(x ₁, x ₂) at the vertices of a regular triangulation T _h composed both of rectangles and triangles is presented. The method assumes that only the interpolant Π _h [u] of u in the finite element space of the linear triangular and bilinear rectangular finite elements from T _h is known. A complete analysis of this method is an extension of the complete analysis concerning the finite element spaces of linear triangular elements from [Dalík J., Averaging of directional derivatives in vertices of nonobtuse regular triangulations, Numer. Math., 2010, 116(4), 619–644]. The second-order approximation of the gradient is extended from the vertices to the whole domain and applied to the a posteriori error estimates of the finite element solutions of the planar elliptic boundary-value problems of second order. Numerical illustrations of the accuracy of the averaging method and of the quality of the a posteriori error estimates are also presented.     

Crank-Nicolson split least-squares procedure for nonlocal reactive flow in porous media

In this paper, we introduce a Crank-Nicolson split least-squares Galerkin finite element procedure for parabolic integro-differential equations, arising in the modeling of nonlocal reactive flows in porous media. By selecting the least-squares functional properly, the procedure can be split into two independent sub-procedures, one of which is for the primitive unknown and the other is for the flux. By carefully choosing projections, we get optimal order H ¹(Ω) and L ²(Ω) norm error estimates for u and sub-optimal (L ²(Ω)) ^d norm error estimate for σ with second-order accuracy in time increment. The numerical examples are given to testify the efficiency of the introduced scheme.     

Superconvergence analysis of a linearized three‐step backward differential formula finite element method for nonlinear Sobolev equation

A linearized three‐step backward differential formula (BDF) Galerkin finite element method (FEM) is developed for nonlinear Sobolev equation with bilinear element. Temporal error and spatial error are discussed through introducing a time‐discrete system. Solutions of the time‐discrete system are bounded in H²‐norm by the temporal error. Superconvergence results of order O(h² + τ³) in H¹‐norm for the original variable are deduced based on the spatial error. Some new tricks are utilized to get higher order of the temporal error and the spatial error. At last, two numerical examples are provided to support the theoretical analysis. Here, h is the subdivision parameter, and τ is the time step.     

On the Attainable Order of Collocation Methods for Delay Differential Equations with Proportional Delay

To analyze the attainable order of m-stage implicit (collocation-based) Runge-Kutta methods for the delay differential equation (DDE) y′(t) = by(qt), 0 < q ≤ 1 with y(0) = 1, and the delay Volterra integral equation (DVIE) y(t) = 1 + $\tfrac{b}{q}\int {_0^{qt} }$ y(s) ds with proportional delay qt, 0 < q ≤ 1, our particular interest lies in the approximations (and their orders) at the first mesh point t = h for the collocation solution v(t) of the DDE and the iterated collocation solution u _it(t) of the DVIE to the solution y(t). Recently, H. Brunner proposed the following open problem: “For m ≤ 3, do there exist collocation points c _i = c _i(q), i = 1, 2,..., m in [0,1] such that the rational approximant v(h)is the (m, m)-Padé approximant to y(h)? If these exist, then |v(h) ? y(h)| = O(h ^2m+1) but what is the collocation polynomial M _m(t; q) = K Π _i=1 ^m (t ? c _i) of v(th), t ∈ [0, 1]?” In this paper, we solve this question affirmatively, and give the related results between the collocation solution v(t) of the DDE and the iterated collocation solution u _it(t) of the DVIE. We also answer to Brunner's second open question in the case that one collocation point is fixed at the right end point of the interval.