Fully parallel code for Monte Carlo simulation

We present a fully parallel version of Monte Carlo simulation of the Ising model using the Metropolis algorithm. In the 3-dimensional version the performance can be enhanced by a factor >20 in 16-bit word processors relative to other multispin codes. This factor could be further increased if implemented in 64-bit word computers.     

能量色散X荧光分析中元素间效应的蒙特卡罗模拟

针对能量色散X荧光分析(EDXRF)技术中元素间效应实验研究的难点问题,采用蒙特卡罗方法对基于Si（PIN）探测器的EDXRF系统建立模拟模型,并对模拟微束软X射线注量谱建立了高斯展宽算法。对Fe-Ni样品进行模拟计算,表明经该算法展宽后的注量谱与实测K系特征X射线谱吻合度较高,并得到了各元素特征X射线归一计数与元素含量关系曲线,结果表明该方法可自行校正EDXRF中元素间效应,获得准确的元素K系特征X射线谱理论强度。     

能量色散X荧光分析中元素间效应的蒙特卡罗模拟

针对能量色散X荧光分析(EDXRF)技术中元素间效应实验研究的难点问题,采用蒙特卡罗方法对基于Si（PIN）探测器的EDXRF系统建立模拟模型,并对模拟微束软X射线注量谱建立了高斯展宽算法。对Fe-Ni样品进行模拟计算,表明经该算法展宽后的注量谱与实测K系特征X射线谱吻合度较高,并得到了各元素特征X射线归一计数与元素含量关系曲线,结果表明该方法可自行校正EDXRF中元素间效应,获得准确的元素K系特征X射线谱理论强度。     

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

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

X射线能谱测量的蒙特卡罗成像模拟

针对高能强流电子束轰击高Z靶产生的X射线的能谱测量问题,采用蒙特卡罗方法进行成像模拟研究。高能X射线能谱通常由对X射线经过衰减体的直穿透射率曲线进行解谱获得。设计了带多准直孔的截锥体模型,在单次模拟成像中获得完整的衰减透射率曲线,有效避免了散射光子对透射率曲线以及X射线能谱重建的影响。成像面采用非均匀划分网格计数,将大部分探测计数点集中于各准直孔出口处,保证关注区域模拟成像的精细度,同时减少总体计算时间。     

Monte Carlo simulation code for confocal 3D micro‐beam X‐ray fluorescence analysis of stratified materials

Stratified materials are of great importance for many branches of modern industry, e.g. electronics or optics and for biomedical applications. Examination of chemical composition of individual layers and determination of their thickness helps to get information on their properties and function. A confocal 3D micro X‐ray fluorescence (3D µXRF) spectroscopy is an analytical method giving the possibility to investigate 3D distribution of chemical elements in a sample with spatial resolution in the micrometer regime in a non‐destructive way. Thin foils of Ti, Cu and Au, a bulk sample of Cu and a three‐layered sandwich sample, made of two thin Fe/Ni alloy foils, separated by polypropylene, were used as test samples. A Monte Carlo (MC) simulation code for the determination of elemental concentrations and thickness of individual layers in stratified materials with the use of confocal 3D µXRF spectroscopy was developed. The X‐ray intensity profiles versus the depth below surface, obtained from 3D µXRF experiments, MC simulation and an analytical approach were compared. Correlation coefficients between experimental versus simulated, and experimental versus analytical model X‐ray profiles were calculated. The correlation coefficients were comparable for both methods and exceeded 99%. The experimental X‐ray intensity profiles were deconvoluted with iterative MC simulation and by using analytical expression. The MC method produced slightly more accurate elemental concentrations and thickness of successive layers as compared to the results of the analytical approach. This MC code is a robust tool for simulation of scanning confocal 3D µXRF experiments on stratified materials and for quantitative interpretation of experimental results. Copyright © 2011 John Wiley & Sons, Ltd.     

A new Monte Carlo code for absorption simulation of laser-skin tissue interaction

