Numerical integration of ordinary differential equations on manifolds

Summary This paper is concerned with the problem of developing numerical integration algorithms for differential equations that, when viewed as equations in some Euclidean space, naturally evolve on some embedded submanifold. It is desired to construct algorithms whose iterates also evolve on the same manifold. These algorithms can therefore be viewed as integrating ordinary differential equations on manifolds. The basic method “decouples” the computation of flows on the submanifold from the numerical integration process. It is shown that two classes of single-step and multistep algorithms can be posed and analyzed theoretically, using the concept of “freezing” the coefficients of differential operators obtained from the defining vector field. Explicit third-order algorithms are derived, with additional equations augmenting those of their classical counterparts, obtained from “obstructions” defined by nonvanishing Lie brackets.     

吴方法和带参数杨—Baxter方程的求解技术

回顾用吴文俊先生提出的消元方法求解带参数的杨—Ｂａｘｔｅｒ方程的技巧和主要结果。     

四角晶相HfO2(001)表面原子和电子结构研究

采用基于第一性原理的密度泛函理论研究了四角晶相二氧化铪(t-HfO₂)体相及 其(001)表面的原子几何与电子结构.理论计算结果表明，t-HfO₂(001)表面不会 产生重构现象.与体相电子结构相比, t-HfO₂(001)表面态密度明显高于体相态 密度.其次，表面原子的态密度更靠近费米能级(E_F)，价带往低能量处移动，并 有表面态产生.计算结果表明了t-HfO₂表面禁带宽度明显低于体相的禁带宽度. t-HfO₂(001)的表面态产生以及表面禁带宽度减小是由于Hf原子与O原子的配位 数减少，表面原子周围的环境发生变化而引起的. 关键词： 密度泛函理论 2(001)')" href="#">t-HfO₂(001) 表面电子结构     

铁分子Fe2的自旋极化效应

采用密度泛函方法(B3P86)对 Fe_2分子结构进行了优化.计算结果中未观察到自旋污染,基态波函数与高态波函数并未混杂,结果表明,Fe_2中有8个未配对电子,这些电子空间分布不同和自旋平行产生的自旋极化效应,使 Fe_2能量最低.计算结果表明,Fe_2分子的基态是~9∑_g~ ,并非~7Δ_u,进而表明 Fe_2的自旋平行效应比电子自旋配对效应要强.计算得到该分子基态的二阶、三阶和四阶力常数分别为1.4115×10~(-2)aJ/nm~2、-37.1751×10~3aJ/nm~3和 98.7596×10~4aJ/nm~4；光谱数据ω_eχ_e、B_e、α_e分别为0.3522、0.0345、 0.4963×10~(-4)cm~(-1)；离解能为3.5522eV,平衡键长为0.2137nm,振动频率为292.914cm~(-1)；并得到了 Murrel-Sorbie 函数.     

人体中的伯努利方程

本文从设计实验入手,形象诠释伯努利原理,并从物理学角度出发,根据人体体液流动的实际情况,运用伯努利原理及经典的伯努利方程,从人体血液循环、房水循环的压力形成和改变方面,把物理学的基本理论运用于分析人体体液压力的变化。     

SINGULAR PERTURBATION OF BOUNDARY VALUE PROBLEM FOR NONLINEAR HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS INVOLVING TWO SMALL PARAMETERS

The singularly perturbed boundary value problem for nonlinear higher order ordinary differential equation involving two small parameters has been considered. Under appropriate assumptions, for the three cases:ε/μ2→0(μ→0),μ2/ε→0 (ε→0) andε=μ2, the uniformly valid asymptotic solution is obtained by using the expansion method of two small parameters and the theory of differential inequality.     

Attractors for a Three-Dimensional Thermo-Mechanical Model of Shape Memory Alloys

Abstract In this note, we consider a Frémond model of shape memory alloys. Let us imagine a piece of a shape memory alloy which is fixed on one part of its boundary, and assume that forcing terms, e.g., heat sources and external stress on the remaining part of its boundary, converge to some time-independent functions, in appropriate senses, as time goes to infinity. Under the above assumption, we shall discuss the asymptotic stability for the dynamical system from the viewpoint of the global attractor. More precisely, we generalize the paper [12] dealing with the one-dimensional case. First, we show the existence of the global attractor for the limiting autonomous dynamical system; then we characterize the asymptotic stability for the non-autonomous case by the limiting global attractor. * Project supported by the MIUR-COFIN 2004 research program on “Mathematical Modelling and Analysis of Free Boundary Problems”.     

355nm Nd∶YAG激光在H_2中的高效一级斯托克斯转换

对脉冲Nd∶YAG激光(355 nm)在H2和H2∶He-Ar混合气体中的受激拉曼散射(SRS)进行了研究。在0.5 MPa的氢气中,同时测量到从二级反斯托克斯到三级斯托克斯的多波长输出,其总转化效率达88%;而高压下只剩下一级和二级斯托克斯输出,其中二级斯托克斯最大能量转化效率达44%(对应量子效率为63%)。由于高级斯托克斯的竞争,纯氢气中一级斯托克斯的最大能量转换效率不超过43%。通过向3 MPa氢气中掺入2 MPaAr气后,很好地抑制了二级斯托克斯的产生,从而获得了能量转换效率高达71%(对应量子效率为83%)的一级斯托克斯输出。对四波混频和级联受激拉曼散射在氢气多级斯托克斯产生中的作用以及惰性气体对它们的影响进行了讨论。     

一类非线性Schr(o)dinger方程的守恒差分法与Fourier谱方法

考察了一类带导数项的非线性Schrodinger方程的周期边值问题，提出了一种守恒的差分格式，在空间方向上采用Fourier谱方法，证明了格式的稳定性和收敛性．数值试验得到了与理论分析一致的结果．     

A rational approximation and its applications to nonlinear partial differential equations on the whole line

An orthogonal system of rational functions is discussed. Some inverse inequalities, imbedding inequalities and approximation results are obtained. Two model problems are considered. The stabilities and convergences of proposed rational spectral schemes and rational pseudospectral schemes are proved. The techniques used in this paper are also applicable to other problems on the whole line. Numerical results show the efficiency of this approach.