用格子Boltzmann模型模拟非等温流场

根据微观和宏观之间的质量、动量、能量守恒准则和在原格子Boltzmann模型基础上,建立了几个新的格子Boltzmann模型,使得在外力场中的格子Boltzmann模型得到进一步完善.通过还原宏观流体力学方程,捕捉到了浮力强迫系数与Grashof数之间的关系.所得动量方程和Navier Stokes方程相比,在黏性输运项上有明显的改进,说明黏性应力不但与流体的速度梯度和流体的压缩性有关,而且还与非定常的内能梯度和动量通量有关.该模型对非等温流场的数值结果证明了其具有很好的数值稳定性和适用性. 关键词： Boltzmann模型 平衡分布函数 流体力学方程     

用基于Daubechies尺度函数的时域多分辨分析计算电磁散射

利用Daubechies尺度函数的紧支撑性，将一种新的时域多分辨分析(MRTD)算法应用到三维目标的电磁散射中并分析了其色散特性，使用MRTD/FDTD接口技术解决了连接和吸收边界的处理.算例表明，在不牺牲精度的情况下，该方法比传统的FDTD需要的网格数少，计算速度快.最后，计算了二维PBG结构的反射特性. 关键词： 时域多分辨分析 Daubechies尺度函数 电磁散射     

无穷积分余项与其被积函数比值的收敛阶

对一类函数的无穷积分余项与该函数的比值得到当x趋于无穷大时的收敛阶,这类函数是幂函数与指数函数的乘积函数,并将其应用到Mittag-Leffler函数.同时考虑了对应的级数情形.     

O2分子激发态a1Δg和b1Σ+g态的分析势能函数

以aug-cc-pVQZ,cc-pV5Z,6-311++g（d,p）和6-311++g（3df,3pd）等为基函数,采用多组态相互作用（MRCI）方法对O₂分子最低的两个激发态¹Δ_g和¹Σ⁺_g的平衡结构、谐振频率和势能曲线进行了计算.并选用Murrell-Sorbie势能函数对曲线进行拟合,利用拟合的参数值计算出了力常数和光谱数据.结果表明计算值与实验值符合较好. 关键词： MRCI 势能函数 力常数 光谱数据     

非线性波方程求解的新方法

从Legendre椭圆积分和Jacobi椭圆函数的定义出发，得到了新的变换，并把它用于非线性演化方程的求解.用三个具体的例子，如非线性Klein-Gordon方程、Boussinesq方程和耦合的mKdV方程组，说明了具体的求解步骤.比较方便地得到非线性演化方程或方程组的新解析解，如周期解、孤子解等. 关键词： Jacobi椭圆函数 非线性方程 周期解 孤子解     

激光辐照下皮肤组织光热响应有限元分析

为了更好地选择临床激光医疗曝光参量,采用有限元数值计算方法,模拟了脉冲激光与连续激光对人皮肤组织的光热作用及导致的温度变化效应,比较了两者的不同,得到了热响应时间及热弛豫时间与组织深度的关系,即组织越深(0～60 μm),其热响应时间(0～4 ms)与热弛豫时间(0.4～12.1 ms)越长;分析了激光脉宽长短对组织升温的影响;建立了评价脉冲间热损失的评价函数δ,并以此对脉冲间隔的选取作了探讨.     

质点从可自由移动的凹曲面上滑下的速度、加速度和时间

对于一个质点从可自由移动的光滑的凹曲面上滑下的速度、加速度和时间作了一般性讨论,建立了相应的公式.并对质点的运动平面与凹曲面的交线分别为直线、抛物线和椭圆等情况进行了具体的计算.     

Computer simulation of some dynamical properties of the Lorentz gas

We carried out molecular dynamics simulations of a Lorentz gas, consisting of a lone hydrogen molecule moving in a sea of stationary argon atoms. A Lennard-Jones form was assumed for the H₂-Ar potential. The calculations were performed at a reduced temperatureK ^* =kT/_{H 2–Ar} = 4.64 and at reduced densities ^*= _{Ar Ar} ³ in the range 0.074–0.414. The placement of Ar atoms was assumed to be random rather than dictated by equilibrium considerations. We followed the trajectories of many H₂ molecules, each of which is assigned in turn a velocity given by the Maxwell-Boltzmann distribution at the temperature of the simulation. Solving the equations of motion classically, we obtained the translational part of the incoherent dynamic structure factor for the H₂ molecule,S _tr(q, ). This was convoluted with the rotational structure factorS _rot(q, ) calculated assuming unhindered rotation to obtain the total structure factorS(q, ). Our results agree well with experimental data on this function obtained by Egelstaffet al. At the highest density ( ^*=0.414) we studied the dependence ofS(q, ) on system size (number of Ar atoms), number of H₂ molecules for which trajectories are generated, and the length of time over which these trajectories are followed.     

负氢离子源中多峰磁场形态数值模拟

 基于磁荷模型,在有限差分法基础上,推导了3维永磁体数值算法,结合全3维PIC/MCC粒子模拟算法,数值分析比较了外置型和“伞状”型两种特殊形态多峰磁场对电子能量沉积的影响。模拟结果表明:二者都能对粒子起到约束作用且能过滤引出负氢离子,同时二者电子能量分布规律基本一致,都呈现了双电子能态,符合等离子体放电基本机理;外置型过滤磁场对粒子的约束能力更强,能明显产生更多的粒子,总粒子数大约是“伞状”型过滤磁场情况下的4倍;“伞状”型过滤磁场能有效地抑制由磁场不均匀所引起的电子漂移,使产生的负氢离子空间分布更均匀。模拟结果与国外实验结果基本一致。     

Calculation of the reactor neutron time of flight spectrum by convolution technique

It is a very complex and time-consuming process to simulate the nuclear reactor neutron spectrum from the reactor core to the export channel by applying a Monte Carlo program. This paper presents a new method to calculate the neutron spectrum by using the convolution technique which considers the channel transportation as a linear system and the transportation scattering as the response function. It also applies Monte Carlo Neutron and Photon Transport Code （MCNP） to simulate the response function numerically. With the application of convolution technique to calculate the spectrum distribution from the core to the channel, the process is then much more convenient only with the simple numerical integral numeration. This saves computer time and reduces some trouble in re-writing of the MCNP program.