Investigations with the finite element method on the Cyber 205

We are trying to investigate systematically the application of the finite element method (FEM) for solving the Schrödinger equation. The present paper is devoted to the calculation of vibrational transition probabilities for the collinear reactive system A + BC (i.e. H+H₂ and their isotopes). The calculations are fully two-dimensional and the results are compared with earlier FEM calculations and conventional basis set expansion methods using the the R-matrix or S-matrix propagation.We made extensive analysis of FEM on the vector-computer Cyber 205 and developed a vector code for the efficient use in two dimensions, so that in the near future applications even in three dimensions will be possible.For the hydrogen exchange reactions we investigated the following isotope combinations: (a) H + H₂, b) H + DH, D + HD and H + MuH (symmetric reaction), (c) D + HH, H + DD and Mu + DD (asymmetric reaction). We calculated the transition probabilities for up to five open vibrational channels and found excellent agreement with known exact values.     

On the applicability of the HSAB principle through the use of improved computational schemes for chemical hardness evaluation

Finite difference schemes, named Compact Finite Difference Schemes with Spectral-like Resolution, have been used for a less crude approximation of the analytical hardness definition as the second-order derivative of the energy with respect to the electron number. The improved computational schemes, at different levels of theory, have been used to calculate global hardness values of some probe bases, traditionally classified as hard and soft on the basis of their chemical behavior, and to investigate the quantitative applicability of the HSAB principle. Exchange acid-base reactions have been used to test the HSAB principle assuming the reaction energies as a measure of the stabilization of product adducts.     

STRUCTURE EVOLUTION OF POLYMER CHAINS FOR NECKING FORMATION IN HIGH-SPEED FIBER SPINNING PROCESS

Finite element method is used to simulate the high-speed melt spinning process,based on the equation system proposed by Doufas et al.Calculation predicts a neck-like deformation,as well as the related profiles of velocity,diameter, temperature,chain orientation,and crystallinity in the fiber spinning process.Considering combined effects on the process such as flow-induced crystallization,viscoelasticity,filament cooling,air drag,inertia,surface tension and gravity,the simulated material flow behaviors are consistent with those observed for semi-crystalline polymers under various spinning conditions.The structure change of polymer coils in the necking region described by the evolution of conformation tensor is also investigated.Based on the relaxation mechanism of macromolecules in flow field different types of morphology change of polymer chains before and in the neck are proposed,giving a complete prospect of structure evolution and crystallization of semi-crystalline polymer in the high speed fiber spinning process.     

Variational reactive scattering calculations: computational optimization strategies

Summary We have developed efficient and accurate techniques for the calculation of quantum mechanical reaction probabilities of atom-diatom exchange reactions in the gas phase, and we have optimized a computer code employing these techniques and applied it sucessfully to several systems. In this paper we consider further strategies for improving the algorithm to allow even more demanding applications. In this context, improvement means that equivalent results can be obtained using fewer computational resources (computer time or storage) or that an equivalent expenditure of resources can yield higher accuracy. The new strategies discussed here lead to improvement in both of these areas. Two areas of special focus in the present paper are (i) the finite difference boundary value method used for calculating distorted wave Green's functions and regular solutions for scattering by the distortion potential and (ii) the choice of the distortion potential itself. Among other results included here is the first application of the outgoing wave or scattered wave variational principle to reactive scattering.     

有限线性扩散区域内EC＇反应催化电流

有限线性扩散区域内ＥＣ＇反应催化电流吴志斌，董绍俊，吴浩青（同济大学化学系，上海，２０００９２）（中国科学院长春应用化学研究所）（复旦大学化学系）关键词有限线性扩散空间，ＥＣ＇反应，催化电流化学修饰电极等有关膜电极及薄层池的电极过程动力学问题均涉及到...     

随机模糊集与随机集

本文研究了三个方面的工作：一是定义了一种模糊集上的可测结构，从而定义了随机模糊集，这些定义都与论域Ｘ上的拓扑结构无关。将通常意义下的集合看成特殊模糊集得到的通常集合上的超可测结构与文（３）中的定义一致；二是给出了随机模糊集、随机集的一些等价条件；三是研究了随机模糊集、随机集的分布与其有限维落影族的关系。     

基于双尺度渐近分析的有限元算法

1．引言正如文山所说，由于复合材料和周期结构的材料系数ail（x）在局部区域内间断且跳跃性很大，加上区域内含有周期性洞穴或裂缝，且周期长度很小．一般而言，直接采用有限元方法进行数值模拟，其计算量大得惊人，甚至难以实现．文山针对这种特征，提出了一种可计算的双尺度渐近分析模式，本文在此基础上给出了相应的有限元算法，它包括：1．周期解在一个基本构造上的有限元计算；2．边界层的有限元计算．同时，给出了相应的误差分析．2．周期解的有限元计算首先考虑下列形式的边值问题；其中把，代E尸（on叫，iii（0关于E—（EI，ZZ…     

束变换环孔激光谐振腔的环束场传输差分计算

 通过使用新的坐标变量和光场的表达形式，推导了新坐标系下环束场传输的偏微分方程，给出了相应差分方程的截断误差和稳定性条件，使得环束场传输计算速度大幅提高；使用二阶精度的Crank-Nicolso差分方法对束变换环孔激光谐振腔的环束场传输进行了计算，并与W.D.Murphy和M.L.Bernabe所给M-B积分近似方法的计算结果进行了比较，证实了该方法无论计算速度和精度都优于M-B积分近似方法。     

脉冲泄流时土壤条件对冲击接地阻抗的影响

 为了探寻电流脉冲泄流时，设备的接地是否有效可靠，用FDTD法分析了脉冲泄流时均一土壤、分层土壤及局部改善土壤的电导率对冲击接地阻抗（TGR）的影响。结果表明：大地电导率对冲击接地阻抗的影响很大；对均一土壤，TGR随大地电导率的增大而减小，大地电导率越大，TGR对接地体的埋深越不敏感。对分层土壤，接地体要埋在大地电导率大的土壤层中，且不要靠近交界面；局部改善大地导电性能对降低TGR效果明显，TGR对改善区域的体积非常敏感，而当电导率大于一定值时，大地电导率的再增加，对TGR的影响不明显。     

Multi‐stage high order semi‐Lagrangian schemes for incompressible flows in Cartesian geometries

Efficient transport algorithms are essential to the numerical resolution of incompressible fluid‐flow problems. Semi‐Lagrangian methods are widely used in grid based methods to achieve this aim. The accuracy of the interpolation strategy then determines the properties of the scheme. We introduce a simple multi‐stage procedure, which can easily be used to increase the order of accuracy of a code based on multilinear interpolations. This approach is an extension of a corrective algorithm introduced by Dupont & Liu (2003, 2007). This multi‐stage procedure can be easily implemented in existing parallel codes using a domain decomposition strategy, as the communication pattern is identical to that of the multilinear scheme. We show how a combination of a forward and backward error correction can provide a third‐order accurate scheme, thus significantly reducing diffusive effects while retaining a non‐dispersive leading error term. Copyright © 2016 John Wiley & Sons, Ltd.