THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

Hilbert空间中涉及一类新的具误差迭代程序的平衡问题的黏性逼近法

In this paper, we introduce an iterative scheme with error by the viscosity approximation method for finding a common element of the set of solutions of an equilibrium problem and the set of fixed points of a nonexpansive mapping in a Hilbert space. A strong convergence theorem is given, which generalizes all the results obtained by S.Takahashi and W.Takahashi in 2007. In addition, some of the methods applied in this paper improve those of S.Takahashi and W.Takahashi.     

ROBUST HIGH ORDER CONVERGENCE OF AN OVERLAPPING SCHWARZ METHOD FOR SINGULARLY PERTURBED SEMILINEAR REACTION-DIFFUSION PROBLEMS

In this article we propose an overlapping Schwarz domain decomposition method for solving a singularly perturbed semilinear reaction-diffusion problem. The solution to this problem exhibits boundary layers of width O（√ε ln（1/√ε）） at both ends of the domain due to the presence of singular perturbation parameter ε. The method splits the domain into three overlapping subdomains, and uses the Numerov or Hermite scheme with a uniform mesh on two boundary layer subdomains and a hybrid scheme with a uniform mesh on the interior subdomain. The numerical approximations obtained from this method are proved to be almost fourth order uniformly convergent （in the maximum norm） with respect to the singular perturbation parameter. Furthermore, it is proved that, for small ε, one iteration is sufficient to achieve almost fourth order uniform convergence. Numerical experiments are given to illustrate the theoretical order of convergence established for the method.     

TWO-LEVEL MULTISCALE FINITE ELEMENT METHODS FOR THE STEADY NAVIER-STOKES PROBLEM

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.     

Design of a new dynamical core for global atmospheric models based on some efficient numerical methods

A careful study on the integral properties of the primitive hydrostatic balance equations for baroclinic atmosphere is carried out, and a new scheme to design the global adiabatic model of atmospheric dynamics is presented. This scheme includes a method of weighted equal-area mesh and a fully discrete finite difference method with quadratic and linear conservations for solving the primitive equation system. Using this scheme, we established a new dynamical core with adjustable high resolution acceptable to the available     

CONSERVATION OF THREE—POINT COMPACT SCHEMES ON SINGLE AND MULTIBLOCK PATCHED GRIDS FOR HYPERBOLIC PROBLEMS

For nonlinear hyperbloic problems,Conservation of the numerical scheme is important for convergence to the correct weak solutions.In this paper the the conservation of the well-known compact scheme up to fourth order of accuracy on a single and uniform grid is studied,and a conservative interface treatment is derived for compact schemes on patched grids .For a pure initial value problem,the compact scheme is shown to be equivalent to a scheme in the usual conservative form .For the case of a mixed initial boundary value problem,the compact scheme is conservative only if the rounding errors are small enough.For a pactched grid interface,a conservative interface condition useful for mesh fefiement and for parallel computation is derived and its order of local accuracy is analyzed.     

集合$\tau$-标准部分的闭性

In this paper,the closeness of the τ-standard part of a set is discussed.Some related propositions of the τ-neighborhood system of a set are given.And then some related conclusions of the τ-monad of a set and the τ-standard part of a set are presented.And based on it,the necessary and sufficient conditions of the enlarged model and the saturated model are showed.Finally,some sufficient conditions that the τ-standard part of a set is closed are proved in the enlarged model and the saturated model.     

The Closeness of the τ-Standard Part of a Set

In this paper,the closeness of the τ-standard part of a set is discussed.Some related propositions of the τ-neighborhood system of a set are given.And then some related conclusions of the τ-monad of a set and the τ-standard part of a set are presented.And based on it,the necessary and sufficient conditions of the enlarged model and the saturated model are showed.Finally,some sufficient conditions that the τ-standard part of a set is closed are proved in the enlarged model and the saturated model.     

一类多元加权最小乖正交多项式

This paper introduces a new notion of weighted least-square orthogonal polynomials in multivariables from the triangular form. Their existence and uniqueness is studied and some methods for their recursive computation are given. As an application, this paper constructs a new family of Pade-type approximates in multi-variables from the triangular form.     

ON THE ERROR ESTIMATES OF A NEW OPERATE SPLITTING METHOD FOR THE NAVIER-STOKES EQUATIONS

In this paper, a new operator splitting scheme is introduced for the numerical solution of the incompressible Navier-Stokes equations. Under some mild regularity assumptions on the PDE solution, the stability of the scheme is presented, and error estimates for the velocity and the pressure of the proposed operator splitting scheme are given.     

Vector subdivision schemes in (L_p(R~5))~r(1≤p≤∞) spaces

