圆柱坐标系下任意方向辐射强度的源项六流法模拟

基于传统热流法,提出一种圆柱坐标系下的源项六流模型(Source Six Flux,SSF),可快速准确地计算参与性介质内任意方向的出射辐射强度.详细介绍SSF模型的基本原理和求解步骤,以圆柱形吸收、散射、发射性介质为例,模拟其沿任意方向的出射辐射强度,并与反向蒙特卡罗法(Backward Monte Carlo,BMC)和二流法(Two Flux Method,TFM)的计算结果进行比较.结果表明,SSF法与BMC法的计算结果吻合较好,计算精度均高于TFM法,但SSF法的计算效率明显优于BMC法.因此,SSF模型是一种适用于计算任意方向辐射强度问题的高效数值模型.     

衬底温度对电子辅助增强化学气相沉积金刚石膜的影响

本工作采用蒙特卡罗(Monte Carlo)方法,对以CH4/H2为源气体的电子助进化学气相沉积(EACVD)金刚石薄膜的气相动力学过程进行了模拟.提出了衬底温度的空间梯度变化模型,研究了衬底温度对EACVD气相过程中的电子群行为以及H2和CH4分解过程的影响.结果表明:电子平均温度随着衬底温度的升高而升高;当气压较低时,分解得到的成膜关键粒子H、CH3数量随衬底温度的增大而减少;而当气压较高时,H、CH3数目随衬底温度的增大而增大;衬底温度主要改变了衬底表面附近的化学反应动力学过程,从而对薄膜质量产生了决定性的影响.     

模体厚度对MC程序计算速度与效率的影响研究

基于蒙特卡罗(MC)程序MCNP4c, DOSXYZnrc, DOSRZnrc 和DPM, 考察了模体厚度变化对MC程序计算速度和MC效率的影响, 发现模体的厚度变化对不同MC程序速度提升的影响差异较大, 其中DOSXYZnrc 提升最多, DOSRZnrc提升最小。 从MC效率方面来讲, DPM和DOSRZnrc具有明显的优势。 就4个程序而言, 薄模体均表现出一定的MC 效率优势。 从MC 效率提升角度来看, MCNP4c, DOSXYZnrc和DOSRZnrc表现出一定的相似性, 而DPM相对低些。 也对DOSRZnrc异常的统计不确定性特征、 速度提升进行了研究, 发现其内部植入的光子强迫方差减小技巧在降低体元统计不确定性的同时, 也削弱了薄模体的速度优势。 但是由于显著影响了方差, 所以总体来讲, 仍然提高了DOSRZnrc薄模体的MC 模拟效率。     

二层Ising模型的短时临界动力学性质

采用动力学Monte Carlo 方法研究了二层Ising模型的临界性质及早期动力学标度行为.结果表明层间耦合不为零时也存在临界点;计算了早期动力学临界指数θ;估计了传统的临界指数1/νz.其结果支持临界线存在的猜想,并表明此模型很可能是一种弱普适模型. 关键词：     

ZnS型薄膜电致发光器件输运过程的解析能带模拟

采用解析能带模型对ZnS型薄膜电致发光器件的输运过程进行了Monte Carlo模拟.利用分段的多项式拟合了ZnS的实际导带结构,计算了能态密度和散射速率.以这些结果为基础,模拟了ZnS中的电子输运过程.通过比较,这一模型既具有非抛物多能谷能带模型运算速度快、使 用方便的优势,又具有与采用全导带模型相近的计算精度.进而,讨论了能带模型对模拟结 果的影响,发现色散关系与能态密度均对模拟结果有较大影响.     

计算复杂边界渗流的蒙特卡罗方法

本文运用不规则游动网格的蒙特卡罗方法求解复杂边界渗流问题。与以往各种数值方法相比,蒙特卡罗方法在计算渗流问题时要更灵活、更方便,可以根据需要独立地计算出渗流区域内任意一点的渗透压力、渗流流速。     

氯化钠水溶液的monte Carlo分子模拟研究

采用MonteCarlo计算机分子模拟方法,研究了氯化钠水溶液在常温和高温情况下粒子的径向分布函数和粒子间各种作用位能的热力学性质。模拟过程采用了NTV正则系综,粒子间的作用能包括离子库仑静电作用,偶极子作用以及色散作用,这些作用构成了电解质溶液的基本框架,模拟结果与微扰理论和平均球近似积分方程理论的预测值进行了比较。     

Electromagnetic Scattering from One-Dimensional Time-Varying Sea Surface with Gaussian-Like Beam Incidence

A widely used iterative technique named the Method of Ordered Multiple Interactions is given for calculating the Gaussian-like beam scattering from a time-varying sea surface with the Pierson-Moskowitz (P-M) spectrum. This is done by solving the magnetic field integral equation for the current induced on an infinite rough surface. Following the discretization of the integral equation, the unknown currents can be determined more rapidly with the LU decomposition. Numerical results are presented with emphasis on the electromagnetic backscattering at low grazing angle incidence. It is shown that the backscattering cross section is proportional to the nearly fourth power of the grazing angle for the plane and beam incidence. This is consistent with the result given in some references. The angular distributions of the backscattering cross section with a different beam waist, surface length, and velocity of the wind are discussed.     

Folding of lattice protein model chains with fixed ends

Using Monte Carlo simulations, we have studied the folding dynamics and thermodynamics of geometricall yconstrained lattice protein model chains. The constraints are realized by fixing one or both terminals of the chains. By comparing the results with that of the free-end chains, we find that the folding behaviours of the end-constralned chains are not completely similar to that of the free-end chains. Both kinds of constraints on the chain ends affect the folding dynamics of the chains: i.e., the folding rate, but not the thermodynamics. The thermodynamic behaviour of the one-end-fixed chains shows less difference from that of the free-end chains, while the thermodynamic behaviour of the two-end-fixed chains has obvious difference from that of the free-end chains. The origin of these differences comes from the differences of the ergodicitv of the chains in the conformational space.     

氢原子电子云的蒙特卡罗模拟

根据量子力学原理,氢原子中电子的几率分布由电子的定态波函数可以唯一地确定.在大多数教科书中,通常是将电子沿径向和角度方向的几率分布分别进行讨论,因此很难给出电子在整个空间几率分布的一个完整图象.我们使用蒙特卡罗随机抽样的方法,通过计算机作图形象而逼真地绘制了氢原子中电子的几率分布图,即所谓的电子云.下面介绍所使用的方法.