A modified Monte Carlo model of speckle tracking of shear wave propagation in scattering media is proposed. The established Monte Carlo model mainly concerns the variations of optical electric field and speckle. The t...A modified Monte Carlo model of speckle tracking of shear wave propagation in scattering media is proposed. The established Monte Carlo model mainly concerns the variations of optical electric field and speckle. The two- dimensional intensity distribution and the time evolution of speckles in different probe locations are obtained. The fluctuation of speckle intensity tracks the acoustic-radiation-force shear wave propagation, and especially the reduction of speckle intensity implies attenuation of shear wave. Then, the shear wave velocity is estimated quantitatively on the basis of the time-to-peak algorithm and linear regression processing. The results reveal that a smaller sampling interval yields higher estimation precision and the shear wave velocity is estimated more efficiently by using speckle intensity difference than by using speckle contrast difference according to the estimation error. Hence, the shear wave velocity is estimated to be 2.25 m/s with relatively high accuracy for the estimation error reaches the minimum （0.071）.展开更多
To sharpen the imaging of structures, it is vital to develop a convenient and efficient quantitative algorithm of the optical coherence tomography (OCT) sampling. In this paper a new Monte Carlo model is set up and ho...To sharpen the imaging of structures, it is vital to develop a convenient and efficient quantitative algorithm of the optical coherence tomography (OCT) sampling. In this paper a new Monte Carlo model is set up and how light propagates in bio-tissue is analyzed in virtue of mathematics and physics equations. The relations,in which light intensity of Class 1 and Class 2 light with different wavelengths changes with their permeation depth,and in which Class 1 light intensity (signal light intensity) changes with the probing depth, and in which angularly resolved diffuse reflectance and diffuse transmittance change with the exiting angle, are studied. The results show that Monte Carlo simulation results are consistent with the theory data.展开更多
To validate the ability of full configuration interaction quantum Monte Carlo （FCIQMC） for studying the 2D Hubbard model near half-filling regime, the ground state energies of a 4×44×4 square lattice syste...To validate the ability of full configuration interaction quantum Monte Carlo （FCIQMC） for studying the 2D Hubbard model near half-filling regime, the ground state energies of a 4×44×4 square lattice system with various interaction strengths are calculated. It is found that the calculated results are in good agreement with those obtained by exact diagonalization （i.e., the exact values for a given basis set） when the population of psi particles （psips） is higher than the critical population required to correctly sample the ground state wave function. In addition, the variations of the average computational time per 20 Monte Carlo cycles with the coupling strength and the number of processors are also analyzed. The calculated results show that the computational efficiency of an FCIQMC calculation is mainly affected by the total population of psips and the communication between processors. These results can provide useful references for understanding the FCIQMC algorithm, studying the ground state properties of the 2D Hubbard model for the larger system size by the FCIQMC method and using a computational budget as effectively as possible.展开更多
A random walk Monte Carlo （RWMC） simulation model of catalytic particle was established on the basis of the structures of bismuth molybdate catalysts and mechanisms of catalytic reactions with propylene selective ox...A random walk Monte Carlo （RWMC） simulation model of catalytic particle was established on the basis of the structures of bismuth molybdate catalysts and mechanisms of catalytic reactions with propylene selective oxidation and ammoxidation. The simulation results show that rationality of the RWMC model is proved on the basis of pulse experimental data. One of the most remarkable factors affecting catalytic behavior is the transfer of bulk lattice oxygen, which decides the rate of ammonia-consuming and propylene-consuming. The selectivity of main products reaches the maximum after the reduction of catalysts to a certain degree. It is inferred that catalytic performance improves greatly if the ratio of capacity for dehydrogenation from adsorbed propylene molecule on catalytically active site of molybdenum metal-imido group （Mo=NH） to that on catalytically active site of molybdenum metal-oxo group （Mo=O） becomes much higher.展开更多
Numerical simulations by means of the Monte Carlo Potts model have been provided to simulate grain structures in two-phase polycrystalline materials. The topological features in the simulated microstructure analyzed f...Numerical simulations by means of the Monte Carlo Potts model have been provided to simulate grain structures in two-phase polycrystalline materials. The topological features in the simulated microstructure analyzed for different diffusion mechanisms over a broad range of volume fractions for both phases. The topological properties include the average number of sides, grain topology distribution </span><span style="font-family:Verdana;">and</span><span style="font-family:Verdana;"> the topological size relation function. It is found that the average number of sides depends proportionally on the volume fraction. It increases as the </span><span style="font-family:Verdana;">volumes</span><span style="font-family:Verdana;"> fraction increases and vice versa. Moreover, it is shown that the grain topology distribution in the self-similar growth regime can be described by </span><span style="font-family:Verdana;">time</span><span style="font-family:Verdana;"> unchanged function of the relative grain size. Additionally, topological size function in the simulated microstructure can be evaluated by a quadratic function.展开更多
The magnetic moment of 2＋1 state for 10Be are calculated and investigated in terms of single particle orbits for protons and neutrons under the framework of ab initio Monte Carlo shell model method in an emax=3 model...The magnetic moment of 2＋1 state for 10Be are calculated and investigated in terms of single particle orbits for protons and neutrons under the framework of ab initio Monte Carlo shell model method in an emax=3 model space. The reduced matrix elements of orbital and spin angular momentum are evaluated. It is found that the orientations of orbital angular momentum in different single particle orbits are consistent. Conversely, the orientations of spin in different single particle orbits tend to be chaotic. The nuclear magnetic moment of 2＋1 state for 10Be is obtained as 1.006 ,UN and is discussed in regards to the contribution of orbital and spin angular momentum both for protons and neutrons. The corresponding g-factor is also given.展开更多
The coupling between the Monte Carlo （MC） method and geometrical optics to improve accuracy is investigated. The results obtained show improved agreement with previous experimental data, demonstrating that the MC me...The coupling between the Monte Carlo （MC） method and geometrical optics to improve accuracy is investigated. The results obtained show improved agreement with previous experimental data, demonstrating that the MC method, when coupled with simple geometrical optics, can simulate multiple scattering with enhanced fidelity.展开更多
In this study,a comprehensive analysis of jets and underlying events as a function of charged particle multiplicity in proton-proton(pp)collisions at a center-of-mass energy of √s=7 TeV is conducted.Various Monte Car...In this study,a comprehensive analysis of jets and underlying events as a function of charged particle multiplicity in proton-proton(pp)collisions at a center-of-mass energy of √s=7 TeV is conducted.Various Monte Carlo(MC)event generators,including Pythia8.308,EPOS 1.99,EPOSLHC,EPOS4_(Hydro),and EPOS4_(noHydro),are employed to predict particle production.The predictions from these models are compared with experimental data from the CMS collaboration.The charged particles are categorized into those associated with underlying events and those linked to jets,and the analysis is restricted to charged particles with|η|<2.4 and p_T>0.25 GeV/c.By comparing the MC predictions with CMS data,we find that EPOS4_(Hydro),EPOSLHC,and Pythia8 consistently reproduce the experimental results for all charged particles,underlying events,intrajets,and leading charged particles.For charged jet rates with p_T^(ch.jet)>5 GeV/c,EPOS4_(Hydro)and Pythia8 perform exceptionally well.In the case of charged jet rates with p_T^(ch.jet)→30 GeV/c,EPOSLHC reproduces satisfactorily good results,whereas EPOS4 Hydro exhibits good agreement with the data at higher charged particle multiplicities compared to the other models.This can be attributed to the conversion of energy into flow when"Hydro=on"leading to an increase in multiplicity.The EPOSLHC model describes the data better owing to the new collective flow effects,correlated flow treatment,and parameterization compared to EPOS 1.99.However,the examination of the jet p_T spectrum and normalized charged p_T density reveals that EPOS4_(Hydro),EPOS4_(noHydro),and EPOSLHC exhibit good agreement with the experimental results,whereas Pythia8 and EPOS 1.99 do not perform as well owing to the lack of correlated flow treatment.展开更多
The unitary correlation operator method （UCOM） and the similarity renormalization group theory （SRG） are compared and discussed in the framework of the no-core Monte Carlo shell model （MCSM） calculations for ^3H...The unitary correlation operator method （UCOM） and the similarity renormalization group theory （SRG） are compared and discussed in the framework of the no-core Monte Carlo shell model （MCSM） calculations for ^3H and ^4He. The treatment of spurious center-of-mass motion by Lawson＇s prescription is performed in the MCSM calculations. These results with both transformed interactions show good suppression of spurious center-of-mass motion with proper Lawson＇s prescription parameter βc.m. values. The UCOM potentials obtain faster convergence of total energy for the ground state than that of SRG potentials in the MCSM calculations, which differs from the cases in the no-core shell model calculations （NCSM）. These differences are discussed and analyzed in terms of the truncation scheme in the MCSM and NCSM, as well as the properties of the potentials of SRG and UCOM.展开更多
Geometrical optics and the Monte Carlo method are very flexible in dealing with the interaction of light with non-spherical particles, but usually diffraction is not considered. To cover this gap, the Heisenberg Uncer...Geometrical optics and the Monte Carlo method are very flexible in dealing with the interaction of light with non-spherical particles, but usually diffraction is not considered. To cover this gap, the Heisenberg Uncertainty Monte Carlo （HUMC） model is applied to calculate separately the diffraction of a ray or a photon. In this paper, we report an improvement of the HUMC model by specifying the phase of the photon subject to the Fraunhofer diffraction condition. After validating the model by comparing its results with analytical results for apertures of simple shapes, the HUMC model is then applied in simulations of Fraunhofer diffraction by apertures of complex shapes, such as those composed of one or two elliptical openings. We have shown that the diffracted intensity distributions of simple apertures obtained by the HUMC model are in good agreement with the results calculated from analytical expressions. The simulations of diffraction by apertures composed of two square or elliptical openings prove that the HUMC model is a powerful and flexible too] for predicting the Fraunhofer diffraction by a complex optical system.展开更多
A study has just been carried out on hot electron effects in GaAs/Al0.3Ga0.7As potential well barrier(PWB)diodes using both Monte Carlo(MC)and drift-diffusion(DD)models of charge transport.We show the operation and be...A study has just been carried out on hot electron effects in GaAs/Al0.3Ga0.7As potential well barrier(PWB)diodes using both Monte Carlo(MC)and drift-diffusion(DD)models of charge transport.We show the operation and behaviour of the diode in terms of electric field,mean electron velocity and potential,mean energy of electrons andΓ-valley population.The MC model predicts lower currents flowing through the diode due to back scattering at anode(collector)and carrier heating at higher bias.At a bias of 1.0 V,the current density obtained from experimental result,MC and DD simulation models are 1.35,1.12 and1.77μA/μm^2 respectively.The reduction in current over conventional model,is compensated to a certain extent because less charge settles in the potential well and so the barrier is slightly reduced.The DD model results in higher currents under the same bias and conditions.However,at very low bias specifically,up to 0.3 V without any carrier heating effects,the DD and MC models look pretty similar as experimental results.The significant differences observed in the I-V characteristics of the DD and MC models at higher biases confirm the importance of energy transport when considering these devices.展开更多
The Monte Carlo probability(MCP)model,which has been used for official tropical cyclone(TC)warnings to the public by the United States’National Hurricane Center(NHC),can estimate the probability of wind speed in the ...The Monte Carlo probability(MCP)model,which has been used for official tropical cyclone(TC)warnings to the public by the United States’National Hurricane Center(NHC),can estimate the probability of wind speed in the vicinity of a TC during the forecast period.It has been successful in the operational environment for many years.However,due to its strong dependence on a given forecast track(e.g.,forecast from the NCEP Global Forecast System),the MCP model may generate a poor probability map for TCs near landfall.In this study,we proposed and tested a modified MCP method for TC forecasts near landfall.We first adjusted the MCP model by adding limits to the direction angle and motion distance to deal with the substantial change in TC moving direction and the low wind speeds during landfall.Then,we combined ensemble probability maps generated from ECMWF,United Kingdom Met Office(UKMO),and NCEP ensemble forecasts,obtained from The International Grand Global Ensemble(TIGGE),into the MCP model to configure a modified MCP model.Wind speed probability maps for the 0-120-h forecast from both the original and modified MCP models are compared.It is found that the modified MCP model can provide a better wind speed probability map during landfall,especially at wind speeds of 20-64 kt near TC landfall.The results from this study prove the benefits of combining the MCP model with ensemble forecasting in potential applications for improved TC forecasts.展开更多
The grain size in the heat-affected zone （ HAZ） of workpieces is an important factor for evaluating the welding quality. The Monte Carlo technique of grain growth in the heat-affected zone of alloy is widely used. A...The grain size in the heat-affected zone （ HAZ） of workpieces is an important factor for evaluating the welding quality. The Monte Carlo technique of grain growth in the heat-affected zone of alloy is widely used. A finite element model （FEM） was employed to simulate the temperature field of gas metal arc welding （GMAW） ; the data of thermal cycles were used in the Monte Carlo model to calculate the grain sizes in different welding heat inputs. The equation for the relationship between Monte Carlo Step and real time has been estimated using regression analysis. Then mathematics model of the grain growth could be worked out by both the kinetic model and the Monte Carlo model. The dynamic process of grain growth was simulated by the result of the Monte Carlo model. The experimental result was used to prove the validity of this method in simulating of microstructure.展开更多
The shift in the percolation threshold of compressed composites was studied by a 3D continuum percolation model. A Monte Carlo (MC) method was employed in the simulations. The percolation threshold was found to rise w...The shift in the percolation threshold of compressed composites was studied by a 3D continuum percolation model. A Monte Carlo (MC) method was employed in the simulations. The percolation threshold was found to rise with the compression strain, which captures the basic trend in compression-induced conductivity variation from the experiments. Both fiber bending and texture formation contribute to the percolation threshold. The results suggest that fillers with a high aspect ratio are more desirable for sensor and electrical switch applications.展开更多
The longitudinal tensile properties of SiCf/Ti-6Al-4V composites with different fiber volume fractions were simulated by the Monte Carlo 2-D finite element model. The random distribution of fiber strength was expresse...The longitudinal tensile properties of SiCf/Ti-6Al-4V composites with different fiber volume fractions were simulated by the Monte Carlo 2-D finite element model. The random distribution of fiber strength was expressed by the two-parameter Weibull function. Meanwhile, contact elements and birth-death elements were used to describe the interfacial sliding process after debonding and fiber breakage（or matrix cracking） respectively, which was realized by subroutine complied in ANSYS-APDL（ANSYS Parametric Design Language）. The experimental results show that the yield stress and ultimate tensile strength of SiCf/Ti-6Al-4V composites increase with increasing fiber volume fraction, while the corresponding strain of them is just on the contrary. In addition, almost the same failure mode is obtained in SiCf/Ti-6Al-4V composites with various fiber volume fractions when the interfacial shear strength is fixed. Finally, the tensile strength predicted by finite element analysis is compared with that predicted by Global load-sharing model, Local load-sharing model and conventional rule of mixtures, thus drawing the conclusion that Local load-sharing model is very perfect for the prediction of the ultimate tensile strength.展开更多
Aiming at the assembly accuracy of a large aircraft transport jig, the effect of component error and the error of work-piece surface on the work-piece position and orientation in the 3-2-1 fixturing scheme is studied ...Aiming at the assembly accuracy of a large aircraft transport jig, the effect of component error and the error of work-piece surface on the work-piece position and orientation in the 3-2-1 fixturing scheme is studied with the object pose space description method. The error mapping model between the connecting part of the front frame rack and its support base is modeled using the homogeneous transformation matrix(HTM) method. The probabilistic error is simulated using the Monte Carlo method. The measurement experiment was conducted by the laser tracker to verify the effectiveness of the approach, and the approach has been successfully applied to the production of transport jig.展开更多
This manuscript presents a stochastic model updating method, taking both uncertainties in models and variability in testing into account. The updated finite element（FE） models obtained through the proposed technique...This manuscript presents a stochastic model updating method, taking both uncertainties in models and variability in testing into account. The updated finite element（FE） models obtained through the proposed technique can aid in the analysis and design of structural systems. The authors developed a stochastic model updating method integrating distance discrimination analysis（DDA） and advanced Monte Carlo（MC） technique to（1） enable more efficient MC by using a response surface model,（2） calibrate parameters with an iterative test-analysis correlation based upon DDA, and（3） utilize and compare different distance functions as correlation metrics. Using DDA, the influence of distance functions on model updating results is analyzed. The proposed stochastic method makes it possible to obtain a precise model updating outcome with acceptable calculation cost. The stochastic method is demonstrated on a helicopter case study updated using both Euclidian and Mahalanobis distance metrics. It is observed that the selected distance function influences the iterative calibration process and thus, the calibration outcome, indicating that an integration of different metrics might yield improved results.展开更多
We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we int...We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we introduced multibody-interaction parameters that establish correlations between the attaching and detaching rate constants and the resulting thickness and width of the crystalline lamella.Using MATLAB and Monte Carlo method,we followed the evolution of the secondary nuclei as a function of various multibody-interaction parameters.We identified three different growth progressions of the crystal：（i） Widening,（ii） thickening and（iii） simultaneously thickening and widening of lamellar crystals,controlled by the corresponding kinetic parameters.展开更多
In recent years, private sectors are encouraged to take an active part in franchising of urban infrastructure investments and operations, which promotes the rapid development of public-private partnership(PPP)in infra...In recent years, private sectors are encouraged to take an active part in franchising of urban infrastructure investments and operations, which promotes the rapid development of public-private partnership(PPP)in infrastructure and public service supply. Value for money(VFM) assessment has been officially proposed to provide a reference for selection of projects planning to adopt PPP. Based on the bottlenecks of VFM application in China and the uncertainties for urban infrastructure PPP projects, a discounted cash flow(DCF) model is established for VFM of infrastructure PPP projects. Then, a Monte Carlo simulation model is established on the basis of uncertainty factors for VFM. Through the analysis of Huai'an trams PPP project, coping strategies of uncertainties for VFM are put forward. Findings of the research may propel the establishment of a complete VFM evaluation system for PPP projects. Key instructional functions of VFM during the process of decision-making can be brought into full play and PPP may develop orderly.展开更多
基金Supported by the National Key Scientific Instrument and Equipment Development Projects of China under Grant No 81127901the National Natural Science Foundation of China under Grant Nos 61372017 and 30970828
文摘A modified Monte Carlo model of speckle tracking of shear wave propagation in scattering media is proposed. The established Monte Carlo model mainly concerns the variations of optical electric field and speckle. The two- dimensional intensity distribution and the time evolution of speckles in different probe locations are obtained. The fluctuation of speckle intensity tracks the acoustic-radiation-force shear wave propagation, and especially the reduction of speckle intensity implies attenuation of shear wave. Then, the shear wave velocity is estimated quantitatively on the basis of the time-to-peak algorithm and linear regression processing. The results reveal that a smaller sampling interval yields higher estimation precision and the shear wave velocity is estimated more efficiently by using speckle intensity difference than by using speckle contrast difference according to the estimation error. Hence, the shear wave velocity is estimated to be 2.25 m/s with relatively high accuracy for the estimation error reaches the minimum （0.071）.
文摘To sharpen the imaging of structures, it is vital to develop a convenient and efficient quantitative algorithm of the optical coherence tomography (OCT) sampling. In this paper a new Monte Carlo model is set up and how light propagates in bio-tissue is analyzed in virtue of mathematics and physics equations. The relations,in which light intensity of Class 1 and Class 2 light with different wavelengths changes with their permeation depth,and in which Class 1 light intensity (signal light intensity) changes with the probing depth, and in which angularly resolved diffuse reflectance and diffuse transmittance change with the exiting angle, are studied. The results show that Monte Carlo simulation results are consistent with the theory data.
基金Supported by the Natural Science Foundation for Colleges and Universities of Jiangsu Province under Grant No 16KJB140008the National Natural Science Foundation of China under Grant Nos 11447204 and 11647164+1 种基金the Natural Science Foundation of Jiangsu Province under Grant No BK20151079the Scientific Research Foundation of Nanjing Xiaozhuang University under Grant No 2015NXY34
文摘To validate the ability of full configuration interaction quantum Monte Carlo （FCIQMC） for studying the 2D Hubbard model near half-filling regime, the ground state energies of a 4×44×4 square lattice system with various interaction strengths are calculated. It is found that the calculated results are in good agreement with those obtained by exact diagonalization （i.e., the exact values for a given basis set） when the population of psi particles （psips） is higher than the critical population required to correctly sample the ground state wave function. In addition, the variations of the average computational time per 20 Monte Carlo cycles with the coupling strength and the number of processors are also analyzed. The calculated results show that the computational efficiency of an FCIQMC calculation is mainly affected by the total population of psips and the communication between processors. These results can provide useful references for understanding the FCIQMC algorithm, studying the ground state properties of the 2D Hubbard model for the larger system size by the FCIQMC method and using a computational budget as effectively as possible.
基金国家自然科学基金，the Fundamental Research Foundation of SINOPEC
文摘A random walk Monte Carlo （RWMC） simulation model of catalytic particle was established on the basis of the structures of bismuth molybdate catalysts and mechanisms of catalytic reactions with propylene selective oxidation and ammoxidation. The simulation results show that rationality of the RWMC model is proved on the basis of pulse experimental data. One of the most remarkable factors affecting catalytic behavior is the transfer of bulk lattice oxygen, which decides the rate of ammonia-consuming and propylene-consuming. The selectivity of main products reaches the maximum after the reduction of catalysts to a certain degree. It is inferred that catalytic performance improves greatly if the ratio of capacity for dehydrogenation from adsorbed propylene molecule on catalytically active site of molybdenum metal-imido group （Mo=NH） to that on catalytically active site of molybdenum metal-oxo group （Mo=O） becomes much higher.
文摘Numerical simulations by means of the Monte Carlo Potts model have been provided to simulate grain structures in two-phase polycrystalline materials. The topological features in the simulated microstructure analyzed for different diffusion mechanisms over a broad range of volume fractions for both phases. The topological properties include the average number of sides, grain topology distribution </span><span style="font-family:Verdana;">and</span><span style="font-family:Verdana;"> the topological size relation function. It is found that the average number of sides depends proportionally on the volume fraction. It increases as the </span><span style="font-family:Verdana;">volumes</span><span style="font-family:Verdana;"> fraction increases and vice versa. Moreover, it is shown that the grain topology distribution in the self-similar growth regime can be described by </span><span style="font-family:Verdana;">time</span><span style="font-family:Verdana;"> unchanged function of the relative grain size. Additionally, topological size function in the simulated microstructure can be evaluated by a quadratic function.
基金supported by the Fundamental Research Funds for the Central Universities (JUSRP1035)NSFC (Grant Nos. 11205068 and 11305077)CPSC (Grant No. 2012M520667)
文摘The magnetic moment of 2＋1 state for 10Be are calculated and investigated in terms of single particle orbits for protons and neutrons under the framework of ab initio Monte Carlo shell model method in an emax=3 model space. The reduced matrix elements of orbital and spin angular momentum are evaluated. It is found that the orientations of orbital angular momentum in different single particle orbits are consistent. Conversely, the orientations of spin in different single particle orbits tend to be chaotic. The nuclear magnetic moment of 2＋1 state for 10Be is obtained as 1.006 ,UN and is discussed in regards to the contribution of orbital and spin angular momentum both for protons and neutrons. The corresponding g-factor is also given.
文摘The coupling between the Monte Carlo （MC） method and geometrical optics to improve accuracy is investigated. The results obtained show improved agreement with previous experimental data, demonstrating that the MC method, when coupled with simple geometrical optics, can simulate multiple scattering with enhanced fidelity.
基金Supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project number(PNURSP2024R106)Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia。
文摘In this study,a comprehensive analysis of jets and underlying events as a function of charged particle multiplicity in proton-proton(pp)collisions at a center-of-mass energy of √s=7 TeV is conducted.Various Monte Carlo(MC)event generators,including Pythia8.308,EPOS 1.99,EPOSLHC,EPOS4_(Hydro),and EPOS4_(noHydro),are employed to predict particle production.The predictions from these models are compared with experimental data from the CMS collaboration.The charged particles are categorized into those associated with underlying events and those linked to jets,and the analysis is restricted to charged particles with|η|<2.4 and p_T>0.25 GeV/c.By comparing the MC predictions with CMS data,we find that EPOS4_(Hydro),EPOSLHC,and Pythia8 consistently reproduce the experimental results for all charged particles,underlying events,intrajets,and leading charged particles.For charged jet rates with p_T^(ch.jet)>5 GeV/c,EPOS4_(Hydro)and Pythia8 perform exceptionally well.In the case of charged jet rates with p_T^(ch.jet)→30 GeV/c,EPOSLHC reproduces satisfactorily good results,whereas EPOS4 Hydro exhibits good agreement with the data at higher charged particle multiplicities compared to the other models.This can be attributed to the conversion of energy into flow when"Hydro=on"leading to an increase in multiplicity.The EPOSLHC model describes the data better owing to the new collective flow effects,correlated flow treatment,and parameterization compared to EPOS 1.99.However,the examination of the jet p_T spectrum and normalized charged p_T density reveals that EPOS4_(Hydro),EPOS4_(noHydro),and EPOSLHC exhibit good agreement with the experimental results,whereas Pythia8 and EPOS 1.99 do not perform as well owing to the lack of correlated flow treatment.
基金Supported by Fundamental Research Funds for the Central Universities(JUSRP1035)National Natural Science Foundation of China(11305077)
文摘The unitary correlation operator method （UCOM） and the similarity renormalization group theory （SRG） are compared and discussed in the framework of the no-core Monte Carlo shell model （MCSM） calculations for ^3H and ^4He. The treatment of spurious center-of-mass motion by Lawson＇s prescription is performed in the MCSM calculations. These results with both transformed interactions show good suppression of spurious center-of-mass motion with proper Lawson＇s prescription parameter βc.m. values. The UCOM potentials obtain faster convergence of total energy for the ground state than that of SRG potentials in the MCSM calculations, which differs from the cases in the no-core shell model calculations （NCSM）. These differences are discussed and analyzed in terms of the truncation scheme in the MCSM and NCSM, as well as the properties of the potentials of SRG and UCOM.
文摘Geometrical optics and the Monte Carlo method are very flexible in dealing with the interaction of light with non-spherical particles, but usually diffraction is not considered. To cover this gap, the Heisenberg Uncertainty Monte Carlo （HUMC） model is applied to calculate separately the diffraction of a ray or a photon. In this paper, we report an improvement of the HUMC model by specifying the phase of the photon subject to the Fraunhofer diffraction condition. After validating the model by comparing its results with analytical results for apertures of simple shapes, the HUMC model is then applied in simulations of Fraunhofer diffraction by apertures of complex shapes, such as those composed of one or two elliptical openings. We have shown that the diffracted intensity distributions of simple apertures obtained by the HUMC model are in good agreement with the results calculated from analytical expressions. The simulations of diffraction by apertures composed of two square or elliptical openings prove that the HUMC model is a powerful and flexible too] for predicting the Fraunhofer diffraction by a complex optical system.
基金supported by the College of Physical Sciences, University of Aberdeen, UK
文摘A study has just been carried out on hot electron effects in GaAs/Al0.3Ga0.7As potential well barrier(PWB)diodes using both Monte Carlo(MC)and drift-diffusion(DD)models of charge transport.We show the operation and behaviour of the diode in terms of electric field,mean electron velocity and potential,mean energy of electrons andΓ-valley population.The MC model predicts lower currents flowing through the diode due to back scattering at anode(collector)and carrier heating at higher bias.At a bias of 1.0 V,the current density obtained from experimental result,MC and DD simulation models are 1.35,1.12 and1.77μA/μm^2 respectively.The reduction in current over conventional model,is compensated to a certain extent because less charge settles in the potential well and so the barrier is slightly reduced.The DD model results in higher currents under the same bias and conditions.However,at very low bias specifically,up to 0.3 V without any carrier heating effects,the DD and MC models look pretty similar as experimental results.The significant differences observed in the I-V characteristics of the DD and MC models at higher biases confirm the importance of energy transport when considering these devices.
基金supported by the US National Science Foundation(ECCS-1839833 and OAC-2004658)。
文摘The Monte Carlo probability(MCP)model,which has been used for official tropical cyclone(TC)warnings to the public by the United States’National Hurricane Center(NHC),can estimate the probability of wind speed in the vicinity of a TC during the forecast period.It has been successful in the operational environment for many years.However,due to its strong dependence on a given forecast track(e.g.,forecast from the NCEP Global Forecast System),the MCP model may generate a poor probability map for TCs near landfall.In this study,we proposed and tested a modified MCP method for TC forecasts near landfall.We first adjusted the MCP model by adding limits to the direction angle and motion distance to deal with the substantial change in TC moving direction and the low wind speeds during landfall.Then,we combined ensemble probability maps generated from ECMWF,United Kingdom Met Office(UKMO),and NCEP ensemble forecasts,obtained from The International Grand Global Ensemble(TIGGE),into the MCP model to configure a modified MCP model.Wind speed probability maps for the 0-120-h forecast from both the original and modified MCP models are compared.It is found that the modified MCP model can provide a better wind speed probability map during landfall,especially at wind speeds of 20-64 kt near TC landfall.The results from this study prove the benefits of combining the MCP model with ensemble forecasting in potential applications for improved TC forecasts.
文摘The grain size in the heat-affected zone （ HAZ） of workpieces is an important factor for evaluating the welding quality. The Monte Carlo technique of grain growth in the heat-affected zone of alloy is widely used. A finite element model （FEM） was employed to simulate the temperature field of gas metal arc welding （GMAW） ; the data of thermal cycles were used in the Monte Carlo model to calculate the grain sizes in different welding heat inputs. The equation for the relationship between Monte Carlo Step and real time has been estimated using regression analysis. Then mathematics model of the grain growth could be worked out by both the kinetic model and the Monte Carlo model. The dynamic process of grain growth was simulated by the result of the Monte Carlo model. The experimental result was used to prove the validity of this method in simulating of microstructure.
基金Project supported by the National Natural Science Foundation of China (No 10832009)the National Basic Research Program (973) of China (No 2004CB619304)the Science Foundation of Chinese University (No 2009QNA4034)
文摘The shift in the percolation threshold of compressed composites was studied by a 3D continuum percolation model. A Monte Carlo (MC) method was employed in the simulations. The percolation threshold was found to rise with the compression strain, which captures the basic trend in compression-induced conductivity variation from the experiments. Both fiber bending and texture formation contribute to the percolation threshold. The results suggest that fillers with a high aspect ratio are more desirable for sensor and electrical switch applications.
基金Funded by the National Natural Science Foundation of China(51271147)
文摘The longitudinal tensile properties of SiCf/Ti-6Al-4V composites with different fiber volume fractions were simulated by the Monte Carlo 2-D finite element model. The random distribution of fiber strength was expressed by the two-parameter Weibull function. Meanwhile, contact elements and birth-death elements were used to describe the interfacial sliding process after debonding and fiber breakage（or matrix cracking） respectively, which was realized by subroutine complied in ANSYS-APDL（ANSYS Parametric Design Language）. The experimental results show that the yield stress and ultimate tensile strength of SiCf/Ti-6Al-4V composites increase with increasing fiber volume fraction, while the corresponding strain of them is just on the contrary. In addition, almost the same failure mode is obtained in SiCf/Ti-6Al-4V composites with various fiber volume fractions when the interfacial shear strength is fixed. Finally, the tensile strength predicted by finite element analysis is compared with that predicted by Global load-sharing model, Local load-sharing model and conventional rule of mixtures, thus drawing the conclusion that Local load-sharing model is very perfect for the prediction of the ultimate tensile strength.
基金Supported by National Key Technology Research and Development Program of China(No.2012BAF01B07)
文摘Aiming at the assembly accuracy of a large aircraft transport jig, the effect of component error and the error of work-piece surface on the work-piece position and orientation in the 3-2-1 fixturing scheme is studied with the object pose space description method. The error mapping model between the connecting part of the front frame rack and its support base is modeled using the homogeneous transformation matrix(HTM) method. The probabilistic error is simulated using the Monte Carlo method. The measurement experiment was conducted by the laser tracker to verify the effectiveness of the approach, and the approach has been successfully applied to the production of transport jig.
基金supported by the National Natural Science Foundation of China (No. 10972019)the Innovation Foundation of BUAA for Ph.D. Graduates of China, and the China Scholarship Council
文摘This manuscript presents a stochastic model updating method, taking both uncertainties in models and variability in testing into account. The updated finite element（FE） models obtained through the proposed technique can aid in the analysis and design of structural systems. The authors developed a stochastic model updating method integrating distance discrimination analysis（DDA） and advanced Monte Carlo（MC） technique to（1） enable more efficient MC by using a response surface model,（2） calibrate parameters with an iterative test-analysis correlation based upon DDA, and（3） utilize and compare different distance functions as correlation metrics. Using DDA, the influence of distance functions on model updating results is analyzed. The proposed stochastic method makes it possible to obtain a precise model updating outcome with acceptable calculation cost. The stochastic method is demonstrated on a helicopter case study updated using both Euclidian and Mahalanobis distance metrics. It is observed that the selected distance function influences the iterative calibration process and thus, the calibration outcome, indicating that an integration of different metrics might yield improved results.
基金financially supported by the National Natural Science Foundation of China(No.21374054)the Sino-German Center for Research Promotion
文摘We present simulations of the mechanism of secondary nucleation of polymer crystallization,based on a new model accounting for the microscopic kinetics of attaching and detaching.As the key feature of the model,we introduced multibody-interaction parameters that establish correlations between the attaching and detaching rate constants and the resulting thickness and width of the crystalline lamella.Using MATLAB and Monte Carlo method,we followed the evolution of the secondary nuclei as a function of various multibody-interaction parameters.We identified three different growth progressions of the crystal：（i） Widening,（ii） thickening and（iii） simultaneously thickening and widening of lamellar crystals,controlled by the corresponding kinetic parameters.
基金the Research Project of China Institute of Urban Governance of SJTU(No.16JCCS09)
文摘In recent years, private sectors are encouraged to take an active part in franchising of urban infrastructure investments and operations, which promotes the rapid development of public-private partnership(PPP)in infrastructure and public service supply. Value for money(VFM) assessment has been officially proposed to provide a reference for selection of projects planning to adopt PPP. Based on the bottlenecks of VFM application in China and the uncertainties for urban infrastructure PPP projects, a discounted cash flow(DCF) model is established for VFM of infrastructure PPP projects. Then, a Monte Carlo simulation model is established on the basis of uncertainty factors for VFM. Through the analysis of Huai'an trams PPP project, coping strategies of uncertainties for VFM are put forward. Findings of the research may propel the establishment of a complete VFM evaluation system for PPP projects. Key instructional functions of VFM during the process of decision-making can be brought into full play and PPP may develop orderly.
