轴流泵端壁区域流动三维粘性数值计算

本文采用通用流场分析软件FLUENT,基于N-S方程,选用RNG κ—ε湍流模型与SIMPLEC算法,对轴流泵叶轮内部(及端壁间隙)流动进行了三维粘性数值计算。通过实验验证,表明数值计算结果和实测数据吻合较好。详细分析了叶顶附近流场形态以及叶轮出口轴向、周向速度分布,进行了性能预估。并给出了实验装置图。     

相对论二极管电子枪设计的分析解与实验的比较

本文对轴对称场二极管区电子运动所遵循的近似轨迹方程进行数值求解,以此做为设计相对论二极管电子枪的依据,与实验结果相比较,符合较好。     

多排叶片内部流动以及叶片排与排气系统相互作用研究

本文采用高分辨率格式和多块多网格方法求解雷诺平均Navier-Stokes方程,数值模拟多排叶片内的三维粘性流动;使用数值激盘模拟叶片排的存在,研究叶片排与非轴对称排气部件之间的相互作用以及复杂的内部流动。文中描述了数值方法,给出了NASA透平导叶和单级透平内部流场的数值结果及其与相应实验结果的对比,也给出了多级透平内部流场的数值结果,以及透平与排气部件之间的耦合流场的结果。     

BEPCⅡ束流相关本底的蒙特卡罗模拟研究(Ⅱ)—束流-气体相互作用

介绍了用蒙特卡罗方法研究正常对撞模式下束流本底中的束流-气体相互作用的原理与效果.使用通用的模拟工具,研究了BEPCⅡ的束流-气体本底对设计中的BESⅢ的影响.结果显示,在储存环上设置挡板后,探测器的本底在允许的范围之内.模拟结果对BEPCⅡ-BESⅢ的设计建造具有指导意义.     

Speed of adaptation in structured populations

Adaptation of populations takes place with the occurrence and subsequent fixation of mutations that confer some selective advantage to the individuals which acquire it. For this reason, the study of the process of fixation of advantageous mutations has a long history in the population genetics literature. Particularly, the previous investigations aimed to find out the main evolutionary forces affecting the strength of natural selection in the populations. In the current work, we investigate the dynamics of fixation of beneficial mutations in a subdivided population. The subpopulations (demes) can exchange migrants among their neighbors, in a migration network which is assumed to have either a random graph or a scale-free topology. We have observed that the migration rate drastically affects the dynamics of mutation fixation, despite of the fact that the probability of fixation is invariant on the migration rate, accordingly to Maruyama's conjecture. In addition, we have noticed a topological dependence of the adaptive evolution of the population when clonal interference becomes effective.     

热动力学问题的数值计算

 对含热传导的流体动力学方程组，用有限元方法进行数值求解。采用傅里叶热传导计算热流、用热流连续条件计算单元间接触面的温度、用三角形传热法计算体单元表面的热流，考虑各向同性弹塑性流体材料模型、三项式物态方程和导热系数与状态的相关性，给出了傅里叶热传导、接触传热、热应力应变效应、以及混合物冲击压缩特性等算例。对混合物冲击温度的数值模拟表明，小颗粒混合物在冲击压缩过程中，颗粒间的温度有差别、稍有波动，并随时间趋向于一致，以至热平衡。     

The influence of temperature gradients on the kinetic boundary layer problem for a condensing droplet

We extend an earlier method for solving kinetic boundary layer problems to the case of particles moving in aspatially inhomogeneous background. The method is developed for a gas mixture containing a supersaturated vapor and a light carrier gas from which a small droplet condenses. The release of heat of condensation causes a temperature difference between droplet and gas in the quasistationary state; the kinetic equation describing the vapor is the stationary Klein-Kramers equation for Brownian particles diffusing in a temperature gradient. By means of an expansion in Burnett functions, this equation is transformed into a set of coupled algebrodifferential equations. By numerical integration we construct fundamental solutions of this equation that are subsequently combined linearly to fulfill appropriate mesoscopic boundary conditions for particles leaving the droplet surface. In view of the intrinsic numerical instability of the system of equations, a novel procedure is developed to remove the admixture of fast growing solutions to the solutions of interest. The procedure is tested for a few model problems and then applied to a slightly simplified condensation problem with parameters corresponding to the condensation of mercury in a background of neon. The effects of thermal gradients and thermodiffusion on the growth rate of the droplet are small (of the order of 1%), but well outside of the margin of error of the method.     

自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。     

Thermal helium clusters at 3.2 Kelvin in classical and semiclassical simulations

The thermodynamic stability of⁴He_4–13 at 3.2 K is investigated with the classical Monte Carlo method, with the semiclassical path-integral Monte Carlo (PIMC) method, and with the semiclassical all-order many-body method. In the all-order many-body simulation the dipole-dipole approximation including short-range correction is used. The resulting stability plots are discussed and related to recent TOF experiments by Stephens and King. It is found that with classical Monte Carlo of course the characteristics of the measured mass spectrum cannot be resolved. With PIMC, switching on more and more quantum mechanics. by raising the number of virtual time steps results in more structure in the stability plot, but this did not lead to sufficient agreement with the TOF experiment. Only the all-order many-body method resolved the characteristic structures of the measured mass spectrum, including magic numbers. The result shows the influence of quantum statistics and quantum mechanics on the stability of small neutral helium clusters.     

A parallel inexact newton method for stochastic programs with recourse

A parallel inexact Newton method with a line search is proposed for two-stage quadratic stochastic programs with recourse. A lattice rule is used for the numerical evaluation of multi-dimensional integrals, and a parallel iterative method is used to solve the quadratic programming subproblems. Although the objective only has a locally Lipschitz gradient, global convergence and local superlinear convergence of the method are established. Furthermore, the method provides an error estimate which does not require much extra computation. The performance of the method is illustrated on a CM5 parallel computer.This work was supported by the Australian Research Council and the numerical experiments were done on the Sydney Regional Centre for Parallel Computing CM5.