Stokes型积分-微分方程Q_2-P_1元的超收敛分析

本文研究Q2-P1混合元对Stokes型积分-微分方程的有限元方法.利用积分恒等式技巧给出了关于流体速度u和压力p的误差估计,特别是在压力p的误差中去掉了影响解的稳定性的1因子t-2,改善了以往文献的结果.同时,通过构造适当的插值后处理算子得到了整体超收敛结果.     

The optimal convergence of the h–p version of the boundary element method with quasiuniform meshes for elliptic problems on polygonal domains

In the framework of the Jacobi-weighted Besov spaces, we analyze the lower and upper bounds of errors in the h–p version of boundary element solutions on quasiuniform meshes for elliptic problems on polygons. Both lower bound and upper bound are optimal in h and p, and they are of the same order. The optimal convergence of the h–p version of boundary element method with quasiuniform meshes is proved, which includes the optimal rates for h version with quasiuniform meshes and the p version with quasiuniform degrees as two special cases. Dedicated to Professor Charles Micchelli on the occasion of his sixtieth birthday Mathematics subject classification (2000) 65N38. Benqi Guo: The work of this author was supported by NSERC of Canada under Grant OGP0046726 and was complete during visiting Newton Institute for Mathematical Sciences, Cambridge University for participating in special program “Computational Challenges in PDEs” in 2003. Norbert Heuer: This author is supported by Fondecyt project No. 1010220 and by the FONDAP Program (Chile) on Numerical Analysis. Current address: Mathematical Sciences, Brunel University, Uxbridge, U.K.     

Superconvergence of H 1-Galerkin Mixed Finite Element Methods for Second-Order Elliptic Equations

We derive superconvergence result for H ¹-Galerkin mixed finite element method for second-order elliptic equations over rectangular partitions. Compared to standard mixed finite element procedure, the method is not subject to the Ladyzhenskaya–Bab?ska–Brezzi (LBB) condition and the approximating finite element spaces are allowed to be of different polynomial degrees. Superconvergence estimate of order 𝒪(h ^k+3), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas elements, is established for the flux via a postprocessing technique.     

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.     

Superconvergence of H 1-Galerkin mixed finite element methods for parabolic problems

In this article, we study the semidiscrete H ¹-Galerkin mixed finite element method for parabolic problems over rectangular partitions. The well-known optimal order error estimate in the L ²-norm for the flux is of order 𝒪(h ^k+1) (SIAM J. Numer. Anal. 35 (2), (1998), pp. 712–727), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas element. We derive a superconvergence estimate of order 𝒪(h ^k+3) between the H ¹-Galerkin mixed finite element approximation and an appropriately defined local projection of the flux variable when k ≥ 1. A the new approximate solution for the flux with superconvergence of order 𝒪(h ^k+3) is realized via a postprocessing technique using local projection methods.     

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).     

An hp‐adaptive Newton‐Galerkin finite element procedure for semilinear boundary value problems

In this paper, we develop an hp‐adaptive procedure for the numerical solution of general, semilinear elliptic boundary value problems in 1d, with possible singular perturbations. Our approach combines both a prediction‐type adaptive Newton method and an hp‐version adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully hp‐adaptive Newton–Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for various examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite volume method based on stabilized finite elements for the nonstationary Navier–Stokes problem

A finite volume method based on stabilized finite element for the two‐dimensional nonstationary Navier–Stokes equations is investigated in this work. As in stabilized finite element method, macroelement condition is introduced for constructing the local stabilized formulation of the nonstationary Navier–Stokes equations. Moreover, for P₁ ? P₀ element, the H¹ error estimate of optimal order for finite volume solution (u_h,p_h) is analyzed. And, a uniform H¹ error estimate of optimal order for finite volume solution (u_h, p_h) is also obtained if the uniqueness condition is satisfied. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On approximating the solution of the non-stationary Stokes equations using the cell discretization algorithm

The cell discretization algorithm, a nonconforming extension of the finite element method, is used to obtain approximations to the velocity and pressure satisfying the nonstationary Stokes equations. Error estimates show convergence of the approximations. An implementation using polynomial bases is described that permits the use of the continuous approximations of the h–p finite element method and exactly satisfies the solenoidal requirement. We express the error estimates in terms of the diameter h of a cell and the degree p of the approximation on each cell. Results of an experiment with p10 are presented that confirm the theoretical estimates.     

On the construction of stable curvilinear pp version elements for mixed formulations of elasticity and Stokes flow

Summary. The use of mixed finite element methods is well-established in the numerical approximation of the problem of nearly incompressible elasticity, and its limit, Stokes flow. The question of stability over curved elements for such methods is of particular significance in the p version, where, since the element size remains fixed, exact representation of the curved boundary by (large) elements is often used. We identify a mixed element which we show to be optimally stable in both p and h refinement over curvilinear meshes. We prove optimal p version (up to ) and h version (p = 2, 3) convergence for our element, and illustrate its optimality through numerical experiments. Received August 25, 1998 / Revised version received February 16, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Two-Grid hp-Version DGFEMs for Strongly Monotone Second-Order Quasilinear Elliptic PDEs

