估计死亡率分布的一个最大熵模型

本文提出了一种估计死亡率分布的新模型一最大熵模型。该模型直接从样本信息出发，不需要对待估分布的概率密度函数或先验分布作任何假定，从而克服了极大似然估计和贝叶斯估计的不足。而且通过两个例子的计算结果，表明该方法与样本数据的拟合效果要好于其它两种方法。     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Calculation of Normal Vibrations of Chlorin by the Method of Density Functionals

Using the DFT/B3LYP method with a 6-31G(d) basis set, the structure, normal vibration frequencies, and the absolute band intensities in the IR spectra of the chlorin molecule and its four symmetric isotopomers have been calculated. Scaling of the force field by the Pulay method in independent and natural coordinates has been carried out. A method for obtaining effective force fields without using experimental data on the fundamental vibration frequencies is proposed. By comparing the vibration modes and constructing special matrices, complete assignment of the fundamental frequencies of porphin and chlorin has been carried out. It has been shown that the majority of porphin macroring vibrations upon pyrrolenine ring hydrogenation are frequency-characteristic and only 12 vibrations change considerably. A frequency correlation with regard for the mode transition between chlorin and all its isotopomers under consideration has been established. Comparative analysis of the force fields of porphin and chlorin in dependent natural coordinates has revealed the unique nonlocal character of the change in force constants of the macroring upon hydrogenation of one pyrrolenine ring. Modeling of the IR spectra of chlorin and its isotopomers has been performed. Assignment and interpretation of the normal vibrations of the molecules under consideration have been carried out.     

对流扩散方程的一种新型差分格式

对流扩散方程可以描述众多的物理化学现象，因而对其寻求稳定的，实用的数值解法有着重要的现实意义。本文针对形式较一般的一维非定常对流扩散方程，构造了对角元严格占优的Crank-Nicholson差分格式，然后对其分别用分离变量的方法以及能量估计的方法作了稳定性的分析，最后给出了数值试验的结果，数值结果表明本文构造的格式能够较好的处理经典的Crank-Nicholson格式所不能处理的对流项系数较大的对流扩散方程，并具有较好的精度。     

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.     

Essay: The Quantum Gravitational Black Hole Is Neither Black Nor White

Understanding the end state of black hole evaporation, the microscopic origin of black hole entropy, the information loss paradox, and the nature of the singularity arising in gravitational collapse - these are outstanding challenges for any candidate quantum theory of gravity. Recently, a midisuperspace model of quantum gravitational collapse has been solved using a lattice regularization scheme. It is shown that the mass of an eternal black hole follows the Bekenstein spectrum, and a related argument provides a fairly accurate estimate of the entropy. The solution also describes a quantized mass-energy distribution around a central black hole, which in the WKB approximation, is precisely Hawking radiation. The leading quantum gravitational correction makes the spectrum non-thermal, thus providing a plausible resolution of the information loss problem.     

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过程的影响,并与相应的离体实验结果进行了对比。数值模拟结果表明,组织的热物性及血液灌注率的动态变化对于热损伤体积的变化具有重要的影响。因此在激光诱导间质肿瘤热疗的数值模拟中应该考虑热物性及血液灌注率的动态变化以期为临床治疗方案的制定提供更为准确的依据。