Particle flow code simulation of intact and fissured granitic rock samples

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.

Numerical simulation of predicting and reducing solid particle erosion of solid-liquid two-phase flow in a choke

Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.

Numerical simulation of flow field and residence time of nanoparticles in a 1000-ton industrial multi-jet combustion reactor

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.

DEM simulation of particle flow on a single deck banana screen

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.

Computational modeling of free-surface slurry flow problems using particle simulation method

The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.

Evaluation of plow system performance in long-wall mining method using particle flow code

The underground or open-pit methods are used for the extraction of mineral resources,each of which is divided into different categories.Coal is one of the mineral resources,which is exploited either by the surface or the underground methods.The long-wall mining is one of the methods for the underground coal mining.In this method,which is a mechanized one,some machines such as the shearer or plow are used for the mining.The coal mine in Parvadeh,Tabas is a mechanized mine that is extracted by the long-wall mining.The modeling with particle flow code software was used in this mine for the evaluation of plow performance using the coal specifications.In this regard,the sample was first calibrated by sampling from the Parvadeh coal mine and performing the uniaxial and Brazilian tests on the model.Then,the modeling was done by constructing the model and using the variables such as the clearance angle and the linear velocity of the plow.After making 28 models for the plow,the best model of the plow was selected based on the maximum force applied to the machine in the X direction.Finally,the results of this study showed that the best plow performance is for a model with the clearance angle of zero and the linear velocity of 9 mm/min,and the maximum force applied to this model is equal to 39,000 kN in the X direction.

NUMERICAL SIMULATION OF ORIENTATION DISTRIBUTION FUNCTION OF CYLINDRICAL PARTICLE SUSPENSIONS

The orientation distribution function of cylindrical particle suspensions was deduced and numerically simulated, and an application was taken in a wedge-shaped flow field. The relationship between the orientation distribution function and particle orientation angles was obtained. The results show that comparing with the most probable angle distribution which comes to being in short time, the distribution of the steady state doesn't vary much. in range; the main difference is the anti-clockwise rotation in the right and upper field, that is, particles rotate more at the points where the velocity gradients are larger. The most probable orientations are close to the direction of local streamlines. In the direction of streamlines, with poleradius decreasing, the most probable angles increase, but the angles between their orientations and the local streamlines decrease.

An improved CFD model of gas flow and particle interception in a fiber material

An improved CFD model of gas flow and particle interception in a fiber material which fiber size is Y-shape was developed in this work. The porous medium model was used to build the model of the whole size of fiber filter medium. Mixture model was adopted. The algorithm of particle interception in the whole size of fiber filter medium was derived and UDF(User Defined Function) that described kinds of particle filtering mechanisms in filter fibrous media was added to the Fluent default conservation equation as source term for simulation. The inertial resistance of the filter was taken into consideration, which provided a more precise measurement of the smoke flow and the particle interception in the filter under higher smoke speed conditions. The commercial software, Fluent 6.3, was used to simulate the smoke flow and particle interception in the filter in a single suction. The velocity and pressure profiles of smoke or nicotine particle in the filter, as well as nicotine particle volume fraction profile were well simulated. Finally, the comparisons of nicotine particle filtration efficiency between Fluent simulation results in this work and experimental results, as well as the model prediction in the literature were made to validate the simulation model. The comparisons showed that the particle entrapment model from simulation results was in good agreement with that from the experimental results. In addition, the Fluent simulation results are closer to reality both at the beginning and the end of the smoke process comparing with the model predicted results in the literature.

Numerical investigation of particle deposition on converging slot-hole film-cooled wall

Numerical research on the dilute particles movement and deposition characteristics in the vicinity of converging slot-hole(console) was carried out, and the effect of hole shape on the particle deposition characteristics was investigated. The EI-Batsh deposition model was used to predict the particle deposition characteristics. The results show that the console hole has an obvious advantage in reducing particle deposition in comparison with cylindrical hole, especially under higher blowing ratio. The coolant jet from console holes can cover the wall well. Furthermore, the rotation direction of vortices near console hole is contrary to that near cylindrical hole. For console holes, particle deposition mainly takes place in the upstream area of the holes.

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

Characteristics of turbulence transport for momentum and heat in particle-laden turbulent vertical channel flows

The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, with a focus on the contributions of particle feedback effect to momentum and heat transfer of turbulence. We take into account the effects of particles on flow drag and Nusselt number and explore the possibility of drag reduction in conjunction with heat transfer enhancement in particle-laden turbulent flows. The effects of particles on momentum and heat transfer are analyzed, and the possibility of drag reduction in conjunction with heat transfer enhancement for the prototypical case of particle-laden turbulent channel flows is addressed. We present results of turbulence modification and heat transfer in turbulent particle-laden channel flow, which shows the heat transfer reduction when large inertial particles with low specific heat capacity are added to the flow. However, we also found an enhancement of the heat transfer and a small reduction of the flow drag when particles with high specific heat capacity are involved. The present results show that particles, which are active agents, interact not only with the velocity field, but also the temperature field and can cause a dissimilarity in momentum and heat transport. This demonstrates that the possibility to increase heat transfer and suppress friction drag can be achieved with addition of particles with different thermal properties.

