Recent progress on the discrete element method simulations for powder transport systems: A review

Powder transport systems are ubiquitous in various industries, where they can encounter single powder flow, two-phase flow with solids carried by gas or liquid, and gas–solid–liquid three-phase flow. System geometry, operating conditions, and particle properties have significant impacts on the flow behavior, making it difficult to achieve good transportation of granular materials. Compared to experimental trials and theoretical studies, the numerical approach provides unparalleled advantages over the investigation and prediction of detailed flow behavior, of which the discrete element method (DEM) can precisely capture complex particle-scale information and attract a plethora of research interests. This is the first study to review recent progress in the DEM and coupled DEM with computational fluid dynamics for extensive powder transport systems, including single-particle, gas–solid/solid–liquid, and gas–solid–liquid flows. Some important aspects (i.e., powder electrification during pneumatic conveying, pipe bend erosion, non-spherical particle transport) that have not been well summarized previously are given special attention, as is the application in some new-rising fields (ocean mining, hydraulic fracturing, and gas/oil production). Studies involving important large-scale computation methods, such as the coarse grained DEM, graphical processing unit-based technique, and periodic boundary condition, are also introduced to provide insight for industrial application. This review study conducts a comprehensive survey of the DEM studies in powder transport systems.     

离子型稀土矿山淋洗及尾水富集试验研究

对浸矿后离子型稀土原地浸矿场采用清水进行淋洗，在184天的清水淋洗过程中，尾水氨氮值从最开始的507mg/L，降低至140mg/L，淋洗尾水pH4.52~3.10。淋洗尾水采用两级反渗透膜分离，既回收有价资源稀土，又能使出水氨氮达标。结果表明，产水氨氮浓度稳定低于15mg/L，对稀土的截留率高于98.25%，浓水中稀土离子平均浓度313.4mg/L，可进一步回收稀土资源。     

地应力条件下的凿岩爆破数值模拟

为了研究地应力对凿岩爆破的影响,采用DDA方法模拟爆炸应力波作用下考虑地应力条件时的单孔和多孔凿岩爆破破岩过程。模拟发现,随着初始地应力水平的增加,裂纹扩展半径和破岩区域面积减小,裂纹发育主方向趋于地应力的最大主应力方向,初始地应力对裂纹的抑制和引导作用明显;初始地应力水平的增加,对拉伸裂纹的抑制作用更为显著,从而降低了拉伸破坏对爆破破岩的贡献。模拟也表明,在初始地应力存在的条件下,通过对爆破载荷和炮孔布置进行针对性的优化,可以克服地应力带来的影响,并取得预期的爆破效果。本研究对地应力条件下的凿岩爆破工程具有理论和参考意义。     

基于有限元法的石油射孔枪弹间干扰研究

石油射孔弹顺序起爆后,弹间冲击波相互作用是一个高度非线性的力学过程。基于有限元法对这一过程进行显式动力分析,由计算结果得知,高孔密射孔枪弹间干扰产生的主要原因是下位弹在不对称的压力场起爆。设计的不同孔密和不同爆速的导爆索两种情形的爆轰试验结果表明,弹间间距越小(即孔密的增加),弹间干扰越严重;随着导爆索爆速的降低,下位弹射流头部质量集中的"中间"偏转角度增大,上位弹对其干扰越明显。以上结论可用于高孔密射孔枪的优化设计。     

Direct shear behavior of gravel-rubber mixtures: Discrete element modeling and microscopic investigations

