DOMAIN DECOMPOSITION METHOD FOR PARABOLIC EQUATION ON GENERAL DOMAIN

The nonoverlapping domain decomposition method for parabolic partialdifferential equation on general domain is considered. A kind of domain decompositionthat uses the finite element procedure is given. The problem over the domains can beimplemented on parallel computer. Convergence analysis is also presented.

FULL DISCRETE NONLINEAR GALERKIN METHOD FOR THE NAVIER-STOKES EQUATIONS

This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.

CHARACTERISTICS DIFFERENCE METHOD FOR INCOMPRESSIBLE MISCIBLE FLOW IN POROUS MEDIA

Abstract Two phase,incompressible miscible flow in porous media is governed by a system of nonlinear partial differential equations.The pressure equation,which is elliptic in appearance,is discretized by a standard five points difference method.The concentration equation is treated by an implicit finite difference method that applies a form of the method of characteristics to the transport terms.A class of biquadratic interpolation is introduced for the method of characteristics.Convergence rate is proved to be O(Δt+h 2).

一类Cauchy型奇异积分方程的线性有限元逼近

考虑了一类Ｃａｕｃｈｙ型奇异积分方程的有限元解，证明了逼近解的误差估计．

CONVERGENCE OF A MIXED FINITE ELEMENT FOR THE STOKES PROBLEM ON ANISOTROPIC MESHES

The main aim of this paper is to study the convergence properties of a low order mixed finite element for the Stokes problem under anisotropic meshes. We discuss the anisotropic convergence and superconvergence independent of the aspect ratio. Without the shape regularity assumption and inverse assumption on the meshes, the optimal error estimates and natural superconvergence at central points are obtained. The global superconvergence for the gradient of the velocity and the pressure is derived with the aid of a suitable postprocessing method. Furthermore, we develop a simple method to obtain the superclose properties which improves the results of the previous works .

Analysis of extrapolation cascadic multigrid method(EXCMG)

Based on an asymptotic expansion of finite element,a new extrapolation formula and extrapolation cascadic multigrid method(EXCMG)are proposed,in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid.In the case of triple grids,the error of the new initial value is analyzed in detail.A larger scale computation is completed in PC.

Finite element formulation based on proper orthogonal decomposition for parabolic equations

A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.

Finite difference scheme based on proper orthogonal decomposition for the nonstationary Navier-Stokes equations

The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.

Interpolation theory of anisotropic finite elements and applications

Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.

The order-preserving convergence for spectral approximation of self-adjoint completely continuous operators

This paper discusses the order-preserving convergence for spectral approximation of the self-adjoint completely continuous operator T.Under the condition that the approximate operator Th converges to T in norm,it is proven that the k-th eigenvalue of Th converges to the k-th eigenvalue of T.(We sorted the positive eigenvalues in decreasing order and negative eigenvalues in increasing order.) Then we apply this result to conforming elements,nonconforming elements and mixed elements of self-adjoint elliptic differential operators eigenvalue problems,and prove that the k-th approximate eigenvalue obtained by these methods converges to the k-th exact eigenvalue.

Analysis of search-extension method for finding multiple solutions of nonlinear problem

For numerical computations of multiple solutions of the nonlinear elliptic problemΔu+ f（u）=0 inΩ, u=0 onΓ, a search-extension method （SEM） was proposed and systematically studied by the authors. This paper shall complete its theoretical analysis. It is assumed that the nonlinearity is non-convex and its solution is isolated, under some conditions the corresponding linearized problem has a unique solution. By use of the compactness of the solution family and the contradiction argument, in general conditions, the high order regularity of the solution u∈H1+α,α】0 is proved. Assume that some initial value searched by suitably many eigenbases is already fallen into the neighborhood of the isolated solution, then the optimal error estimates of its nonlinear finite element approximation are shown by the duality argument and continuation method.

Cascadic multigrid methods for parabolic problems

In this paper,we consider the cascadic multigrid method for a parabolic type equation.Backward Euler approximation in time and linear finite element approximation in space are employed.A stability result is established under some conditions on the smoother.Using new and sharper estimates for the smoothers that reflect the precise dependence on the time step and the spatial mesh parameter,these conditions are verified for a number of popular smoothers.Optimal error bound sare derived for both smooth and non-smooth data.Iteration strategies guaranteeing both the optimal accuracy and the optimal complexity are presented.

Residual-based a posteriori error estimates of nonconforming finite element method for elliptic problems with Dirac delta source terms

Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.

Generalized Rayleigh quotient and finite element two-grid discretization schemes

This study discusses generalized Rayleigh quotient and high efficiency finite element discretization schemes. Some results are as follows: 1) Rayleigh quotient accelerate technique is extended to nonselfadjoint problems. Generalized Rayleigh quotients of operator form and weak form are defined and the basic relationship between approximate eigenfunction and its generalized Rayleigh quotient is established. 2) New error estimates are obtained by replacing the ascent of exact eigenvalue with the ascent of finite element approximate eigenvalue. 3) Based on the work of Xu Jinchao and Zhou Aihui, finite element two-grid discretization schemes are established to solve nonselfadjoint elliptic differential operator eigenvalue problems and these schemes are used in both conforming finite element and non-conforming finite element. Besides, the efficiency of the schemes is proved by both theoretical analysis and numerical experiments. 4) Iterated Galerkin method, interpolated correction method and gradient recovery for selfadjoint elliptic differential operator eigenvalue problems are extended to nonselfadjoint elliptic differential operator eigenvalue problems.

Superconvergence of tricubic block finite elements

In this paper, we first introduce interpolation operator of projection type in three dimen- sions, from which we derive weak estimates for tricubic block finite elements. Then using the estimate for the W 2, 1-seminorm of the discrete derivative Green's function and the weak estimates, we show that the tricubic block finite element solution uh and the tricubic interpolant of projection type Πh3u have superclose gradient in the pointwise sense of the L∞-norm. Finally, this supercloseness is applied to superconvergence analysis, and the global superconvergence of the finite element approximation is derived.

Further results on ultraconvergence derivative recovery for odd-order rectangular finite elements

For rectangular finite element, we give a superconvergence method by SPR technique based on the generalization of a new ultraconvergence record and the sharp Green function estimates, by which we prove that the derivative has ultra-convergence of order O(hk+3) (k 3 being odd) and displacement has order of O(hk+4) (k 4 being even) at the locally symmetry points.

Two lower order nonconforming quadrilateral elements for the Reissner-Mindlin plate

This paper generalizes two nonconforming rectangular elements of the Reissner-Mindlin plate to the quadrilateral mesh. The first quadrilateral element uses the usual conforming bilinear element to approximate both components of the rotation, and the modified nonconforming rotated Q1 element enriched with the intersected term on each element to approximate the displacement, whereas the second one uses the enriched modified nonconforming rotated Q1 element to approximate both the rotation and the displacement. Both elements employ a more complicated shear force space to overcome the shear force locking, which will be described in detail in the introduction. We prove that both methods converge at optimal rates uniformly in the plate thickness t and the mesh distortion parameter in both the H1-and the L2-norms, and consequently they are locking free.

Uniqueness and stability in an inverse problem for a Poisson's equation

Consider the Poisson’s equation（?）″（x）=-ev-（?）+e（?）-v-N（x）with the Diriehlet boundary data,and we mainly investigate the inverse problem of determining the unknown function N（x）from a parameter function family.Some uniqueness and stability results in the inverse problem are obtained.