Numerical simulation of slurry jets using mixture model

Slurry jets in a static uniform environment were simulated with a two-phase mixture model in which flow-particle interactions were considered. A standard k-e turbulence model was chosen to close the governing equations. The computational results were in agreement with previous laboratory measurements. The characteristics of the two-phase flow field and the influences of hydraulic and geometric parameters on the distribution of the slurry jets were analyzed on the basis of the computational results. The calculated results reveal that if the initial velocity of the slurry jet is high, the jet spreads less in the radial direction. When the slurry jet is less influenced by the ambient fluid （when the Stokes number St is relatively large）, the turbulent kinetic energy k and turbulent dissipation rate e, which are relatively concentrated around the jet axis, decrease more rapidly after the slurry jet passes through the nozzle. For different values of St, the radial distributions of streamwise velocity and particle volume fraction are both self-similar and fit a Gaussian profile after the slurry jet fully develops. The decay rate of the particle velocity is lower than that of water velocity along the jet axis, and the axial distributions of the centerline particle streamwise velocity are self-similar along the jet axis. The pattern of particle dispersion depends on the Stokes number St. When St = 0.39, the panicle dispersion along the radial direction is considerable, and the relative velocity is very low due to the low dynamic response time. When St = 3.08, the dispersion of particles along the radial direction is very little, and most of the particles have high relative velocities along the streamwise direction.

Effects of silt particle on cavitation flow in centrifugal pump

Based on numerical method, effects of silt particle with certain silt mean diameter and silt concentration on the evolution of cavitation in a centrifugal pump were studied. Silt mean diameter 0.005 mm and silt concentration 1.0% were adopted in numerical simulations. Cavitation flow in a flat- nosed cylinder was simulated to validate the designed algorithm. Cavitaton flows of water and silt-laden water were simulated and compared. The results indicate that the silt particles promote the evolution of cavitation. At the outlet pressure of 6.0×10^5 Pa, cavitation bubbles do not exist in the water flow, but a few cavitation bubbles appeare in the silt-laden water flow, demonstrating the silt particles induce the formation of cavitation bubbles. At the outlet pressure of 5.29×10^5 Pa, the vapor volume fraction in the silt-laden water flow is much larger than that in the water flow, indicating that the silt particles enhance the evolution of cavitation. The properties of silt particle, static pressure, flow field structure, turbulent kinetic energy and density difference have a close relationship with the evolution of cavitation.

Effects of finite-size neutrally buoyant particles on the turbulent channel flow at a Reynolds number of 395

A direct-forcing fictitious domain(DFFD) method is used to perform fully resolved numerical simulations of turbulent channel flows laden with large neutrally buoyant particles. The effects of the particles on the turbulence(including the mean velocity,the root mean square(RMS) of the velocity fluctuation, the probability density function(PDF) of the velocity, and the vortex structures) at a friction Reynolds number of 395 are investigated. The results show that the drag-reduction effect caused by finite-size spherical particles at low particle volumes is negligibly small. The particle effects on the RMS velocities at Re_τ = 395 are significantly smaller than those at Re_τ = 180, despite qualitatively the same effects, i.e., the presence of particles decreases the maximum streamwise RMS velocity near the wall via weakening the large-scale streamwise vortices,and increases the transverse and spanwise RMS velocities in the vicinity of the wall by inducing smaller-scale vortices. The effects of the particles on the PDFs of the fluid fluctuating velocities normalized with the RMS velocities are small, regardless of the particle size, the particle volume fraction, and the Reynolds number.

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation （DNS） of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed. The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers （Pτ = 0.71, 1.5, and 3.0） and a shear Reynolds number （Reτ = 180）. Some typical thermal statistics, including normalized mean temperature and their fluctuations, turbulent heat fluxes, Nusselt number and so on, are analyzed. The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number. Two reasons can explain this. First, the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases, and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field. Second, the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number, and the thermal feedback of the particles to turbulence becomes weak.

Dynamics of self-organizing single-line particle trains in the channel flow of a power-law fluid

The formation of self-organizing single-line particle train in a channel flow of a power-law fluid is studied using the lattice Boltzmann method with power-law index 0.6≤n≤1.2,particle volume concentration 0.8%≤Φ≤6.4%,Reynolds number 10≤Re≤100,and blockage ratio 0.2≤k≤0.4.The numerical method is validated by comparing the present results with the previous ones.The effect n,Φ,Re and k on the interparticle spacing and parallelism of particle train is discussed.The results showed that the randomly distributed particles would migrate towards the vicinity of the equilibrium position and form the ordered particle train in the power-law fluid.The equilibrium position of particles is closer to the channel centerline in the shear-thickening fluid than that in the Newtonian fluid and shear-thinning fluid.The particles are not perfectly parallel in the equilibrium position,hence IH is used to describe the inclination of the line linking the equilibrium position of each particle.When self-organizing single-line particle train is formed,the particle train has a better parallelism and hence benefit for particle focusing in the shearthickening fluid at highΦ,low Re and small k.Meanwhile,the interparticle spacing is the largest and hence benefit for particle separation in the shear-thinning fluid at lowΦ,low Re and small k.

