Numerical Approximation and Error Analysis for the Timoshenko Beam Equations with Boundary Feedback

In this paper,the numerical approximation of a Timoshenko beam with bound- ary feedback is considered.We derived a linearized three-level difference scheme on uniform meshes by the method of reduction of order for a Timoshenko beam with boundary feedback.It is proved that the scheme is uniquely solvable,unconditionally stable and second order convergent in L_∞norm by using the discrete energy method. A numerical example is presented to verify the theoretical results.     

有限分形介质中分数阶反应扩散方程及其解析解

建立了有限分形介质中具有吸附效应的分数阶反应扩散积分方程.利用Lap lace变换、广义有限H ankel变换及其相应的逆变换得到了以M ittag-Leffler函数为主要形式的解析解,并研究了解的渐近性态.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome (OSAHS), three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed three-dimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosing diseases related to anatomical structure and function of the upper airway. The project supported by the National Natural Science Foundation of China (10672036, 10472025 and 10421002), the Natural Science Foundation of Liaoning Province (20032109). English text was polished by Yunming Chen.     

����Լ���������������MinimaxԤ��

本文研究了独立随机变量之和的绝对矩的几个性质, 其中包括$\ep|X+Y|-\ep|X-Y|$的表达式, 这里$X$和$Y$是相互独立的随机变量.     

四阶方程两点边值问题Hermite有限元解的渐近展式与外推

1引言有限元解的渐近展式是提高微分方程数值解精度的重要工具,比如亏量校正和外推就是建立在有限元解的渐近展式的基础之上．许多作者对此进行了大量的研究(见[1]-[4]),特别是文[1],提出了在研究有限元解的渐近展式中十分有用的能量嵌入技巧．本文利用能量嵌入定理得到了四阶方程两点边值问题Hermite有限元解及其二阶平均导数的渐近展式,进一步我们还讨论了它们的Richardson外推公式．考虑四阶方程两点边值问题     

Sparse generalized Fourier transforms

Block-diagonalization of sparse equivariant discretization matrices is studied. Such matrices typically arise when partial differential equations that evolve in symmetric geometries are discretized via the finite element method or via finite differences. By considering sparse equivariant matrices as equivariant graphs, we identify a condition for when block-diagonalization via a sparse variant of a generalized Fourier transform (GFT) becomes particularly simple and fast. Characterizations for finite element triangulations of a symmetric domain are given, and formulas for assembling the block-diagonalized matrix directly are presented. It is emphasized that the GFT preserves symmetric (Hermitian) properties of an equivariant matrix. By simulating the heat equation at the surface of a sphere discretized by an icosahedral grid, it is demonstrated that the block-diagonalization is beneficial. The gain is significant for a direct method, and modest for an iterative method. A comparison with a block-diagonalization approach based upon the continuous formulation is made. It is found that the sparse GFT method is an appropriate way to discretize the resulting continuous subsystems, since the spectrum and the symmetry are preserved. AMS subject classification (2000)  43A30, 65T99, 20B25     

有限长FGM圆筒的热/机应力的解析解

利用微分方程的级数求解方法,分析了两端简支的有限长功能梯度圆筒的轴对称稳态热弹性问题,推导出了稳态温度场与应力场的解析解。分析中采用指数函数模型来描述FGM圆筒中材料性能在厚度方向的连续变化,同时忽略温度对材料性能的影响。另外,论文以金属钼和多铝红柱石制成的功能梯度圆筒为例,给出了稳态温度场和应力场的数值结果。     

Error estimate for the approximation of nonlinear conservation laws on bounded domains by the finite volume method

In this paper we derive a priori and a posteriori error estimates for cell centered finite volume approximations of nonlinear conservation laws on polygonal bounded domains. Numerical experiments show the applicability of the a posteriori result for the derivation of local adaptive solution strategies.

     

3004铝合金“反常”锯齿屈服现象的研究

在应变速率为5.56×10^-5s^-1—5.56×10^-3s^-1的范围内，在不同温度下(从223K至773K)，对3004铝合金进行系列拉伸试验，探索其锯齿屈服规律；通过激活能的计算、内耗研究、微观组织观察和能谱分析，探讨锯齿屈服的机理与物理本质.结果表明，3004铝合金在形变过程中会出现动态应变时效现象；发现了一种“反常”的锯齿屈服现象：在出现锯齿屈服的温区内，存在锯齿屈服临界应变量转变温度T_t， 关键词： 动态应变时效 锯齿屈服 铝合金 内耗