辐射传递蒙特卡洛法精度分析及数值试验

本文建立了蒙特卡洛法模拟散射参与性介质内辐射传递计算模型。对蒙特卡洛法的计算精度及运行时间进行了较为详细的分析,提出了几种判断计算精度的方法。同时,借助蒙特卡洛法模拟辐射传递过程,进行数值“辐射实验”。利用该“实验结果”进行了物性反问题研究。在已知光学厚度的前提下,得到散射反照率与后半球辐射热流之间的单值函数关系。     

Numerical Prediction of the Strain State of Thermoset Matrix/Mineral Filler Composite Plates

The numerical prediction of the fields of inelastic strains (the linear invariant of the tensor of inelastic strains) in thermoset polyester/marble filler composite plates is discussed. A uniformly distributed load is applied to the plates, which lie on a steel base. The strain fields are predicted by means of the boundary element method by using creep test data for the composites and the polyester matrix itself. Identical creep tests were performed for two ages of the materials (1 month and 13 years), which allowed evaluating the aging effect. The study is carried out in two stages. At the first stage, the application of the generalized Maxwell-Gurevich equation to the thermoset matrix/mineral filler composite is demonstrated. The model parameters determined from the experimental creep data is used for the second stage, where the state of inelastic strains in the plates is predicted by applying the boundary element method. The influence of composite formulation (filler content) and physical aging of the polyester matrix on the state of inelastic strains in the plates is shown.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 145–156, March–April, 2005.     

基于MMD算法的HMM参数重估计的矩阵表示

本文提出的MMD算法用于提高模型区别错误信息和正确信息的能力．利用该算法在对模型的参数进行重估计时．涉及到复杂的目标函数的梯度运算．击运用矩阵运算使得梯度运算变得简单明了，因此本文给出了MMD算法下的HMM参数重估计的矩阵表示形式并给出了证明．     

An Additive Problem with Primes in Arithmetic Progressions

In this paper, we extend a classical result of Hua to arithmetic progressions with large moduli. The result implies the Linnik Theorem on the least prime in an arithmetic progression.     

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.     

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.     

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

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

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法的改进越显著。