ADAPTIVE FINITE ELEMENT APPROXIMATION FOR A CLASS OF PARAMETER ESTIMATION PROBLEMS

In this paper, we study adaptive finite element discretisation schemes for a class of parameter estimation problem. We propose to efficient algorithms for the estimation problem use adaptive multi-meshes in developing We derive equivalent a posteriori error estimators for both the state and the control approximation, which particularly suit an adaptive multi-mesh finite element scheme. The error estimators are then implemented and tested with promising numerical results.     

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.     

Uniform analysis of a stabilized hybrid finite element method for Reissner-Mindlin plates

This paper presents a low order stabilized hybrid quadrilateral finite element method for ReissnerMindlin plates based on Hellinger-Reissner variational principle,which includes variables of displacements,shear stresses and bending moments.The approach uses continuous piecewise isoparametric bilinear interpolations for the approximations of the transverse displacement and rotation.The stabilization achieved by adding a stabilization term of least-squares to the original hybrid scheme,allows independent approximations of the stresses and moments.The stress approximation adopts a piecewise independent 4-parameter mode satisfying an accuracy-enhanced condition.The approximation of moments employs a piecewise-independent 5-parameter mode.This method can be viewed as a stabilized version of the hybrid finite element scheme proposed in [Carstensen C,Xie X,Yu G,et al.A priori and a posteriori analysis for a locking-free low order quadrilateral hybrid finite element for Reissner-Mindlin plates.Comput Methods Appl Mech Engrg,2011,200:1161-1175],where the approximations of stresses and moments are required to satisfy an equilibrium criterion.A priori error analysis shows that the method is uniform with respect to the plate thickness t.Numerical experiments confirm the theoretical results.     

MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

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.     

CONSTRAINT-PRESERVING ENERGY-STABLE SCHEME FOR THE 2D SIMPLIFIED ERICKSEN-LESLIE SYSTEM

Here we consider the numerical approximations of the 2D simplified Ericksen-Leslie system.We first rewrite the system and get a new system.For the new system,we propose an easy-to-implement time discretization scheme which preserves the sphere constraint at each node,enjoys a discrete energy law,and leads to linear and decoupled elliptic equations to be solved at each time step.A discrete maximum principle of the schemc in the finite element form is also proved.Some numerical simulations are performed to validate the scheme and simulate the dynamic motion of liquid crystals.     

FINITE ELEMENT METHODS FOR THE NAVIER-STOKES EQUATIONS BY H（div） ELEMENTS

We derived and analyzed a new numerical scheme for the Navier-Stokes equations by using H（div） conforming finite elements. A great deal of effort was given to an establishment of some Sobolev-type inequalities for piecewise smooth functions. In particular, the newly derived Sobolev inequalities were employed to provide a mathematical theory for the H（div） finite element scheme. For example, it was proved that the new finite element scheme has solutions which admit a certain boundedness in terms of the input data. A solution uniqueness was also possible when the input data satisfies a certain smallness condition. Optimal-order error estimates for the corresponding finite element solutions were established in various Sobolev norms. The finite element solutions from the new scheme feature a full satisfaction of the continuity equation which is highly demanded in scientific computing.     

ON LOCKING-FREE FINITE ELEMENT SCHEMES FOR THREE-DIMENSIONAL ELASTICITY

In the present paper,the authors discuss the locking phenomenon of the finite elementmethod for three-dimensional elasticity as the Laméconstant λ→∞.Three kinds of finiteelements are proposed and analyzed to approximate the three-dimensional elasticity withpure displacement boundary condition.Optimal order error estimates which are uniformwith respect to λ∈(0,+∞)are obtained for three schemes.Furthermore,numerical resultsare presented to show that,our schemes are locking-free and and the trilinear conformingfinite element scheme is locking.     

THE REDUCED BASIS TECHNIQUE AS A COARSE SOLVER FOR PARAREAL IN TIME SIMULATIONS

