基于杂交应力元法的颗粒增强复合材料Monte⁃Carlo模拟

杂交应力元假设的高阶应力场可以用较疏的网格获得较高的计算精度.采用四叉树网格离散非均质计算域,四叉树杂交应力单元悬挂节点的位移协调条件自动满足,且得益于单元类型数量有限,单元刚度矩阵可以预计算,以便在实际计算时直接读取调用,大幅提高了计算效率.考虑夹杂的随机性对颗粒增强复合材料力学性能的影响,采用均匀化方法和Monte⁃Carlo方法,研究了随机夹杂的体积比、数量、长宽比对材料均质等效模量的影响,结果表明,复合材料的等效弹性模量随夹杂体积比、数量、长宽比的增大而增大,且对体积比最敏感.

基于Monte-Carlo仿真的恒应力加速试验方案模拟评价

为了保证可靠性加速试验方案的正确性和可行性,提出基于仿真的模拟评价方法。以最具代表性的Weibull寿命分布产品的恒应力加速寿命试验为例,以代表试验结果准确性和稳定性的指标——中位寿命估计值渐进方差的均值和标准差作为评价指标,通过Monte-Carlo仿真模拟生成失效数据,对不同试验方案的统计精度和稳定性进行评价,为加速寿命试验优化设计的工程实现提供了思路和实例验证。

非对称循环弯曲变形再结晶组织的Monte Carlo模拟

利用Monte Carlo方法模拟了非对称循环弯曲变形工艺下的静态再结晶过程,分析了循环道次、退火温度等工艺参数对板材表面及厚度方向上组织演变和晶粒尺寸分布的影响,并基于模拟数据,建立了静态再结晶晶粒尺寸预测模型。结果表明,在450~750℃的退火温度范围内,随着循环道次增加、退火温度升高以及与板材表面距离的减小,静态再结晶晶粒尺寸显著减小。当达到7个循环道次后,板材表面的晶粒尺寸趋于稳定;而在11个循环道次后,厚度方向的晶粒尺寸均匀化,并细化为等轴晶。此外,通过将模拟结果和金相实验得到的再结晶晶粒尺寸进行对比,模型的预测结果和实验数据吻合度高,所建模型的有效性和准确性得到验证。

基于Monte Carlo随机模拟算法的电动汽车无序充电负荷计算模型

为了有效降低电动汽车无序充电负荷对于电力系统稳定性的影响,研究了基于Monte Carlo随机模拟算法的计算模型。综合考虑用户充电行为各参数的概率密度,使得构建的模型更精确。分析了不同的充电策略对于充电负荷曲线的影响,并利用综合策略模型使得用户的充电成本更低,利用基于Levy飞行策略的改进麻雀算法对综合充电策略的权重系数进行寻优计算,以使得用户的成本和电网稳定均能获得不错的效果。在此基础上,提出利用双层优化模型对模型进一步优化,综合考虑用户满意度和电力波动的情况,结果表明改进后的模型,可以节省充电成本15.9%,降低电力系统负荷峰谷差约10.1%。

Monte Carlo方法在移动通信可靠性分析中的应用

文章主要探讨Monte Carlo方法在移动通信可靠性分析中的应用,并引入重要抽样技术对Monte Carlo方法进行优化。首先,文章介绍Monte Carlo方法的基本原理及其在移动通信可靠性分析中的应用。其次,引入重要抽样技术对该方法进行优化。最后,利用优化网络工程工具(Optimized Network Engineering Tools,OPNET)模拟平台构建相应的移动通信网络模型,并对比Monte Carlo方法和本文研究的优化方法在不同样本数量下的性能差异。结果表明,重要抽样技术的引入显著提高了Monte Carlo方法在移动通信可靠性分析中的准确性和效率,减少了方差,验证了该方法在移动通信可靠性分析中的有效性和实用性。

Exhaustive review of acceleration strategies for Monte Carlo simulations in photon transit

