含有密集随机分布内核的椭球粒子光散射特性研究

给出了一种结合射线追踪和蒙特卡罗方法计算含核粒子光散射的方法,内核粒子可以为稀疏分布也可为浓密分布.粒子外边界的反射和折射由射线追踪方法计算,而粒子内部的多次散射过程由蒙特卡罗方法模拟;当内核粒子为浓密随机分布时,其单次散射特性由基于静态结构因子（static structure factor）的浓密介质光散射理论计算.最后讨论了含核椭球粒子模型的单次散射特性. 关键词： 射线追踪技术 蒙特卡洛方法 光散射 椭球粒子     

蒙特卡罗粒子输运程序JMCT研制

介绍了具有自主知识产权的蒙特卡罗(MC)粒子输运程序JMCT的初步研制成果。JMCT基于三维组合几何支撑软件框架JCOGIN,采用模块化,分成多层管理结构,可处理多群碰撞和连续能量碰撞,可进行粒子并行和区域分解并行两种并行方法,并具有良好的可扩展性和高速通信技术,同时配有可视化建模前处理软件。介绍了JMCT采用的多群物质碰撞机制,展示了程序模拟计算测试模型的结果,与MCNP程序计算结果一致。JMCT串行计算速度相比MCNP提高了约3 倍;在20 480个处理器核上模拟2109样本,并行效率可达70%。     

蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用

物理建模的精细化和三维模拟给蒙特卡罗粒子输运计算的规模成千上万倍的增加,甚至超过单核内存的最大规模,仅仅依靠传统的粒子并行蒙特卡罗计算无法实现对模型模拟,区域分解并行是可能的解决方法之一.然而区域分解带来了粒子在各区域间进行迁移,导致现有的伪随机数应用方式无法确保串行计算和并行计算的结果一致.针对这种现象,本文提出赋予粒子随机数属性和动态派生次级粒子随机数的技巧来确保区域分解并行计算的串并行结果一致.     

激光在不同类型气溶胶中传输特性研究

激光在大气中的传输衰减特性是激光工程应用中需要考虑的一个重要问题.本文针对常用的1.06 μm和10.6μm激光,基于Mie散射理论计算了气溶胶粒子的单次散射参量;对于激光在气溶胶中多次散射传输衰减,建立了蒙特卡罗模拟计算模型,利用Matlab语言编制了相应的计算程序,计算分析了两种波长的激光分别在沙尘性、水溶性、海洋性和煤烟性四种不同类型气溶胶中透过率与传播距离、能见度的关系,并将蒙特卡罗方法和单次散射的计算结果进行了比较.结果表明,当能见度较低、气溶胶粒子反照率较高时,单次散射计算存在很大的误差,用蒙特卡罗方法更能揭示多重散射现象;煤烟性气溶胶对1.06 μm激光的传输衰减影响最大,沙尘性气溶胶对10.6 μm激光的传输衰减影响最大.     

基于负氢离子源的全三维PIC/MCC模拟算法研究

在分析负氢离子源中等离子体物理机理基础下, 研究并优化粒子模拟算法, 设计高效的粒子存储方法. 研究并运用粒子碰撞蒙特卡罗方法, 考虑等离子体势以及带电粒子间库仑碰撞, 研制了全三维粒子模拟/蒙特卡罗算法(PIC/MCC). 采用磁荷模型, 运用时域有限差分方法计算多峰磁场, 并结合国外负氢离子源JT-60U, 考虑负氢离子源中主要反应, 对全三维PIC/MCC模拟算法模拟验证.     

组合几何蒙特卡罗粒子输运支撑软件框架JCOGIN介绍

介绍了JCOGIN支撑软件框架层次式、模块化的体系结构及其核心数据结构。这些使得JCOGIN软件框架能够支撑大规模复杂几何的蒙特卡罗（MC）粒子输运并行计算。同时介绍了在JCOGIN支撑软件框架下研发的JMCT程序,并进行了性能测试,在2万处理器核上模拟20亿粒子,并行效率达到70%。     

