THE OPTIMAL ORDER ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS TO HYPERBOLIC PROBLEMS

In this paper, the linear finite element approximation to the positive and symmetric,linear hyperbolic systems is analyzed and an O(h^2) order error estimate is established under the conditions of strongly regular triangulation and the H^3-regularity for the exact solutions. The convergence analysis is based on some superclose estimates derived in this paper. Our method and result here are also applicable to general hyperbolic problems.Finally, we discuss the linearized shallow water system of equations.     

THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

A NEW APPROXIMATE INERTIAL MANIFOLD AND ASSOCIATED ALGORITHM

In this article the authors propose a new approximate inertial manifold(AIM) to the Navier-Stokes equations. The solutions are in the neighborhoods of this AIM with thickness δ=o(h2k 1-ε). The article aims to investigate a two grids finite element approximation based on it and give error estimates of the approximate solution where (h, h*) and (k, m) are coarse and fine meshes and degree of finite element subspaces, respectively. These results are much better than Standard Galerkin(SG) and nonlinear Galerkin (NG) methods. For example, for 2D NS eqs and linear element, let uh,uh,u* be the SG, NG and their approximate solutions respectively, then and h*≈h2 for NG, h*≈h3/2 for theirs.     

积分微分方程的一个高精度导数恢复公式

A highly accurate derivative recovery formula is presented for the korder finite element approximations to integro-differential equations. This formula possesses O(h^k l) order of superconvergence on the whole domain in L∞ norm and O(h^2k) order of ultraconvergence at the mesh points, respectively, and the same type of results are demonstrated for the semi-discrete finite element approximates solutions of Sobolev equations.     

Superconvergence of Continuous Finite Elements with Interpolated Coeffcients for Initial Value Problems of Nonlinear Ordinary Differential Equation

In this paper, n-degree continuous finite element method with interpolated coefficients for nonlinear initial value problem of ordinary differential equation is introduced and analyzed. An optimal superconvergence u-uh = O(hn+2), n ≥ 2, at (n + 1)-order Lobatto points in each element respectively is proved. Finally the theoretical results are tested by a numerical example.     

THE DUAL MIXED METHOD FOR AN UNILATERAL PROBLEM

In this paper, the dual mixed method for an unilateral problem, which is the simplified modelling of scalar function for the friction-free contact problem, is considered. The dual mixed problem is introduced, the existence and uniqeness of the solution of the problem are presented, and error bounds O(h^3/4 ) and O(h^3/2 ) are obtained for the dual mixed finite element approximations of Raviart-Thomas elements for k= 0 and k= 1 respectively.     

SUPERCONVERGENCE OF LEAST-SQUARES MIXED FINITE ELEMENT FOR SECOND-ORDER ELLIPTIC PROBLEMS

In this paper the least-squares mixed finite element is considered for solving secondorder elliptic problems in two dimensional domains. The primary solution u and the flux er are approximated using finite element spaces consisting of piecewise polynomials of degree k and r respectively. Based on interpolation operators and an auxiliary projection,superconvergent H^1-error estimates of both the primary solution approximation uh and the flux approximation σh are obtained under the standard quasi-uniform assumption on finite element partition. The superconvergence indicates an accuracy of O(h^r 2) for the least-squares mixed finite element approximation if Raviart-Thomas or Brezzi-DouglasFortin-Marini elements of order r are employed with optimal error estimate of O(h^r l).     

NODAL O（h^4）-SUPERCONVERGENCE IN 3D BY AVERAGING PIECEWISE LINEAR,BILINEAR, AND TRILINEAR FE APPROXIMATIONS

We construct and analyse a nodal O（h^4）-superconvergent FE scheme for approximating the Poisson equation with homogeneous boundary conditions in three-dimensional domains by means of piecewise trilinear functions. The scheme is based on averaging the equations that arise from FE approximations on uniform cubic, tetrahedral, and prismatic partitions. This approach presents a three-dimensional generalization of a two-dimensional averaging of linear and bilinear elements which also exhibits nodal O（h^4）-superconvergence （ultracon- vergence）. The obtained superconvergence result is illustrated by two numerical examples.     

A two-grid method with expanded mixed element for nonlinear reaction-diffusion equations

Expanded mixed finite element approximation of nonlinear reaction-diffusion equations is discussed. The equations considered here are used to model the hydrologic and bio-geochemical phenomena. To linearize the mixed-method equations, we use a two-grid method involving a small nonlinear system on a coarse gird of size H and a linear system on a fine grid of size h. Error estimates are derived which demonstrate that the error is O(△t + h k+1 + H 2k+2 d/2 ) (k ≥ 1), where k is the degree of the approximating space for the primary variable and d is the spatial dimension. The above estimates are useful for determining an appropriate H for the coarse grid problems.     

SURFACE FINITE ELEMENTS FOR

In this article we define a surface finite element method （SFEM） for the numerical solution of parabolic partial differential equations on hypersurfaces F in R^n＋1. The key idea is based on the approximation of F by a polyhedral surface Гh consisting of a union of simplices （triangles for n = 2, intervals for n = 1） with vertices on Г. A finite element space of functions is then defined by taking the continuous functions on Гh which are linear affine on each simplex of the polygonal surface. We use surface gradients to define weak forms of elliptic operators and naturally generate weak formulations of elliptic and parabolic equations on Г. Our finite element method is applied to weak forms of the equations. The computation of the mass and element stiffness matrices are simple and straightforward. We give an example of error bounds in the case of semi-discretization in space for a fourth order linear problem. Numerical experiments are described for several linear and nonlinear partial differential equations. In particular the power of the method is demorrstrated by employing it to solve highly nonlinear second and fourth order problems such as surface Allen-Cahn and Cahn-Hilliard equations and surface level set equations for geodesic mean curvature flow.     

