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 predictio...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.展开更多
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...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.展开更多
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 supercondu...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.展开更多
Haze concentration prediction,especially PM2.5,has always been a significant focus of air quality research,which is necessary to start a deep study.Aimed at predicting the monthly average concentration of PM2.5 in Bei...Haze concentration prediction,especially PM2.5,has always been a significant focus of air quality research,which is necessary to start a deep study.Aimed at predicting the monthly average concentration of PM2.5 in Beijing,a novel method based on Monte Carlo model is conducted.In order to fully exploit the value of PM2.5 data,we take logarithmic processing of the original PM2.5 data and propose two different scales of the daily concentration and the daily chain development speed of PM2.5 respectively.The results show that these data are both approximately normal distribution.On the basis of the results,a Monte Carlo method can be applied to establish a probability model of normal distribution based on two different variables and random sampling numbers can also be generated by computer.Through a large number of simulation experiments,the average monthly concentration of PM2.5 in Beijing and the general trend of PM2.5 can be obtained.By comparing the errors between the real data and the predicted data,the Monte Carlo method is reliable in predicting the PM2.5 monthly mean concentration in the area.This study also provides a feasible method that may be applied in other studies to predict other pollutants with large scale time series data.展开更多
Calculation results of the Monte Carlo method of the average energy of the electrostatic interaction between the quarks are presented to the neutron and proton. The proposed model of the distribution of quarks in prot...Calculation results of the Monte Carlo method of the average energy of the electrostatic interaction between the quarks are presented to the neutron and proton. The proposed model of the distribution of quarks in protons and neutrons is possible to assess the area which included a strong (gluon) interaction. Given the fact that the probability of finding a quark in the field with strong interaction is less than one, there is a good agreement between the experimental and calculated values of the mass difference between the neutron and the proton.展开更多
Solution of the system stochastic differential equations in multi dimensional case using Monte Carlo method had many useful features in compare with the other computational methods. One of them is the solution of boun...Solution of the system stochastic differential equations in multi dimensional case using Monte Carlo method had many useful features in compare with the other computational methods. One of them is the solution of boundary value problems to be found at just one point, if required (with associated saving in computation), whereas deterministic methods necessarily find the solution at large number of points simultaneously. This property can be particularly useful in problems such option pricing, where the value of an option is required only at the time of striking, and for the state of the market at that time. In this work we consider a European multi-asset options which mathematically described by the system of stochastic differential equations. We will apply Monte Carlo method for the solution of that system which is the price of Multi-asset rainbow options.展开更多
A random simulation method was used for treatment of systems of Volterra integral equations of the second kind. Firstly, a linear algebra system was obtained by discretization using quadrature formula. Secondly, this ...A random simulation method was used for treatment of systems of Volterra integral equations of the second kind. Firstly, a linear algebra system was obtained by discretization using quadrature formula. Secondly, this algebra system was solved by using relaxed Monte Carlo method with importance sampling and numerical approximation solutions of the integral equations system were achieved. It is theoretically proved that the validity of relaxed Monte Carlo method is based on importance sampling to solve the integral equations system. Finally, some numerical examples from literatures are given to show the efficiency of the method.展开更多
Delayed neutron loss is an important parameter in the safety analysis of molten salt reactors. In this study,to obtain the effective delayed neutron fraction under flow condition, a delayed neutron precursor transport...Delayed neutron loss is an important parameter in the safety analysis of molten salt reactors. In this study,to obtain the effective delayed neutron fraction under flow condition, a delayed neutron precursor transport was implemented in the Monte Carlo code MCNP. The moltensalt reactor experiment(MSRE) model was used to analyze the reliability of this method. The obtained flow losses of reactivity for 235 U and 233 U fuels in the MSRE are223 pcm and 100.8 pcm, respectively, which are in good agreement with the experimental values(212 pcm and100.5 pcm, respectively). Then, six groups of effective delayed neutron fractions in a small molten salt reactor were calculated under different mass flow rates. The flow loss of reactivity at full power operation is approximately105.6 pcm, which is significantly lower than that of the MSRE due to the longer residence time inside the active core. The sensitivity of the reactivity loss to other factors,such as the residence time inside or outside the core and flow distribution, was evaluated as well. As a conclusion,the sensitivity of the reactivity loss to the residence time inside the core is greater than to other parameters.展开更多
Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product specification(GPS) requires the measurement uncertainty characterizing the reliability of the results...Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product specification(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.展开更多
To improve the precisions of pose error analysis for 6-dof parallel kinematic mechanism( PKM)during assembly quality control,a Sobol sequence based on Quasi Monte Carlo( QMC) method is introduced and implemented in po...To improve the precisions of pose error analysis for 6-dof parallel kinematic mechanism( PKM)during assembly quality control,a Sobol sequence based on Quasi Monte Carlo( QMC) method is introduced and implemented in pose accuracy analysis for the PKM in this paper. The Sobol sequence based on Quasi Monte Carlo with the regularity and uniformity of samples in high dimensions,can prevail traditional Monte Carlo method with up to 98. 59% and 98. 25% enhancement for computational precision of pose error statistics.Then a PKM tolerance design system integrating this method is developed and with it pose error distributions of the PKM within a prescribed workspace are finally obtained and analyzed.展开更多
The performance in vibration environment of switching apparatus containing mechanical contact is an important element when judging the apparatus’s reliability. A piecewise linear two-degrees-of-freedom mathematical m...The performance in vibration environment of switching apparatus containing mechanical contact is an important element when judging the apparatus’s reliability. A piecewise linear two-degrees-of-freedom mathematical model considering contact loss was built in this work, and the vibration performance of the model under random external Gaussian white noise excitation was investigated by using Monte Carlo simulation in Matlab/Simulink. Simulation showed that the spectral content and statistical characters of the contact force coincided strongly with reality. The random vibration character of the contact system was solved using time (numerical) domain simulation in this paper. Conclusions reached here are of great importance for reliability design of switching apparatus.展开更多
Based on the assumption of the plain-strain problem, various optimization or random search methods have been developed for locating the critical slip surfaces in slope-stability analysis, but none of such methods is a...Based on the assumption of the plain-strain problem, various optimization or random search methods have been developed for locating the critical slip surfaces in slope-stability analysis, but none of such methods is applicable to the 3D case. In this paper, a simple Monte Carlo random simulation method is proposed to identify the 3D critical slip surface. Assuming the initial slip to be the lower part of a slip ellipsoid, the 3D critical slip surface is located by means of a minimized 3D safety factor. A column-based 3D slope stability analysis model is used to calculate this factor. In this study, some practical cases of known minimum safety factors and critical slip surfaces in 2D analysis are extended to 3D slope problems to locate the critical slip surfaces. Compared with the 2D result, the resulting 3D critical slip surface has no apparent difference in terms of only cross section, but the associated 3D safety factor is definitely higher.展开更多
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 M...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 efficiency and stability of inversion are both taken into consideration 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.展开更多
An improved Monte Carlo method was used to simulate the motion of electrons in c-C_4F_8 and SF_6 gas mixtures for pulsed townsend discharge. The electron swarm parameters such as effective ionization coefficient, [AK...An improved Monte Carlo method was used to simulate the motion of electrons in c-C_4F_8 and SF_6 gas mixtures for pulsed townsend discharge. The electron swarm parameters such as effective ionization coefficient, [AKα-] and drift velocity over the E/N range from 280~700 Td(1Td=10-21 V·m2) were calculated by employing a set of cross sections available in literature. From the variation cure of [AKα-] with SF_6 partial pressure p, the limiting field (E/N)_lim of gas mixture at different gas content was determined. It is found that the limiting field of c-C_4F_8 and SF_6 gas mixture is higher than that of pure SF_6 at any SF_6 mixture ratio. Simulation results show excellent agreement with experiment data available in previous literature.展开更多
A new reliability allocation model has been built for engine system,which is a repairable system,and consists of a large number of mechanical components.The cost and reliability are taken as objective function and con...A new reliability allocation model has been built for engine system,which is a repairable system,and consists of a large number of mechanical components.The cost and reliability are taken as objective function and constraint condition respectively.The parameters of components lifetime distribution are given as decision variables,and the component lifetimes are assumed to follow that Weibull distribution.The allocation is separated into two steps to reduce calculated amount of one allocation.Genetic algorithm and Monte Carlo method are applied to solve distribution parameters and system cost separately.展开更多
Three-dimensional normal grain growth was appropriately simulated using a Potts model Monte Carlo algorithm. The quasi-stationary grain size distribution obtained from simulation agreed well with the experimental resu...Three-dimensional normal grain growth was appropriately simulated using a Potts model Monte Carlo algorithm. The quasi-stationary grain size distribution obtained from simulation agreed well with the experimental result of pure iron. The Weibull function with a parameter β=2.77 and the Yu-Liu function with a parameter ν =2.71 fit the quasi-stationary grain size distribution well. The grain volume distribution is a function that decreased exponentially with increasing grain volume. The distribution of boundary area of grains has a peak at S/展开更多
Interface roughness strongly influences the performance of germanium metal-organic-semiconductor field effect transistors(MOSFETs).In this paper,a 2D full-band Monte Carlo simulator is used to study the impact of inte...Interface roughness strongly influences the performance of germanium metal-organic-semiconductor field effect transistors(MOSFETs).In this paper,a 2D full-band Monte Carlo simulator is used to study the impact of interface roughness scattering on electron and hole transport properties in long-and short-channel Ge MOSFETs inversion layers.The carrier effective mobility in the channel of Ge MOSFETs and the in non-equilibrium transport properties are investigated.Results show that both electron and hole mobility are strongly influenced by interface roughness scattering.The output curves for 50 nm channel-length double gate n and p Ge MOSFET show that the drive currents of n-and p-Ge MOSFETs have significant improvement compared with that of Si n-and p-MOSFETs with smooth interface between channel and gate dielectric.The 82% and 96% drive current enhancement are obtained for the n-and p-MOSFETs with the completely smooth interface.However,the enhancement decreases sharply with the increase of interface roughness.With the very rough interface,the drive currents of Ge MOSFETs are even less than that of Si MOSFETs.Moreover,the significant velocity overshoot also has been found in Ge MOSFETs.展开更多
CH3SiCl3 (MTS)-H2-Ar system has been applied to prepare SiC film with chemical vapor deposition (CVD) method in this paper. For three facets of SiC film, some significant influence on growth rate, surface roughness, t...CH3SiCl3 (MTS)-H2-Ar system has been applied to prepare SiC film with chemical vapor deposition (CVD) method in this paper. For three facets of SiC film, some significant influence on growth rate, surface roughness, thickness and relative density brought by MTS consistency has been mainly discussed with kinetic monte carlo (KMC) method. The simulation results show that there is a certain scale for mol ratio of H2 to MTS (H2/MTS) with different deposition temperature. When MTS consistency increases, growth rate and surface roughness of three facets all increase, which manifests approximate linearity relationship. Thickness of three facets also increases while increasing trend of three facets thickness is different obviously. Although relative density of three facets all increases, increasing trend shows a little difference with MTS consistency increasing.展开更多
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 supplemen...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.展开更多
基金Project supported by the Hefei National Research Center for Physical Sciences at the Microscale (Grant No.KF2021002)the Natural Science Foundation of Shanxi Province,China (Grant Nos.202303021221029 and 202103021224051)+2 种基金the National Natural Science Foundation of China (Grant Nos.11975024,12047503,and 12275263)the Anhui Provincial Supporting Program for Excellent Young Talents in Colleges and Universities (Grant No.gxyq ZD2019023)the National Key Research and Development Program of China (Grant No.2018YFA0306501)。
文摘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.
基金Supported by Basic Research and Development Plan of China (973 Program,Grant Nos.2011CB013401,2011CB013402)Special Fundamental Research Funds for Central Universities of China(Grant No.DUT14QY21)
文摘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.
基金support and help in this research.This work was supported by Users with Excellence Program of Hefei Science Center CAS(No.2020HSC-UE012)Comprehensive Research Facility for Fusion Technology Program of China(No.2018-000052-73-01-001228)National Natural Science Foundation of China(No.11605241)。
文摘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.
文摘Haze concentration prediction,especially PM2.5,has always been a significant focus of air quality research,which is necessary to start a deep study.Aimed at predicting the monthly average concentration of PM2.5 in Beijing,a novel method based on Monte Carlo model is conducted.In order to fully exploit the value of PM2.5 data,we take logarithmic processing of the original PM2.5 data and propose two different scales of the daily concentration and the daily chain development speed of PM2.5 respectively.The results show that these data are both approximately normal distribution.On the basis of the results,a Monte Carlo method can be applied to establish a probability model of normal distribution based on two different variables and random sampling numbers can also be generated by computer.Through a large number of simulation experiments,the average monthly concentration of PM2.5 in Beijing and the general trend of PM2.5 can be obtained.By comparing the errors between the real data and the predicted data,the Monte Carlo method is reliable in predicting the PM2.5 monthly mean concentration in the area.This study also provides a feasible method that may be applied in other studies to predict other pollutants with large scale time series data.
文摘Calculation results of the Monte Carlo method of the average energy of the electrostatic interaction between the quarks are presented to the neutron and proton. The proposed model of the distribution of quarks in protons and neutrons is possible to assess the area which included a strong (gluon) interaction. Given the fact that the probability of finding a quark in the field with strong interaction is less than one, there is a good agreement between the experimental and calculated values of the mass difference between the neutron and the proton.
文摘Solution of the system stochastic differential equations in multi dimensional case using Monte Carlo method had many useful features in compare with the other computational methods. One of them is the solution of boundary value problems to be found at just one point, if required (with associated saving in computation), whereas deterministic methods necessarily find the solution at large number of points simultaneously. This property can be particularly useful in problems such option pricing, where the value of an option is required only at the time of striking, and for the state of the market at that time. In this work we consider a European multi-asset options which mathematically described by the system of stochastic differential equations. We will apply Monte Carlo method for the solution of that system which is the price of Multi-asset rainbow options.
文摘A random simulation method was used for treatment of systems of Volterra integral equations of the second kind. Firstly, a linear algebra system was obtained by discretization using quadrature formula. Secondly, this algebra system was solved by using relaxed Monte Carlo method with importance sampling and numerical approximation solutions of the integral equations system were achieved. It is theoretically proved that the validity of relaxed Monte Carlo method is based on importance sampling to solve the integral equations system. Finally, some numerical examples from literatures are given to show the efficiency of the method.
基金supported by the Chinese TMSR Strategic Pioneer Science and Technology Project(No.XDA02010000)the Frontier Science Key Program of the Chinese Academy of Sciences(No.QYZDY-SSW-JSC016)
文摘Delayed neutron loss is an important parameter in the safety analysis of molten salt reactors. In this study,to obtain the effective delayed neutron fraction under flow condition, a delayed neutron precursor transport was implemented in the Monte Carlo code MCNP. The moltensalt reactor experiment(MSRE) model was used to analyze the reliability of this method. The obtained flow losses of reactivity for 235 U and 233 U fuels in the MSRE are223 pcm and 100.8 pcm, respectively, which are in good agreement with the experimental values(212 pcm and100.5 pcm, respectively). Then, six groups of effective delayed neutron fractions in a small molten salt reactor were calculated under different mass flow rates. The flow loss of reactivity at full power operation is approximately105.6 pcm, which is significantly lower than that of the MSRE due to the longer residence time inside the active core. The sensitivity of the reactivity loss to other factors,such as the residence time inside or outside the core and flow distribution, was evaluated as well. As a conclusion,the sensitivity of the reactivity loss to the residence time inside the core is greater than to other parameters.
基金Project supported by the Special Foundation for State Major Basic Research Program of China (Grant No G2000035602) and the National Natural Science Foundation of China (Grant No 90307006).
基金supported by National Natural Science Foundation of China (Grant No. 51075198)Jiangsu Provincial Natural Science Foundation of China (Grant No. BK2010479)+2 种基金Innovation Research of Nanjing Institute of Technology, China (Grant No. CKJ20100008)Jiangsu Provincial Foundation of 333 Talents Engineering of ChinaJiangsu Provincial Foundation of Six Talented Peak of China
文摘Straightness error is an important parameter in measuring high-precision shafts. New generation geometrical product specification(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.
基金Sponsored by the National Defense Basic Scientific Research Program(Grant No.A0320110019)the Shanghai Science and Technology Innovation Action Plan(Grant No.11DZ1120800)
文摘To improve the precisions of pose error analysis for 6-dof parallel kinematic mechanism( PKM)during assembly quality control,a Sobol sequence based on Quasi Monte Carlo( QMC) method is introduced and implemented in pose accuracy analysis for the PKM in this paper. The Sobol sequence based on Quasi Monte Carlo with the regularity and uniformity of samples in high dimensions,can prevail traditional Monte Carlo method with up to 98. 59% and 98. 25% enhancement for computational precision of pose error statistics.Then a PKM tolerance design system integrating this method is developed and with it pose error distributions of the PKM within a prescribed workspace are finally obtained and analyzed.
基金Project (No. FEBQ24409102) supported by the Space Technology Innovation Fund, China
文摘The performance in vibration environment of switching apparatus containing mechanical contact is an important element when judging the apparatus’s reliability. A piecewise linear two-degrees-of-freedom mathematical model considering contact loss was built in this work, and the vibration performance of the model under random external Gaussian white noise excitation was investigated by using Monte Carlo simulation in Matlab/Simulink. Simulation showed that the spectral content and statistical characters of the contact force coincided strongly with reality. The random vibration character of the contact system was solved using time (numerical) domain simulation in this paper. Conclusions reached here are of great importance for reliability design of switching apparatus.
文摘Based on the assumption of the plain-strain problem, various optimization or random search methods have been developed for locating the critical slip surfaces in slope-stability analysis, but none of such methods is applicable to the 3D case. In this paper, a simple Monte Carlo random simulation method is proposed to identify the 3D critical slip surface. Assuming the initial slip to be the lower part of a slip ellipsoid, the 3D critical slip surface is located by means of a minimized 3D safety factor. A column-based 3D slope stability analysis model is used to calculate this factor. In this study, some practical cases of known minimum safety factors and critical slip surfaces in 2D analysis are extended to 3D slope problems to locate the critical slip surfaces. Compared with the 2D result, the resulting 3D critical slip surface has no apparent difference in terms of only cross section, but the associated 3D safety factor is definitely higher.
基金the sponsorship of the National Basic Research Program of China (973 Program,2013CB228604,2014CB239201)the National Oil and Gas Major Projects of China (2011ZX05014-001-010HZ,2011ZX05014-001-006-XY570) for their funding of this research
文摘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 efficiency and stability of inversion are both taken into consideration 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.
文摘An improved Monte Carlo method was used to simulate the motion of electrons in c-C_4F_8 and SF_6 gas mixtures for pulsed townsend discharge. The electron swarm parameters such as effective ionization coefficient, [AKα-] and drift velocity over the E/N range from 280~700 Td(1Td=10-21 V·m2) were calculated by employing a set of cross sections available in literature. From the variation cure of [AKα-] with SF_6 partial pressure p, the limiting field (E/N)_lim of gas mixture at different gas content was determined. It is found that the limiting field of c-C_4F_8 and SF_6 gas mixture is higher than that of pure SF_6 at any SF_6 mixture ratio. Simulation results show excellent agreement with experiment data available in previous literature.
文摘A new reliability allocation model has been built for engine system,which is a repairable system,and consists of a large number of mechanical components.The cost and reliability are taken as objective function and constraint condition respectively.The parameters of components lifetime distribution are given as decision variables,and the component lifetimes are assumed to follow that Weibull distribution.The allocation is separated into two steps to reduce calculated amount of one allocation.Genetic algorithm and Monte Carlo method are applied to solve distribution parameters and system cost separately.
基金supported by the National Natural Science Foundation of China (No.50671010)
文摘Three-dimensional normal grain growth was appropriately simulated using a Potts model Monte Carlo algorithm. The quasi-stationary grain size distribution obtained from simulation agreed well with the experimental result of pure iron. The Weibull function with a parameter β=2.77 and the Yu-Liu function with a parameter ν =2.71 fit the quasi-stationary grain size distribution well. The grain volume distribution is a function that decreased exponentially with increasing grain volume. The distribution of boundary area of grains has a peak at S/
基金Project supported by the National Natural Science Foundation of China (Grant No. 60606013)the National Basic Research Program of China (Grant No. 2006CB302705)
文摘Interface roughness strongly influences the performance of germanium metal-organic-semiconductor field effect transistors(MOSFETs).In this paper,a 2D full-band Monte Carlo simulator is used to study the impact of interface roughness scattering on electron and hole transport properties in long-and short-channel Ge MOSFETs inversion layers.The carrier effective mobility in the channel of Ge MOSFETs and the in non-equilibrium transport properties are investigated.Results show that both electron and hole mobility are strongly influenced by interface roughness scattering.The output curves for 50 nm channel-length double gate n and p Ge MOSFET show that the drive currents of n-and p-Ge MOSFETs have significant improvement compared with that of Si n-and p-MOSFETs with smooth interface between channel and gate dielectric.The 82% and 96% drive current enhancement are obtained for the n-and p-MOSFETs with the completely smooth interface.However,the enhancement decreases sharply with the increase of interface roughness.With the very rough interface,the drive currents of Ge MOSFETs are even less than that of Si MOSFETs.Moreover,the significant velocity overshoot also has been found in Ge MOSFETs.
基金National Natural Science Foundation of China(No.50871086)the Research Fund of State Key Laboratory of Solidification Processing (NPU), China(No.KP200906)the President Research Foundation of Xi'an Technological University(No.XAGDXJJ1008)
文摘CH3SiCl3 (MTS)-H2-Ar system has been applied to prepare SiC film with chemical vapor deposition (CVD) method in this paper. For three facets of SiC film, some significant influence on growth rate, surface roughness, thickness and relative density brought by MTS consistency has been mainly discussed with kinetic monte carlo (KMC) method. The simulation results show that there is a certain scale for mol ratio of H2 to MTS (H2/MTS) with different deposition temperature. When MTS consistency increases, growth rate and surface roughness of three facets all increase, which manifests approximate linearity relationship. Thickness of three facets also increases while increasing trend of three facets thickness is different obviously. Although relative density of three facets all increases, increasing trend shows a little difference with MTS consistency increasing.
基金Project(51318010402)supported by General Armament Department Pre-Research Program of China
文摘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.