Monte Carlo method for evaluation of surface emission rate measurement uncertainty

The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.

Emergent topological ordered phase for the Ising-XY model revealed by cluster-updating Monte Carlo method

The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.

Application of the Monte Carlo Method in Simulative Generation Test Problems Simulative Generation Test Problems

An approach for generating test problems by a computer using the Monte Carlo method based upon user-given characterizations is described.A single point X~* is prespocified by the user to be a solution of the test problems.The approach is flex- ible enough to specify function values,gradients,Hesse,degeneracy degree and ill- conditioned degree at the point X~*.Other numerical features such as indefiniteness, convexity are also under user's control.

Application of Monte Carlo method in rudder control precision

After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting control moment error is carried out respectively by means of Monte Carlo method and the distribution of impact points of rudder control rocket is counted from all the successful subsamples. In the case of adding interference errors associated with rudder angle error and starting time error, the simulation analysis of impact point dispersion is done and its lateral and longitudinal correction abilities at different targeting angles are simulated to identify the effects of these factors on characteristics and control precision of the rudder control rocket, which provides the relevant reference for high-precision design of rudder control system.

Study of the Conductance Characteristic of Doped Polymer by Monte Carlo Method

The conduct mechanism of the doped polymer is considered. In an asymmetrysystem composed of high polymer and doping conductive matte, chain or congeries framework will beformed between the conductive particles to improve the conductance characteristic. In thisprocession, the conductive particles interact to each other. In this paper, we describe theconductance of the doped polymer by Monte Carlo method. The results accord with the experimentsquite well. It can be concluded that there is an evident change of doped polymer from nonconductorto metal.

Monte Carlo Method for the Uncertainty Evaluation of Spatial Straightness Error Based on New Generation Geometrical Product Specification

Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product speeifieation（GPS） requires the measurement uncertainty characterizing the reliability of the results should be given together when the measurement result is given. Nowadays most researches on straightness focus on error calculation and only several research projects evaluate the measurement uncertainty based on ＂The Guide to the Expression of Uncertainty in Measurement（GUM）＂. In order to compute spatial straightness error（SSE） accurately and rapidly and overcome the limitations of GUM, a quasi particle swarm optimization（QPSO） is proposed to solve the minimum zone SSE and Monte Carlo Method（MCM） is developed to estimate the measurement uncertainty. The mathematical model of minimum zone SSE is formulated. In QPSO quasi-random sequences are applied to the generation of the initial position and velocity of particles and their velocities are modified by the constriction factor approach. The flow of measurement uncertainty evaluation based on MCM is proposed, where the heart is repeatedly sampling from the probability density function（PDF） for every input quantity and evaluating the model in each case. The minimum zone SSE of a shaft measured on a Coordinate Measuring Machine（CMM） is calculated by QPSO and the measurement uncertainty is evaluated by MCM on the basis of analyzing the uncertainty contributors. The results show that the uncertainty directly influences the product judgment result. Therefore it is scientific and reasonable to consider the influence of the uncertainty in judging whether the parts are accepted or rejected, especially for those located in the uncertainty zone. The proposed method is especially suitable when the PDF of the measurand cannot adequately be approximated by a Gaussian distribution or a scaled and shifted t-distribution and the measurement model is non-linear.

Quantum Boltzmann equation solved by Monte Carlo method for nano-scale semiconductor devices simulation

A two-dimensional （2D） full band self-consistent ensemble Monte Carlo （MC） method for solving the quantum Boltzmann equation, including collision broadening and quantum potential corrections, is developed to extend the MC method to the study of nano-scale semiconductor devices with obvious quantum mechanical （QM） effects. The quantum effects both in real space and momentum space in nano-scale semiconductor devices can be simulated. The effective mobility in the inversion layer of n and p channel MOSFET is simulated and compared with experimental data to verify this method. With this method 50nm ultra thin body silicon on insulator MOSFET are simulated. Results indicate that this method can be used to simulate the 2D QM effects in semiconductor devices including tunnelling effect.

