Approximation of photonuclear interaction cross-sections

Photonuclear interaction cross-sections from the GEANT4 database are approximated for all nuclei and all energies (from the hadron production threshold to about 40 TeV). The approximation methods in the giant-dipole resonance region, nucleon resonance region, and high-energy region are improved with respect to existing approximations. As an application of the approximation for photonuclear cross-sections, an improved method of calculating electronuclear cross-sections is developed. The interaction cross-section of virtual photons with nuclei at high Q² are approximated and a simple algorithm for describing the electronuclear reactions, including high-Q² scattering, is proposed. Received: 22 February 2002 / Accepted: 6 May 2002     

Tomographic reconstruction using heuristic Monte Carlo methods

A tomographic reconstruction method based on Monte Carlo random searching guided by the information contained in the projections of radiographed objects is presented. In order to solve the optimization problem, a multiscale algorithm is proposed to reduce computation. The reconstruction is performed in a coarse-to-fine multigrid scale that initializes each resolution level with the reconstruction of the previous coarser level, which substantially improves the performance. The method was applied to a real case reconstructing the internal structure of a small metallic object with internal components, showing excellent results.     

A Multiobjective Hybrid Metaheuristic Approach for GIS-based Spatial Zoning Model

This paper presents a multiobjective hybrid metaheuristic approach for an intelligent spatial zoning model in order to draw territory line for geographical or spatial zone for the purpose of space control. The model employs a Geographic Information System (GIS) and uses multiobjective combinatorial optimization techniques as its components. The proposed hybrid metaheuristic consists of the symbiosis between tabu search and scatter search method and it is used heuristically to generate non-dominated alternatives. The approach works with a set of current solution, which through manipulation of weights are optimized towards the non-dominated frontier while at the same time, seek to disperse over the frontier by a strategic oscillation concept. The general procedure and its algorithms are given as well as its implementation in the GIS environment. The computation has resulted in tremendous improvements in spatial zoning.     

激光诱导间质肿瘤热疗的数值模拟和实验研究

本文在考虑生物组织物性动态变化的情况下建立了激光诱导间质肿瘤热疗(LITT)的物理数学模型,采用MonteCarlo方法数值模拟了LITT中激光能量在生物组织内的传输过程,基于Pennes生物传热方程和Arrhenius方程数值求解了组织内的温度分布和热损伤体积的变化,分析了热物性及血液灌注率的动态变化对LITT过程的影响,并与相应的离体实验结果进行了对比。数值模拟结果表明,组织的热物性及血液灌注率的动态变化对于热损伤体积的变化具有重要的影响。因此在激光诱导间质肿瘤热疗的数值模拟中应该考虑热物性及血液灌注率的动态变化以期为临床治疗方案的制定提供更为准确的依据。     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

晶格失配对异质外延超薄膜生长中成核特性的影响

利用动力学蒙特卡罗方法模拟了异质外延超薄膜生长中的成核过程.研究了薄膜与衬底的晶格失配对超薄膜生长中成核密度、平均核尺寸、标度关系及生长模式的影响.结果发现产生压（张）应变的晶格负（正）失配使生长过程更早（迟）从成核区进入过渡区，失配越大，这一效应越明显.在相同的沉积条件下，负失配导致超薄膜形成较低的成核密度与较大的平均核尺寸，而正失配则相反.成核密度满足标度关系Ns≈(F/D)χ，随着失配度从-0.04增加到0.02，标度系数χ从0.37逐渐减小到0.33，对应超薄膜生长过程从包含二聚体扩散模式转变到无 关键词： 薄膜生长 成核 晶格失配 蒙特卡罗模拟     

闪烁光纤在γ射线辐照下部分特性的蒙特卡罗模拟

 闪烁光纤在射线成像方面的应用越来越广泛。为了进一步了解闪烁光纤在射线辐照下的基本特性，基于蒙特卡罗方法，利用计算机模拟分析了γ射线在闪烁光纤中的好事例率与光纤长度及射线能量的关系，能量沉积效率与光纤长度及射线能量的关系。此项工作对于闪烁光纤阵列在射线成像，剂量场测量等方面的研究很有价值。     

Studies on the Adsorption of Asymmetrical Triblock Copolymers by Scheutjens-Fleer Theory and Monte Carlo Simulation

The adsorption of asymmetrical triblock copolymers from a non-selective solvent on solid surface has been studied by using Scheutjens-Fleer mean-field theory and Monte Carlo simulation method on lattice model. The main aim of this paper is to provide detailed computer simulation data, taking A8-kB20Ak as a key example, to study the influence of the structure of copolymer on adsorption behavior and make a comparison between MC and SF results. The simulated results show that the size distribution of various configurations and density-profile are dependent on molecular structure and adsorption energy. The molecular structure will lead to diversity of adsorption behavior. This discrepancy between different structures would be enlarged for the surface coverage and adsorption amount with increasing of the adsorption energy. The surface coverage and the adsorption amount as well as the bound fraction will become larger as symmetry of the molecular structure becomes gradually worse. The adsorption layer becomes thicker with increasing of symmetry of the molecule when adsorption energy is smaller but it becomes thinner when adsorption energy is higher. It is shown that SF theory can reproduce the adsorption behavior of asymmetrical triblock copolymers. However, systematic discrepancy between the theory and simulation still exists.The approximations inherited in the mean-filed theory such as random mixing and the allowance of direct back folding may be responsible for those deviations.     

Vibrational control of crystal growth from liquid phase

Modeling and numerical simulations of the convective flows induced by the vibration of the monocrystal during crystal growth have been performed for two configurations simulating the Cz and FZ methods. This permitted to emphasize the role of different vibrational mechanisms in the formation of the average flows. It is shown that an appropriate combination of these mechanisms can be used to counteract the usual convective flows (buoyancy- and/or thermocapillary-driven) inherent to crystal growth processes from the liquid phase. While vibrational convection is rather complex due to these identified mechanisms, the new modeling used in the present paper opens up very promising perspectives to efficiently control heat and mass transfer during real industrial applications of crystal growth from the liquid phase.     

一种改进的计算探测器校正因子的相关抽样方法

对于小尺寸探测器处于大块介质的情形, 在探测器的校正因子的Monte Carlo模拟中, 存在两个难题: 一是由于探测器尺寸很小粒子难以到达探测器并发生碰撞; 二是两个随机变量比值难以达到要求精度. 本文使用经过改进的粒子碰撞自动重要抽样方法, 再结合相关抽样方法, 解决了这两个难题, 并在MCNP-4C程序平台上加以实现. 除了粒子碰撞自动重要抽样以外, 还选用了其他3种方法: 直接模拟、区域分裂、强迫碰撞+Dxtran球分别与相关抽样方法结合, 对一个简化的探测器校正因子计算模型进行了计算. 实际计算结果表明, 相关抽样方法无论与哪种方法结合, 都起到了提高相关量计算效率的作用; 而它与粒子碰撞自动重要抽样结合, 比其他方法具有明显的优越性.