Monte Carlo simulation techniques have become the quintessence and a pivotal nexus of inquiry in the realm of simulating photon movement within biological fabrics.Through the stochastic sampling of tissue archetypes delineated by explicit optical characteristics,Monte Carlo simulations possess the theoretical capacity to render unparalleled accuracy in the depiction of exceedingly intricate phenomena.Nonetheless,the quintessential challenge associated with Monte Carlo simulation methodologies resides in their extended computational duration,which significantly impedes the refinement of their precision.Consequently,this discourse is specifically dedicated to exploring innovations in strategies and technologies aimed at expediting Monte Carlo simulations.It delves into the foundational concepts of various acceleration tactics,evaluates these strategies concerning their speed,accuracy,and practicality,and amalgamates a comprehensive overview and critique of acceleration methodologies for Monte Carlo simulations.Ultimately,the discourse envisages prospective trajectories for the employment of Monte Carlo techniques within the domain of tissue optics.

Research on a Monte Carlo global variance reduction method based on an automatic importance sampling method

Global variance reduction is a bottleneck in Monte Carlo shielding calculations.The global variance reduction problem requires that the statistical error of the entire space is uniform.This study proposed a grid-AIS method for the global variance reduction problem based on the AIS method,which was implemented in the Monte Carlo program MCShield.The proposed method was validated using the VENUS-Ⅲ international benchmark problem and a self-shielding calculation example.The results from the VENUS-Ⅲ benchmark problem showed that the grid-AIS method achieved a significant reduction in the variance of the statistical errors of the MESH grids,decreasing from 1.08×10^(-2) to 3.84×10^(-3),representing a 64.00% reduction.This demonstrates that the grid-AIS method is effective in addressing global issues.The results of the selfshielding calculation demonstrate that the grid-AIS method produced accurate computational results.Moreover,the grid-AIS method exhibited a computational efficiency approximately one order of magnitude higher than that of the AIS method and approximately two orders of magnitude higher than that of the conventional Monte Carlo method.

基于Monte Carlo的风力机叶片强度可靠性分析

随着风力机应用数量增多,以及应用规模扩大,风力机叶片可靠度的保证尤为关键。利用Monte Carlo客观分析叶片强度可靠性,为叶片安全状态长效维持提供依据,实现风力机的平稳运转。该文在介绍Monte Carlo的基础上,分析风力机类型及其叶片状态和叶片参数,并以有限元模型构建的方式分析风力机叶片载荷及强度,最后探究基于Monte Carlo的风力机叶片可靠性。

Fe-Cr二元合金微观组织演化的质量密度场耦合动力学Monte-Carlo模拟研究

本文建立了一种全新的将动力学Monte-Carlo粒子模拟与基于归一化Gauss函数基组的质量密度场空间粗粒化模型耦合的杂化模拟算法.采用该杂化模拟算法,系统对比研究了4种Cr原子含量分别为12.8%,20.0%,30.0%和40.0%的Fe-Cr合金中Cr相在温度为673 K下的时效析出动力学机制,及其时效不同阶段微观组织形貌的演变规律.研究得出Fe-Cr(12.8%)合金富Cr相时效组织形貌呈现孤立颗粒状空间分布形态,时效机制属于形核-长大(NG)机制;对于Fe-Cr(30.0%)和Fe-Cr(40.0%),富Cr相时效形貌在形核-生长及熟化阶段均呈现为三维蠕虫状空间分布特征,时效机制属于条幅分解(SD)机制;对于Fe-Cr(20.0%)合金,其富Cr相组织演化特征介于NG和SD机制之间.研究进一步发现Cr原子短程序参量可用来分析富Cr相形核-生长阶段Fe-Cr合金原子尺度结构的演变,但对于时效熟化阶段微观结构组织变化不敏感.基于空间粗粒化后Fe-Cr合金微观组织形貌,进一步分析了4种Cr原子含量下Fe-Cr合金相变动力学参数如富Cr相体积分数、平均粒径及相颗粒数密度随时效时间演变.本文建立的质量密度场耦合动力学Monte-Carlo模拟方法,为开发多尺度算法模拟合金时效动力学机制及微观组织形貌演变提供了新的思路和研究基础.

