This paper proposes an active learning accelerated Monte-Carlo simulation method based on the modified K-nearest neighbors algorithm.The core idea of the proposed method is to judge whether or not the output of a rand...This paper proposes an active learning accelerated Monte-Carlo simulation method based on the modified K-nearest neighbors algorithm.The core idea of the proposed method is to judge whether or not the output of a random input point can be postulated through a classifier implemented through the modified K-nearest neighbors algorithm.Compared to other active learning methods resorting to experimental designs,the proposed method is characterized by employing Monte-Carlo simulation for sampling inputs and saving a large portion of the actual evaluations of outputs through an accurate classification,which is applicable for most structural reliability estimation problems.Moreover,the validity,efficiency,and accuracy of the proposed method are demonstrated numerically.In addition,the optimal value of K that maximizes the computational efficiency is studied.Finally,the proposed method is applied to the reliability estimation of the carbon fiber reinforced silicon carbide composite specimens subjected to random displacements,which further validates its practicability.展开更多
As one of the most serious natural disasters,many typhoons affect southeastern China every year.Taking Shenzhen,a coastal city in southeast China as an example,we employed a Monte-Carlo simulation to generate a large ...As one of the most serious natural disasters,many typhoons affect southeastern China every year.Taking Shenzhen,a coastal city in southeast China as an example,we employed a Monte-Carlo simulation to generate a large number of virtual typhoons for wind hazard analysis.By analyzing 67-year historical typhoons data from 1949 to 2015 using the Best Track Dataset for Tropical Cyclones over the Western North Pacific recorded by the Shanghai Typhoon Institute,China Meteorological Administration(CMASTI),typhoon characteristic parameters were extracted and optimal statistical distributions established for the parameters in relation to Shenzhen.We employed the Monte-Carlo method to sample each distribution to generate the characteristic parameters of virtual typhoons.In addition,the Yah Meng(YM)wind field model was introduced,and the sensitivity of the YM model to several parameters discussed.Using the YM wind field model,extreme wind speeds were extracted from the virtual typhoons.The extreme wind speeds for different return periods were predicted and compared with the current structural code to provide improved wind load information for wind-resistant structural design.展开更多
The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By thi...The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.展开更多
A modified Monte-Carlo(MC) method to simulate the regular growth of binary eutectic alloys is presented. It is found that the growth rate has a linear dependence on the chemical potential difference between the solid ...A modified Monte-Carlo(MC) method to simulate the regular growth of binary eutectic alloys is presented. It is found that the growth rate has a linear dependence on the chemical potential difference between the solid and liquid; the relation between the lamellar spacing A and growth rate R accords well with the prediction of Jackson-Hunt(JH) theory unless the growth rate is very low.展开更多
Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with pr...Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with practical values. The migrations of adatoms in coating's growth are different from that of solidification of amorphous materials. In some cases, the migrated adatoms in the process of growth of amorphous coating are not enough to occupy all vacancies and traps, so the amorphous coating is micro-porous. The immovable probability k and the largest migration distance of adatoms, which lie on the electroless bath components, affect the PDF, volume density and microporosity remarkably.展开更多
Usually, there are several methods, e.g. experiment, interpolation experiment-based, analytic function, and Monte-Carlo simulation, to calculate the response functions in LaBr3(Ce) detectors. In logging applications...Usually, there are several methods, e.g. experiment, interpolation experiment-based, analytic function, and Monte-Carlo simulation, to calculate the response functions in LaBr3(Ce) detectors. In logging applications, the experiment-based methods cannot be adopted because of their limitations. Analytic function has the advantage of fast calculating speed, but it is very difficult to take into account many effects that occur in practical applications. On the contrary, Monte-Carlo simulation can deal with physical and geometric configurations very tactfully. It has a distinct advantage for calculating the functions with complex configurations in borehole. A new application of LaBr3(Ce) detector is in natural gamma-rays borehole spectrometer for uranium well logging. Calculation of response functions must consider a series of physical and geometric factors under complex logging conditions, including earth formations and its relevant parameters, different energies, material and thickness of the casings, the fluid between the two tubes, and relative position of the LaBr3(Ce) crystal to steel ingot at the front of logging tube. The present work establishes Monte-Carlo simulation models for the above-mentioned situations, and then performs calculations for main gamma-rays from natural radio-elements series. The response functions can offer experimental directions for the design of borehole detection system, and provide technique basis and basic data for spectral analysis of natural gamma-rays, and for sonrceless calibration in uranium quantitative interpretation.展开更多
For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo(MC) techniques are considered to be more effective than other algorithms, such as finite element method or f...For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo(MC) techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.展开更多
The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-se...The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-section dependence on ion velocity and random generation of ion velocities and scattering angles after collisions. The developed algorithm of simulation of interval between collisions takes into account the change of relative velocity of ion-atom pair as well as the change of cross-section of collision and atomic concentration. At the same time, unlike the widely used “null-collision” method, both the probability of collision and change of particles’ state which determines this probability are taken into consideration for each particle independently in time. The simulation results according to the techniques proposed are found to be close to the theoretical values of ion drift velocities. It is revealed that the “null-collision” method results in exceeding of drift velocity in strong and intermediate fields. At the same time the proposed method of accumulation of probability under the same conditions gives values close to theoretical ones. In weak fields calculated values of drift velocity in both methods exceed theoretical values to some small extent.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.12002246 and No.52178301)Knowledge Innovation Program of Wuhan(Grant No.2022010801020357)+2 种基金the Science Research Foundation of Wuhan Institute of Technology(Grant No.K2021030)2020 annual Open Fund of Failure Mechanics&Engineering Disaster Prevention and Mitigation,Key Laboratory of Sichuan Province(Sichuan University)(Grant No.2020JDS0022)Open Research Fund Program of Hubei Provincial Key Laboratory of Chemical Equipment Intensification and Intrinsic Safety(Grant No.2019KA03)。
文摘This paper proposes an active learning accelerated Monte-Carlo simulation method based on the modified K-nearest neighbors algorithm.The core idea of the proposed method is to judge whether or not the output of a random input point can be postulated through a classifier implemented through the modified K-nearest neighbors algorithm.Compared to other active learning methods resorting to experimental designs,the proposed method is characterized by employing Monte-Carlo simulation for sampling inputs and saving a large portion of the actual evaluations of outputs through an accurate classification,which is applicable for most structural reliability estimation problems.Moreover,the validity,efficiency,and accuracy of the proposed method are demonstrated numerically.In addition,the optimal value of K that maximizes the computational efficiency is studied.Finally,the proposed method is applied to the reliability estimation of the carbon fiber reinforced silicon carbide composite specimens subjected to random displacements,which further validates its practicability.
基金Supported by the National Key Research and Development Program of China(Nos.2016YFC1402004,2016YFC1402000,2018YFC1407003)the National Natural Science Foundation of China(Nos.U1706216,U1606402,41421005)
文摘As one of the most serious natural disasters,many typhoons affect southeastern China every year.Taking Shenzhen,a coastal city in southeast China as an example,we employed a Monte-Carlo simulation to generate a large number of virtual typhoons for wind hazard analysis.By analyzing 67-year historical typhoons data from 1949 to 2015 using the Best Track Dataset for Tropical Cyclones over the Western North Pacific recorded by the Shanghai Typhoon Institute,China Meteorological Administration(CMASTI),typhoon characteristic parameters were extracted and optimal statistical distributions established for the parameters in relation to Shenzhen.We employed the Monte-Carlo method to sample each distribution to generate the characteristic parameters of virtual typhoons.In addition,the Yah Meng(YM)wind field model was introduced,and the sensitivity of the YM model to several parameters discussed.Using the YM wind field model,extreme wind speeds were extracted from the virtual typhoons.The extreme wind speeds for different return periods were predicted and compared with the current structural code to provide improved wind load information for wind-resistant structural design.
文摘The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.
基金The authors wish to acknowledge the financial support from the Post-doctoral Foundation of Laboratory of Solid State Microstruc-ture of Nanjing University, and the National Natural Science Foun-dation of China. (No.10074028)
文摘A modified Monte-Carlo(MC) method to simulate the regular growth of binary eutectic alloys is presented. It is found that the growth rate has a linear dependence on the chemical potential difference between the solid and liquid; the relation between the lamellar spacing A and growth rate R accords well with the prediction of Jackson-Hunt(JH) theory unless the growth rate is very low.
文摘Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with practical values. The migrations of adatoms in coating's growth are different from that of solidification of amorphous materials. In some cases, the migrated adatoms in the process of growth of amorphous coating are not enough to occupy all vacancies and traps, so the amorphous coating is micro-porous. The immovable probability k and the largest migration distance of adatoms, which lie on the electroless bath components, affect the PDF, volume density and microporosity remarkably.
基金supported by Key Laboratory of Radioactive Geology and Exploration Technology Fundamental Science for National Defense(No. 2011RGET04)East China Institute of Technology, and National Natural Science Foundation of China (No. 41074078)
文摘Usually, there are several methods, e.g. experiment, interpolation experiment-based, analytic function, and Monte-Carlo simulation, to calculate the response functions in LaBr3(Ce) detectors. In logging applications, the experiment-based methods cannot be adopted because of their limitations. Analytic function has the advantage of fast calculating speed, but it is very difficult to take into account many effects that occur in practical applications. On the contrary, Monte-Carlo simulation can deal with physical and geometric configurations very tactfully. It has a distinct advantage for calculating the functions with complex configurations in borehole. A new application of LaBr3(Ce) detector is in natural gamma-rays borehole spectrometer for uranium well logging. Calculation of response functions must consider a series of physical and geometric factors under complex logging conditions, including earth formations and its relevant parameters, different energies, material and thickness of the casings, the fluid between the two tubes, and relative position of the LaBr3(Ce) crystal to steel ingot at the front of logging tube. The present work establishes Monte-Carlo simulation models for the above-mentioned situations, and then performs calculations for main gamma-rays from natural radio-elements series. The response functions can offer experimental directions for the design of borehole detection system, and provide technique basis and basic data for spectral analysis of natural gamma-rays, and for sonrceless calibration in uranium quantitative interpretation.
基金supported in part by the Higher Education Commission of Pakistan under PPCR programsupported by the National Magnetic Confinement Fusion Program under Grant No.2013GB104004Fundamental Research Fund for Chinese Central Universities
文摘For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo(MC) techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.
文摘The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-section dependence on ion velocity and random generation of ion velocities and scattering angles after collisions. The developed algorithm of simulation of interval between collisions takes into account the change of relative velocity of ion-atom pair as well as the change of cross-section of collision and atomic concentration. At the same time, unlike the widely used “null-collision” method, both the probability of collision and change of particles’ state which determines this probability are taken into consideration for each particle independently in time. The simulation results according to the techniques proposed are found to be close to the theoretical values of ion drift velocities. It is revealed that the “null-collision” method results in exceeding of drift velocity in strong and intermediate fields. At the same time the proposed method of accumulation of probability under the same conditions gives values close to theoretical ones. In weak fields calculated values of drift velocity in both methods exceed theoretical values to some small extent.