In this article we develop the a priori error analysis of so-called two-grid hp-version discontinuous Galerkin finite element methods for the numerical approximation of strongly monotone second-order quasilinear partial differential equations. In this setting, the fully nonlinear problem is first approximated on a coarse finite element space V(𝒯_H, P ). The resulting ‘coarse’ numerical solution is then exploited to provide the necessary data needed to linearize the underlying discretization on the finer space V(𝒯_h, p ); thereby, only a linear system of equations is solved on the richer space V(𝒯_h, p ). Numerical experiments confirming the theoretical results are presented. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Error estimates for the combinedh andp versions of the finite element method

Summary In theh-version of the finite element method, convergence is achieved by refining the mesh while keeping the degree of the elements fixed. On the other hand, thep-version keeps the mesh fixed and increases the degree of the elements. In this paper, we prove estimates showing the simultaneous dependence of the order of approximation on both the element degrees and the mesh. In addition, it is shown that a proper design of the mesh and distribution of element degrees lead to a better than polynomial rate of convergence with respect to the number of degrees of freedom, even in the presence of corner singularities. Numerical results comparing theh-version,p-version, and combinedh-p-version for a one dimensional problem are presented.     

An iterative and parallel solver based on domain decomposition for the h-p version of the finite element method

In this paper, we propose a new interative and parallel solver, based on domain decomposition, for the h-p version of the finite element method in two dimensions. It improves our previous work in two aspects: (1) A subdomain may contain several super-elements of the coarse mesh, thus can be of arbitrary shape and size. This makes the solver more efficient and more flexible in computational practice. (2) The p-version components (i.e., the high order side and internal modes) in every element are treated separately, which results in better parallelism.     

Congruences for the class numbers of real cyclic sextic number fields

LetK ₆ be a real cyclic sextic number field, andK ₂,K ₃ its quadratic and cubic subfield. Leth(L) denote the ideal class number of fieldL. Seven congruences forh ^- =h (K ₆)/(h(K ₂)h(K ₃)) are obtained. In particular, when the conductorf ₆ ofK ₆ is a primep, , whereC is an explicitly given constant, andB _n is the Bernoulli number. These results on real cyclic sextic fields are an extension of the results on quadratic and cyclic quartic fields. Project supported by the National Natural Science Foundation of China (Grant No. 19771052).     

Applications of discrete maximal <Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript> regularity for finite element operators

In this paper, we present applications of discrete maximal L _p regularity for finite element operators. More precisely, we show error estimates of order h ² for linear and certain semilinear problems in various L _p (Ω)-norms. Discrete maximal regularity allows us to prove error estimates in a very easy and efficient way. Moreover, we also develop interpolation theory for (fractional powers of) finite element operators and extend the results on discrete maximal L _p regularity formerly proved by the author. The author was supported by the DFG-Graduiertenkolleg 853.     

Two-parameter estimates for joint spectral projections on complex spheres

We prove sharp two-parameter estimates for the L ^p -L ² norm, 1 ≤ p ≤ 2, of the joint spectral projectors associated to the Laplace–Beltrami operator and to the Kohn Laplacian on the unit sphere S ^2n-1 in . Then, by using the notion of contraction of Lie groups, we deduce the estimates recently obtained by H. Koch and F. Ricci for joint spectral projections on the reduced Heisenberg group h ¹.      

Additive Schwarz algorithms for the p version of the Galerkin boundary element method

Summary. We study some additive Schwarz algorithms for the version Galerkin boundary element method applied to some weakly singular and hypersingular integral equations of the first kind. Both non-overlapping and overlapping methods are considered. We prove that the condition numbers of the additive Schwarz operators grow at most as independently of h, where p is the degree of the polynomials used in the Galerkin boundary element schemes and h is the mesh size. Thus we show that additive Schwarz methods, which were originally designed for finite element discretisation of differential equations, are also efficient preconditioners for some boundary integral operators, which are non-local operators. Received June 15, 1997 / Revised version received July 7, 1998 / Published online February 17, 2000     

The approximation theory for the p-version finite element method and application to non-linear elliptic PDEs

Approximation theoretic results are obtained for approximation using continuous piecewise polynomials of degree p on meshes of triangular and quadrilateral elements. Estimates for the rate of convergence in Sobolev spaces , are given. The results are applied to estimate the rate of convergence when the p-version finite element method is used to approximate the -Laplacian. It is shown that the rate of convergence of the p-version is always at least that of the h-version (measured in terms of number of degrees of freedom used). If the solution is very smooth then the p-version attains an exponential rate of convergence. If the solution has certain types of singularity, the rate of convergence of the p-version is twice that of the h-version. The analysis generalises the work of Babuska and others to the case . In addition, the approximation theoretic results find immediate application for some types of spectral and spectral element methods. Received August 2, 1995 / Revised version received January 26, 1998