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.     

准移动闭塞系统的元胞自动机模型及列车延迟传播规律的研究

在NaSch模型的基础上，针对铁路交通的特点提出一种用于模拟准移动闭塞系统的元胞自动机模型.应用该模型模拟了准移动闭塞系统列车延迟传播的交通现象，分析了准移动闭塞系统中的轨道定位单元长度、发车时间间隔、初始延迟时间等因素对列车延迟传播的影响. 关键词： 元胞自动机 交通流 准移动闭塞     

识别过渡金属离子的1,8-萘二甲酰亚胺多胺衍生物荧光传感器的研究

设计合成了用以检测过渡金属离子的荧光化学敏感器体系,它们是由1,8-萘二酰亚胺为荧光团,多胺衍生物为金属离子受体组成．在室温下对其光物理性质的研究中发现,在没有加入过渡金属离子时,由于体系内存在有效的光诱导电子转移过程使得荧光团的荧光被淬灭．加入过渡金属离子后,金属离子受体中的氮原子和过渡金属离子之间的配位作用阻断了光诱导电子转移过程,体系的荧光增强．不同的金属离子受体表现出了和过渡金属离子不同的配位识别能力,并且通过荧光的变化传递出受体-金属离子作用的信息．     

水斗非定常自由水膜流三维贴体数值模拟

本研究采用水斗三维非正交贴体坐标系进行了非定常自由水膜流动的数值解析。对不规则水斗内表面采用三维非正交贴体坐标系下离散点进行拟合,推导了曲面离散点的法向矢量和曲面微元面高斯曲率、平均曲率等几何特征量的计算公式,进而导出流体粒子在运动方向上曲率计算式。在水斗三维贴体坐标系中,还推导了流体粒子在水斗曲面上的运动控制方程。最后对某水轮机水斗内表面非定常自由水膜流进行了数值模拟,得到其非定常水膜流态分布。     

超强脉冲激光在低密度等离子体中的相对论自导引效应

分析了相对论条件下激光超短脉冲在等离子体中的传输特性 ,在傍轴近似和慢变振幅近似条件下 ,推导了折射率、电子密度、静电场以及电子空腔尺度的表达式。当激光功率超过产生自导引阈值功率时 ,激光束斑沿着传输光轴方向振荡。在有质动力产生的压力非常强时 ,聚焦光束中央部分的电子被全部排开形成电子空腔。给出了电子空腔的尺寸以及在出现电子空腔时的处理方法。在超过形成电子空腔的阈值功率 (Pc≈ 2 .5TW )时 ,空腔的尺度几乎与激光功率无关 ,这意味着电子空腔阻止了激光脉冲的进一步聚焦。     

Surface structures formed by individual adsorption and coadsorption of Mn and Bi on Cu(0 0 1), studied by LEED

We have studied the individual adsorption of Mn and Bi, and their coadsorption on Cu(0 0 1) by low-energy electron diffraction (LEED). For Mn, we have determined the c(2 × 2) structure formed at 300 K, whose structure had been determined by several methods. We reconfirmed by a tensor LEED analysis that it is a substitutional structure and that a previously reported large corrugation (0.30 Å) between substitutional Mn and remaining surface Cu atoms coincides perfectly with the present value. In the individual adsorption of Bi, we have found a c(4 × 2) structure, which is formed by cooling below ∼250 K a surface prepared by Bi deposition of ∼0.25 ML coverage at 300 K where streaky half-order LEED spots appear. The c(4 × 2) structure has been determined by the tensor LEED analysis at 130 K and it is a substitutional structure. In the coadsorption, we found a c(6 × 4) structure, which has been determined by the tensor LEED analysis. It is very similar to the previously determined structure of the c(6 × 4) formed by coadsorption of Mg and Bi, and embedded MnBi₄ clusters are arranged in the top Cu layer instead of MgBi₄. Large lateral displacements of Bi atoms in the c(6 × 4)-(Mn + Bi) suggest that the Mn atoms undergo the size-enhancement caused by their large magnetic moment.     

Roles of SiH4 and SiF4 in growth and structural changes of poly-Si films

The structural properties of polycrystalline silicon films, prepared by plasma enhanced chemical vapor deposition system, with different flow rates of SiH₄/SiF₄ mixtures at 300 °C were investigated. This study indicates that the low hydrogen coverage on the growing surface, under optimum fluorine radicals, will be leaded to an improvement of crystallized area as compared with case of high hydrogen coverage surface. Moreover, the studies of the role of SiH₄ and SiF₄ radicals show that the SiH₄ radicals are important in the nucleation and growth of grains. However, SiF₄ radicals are effective in the structural change of grain boundaries regions and by this way, in the present system, establish the growth of grains under the dominant 〈1 1 0〉 direction. The stress investigation indicates that addition of high flow rate of SiF₄ in amorphous film, results in the nearly stress free films. Finally, we found that the changes in g-value reflect the changes in the intrinsic compressive and tensile stress in the both polycrystalline and amorphous silicon films.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulation of the Action of a Permeable Barrier on an Ideal Incompressible Channel Flow

Plane ideal incompressible flow in a rectangular channel partitioned by a thin permeable barrier (lattice) is considered. In flowing through the lattice the stream suddenly (jumpwise) changes direction and loses energy. The flow is assumed to be vortical; the vorticity is discontinuous on the lattice. A mathematical formulation of the problem for the stream function is proposed in the form of a nonlinear elliptic equation with coefficients discontinuous on the lattice line. A numerical solution is constructed using the finite-element iteration method. The results of the numerical simulation show how the flow velocity profile in the channel can be controlled by means of permeable barriers.