The Monte Carlo(MC)method offers significant advantages in handling complex geometries and physical processes in particle transport problems and has become a widely used approach in reactor physics analysis,radiation ...The Monte Carlo(MC)method offers significant advantages in handling complex geometries and physical processes in particle transport problems and has become a widely used approach in reactor physics analysis,radiation shielding design,and medical physics.However,with the rapid advancement of new nuclear energy systems,the Monte Carlo method faces challenges in efficiency,accuracy,and adaptability,limiting its effectiveness in meeting modern design requirements.Overcoming technical obstacles related to high-fidelity coupling,high-resolution computation,and intelligent design is essential for using the Monte Carlo method as a reliable tool in numerical analysis for these new nuclear energy systems.To address these challenges,the Nuclear Energy and Application Laboratory(NEAL)team at the University of South China developed a multifunctional and generalized intelligent code platform called MagicMC,based on the Monte Carlo particle transport method.MagicMC is a developing tool dedicated to nuclear applications,incorporating intelligent methodologies.It consists of two primary components:a basic unit and a functional unit.The basic unit,which functions similarly to a standard Monte Carlo particle transport code,includes seven modules:geometry,source,transport,database,tally,output,and auxiliary.The functional unit builds on the basic unit by adding functional modules to address complex and diverse applications in nuclear analysis.MagicMC introduces a dynamic Monte Carlo particle transport algorithm to address time-space particle transport problems within emerging nuclear energy systems and incorporates a CPU-GPU heterogeneous parallel framework to enable high-efficiency,high-resolution simulations for large-scale computational problems.Anticipating future trends in intelligent design,MagicMC integrates several advanced features,including CAD-based geometry modeling,global variance reduction methods,multi-objective shielding optimization,high-resolution activation analysis,multi-physics coupling,and radiation therapy.In this paper,various numerical benchmarks-spanning reactor transient simulations,material activation analysis,radiation shielding optimization,and medical dosimetry analysis-are presented to validate MagicMC.The numerical results demonstrate MagicMC's efficiency,accuracy,and reliability in these preliminary applications,underscoring its potential to support technological advancements in developing high-fidelity,high-resolution,and high-intelligence MC-based tools for advanced nuclear applications.展开更多
Scramjet is the most promising propulsion system for Air-breathing Hypersonic Vehicle(AHV),and the Infrared(IR)radiation it emits is critical for early warning,detection,and identification of such weapons.This work pr...Scramjet is the most promising propulsion system for Air-breathing Hypersonic Vehicle(AHV),and the Infrared(IR)radiation it emits is critical for early warning,detection,and identification of such weapons.This work proposes an Adaptive Reverse Monte Carlo(ARMC)method and develops an analytical model for the IR radiation of scramjet considering gaseous kerosene and hydrogen fueled conditions.The evaluation studies show that at a global equivalence ratio of 0.8,the IR radiation from hydrogen-fueled plume is predominantly from H_(2)O and spectral peak is 1.53 kW·Sr^(-1)·μm^(-1)at the 2.7μm band,while the kerosene-fueled plume exhibits a spectral intensity approaching 7.0 kW·Sr^(-1)·μm^(-1)at the 4.3μm band.At the backward detection angle,both types of scramjets exhibit spectral peaks within the 1.3-1.4μm band,with intensities around10 kW·Sr^(-1)·μm^(-1).The integral radiation intensity of hydrogen-fueled scramjet is generally higher than kerosene-fueled scramjet,particularly in 1-3μm band.Meanwhile,at wide detection angles,the solid walls become the predominant radiation source.The radiation intensity is highest in1-3μm and weakest in 8-14μm band,with values of 21.5 kW·Sr^(-1)and 0.57 kW·Sr^(-1)at the backward detection angles,respectively.Significant variations in the radiation contributions from gases and solids are observed across different bands under the two fuel conditions,especially within 3-5μm band.This research provides valuable insights into the IR radiation characteristics of scramjets,which can aid in the development of IR detection systems for AHV.展开更多
Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patte...Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.展开更多
In contrast to conventional reservoirs,tight formations have more complex pore structures and significant boundary layer effect,making it difficult to determine the effective permeability.To address this,this paper fi...In contrast to conventional reservoirs,tight formations have more complex pore structures and significant boundary layer effect,making it difficult to determine the effective permeability.To address this,this paper first proposes a semi-empirical model for calculating boundary layer thickness based on dimensional analysis,using published experimental data on microcapillary flow.Furthermore,considering the non-uniform distribution of fluid viscosity in the flow channels of tight reservoirs,a theoretical model for boundary layer thickness is established based on fractal theory,and permeability predictions are conducted through Monte Carlo simulations.Finally,sensitivity analyses of various influencing parameters are performed.The results show that,compared to other fractal-based analytical models,the proposed permeability probabilistic model integrates parameters affecting fluid flow with random numbers,reflecting both the fractal and randomness characteristics of capillary size distribution.The computational results exhibit the highest consistency with experimental data.Among the factors affecting the boundary layer,in addition to certain conventional physical and mechanical parameters,different microstructure parameters significantly influence the boundary layer as well.A higher tortuosity fractal dimension results in a thicker boundary layer,while increases in pore fractal dimension,porosity,and maximum capillary size help mitigate the boundary layer effect.It is also observed that the permeability of large pores exhibits greater sensitivity to changes in various influencing parameters.Considering micro-scale flow effects,the proposed model enhances the understanding of the physical mechanisms of fluid transport in dense porous media.展开更多
We extend a semiclassical numerical method, bosonic auxiliary-field Monte Carlo, to quantum spin systems. This method breaks the lattice into clusters, solves each cluster precisely and couples them with classical aux...We extend a semiclassical numerical method, bosonic auxiliary-field Monte Carlo, to quantum spin systems. This method breaks the lattice into clusters, solves each cluster precisely and couples them with classical auxiliary fields through classical Monte Carlo simulation. We test the method with antiferromagnetic spin models in one-dimensional chains, square lattices and triangular lattices, and obtain reasonable results at finite temperatures. This algorithm builds a bridge between classical Monte Carlo method and quantum methods. The algorithm can be improved with either progress in classical Monte Carlo sampling or the development of quantum solvers, and can also be further applied to systems with different lattices or interactions.展开更多
针对多部位损伤(Multiple Site Damage,MSD)结构安全性评估问题,通过Monte-Carlo方法对MSD结构的失效概率进行预测和分析。首先,基于多孔铝板的多裂纹萌生试验,得出裂纹萌生寿命服从对数正态分布,为多裂纹萌生分析提供支持;通过多孔铝...针对多部位损伤(Multiple Site Damage,MSD)结构安全性评估问题,通过Monte-Carlo方法对MSD结构的失效概率进行预测和分析。首先,基于多孔铝板的多裂纹萌生试验,得出裂纹萌生寿命服从对数正态分布,为多裂纹萌生分析提供支持;通过多孔铝板的剩余强度试验得到铆钉孔直径、铆钉孔间距和裂纹萌生位置对结构剩余强度均有一定影响。其次,通过对裂纹萌生寿命分布进行随机抽样生成初始裂纹并使用组合法结合Paris公式,实现多裂纹随机扩展的模拟;在试验数据基础上,对传统的Irwin塑性区连通准则进行改进,发现改进的Irwin塑性区连通准则在孔间距大于10mm时的误差大大降低,并结合净截面屈服准则以获得更好的剩余强度预测结果;将随机性的裂纹萌生和扩展过程与确定性的剩余强度预测方法相结合,建立基于Monte-Carlo方法的MSD结构的失效概率预测模型。最后,通过算例分析,该模型能够得到MSD结构的失效概率曲线,实现结构安全性评估。结果表明MSD结构的失效概率会在短时间内迅速增加,需要在裂纹萌生寿命附近进行限制。展开更多
Considering the tunneling effect and the Schottky effect,the metal semiconductor contact is simulated by using self consistent ensemble Monte Carlo method.Under different biases or at different barrier heights,the i...Considering the tunneling effect and the Schottky effect,the metal semiconductor contact is simulated by using self consistent ensemble Monte Carlo method.Under different biases or at different barrier heights,the investigation into the tunneling current indicates that the tunneling effect is of great importance under reverse biases.The Schottky barrier diode current due to Schottky effect is in agreement with the theoretical one.The barrier lowering is found a profound effect on the current transport at the metal semiconductor interface.展开更多
An approach for generating test problems by a computer using the Monte Carlo method based upon user-given characterizations is described.A single point X~* is prespocified by the user to be a solution of the test prob...An approach for generating test problems by a computer using the Monte Carlo method based upon user-given characterizations is described.A single point X~* is prespocified by the user to be a solution of the test problems.The approach is flex- ible enough to specify function values,gradients,Hesse,degeneracy degree and ill- conditioned degree at the point X~*.Other numerical features such as indefiniteness, convexity are also under user's control.展开更多
After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting cont...After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting control moment error is carried out respectively by means of Monte Carlo method and the distribution of impact points of rudder control rocket is counted from all the successful subsamples. In the case of adding interference errors associated with rudder angle error and starting time error, the simulation analysis of impact point dispersion is done and its lateral and longitudinal correction abilities at different targeting angles are simulated to identify the effects of these factors on characteristics and control precision of the rudder control rocket, which provides the relevant reference for high-precision design of rudder control system.展开更多
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.展开更多
The conduct mechanism of the doped polymer is considered. In an asymmetrysystem composed of high polymer and doping conductive matte, chain or congeries framework will beformed between the conductive particles to impr...The conduct mechanism of the doped polymer is considered. In an asymmetrysystem composed of high polymer and doping conductive matte, chain or congeries framework will beformed between the conductive particles to improve the conductance characteristic. In thisprocession, the conductive particles interact to each other. In this paper, we describe theconductance of the doped polymer by Monte Carlo method. The results accord with the experimentsquite well. It can be concluded that there is an evident change of doped polymer from nonconductorto metal.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China(Nos.12475174 and U2267207)YueLuShan Center Industrial Innovation(No.2024YCII0108)+2 种基金Natural Science Foundation of Hunan Province(No.2022JJ40345)Science and Technology Innovation Project of Hengyang(No.202250045336)the Project of State Key Laboratory of Radiation Medicine and Protection,Soochow University(No.GZK12023031)。
文摘The Monte Carlo(MC)method offers significant advantages in handling complex geometries and physical processes in particle transport problems and has become a widely used approach in reactor physics analysis,radiation shielding design,and medical physics.However,with the rapid advancement of new nuclear energy systems,the Monte Carlo method faces challenges in efficiency,accuracy,and adaptability,limiting its effectiveness in meeting modern design requirements.Overcoming technical obstacles related to high-fidelity coupling,high-resolution computation,and intelligent design is essential for using the Monte Carlo method as a reliable tool in numerical analysis for these new nuclear energy systems.To address these challenges,the Nuclear Energy and Application Laboratory(NEAL)team at the University of South China developed a multifunctional and generalized intelligent code platform called MagicMC,based on the Monte Carlo particle transport method.MagicMC is a developing tool dedicated to nuclear applications,incorporating intelligent methodologies.It consists of two primary components:a basic unit and a functional unit.The basic unit,which functions similarly to a standard Monte Carlo particle transport code,includes seven modules:geometry,source,transport,database,tally,output,and auxiliary.The functional unit builds on the basic unit by adding functional modules to address complex and diverse applications in nuclear analysis.MagicMC introduces a dynamic Monte Carlo particle transport algorithm to address time-space particle transport problems within emerging nuclear energy systems and incorporates a CPU-GPU heterogeneous parallel framework to enable high-efficiency,high-resolution simulations for large-scale computational problems.Anticipating future trends in intelligent design,MagicMC integrates several advanced features,including CAD-based geometry modeling,global variance reduction methods,multi-objective shielding optimization,high-resolution activation analysis,multi-physics coupling,and radiation therapy.In this paper,various numerical benchmarks-spanning reactor transient simulations,material activation analysis,radiation shielding optimization,and medical dosimetry analysis-are presented to validate MagicMC.The numerical results demonstrate MagicMC's efficiency,accuracy,and reliability in these preliminary applications,underscoring its potential to support technological advancements in developing high-fidelity,high-resolution,and high-intelligence MC-based tools for advanced nuclear applications.
基金supported by the National Natural Science Foundation of China(No.12102356)。
文摘Scramjet is the most promising propulsion system for Air-breathing Hypersonic Vehicle(AHV),and the Infrared(IR)radiation it emits is critical for early warning,detection,and identification of such weapons.This work proposes an Adaptive Reverse Monte Carlo(ARMC)method and develops an analytical model for the IR radiation of scramjet considering gaseous kerosene and hydrogen fueled conditions.The evaluation studies show that at a global equivalence ratio of 0.8,the IR radiation from hydrogen-fueled plume is predominantly from H_(2)O and spectral peak is 1.53 kW·Sr^(-1)·μm^(-1)at the 2.7μm band,while the kerosene-fueled plume exhibits a spectral intensity approaching 7.0 kW·Sr^(-1)·μm^(-1)at the 4.3μm band.At the backward detection angle,both types of scramjets exhibit spectral peaks within the 1.3-1.4μm band,with intensities around10 kW·Sr^(-1)·μm^(-1).The integral radiation intensity of hydrogen-fueled scramjet is generally higher than kerosene-fueled scramjet,particularly in 1-3μm band.Meanwhile,at wide detection angles,the solid walls become the predominant radiation source.The radiation intensity is highest in1-3μm and weakest in 8-14μm band,with values of 21.5 kW·Sr^(-1)and 0.57 kW·Sr^(-1)at the backward detection angles,respectively.Significant variations in the radiation contributions from gases and solids are observed across different bands under the two fuel conditions,especially within 3-5μm band.This research provides valuable insights into the IR radiation characteristics of scramjets,which can aid in the development of IR detection systems for AHV.
基金the project PID2022-139202OB-I00Neural Networks and Optimization Techniques for the Design and Safe Maintenance of Transportation Infrastructures:Volcanic Rock Geotechnics and Slope Stability(IA-Pyroslope),funded by the Spanish State Research Agency of the Ministry of Science,Innovation and Universities of Spain and the European Regional Development Fund,MCIN/AEI/10.13039/501100011033/FEDER,EU。
文摘Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.
基金supported by the Hebei Provincial Natural Science Foundation of China(No.D2023402012)the Major Science and Technology Project of China National Petroleum Corporation(No.2024DJ87).
文摘In contrast to conventional reservoirs,tight formations have more complex pore structures and significant boundary layer effect,making it difficult to determine the effective permeability.To address this,this paper first proposes a semi-empirical model for calculating boundary layer thickness based on dimensional analysis,using published experimental data on microcapillary flow.Furthermore,considering the non-uniform distribution of fluid viscosity in the flow channels of tight reservoirs,a theoretical model for boundary layer thickness is established based on fractal theory,and permeability predictions are conducted through Monte Carlo simulations.Finally,sensitivity analyses of various influencing parameters are performed.The results show that,compared to other fractal-based analytical models,the proposed permeability probabilistic model integrates parameters affecting fluid flow with random numbers,reflecting both the fractal and randomness characteristics of capillary size distribution.The computational results exhibit the highest consistency with experimental data.Among the factors affecting the boundary layer,in addition to certain conventional physical and mechanical parameters,different microstructure parameters significantly influence the boundary layer as well.A higher tortuosity fractal dimension results in a thicker boundary layer,while increases in pore fractal dimension,porosity,and maximum capillary size help mitigate the boundary layer effect.It is also observed that the permeability of large pores exhibits greater sensitivity to changes in various influencing parameters.Considering micro-scale flow effects,the proposed model enhances the understanding of the physical mechanisms of fluid transport in dense porous media.
基金supports by the National Key Research and Development Program of China (Grant No.2024YFA1409200)the National Natural Science Foundation of China (Grant Nos.12222412 and 12047503)+1 种基金CAS Project for Young Scientists in Basic Research (Grant No.YSBR-057)supports by the National Natural Science Foundation of China (Grant No.12374144)。
文摘We extend a semiclassical numerical method, bosonic auxiliary-field Monte Carlo, to quantum spin systems. This method breaks the lattice into clusters, solves each cluster precisely and couples them with classical auxiliary fields through classical Monte Carlo simulation. We test the method with antiferromagnetic spin models in one-dimensional chains, square lattices and triangular lattices, and obtain reasonable results at finite temperatures. This algorithm builds a bridge between classical Monte Carlo method and quantum methods. The algorithm can be improved with either progress in classical Monte Carlo sampling or the development of quantum solvers, and can also be further applied to systems with different lattices or interactions.
文摘针对多部位损伤(Multiple Site Damage,MSD)结构安全性评估问题,通过Monte-Carlo方法对MSD结构的失效概率进行预测和分析。首先,基于多孔铝板的多裂纹萌生试验,得出裂纹萌生寿命服从对数正态分布,为多裂纹萌生分析提供支持;通过多孔铝板的剩余强度试验得到铆钉孔直径、铆钉孔间距和裂纹萌生位置对结构剩余强度均有一定影响。其次,通过对裂纹萌生寿命分布进行随机抽样生成初始裂纹并使用组合法结合Paris公式,实现多裂纹随机扩展的模拟;在试验数据基础上,对传统的Irwin塑性区连通准则进行改进,发现改进的Irwin塑性区连通准则在孔间距大于10mm时的误差大大降低,并结合净截面屈服准则以获得更好的剩余强度预测结果;将随机性的裂纹萌生和扩展过程与确定性的剩余强度预测方法相结合,建立基于Monte-Carlo方法的MSD结构的失效概率预测模型。最后,通过算例分析,该模型能够得到MSD结构的失效概率曲线,实现结构安全性评估。结果表明MSD结构的失效概率会在短时间内迅速增加,需要在裂纹萌生寿命附近进行限制。
文摘Considering the tunneling effect and the Schottky effect,the metal semiconductor contact is simulated by using self consistent ensemble Monte Carlo method.Under different biases or at different barrier heights,the investigation into the tunneling current indicates that the tunneling effect is of great importance under reverse biases.The Schottky barrier diode current due to Schottky effect is in agreement with the theoretical one.The barrier lowering is found a profound effect on the current transport at the metal semiconductor interface.
文摘An approach for generating test problems by a computer using the Monte Carlo method based upon user-given characterizations is described.A single point X~* is prespocified by the user to be a solution of the test problems.The approach is flex- ible enough to specify function values,gradients,Hesse,degeneracy degree and ill- conditioned degree at the point X~*.Other numerical features such as indefiniteness, convexity are also under user's control.
文摘After the trajectory simulation model of rudder control rocket with six degrees of freedom is established by Matlab/ Simulink, the simulated targeting of rudder control rocket with rudder angle error and starting control moment error is carried out respectively by means of Monte Carlo method and the distribution of impact points of rudder control rocket is counted from all the successful subsamples. In the case of adding interference errors associated with rudder angle error and starting time error, the simulation analysis of impact point dispersion is done and its lateral and longitudinal correction abilities at different targeting angles are simulated to identify the effects of these factors on characteristics and control precision of the rudder control rocket, which provides the relevant reference for high-precision design of rudder control system.
基金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.
文摘The conduct mechanism of the doped polymer is considered. In an asymmetrysystem composed of high polymer and doping conductive matte, chain or congeries framework will beformed between the conductive particles to improve the conductance characteristic. In thisprocession, the conductive particles interact to each other. In this paper, we describe theconductance of the doped polymer by Monte Carlo method. The results accord with the experimentsquite well. It can be concluded that there is an evident change of doped polymer from nonconductorto metal.
基金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.
基金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）.
文摘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.
文摘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.
