The method of lines solution of discrete ordinates method for radiative heat transfer in cylindrical enclosures

A radiation code based on method of lines solution of discrete ordinates method for radiative heat transfer in axisymmetric cylindrical enclosures containing absorbing-emitting medium was developed and tested for predictive accuracy by applying it to (i) test problems with black and grey walls (ii) a gas turbine combustor simulator enclosing a non-homogeneous absorbing-emitting medium and benchmarking its steady-state predictions against exact solutions and measurements. Comparisons show that it provides accurate solutions for radiative heat fluxes and can be used with confidence in conjunction with CFD codes based on the same approach.     

Comment on: “Confined states in two-dimensional flat elliptic quantum dots and elliptic quantum wires” [Physica E 11 (2001) 345]

We argue that the two-dimensional elliptic quantum dot problem with finite barrier cannot be exactly solved, contrary to a recent assertion (van den Broek and Peeters, Physica E 11 (2001) 345. We also prove it explicitly by numerically calculating the correct energy spectrum.     

重新初始化的LS方法跟踪二维可压缩多介质流界面运动

 用非耦合求解方法计算Level Set函数方程与流体力学方程组，应用重新初始化的Level Set函数确保距离函数性质，流体力学方程组采用二阶精度多介质流波传播差分格式计算，重新初始化方程采用五阶WENO格式计算。并给出了二维可压缩多介质流界面运动的计算结果。     

Density-Functional Theory Calculations of Normal Modes and Raman Intensities for Tetraazaporphin

The structure, harmonic frequencies, and nonresonance Raman intensities for porphin, tetraazaporphin (TAP), and three of its isotopomers are calculated by the density-functional theory of B3LYP/6-31G(d). Scaling of force constants for porphin in nonredundant natural coordinates is performed. The scaling factors obtained were used to predict the force field and normal modes of TAP and three of its isotopomers. Two alternative methods are used to carry out reliable assignment of the TAP frequencies: wavenumber-linear scaling method and frequency-shift method. There is good agreement between the frequencies predicted within the framework of the three methods used. The conservativeness of the out-of-plane B _2g - and B _3g -modes for porphin and TAP is examined. The Raman spectrum for TAP is simulated. A refinement of the assignment of the experimental frequencies for TAP of even symmetry types on the basis of the calculations performed is made.     

激光诱导间质肿瘤热疗的数值模拟和实验研究

本文在考虑生物组织物性动态变化的情况下建立了激光诱导间质肿瘤热疗(LITT)的物理数学模型,采用MonteCarlo方法数值模拟了LITT中激光能量在生物组织内的传输过程,基于Pennes生物传热方程和Arrhenius方程数值求解了组织内的温度分布和热损伤体积的变化,分析了热物性及血液灌注率的动态变化对LITT过程的影响,并与相应的离体实验结果进行了对比。数值模拟结果表明,组织的热物性及血液灌注率的动态变化对于热损伤体积的变化具有重要的影响。因此在激光诱导间质肿瘤热疗的数值模拟中应该考虑热物性及血液灌注率的动态变化以期为临床治疗方案的制定提供更为准确的依据。     

线性均衡约束最优化的一个广义投影强次可行方向法

本文讨论带线性均衡约束最优化问题,首先利用摄动技术和一个互补函数将问题等价转化为一般约束最优化问题,然后结合广义投影技术和强次可行方向法思想,建立了问题的一个新算法.算法在迭代过程中保证搜索方向不为零,从而使得每次迭代只需计算一次广义投影.在适当的条件下,证明了算法的全局收敛性,并对算法进行了初步的数值试验.     

三苯甲基自由基及其二聚体稳定性的理论研究

采用量子化学AM1 MO方法优化了三苯甲基自由基及其两种二聚体的构型，并计算了反应物和二聚物的生成焓、活化能以及它们的电荷布居，计算结果表明，三苯甲基自由基是一种较稳定的自由基，在这个自由基中存在明显的共轭效应，分子中电子云呈平均化分布．它的二聚体1，4—环己二烯衍生物分子中的5个苯基排布在环己二烯环平面的两端，明显降低了苯基之间的空间排斥力，使二聚物—环己二烯衍生物分子的稳定性大大高于六苯乙烷．因此，三苯甲基自由基的二聚体是1，4—环己二烯衍生物而不是六苯基乙烷，得到了与实验事实一致的结果。     

Detecting instabilities in flows of viscoelastic fluids

There is a growing interest in developing numerical tools to investigate the onset of physical instabilities observed in experiments involving viscoelastic flows, which is a difficult and challenging task as the simulations are very sensitive to numerical instabilities. Following a recent linear stability analysis carried out in order to better understand qualitatively the origin of numerical instabilities occurring in the simulation of flows viscoelastic fluids, the present paper considers a possible extension for more complex flows. This promising method could be applied to track instabilities in complex (i.e. essentially non‐parallel) flows. In addition, results related to transient growth mechanism indicate that it might be responsible for the development of numerical instabilities in the simulation of viscoelastic fluids. Copyright © 2003 John Wiley & Sons, Ltd.     

应用自动微分的Newton-PCG算法

一类新的使用符号微分的Newton-PCG型算法在文献[1]和[2]被导出来了。本文建立和研究应用自动微分的相应的Newton-PCG算法，理论分析和数值实验结果显示应用自动微分之后，目标函数的维数或复杂性越大，Newton-PCG算法对Newton法的改进越显著。     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.