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)的一些物理数据。     

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.     

曙光一号自由电子激光器的理论计算

系统总结曙光一号自由电子激光器理论计算的主要结果:包括曙光一号装置主要参数的选取和理解;磁场失谐曲线的计算;常参数摇摆器和变参数摇摆器的主要结果;高阶波导模的贡献;电子束参数扰动对激光性能的影响;空间电荷效应等。计算结果表明,常参数摇摆器激光输出功率可达80MW,效率约50％;变参数摇摆器激光输出功率可达250MW左右,效率约16％。     

Numerical solution of eigenvalue problems for Hamiltonian systems

This paper discusses the numerical solution of eigenvalue problems for Hamiltonian systems of ordinary differential equations. Two new codes are presented which incorporate the algorithms described here; to the best of the author’s knowledge, these are the first codes capable of solving numerically such general eigenvalue problems. One of these implements a new new method of solving a differential equation whose solution is a unitary matrix. Both codes are fully documented and are written inPfort-verifiedFortran 77, and will be available in netlib/aicm/sl11f and netlib/aicm/sl12f.     

Numerical propagation of partially coherent laser beams through optical systems

The propagation of partially coherent beams through optical systems is computed numerically in one transverse dimension. The optical system is divided into different elementary segments, through which the propagation of light can be calculated by appropriate operators, working on the coherence function or the Wigner distribution function respectively. For the necessary changes between these two functions describing the partially coherent beams, the use of the remarkable z-transform seems to be an advantage. With this algorithm the grid and the resolution in the spatial frequency domain can be arbitrarily chosen in contrast to the usual Fourier transform, the influence of phase aberrations on the focusability of Gauss-Schell model beams is discussed as an application example of the numerical model. With the help of this tool, practical beam guiding systems can be simulated for use with multimode laser radiation.     

斜孔气膜冷却数值模拟分析

本文通过数值模拟分析圆柱孔和扩散孔的单斜孔气膜冷却特性,考察复合角、孔型和吹风比对流场和气膜冷却效果的影响。结果表明,复合角的引入使气膜侧向分布更宽,但冷却效果沿主流方向衰减更快。适中吹风比得到的气膜能更有效的保护壁面。在相同吹风比和复合角条件下,扩散孔的气膜冷却效率比圆柱孔更好,且冷却更为均匀持久。     

Numerical Computation of Figure-eight Solution for 3-body Problems

The main goal of this paper is to compute the Figure-eight solutions for the planar Newtonian 3-body problem with equal masses by finding the critical points of the functional associated with the motion equations of 3-body in plane R2. The algorithm adopted here is the steepest descent method, which is simple but very valid for our problem.     

含腔电大尺寸导体目标电磁散射的一体化数值模拟

应用等效原理，通过引入口面上等效磁流将含腔导电目标电磁散射简化为腔内、外两个等效 问题. 腔内问题分段求解并应用级联法获得口面等效导纳矩阵；腔内外的耦合关系应用近似 边界元方法描述并由此获得口面等效磁流；最后，这一具有混合源的腔体内外一体化散射问 题则应用所提出的广义混合场积分方程方法建立电磁模型，并用多层快速多极子方法实现高 效数值求解. 实例计算结果与测试结果具有很好的一致性. 关键词： 含腔目标 电磁散射 混合场积分方程 数值分析     

Runge–Kutta Methods for Itô Stochastic Differential Equations with Scalar Noise

A general class of stochastic Runge–Kutta methods for Itô stochastic differential equation systems w.r.t. a one-dimensional Wiener process is introduced. The colored rooted tree analysis is applied to derive conditions for the coefficients of the stochastic Runge–Kutta method assuring convergence in the weak sense with a prescribed order. Some coefficients for new stochastic Runge–Kutta schemes of order two are calculated explicitly and a simulation study reveals their good performance.