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.展开更多
现有安全稳定控制系统(简称稳控系统)的可靠性评估方法本质上属于静态建模,由于未能体现系统内各装置老化和检修等动态过程,在一定程度上影响了评估结果的准确性。为此,文中提出一种基于马尔可夫链蒙特卡洛(Markov chain Monte Carlo,MC...现有安全稳定控制系统(简称稳控系统)的可靠性评估方法本质上属于静态建模,由于未能体现系统内各装置老化和检修等动态过程,在一定程度上影响了评估结果的准确性。为此,文中提出一种基于马尔可夫链蒙特卡洛(Markov chain Monte Carlo,MCMC)的稳控系统动态可靠性评估方法。首先针对失效过程,构建四状态非齐次马尔可夫模型来模拟装置老化过程,并给出各状态评判方法;其次针对修复过程,分析不同检修策略对装置状态转移的影响以体现状态检修的差异性;最后考虑稳控装置状态转移过程的时序或条件相关性,对稳控系统可靠性进行动态建模。以实际稳控系统为例,仿真对比不同检修策略下的可靠性,并对模型参数进行灵敏度分析。评估结果表明,该方法可以求解稳控系统的时变可用度,用于指导稳控装置现场合理检修。展开更多
Background Many phenotypes in animal breeding are derived from incomplete measures,especially if they are challenging or expensive to measure precisely.Examples include time-dependent traits such as reproductive statu...Background Many phenotypes in animal breeding are derived from incomplete measures,especially if they are challenging or expensive to measure precisely.Examples include time-dependent traits such as reproductive status,or lifespan.Incomplete measures for these traits result in phenotypes that are subject to left-,interval-and rightcensoring,where phenotypes are only known to fall below an upper bound,between a lower and upper bound,or above a lower bound respectively.Here we compare three methods for deriving phenotypes from incomplete data using age at first elevation(>1 ng/mL)in blood plasma progesterone(AGEP4),which generally coincides with onset of puberty,as an example trait.Methods We produced AGEP4 phenotypes from three blood samples collected at about 30-day intervals from approximately 5,000 Holstein–Friesian or Holstein–Friesian×Jersey cross-bred dairy heifers managed in 54 seasonal-calving,pasture-based herds in New Zealand.We used these actual data to simulate 7 different visit scenarios,increasing the extent of censoring by disregarding data from one or two of the three visits.Three methods for deriving phenotypes from these data were explored:1)ordinal categorical variables which were analysed using categorical threshold analysis;2)continuous variables,with a penalty of 31 d assigned to right-censored phenotypes;and 3)continuous variables,sampled from within a lower and upper bound using a data augmentation approach.Results Credibility intervals for heritability estimations overlapped across all methods and visit scenarios,but estimated heritabilities tended to be higher when left censoring was reduced.For sires with at least 5 daughters,the correlations between estimated breeding values(EBVs)from our three-visit scenario and each reduced data scenario varied by method,ranging from 0.65 to 0.95.The estimated breed effects also varied by method,but breed differences were smaller as phenotype censoring increased.Conclusion Our results indicate that using some methods,phenotypes derived from one observation per offspring for a time-dependent trait such as AGEP4 may provide comparable sire rankings to three observations per offspring.This has implications for the design of large-scale phenotyping initiatives where animal breeders aim to estimate variance parameters and estimated breeding values(EBVs)for phenotypes that are challenging to measure or prohibitively expensive.展开更多
食物网结构特征和能量流动的研究,对于维持海洋生态系统结构和功能的稳定具有重要意义,有助于深入理解海洋生态系统的复杂过程。本研究基于2019-2021年在江苏近海北部海域开展的季节性渔业资源底拖网调查数据,通过构建基于蒙特卡罗马尔...食物网结构特征和能量流动的研究,对于维持海洋生态系统结构和功能的稳定具有重要意义,有助于深入理解海洋生态系统的复杂过程。本研究基于2019-2021年在江苏近海北部海域开展的季节性渔业资源底拖网调查数据,通过构建基于蒙特卡罗马尔科夫链算法的逆线性模型(Linear Inverse Models using a Monte Carlo Method Coupled with Markov Chain, LIM-MCMC),结合生态网络分析(Ecological Network Analysis,ENA)的方法,分析了该海域生态系统状态和食物网能量流动特征,旨在为江苏近海北部海域食物网营养动力学研究提供参考依据。结果表明,该海域生态系统共包含299条能量流动路径,能量流动分布整体呈典型的金字塔结构,各功能群呼吸消耗和流入有机碎屑的能量保持同步性。通过与其他海域比较发现,江苏近海北部海域生态系统的连接指数(Connectance,C)和系统杂食指数(System Omnivory Index,SOI)分别为0.40和0.22,处于较高水平,表明该生态系统不同营养级间的营养联系较为紧密,食物网结构相对复杂,能够在较大程度上抵御外界扰动。总初级生产力/总呼吸(Total Primary Production/Total Respiration,TPP/TR)和Finn’s循环指数(Finn’s Cycling Index,FCI)分别为1.05和5.76%,表明该生态系统对能量利用效率较高。此外,约束效率(Constraint Efficiency,CE)、发展程度(Extent of Development,AC)、协同效应指数(Synergism Index,b/c)和主导间接效应(Dominance Indirect Effects,i/d)也表明该生态系统具有较高的系统发展程度、再生潜力和系统发展空间。本研究将有助于为江苏近海北部海域生态系统的修复和渔业资源的可持续利用提供理论基础,为实施基于生态系统的渔业管理提供科学依据。展开更多
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 resul...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.展开更多
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 exten...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.展开更多
Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random error...Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system's dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error(DTE) fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system's dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.展开更多
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 t...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.展开更多
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 developi...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.展开更多
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 sam...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.展开更多
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 accelerati...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.展开更多
The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detect...The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detector distance in a high-purity germanium detector system was studied using a Geant4 Monte Carlo simulation.The effect of the deadlayer in the germanium crystal was estimated by reference to experiments on 241 Am and the relative efficiency of the detector.The peak efficiency dependence on the source-to-detector distance was compared between the simulation and measurements.The modified sum-peak method is discussed with respect to these peak efficiencies.展开更多
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 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.展开更多
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...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.展开更多
An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face ...An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face is numerically generated with an ocean waveheight spectrum for a given wind velocity. A two-layered medium is then generated by adding a thin oil layer on the simulated rough sea surface. The electric fields backscattered from the sea surface with two-layered medium are computed with the method of moments (MoM), and the backscattering coefficients are statistically obtained with N independent samples for each oil-spilled surface using the Monte-Carlo technique for various conditions of surface roughness, oil-layer thickness, frequency, polarization and incidence angle. The numerical simulation results are compared with theoretical models for clean sea surfaces and SAR images of an off-spilled sea surface caused by the Hebei (Hebei province, China) Spirit oil tanker in 2007. Further, conditions for better oil spill extraction are sought by the numerical simulation on the effects of wind speed and oil-layer thickness at different inci- dence angles on the backscattering coefficients.展开更多
基金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.
文摘现有安全稳定控制系统(简称稳控系统)的可靠性评估方法本质上属于静态建模,由于未能体现系统内各装置老化和检修等动态过程,在一定程度上影响了评估结果的准确性。为此,文中提出一种基于马尔可夫链蒙特卡洛(Markov chain Monte Carlo,MCMC)的稳控系统动态可靠性评估方法。首先针对失效过程,构建四状态非齐次马尔可夫模型来模拟装置老化过程,并给出各状态评判方法;其次针对修复过程,分析不同检修策略对装置状态转移的影响以体现状态检修的差异性;最后考虑稳控装置状态转移过程的时序或条件相关性,对稳控系统可靠性进行动态建模。以实际稳控系统为例,仿真对比不同检修策略下的可靠性,并对模型参数进行灵敏度分析。评估结果表明,该方法可以求解稳控系统的时变可用度,用于指导稳控装置现场合理检修。
基金funded by New Zealand dairy farmers through Dairy NZ Inc. and by the New Zealand Ministry of Business,Innovation and Employment (DRCX1302)support was kindly received from LIC and CRV
文摘Background Many phenotypes in animal breeding are derived from incomplete measures,especially if they are challenging or expensive to measure precisely.Examples include time-dependent traits such as reproductive status,or lifespan.Incomplete measures for these traits result in phenotypes that are subject to left-,interval-and rightcensoring,where phenotypes are only known to fall below an upper bound,between a lower and upper bound,or above a lower bound respectively.Here we compare three methods for deriving phenotypes from incomplete data using age at first elevation(>1 ng/mL)in blood plasma progesterone(AGEP4),which generally coincides with onset of puberty,as an example trait.Methods We produced AGEP4 phenotypes from three blood samples collected at about 30-day intervals from approximately 5,000 Holstein–Friesian or Holstein–Friesian×Jersey cross-bred dairy heifers managed in 54 seasonal-calving,pasture-based herds in New Zealand.We used these actual data to simulate 7 different visit scenarios,increasing the extent of censoring by disregarding data from one or two of the three visits.Three methods for deriving phenotypes from these data were explored:1)ordinal categorical variables which were analysed using categorical threshold analysis;2)continuous variables,with a penalty of 31 d assigned to right-censored phenotypes;and 3)continuous variables,sampled from within a lower and upper bound using a data augmentation approach.Results Credibility intervals for heritability estimations overlapped across all methods and visit scenarios,but estimated heritabilities tended to be higher when left censoring was reduced.For sires with at least 5 daughters,the correlations between estimated breeding values(EBVs)from our three-visit scenario and each reduced data scenario varied by method,ranging from 0.65 to 0.95.The estimated breed effects also varied by method,but breed differences were smaller as phenotype censoring increased.Conclusion Our results indicate that using some methods,phenotypes derived from one observation per offspring for a time-dependent trait such as AGEP4 may provide comparable sire rankings to three observations per offspring.This has implications for the design of large-scale phenotyping initiatives where animal breeders aim to estimate variance parameters and estimated breeding values(EBVs)for phenotypes that are challenging to measure or prohibitively expensive.
文摘食物网结构特征和能量流动的研究,对于维持海洋生态系统结构和功能的稳定具有重要意义,有助于深入理解海洋生态系统的复杂过程。本研究基于2019-2021年在江苏近海北部海域开展的季节性渔业资源底拖网调查数据,通过构建基于蒙特卡罗马尔科夫链算法的逆线性模型(Linear Inverse Models using a Monte Carlo Method Coupled with Markov Chain, LIM-MCMC),结合生态网络分析(Ecological Network Analysis,ENA)的方法,分析了该海域生态系统状态和食物网能量流动特征,旨在为江苏近海北部海域食物网营养动力学研究提供参考依据。结果表明,该海域生态系统共包含299条能量流动路径,能量流动分布整体呈典型的金字塔结构,各功能群呼吸消耗和流入有机碎屑的能量保持同步性。通过与其他海域比较发现,江苏近海北部海域生态系统的连接指数(Connectance,C)和系统杂食指数(System Omnivory Index,SOI)分别为0.40和0.22,处于较高水平,表明该生态系统不同营养级间的营养联系较为紧密,食物网结构相对复杂,能够在较大程度上抵御外界扰动。总初级生产力/总呼吸(Total Primary Production/Total Respiration,TPP/TR)和Finn’s循环指数(Finn’s Cycling Index,FCI)分别为1.05和5.76%,表明该生态系统对能量利用效率较高。此外,约束效率(Constraint Efficiency,CE)、发展程度(Extent of Development,AC)、协同效应指数(Synergism Index,b/c)和主导间接效应(Dominance Indirect Effects,i/d)也表明该生态系统具有较高的系统发展程度、再生潜力和系统发展空间。本研究将有助于为江苏近海北部海域生态系统的修复和渔业资源的可持续利用提供理论基础,为实施基于生态系统的渔业管理提供科学依据。
基金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 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.
基金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).
文摘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.
基金Supported by National Natural Science Foundation of China(Grant No.51375013)Anhui Provincial Natural Science Foundation of China(Grant No.1208085ME64)
文摘Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system's dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error(DTE) fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system's dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.
文摘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.
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.11775126,11545013,11605101)the Young Elite Scientists Sponsorship Program by CAST(No.2016QNRC001)+1 种基金Science Challenge Project by MIIT of China(No.TZ2018001)Tsinghua University,Initiative Scientific Research Program。
文摘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.
文摘The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detector distance in a high-purity germanium detector system was studied using a Geant4 Monte Carlo simulation.The effect of the deadlayer in the germanium crystal was estimated by reference to experiments on 241 Am and the relative efficiency of the detector.The peak efficiency dependence on the source-to-detector distance was compared between the simulation and measurements.The modified sum-peak method is discussed with respect to these peak efficiencies.
基金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 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.
基金Supported by the National Natural Science Foundation of China under Grant No 61572529
文摘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.
基金The Project "Development of Korea Operational Oceanographic System (PM57041)" funded by the Ministry of Land, Transport and Maritime Affairs of Korean Governmentthe Project "Cooperation on the Development of Basic Technologies for the Yellow Sea and East China Sea Operational Oceanographic System (YOOS)" funded by CKJORC and the Basic Research Projects (PE98731, PG47770 and PE98732) of the Korea Institute Ocean Science and Technologysupport by the PASCO Corporation,Japan is also apreciated
文摘An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face is numerically generated with an ocean waveheight spectrum for a given wind velocity. A two-layered medium is then generated by adding a thin oil layer on the simulated rough sea surface. The electric fields backscattered from the sea surface with two-layered medium are computed with the method of moments (MoM), and the backscattering coefficients are statistically obtained with N independent samples for each oil-spilled surface using the Monte-Carlo technique for various conditions of surface roughness, oil-layer thickness, frequency, polarization and incidence angle. The numerical simulation results are compared with theoretical models for clean sea surfaces and SAR images of an off-spilled sea surface caused by the Hebei (Hebei province, China) Spirit oil tanker in 2007. Further, conditions for better oil spill extraction are sought by the numerical simulation on the effects of wind speed and oil-layer thickness at different inci- dence angles on the backscattering coefficients.