Machine learning-enhanced Monte Carlo and subset simulations for advanced risk assessment in transportation infrastructure

The maintenance of safety and dependability in rail and road embankments is of utmost importance in order to facilitate the smooth operation of transportation networks.This study introduces a comprehensive methodology for soil slope stability evaluation,employing Monte Carlo Simulation(MCS)and Subset Simulation(SS)with the"UPSS 3.0 Add-in"in MS-Excel.Focused on an 11.693-meter embankment with a soil slope(inclination ratio of 2H:1V),the investigation considers earthquake coefficients(kh)and pore water pressure ratios(ru)following Indian zoning requirements.The chance of slope failure showed a considerable increase as the Coefficient of Variation(COV),seismic coefficients(kh),and pore water pressure ratios(ru)experienced an escalation.The SS approach showed exceptional efficacy in calculating odds of failure that are notably low.Within computational modeling,the study optimized the worst-case scenario using ANFIS-GA,ANFIS-GWO,ANFIS-PSO,and ANFIS-BBO models.The ANFIS-PSO model exhibits exceptional accuracy(training R2=0.9011,RMSE=0.0549;testing R2=0.8968,RMSE=0.0615),emerging as the most promising.This study highlights the significance of conducting thorough risk assessments and offers practical insights into evaluating and improving the stability of soil slopes in transportation infrastructure.These findings contribute to the enhancement of safety and reliability in real-world situations.

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

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

Chord length sampling correction analysis for dispersion fuel in Monte Carlo simulation

Dispersion fuels,knowned for their excellent safety performance,are widely used in advanced reactors,such as hightemperature gas-cooled reactors.Compared with deterministic methods,the Monte Carlo method has more advantages in the geometric modeling of stochastic media.The explicit modeling method has high computational accuracy and high computational cost.The chord length sampling(CLS)method can improve computational efficiency by sampling the chord length during neutron transport using the matrix chord length?s probability density function.This study shows that the excluded-volume effect in realistic stochastic media can introduce certain deviations into the CLS.A chord length correction approach is proposed to obtain the chord length correction factor by developing the Particle code based on equivalent transmission probability.Through numerical analysis against reference solutions from explicit modeling in the RMC code,it was demonstrated that CLS with the proposed correction method provides good accuracy for addressing the excludedvolume effect in realistic infinite stochastic media.

Monte Carlo method for evaluation of surface emission rate measurement uncertainty

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

基于二次电子发射的Monte Carlo仿真模拟

二次电子发射特性在现代先进真空电子器件中有着广泛的应用,准确测量二次电子的产额数据对进一步确定二次电子发射特性至关重要。然而,目前的实验条件无法保证高真空的实验环境和高精准的测量技术,导致得出的实验数据往往存在着不可避免的误差。因此,利用Monte Carlo计算方法模拟电子束与金属固体的相互作用过程已然成为了一种较为实用的解题方法。此方法是一种基于随机模拟的数值计算方法,通过重复随机采样,根据采样结果计算目标量的期望值,进而得到数值解,它在研究电子运输问题上具有极强的解题能量。模型采用经验Mott散射截面公式来快速描述电子的弹性散射过程,利用介电函数模型来描述电子的非弹性散射过程,从而快速计算出二次电子在金属内部的散射运动轨迹,通过定量分析出二次电子发射产额的影响因素。研究发现,二次电子产额随着入射粒子能量的增加而增加,在一定范围内产额随能量增加的趋势比较明显,但达到一定能量后,产额将趋于饱和。此外,入射角度对二次电子的产额分布也有显著影响。

Health risk assessment of trace metal(loid)s in agricultural soils based on Monte Carlo simulation coupled with positive matrix factorization model in Chongqing, southwest China