用蒙特卡罗方法计算光通过薄层随机分布粒子的散射

本文利用蒙特卡罗方法计算准直光束通过薄层随机分布粒子散射的透射和反射光强,并和输运理论的扩散近似结果做了比较。当散射接近各向同性时两者符合良好。当散射明显地成为各向异性时,蒙特卡罗方法的结果是合理的,而输运理论的扩散近似失效。     

闪光X射线照相蒙特卡罗程序的并行化

研究了闪光X射线辐射照像蒙特卡罗程序(FXRMC)在MPI平台下的并行计算实现,给出了实现过程中并行随机数的产生方法。对并行程序的测试结果表明:并行程序与串行程序结果一致,加速比比较理想,呈线性增长,并行效率在16个处理器上可达80%以上。算例的结果说明了并行化可有效地解决程序计算散射技术性能时的耗时问题,从而有效化解FXRMC耗时和大规模计算的困难,提高了FXRMC程序的计算规模和计算速度,达到了研究要求。 (Institute of Fluid Physics, CAEP, P. O. Box 919-105, Mianyang 621900, China)     

基于粒子模拟和并行遗传算法的高功率微波源优化设计

提出了基于粒子模拟和并行遗传算法的高功率微波源优化设计方法, 以全电磁粒子模拟软件UNIPIC模拟的高功率微波器件输出功率作为适应度函数, 采用浮点数编码的遗传算法对高功率微波源器件进行优化. 采用该算法, 对相对论返波管的布拉格反射器位置以及高度进行了浮点数编码,然后在巨型机上进行参数的全局优化, 获得了该返波管布拉格反射器的全局最优参数. 关键词： 并行遗传算法 相对论返波管 粒子模拟 高功率微波源     

三维光滑粒子流体动力学并行计算程序

介绍了基于消息传递并行程序设计平台研制的三维光滑粒子流体动力学并行程序CSPH3D.包括计算格式、并行方案、并行程序逻辑,以及加快邻域粒子搜索的处理方法.对三维微喷射和斜侵彻的计算表明:CSPH3D程序可以较好地计算这类问题.并且程序具有较高的并行效率.对于粒子总数为1527402的微喷射算例和粒子总数为1454225的斜侵彻算例,使用100个处理器时,并行效率可以达到80%.     

A boundary-based net-exchange Monte Carlo method for absorbing and scattering thick media

A boundary-based net-exchange Monte Carlo method was introduced (JQSRT 74 (2001) 563) that allows to bypass the difficulties encountered by standard Monte Carlo algorithms in the limit of optically thick absorption (and/or for quasi-isothermal configurations). With the present paper, this method is extended to scattering media. Developments are fully 3D, but illustrations are presented for plane parallel configuration. Compared to standard Monte Carlo algorithms, convergence qualities have been improved over a wide range of absorption and scattering optical thicknesses. The proposed algorithm still encounters a convergence difficulty in the case of optically thick, highly scattering media.     

纺织纤维光散射特性的模拟与实验

纺织纤维的光散射特性在纺织材料微观结构、光学性质以及无损检测中至关重要。利用纺织纤维各向异性的结构特点改进现有的角谱展开法,得到纺织纤维对倾斜入射平面波的散射以及表征其光散射特性的Mueller矩阵。采用蒙特卡罗算法模拟了平行棉纤维束对倾斜入射偏振光的多次散射,并对其计算结果进行了实验验证。结果表明理论计算和实验测量的散射光斑在形状特征、光强分布上均一致,从而验证了蒙特卡洛方法模拟纺织纤维光散射的正确性。结果同时也表明蒙特卡罗方法在纺织材料光传播特性的理论研究中具有重要的指导意义。     

Analysis of communication costs for domain decomposed Monte Carlo methods in nuclear reactor analysis

A domain decomposed Monte Carlo communication kernel is used to carry out performance tests to establish the feasibility of using Monte Carlo techniques for practical Light Water Reactor (LWR) core analyses. The results of the prototype code are interpreted in the context of simplified performance models which elucidate key scaling regimes of the parallel algorithm.     