Forced propagation method for Monte Carlo fission source convergence acceleration in the RMC

In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.

Elbow precision machining technology by abrasive flow based on direct Monte Carlo method

The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.

Comparison of GUF and Monte Carlo methods to evaluate task-specific uncertainty in laser tracker measurement

Measurement uncertainty plays an important role in laser tracking measurement analyses. In the present work, the guides to the expression of uncertainty in measurement(GUM) uncertainty framework(GUF) and its supplement, the Monte Carlo method, were used to estimate the uncertainty of task-specific laser tracker measurements. First, the sources of error in laser tracker measurement were analyzed in detail, including instruments, measuring network fusion, measurement strategies, measurement process factors(such as the operator), measurement environment, and task-specific data processing. Second, the GUM and Monte Carlo methods and their application to laser tracker measurement were presented. Finally, a case study involving the uncertainty estimation of a cylindricity measurement process using the GUF and Monte Carlo methods was illustrated. The expanded uncertainty results(at 95% confidence levels) obtained with the Monte Carlo method are 0.069 mm(least-squares criterion) and 0.062 mm(minimum zone criterion), respectively, while with the GUM uncertainty framework, none but the result of least-squares criterion can be got, which is 0.071 mm. Thus, the GUM uncertainty framework slightly underestimates the overall uncertainty by 10%. The results demonstrate that the two methods have different characteristics in task-specific uncertainty evaluations of laser tracker measurements. The results indicate that the Monte Carlo method is a practical tool for applying the principle of propagation of distributions and does not depend on the assumptions and limitations required by the law of propagation of uncertainties(GUF). These features of the Monte Carlo method reduce the risk of an unreliable measurement of uncertainty estimation, particularly in cases of complicated measurement models, without the need to evaluate partial derivatives. In addition, the impact of sampling strategy and evaluation method on the uncertainty of the measurement results can also be taken into account with Monte Carlo method, which plays a guiding role in measurement planning.

Monte Carlo Analytic Hierarchy Process (MAHP) approach to selection of optimum mining method

One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.

IMPROVED STRATIFIED SAMPLING MONTE CARLO METHOD TO ANALYZE RELIABILITY OF STRUCTURAL SYSTEMS

A stratified sampling Monte Carlo method to analyze the reliability of structural systems is presented. Introducing a small exploratory simulation, this method overcomes the difficulties for getting the systematic sampling probability of all the strata. Several useful and efficient stratification methods are given and the strategies of stratification and simulation are studied. A general conclusion has been presented corresponding to actual engineering structures. The strict theoretical proof has been given,and it is especially effective to solve probabilistic integration. Statistic error of evaluating failure probability is reduced obviously. Especially in highly non-linear and nonreonvex problems, it is more accurate than other methods. Compared with other variance reduction techniques, this method can obtain a more obvious variance reduction and an increased sampling efficiency. Moreover, without strict limiting condition, it is convenient to use. This method is especially suitable to solve the reliability problem of structural systems with multiple failure modes and highly non-linear safety margin equations.

Application of local Monte Carlo method in neutronics calculation of EAST radial neutron camera

The Local Monte Carlo(LMC)method is used to solve the problems of deep penetration and long time in the neutronics calculation of the radial neutron camera(RNC)diagnostic system on the experimental advanced superconducting tokamak(EAST),and the radiation distribution of the RNC and the neutron flux at the detector positions of each channel are obtained.Compared with the results calculated by the global variance reduction method,it is shown that the LMC calculation is reliable within the reasonable error range.The calculation process of LMC is analyzed in detail,and the transport process of radiation particles is simulated in two steps.In the first step,an integrated neutronics model considering the complex window environment and a neutron source model based on EAST plasma shape are used to support the calculation.The particle information on the equivalent surface is analyzed to evaluate the rationality of settings of equivalent surface source and boundary.Based on the characteristic that only a local geometric model is needed in the second step,it is shown that the LMC is an advantageous calculation method for the nuclear shielding design of tokamak diagnostic systems.

Security Simulation of Continuous-Variable Quantum Key Distribution over Air-to-Water Channel Using Monte Carlo Method

