Parallel Galerkin domain decomposition procedures for wave equation

Parallel Galerkin domain decomposition procedures for wave equation are given. These procedures use implicit method in the sub-domains and simple explicit flux calculation on the inter-boundaries of sub-domains by integral mean method or extrapolation method. Thus, the parallelism can be achieved by these procedures. The explicit nature of the flux prediction induces a time step constraint that is necessary to preserve the stability. L²-norm error estimates are derived for these procedures. Experimental results are presented to confirm the theoretical results.     

Parallel domain decomposition procedures of improved D-D type for parabolic problems

Two parallel domain decomposition procedures for solving initial-boundary value problems of parabolic partial differential equations are proposed. One is the extended D-D type algorithm, which extends the explicit/implicit conservative Galerkin domain decomposition procedures, given in [5], from a rectangle domain and its decomposition that consisted of a stripe of sub-rectangles into a general domain and its general decomposition with a net-like structure. An almost optimal error estimate, without the factor H^−1/2 given in Dawson-Dupont’s error estimate, is proved. Another is the parallel domain decomposition algorithm of improved D-D type, in which an additional term is introduced to produce an approximation of an optimal error accuracy in L²-norm.     

A posteriori error estimates for mixed finite element approximations of elliptic problems

We derive residual based a posteriori error estimates of the flux in L ²-norm for a general class of mixed methods for elliptic problems. The estimate is applicable to standard mixed methods such as the Raviart–Thomas–Nedelec and Brezzi–Douglas–Marini elements, as well as stabilized methods such as the Galerkin-Least squares method. The element residual in the estimate employs an elementwise computable postprocessed approximation of the displacement which gives optimal order.     

Nonconforming mixed finite element approximation to the stationary Navier-Stokes equations on anisotropic meshes

The main aim of this paper is to study the error estimates of a rectangular nonconforming finite element for the stationary Navier-Stokes equations under anisotropic meshes. That is, the nonconforming rectangular element is taken as approximation space for the velocity and the piecewise constant element for the pressure. The convergence analysis is presented and the optimal error estimates both in a broken H¹-norm for the velocity and in an L²-norm for the pressure are derived on anisotropic meshes.     

Convergence and superconvergence analysis of an anisotropic nonconforming finite element methods for singularly perturbed reaction-diffusion problems

The numerical approximation by a lower order anisotropic nonconforming finite element on appropriately graded meshes are considered for solving singular perturbation problems. The quasi-optimal order error estimates are proved in the ε-weighted H¹-norm valid uniformly, up to a logarithmic factor, in the singular perturbation parameter. By using the interpolation postprocessing technique, the global superconvergent error estimates in ε-weighted H¹-norm are obtained. Numerical experiments are given to demonstrate validity of 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.     

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.     

An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method

An approximation scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using two methods. Standard mixed finite element is used for the Darcy velocity equation. A characteristics-mixed finite element method is presented for the concentration equation. Characteristic approximation is applied to handle the convection part of the concentration equation, and a lowest-order mixed finite element spatial approximation is adopted to deal with the diffusion part. Thus, the scalar unknown concentration and the diffusive flux can be approximated simultaneously. In order to derive the optimal L²$L^2$-norm error estimates, a post-processing step is included in the approximation to the scalar unknown concentration. This scheme conserves mass globally; in fact, on the discrete level, fluid is transported along the approximate characteristics. Numerical experiments are presented finally to validate the theoretical analysis.     

A coupling method based on new MFE and FE for fourth-order parabolic equation

In this article, a coupling method of new mixed finite element (MFE) and finite element (FE) is proposed and analyzed for fourth-order parabolic partial differential equation. First, the fourth-order parabolic equation is split into the coupled system of second-order equations. Then, an equation is solved by finite element method, the other equation is approximated by the new mixed finite element method, whose flux belongs to the square integrable space replacing the classical H(div;Ω) space. The stability for fully discrete scheme is derived, and both semi-discrete and fully discrete error estimates are obtained. Moreover, the optimal a priori error estimates in L ² and H ¹-norm for both the scalar unknown u and the diffusion term γ and a priori error estimate in (L ²)²-norm for its flux σ are derived. Finally, some numerical results are provided to validate our theoretical analysis.     

A two-grid method for finite volume element approximations of second-order nonlinear hyperbolic equations

The two-grid method is studied for solving a two-dimensional second-order nonlinear hyperbolic equation using finite volume element method. The method is based on two different finite element spaces defined on one coarse grid with grid size H and one fine grid with grid size h, respectively. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. A prior error estimate in the H¹-norm is proved to be O(h+H³|lnH|) for the two-grid semidiscrete finite volume element method. With these proposed techniques, solving such a large class of second-order nonlinear hyperbolic equations will not be much more difficult than solving one single linearized equation. Finally, a numerical example is presented to validate the usefulness and efficiency of the method.     

A finite difference analysis of a streamline diffusion method on a Shishkin mesh

We consider streamline diffusion finite element methods applied to a singularly perturbed convection–diffusion two‐point boundary value problem whose solution has a single boundary layer. To analyse the convergence of these methods, we rewrite them as finite difference schemes. We first consider arbitrary meshes, then, in analysing the scheme on a Shishkin mesh, we consider two formulations on the fine part of the mesh: the usual streamline diffusion upwinding and the standard Galerkin method. The error estimates are given in the discrete L _∞ norm; in particular we give the first analysis that shows precisely how the error depends on the user-chosen parameter _τ0 specifying the mesh. When _τ0 is too small, the error becomes O(1), but for _τ0 above a certain threshold value, the error is small and increases either linearly or quadratically as a function of . Numerical tests support our theoretical results. This revised version was published online in August 2006 with corrections to the Cover Date.     

