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

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

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

HgCl2/I2复合催化剂协同催化合成异丙醇钇的研究

以金属钇和异丙醇为原料，以HgCl2/I2为复合催化剂，通过对金属钇的机械加工以增加其比表面，并将异丙醇脱水使其含水量降低至0．05％，体系在82℃回流5h，经过滤、减压蒸馏，得到了白色海绵状异丙醇钇，其产率高达83％，合成时间比文献报道的缩短了19h，产率提高了8％。文章确定了催化剂的最佳用量为20gY加入60mg HgCl2/I2，研究了合成产率与HgCl2/I2催化剂和HgCl2催化剂的依赖关系及异丙醇中含水量对合成产率的影响，并对HgCl2/I2的催化作用机理进行了初步探讨。     

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.     

Remote stereocontrol by the sulfinyl group: Mukaiyama aldol reactions of (S)-2-[2-(p-tolylsulfinyl)phenyl]acetaldehyde in the asymmetric synthesis of β-hydroxyacids and 1,3-diols

(S)-2-[2-(p-Tolylsulfinyl)phenyl]acetaldehyde reacts with different O-silylated ketenethioacetals in the presence of Yb(OTf)₃ yielding β-hydroxythioesters in high yields and diastereoselectivities. The obtained compounds were readily transformed into β-hydroxyacids and their corresponding diols. These Mukaiyama aldol reactions are a direct evidence of the ability of the sulfinyl group to control 1,5- and 1,6-asymmetric induction processes.     

A study of structural, compositional and optical properties of spray-deposited non-stoichiometric (Zn,Fe)S thin films

Non-stoichiometric ternary chalcogenides (Zn,Fe)S were prepared in the film form by pyrolytic spray deposition technique, using air/nitrogen as the carrier gas. The precursor solution comprised of ZnCl₂, FeCl₂ and thiourea. The depositions were carried out under optimum conditions of experimental parameters viz. carrier gas (air/nitrogen) flow rate, concentration of precursor constituents, nozzle substrate distance and temperature of quartz substrate. The deposited thin films were later sintered in argon at 1073 K for 120 min.The structural, compositional and optical properties of the sintered thin films were studied. X-ray diffraction studies of the thin films indicated the presence of (Zn,Fe)S solid solution with prominent cubic sphalerite phase while surface morphology as determined by scanning electron microscopy (SEM) revealed a granular structure.The chemical composition of the resulting thin films as analyzed by energy dispersive X-ray analysis (EDAX) reflected the composition of the precursor solutions from which the depositions were carried out with Fe at% values ranging from 0.4 up to 33.SEM micrographs of thin films reveal that the grain sizes of the thin films prepared using air as carrier gas and N₂ as carrier gas are in the vicinity of 300 and 150 nm, respectively.The diffuse transmittance measurements for thin films, as a function of wavelength reveal the dependence of direct optical band gap on Fe content and type of phase.     

Convergent Flows of Molten Polymers Modeled by Generalized Second-Grade Fluids of Power-Law Type

The molding processes of polymer melts involve geometrically complex dies. Such dies are usually tapered or streamlined to achieve a maximum output rate under conditions of laminar flow. The model of a generalized second-grade fluid of power-law type is used and the results obtained are illustrated by examples of convergent flows in conical and wedge-shaped dies.