Considering the ocean water＇s optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable （CV） quantum key distribution （QKD） based Mr-to-water channel. The effects of the absorp- tion and scattering on the transmittance of underwater quantum channel and the maximum secure transmission distance are studied. Considering the roughness of the sea surface, we simulate the performance bounds of CV QKD with different wind speeds using the Monte Carlo method. The results show that even if the secret key rate gradually reduces as the wind speed increases, the maximum transmission distance will not be affected obviously. Compared to the works regarding short-distance underwater optical communication, our research represents a significant step towards establishing secure communication between air platform and submarine vehicle.

Seismic fluid identification using a nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo method

Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.

SEMI-BLIND CHANNEL ESTIMATION OF MULTIPLE-INPUT/MULTIPLE-OUTPUT SYSTEMS BASED ON MARKOV CHAIN MONTE CARLO METHODS

This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and the transmitted signals. The deduced algorithms can work well under circumstances of low Signal-to-Noise Ratio (SNR). Simulation results are presented to demonstrate their effectiveness.

Sequential Monte Carlo Method Toward Online RUL Assessment with Applications

Online assessment of remaining useful life(RUL) of a system or device has been widely studied for performance reliability, production safety, system conditional maintenance, and decision in remanufacturing engineering. However,there is no consistency framework to solve the RUL recursive estimation for the complex degenerate systems/device.In this paper, state space model(SSM) with Bayesian online estimation expounded from Markov chain Monte Carlo(MCMC) to Sequential Monte Carlo(SMC) algorithm is presented in order to derive the optimal Bayesian estimation.In the context of nonlinear & non-Gaussian dynamic systems, SMC(also named particle filter, PF) is quite capable of performing filtering and RUL assessment recursively. The underlying deterioration of a system/device is seen as a stochastic process with continuous, nonreversible degrading. The state of the deterioration tendency is filtered and predicted with updating observations through the SMC procedure. The corresponding remaining useful life of the system/device is estimated based on the state degradation and a predefined threshold of the failure with two-sided criterion. The paper presents an application on a milling machine for cutter tool RUL assessment by applying the above proposed methodology. The example shows the promising results and the effectiveness of SSM and SMC online assessment of RUL.

COMPUTERIZED SIMULATION OF THERMODYNAMICS FOR MOLTEN LiF-KF SOLUTION BY MONTE CARLO METHOD

The structure and behaviour of LiF-KF solution,as a typical common-anion system,has been simulated by Monte Carlo method.The calculation of partial radial distribution function of ions,heat of mixing and potential energy distribution shows that the average distance be- tween Li^+ and F^- ions will significantly narrow after mixing of molten LiF and KF.This is very similar to the lean-on-one-side effect in molten LiF-KCl solution.The calculated heat of mixing is in fair agreement with the measured one.The dominant source of the energy of mixing may be that the decrease of the repulsion energy between cations,the decrease of the attraction energy between cations and anions,and the decrease of the repulsion energy be- tween anions.

Monte Carlo Method and Its Variance Reduction Techniques for On-stream Coal Slurry Parameter Determination

The special purpose Monte Carlo program McMesh was used to study neutron transport in coal slurries for on stream determination of the slurry parameters. McMesh uses the mesh weight window method as the major variance reduction technique with other methods such as exponential transforms and correlated sampling included as options. There was good agreement between the calculated results from McMesh and from MCNP, a general Monte Carlo program, but McMesh was more efficient and more convenient.

Speckle intensity images of target based on Monte Carlo method

Speckle intensity in the detector plane is deduced in the free-space optical system and imaging system based on Van Cittert-Zemike theorem. The speckle intensity images of plane target and conical target are obtained by using the Monte Carlo method and measured experimentally. The results show that when the range extent of target is smaller, the speckle size along the same direction become longer, and the speckle size increase with increasing incident light wavelengths. The speckle size increases and the speckle intensity images of target is closer to the actual object when the aperture scale augments. These findings are useful to access the target information by speckle in laser radar systems.