On a class of self-similar processes with stationary increments in higher order Wiener chaoses

We study a class of self-similar processes with stationary increments belonging to higher order Wiener chaoses which are similar to Hermite processes. We obtain an almost sure wavelet-like expansion of these processes. This allows us to compute the pointwise and local Hölder regularity of sample paths and to analyse their behaviour at infinity. We also provide some results on the Hausdorff dimension of the range and graphs of multidimensional anisotropic self-similar processes with stationary increments defined by multiple Wiener–Itô integrals.     

A REPRESENTATION FORMULA RELATED TO SCHR(O)DINGER OPERATORS

Schr(o)dinger operator is a central subject in the mathematical study of quantum mechanics.Consider the Schrodinger operator H = -△ V on R, where △ = d2/dx2 and the potential function V is real valued. In Fourier analysis, it is well-known that a square integrable function admits an expansion with exponentials as eigenfunctions of -△. A natural conjecture is that an L2 function admits a similar expansion in terms of "eigenfunctions" of H, a perturbation of the Laplacian (see [7], Ch. Ⅺ and the notes), under certain condition on V.     

On a class of Riemann surfaces

It is considered the class of Riemann surfaces with dimT1 = 0, where T1 is a subclass of exact harmonic forms which is one of the factors in the orthogonal decomposition of the spaceΩH of harmonic forms of the surface, namely The surfaces in the class OHD and the class of planar surfaces satisfy dimT1 = 0. A.Pfluger posed the question whether there might exist other surfaces outside those two classes. Here it is shown that in the case of finite genus g, we should look for a surface S with dimT1 = 0 among the surfaces of the form Sg\K , where Sg is a closed surface of genus g and K a compact set of positive harmonic measure with perfect components and very irregular boundary.     

Qualitative Analysis of Bifurcating Solutions in the Lengyel-Epstein Model

One of the most fundamental problems in theoretical biology is to explain the mechanisms by which patterns and forms are created in the＇living world. In his seminal paper ＂The Chemical Basis of Morphogenesis＂, Turing showed that a system of coupled reaction-diffusion equations can be used to describe patterns and forms in biological systems. However, the first experimental evidence to the Turing patterns was observed by De Kepper and her associates（1990） on the CIMA reaction in an open unstirred reactor, almost 40 years after Turing＇s prediction. Lengyel and Epstein characterized this famous experiment using a system of reaction-diffusion equations. The Lengyel-Epstein model is in the form as follows     

(0,1 ;0)-INTERPOLATION ON INFINITE INTERVAL (-∞, +∞)

In this paper, we study the explicit representation and convergence of (0, 1; 0)-interpolation on infinite interval, which means to determine a polynomial of degree ≤ 3n - 2 when the function values are prescribed at two set of points namely the zeros of Hn(x) and H′n(x) and the first derivatives at the zeros of H′n(x).     

CONTENTS OF VOLUME 30

正~~     

Instructions for contributions

正Applied Mathematics-A Journal of Chinese Universities,Series B(Appl.Math.J.Chinese Univ.,Ser.B)is a comprehensive applied mathematics journal jointly sponsored by Zhejiang University,China Society for Industrial and Applied Mathematics,and Springer-Verlag.It is a quarterly journal with     

Journal of Mathematical Research with Applications Guide for Authors

正Journal overview:Journal of Mathematical Research with Applications(JMRA),formerly Journal of Mathematical Research and Exposition(JMRE)created in 1981,one of the transactions of China Society for Industrial and Applied Mathematics,is a home for original research papers of the highest quality in all areas of mathematics with applications.The target audience comprises:pure and applied mathematicians,graduate students in broad fields of sciences and technology,scientists and engineers interested in mathematics.     

Staircase transportation problems with superadditive rewards and cumulative capacities

A cumulative-capacitated transportation problem is studied. The supply nodes and demand nodes are each chains. Shipments from a supply node to a demand node are possible only if the pair lies in a sublattice, or equivalently, in a staircase disjoint union of rectangles, of the product of the two chains. There are (lattice) superadditive upper bounds on the cumulative flows in all leading subrectangles of each rectangle. It is shown that there is a greatest cumulative flow formed by the natural generalization of the South-West Corner Rule that respects cumulative-flow capacities; it has maximum reward when the rewards are (lattice) superadditive; it is integer if the supplies, demands and capacities are integer; and it can be calculated myopically in linear time. The result is specialized to earlier work of Hoeffding (1940), Fréchet (1951), Lorentz (1953), Hoffman (1963) and Barnes and Hoffman (1985). Applications are given to extreme constrained bivariate distributions, optimal distribution with limited one-way product substitution and, generalizing results of Derman and Klein (1958), optimal sales with age-dependent rewards and capacities.To our friend, Philip Wolfe, with admiration and affection, on the occasion of his 65th birthday.Research was supported respectively by the IBM T.J. Watson and IBM Almaden Research Centers and is a minor revision of the IBM Research Report [6].