X射线能谱处理中NaI(Tl)闪烁体边界效应矩阵的蒙特-卡罗计算

用 NaI(Tl)闪烁谱仪测量 X 射线机产生的 X 射线能谱是一个较为方便的方法。在此方法的应用中应当考虑到闪烁体边界效应的影响。本文报道了在解谱中表征这一影响的两个矩阵的蒙特-卡罗计算方法和计算结果,并对计算结果作了简要讨论。

Study on slope reliability method Response Surface Method based on Morgenstern-Price Method

The combination of Morgenstern-Price Method （MPM） with Response Surface Method （RSM） improves the eonventional RSM by replacing Finite-Element Method （FEM） with MPM, which is a type of reliability analysis method-RSM based on MPM, and makes the computing process simpler. The stability of the slope was re-evaluated by use of Monte Carlo method （MCM） and the result contrasted with the fore is the conclusion. The reliability analysis of the artificial slope in Donggang Power Station was carried on, and obtained the reliable charaeteristics, which provided basis for the project.

Reduced Expanding Load Method for Simulation-Based Structural System Reliability Analysis

The current situation and difficulties of the structural system reliability analysis are mentioned. Then on the basis of Monte Carlo method and computer simulation, a new analysis method reduced expanding load method (RELM) is presented, which can be used to solve structural reliability problems effectively and conveniently. In this method, the uncertainties of loads, structural material properties and dimensions can be fully considered. If the statistic parameters of stochastic variables are known, by using this method, the probability of failure can be estimated rather accurately. In contrast with traditional approaches,RELM method gives a much better understanding of structural failure frequency and its reliability index β is more meaningful.To illustrate this new idea, a specific example is given.

Asymptotic sum rate of MMSE receivers in multi-user distributed antenna systems

In this paper, the asymptotic sum rate of a multi-user distributed antenna system （DAS） is analyzed. To mitigate inter-user interference, minimum mean squared error （MMSE） receivers are utilized to cooperatively process received signals in the uplink. It shows that inter-user interference is efficiently mitigated and the uplink sum rate of a multi-user DAS is greatly improved by adopting MMSE receivers. For very large number of users and remote antennas, the asymptotic uplink sum rate of MMSE receivers is derived by using virtue of the random matrix theory, which can be The approximation is verified to be quite accurate by Monte Carlo simply calculated in an iterative way simulations.

Monte Carlo Hamiltonian： Inverse Potential

The Monte Carlo Hamiltonian method developed recently allows to investigate the ground state and low-lying excited states of a quantum system,using Monte Carlo(MC)algorithm with importance sampling.However,conventional MC algorithm has some difficulties when applied to inverse potentials.We propose to use effective potential and extrapolation method to solve the problem.We present examples from the hydrogen system.

Environmental Capacity of Petroleum Hydrocarbon Pollutants in Jiaozhou Bay, China: Modeling and Calculation

An environmental capacity model for the petroleum hydrocarbon pollutions （PHs） in Jiaozhou Bay is constructed based on field surveys, mesocosm, and parallel laboratory experiments. Simulated results of PHs seasonal successions in 2003 match the field surveys of Jiaozhou Bay resaonably well with a highest value in July. The Monte Carlo analysis confirms that the variation of PHs concentration significantly correlates with the river input. The water body in the bay is reasonably subjected to self-purification processes, such as volatilization to the atmosphere, biodegradation by microorganism, and transport to the Yellow Sea by water exchange. The environmental capacity of PHs in Jiaozhou Bay is 1500 tons per year IF the seawater quality criterion （Grade Ⅰ/Ⅱ, 0.05 mgLl） in the region is to be satisfied. The contribution to self-purification by volatilization, biodegradation, and transport to the Yellow Sea accounts for 48%, 28%, and 23%, respectively, which make these three processes the main ways of PHs purification in Jiaozhou Bay.

A New Alternative Method for Checking Total Stations