In this paper, we extend the reduced basis methods for parameter dependent problems to the parareal in time algorithm introduced by Lions et al. [12] and solve a nonlinear evolutionary parabolic partial differential equation. The fine solver is based on the finite element method or spectral element method in space and a semi-implicit Runge-Kutta scheme in time. The coarse solver is based on a semi-implicit scheme in time and the reduced basis approximation in space. Of[line-online procedures are developed, and it is proved that the computational complexity of the on-line stage depends only on the dimension of the reduced basis space （typically small）. Parareal in time algorithms based on a multi-grids finite element method and a multi-degrees finite element method are also presented. Some numerical results are reported.     

FINITE ELEMENT METHOD FOR INCOMPRESSIBLE MISCIBLE DISPLACEMENT IN POROUS MEDIA

Under the assumptions of nonlinear finite element and △t=o(h),Ewing and Wheeler discussed a Galerkin method for the single phase incompressible miscible displacement of one fluid by another in porous media.In the present paper we give a finite element scheme which weakens the △t=o(h)-restriction to △t=o(h~ε),0<ε≤1/2.Furthermore,this scheme is suitable for both linear element and nonlinear element.We also derive the optimal approximation estimates for concentration c,its gradient ▽c and the gradient ▽p of the pressure p.     

一类矩阵的AOR迭代收敛性分析及其与SOR迭代的比较

1 引言 许多实际问题最后常归结为解一个或一些矩阵的线性代数方程组Ax=b （1.1）这里讨论A为（1，1）相容次序矩阵的情形。     

双模范畴间同构态射的构造

本文研究了双模范畴上的同构态射.利用代数上常用的构造方法,给出了模范畴R(C:T)和L(C:T)的定义及双模范畴TMC和TMCM之间的同构映射和证明过程,将代数上双模范畴的一些重要结论进行了推广.     

关于“停走”生成器概率模型中的符合率问题

本文首先建立了“停走”生成器辅出序列的概率模型,给出了“停走”生成器输出序列与其线性移位寄存器序列之间的符合率的计算公式。     

半序集的碰撞数与分层深度贪婪算法

设Ｐ＝（Ｘ，≤）是一个半序集，本文在关于碰撞数的深度贪婪算法的基础上，直接证明了对任意的Ｐ存在一个最优的ＤＬＧ扩张，给出了ＤＬＧ半序集的定义，并证明了半序集Ｐ是ＤＬＧ半序集的一个充分条件，最后给出了ＤＬＧ扩张算法。     

非齐次空间上的奇异积分交换子的有界性

本文在非齐次空间上给出了交换子[b,T](f)=bTf(x)-T(bf)(x)在b(x)是Lipschitz函数时的 Lp(p>1)有界性.     

THE SOLUTION TO THE SK-MODEL

IIntroductlonAs one ofwell－kn。mean ield models for spin glasses,the SK(Sherrin红on－Kirkpatri山)model has been studied by many authors恤叫2]nd[81,andthe references therein)．Particu－larl儿丁劝a以andls]repm眈* some quite lmerestingresults on It in his one－hour Invited talk tthe International Congress ofMathem航icians held t Berlin in August,ig98．In mathematical terms;the SK－Model Is the study of a cert。n random measure on Z。:＝{一1;1}”for a natural。mber N．Z。Is called configu…     

三项指数和四次均值的精确计算公式(2)

本文进一步深入研究了三项指数和四次均值的计算问题.运用指数和的相关性质并结合求解同余方程组的方法与技巧,利用两种不同的方法获得了两个精确的均值计算公式,揭示了三项指数和的计算与同余方程组解的个数之间的本质联系,推广了已有的结果.     

非齐次空间上的奇异积分交换子的有界性

本文在非齐次空间上给出了交换子[b,T](f)=bTf(x)-T(bf)(x)在b(x)是Lipschitz函数时的Lp(p＞1)有界性.     

用离散纵标法计算平板几何反应堆厚度临界尺度的合理性~

自从1944年 chandrasekhar 在辐射迁移现象计算中使用离散纵标法之后,该方法在核反应堆实际计算中有了广泛的应用,因而引起了许多数学工作者的关心.他们去研究和证明该方法的合理性,并已得到很多结果(如[2—10]).本文的目的是证明用离散纵标法计算平板几何反应堆关于厚度的临界尺度本征值的合理性.这里我们讨论介质体     

关于TLS和LS解的扰动分析

1．引言本文采用卜］的记号．最小二乘（LS）和总体最小二乘（TLS）是科学计算中的两种重要方法．尤是TLS，近来已有多篇论文讨论［1－6，8－16］．奇异值分解（SVD）和CS分解是研究TLS和LS的重要工具．令ACm，BCm,C=(A,B)，A和C的SVD分别为（1．1）（1．2）其中P51为某个正整数，U，U，V，V均为西矩阵，UI，UI，VI，VI为上述矩阵的前P列，z1一山。g（。1，…，内），】2＝di。g（内十l，…，。小】1＝dl。g（61；…，站，】2二diag（4＋1；…，dk），。l三··2。120和dl三…三d。20分别为C和A的奇异值，Z＝mhfm．n十以…