In this paper, a newly developed 3-dimentional discrete element model (DEM) for gravel-rubber mixtures (GRMs), namely DEM4GRM, that is capable of accurately describing the macro-scale shear response (from small to large deformation) of GRMs in a direct shear box apparatus is presented. Rigid gravel grains are modelled as simple multi-shape clumps, while soft rubber particles are modeled by using deformable 35-ball body-centered-cubic clusters. Mixtures are prepared with different volumetric rubber content (VRC) at 0, 10, 25, 40 and 100%, statically compressed under 30, 60 and 100 kPa vertical stress and then sheared, by closely simulating a reference laboratory test procedure. The variation of micro-scale factors such as fabric, normal and tangential force anisotropy is carefully examined throughout the shearing process and described by means of novel micro-mechanical relationships valid for GRMs. Moreover, strong-force chains are scrutinized to identify the transition from rigid to soft granular skeleton and gain insights on the load transfer and deformation mechanisms of GRMs. It is shown that the development of the fabric and force anisotropy during shearing is closely related to the macro-scale shear strength of GRMs, and strongly depends on the VRC. Besides, strong-force chains appear to be primarily formed by gravel-gravel contacts (resulting in a rigid-like mechanical behavior) up to VRC = 30%, and by rubber-rubber contacts (causing a soft-like mechanical response) beyond VRC = 60%. Alternatively, at 30% < VRC < 60%, gravel-rubber contacts are predominant in the strong-force network and an intermediate mechanical behavior is observed. This is consistent with the behavioral trends observed in the macro- and micro-mechanical responses.     

In-situ measurement of mineral phase transition and kinetics in roasting process of Bayan Obo mixed rare earth concentrate by sodium carbonate

In this study, the Bayan Obo rare earth concentrates mixed with Na₂CO₃ were used for roasting research. The phase change process of each firing stage was analyzed. The kinetic mechanism model of the continuous heating process was calculated. This study aims to recover valuable elements and optimize the production process to provide a certain theoretical basis. Using X-ray diffraction (XRD), Fourier infrared spectroscopy, scanning electron microscopy with energy dispersive spectrometry, the reaction process and the existence of mineral phases were analyzed. The variable temperature XRD and thermogravimetric method were used to calculate the roasting kinetics. The phase transition results show that carbonate-like substances first decompose into fine mineral particles, and CaO, MgO, and SiO₂ react to form silicates, causing hardening. Further, REPO₄ and NaF can directly generate CeF₃ and CeF₄ at high temperatures, and a part of CeF₄ and NaF forms a solid solution substance Na₃CeF₇. Rare earth oxides calcined at a high temperature of 750 °C were separated to produce Ce_0.6Nd_0.4O_1.8, Ce₄O₇, and LaPrO_3+x. Then, BaSO₄, Na₂CO₃, and Fe₂O₃ react to form barium ferrite BaFe₁₂O₁₉; the kinetic calculation results show that during the continuous heating process, the apparent activation energy E reaches the minimum in the entire reaction stage in the temperature range of 440–524 °C, and the reaction order n reaches the maximum, which indicates that the decomposition product REFO significantly impacts the reaction system and reduces the activation energy. The mechanism function is F(α) = [?ln (1?α)]^1/3. The reaction order n reaches the minimum in the temperature range of 680–757 °C, and the apparent activation energy E is large. The difficulty of the reaction increases during the final stage. The reaction mechanism function is F(α) = [1?(1?α)^1/3]². Observing the entire reaction stage, the step of controlling the reaction rate changes from random nucleation to three-dimensional diffusion (spherical symmetry).     

Synthesis,microstructures, and corrosion behaviors of multi-components rare-earth silicates

Multi-components and equimolar rare earth monosilicates, (Y_1/3Dy_1/3Er_1/3)₂SiO₅, (Y_1/3Dy_1/3Lu_1/3)₂SiO₅, (Y_1/4Dy_1/4Ho_1/4Er_1/4)₂SiO₅ and (Yb_1/4Dy_1/4Ho_1/4Er_1/4)₂SiO₅, were prepared by solid-state reactions and the following hot-pressing. Dense microstructures with uniform elemental distributions were obtained for all samples. These investigated multi-components monosilicates exhibit low thermal conductivities and similar coefficients of thermal expansion with SiC. Moreover, they exhibit high corrosion resistances in 1400 °C water vapor, especially, four-components (Y_1/4Dy_1/4Ho_1/4Er_1/4)₂SiO₅ and (Yb_1/4Dy_1/4Ho_1/4Er_1/4)₂SiO₅ experienced almost invariable weights after small weight losses during the initial 0.5 h. All those results indicate that multi-components rare earth monosilicates are promising candidates of environmental barrier coatings for SiC/SiC composites.