Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

聚丙烯孔洞铁电驻极体膜的电极化及其电荷动态特性

根据栅控恒压电晕充电组合反极性电晕补偿充电法的实验结果计算出铁电驻极体的极化强度.结果说明，伴随着薄膜内孔洞气体的Paschen击穿，该铁电体的极化强度随栅压增加而显著上升.利用上述充电方法和热刺激放电(TSD)谱的分析讨论了这类空间电荷型宏观电偶极子，及与其补偿的空间电荷热退极化的电荷动态特性；阐明了这两类俘获电荷的能阱分布，即构成宏观电偶极子的位于孔洞上下介质层内的等值异号空间电荷分别被俘获在深、浅两种能值陷阱内，而位于薄膜表面层的注入空间电荷则被俘获在中等能值陷阱中. 关键词： 反极性电晕补偿充电法 铁电驻极体 充电电流 热刺激放电     

如何提高课堂教学效果

就如何提高课堂教学效果做了一些探讨。     

A method for calculating invariant subspaces of symmetric hyperbolic equations

An algorithm is constructed for calculating invariant subspaces of symmetric hyperbolic systems arising in electromagnetic, acoustic, and elasticity problems. Discrete approximations are calculated for subspaces that correspond to minimal eigenvalues and smooth eigenfunctions. Difficulties related to the presence of an infinite-dimensional kernel in the differential operator are successfully handled. The efficiency of the algorithm is demonstrated using acoustics equations.     

Comparison of growth stress measurements with modelling in thin iron oxide films

High temperature oxidation of metals leads to residual stresses both in the metal and in the growing oxide. In this work, the evolution of this residual stresses is theoretically predicted in the growing oxide layers. The origin of these stresses is based on a microstructural model. Using experimental results providing from the oxidation kinetics, and an analysis proposed to describe the growth strain occurring in the thin layers, a set of equations is established allowing determining the stresses evolution with oxidation time. Then, the model is compared with experimental results obtained on both α-Fe and phosphated α-Fe, oxidised at different temperatures. Numerical data are extracted from experiments either with an asymptotic formulation or with an inverse method. These two methods give good agreement with experiments and allow extracting the model parameters.     

有限时滞微分系统依照双测度的实际稳定性

本文对具有有限时滞的泛函微分方程建立了关于依照两种测度的实际稳定性的Razumikhin型判定定理，其中未采用通常的辅助函数，且可运用多个含有状态变量x的部分变元的Lyapunov函数，得出部分变远实际稳定性的判定定理，从而改进了已有的结果。     

关于圈的完全正矩阵实现

本文给出了一个关联图为圈的非负、半正定矩阵A为完全正的一个充要条件.我们还证明了这样的矩阵A（当A为完全正时）的分解指数即为A的阶数.     

一维光学格子孤子的传输特性及控制研究

利用解析和数值方法研究了在具有横向折射率周期性调制的克尔型非线性介质中光学格子孤子的传输，得到了孤子参数的演化方程以及格子孤子的形成和稳定传输的条件.结果表明：当光束的入射角小于某临界角度时，光束可被类似波导形式的路径俘获而稳定传输，该临界角随折射率调制周期、调制深度的增加而增大，且光束越窄临界值越大.此外，线性空间啁啾虽然对光束传输的中心位置没有任何影响，但会导致光束发散从而破坏格子孤子的形成和稳定传输，对此提出了采用特定功率取值来补偿啁啾作用从而形成格子孤子的方案. 关键词： 光孤子 光学格子 光传输 矩方法     

宽带膜厚实时监控过程中膜层折射率的确定方法

为了准确计算出镀膜过程中每层膜的折射率，介绍了实时监控过程中确定膜层折射率的2种方法：一种是由实测的透射比光谱直接反算出膜层的折射率；另一种是用最小二乘法的优化算法实时拟合折射率。试验结果表明：在线反算适合单点监控，所得折射率误差小于2%。然而在实际镀膜过程中,由于宽带内膜层参数误差较大,一般大于25%。为此，采用最小二乘法拟合，即在整个宽光谱范围内采集每个波长点的信息，所得结果误差很小，一般都在2％～5%之间,有时可达到10%,在很大程度上提高了实际镀膜时膜厚监控的精度。     

A Modified Quasi-Newton Method for Structured Optimization with Partial Information on the Hessian

This paper develops a modified quasi-Newton method for structured unconstrained optimization with partial information on the Hessian, based on a better approximation to the Hessian in current search direction. The new approximation is decided by both function values and gradients at the last two iterations unlike the original one which only uses the gradients at the last two iterations. The modified method owns local and superlinear convergence. Numerical experiments show that the proposed method is encouraging comparing with the methods proposed in [4] for structured unconstrained optimization Presented at the 6th International Conference on Optimization: Techniques and Applications, Ballarat, Australia, December 9–11, 2004