This study aimed to investigate the pollution characteristics, source apportionment, and health risks associated with trace metal(loid)s(TMs) in the major agricultural producing areas in Chongqing, China. We analyzed the source apportionment and assessed the health risk of TMs in agricultural soils by using positive matrix factorization(PMF) model and health risk assessment(HRA) model based on Monte Carlo simulation. Meanwhile, we combined PMF and HRA models to explore the health risks of TMs in agricultural soils by different pollution sources to determine the priority control factors. Results showed that the average contents of cadmium(Cd), arsenic (As), lead(Pb), chromium(Cr), copper(Cu), nickel(Ni), and zinc(Zn) in the soil were found to be 0.26, 5.93, 27.14, 61.32, 23.81, 32.45, and 78.65 mg/kg, respectively. Spatial analysis and source apportionment analysis revealed that urban and industrial sources, agricultural sources, and natural sources accounted for 33.0%, 27.7%, and 39.3% of TM accumulation in the soil, respectively. In the HRA model based on Monte Carlo simulation, noncarcinogenic risks were deemed negligible(hazard index <1), the carcinogenic risks were at acceptable level(10^(-6)<total carcinogenic risk ≤ 10^(-4)), with higher risks observed for children compared to adults. The relationship between TMs, their sources, and health risks indicated that urban and industrial sources were primarily associated with As, contributing to 75.1% of carcinogenic risks and 55.7% of non-carcinogenic risks, making them the primary control factors. Meanwhile, agricultural sources were primarily linked to Cd and Pb, contributing to 13.1% of carcinogenic risks and 21.8% of non-carcinogenic risks, designating them as secondary control factors.

Optimization of the Use of Spherical Targets for Point Cloud Registration Using Monte Carlo Simulation

Registrations based on the manual placement of spherical targets are still being employed by many professionals in the industry.However,the placement of those targets usually relies solely on personal experience without scientific evidence supported by numerical analysis.This paper presents a comprehensive investigation,based on Monte Carlo simulation,into determining the optimal number and positions for efficient target placement in typical scenes consisting of a pair of facades.It demonstrates new check-up statistical rules and geometrical constraints that can effectively extract and analyze massive simulations of unregistered point clouds and their corresponding registrations.More than 6×10^(7) sets of the registrations were simulated,whereas more than IOO registrations with real data were used to verify the results of simulation.The results indicated that using five spherical targets is the best choice for the registration of a large typical registration site consisting of two vertical facades and a ground,when there is only a box set of spherical targets available.As a result,the users can avoid placing extra targets to achieve insignificant improvements in registration accuracy.The results also suggest that the higher registration accuracy can be obtained when the ratio between the facade-to-target distance and target-to-scanner distance is approximately 3:2.Therefore,the targets should be placed closer to the scanner rather than in the middle between the facades and the scanner,contradicting to the traditional thought. Besides,the results reveal that the accuracy can be increased by setting the largest projected triangular area of the targets to be large.

Monte Carlo模拟在尾矿土壤铀及伴生重金属污染风险评估中的应用

尾矿及其周边土壤中的铀及重金属污染对生态系统安全产生严重威胁,并对人类健康造成潜在风险。本研究以美国环境保护署提出的健康风险评估模型为基础,借助蒙特卡洛模拟(Monte Carlo Simulation)的确定性和概率性模型及敏感性分析的综合模拟方法,通过对多种暴露途径的综合估算,评估了中国南部某铀尾矿及周边区域土壤对成年男性、成年女性及儿童的致癌和非致癌风险。结果表明,研究区域内存在儿童及成年女性的累积非致癌健康风险和致癌健康风险。敏感性分析结果表明,优先控制砷、铅、铀、锑、钒可有效降低非致癌风险;优先控制砷、铬和镉可有效降低致癌风险。

Enhancing Precision in Radiotherapy Delivery: Validating Monte Carlo Simulation Models for 6 MV Elekta Synergy Agility LINAC Photon Beam Using Two Models of the GAMOS Code

