Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

Highlights and Prospects of Hyperspherical Approach in Few-Body Quantum-Mechanical Systems

The binding energy of excitonium negative ion for ground ¹S-state in bulk conductors: Ge, Si, CdSe and for green and yellow lines of Cu₂O in hyperspherical coordinate method are found. Angular and radial correlations between electrons are taken into account by the channel functions, which are the eigen-functions of operator on the surface of sphere in six-dimensional sphere. The calculation of energies have been done using the adiabatic and Born-Oppenheimer approximations. In Born-Oppenheimer approximation is enough to give a binding energy with only 1.2% error, in adiabatic approximation this error drops to 0.7 %.     

解线性最小二乘问题的一个新并行算法

讨论了求解无约束线性最小二乘问题的一种并行单纯形法以及对它的改进算法并行共轭梯度—单纯形法 .算法本身具有很强的并行机制 ,能够充分地发挥并行机快速省时的特点 .本文也对算法做了理论分析 ,对算法的收敛性给予了证明 (在二维情形下 ) .最后做了数值实验 (由于软硬件条件的限制 ,并行算法未能在并行计算机上实现 ,鉴于这种情况 ,我们所做的数值实验均是在串行机上完成的 )     

Oscillation theorems for differential equations with a nonlinear damping term

In this paper, we are concerned with a class of nonlinear second-order differential equations with a nonlinear damping term. Passage to more general class of equations allows us to remove a restrictive condition usually imposed on the nonlinearity, and, as a consequence, our results apply to wider classes of nonlinear differential equations. Two illustrative examples are considered.     

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

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

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