Influence of coherent structures in the gas-particle circular cylinder wake flow

To investigate the influence of coherent structures in the gas-particle wake flow, direct numerical simulation (DNS) method was adopted to compute a two-dimensional particle laden wake flow. A high accuracy spectral element method (SEM) was employed to simulate the gas flow field and a Lagrangian approach was used to compute the particles movement. Numerical results showed that at the same Stokes numbers, particles would be greatly impacted by the development of the coherent structure. But with different Stokes numbers, it can be seen that the large-scale vortex structures would influence the particle flow differently. While under different Reynolds numbers (150 and 200), there are no great changes in the particle laden flow.

Effects of Particle Concentration on the Dynamics of a Single-Channel Sewage Pump under Low-Flow-Rate Conditions

Single-channel sewage pumps are generally used to transport solid-liquid two-phase media consisting of a fluid and solid particles due to the good non-clogging property of such devices.However,the non-axisymmetric structure of the impeller of this type of pumps generally induces flow asymmetry,oscillatory outflow during operations,and hydraulic imbalance.In severe cases,these effects can jeopardize the safety and stability of the overall pump.In the present study,such a problem is investigated in the framework of a Mixture multiphase flow method coupled with a RNG turbulence model used to determine the structure of the flow field and the related motion of transported particles.It is shown that under different inlet particle concentrations,the flow field in the pump exhibits periodic variations of the pressure.The volume fraction of solid particles at the trailing edge of the suction surface of the blade is the largest,and solid particles tend to be concentrated at the outer edge of the pump body.With a rise in import particle content,the pressure and volume fraction of particles in the sewage pump also increase;for a fixed inlet particle concentration,the pressure pulsation amplitude increases with an increase in the flow rate.In addition,under small flow conditions,as the inlet particle concentration increases,the flow field leaving the sewage pump diaphragm near the outlet of the volute becomes more turbulent,and even a secondary back-flow vortex appears.

STUDY ON MESO-MECHANICAL SIMULATION OF GROUTING FOR UNFAVORABLE GEOLOGICAL CONDITIONS

The numerical simulation program of PFC2D(Particle Flow Code in 2 Dimension)particle flow based on the flow-solid coupling principle and,on its built-in FISHTANK function library and FISH language,defines the flow equation and pressure equation of fluid domain respectively,and carries out numerical simulation calculations on the diffusion process and,on the morphology and particle displacement of slurry during the slurry injection process.By adjusting the parameters of hist,n_bond,s_bond and measure in the PFC command flow,the tracking of granular body displacement changes is achieved,and the mesoscopic mechanism such as the diffusion law of soil slurry at different depths and the change of formation porosity is revealed.The numerical calculations show that:the grouting pressure has a significant effect on the alteration and destruction of the formation structure,and the fracturing effect becomes gradually worse with increasing adhesive strength,while the porosity increases significantly with increasing grouting pressure.Based on the elastic-plastic theory of the Mohr-Colomb criterion to theoretically derive the stress field of the soil around the borehole,it is pointed out that the mechanical mechanism of annular tension and radial compression is the fundamental reason for the appearance of fracturing grouting action mode.The increase of slurry viscosity is beneficial to improve the grouting effect of fracturing-compacting grouting,while the increase of friction coefficient has little effect on the grouting effect.The comparative analysis of the laboratory tests shows that the PFC2D simulation of the grouting process is feasible.

隧道围岩-衬砌接触面剪切特性模拟研究

探究隧道围岩-衬砌接触面在不同法向应力及不同粗糙度条件下的剪切特性规律,对进一步研究隧道围岩与衬砌间的相互作用及岩体稳定性具有重要的理论和工程意义。为考虑粗糙度条件,利用分形理论,通过改变分形公式中的分形维数D及分形粗糙度G生成5种不同接触面,基于二维颗粒流程序(PFC^(2D)),建立接触面数值剪切模型,尺寸为150 mm×150 mm。开展同一接触面在不同法向应力下,数值模拟剪切试验及不同接触面在同一法向应力下的数值模拟试验;将数值剪切试验与室内剪切试验获得的剪切应力-剪切位移曲线进行对比,验证数值剪切模拟的可靠性。研究结果表明:(1)数值模型剪切破坏产生的裂隙主要集中于中间接触面处,破坏后的应力主要集中于接触面粗糙度最大处,破坏裂纹可分为剪切裂纹和张拉裂纹,其中张拉裂纹占主要优势;(2)剪切应力-剪切位移过程曲线存在爬坡、啃断、滑移三大阶段,具有明显的峰值强度和残余强度,法向位移随着剪切位移的增加而增大,产生剪胀效应;(3)随着法向应力增加,剪切和张拉裂纹数量增加,破坏现象越明显,峰值强度和残余强度增大,剪胀效应减小;(4)随着粗糙起伏度增大,即分形维数D的增大或者分形粗糙度G的减小,破坏裂纹数量增加,其中张拉裂纹增加较明显,峰值强度和残余强度增大,剪胀效应增大。