The most crucial requirement in radiation therapy treatment planning is a fast and accurate treatment planning system that minimizes damage to healthy tissues surrounding cancer cells. The use of Monte Carlo toolkits has become indispensable for research aimed at precisely determining the dose in radiotherapy. Among the numerous algorithms developed in recent years, the GAMOS code, which utilizes the Geant4 toolkit for Monte Carlo simula-tions, incorporates various electromagnetic physics models and multiple scattering models for simulating particle interactions with matter. This makes it a valuable tool for dose calculations in medical applications and throughout the patient’s volume. The aim of this present work aims to vali-date the GAMOS code for the simulation of a 6 MV photon-beam output from the Elekta Synergy Agility linear accelerator. The simulation involves mod-eling the major components of the accelerator head and the interactions of the radiation beam with a homogeneous water phantom and particle information was collected following the modeling of the phase space. This space was po-sitioned under the X and Y jaws, utilizing three electromagnetic physics mod-els of the GAMOS code: Standard, Penelope, and Low-Energy, along with three multiple scattering models: Goudsmit-Saunderson, Urban, and Wentzel-VI. The obtained phase space file was used as a particle source to simulate dose distributions (depth-dose and dose profile) for field sizes of 5 × 5 cm2 and 10 × 10 cm2 at depths of 10 cm and 20 cm in a water phantom, with a source-surface distance (SSD) of 90 cm from the target. We compared the three electromagnetic physics models and the three multiple scattering mod-els of the GAMOS code to experimental results. Validation of our results was performed using the gamma index, with an acceptability criterion of 3% for the dose difference (DD) and 3 mm for the distance-to-agreement (DTA). We achieved agreements of 94% and 96%, respectively, between simulation and experimentation for the three electromagnetic physics models and three mul-tiple scattering models, for field sizes of 5 × 5 cm2 and 10 × 10 cm2 for depth-dose curves. For dose profile curves, a good agreement of 100% was found between simulation and experimentation for the three electromagnetic physics models, as well as for the three multiple scattering models for a field size of 5 × 5 cm2 at 10 cm and 20 cm depths. For a field size of 10 × 10 cm2, the Penelope model dominated with 98% for 10 cm, along with the three multiple scattering models. The Penelope model and the Standard model, along with the three multiple scattering models, dominated with 100% for 20 cm. Our study, which compared these different GAMOS code models, can be crucial for enhancing the accuracy and quality of radiotherapy, contributing to more effective patient treatment. Our research compares various electro-magnetic physics models and multiple scattering models with experimental measurements, enabling us to choose the models that produce the most reli-able results, thereby directly impacting the quality of simulations. This en-hances confidence in using these models for treatment planning. Our re-search consistently contributes to the progress of Monte Carlo simulation techniques in radiation therapy, enriching the scientific literature.

Monte Carlo Simulations of Doping Properties of a Spin-3/2 Ising Nanotube

The effect of spin-1 impurities doping on the magnetic properties of a spin-3/2 Ising nanotube is investigated using Monte Carlo simulations within the Blume-Emery-Griffiths model in the presence of an external magnetic field. The thermal behaviors of the order parameters and different macroscopic instabilities as well as the hysteretic behavior of the material are examined in great detail as a function of the dopant density. It is found that the impurities concentration affects all the system magnetic properties generating for some specific values, compensation points and multi-cycle hysteresis. Doping conditions where the saturation/remanent magnetization and coercive field of the investigated material can be modified for permanent or soft magnets synthesis purpose are discussed.

A hybrid method integrating Green's function Monte Carlo and projected entangled pair states

This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function in GFMC.By leveraging PEPS’s proficiency in capturing quantum state entanglement and GFMC’s efficient parallel architecture,the hybrid method is well-suited for the accurate and efficient treatment of frustrated quantum spin systems.As a benchmark,we applied this approach to study the frustrated J_(1)–J_(2) Heisenberg model on a square lattice with periodic boundary conditions(PBCs).Compared with other numerical methods,our approach integrating PEPS and GFMC shows competitive accuracy in the performance of ground-state energy.This paper provides systematic and comprehensive discussion of the approach of our previous work[Phys.Rev.B 109235133(2024)].