晶体微观结构的数值计算

晶体微观结构是晶体材料在特定物理条件下其多个能量极小平衔态在空间形成的某种微尺度的规则分布．几何非线性的连续介质力学理论可以用能量极小化原理来解释晶体微观结构的形成，并用Young测度来刻画平衡态各变体在空间的概率分布．定性的理解与定量地分析和计算晶体材料的微观结构对于发展和改进高级晶体功能材料，如形状记忆合金、铁电体、磁至伸缩材料等，有重要的意义．本文回顾了近年来晶体微观结构数值计算方面的最新进展．介绍了计算晶体微观结构的几种数值方法及有关的数值分析结果。     

带对流的稳态Stefan问题的有限元逼近

关于带有对流情形的稳态Stefan问题,其中假设液相部分的流动由Stokes方程确定,Canon,DiBennedetto,Knightly曾作过理论研究.本文将讨论其有限元逼近问题,并且得到了在合理正则性假设下的误差估计.     

Stokes问题有限元逼近中求导数的提取法

1.引言 我们知道Poisson方程和平面弹性问题的解的导数的近似值可以通过所谓提取公式得到,而不必对近似解直接求导数.这样我们可以得到具有与近似解本身同阶精度的导数的近似值.这一方法已被用于基于插值误差的后验误差估计及相应的自适应有限元方法中本文将这一方法应用于Stokes问题的有限元逼近,从Stokes方程的解的     

浅水波方程的一种基于特征方向的Galerkin方法

本文给出求解二维浅水波方程组的一种Crank-Nicolson型基于特征方向的Galerkin有限元方法，证明了该方法的一个误差估计结果，并给出了该方法的一个算例．     

一个改进的Reissner-Mindlin矩形元

本文提出了一个改进的Reissner-Mindlin矩形非协调元方法：旋度用连续双线性元逼近，横向位移用旋转矩形非协调元逼近，而作为中间变量的剪切力用增广的分片常数元逼近，我们证明：该方法具有关于板厚一致稳定性和一致最优收敛性。     

定常的热传导-对流问题的Galerkin/Petrov最小二乘混合元方法

1.引言 热传导-对流问题是大气动力学中的一个重要的方程,这个方程组也称为强迫耗散的非线性系统方程组,其较Navier-Stokes方程多了一个未知函数温度场,且温度与速度和压力之间存在着复杂的非线性关系.从热动力学可知,任何运动都会产生热量即有温度,而且温度与速度和压力之间必定互相转化,因此对该非线性系统的研究更具有实际意义.[1]先对     

一类基于变量分离的稳定化混合有限元方法

Based on a seperated model for saddle-point problems, we develop a new sta-bilized mixed finite element method. Such a model consists of two subproblems with respect to the primal and the dual variables, respectively. We show that the new method is coercive and that optimal error bounds hold. As an application,the nearly incompressible elastic problem is analyzed with our method.     

ROBUSTNESS OF AN UPWIND FINITE DIFFERENCE SCHEME FOR SEMILINEAR CONVECTION-DIFFUSION PROBLEMS WITH BOUNDARY TURNING POINTS

We consider a singularly perturbed semilinear convection-diffusion problem with a boundary layer of attractive turning-point type. It is shown that its solution can be decomposed into a regular solution component and a layer component. This decomposi-tion is used to analyse the convergence of an upwinded finite difference scheme on Shishkin meshes.     

THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

In this paper, we investigate the coupling of natural boundary element and finite element methods of exterior initial boundary value problems for hyperbolic equations. The governing equation is first discretized in time, leading to a time-step scheme, where an exterior elliptic problem has to be solved in each time step. Second, a circular artificial boundary FR consisting of a circle of radius R is introduced, the original problem in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formula are obtained in the infinite domainΩ2 outside circle of radius R. The coupled variational formulation is given. Only the function itself, not its normal derivative at artificial boundary ΓR, appears in the variational equation, so that the unknown numbers are reducedand the boundary element stiffness matrix has a few different elements. Such a coupled method is superior to the one based on direct boundary element method. This paper discusses finite element discretization for variational problem and its corresponding numerical technique, and the convergence for the numerical solutions. Finally, the numerical example is presented to illustrate feasibility and efficiency of this method.     

Asymptotic Expansions of the Heat Kernel of the Laplacian for General Annular Bounded Domains with Robin Boundary Conditions：Further Results

The asymptotic expansions of the trace of the heat kernel θ(t)=∑^∞v=1^exp(-tλv) for small positive t,where {λv} are the eigenvalues of the negative Laplacian -△n=-∑^ni=1(D/Dx^1)^2 in R^2(n=2 or 3),are studied for a general annular bounded domain Ω with a smooth inner boundary DΩ1 and a smooth outer boundary DΩ2,where a finite number of piecewise smooth Robin boundary conditions(D/Dnj γh)Ф=0 on the components Гj(j= 1,...,m) of (DΩ1 and on the components Гj (j=k 1,…,m) of of DΩ2 are considered such that DΩl=U^kj=lГj and DΩ2= U^m=k 1Гj and where the coefficients γj(j=1,...,m) are piecewise smooth positive functions. Some applications of θ(t) for an ideal gas enclosed in the general annular bounded domain Ω are given. Further results are also obtained.