高能碰撞各向异性动力学起伏的一个定量判据 *

通过对自仿射随机级联模型的Levy稳定性研究,发现其Levy稳定性指数不仅依赖于动力学起伏参数α而且依赖于统计阶乘矩时相空间的收缩方式.这样就提供了一种定量地判断动力学起伏是否各向异性的方法.另外,还从各向异性动力学起伏的假定出发,成功地解释了NA22实验组关于Levy稳定性破坏的实验结果.     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

表面杂质对NO+CO/Pt(100)反应体系振荡动力学的影响

提出了一个格气模型， 探讨催化表面的惰性杂质对NO+CO/Pt(100)反应体系振荡动力学行为的影响．研究发现表面杂质较少时，反应进程中表面重构能够形成相连通的1*1相和持续的振荡行为．当表面杂质不断增加时，反应进程中1*1相只能形成许多孤立的畴，从而在1*1相形成的随机的局域振荡的空间关联被削弱，因而体系全局的持续振荡行为演变为衰减振荡行为．当表面吸附的CO和NO的扩散速率增加时，局域振荡的空间协同又得到加强，体系又呈现出持续振荡行为．     

基于Fournier Forand体积散射函数的水中激光脉冲后向散射特性分析模型

基于Fournier Forand体积散射函数,建立了一种水中激光脉冲后向散射仿真模型.运用该模型可用Monte Carlo方法模拟光子在水中的传播过程,并得到光波后向散射的冲击响应.将初始激光脉冲与冲击响应进行卷积并求其傅里叶谱,即可得到激光脉冲后向散射信号的时域和频域特征.利用该模型分析了入射为高斯型激光脉冲时,水...     

激光等离子体相互作用中的电磁不稳定性和能量运输

利用相对论电磁粒子模拟程序研究了超强激光与等离子体相互作用中产生的电磁不稳定性的发展演化过程.讨论了电磁不稳定性的发生和非线性饱和过程.给出了不稳定性的线性增长率和各向异性参量之间的函数关系,用Spitzer-Harm理论分析了电子热传导中能量的运输情况,观察到由激光的非等方加热引起的电子纵向加热现象.结果表明,不稳定...     

含葡萄糖的混浊介质的Mueller矩阵模拟与分析

利用矢量Monte Carlo算法结合葡萄糖的旋光特性模拟了含糖混浊介质的背散射Mueller矩阵.数值模拟表明:含糖混浊介质的背散射Mueller矩阵元素中除Mll、M14、M41、M44外的其它元素花样均发生龙卷风状的不同程度的扭曲,其中M24、M34、M42、M43的扭曲程度更为显著.研究发现M24、M34、M4...     

MCNP Dose Calculations in a CT Phantom for Therapeutic External Photon Beam

In this paper, we have addressed the problem of the radiation transport with the Monte Carlo N particle(MCNP) code. This is a general purpose Monte Carlo tool designed to transport neutron, photon and electron in three dimensional geometries. To examine the performance of MCNP5 code in the field of external radiotherapy, we performed the modeling of an Electron Density phantom (EDP) irradiated by photons from 60Co source. The model was used to calculate the Percent Depth Dose (PDD) at different depths in an EDP. One field size for PDD has been examined. A 60Co photons source placed at 80 cm source to surface distance (SSD). The results of calculations were compared to TPS data obtained at National Institute of Oncology of Rabat.     

(CuAu3)n和(Cu3Au)n(n=1～15)团簇基态能量和结构的MOnte Carlo模拟研究

本文采用Monte Carlo方法和Gupta势函数对(CuAu3)n和(Cu3Au)n团簇的基态能量和结构进行了模拟研究,通过计算平均结合能、结合能的一阶差分和二阶差分,分析了团簇的稳定性.结果表明:(Cu3Au)n和(CuAu3)n团簇都为立体结构,都是以二十面体为基础形成的;(Cu3Au)n和(CuAu3)n团簇结构中金原子都有位于团簇表面的倾向;这两类团簇结构的区别在于:在(Cu3Au)n团簇中,铜一金原子混合程度高;而(CuAu3)n团簇中,形成金原子位于表层,铜原子位于中心的层状结构;且当n=3、5、7、9时,(Cu3Au)n和(CuAu3)n团簇在各自的序列中相对稳定性较邻近团簇高,特别是n=7的团簇,相对稳定性最高.     

Hausman检验统计量有效性的Monte Carlo模拟分析

在大样本下,Hausman检验统计量渐进地服从卡方分布,但在有限样本分析中却经常出现负值。对此,文献中并未形成一致看法。本文采用Monte Carlo模拟,从小样本偏误、内生性偏误、序列相关、非平行面板等角度,研究了Hausman检验统计量的小样本性质。结果表明,内生性问题(解释变量与个体效应相关)是导致Hausman统计量出现负值的主要原因。进一步分析表明,修正后的Hausman统计量,以及过度识别检验方法能够很好地克服上述缺陷,且具有很好的有限样本性质。     

基于Monte Carlo模拟比较K近邻和局部线性分位数回归

局部线性分位数回归是目前比较流行的非参数分位数回归,其潜在假定待估函数线性光滑.K近邻分位数回归也是非参数分位数回归的重要组成部分,其具有不需待估函数光滑和不同分位点的回归曲线不相交等优点.通过Monte Carlo模拟,比较了两者的估计,得到当待估函数的跳跃点或突变点比较多时,K近邻分位数回归的拟合效果优于局部线性回归.其中模拟的函数是Blocks、Bumps和HeaviSine的函数,它们分别代表跳跃性、波动性、斜率突变性的函数.