The purpose of this paper is to investigate the refinement equations of the formwhere the vector of functions = (1, … ,r)T is in (LP(R8))T,1 ≤ p ≤∞, α(α),α ∈ Z5, is a finitely supported sequence of r × r matrices called the refinement mask, and M is an s x a integer matrix such that limn→ ∞ M-n = 0. In order to solve the refinement equation mentioned above, we start with a vector of compactly supported functions (0 ∈ (LP(R8))r and use the iteration schemes fn := Qan0,n = 1,2,…, where Qa is the linear operator defined on (Lp(R8))r given byThis iteration scheme is called a subdivision scheme or cascade algorithm. In this paper, we characterize the Lp-convergence of subdivision schemes in terms of the p-norm joint spectral radius of a finite collection of some linear operators determined by the sequence a and the set B restricted to a certain invariant subspace, where the set B is a complete set of representatives of the distinct cosets of the quotient group Z8/MZ8 containing 0.     

参数曲面用插值三角平面片逼近的误差估计

Triangulation, subdivision and intersection of the surface are the basic and common operations in Computer Aided Design and computer graphics. A complicated surface can be approximated by some triangular patches. The key technique of this problem is the estimate of the approximate error, this work improves the previous result and obtains a new result. These results are valuable for improving the triangulation and subdivision algorithm of parametric surface in the design system.     

Design of Finite Element Tools for Coupled Surface and Volume Meshes

Many problems with underlying variational structure involve a coupling of volume with surface effects.A straight-forward approach in a finite element discretiza- tion is to make use of the surface triangulation that is naturally induced by the volume triangulation.In an adaptive method one wants to facilitate"matching"local mesh modifications,i.e.,local refinement and/or coarsening,of volume and surface mesh with standard tools such that the surface grid is always induced by the volume grid. We describe the concepts behind this approach for bisectional refinement and describe new tools incorporated in the finite element toolbox ALBERTA.We also present several important applications of the mesh coupling.     

Calculating the vertex unknowns of nine point scheme on quadrilateral meshes for diffusion equation

In the construction of nine point scheme,both vertex unknowns and cell-centered unknowns are introduced,and the vertex unknowns are usually eliminated by using the interpolation of neighboring cell-centered unknowns,which often leads to lose accuracy.Instead of using interpolation,here we propose a different method of calculating the vertex unknowns of nine point scheme,which are solved independently on a new generated mesh.This new mesh is a Vorono¨i mesh based on the vertexes of primary mesh and some additional points on the interface.The advantage of this method is that it is particularly suitable for solving diffusion problems with discontinuous coeffcients on highly distorted meshes,and it leads to a symmetric positive definite matrix.We prove that the method has first-order convergence on distorted meshes.Numerical experiments show that the method obtains nearly second-order accuracy on distorted meshes.     

NON C~0 NONCONFORMING ELEMENTS FOR ELLIPTIC FOURTH ORDER SINGULAR PERTURBATION PROBLEM

In this paper we give a convergence theorem for non C^0 nonconforming finite element to solve the elliptic fourth order singular perturbation problem. Two such kind of elements, a nine parameter triangular element and a twelve parameter rectangular element both with double set parameters, are presented. The convergence and numerical results of the two elements are given.     

实对称半正定矩阵LDL~T分解的存在性与唯一性及有关问题

The present paper proves the existence and uniqueness (in some given sense) ofthe LDL~T decomposition of real symmetric non-negative definite matrices, where Lis a unit lower triangular matrix with real elements and D is a diagonal matrixwith real elements. The proof is made in a constructive way. By taking advantageof this decomposition, a criterion for the consistency of the linear equation with sucha coefficient matrix and its whole solution set (or the least-squares solution if it isinconsistent) are obtained. Since it involves no row or column permutation, the pro-cess may be combined with any sparse technique on the computer, and hence is ofpractical importance in treating the large scale sparse matrices derived from suchproblems as the structure design by finite elements methods. Finally, the stability ofsuch a decomposition is discussed and a backward error analysis is given.     

The optimal solution of multi-kernel regularization learning

In regularized kernel methods, the solution of a learning problem is found by minimizing a functional consisting of a empirical risk and a regularization term. In this paper, we study the existence of optimal solution of multi-kernel regularization learning. First, we ameliorate a previous conclusion about this problem given by Micchelli and Pontil, and prove that the optimal solution exists whenever the kernel set is a compact set. Second, we consider this problem for Gaussian kernels with variance σ∈(0,∞), and give some conditions under which the optimal solution exists.     

Convergence of subdivision schemes in L p spaces

Abstract. In this paper it is proved that Lp solutions of a refinement equation exist if and only ifthe corresponding subdivision scheme with suitable initial function converges in Lp without anyassumption on the stability of the solutions of the refinement equation. A characterization forconvergence of subdivision scheme is also given in terms of the refinement mask. Thus a com-plete answer to the relation between the existence of Lp solutions of the refinement equation andthe convergence of the corresponding subdivision schemes is given.     

带分数阶Robin边界条件的时间-空间分数阶扩散方程的有限差分方法（英文）

In this paper, an efficient numerical method is proposed to solve the Caputo-Riesz fractional diffusion equation with fractional Robin boundary conditions. We approximate the Riesz space fractional derivatives using the fractional central difference scheme with second-order accurate. A priori estimation of the solution of the numerical scheme is given, and the stability and convergence of the numerical scheme are analyzed.Finally, a numerical example is used to verify the accuracy and efficiency...     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.