This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fis...This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fissured granite samples is recalled.Then,PFC3D is introduced,with focus on the bonded particle models(BPM).After that,we present previous studies where intact rock is simulated by means of flatjoint approaches,and how improved accuracy was gained with the help of parametric studies.Then,models of the pre-fissured rock specimens were generated,including modeled fissures in the form of“smooth joint”type contacts.Finally,triaxial testing simulations of 1 t 2 and 2 t 3 jointed rock specimens were performed.Results show that both elastic behavior and the peak strength levels are closely matched,without any additional fine tuning of micro-mechanical parameters.Concerning the postfailure behavior,models reproduce the trends of decreasing dilation with increasing confinement and plasticity.However,the dilation values simulated are larger than those observed in practice.This is attributed to the difficulty in modeling some phenomena of fissured rock behaviors,such as rock piece corner crushing with dust production and interactions between newly formed shear bands or axial splitting cracks with pre-existing joints.展开更多
Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions(PFC^(2D)) method, and the simulation results agreed with the experimental data. Th...Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions(PFC^(2D)) method, and the simulation results agreed with the experimental data. The particle flow code method can simulate the mechanical properties of the polymer. The triaxial cyclic loading tests of the polymer material under different confining pressures were carried out via PFC^(2D) to analyze its mechanical performance. The PFC^(2D) simulation results show that the value of the elastic modulus of the polymer decreases slowly at first and fluctuated within a narrow range near the value of the peak strength; the cumulative plastic strain increases slowly at first and then increases rapidly; the peak strength and elastic modulus of polymer increase with the confining pressure; the PFC^(2D) method can be used to quantitatively evaluate the damage behavior of the polymer material and estimate the fatigue life of the materials under fatigue load based on the number and the location of micro-cracks. Thus, the PFC^(2D) method is an effective tool to study polymers.展开更多
A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck ...A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.展开更多
Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the...Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the development process were investigated.On the basis of the scanning electron microscopy(SEM)and mechanical test results of CPB,the particle size information of CPB was extracted,and a two-dimensional particle flow code(PFC)model of CPB was established to analyze the evolution rule of mesoscopic parameters during CPB development.The embedded FISH language in PFC was used to develop a program for establishing a PFC model on the basis of the SEM results.The mesoscopic parameters of CPB samples at different curing times,such as coordination number(C_(n)),contact force chain,and rose diagram,were obtained by recording and loading and used to analyze the intrinsic relationship between mesoscopic parameter variations and macroscopic mechanical response during CPB development.It is of considerable significance to establish the physical model of CPB using the PFC to reveal the mesoscopic structure of CPB.展开更多
Based on three-dimensional cellular automata (CA), a new stochastic simulation model to simulate the microstructures and particle flow of talus deposit is proposed. Ill addition, an auto-modeling program CARS is dev...Based on three-dimensional cellular automata (CA), a new stochastic simulation model to simulate the microstructures and particle flow of talus deposit is proposed. Ill addition, an auto-modeling program CARS is developed, with which nunaerical simulations can be conducted conveniently. For the problem of simulating mechanical behaviors of talus deposit, spatial anangement or sphere shapes should be considered. In the new modeling method, four sphere anangement models are developed for the particle flow simulation of talus deposit. Numerical results show that the talus deposit has the mechanical characteristics of typical stress-strain curves, as other rock-like materials. The cohesion of talus deposit decreases with increasing rock content, while the internal friction angle increases with increasing rock contents. Finally, numerical simulation is verified with the results of field test.展开更多
Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carr...Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.展开更多
Vibrating flip-flow screens(VFFS)with stretchable polyurethane sieve mats have been widely used in screening fine-grained materials in recent years.In this work,the discrete element method(DEM)is used to study the scr...Vibrating flip-flow screens(VFFS)with stretchable polyurethane sieve mats have been widely used in screening fine-grained materials in recent years.In this work,the discrete element method(DEM)is used to study the screening process in VFFS to explain particle flow and separation behavior at the particle scale.Unlike traditional vibrating screens,for VFFS,the amplitude response of each point on the elastic sieve mat is different everywhere.This study measures the kinematics of the elastic sieve mat under different conditions such as different stretched lengths and material loads.To establish the elastic sieve mat model in a DEM simulation,the continuous elastic sieve mat is discretized into multiple units,and the displacement signal of each unit tested is analyzed by Fourier series.The Fourier series analysis results of each unit are used as the setting parameters for motion.In this way,the movement of the elastic sieve mat is approximately simulated,and a DEM model of VFFS is produced.Through the simulation,the flow and separation of different-sized particles in VFFS are studied,and the reasonability of the simulation is verified by a pilot-scale screening experiment.The present study demonstrates the potential of the DEM method for the analysis of screening processes in VFFS.展开更多
The thermomechanical coupling of rocks refers to the interaction between the mechanical and thermodynamic behaviors of rocks induced by temperature changes.The study of this coupling interaction is essential for under...The thermomechanical coupling of rocks refers to the interaction between the mechanical and thermodynamic behaviors of rocks induced by temperature changes.The study of this coupling interaction is essential for understanding the mechanical and thermodynamic properties of the surrounding rocks in underground engineering.In this study,an improved temperature-dependent linear parallel bond model is introduced under the framework of a particle flow simulation.A series of numerical thermomechanical coupling tests are then conducted to calibrate the micro-parameters of the proposed model by considering the mechanical behavior of the rock under different thermomechanical loadings.Good agreement between the numerical results and experimental data are obtained,particularly in terms of the compression,tension,and elastic responses of granite.With this improved model,the thermodynamic response and underlying cracking behavior of a deep-buried tunnel under different thermal loading conditions are investigated and discussed in detail.展开更多
Using tunnel boring machines to excavate high-strength intact rock masses is becoming more common.Due to the interactions between disc cutters and rocks,abnormal wear of disc cutters,especially cutter chipping,has bec...Using tunnel boring machines to excavate high-strength intact rock masses is becoming more common.Due to the interactions between disc cutters and rocks,abnormal wear of disc cutters,especially cutter chipping,has become a common phenomenon.Existing research has mainly focused on normal wear of disc cutters without addressing abnormal wear cases.This study used the disc cutter consumption data of a tunnel project in China to investigate the abovementioned problem based on field research.According to the fail-ure patterns and fracture surface characteristics,the cutter chipping patterns were mainly categorized into four types:granule chipping,patch chipping,primary collapse,and secondary collapse.To further simulate the evolution of disc cutter chipping,based on the linear plastic bond model,a new contact model called the modified plastic bond(MPB)model was proposed to solve the metal simulation prob-lem in Particle Flow Code software.To this end,a set of uniaxial tensile and compressive tests were initially conducted to verify the applicability of the MPB model.Then,a series of three-dimensional rock-cutting simulation tests were conducted to reflect the evolu-tionary processes involved in each type of cutter chipping.The cutter chipping mechanism and morphological characteristics were clas-sified and summarized in detail.The results revealed that the cutting speed and penetration growth led to a rising trend in the probability and intensity of the cutter chipping.The presence of initial defects also induced an adverse effect on the service life of the cutter.The results indicated suitable working conditions for the cutter and suggested ways to control tunneling parameters and avoid frequent cutter chipping cases.展开更多
In this work,by establishing a three-dimensional physical model of a 1000-ton industrial multi-jet combustion reactor,a hexahedral structured grid was used to discretize the model.Combined with realizable k–εmodel,e...In this work,by establishing a three-dimensional physical model of a 1000-ton industrial multi-jet combustion reactor,a hexahedral structured grid was used to discretize the model.Combined with realizable k–εmodel,eddy-dissipation-concept,discrete-ordinate radiation model,hydrogen 19-step detailed reaction mechanism,air age user-defined-function,velocity field,temperature field,concentration field and gas arrival time in the reactor were numerically simulated.The Euler–Lagrange method combined with the discrete-phase-model was used to reveal the flow characteristics of particles in the reactor,and based on this,the effects of the reactor aspect ratios,central jet gas velocity and particle size on the flow field characteristics and particle back-mixing degree in the reactor were investigated.The results show that with the decrease of aspect ratio in the combustion reactors,the velocity and temperature attenuation in the reactor are intensified,the vortex phenomenon is aggravated,and the residence time distribution of nanoparticles is more dispersed.With the increase in the central jet gas velocities in reactors,the vortex lengthens along the axis,the turbulence intensity increases,and the residence time of particles decreases.The back-mixing degree and residence time of particles in the reactor also decrease with the increase in particle size.The simulation results can provide reference for the structural regulation of nanoparticles and the structural design of combustion reactor in the process of gas combustion synthesis.展开更多
In a gas/particle two-phase test facility, a three-component particle-dynamics anemometer was used to measure the characteristics of gas/particle two-phase flows in a 29 megawatt (MW) pulverized coal industrial boiler...In a gas/particle two-phase test facility, a three-component particle-dynamics anemometer was used to measure the characteristics of gas/particle two-phase flows in a 29 megawatt (MW) pulverized coal industrial boiler equipped with a new type of swirling pulverized coal burner. The distributions of three-dimensional gas/particle velocity, particle volume flux, and particle size distribution were measured under different working conditions. The mean axial velocity and the particle volume flux in the central region of the burner outlet were found to be negative. This indicated that a central recirculation zone was formed in the center of the burner. In the central recirculation zone, the absolute value of the mean axial velocity and the particle volume flux increased when the external secondary air volume increased. The size of the central reflux zone remained stable when the air volume ratio changed. Along the direction of the jet, the peak value formed by the tertiary air gradually moved toward the center of the burner. This tertiary air was mixed with the peak value formed by the air in the adiabatic combustion chamber after the cross-section of x/d = 0.7. Large particles were concentrated near the wall area, and the particle size in the recirculation zone was small.展开更多
Swirl burner design was optimized by simulating swirl gas\|particle flows with different swirl numbers at the exit of a small\|scale swirl burner for pulverized\|coal furnaces using the k\|ε\|k p model. The pred...Swirl burner design was optimized by simulating swirl gas\|particle flows with different swirl numbers at the exit of a small\|scale swirl burner for pulverized\|coal furnaces using the k\|ε\|k p model. The predicted two\|phase time\|averaged velocities and particle concentration distributions for several different cases were compared to improve the design. The effect of the swirl number on the two\|phase velocities and particle concentration was investigated. The results give the two\|phase axial and tangential time\|averaged and fluctuation velocities and particle concentrations, showing that large recirculation zones of gas and particles forms in the near\|axis region of the burner exit, but the particle concentration in the recirculating zone is very low.展开更多
In order to better understand the mechanical properties of graded crushed rocks (GCRs) and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed fo...In order to better understand the mechanical properties of graded crushed rocks (GCRs) and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio (CBR) test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson's ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson's ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.展开更多
A new type of swirl burner has been developed to stabilize pulverized\|coal combustion by burning different types of coal at different loads and to reduce NO x formation during combustion. The burner uses a device to...A new type of swirl burner has been developed to stabilize pulverized\|coal combustion by burning different types of coal at different loads and to reduce NO x formation during combustion. The burner uses a device to concentrate the coal powder in the primary\|air tube that divides the primary coal\|air into two streams with different pulverized\|coal concentrations. This paper reports the measurement of gas\|particle flows at the exit of the different swirl burners using a 3\|D Phase Doppler Particle Anemometer (PDPA). The effect of different geometrical configurations on the two\|phase flow field is studied. The results that give the two\|phase flow fields and particle concentrations show the superiority of the new swirl burner.展开更多
This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbule...This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbulence effects inside the separator were modeled using the Reynolds stress model. The two phase gas–solid particles flow was modeled using a hybrid Euler–Lagrange approach, which accounts for the four-way coupling between phases. The simulations were performed for three inlet velocities of the gaseous phase and several cyclone mass particle loadings. Moreover, the influences of several submodel parameters on the calculated results were investigated. The obtained results were compared against experimental data collected at the in-house experimental rig. The cyclone pressure drop evaluated numerically underpredicts the measured values. The possible reason of this discrepancies was disused.展开更多
The double blade pump is widely used in sewage treatment industry,however,the research on the internal flow characteristics of the double blade pump with particle image velocimetry(PIV) technology is very little at ...The double blade pump is widely used in sewage treatment industry,however,the research on the internal flow characteristics of the double blade pump with particle image velocimetry(PIV) technology is very little at present.To reveal inner flow characteristics in double blade pump impeller under off-design and design conditions,inner flows in a double blade pump impeller,whose specific speed is 111,are measured under the five off-design conditions and design condition by using 3D PIV test technology.In order to ensure the accuracy of the 3D PIV test,the external trigger synchronization system which makes use of fiber optic and equivalent calibration method are applied.The 3D PIV relative velocity synthesis procedure is compiled by using Visual C++ 2005.Then absolute velocity distribution and relative velocity distribution in the double blade pump impeller are obtained.Test results show that vortex exists in each condition,but the location,size and velocity of vortex core are different.Average absolute velocity value of impeller outlet increases at first,then decreases,and then increases again with increase of flow rate.Again average relative velocity values under 0.4,0.8,and 1.2 design condition are higher than that under 1.0 design condition,while under 0.6 and 1.4 design condition it is lower.Under low flow rate conditions,radial vectors of absolute velocities at impeller outlet and blade inlet near the pump shaft decrease with increase of flow rate,while that of relative velocities at the suction side near the pump shaft decreases.Radial vectors of absolute velocities and relative velocities change slightly under the two large flow rate conditions.The research results can be applied to instruct the hydraulic optimization design of double blade pumps.展开更多
For any study ofa suspension entering a pore, the knowledge of the force and moment exerted on a solute particle in an arbitrary position outside the pore is essential, 'This paper for the first lime presents appr...For any study ofa suspension entering a pore, the knowledge of the force and moment exerted on a solute particle in an arbitrary position outside the pore is essential, 'This paper for the first lime presents approximate analytical expressions (in closed form) of all the twelve force and moment coefficienis for a sphere outsied a circular orifice, on the basis of a number of discrete data computed by Yan et al(1987).These coefficients are then applied to calculate the trajectory and angular velocity of a spherical particle approaching the pore at zero Reynolds number. The trajectory is in excellent agreement with the available experimental results. An analysis of the relative importance of the coefficients shows that the rotation effect cannot be neglected near the pore opening or near the wall, and that the lateral force effect must be taken into account in the neighborhood of the edge of the pore opening. It is due to neglecting these factors that previous theoretical results deviate from the experimental ones near the pore opening. The effects of the ratio of the particle to pore radii as well as the influences of the graritytbuoyance on the particle trajectory, velocity distribution and rotation are discnssed in detail. It is pointed out that in the experiments of neutrally-buoyant suspensions, the restriction on the density of the particle is most demanding for a large particle size.The expressions of forces and moments presenled herein are complete, relatively accurate and convenient, thus providing a good prerequisite for further studies of any problems involving the entrance of particles to a pare.展开更多
The elliptic flow v2, for π±, K±, p and p in Au+Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the ...The elliptic flow v2, for π±, K±, p and p in Au+Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the v2 values for p and p is observed, and the values of v2 splitting are larger compared with π+ and π-, K+ and K-. The difference increases with decreasing the center-of-mass energy. The effect of the quark coalescence mechanism in a multi-phase transport model to the value of elliptic difference △v2 between p and p- has been discussed. The simulation of Au+Au collisions at 14.5 GeV shows that the effect of hadron cascade to △v2 is not obvious, and a larger patton-scattering cross section can lead to a larger △v2.展开更多
The Fenwei Basin, covered by loess, experiences severe ground fissure disasters. These disasters disrupt the continuity of the loess and pose significant threats to engineering construction safety along transportation...The Fenwei Basin, covered by loess, experiences severe ground fissure disasters. These disasters disrupt the continuity of the loess and pose significant threats to engineering construction safety along transportation routes. Nevertheless, the crack characteristics and the influence zone of ground fissures in the loess layer remain inadequately investigated. To effectively prevent and control ground fissure disasters, physical model tests and the PFC(particle flow code) numerical simulation method are used to investigate the crack mechanism of buried ground fissures in the loess layer. The results show that there are two main cracks in the layer profile, which have a Y-shape morphology. As the dip angle of the preset cracks increased from 60° to 90°, the main deformation zone at the surface gradually shifted towards the footwall. The process of crack propagation from depth to surface is divided into five stages. Additionally, the results confirm the accuracy of the width of the rupture zone d2in the footwall calculated by the cantilever beam theory. These findings can offer theoretical guidance for determining the avoidance distance of ground fissures in loess regions, as well as for implementing disaster prevention and corresponding control measures for various stages of buried ground fissure propagation.展开更多
基金The University of Vigo is acknowledged for financing part of the first author’s PhD studiesthe Spanish Ministry of Economy and Competitiveness for funding of the project‘Deepening on the behaviour of rock masses:Scale effects on the stressestrain response of fissured rock samples with particular emphasis on post-failure’,awarded under Contract Reference No.RTI2018-093563-B-I00partially financed by means of European Regional Development Funds from the European Union(EU)。
文摘This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fissured granite samples is recalled.Then,PFC3D is introduced,with focus on the bonded particle models(BPM).After that,we present previous studies where intact rock is simulated by means of flatjoint approaches,and how improved accuracy was gained with the help of parametric studies.Then,models of the pre-fissured rock specimens were generated,including modeled fissures in the form of“smooth joint”type contacts.Finally,triaxial testing simulations of 1 t 2 and 2 t 3 jointed rock specimens were performed.Results show that both elastic behavior and the peak strength levels are closely matched,without any additional fine tuning of micro-mechanical parameters.Concerning the postfailure behavior,models reproduce the trends of decreasing dilation with increasing confinement and plasticity.However,the dilation values simulated are larger than those observed in practice.This is attributed to the difficulty in modeling some phenomena of fissured rock behaviors,such as rock piece corner crushing with dust production and interactions between newly formed shear bands or axial splitting cracks with pre-existing joints.
基金the National Key R&D Program of China(No.2017YFC0405002)
文摘Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions(PFC^(2D)) method, and the simulation results agreed with the experimental data. The particle flow code method can simulate the mechanical properties of the polymer. The triaxial cyclic loading tests of the polymer material under different confining pressures were carried out via PFC^(2D) to analyze its mechanical performance. The PFC^(2D) simulation results show that the value of the elastic modulus of the polymer decreases slowly at first and fluctuated within a narrow range near the value of the peak strength; the cumulative plastic strain increases slowly at first and then increases rapidly; the peak strength and elastic modulus of polymer increase with the confining pressure; the PFC^(2D) method can be used to quantitatively evaluate the damage behavior of the polymer material and estimate the fatigue life of the materials under fatigue load based on the number and the location of micro-cracks. Thus, the PFC^(2D) method is an effective tool to study polymers.
基金financial support from the National Natural Science Foundation of China (No. 51204181)the Research Fund for the Doctoral Program of Higher Education of China (No. 20110095120004)+2 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe Fundamental Research Funds for the Central Universities (Nos. 2011QNA10 and 2010QNB17)the China Postdoctoral Science Foundation (No. 20110491485)
文摘A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.
基金financially supported by the National Natural Science Foundation of China(Nos.51874229,52074212,51674188,51504182,51404191,and 51405381)the Natural Science Basic Research Plan of Shaanxi Province of China(Nos.2015JQ5187,2018JQ5183,and 2018JM5161)+3 种基金the Scientific Research Program funded by the Shaanxi Education Department(No.15JK1466)the China Postdoctoral Science Foundation(No.2015M582685)the Outstanding Youth Science Fund of Xi’an University of Science and Technology(No.2018YQ2-01)supported by the National Research Council of Science&and Technology(NST)grant by the Korea Korean government(MSIP)(No.CRC-16-38502-KICT)。
文摘Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the development process were investigated.On the basis of the scanning electron microscopy(SEM)and mechanical test results of CPB,the particle size information of CPB was extracted,and a two-dimensional particle flow code(PFC)model of CPB was established to analyze the evolution rule of mesoscopic parameters during CPB development.The embedded FISH language in PFC was used to develop a program for establishing a PFC model on the basis of the SEM results.The mesoscopic parameters of CPB samples at different curing times,such as coordination number(C_(n)),contact force chain,and rose diagram,were obtained by recording and loading and used to analyze the intrinsic relationship between mesoscopic parameter variations and macroscopic mechanical response during CPB development.It is of considerable significance to establish the physical model of CPB using the PFC to reveal the mesoscopic structure of CPB.
基金Supported by the National Natural Science Foundation of China(50979030 and 50911130366)
文摘Based on three-dimensional cellular automata (CA), a new stochastic simulation model to simulate the microstructures and particle flow of talus deposit is proposed. Ill addition, an auto-modeling program CARS is developed, with which nunaerical simulations can be conducted conveniently. For the problem of simulating mechanical behaviors of talus deposit, spatial anangement or sphere shapes should be considered. In the new modeling method, four sphere anangement models are developed for the particle flow simulation of talus deposit. Numerical results show that the talus deposit has the mechanical characteristics of typical stress-strain curves, as other rock-like materials. The cohesion of talus deposit decreases with increasing rock content, while the internal friction angle increases with increasing rock contents. Finally, numerical simulation is verified with the results of field test.
基金The project supported by the National Natural Science Foundation of China(50606026 and 50736006)
文摘Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.
基金supported by the Fundamental Research Funds for the Central Universities (grant No.2021YJSHH32)Anhui Province Major Science and Technology Achievements Engineering Research and Development Special Project (grant No.202103c08020007).
文摘Vibrating flip-flow screens(VFFS)with stretchable polyurethane sieve mats have been widely used in screening fine-grained materials in recent years.In this work,the discrete element method(DEM)is used to study the screening process in VFFS to explain particle flow and separation behavior at the particle scale.Unlike traditional vibrating screens,for VFFS,the amplitude response of each point on the elastic sieve mat is different everywhere.This study measures the kinematics of the elastic sieve mat under different conditions such as different stretched lengths and material loads.To establish the elastic sieve mat model in a DEM simulation,the continuous elastic sieve mat is discretized into multiple units,and the displacement signal of each unit tested is analyzed by Fourier series.The Fourier series analysis results of each unit are used as the setting parameters for motion.In this way,the movement of the elastic sieve mat is approximately simulated,and a DEM model of VFFS is produced.Through the simulation,the flow and separation of different-sized particles in VFFS are studied,and the reasonability of the simulation is verified by a pilot-scale screening experiment.The present study demonstrates the potential of the DEM method for the analysis of screening processes in VFFS.
基金supported by the Postgraduate Research&Practice Innovation Program of Jiangsu Province (No.KYCX21_0494)the National Natural Science Foundation of China (Grant Nos.51679071 and 41831278)the Key Laboratory of the Ministry of Education on Safe Mining of Deep Metal Mines (No.DM2019K02).
文摘The thermomechanical coupling of rocks refers to the interaction between the mechanical and thermodynamic behaviors of rocks induced by temperature changes.The study of this coupling interaction is essential for understanding the mechanical and thermodynamic properties of the surrounding rocks in underground engineering.In this study,an improved temperature-dependent linear parallel bond model is introduced under the framework of a particle flow simulation.A series of numerical thermomechanical coupling tests are then conducted to calibrate the micro-parameters of the proposed model by considering the mechanical behavior of the rock under different thermomechanical loadings.Good agreement between the numerical results and experimental data are obtained,particularly in terms of the compression,tension,and elastic responses of granite.With this improved model,the thermodynamic response and underlying cracking behavior of a deep-buried tunnel under different thermal loading conditions are investigated and discussed in detail.
基金supported by the National Natural Science Foundation of China(Grant No.52078377)the Key Field Science and Technology Project of Yunnan Province(Grant No.202002AC080002)supported by the China Atomic Energy Authority(CAEA)through the Geological Disposal Program.
文摘Using tunnel boring machines to excavate high-strength intact rock masses is becoming more common.Due to the interactions between disc cutters and rocks,abnormal wear of disc cutters,especially cutter chipping,has become a common phenomenon.Existing research has mainly focused on normal wear of disc cutters without addressing abnormal wear cases.This study used the disc cutter consumption data of a tunnel project in China to investigate the abovementioned problem based on field research.According to the fail-ure patterns and fracture surface characteristics,the cutter chipping patterns were mainly categorized into four types:granule chipping,patch chipping,primary collapse,and secondary collapse.To further simulate the evolution of disc cutter chipping,based on the linear plastic bond model,a new contact model called the modified plastic bond(MPB)model was proposed to solve the metal simulation prob-lem in Particle Flow Code software.To this end,a set of uniaxial tensile and compressive tests were initially conducted to verify the applicability of the MPB model.Then,a series of three-dimensional rock-cutting simulation tests were conducted to reflect the evolu-tionary processes involved in each type of cutter chipping.The cutter chipping mechanism and morphological characteristics were clas-sified and summarized in detail.The results revealed that the cutting speed and penetration growth led to a rising trend in the probability and intensity of the cutter chipping.The presence of initial defects also induced an adverse effect on the service life of the cutter.The results indicated suitable working conditions for the cutter and suggested ways to control tunneling parameters and avoid frequent cutter chipping cases.
基金supported by the National Natural Science Foundation of China(21978088,91534202,51673063)Shanghai Technology Research Leader(20XD1433600)+4 种基金the Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutes of High Learningthe Basic Research Program of Shanghai(17JC1402300)the Shanghai City Board of education research and innovation projectthe Fundamental Research Funds for the Central Universities(222201718002)provided by Feringa Nobel Prize Scientist Joint Research Center。
文摘In this work,by establishing a three-dimensional physical model of a 1000-ton industrial multi-jet combustion reactor,a hexahedral structured grid was used to discretize the model.Combined with realizable k–εmodel,eddy-dissipation-concept,discrete-ordinate radiation model,hydrogen 19-step detailed reaction mechanism,air age user-defined-function,velocity field,temperature field,concentration field and gas arrival time in the reactor were numerically simulated.The Euler–Lagrange method combined with the discrete-phase-model was used to reveal the flow characteristics of particles in the reactor,and based on this,the effects of the reactor aspect ratios,central jet gas velocity and particle size on the flow field characteristics and particle back-mixing degree in the reactor were investigated.The results show that with the decrease of aspect ratio in the combustion reactors,the velocity and temperature attenuation in the reactor are intensified,the vortex phenomenon is aggravated,and the residence time distribution of nanoparticles is more dispersed.With the increase in the central jet gas velocities in reactors,the vortex lengthens along the axis,the turbulence intensity increases,and the residence time of particles decreases.The back-mixing degree and residence time of particles in the reactor also decrease with the increase in particle size.The simulation results can provide reference for the structural regulation of nanoparticles and the structural design of combustion reactor in the process of gas combustion synthesis.
基金This work was supported by National Key Research and Development Program of China(Grant No.2017YF C0212500).
文摘In a gas/particle two-phase test facility, a three-component particle-dynamics anemometer was used to measure the characteristics of gas/particle two-phase flows in a 29 megawatt (MW) pulverized coal industrial boiler equipped with a new type of swirling pulverized coal burner. The distributions of three-dimensional gas/particle velocity, particle volume flux, and particle size distribution were measured under different working conditions. The mean axial velocity and the particle volume flux in the central region of the burner outlet were found to be negative. This indicated that a central recirculation zone was formed in the center of the burner. In the central recirculation zone, the absolute value of the mean axial velocity and the particle volume flux increased when the external secondary air volume increased. The size of the central reflux zone remained stable when the air volume ratio changed. Along the direction of the jet, the peak value formed by the tertiary air gradually moved toward the center of the burner. This tertiary air was mixed with the peak value formed by the air in the adiabatic combustion chamber after the cross-section of x/d = 0.7. Large particles were concentrated near the wall area, and the particle size in the recirculation zone was small.
文摘Swirl burner design was optimized by simulating swirl gas\|particle flows with different swirl numbers at the exit of a small\|scale swirl burner for pulverized\|coal furnaces using the k\|ε\|k p model. The predicted two\|phase time\|averaged velocities and particle concentration distributions for several different cases were compared to improve the design. The effect of the swirl number on the two\|phase velocities and particle concentration was investigated. The results give the two\|phase axial and tangential time\|averaged and fluctuation velocities and particle concentrations, showing that large recirculation zones of gas and particles forms in the near\|axis region of the burner exit, but the particle concentration in the recirculating zone is very low.
基金supported by the Program for New Century Excellent Talents in University (NCET-08-0749)Fundamental Research Funds for the Central Universities (CHD2012JC054)
文摘In order to better understand the mechanical properties of graded crushed rocks (GCRs) and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio (CBR) test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson's ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson's ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.
文摘A new type of swirl burner has been developed to stabilize pulverized\|coal combustion by burning different types of coal at different loads and to reduce NO x formation during combustion. The burner uses a device to concentrate the coal powder in the primary\|air tube that divides the primary coal\|air into two streams with different pulverized\|coal concentrations. This paper reports the measurement of gas\|particle flows at the exit of the different swirl burners using a 3\|D Phase Doppler Particle Anemometer (PDPA). The effect of different geometrical configurations on the two\|phase flow field is studied. The results that give the two\|phase flow fields and particle concentrations show the superiority of the new swirl burner.
文摘This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbulence effects inside the separator were modeled using the Reynolds stress model. The two phase gas–solid particles flow was modeled using a hybrid Euler–Lagrange approach, which accounts for the four-way coupling between phases. The simulations were performed for three inlet velocities of the gaseous phase and several cyclone mass particle loadings. Moreover, the influences of several submodel parameters on the calculated results were investigated. The obtained results were compared against experimental data collected at the in-house experimental rig. The cyclone pressure drop evaluated numerically underpredicts the measured values. The possible reason of this discrepancies was disused.
基金supported by National Outstanding Young Scientists Founds of China (Grant No. 50825902)National Key Technology R&D Program of China (Grant No. 2008BAF34B15)+2 种基金National Natural Science Foundation of China (Grant No. 51079062)Jiangsu Provincial 333 High-level Personal Training Project of ChinaJiangsu Provincial Six Talent Peak of China
文摘The double blade pump is widely used in sewage treatment industry,however,the research on the internal flow characteristics of the double blade pump with particle image velocimetry(PIV) technology is very little at present.To reveal inner flow characteristics in double blade pump impeller under off-design and design conditions,inner flows in a double blade pump impeller,whose specific speed is 111,are measured under the five off-design conditions and design condition by using 3D PIV test technology.In order to ensure the accuracy of the 3D PIV test,the external trigger synchronization system which makes use of fiber optic and equivalent calibration method are applied.The 3D PIV relative velocity synthesis procedure is compiled by using Visual C++ 2005.Then absolute velocity distribution and relative velocity distribution in the double blade pump impeller are obtained.Test results show that vortex exists in each condition,but the location,size and velocity of vortex core are different.Average absolute velocity value of impeller outlet increases at first,then decreases,and then increases again with increase of flow rate.Again average relative velocity values under 0.4,0.8,and 1.2 design condition are higher than that under 1.0 design condition,while under 0.6 and 1.4 design condition it is lower.Under low flow rate conditions,radial vectors of absolute velocities at impeller outlet and blade inlet near the pump shaft decrease with increase of flow rate,while that of relative velocities at the suction side near the pump shaft decreases.Radial vectors of absolute velocities and relative velocities change slightly under the two large flow rate conditions.The research results can be applied to instruct the hydraulic optimization design of double blade pumps.
基金Project supported by the National Natural Science Foundation of China
文摘For any study ofa suspension entering a pore, the knowledge of the force and moment exerted on a solute particle in an arbitrary position outside the pore is essential, 'This paper for the first lime presents approximate analytical expressions (in closed form) of all the twelve force and moment coefficienis for a sphere outsied a circular orifice, on the basis of a number of discrete data computed by Yan et al(1987).These coefficients are then applied to calculate the trajectory and angular velocity of a spherical particle approaching the pore at zero Reynolds number. The trajectory is in excellent agreement with the available experimental results. An analysis of the relative importance of the coefficients shows that the rotation effect cannot be neglected near the pore opening or near the wall, and that the lateral force effect must be taken into account in the neighborhood of the edge of the pore opening. It is due to neglecting these factors that previous theoretical results deviate from the experimental ones near the pore opening. The effects of the ratio of the particle to pore radii as well as the influences of the graritytbuoyance on the particle trajectory, velocity distribution and rotation are discnssed in detail. It is pointed out that in the experiments of neutrally-buoyant suspensions, the restriction on the density of the particle is most demanding for a large particle size.The expressions of forces and moments presenled herein are complete, relatively accurate and convenient, thus providing a good prerequisite for further studies of any problems involving the entrance of particles to a pare.
基金Supported by the National Natural Science Foundation of China under Grant No U1332125the Program for Innovation Research of Science in Harbin Institute of Technology under Grant No B201408
文摘The elliptic flow v2, for π±, K±, p and p in Au+Au collisions at center-of-mass energies √sNN=7.7, 11.5, 14.5 and 19.6 GeV, is analyzed using a multiphase transport model. A significant difference in the v2 values for p and p is observed, and the values of v2 splitting are larger compared with π+ and π-, K+ and K-. The difference increases with decreasing the center-of-mass energy. The effect of the quark coalescence mechanism in a multi-phase transport model to the value of elliptic difference △v2 between p and p- has been discussed. The simulation of Au+Au collisions at 14.5 GeV shows that the effect of hadron cascade to △v2 is not obvious, and a larger patton-scattering cross section can lead to a larger △v2.
基金supported by National Natural Science Foundation of China (Grant NOs. 41877250, 41272284)the Key Laboratory of Earth Fissures Geological Disaster, Ministry of Natural Resources (No. EFGD20240601)the Natural Science Foundation of Shaanxi Province-General Project (grant number 2023-JC-YB-231)-Suitability Evaluation of Precast Prestressed Underground Comprehensive Pipe Gallery Crossing Active Ground Fissure。
文摘The Fenwei Basin, covered by loess, experiences severe ground fissure disasters. These disasters disrupt the continuity of the loess and pose significant threats to engineering construction safety along transportation routes. Nevertheless, the crack characteristics and the influence zone of ground fissures in the loess layer remain inadequately investigated. To effectively prevent and control ground fissure disasters, physical model tests and the PFC(particle flow code) numerical simulation method are used to investigate the crack mechanism of buried ground fissures in the loess layer. The results show that there are two main cracks in the layer profile, which have a Y-shape morphology. As the dip angle of the preset cracks increased from 60° to 90°, the main deformation zone at the surface gradually shifted towards the footwall. The process of crack propagation from depth to surface is divided into five stages. Additionally, the results confirm the accuracy of the width of the rupture zone d2in the footwall calculated by the cantilever beam theory. These findings can offer theoretical guidance for determining the avoidance distance of ground fissures in loess regions, as well as for implementing disaster prevention and corresponding control measures for various stages of buried ground fissure propagation.