三维中子-光子输运的蒙特卡罗程序MCMG

三维中子-光子输运蒙特卡罗程序MCMG发展了针对物质的碰撞机制,几何块、几何面动态可扩展, 随机数周期进一步扩大到261。可进行多群-连续截面耦合计算,多群散射展开到P5,并考虑了中子上散射,程序配备了通用和专用多群截面库。MCMG模拟取得了与MCNP程序和实验一致的结果,串行计算速度较MCNP快2~4倍,可进行上万处理器核的并行计算。     

用Monte Carlo方法模拟闪锌矿相ZnS电子的输运特性

用Monte Carlo方法模拟闪锌矿相(zinc blende)ZnS电子的输运特性.实验采用的是非抛物线模型计算电子的能带结构,模拟包含了声学声子散射,极性光学声子散射,压电散射,电离杂质散射,能谷间散射以及自散射等散射机理.通过模拟得到了ZnS材料的平均漂移速度、平均电子能量随电场强度变化的曲线图,以及总散射率随电子能量变化图,并将结果与文献报道的模拟结果[1]进行比较得出:本实验方法具有模型简单,计算速度快,获得结果比较准确的优点.     

光子与相对论麦克斯韦分布电子散射的能谱角度谱研究

光子与相对论麦克斯韦分布电子散射的描述及能谱角度谱计算非常复杂且费时.本文提出了一种光子与相对论麦克斯韦速度分布电子散射的蒙特卡罗(MC)模拟方法,该方法能够细致模拟高温等离子体中任意能量光子与任意温度电子的Compton和逆Compton散射问题.对于散射后光子的能谱和角度谱参数,可以根据电子温度抽样若干不同状态的电子,分别模拟其与光子发生散射,可以得到各次散射后的光子能量和偏转角度,取统计平均后的结果即可获得该光子与该温度电子散射的能谱和角度谱分布.根据该方法编写了光子与相对论电子散射MC模拟程序,开展了高温全电离等离子体中光子与相对论电子散射的能谱角度谱计算和分析,分析结果显示:热运动电子将展宽出射光子能谱,且低能光子与高温电子散射后的蓝移现象明显;出射光子的角度谱很复杂,其决定于入射光子能量、出射光子能量及电子温度.基于该方法计算并以数表形式给出的光子-相对论电子散射能谱角度谱数据,可以供辐射输运数值模拟程序使用.     

An efficient,robust, domain-decomposition algorithm for particle Monte Carlo

A previously described algorithm [T.A. Brunner, T.J. Urbatsch, T.M. Evans, N.A. Gentile, Comparison of four parallel algorithms for domain decomposed implicit Monte Carlo, Journal of Computational Physics 212 (2) (2006) 527–539] for doing domain decomposed particle Monte Carlo calculations in the context of thermal radiation transport has been improved. It has been extended to support cases where the number of particles in a time step are unknown at the beginning of the time step. This situation arises when various physical processes, such as neutron transport, can generate additional particles during the time step, or when particle splitting is used for variance reduction. Additionally, several race conditions that existed in the previous algorithm and could cause code hangs have been fixed. This new algorithm is believed to be robust against all race conditions. The parallel scalability of the new algorithm remains excellent.     

反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进

在反应堆pin-by-pin精细建模及蒙特卡罗模拟计算研究中, 由于不同栅元的功率密度差异较大, 导致蒙特卡罗方法临界计算的样本在不同栅元之间的分配不均衡, 由此引起栅元内的各种计数的统计误差差异较大. 为使大部分栅元内计数的统计误差降至一个合理的水平, 单纯增加总样本已不是一个高效的解决方法. 通过在特定临界计算迭代算法的基础上改进并实现均匀裂变源算法的思想, 对大亚湾压水堆pin-by-pin模型取得了具有较高效率的数值结果. 本工作为具有自主知识产权的蒙特卡罗粒子输运模拟软件JMCT最终达到反应堆pin-by-pin模型(包括一系列国际基准模型)的模拟性能要求提供了一个有效的工具.