Today, just by activating a button on the display panel of a total station, one can measure horizontal and vertical angles or calculate distances and coordinates of points. For this reason, it has become necessary to improve the check of a total station in order to respond to its contemporary way of use. The aim of this research is to create a new method for checking total stations using only one check--field, in a different way of the today＇s standards. The basic principle, of the comparison of an instrument＇s measurements to those of another instrument with higher accuracy, is implemented into a mathematical model such as the adjustment of an indoor geodetic network. Thus, the check of the proper function of a total station is applied. An a-priori calculation Of^o of the network＇s adjustment can be calculated by modeling the uncertainties using statistical methods as the random number generators by means of the Monte Carlo method. The prospective ao can then be compared to the a-posteriori ao provided by the adjustment of the measurements which are carried out by using the total station under check.

Calculation model of military aircraft suryiyability to a missile

Since missiles are main threat against aircrafts in air war,a model is proposed for calculating the aircraft survivability to a missile.The hit characteristic of aircraft to a missile is analyzed,and then Monte Carlo method is applied to generate missile detonation location according to its distribution rule.In addition,based on the analysis of fragment trajectory and critical components,the intersection point of these two is determined.Then the kill probability of critical component to a fragment can be calculated,and the aircraft survivability to a missile is obtained accordingly.Finally,the feasibility of the proposed method is demonstrated.Simulation results show that this method captures the basic effects of missile detonation locations on aircraft survivability,which may provide an effective reference to aircraft survivability research.

Random Model and Simulation for Nucleation Stage of Emulsion Polymerization

The method of mathematical model and further computer simulation is an effective way to the theoretical study of emulsion polymerization and the scale-up of the reactors. In this work, Monte Carlo method has been used to simulate the nucleation of emulsion polymerization. The effects of emulsifier concentration [S] and initiator concentration [I] on various parameters such as the number of the particles (N p), the average diameter of the latex particles (D p), monomer conversion (x) and average radical number per particle (n) have been studied. The quantitative equations between [S], [I] and N p are in accord absolutely with the classical theory of Smith-Ewart.

Band Structure and Electron Transport of Bulk Based on an Ensemble Monte Carlo Calculation

The current demand growth of new components capable of operating at high power, high frequency, high temperatures and convergence towards miniaturization has lead to the development of new fields of nanotechnology based on II-VI semiconductor Interest in nanostructure：s based on II-VI semiconductor narrow gap containing mercury （such as super lattices HgTe/CdTe） was due to their advantages over alloys with cadmium telluride Mercury （MCT： HgCdTe）. The ternary alloy is a semiconductor band-gap direct, in that work the main interest is about the ternary compound. The results obtained are very satisfactory, they are compared with experimental results, and are in good agreement. These results are very promising and open new perspectives for the realization of solar cells and applications in the field of sensors.

Neutron Radiography System Collimator Design via Monte Carlo Calculation

Neutron radiography uses the unique interaction probabilities of neutrons to create images of materials. This imaging technique is non-destructive. MCNP Monte Carlo Code has been used to design an optimized neutron radiography system that utilizes 241Am-Be neutron source. Many different arrangements have been simulated to obtain a neutron flux with higher amplitude and more uniform distribution in the collimator outlet, next to image plane. In the final arrangement the specifications of neutron filter, Gamma-ray shield and beam collimator has been determined. Simulations has been carried out for a 5Ci 241Am-Be neutron source. In this case 43.8 n/cm2s thermal neutron flux has been achieved at a distance of 35 cm from neutron source.

Alternative Monte Carlo Approach for General Global Illumination

An alternative Monte Carlo strategy for the computation of global illumination problem was presented.The proposed approach provided a new and optimal way for solving Monte Carlo global illumination based on the zero variance importance sampling procedure. A new importance driven Monte Carlo global illumination algorithm in the framework of the new computing scheme was developed and implemented. Results, which were obtained by rendering test scenes, show that this new framework and the newly derived algorithm are effective and promising.

Simulating Particle Swarm Optimization Algorithm to Estimate Likelihood Function of ARMA（1, 1） Model

This paper present a simulation study of an evolutionary algorithms, Particle Swarm Optimization PSO algorithm to optimize likelihood function of ARMA（1, 1） model, where maximizing likelihood function is equivalent to maximizing its logarithm, so the objective function ＇obj.fun＇ is maximizing log-likelihood function. Monte Carlo method adapted for implementing and designing the experiments of this simulation. This study including a comparison among three versions of PSO algorithm “Constriction coefficient CCPSO, Inertia weight IWPSO, and Fully Informed FIPSO”, the experiments designed by setting different values of model parameters al, bs sample size n, moreover the parameters of PSO algorithms. MSE used as test statistic to measure the efficiency PSO to estimate model. The results show the ability of PSO to estimate ARMA＇ s parameters, and the minimum values of MSE getting for COPSO.