In laser clinical applications, the process of photon absorption and thermal energy diffusion in the target tissue and its surrounding tissue during laser irradiation are crucial. Such information allows the selection of proper operating parameters such as laser power, and exposure time for optimal therapeutic. The Monte Carlo method is a useful tool for studying laser-tissue interaction and simulation of energy absorption in tissue during laser irradiation. We use the principles of this technique and write a new code with MATLAB 6.5, and then validate it against Monte Carlo multi layer (MCML) code. The new code is proved to be with good accuracy. It can be used to calculate the total power bsorbed in the region of interest. This can be combined for heat modelling with other computerized programs.     

ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision

A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations(ISSDE), is developed for solving stochastic differential equations(SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker–Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge–Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability.     

Risk analysis in investment appraisal based on the Monte Carlo simulation technique

This work has been prepared for the purpose of presenting the methodology and uses of the Monte Carlo simulation technique as applied in the evaluation of investment projects to analyze and assess risk. In the deterministic appraisal the basic decision rule for a project is simply to accept or reject the project depending on whether its net present value (NPV) is positive or negative, respectively. Similarly, when choosing among alternative (mutually exclusive) projects, the decision rule is to select the one with the highest NPV, provided that it is positive. Recognizing the fact that the key project variables (as: volume of sales, sales price, costs) are not certain, an appraisal report is usually supplemented to include sensitivity and scenario analysis tests. Both tests are static and rather arbitrary in their nature. During the simulation process, random scenarios are built up using input values for the project's key uncertain variables, which are selected from appropriate probability distributions. The results are collected and analyzed statistically so as to arrive at a probability distribution of the potential outcomes of the project and to estimate various measures of project risk. Received 25 September 2000     

快点火背光照相数值模拟

在激光惯性约束聚变实验研究中,诊断燃料层密度通常利用纳秒激光束辐照金属靶片产生特征X射线作为背光照明光源,对内爆球壳进行瞬态辐射照相,进而推断球壳内爆燃料层的密度。利用蒙特卡罗程序,对快点火X射线背光照相进行了研究,模拟了针孔成像中不同的针孔大小、障板厚度、背光能量对成像质量的影响,对于快点火实验中微米尺寸塑料小球,进行了诊断方案优化。模拟结果表明,综合考虑光通量及成像空间分辨,针孔半径为10 m时成像效果最佳;此时,采用2.7 keV的Mo背光源可以获得最佳的密度分辨。     

快点火背光照相数值模拟

在激光惯性约束聚变实验研究中,诊断燃料层密度通常利用纳秒激光束辐照金属靶片产生特征X射线作为背光照明光源,对内爆球壳进行瞬态辐射照相,进而推断球壳内爆燃料层的密度。利用蒙特卡罗程序,对快点火X射线背光照相进行了研究,模拟了针孔成像中不同的针孔大小、障板厚度、背光能量对成像质量的影响,对于快点火实验中微米尺寸塑料小球,进行了诊断方案优化。模拟结果表明,综合考虑光通量及成像空间分辨,针孔半径为10 m时成像效果最佳;此时,采用2.7 keV的Mo背光源可以获得最佳的密度分辨。     

Reverse electric field Monte Carlo simulation for vector radiative transfer in the atmosphere

In this paper, a reverse electric field Monte Carlo （REMC） method is proposed to study the vector radiation transfer in the atmosphere. The REMC is based on tracing the multiply scattered electric field to simulate the vector transmitted radiance. The reflected intensities with different total optical depth values are obtained, which accord well with the results in the previous research. Stokes vector and the degree of polarization are numerically investigated. The simulation result shows that when the solar zenith angle is determined, the zenith angle of detector has two points, of which the degree of polarization does not change with the ground albedo and the optical depth. The two points change regularly with the solar zenith angle. Moreover, our REMC method can be applied to the vector radiative transfer in the atmosphere-ocean system.     

Z箍缩中子编码诊断系统模拟平台搭建

开发了基于Geant4的Z箍缩中子编码成像系统模拟平台,实现聚变中子编码成像诊断系统各关键部件的完整模拟。获得了低中子产额（约1010量级）下,中子经编码孔编码后在闪烁体阵列中形成的发光分布图像。利用维纳滤波、Richardson-Lucy（RL）及遗传算法（GA）对低中子产额下获得的极低信噪比图像进行重建,并对信噪比、中子产额及重建效果进行了对比研究,结果表明：遗传算法对低信噪比中子编码图像的重建具有很强的鲁棒性;中子编码图像的信噪比与遗传算法重建结果的准确性呈正比。     

Z箍缩中子编码诊断系统模拟平台搭建

开发了基于Geant4的Z箍缩中子编码成像系统模拟平台,实现聚变中子编码成像诊断系统各关键部件的完整模拟。获得了低中子产额（约1010量级）下,中子经编码孔编码后在闪烁体阵列中形成的发光分布图像。利用维纳滤波、Richardson-Lucy（RL）及遗传算法（GA）对低中子产额下获得的极低信噪比图像进行重建,并对信噪比、中子产额及重建效果进行了对比研究,结果表明:遗传算法对低信噪比中子编码图像的重建具有很强的鲁棒性;中子编码图像的信噪比与遗传算法重建结果的准确性呈正比。     

基于蒙特卡罗方法的中子屏蔽材料设计

基于蒙特卡罗粒子输运程序MCNP,设计了一种强度高、密度低、具有优异中子屏蔽性能的新型玻璃纤维/B4C/环氧树脂复合材料,模拟计算了镅-铍（Am-Be）中子源产生中子对该材料的透射率;研究了该材料的中子屏蔽性能与传统屏蔽材料的差异以及不同B4C质量分数对该材料的屏蔽性能影响;根据模拟结果分析了该材料对不同能区中子（慢中子、中能中子、快中子）具有的不同屏蔽性能。研究发现：B4C质量分数为10%的该种新型玻璃纤维/B4C/环氧树脂复合材料的中子屏蔽性能,尤其是慢中子屏蔽性能较传统的含硼聚乙烯和Al-B4C合金材料更为优异;但随着B4C质量分数的增大,屏蔽性能提升不明显。结果验证了蒙特卡罗方法用于中子屏蔽材料优化设计的可行性。     

基于VERA基准题的cosRMC程序验证

为了验证反应堆物理软件和方法的计算能力,美国CASL (Consortium for Advanced Simulation of LWRs) 项目提出了VERA (Virtual Environment for Reactor Application) 堆芯物理基准题。该基准题以Watts Bar初始堆芯为模型,涵盖从二维单栅元到三维全堆芯的燃耗及换料的十个基准问题。针对VERA基准题模型,利用COSINE软件包中的反应堆蒙特卡罗分析程序cosRMC进行临界计算,得到了有效增殖因子、组件功率分布、控制棒微积分价值和反应性系数等结果。通过与基准题中提供的KENO结果对比,两种蒙特卡罗程序的计算结果吻合良好。这表明cosRMC程序具有从组件到堆芯的计算能力,其临界计算精度基本与KENO程序相当。     

FLUKA Monte Carlo simulation code as used for radiation studies in the Alice experiment

The main features of the FLUKA Monte Carlo code, which can deal with transport and interaction of electromagnetic and hadronic particles, are summarized. The physical models embedded in FLUKA are described. Application of the code for the radiation background calculations. Especially, the origin and composition of the intense radiation field to be expected in parts of the ALICE detector for the coming multi-TeV LHC collider is described. It is important to evaluate the risk of radiation damage in detectors and electronics equipment. The text was submitted by the author in English.     

A high-order Monte Carlo algorithm for the direct simulation of Boltzmann equation

A high-order algorithm of the direct simulation Monte Carlo (DSMC) method, H-DSMC, has been developed to simulate rarefied flow regimes. The mth order Taylor series expansion has been employed to obtain a more generalized form of the time discretization for the collision part of Boltzmann equation. In the purposed algorithm, the higher order collision terms are introduced as well as higher order terms in the time step of the probabilistic coefficients. These newly implemented higher order terms improve the accuracy and efficiency of the solution and enhance the convergence rate quite significantly. Comparison between results of the classic DSMC method and the H-DSMC method shows the promising performance of the introduced technique.