A new a posteriori error estimate for convection–reaction–diffusion problems

A new a posteriori error estimate is derived for the stationary convection–reaction–diffusion equation. In order to estimate the approximation error in the usual energy norm, the underlying bilinear form is decomposed into a computable integral and two other terms which can be estimated from above using elementary tools of functional analysis. Two auxiliary parameter-functions are introduced to construct such a splitting and tune the resulting bound. If these functions are chosen in an optimal way, the exact energy norm of the error is recovered, which proves that the estimate is sharp. The presented methodology is completely independent of the numerical technique used to compute the approximate solution. In particular, it is applicable to approximations which fail to satisfy the Galerkin orthogonality, e.g. due to an inconsistent stabilization, flux limiting, low-order quadrature rules, round-off and iteration errors, etc. Moreover, the only constant that appears in the proposed error estimate is global and stems from the Friedrichs–Poincaré inequality. Numerical experiments illustrate the potential of the proposed error estimation technique.     

A locking-free nonconforming triangular element for planar elasticity with pure traction boundary condition

A new nonconforming triangular element for the equations of planar linear elasticity with pure traction boundary conditions is considered. By virtue of construction of the element, the discrete version of Korn’s second inequality is directly proved to be valid. Convergence rate of the finite element methods is uniformly optimal with respect to λ. Error estimates in the energy norm and L²-norm are O(h²) and O(h³), respectively.     

A posteriori error estimates for hp finite element solutions of convex optimal control problems

In this paper, we present a posteriori error analysis for hp finite element approximation of convex optimal control problems. We derive a new quasi-interpolation operator of Clément type and a new quasi-interpolation operator of Scott-Zhang type that preserves homogeneous boundary condition. The Scott-Zhang type quasi-interpolation is suitable for an application in bounding the errors in L²-norm. Then hp a posteriori error estimators are obtained for the coupled state and control approximations. Such estimators can be used to construct reliable adaptive finite elements for the control problems.     

A posteriori error estimation of residual type for anisotropic diffusion-convection-reaction problems

This paper presents a robust a posteriori residual error estimator for diffusion-convection-reaction problems with anisotropic diffusion, approximated by a SUPG finite element method on isotropic or anisotropic meshes in R^d, d=2 or 3. The equivalence between the energy norm of the error and the residual error estimator is proved. Numerical tests confirm the theoretical results.     

Superconvergence analysis of approximate boundary-flux calculations

Summary Certain projection post-processing techniques have been proposed for computing the boundary flux for two-dimensional problems (e.g., see Carey, et al. [5]). In a series of numerical experiments on elliptic problems they observed that these post-processing formulas for approximate fluxes were almost (O(h ²)-accurate for linear triangular elements. In this paper we prove that the computed boundary flux isO(h ² ln 1/h)-accurate in the maximum norm for the partial method of [5]. If the solutionuH ³() then the boundary flux error isO(h ^3/2) in theL ²-norm.     

Error estimates of the lumped mass finite element method for semilinear elliptic problems

We are concerned with the semilinear elliptic problems. We first investigate the L²-error estimate for the lumped mass finite element method. We then use the cascadic multigrid method to solve the corresponding discrete problem. On the basis of the finite element error estimates, we prove the optimality of the proposed multigrid method. We also report some numerical results to support the theory.     

A remark on least-squares mixed element methods for reaction–diffusion problems

In this paper, we propose a least-squares mixed element procedure for a reaction–diffusion problem based on the first-order system. By selecting the least-squares functional properly, the resulting procedure can be split into two independent symmetric positive definite schemes, one of which is for the unknown variable and the other of which is for the unknown flux variable, which lead to the optimal order H¹(Ω)$H^1(\Omega )$ and L²(Ω)$L^2(\Omega )$ norm error estimates for the primal unknown and optimal H(div;Ω)$H(\mathrm{div};\Omega )$ norm error estimate for the unknown flux. Finally, we give some numerical examples.     

The 2D Boussinesq equations with vertical viscosity and vertical diffusivity

This paper aims at the global regularity of classical solutions to the 2D Boussinesq equations with vertical dissipation and vertical thermal diffusion. We prove that the Lr-norm of the vertical velocity v for any 1<r<∞ is globally bounded and that the L^∞-norm of v controls any possible breakdown of classical solutions. In addition, we show that an extra thermal diffusion given by the fractional Laplace δ(−Δ) for δ>0 would guarantee the global regularity of classical solutions.     

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain are given. These procedures use implicit Galerkin method in the subdomains and simple explicit flux calculation on the interdomain boundaries by integral mean method or extrapolation method to predict the inner‐boundary conditions. Thus, the parallelism can be achieved by these procedures. These procedures are conservative both in the subdomains and across interboundaries. The explicit nature of the flux prediction induces a time‐step limitation that is necessary to preserve stability, but this constraint is less severe than that for a fully explicit method. L²‐norm error estimates are derived for these procedures. Compared with the work of Dawson and Dupont [Math Comp 58 (1992), 21–35], these L²‐norm error estimates avoid the loss of H^?1/2 factor. Experimental results are presented to confirm the theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009