Research of Pipeline Pressure Based on the Reliability Theory

The former gas pipeline operating pressure transmission determined mostly according to the determined formula into the corresponding pipeline parameters are obtained, and the parameters of the actual pipeline due to many reasons, such as measurement error, production batch, etc., is not a fixed value. This paper on pipeline integrity established limit state equation, using Monte Carlo method to calculate the gas pipeline in different pressures of reliability, according to the API （American Petroleum Institute） 579 recommended target reliability to determine the operating pressure of the different regions, and design coefficient method, the calculated results are compared, results of calculation reliability is more reasonable and improve the delivery pressure of the area, the pipeline safety management provides the basis.

Effects of source opacity and expansion on two-pion Bose-Einstein correlations

We examine two-pion Bose-Einstein correlations for partially coherent particle-emitting sources within quantum statistical formal- ism, where the sources are treated as classical currents with chaotic and coherent components. The two-pion correlation functions of the partially coherent sources contain a phase which is sensitive to the asymmetry of the source emission function. We investigate the influence of source opacity and expansion in high energy heavy ion collisions on the phase by Monte Carlo calculations. We find that these two physical effects shift the phase from zero. The Gaussian-formula fit results to the simulated two-pion correlation functions indicate that the opaque and expansion effects lead to a smaller interferometry radius Rout and a larger 2 parameter.

The geometric factor of high energy protons detector on FY-3 satellite

Geometric factor is the key parameter for inversion of particle spectrum in space particle detection. Traditional geometric factor is obtained through the method of numerical calculation with the actual structure of the detector as the input condition. The degree of accuracy for data inversion is reduced since traditional geometric factor fails to take into account the physical process of interaction between the particle and substance as well as the influence of factors such as the particle interference between different energy channels on the measurement result. Here we propose an improved geometrical factor calculation method, the concept of which is to conduct actual structural modelling of the detector in the GEANT4 program, consider the process of interaction between the particle and substance, obtain the response function of the detector to particles of different energy channels through the method of Monte Carlo simulation, calculate the influence of contaminated particle on the geometrical factor, and finally get the geometrical factors for different energy channels of the detector. The imrpoved geometrical factor obtained through the method has carried out inversion for the data of high energy protons detector on China's FY-3 satellite, the energy spectrum after which is more in line with the power law distribution recognized by space physics. The comparison with the measured result of POES satellite indicates that the FY-3 satellite data are in good accordance with the satellite data, which shows the method may effectively improve the quality of data inversion.

Efficient implementation of the Barnes-Hut octree algorithm for Monte Carlo simulations of charged systems

Computer simulation with Monte Carlo is an important tool to investigate the function and equilibrium properties of many biological and soft matter materials solvable in solvents.The appropriate treatment of long-range electrostatic interaction is essential for these charged systems,but remains a challenging problem for large-scale simulations.We develop an efficient Barnes-Hut treecode algorithm for electrostatic evaluation in Monte Carlo simulations of Coulomb many-body systems.The algorithm is based on a divide-and-conquer strategy and fast update of the octree data structure in each trial move through a local adjustment procedure.We test the accuracy of the tree algorithm,and use it to perform computer simulations of electric double layer near a spherical interface.It is shown that the computational cost of the Monte Carlo method with treecode acceleration scales as log N in each move.For a typical system with ten thousand particles,by using the new algorithm,the speed has been improved by two orders of magnitude from the direct summation.

Progress in mesoscopic modeling of microstructure evolution in steels

The mesoscopic modeling developed rapidly in the past three decades is a promising tool for predicting and understanding the microstructure evolution at grain scale.In this paper,the recent development of mesoscopic modeling and its application to microstructure evolution in steels is reviewed.Firstly,some representative computational models are briefly introduced,e.g.,the phase field model,the cellular automaton model and the Monte Carlo model.Then,the emphasis is put on the application of mesoscopic modeling of the complex features of microstructure evolution,including solidification,solid-state phase transformation,recrystallization and grain growth.Finally,some issues in the present mesoscopic modeling and its perspective are discussed.