In public health,simulation modeling stands as an invaluable asset,enabling the evaluation of new systems without their physical implementation,experimentation with existing systems without operational adjustments,and...In public health,simulation modeling stands as an invaluable asset,enabling the evaluation of new systems without their physical implementation,experimentation with existing systems without operational adjustments,and testing system limits without real-world repercussions.In simulation modeling,the Monte Carlo method emerges as a powerful yet underutilized tool.Although the Monte Carlo method has not yet gained widespread prominence in healthcare,its technological capabilities hold promise for substantial cost reduction and risk mitigation.In this review article,we aimed to explore the transformative potential of the Monte Carlo method in healthcare contexts.We underscore the significance of experiential insights derived from simulated experimentation,especially in resource-constrained scenarios where time,financial constraints,and limited resources necessitate innovative and efficient approaches.As public health faces increasing challenges,incorporating the Monte Carlo method presents an opportunity for enhanced system construction,analysis,and evaluation.展开更多
Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random error...Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system's dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error(DTE) fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system's dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.展开更多
This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced ...This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.展开更多
The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detect...The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detector distance in a high-purity germanium detector system was studied using a Geant4 Monte Carlo simulation.The effect of the deadlayer in the germanium crystal was estimated by reference to experiments on 241 Am and the relative efficiency of the detector.The peak efficiency dependence on the source-to-detector distance was compared between the simulation and measurements.The modified sum-peak method is discussed with respect to these peak efficiencies.展开更多
Considering the ocean water's optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable (CV) quantum key distribution (QKD) based Mr-to-water channel. The effects of...Considering the ocean water's optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable (CV) quantum key distribution (QKD) based Mr-to-water channel. The effects of the absorp- tion and scattering on the transmittance of underwater quantum channel and the maximum secure transmission distance are studied. Considering the roughness of the sea surface, we simulate the performance bounds of CV QKD with different wind speeds using the Monte Carlo method. The results show that even if the secret key rate gradually reduces as the wind speed increases, the maximum transmission distance will not be affected obviously. Compared to the works regarding short-distance underwater optical communication, our research represents a significant step towards establishing secure communication between air platform and submarine vehicle.展开更多
The structure and behaviour of LiF-KF solution,as a typical common-anion system,has been simulated by Monte Carlo method.The calculation of partial radial distribution function of ions,heat of mixing and potential ene...The structure and behaviour of LiF-KF solution,as a typical common-anion system,has been simulated by Monte Carlo method.The calculation of partial radial distribution function of ions,heat of mixing and potential energy distribution shows that the average distance be- tween Li^+ and F^- ions will significantly narrow after mixing of molten LiF and KF.This is very similar to the lean-on-one-side effect in molten LiF-KCl solution.The calculated heat of mixing is in fair agreement with the measured one.The dominant source of the energy of mixing may be that the decrease of the repulsion energy between cations,the decrease of the attraction energy between cations and anions,and the decrease of the repulsion energy be- tween anions.展开更多
Pricing the American put options requires solving an optimal stopping problem and therefore is a challenge for the setting up of simulation parameters.This paper uses least square Monte Carlo(LSMC) simulation to pri...Pricing the American put options requires solving an optimal stopping problem and therefore is a challenge for the setting up of simulation parameters.This paper uses least square Monte Carlo(LSMC) simulation to price the American put options and output the optimal simulation steps and number of Hermite basis functions.The results suggest:with different time cost and error tolerance,investors can choose the optimal simulation steps and number of basis function individually to price American put options numerically.Generally,with the pre-limitation in the section "least square Monte Carlo simulation",a number of basis equals 4,15 000 simulation steps for Hermite basis function appear to be sufficient for the method.展开更多
Artificial earthquake catalogue simulation is one of the ways to effectively improve the incompleteness of the existing earthquake catalogue,the scarcity of large earthquake records and the improvement of seismologica...Artificial earthquake catalogue simulation is one of the ways to effectively improve the incompleteness of the existing earthquake catalogue,the scarcity of large earthquake records and the improvement of seismological research. Based on the Poisson distribution model of seismic activity and the Gutenberg-Richter magnitude-frequency relationship,the Monte Carlo method which can describe the characteristics of the stochastic nature and the physical experiment process is used. This paper simulates the future seismic catalogues of the Fenhe-Weihe seismic belt of different durations and conducts statistical tests on them.The analysis shows that the simulation catalogue meets the set seismic activity parameters and meets the Poisson distribution hypothesis,which can obtain a better simulated earthquake catalogues that meets the seismic activity characteristics. According to the simulated earthquake catalogues,future earthquake trends in this region are analyzed to provide reference for seismic hazard analysis.展开更多
Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key ...Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key to measure the floating reference point(FRP)precisely for the application of FRT.Monte Carlo(MC)simulation has been introduced to quantitatively in-vestigate the effects of positioning errors and light source drifts on measuring FRP.In this article,thinning and calculating method(TCM)is proposed to quantify the positioning error.Mean-while,the normalization process(NP)is developed to significantly reduce the error induced by light source drift.The results according to TCM show that 7 purm deviations in positioning can generate about 10.63%relative error in FRP.It is more noticeable that 1%fluctuation in light source intensity may lead to 12.21%relative errors.Gratifyingly,the proposed NP model can effectively reduce the error caused by light source drift.Therefore,the measurement system for FRPs must meet that the positioning error is less than 7 purm,and the light source drift is kept within 1%.Furthermore,an improvement for measurement system is proposed in order to take advantage of the NP model.展开更多
A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two exam...A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided: 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.展开更多
The microstructures and their kinetics of normal grain growth are simulated using different Monte Carlo (MC) algorithms. Compared with the relative figures and the theoretical normal grain growth exponents of n =0.5...The microstructures and their kinetics of normal grain growth are simulated using different Monte Carlo (MC) algorithms. Compared with the relative figures and the theoretical normal grain growth exponents of n =0.5, the effects of some factors of MC algorithm, i.e. the lattice types, the methods of selecting lattice sites, and the neighbors selection for energy calculations, on the simulation results of grain growth are studied. Two methods of regression were compared, and the three-parameter nonlinear regression is much more suitable for fitting the grain growth kinetics. A better model with appropriate factors included triangular lattice, the attempted site randomly selected, and the first and second nearest neighbors for energy calculations is obtained.展开更多
A Monte Carlo algorithm has been developed by the authors to simulate the chemical vapor deposition (CVD) processes of diamond films. The method considers both the diffusion and the incorporation of the growth radical...A Monte Carlo algorithm has been developed by the authors to simulate the chemical vapor deposition (CVD) processes of diamond films. The method considers both the diffusion and the incorporation of the growth radicals on the growing surface in simulating the evolution of the morphology and microstructure. The calculation of configuration energy is used to determine the orientation of adsorbed growth radicals. The effect of processing variables such as nucleation density and substrate temperature on the morphology and microstructure is discussed. It is found that competitive characteristic and coarsening effect exist in the simulation results, which agree with the experimental observations.展开更多
One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developi...One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.展开更多
The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and ...The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.展开更多
Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using exten...Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using extended continuum configurational-bias (ECCB) method. It is shown that the enrichment of beads near surfaces is happened at high densities due to the bulk packing effect, on the contrary, the depletion is revealed at low densities owing to the configurational entropic contribution. Comparisons with those calculated by density functional theory presented by Cai et al. indicate that the agreement between simulations and predictions is good. Compressibility factors of bulk HSCFs calculated using volume fractions at surfaces were also used to test the reliability of various equations of state of HSCFs by different authors.展开更多
基金Supported by the European Union-NextGenerationEU,through the National Recovery and Resilience Plan of the Republic of Bulgaria,No.BG-RRP-2.004-0008.
文摘In public health,simulation modeling stands as an invaluable asset,enabling the evaluation of new systems without their physical implementation,experimentation with existing systems without operational adjustments,and testing system limits without real-world repercussions.In simulation modeling,the Monte Carlo method emerges as a powerful yet underutilized tool.Although the Monte Carlo method has not yet gained widespread prominence in healthcare,its technological capabilities hold promise for substantial cost reduction and risk mitigation.In this review article,we aimed to explore the transformative potential of the Monte Carlo method in healthcare contexts.We underscore the significance of experiential insights derived from simulated experimentation,especially in resource-constrained scenarios where time,financial constraints,and limited resources necessitate innovative and efficient approaches.As public health faces increasing challenges,incorporating the Monte Carlo method presents an opportunity for enhanced system construction,analysis,and evaluation.
基金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 authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(Nos.51904202,11902212,11901578).
文摘This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.
文摘The modified sum-peak method estimates radioactivity by using only the peak and the sum-peak count rates.To verify the modified sum-peak method,the dependence of the full energy peak efficiency on the source-to-detector distance in a high-purity germanium detector system was studied using a Geant4 Monte Carlo simulation.The effect of the deadlayer in the germanium crystal was estimated by reference to experiments on 241 Am and the relative efficiency of the detector.The peak efficiency dependence on the source-to-detector distance was compared between the simulation and measurements.The modified sum-peak method is discussed with respect to these peak efficiencies.
基金Supported by the National Natural Science Foundation of China under Grant No 61572529
文摘Considering the ocean water's optical attenuation and the roughness of the sea surface, we analyze the security of continuous-variable (CV) quantum key distribution (QKD) based Mr-to-water channel. The effects of the absorp- tion and scattering on the transmittance of underwater quantum channel and the maximum secure transmission distance are studied. Considering the roughness of the sea surface, we simulate the performance bounds of CV QKD with different wind speeds using the Monte Carlo method. The results show that even if the secret key rate gradually reduces as the wind speed increases, the maximum transmission distance will not be affected obviously. Compared to the works regarding short-distance underwater optical communication, our research represents a significant step towards establishing secure communication between air platform and submarine vehicle.
文摘The structure and behaviour of LiF-KF solution,as a typical common-anion system,has been simulated by Monte Carlo method.The calculation of partial radial distribution function of ions,heat of mixing and potential energy distribution shows that the average distance be- tween Li^+ and F^- ions will significantly narrow after mixing of molten LiF and KF.This is very similar to the lean-on-one-side effect in molten LiF-KCl solution.The calculated heat of mixing is in fair agreement with the measured one.The dominant source of the energy of mixing may be that the decrease of the repulsion energy between cations,the decrease of the attraction energy between cations and anions,and the decrease of the repulsion energy be- tween anions.
文摘Pricing the American put options requires solving an optimal stopping problem and therefore is a challenge for the setting up of simulation parameters.This paper uses least square Monte Carlo(LSMC) simulation to price the American put options and output the optimal simulation steps and number of Hermite basis functions.The results suggest:with different time cost and error tolerance,investors can choose the optimal simulation steps and number of basis function individually to price American put options numerically.Generally,with the pre-limitation in the section "least square Monte Carlo simulation",a number of basis equals 4,15 000 simulation steps for Hermite basis function appear to be sufficient for the method.
基金supported by the National Key R&D Program of China(No.2017YFB0504104)
文摘Artificial earthquake catalogue simulation is one of the ways to effectively improve the incompleteness of the existing earthquake catalogue,the scarcity of large earthquake records and the improvement of seismological research. Based on the Poisson distribution model of seismic activity and the Gutenberg-Richter magnitude-frequency relationship,the Monte Carlo method which can describe the characteristics of the stochastic nature and the physical experiment process is used. This paper simulates the future seismic catalogues of the Fenhe-Weihe seismic belt of different durations and conducts statistical tests on them.The analysis shows that the simulation catalogue meets the set seismic activity parameters and meets the Poisson distribution hypothesis,which can obtain a better simulated earthquake catalogues that meets the seismic activity characteristics. According to the simulated earthquake catalogues,future earthquake trends in this region are analyzed to provide reference for seismic hazard analysis.
基金the National High Technology Research and Development Program of China(863 Program:2012AA022602)the 111 Project(B07014)and Tianjin Natural Science Foundation(No.16JCZDJC31200).
文摘Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key to measure the floating reference point(FRP)precisely for the application of FRT.Monte Carlo(MC)simulation has been introduced to quantitatively in-vestigate the effects of positioning errors and light source drifts on measuring FRP.In this article,thinning and calculating method(TCM)is proposed to quantify the positioning error.Mean-while,the normalization process(NP)is developed to significantly reduce the error induced by light source drift.The results according to TCM show that 7 purm deviations in positioning can generate about 10.63%relative error in FRP.It is more noticeable that 1%fluctuation in light source intensity may lead to 12.21%relative errors.Gratifyingly,the proposed NP model can effectively reduce the error caused by light source drift.Therefore,the measurement system for FRPs must meet that the positioning error is less than 7 purm,and the light source drift is kept within 1%.Furthermore,an improvement for measurement system is proposed in order to take advantage of the NP model.
基金Supported by the National Natural Science Foundation of China (60672104,10675013)the Na-tional Basic Research Program of China ("973"Program)(2006CB705705)+1 种基金the 10th Five-Year Plan of the Ministry of Science and Technology of China(2001BA706B-05)the Joint Research Foundation of Beijing Municipal Commissionof Education~~
文摘A method for designing an X-ray flatness filter for medical electron linac is developed. It is used in the optimization process in the electron beam radiation system. Monte Carlo simulation method is used and two examples of real radiation system optimization processes for China-made medical electron linac are provided: 15 MV X- ray system of BJ-20 linac, and 12 MeV electron system of BJ-14. Results are verified by using the traditional method.
基金the International Science & Technology Cooperation Project of Shandong Province(2006)the Natural Science Foundation of Shandong Province(Y2007F06).
文摘The microstructures and their kinetics of normal grain growth are simulated using different Monte Carlo (MC) algorithms. Compared with the relative figures and the theoretical normal grain growth exponents of n =0.5, the effects of some factors of MC algorithm, i.e. the lattice types, the methods of selecting lattice sites, and the neighbors selection for energy calculations, on the simulation results of grain growth are studied. Two methods of regression were compared, and the three-parameter nonlinear regression is much more suitable for fitting the grain growth kinetics. A better model with appropriate factors included triangular lattice, the attempted site randomly selected, and the first and second nearest neighbors for energy calculations is obtained.
文摘A Monte Carlo algorithm has been developed by the authors to simulate the chemical vapor deposition (CVD) processes of diamond films. The method considers both the diffusion and the incorporation of the growth radicals on the growing surface in simulating the evolution of the morphology and microstructure. The calculation of configuration energy is used to determine the orientation of adsorbed growth radicals. The effect of processing variables such as nucleation density and substrate temperature on the morphology and microstructure is discussed. It is found that competitive characteristic and coarsening effect exist in the simulation results, which agree with the experimental observations.
文摘One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.
基金Projects(51206011,U1937201)supported by the National Natural Science Foundation of ChinaProject(20200301040RQ)supported by the Science and Technology Development Program of Jilin Province,China+1 种基金Project(JJKH20190541KJ)supported by the Education Department of Jilin Province,ChinaProject(18DY017)supported by Changchun Science and Technology Program of Changchun City,China。
文摘The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.
基金Supported by the National Science Foundation of China (No. 29736170, No. 20025618) and the Doctoral Research Foundation by Ministry of Education of China (No. 1999025103). Additional support provided by the Visiting Researcher Foundation of University La
文摘Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using extended continuum configurational-bias (ECCB) method. It is shown that the enrichment of beads near surfaces is happened at high densities due to the bulk packing effect, on the contrary, the depletion is revealed at low densities owing to the configurational entropic contribution. Comparisons with those calculated by density functional theory presented by Cai et al. indicate that the agreement between simulations and predictions is good. Compressibility factors of bulk HSCFs calculated using volume fractions at surfaces were also used to test the reliability of various equations of state of HSCFs by different authors.