In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. ...In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.展开更多
Although conventional coal mine designs are conservative regarding pillar strength,local failures such as roof-falls and pillar bursts still affect mine safety and operations.Previous studies have identified that disc...Although conventional coal mine designs are conservative regarding pillar strength,local failures such as roof-falls and pillar bursts still affect mine safety and operations.Previous studies have identified that discontinuous,layered roof materials have some self-supporting capacity.This research is a preliminary step towards understanding these mechanics in coal-measure rocks.Although others have considered broad conceptual models and simplified analogs for mine roof behavior,this study presents a unique numerical model that more completely represents in-situ roof conditions.The discrete element method(DEM)is utilized to conduct a parametric analysis considering a range of in-situ stress ratios,material properties,and joint networks to determine the parameters controlling the stability of single-entries modeled in two-dimensions.Model results are compared to empirical observations of roof-support effectiveness(ARBS)in the context of the coal mine roof rating(CMRR)system.Results such as immediate roof displacement,overall stability,and statistical relationships between model parameters and outcomes are presented herein.Potential practical applications of this line of research include:(1)roof-support optimization for a range of coal-measure rocks,(2)establishment of a relationship between roof stability and pillar stress,and(3)determination of which parameters are most critical to roof stability and therefore require concentrated evaluation.展开更多
This study investigated the correlations between mechanical properties and mineralogy of granite using the digital image processing(DIP) and discrete element method(DEM). The results showed that the X-ray diffraction(...This study investigated the correlations between mechanical properties and mineralogy of granite using the digital image processing(DIP) and discrete element method(DEM). The results showed that the X-ray diffraction(XRD)-based DIP method effectively analyzed the mineral composition contents and spatial distributions of granite. During the particle flow code(PFC2D) model calibration phase, the numerical simulation exhibited that the uniaxial compressive strength(UCS) value, elastic modulus(E), and failure pattern of the granite specimen in the UCS test were comparable to the experiment. By establishing 351 sets of numerical models and exploring the impacts of mineral composition on the mechanical properties of granite, it indicated that there was no negative correlation between quartz and feldspar for UCS, tensile strength(σ_(t)), and E. In contrast, mica had a significant negative correlation for UCS, σ_(t), and E. The presence of quartz increased the brittleness of granite, whereas the presence of mica and feldspar increased its ductility in UCS and direct tensile strength(DTS) tests. Varying contents of major mineral compositions in granite showed minor influence on the number of cracks in both UCS and DTS tests.展开更多
The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in...The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.展开更多
A series of accidents caused by crowds within the last decades evoked a lot of scientific interest in modeling the movement of pedestrian crowds. Based on the discrete element method, a granular dynamic model, in whic...A series of accidents caused by crowds within the last decades evoked a lot of scientific interest in modeling the movement of pedestrian crowds. Based on the discrete element method, a granular dynamic model, in which the human body is simplified as a self-driven sphere, is proposed to simulate the characteristics of crowd flow through an exit. In this model, the repulsive force among people is considered to have an anisotropic feature, and the physical contact force due to body deformation is quantified by the Hertz contact model. The movement of the human body is simulated by applying the second Newton＇s law. The crowd flow through an exit at different desired velocities is studied and simulation results indicated that crowd flow exhibits three distinct states, i.e., smooth state, transition state and phase separation state. In the simulation, the clogging phenomenon occurs more easily when the desired velocity is high and the exit may as a result be totally blocked at a desired velocity of 1.6 m/s or above, leading to faster-to-frozen effect.展开更多
Mine overburden dumps have posed significant safety issues in the operations of various unit operations of open pit min-ing especially the external dumps. The external dumps are composed of a mixture of fragmented roc...Mine overburden dumps have posed significant safety issues in the operations of various unit operations of open pit min-ing especially the external dumps. The external dumps are composed of a mixture of fragmented rocks and loose soil. Their charac-teristic is comparable to heavily discontinuous solid mass. The conventional approach of limit equilibrium methods provide safety factors for the slope but nothing about the stress-strain characteristics of the large dump mass. The designs of dump location and their respective geometry are integrated for the know-how of the stability characteristics of these dumps. The discrete element method uses a circular disk to represent the granular solid mass and their interactions are described by the Newton’s third law of motion. The displacement is described by the sliding of the circular disk. This work is focused on the modeling efficiency of the discrete element methods to represent the behaviour of mine dump masses with the specified joint plane for the limit equilibrium method. The advantage of the work lies on the ease of information retrieval at any point at the dump mass concerning the stress and strain histories, displacement, failures etc. which when integrated produces a better understanding of the stability of the dump masses.展开更多
The pebble-bed reactor is one of the most promising designs for the nuclear energy industry. In this paper,a discrete element method-computational fluid dynamics(DEM-CFD) approach that includes thermal conduction, rad...The pebble-bed reactor is one of the most promising designs for the nuclear energy industry. In this paper,a discrete element method-computational fluid dynamics(DEM-CFD) approach that includes thermal conduction, radiation, and natural convection mechanisms was proposed to simulate the thermal-fluid phenomena after the failure of forced circulation cooling system in a pebble-bed core. The whole large-scale packed bed was created using the DEM technique, and the calculated radial porosity of the bed was validated with empirical correlations reported by researchers. To reduce computational costs, a segment of the bed was extracted, which served as a good representative of the large-scale packed bed for CFD calculation. The temperature distributions simulated with two different fluids in this DEM-CFD approach were in good agreement with SANA experimental data. The influence of the natural convection mechanism on heat transfer must be taken into account for coolants with strong convective capacity. The proposed DEM-CFD methodology offers a computationally efficient and widely applied method for understanding the heat transfer process in a pebble-bed core. The method can also be easily extended to assess the passive safety features of newly designed fluoride-salt-cooled pebble-bed reactors.展开更多
To simulate ballast performance accurately and efficiently,the input in discrete element models should be carefully selected,including the contact model and applied particle shape.To study the effects of the contact m...To simulate ballast performance accurately and efficiently,the input in discrete element models should be carefully selected,including the contact model and applied particle shape.To study the effects of the contact model and applied particle shape on the ballast performance(shear strength and deformation),the direct shear test(DST)model and the large-scale process simulation test(LPST)model were developed on the basis of two types of contact models,namely the rolling resistance linear(RRL)model and the linear contact(LC)model.Particle shapes are differentiated by clumps.A clump is a sphere assembly for one ballast particle.The results show that compared with the typical LC model,the RRL method is more efficient and realistic to predict shear strength results of ballast assemblies in DSTs.In addition,the RRL contact model can also provide accurate vertical and lateral ballast deformation under the cyclic loading in LPSTs.展开更多
The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on th...The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.展开更多
With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and re...With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.展开更多
船舶在冰区海域中航行会受到冰水环境阻力的影响,是冰区船舶快速性研究中的重要影响因素。为合理分析冰区船舶的快速性能,该文采用基于离散元(discrete element method,DEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)...船舶在冰区海域中航行会受到冰水环境阻力的影响,是冰区船舶快速性研究中的重要影响因素。为合理分析冰区船舶的快速性能,该文采用基于离散元(discrete element method,DEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法的流固耦合模型模拟船舶冰区航行过程,获得不同航速下的船舶阻力和推进力,进一步计算出螺旋桨的推力、扭矩以及定速航行所需的螺旋桨转速等参数。为研究船体结构、海冰与海水之间的流固耦合作用,文中通过SPH粒子与固定粒子边界相对运动的拟合项直接计算固体与流体之间的相互作用力,建立船体结构、海冰与海水耦合的DEM-SPH模型,并基于该模型分别对船舶在冰区的航行阻力和推进力进行模拟,通过拟合的方式匹配航行阻力和推进力,并考虑尾部流场导致的船体阻力增额,从而预报船舶在特定航速下实现自航所需的螺旋桨转速。此外,文中还模拟了DTMB 5415船模在浮冰区和层冰区中航行的阻力和不同螺旋桨转速下的推力,对船模在不同工况下实现特定航速航行所需的螺旋桨转速进行了预报。计算结果表明:DEM-SPH耦合模型对船-冰、桨-冰作用中的流固耦合过程模拟效果出色,可完整描述船体及尾部伴流场对海冰的拖曳作用;通过文中所述阻力-推力模拟算例及强制力的拟合分析,所形成的基于数值模拟方法的船舶自航下螺旋桨转速预报,可为进一步的试验验证和工程应用推广奠定基础。展开更多
文摘In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.
基金sponsored by the Alpha Foundation for the Improvement of Mine Safety and Health, Inc. (Alpha Foundation)the funding provided for this project by the Alpha Foundationpartially funded by the National Institute of Occupational Health and Science (NIOSH) under Grant Number 200-2016-90154.
文摘Although conventional coal mine designs are conservative regarding pillar strength,local failures such as roof-falls and pillar bursts still affect mine safety and operations.Previous studies have identified that discontinuous,layered roof materials have some self-supporting capacity.This research is a preliminary step towards understanding these mechanics in coal-measure rocks.Although others have considered broad conceptual models and simplified analogs for mine roof behavior,this study presents a unique numerical model that more completely represents in-situ roof conditions.The discrete element method(DEM)is utilized to conduct a parametric analysis considering a range of in-situ stress ratios,material properties,and joint networks to determine the parameters controlling the stability of single-entries modeled in two-dimensions.Model results are compared to empirical observations of roof-support effectiveness(ARBS)in the context of the coal mine roof rating(CMRR)system.Results such as immediate roof displacement,overall stability,and statistical relationships between model parameters and outcomes are presented herein.Potential practical applications of this line of research include:(1)roof-support optimization for a range of coal-measure rocks,(2)establishment of a relationship between roof stability and pillar stress,and(3)determination of which parameters are most critical to roof stability and therefore require concentrated evaluation.
基金This research was supported by the Department of Mining Engineering at the University of Utah.In addition,the lead author wishes to acknowledge the financial support received from the Talent Introduction Project,part of the Elite Program of Shandong University of Science and Technology(No.0104060540171).
文摘This study investigated the correlations between mechanical properties and mineralogy of granite using the digital image processing(DIP) and discrete element method(DEM). The results showed that the X-ray diffraction(XRD)-based DIP method effectively analyzed the mineral composition contents and spatial distributions of granite. During the particle flow code(PFC2D) model calibration phase, the numerical simulation exhibited that the uniaxial compressive strength(UCS) value, elastic modulus(E), and failure pattern of the granite specimen in the UCS test were comparable to the experiment. By establishing 351 sets of numerical models and exploring the impacts of mineral composition on the mechanical properties of granite, it indicated that there was no negative correlation between quartz and feldspar for UCS, tensile strength(σ_(t)), and E. In contrast, mica had a significant negative correlation for UCS, σ_(t), and E. The presence of quartz increased the brittleness of granite, whereas the presence of mica and feldspar increased its ductility in UCS and direct tensile strength(DTS) tests. Varying contents of major mineral compositions in granite showed minor influence on the number of cracks in both UCS and DTS tests.
基金supported by the National Natural Science Foundation of China(No.11975097)the Fundamental Research Funds for the Central Universities(No.2019MS038).
文摘The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.71473207,51178445,and 71103148)the Research Grant Council,Government of Hong Kong,China(Grant No.City U119011)the Fundamental Research Funds for the Central Universities,China(Grant Nos.2682014CX103 and 2682014RC05)
文摘A series of accidents caused by crowds within the last decades evoked a lot of scientific interest in modeling the movement of pedestrian crowds. Based on the discrete element method, a granular dynamic model, in which the human body is simplified as a self-driven sphere, is proposed to simulate the characteristics of crowd flow through an exit. In this model, the repulsive force among people is considered to have an anisotropic feature, and the physical contact force due to body deformation is quantified by the Hertz contact model. The movement of the human body is simulated by applying the second Newton＇s law. The crowd flow through an exit at different desired velocities is studied and simulation results indicated that crowd flow exhibits three distinct states, i.e., smooth state, transition state and phase separation state. In the simulation, the clogging phenomenon occurs more easily when the desired velocity is high and the exit may as a result be totally blocked at a desired velocity of 1.6 m/s or above, leading to faster-to-frozen effect.
文摘Mine overburden dumps have posed significant safety issues in the operations of various unit operations of open pit min-ing especially the external dumps. The external dumps are composed of a mixture of fragmented rocks and loose soil. Their charac-teristic is comparable to heavily discontinuous solid mass. The conventional approach of limit equilibrium methods provide safety factors for the slope but nothing about the stress-strain characteristics of the large dump mass. The designs of dump location and their respective geometry are integrated for the know-how of the stability characteristics of these dumps. The discrete element method uses a circular disk to represent the granular solid mass and their interactions are described by the Newton’s third law of motion. The displacement is described by the sliding of the circular disk. This work is focused on the modeling efficiency of the discrete element methods to represent the behaviour of mine dump masses with the specified joint plane for the limit equilibrium method. The advantage of the work lies on the ease of information retrieval at any point at the dump mass concerning the stress and strain histories, displacement, failures etc. which when integrated produces a better understanding of the stability of the dump masses.
基金supported by the Chinese TMSR Strategic Pioneer Science and Technology Project(No.XDA02010000)the Frontier Science Key Program of the Chinese Academy of Sciences(No.QYZDY-SSW-JSC016)
文摘The pebble-bed reactor is one of the most promising designs for the nuclear energy industry. In this paper,a discrete element method-computational fluid dynamics(DEM-CFD) approach that includes thermal conduction, radiation, and natural convection mechanisms was proposed to simulate the thermal-fluid phenomena after the failure of forced circulation cooling system in a pebble-bed core. The whole large-scale packed bed was created using the DEM technique, and the calculated radial porosity of the bed was validated with empirical correlations reported by researchers. To reduce computational costs, a segment of the bed was extracted, which served as a good representative of the large-scale packed bed for CFD calculation. The temperature distributions simulated with two different fluids in this DEM-CFD approach were in good agreement with SANA experimental data. The influence of the natural convection mechanism on heat transfer must be taken into account for coolants with strong convective capacity. The proposed DEM-CFD methodology offers a computationally efficient and widely applied method for understanding the heat transfer process in a pebble-bed core. The method can also be easily extended to assess the passive safety features of newly designed fluoride-salt-cooled pebble-bed reactors.
基金by the China Scholarship Council and the Natural Science Foundation of China(Grant No.51578469)We also would like to acknowledge the support of the Chinese Program of Introducing Talents of Discipline to Universities(111 Project,Grant No.B16041)。
文摘To simulate ballast performance accurately and efficiently,the input in discrete element models should be carefully selected,including the contact model and applied particle shape.To study the effects of the contact model and applied particle shape on the ballast performance(shear strength and deformation),the direct shear test(DST)model and the large-scale process simulation test(LPST)model were developed on the basis of two types of contact models,namely the rolling resistance linear(RRL)model and the linear contact(LC)model.Particle shapes are differentiated by clumps.A clump is a sphere assembly for one ballast particle.The results show that compared with the typical LC model,the RRL method is more efficient and realistic to predict shear strength results of ballast assemblies in DSTs.In addition,the RRL contact model can also provide accurate vertical and lateral ballast deformation under the cyclic loading in LPSTs.
基金This work was supported by the National Key R&D Program‘Transportation Infrastructure’project(No.2022YFB2603400).
文摘The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.
文摘With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.
