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.     

一种“KooN”结构的马尔可夫简并建模和简化计算的方法

IEC61508等标准提出了几种计算安全仪表系统需求时平均失效概率的方法,但是,对于具有冗余配置的复杂系统,随着组件数量的增加,系统的中间状态数量快速增长,用户难以构建马尔可夫模型,即便借助计算机来建模运算也较为耗时。提出了一种同型“K oo N”简并状态的马尔可夫建模的通用方法,首先是根据降级状态进行判断,将符合条件的状态进行简并,然后对标记为危险失效状态的概率进行计算。通过严格的理论推导,该简并状态方法可以在不损失精度的前提下简化马尔可夫建模。     

Combined effect of channel to rib width ratio and gas diffusion layer deformation on high temperature – Polymer electrolyte membrane fuel cell performance

The present study investigates the combined influence of Channel to Rib Width (CRW) ratio and clamping pressure on the structure and performance of High Temperature-Polymer Electrolyte Membrane Fuel Cell (HT-PEMFC) using a three-dimensional numerical model developed previously. It also considers the impact of interfacial contact resistance between the Gas Diffusion Layer (GDL) and Bipolar Plate (BPP). The structural analysis of the single straight channel HT-PEMFC geometry shows that the von-Mises stress greatly increases in the GDL under the ribs as the CRW ratio increases resulting in considerably high deformation. The cell performance analysis depicts the significance of ohmic resistance and concentration polarization for different CRW ratios, particularly at higher operating current densities. However, in low to medium current density regions, the CRW ratio has little influence on cell performance. A substantial impact on the species, overpotential, and current distributions is observed. The findings also reveal that the CRW ratio significantly affects the temperature distribution in the cell.     

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

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

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.     

Acidification of the electrolyte during the galvanic corrosion of AA7075: A numerical and experimental study

The electrochemical interactions between aluminum alloy 7075 and low-carbon steels under gelled electrolytes were studied. Such electrolytes provided the opportunity to investigate both thick and thin electrolyte systems. The electrolyte was chemically modified to visually track the acidic fronts during the anodic reaction and the subsequent hydrolysis process. Two mathematical models were validated for both thick and ultrathin electrolytes. The acidification of thick electrolytes was extended some millimeters beyond the aluminum alloy surface, whereas the acidic front was localized next to the metallic joint using ultrathin electrolytes. The combination of both numerical and experimental results allows proving (and explaining why) that the acidification process is more aggressive under dilute than under concentrated electrolytes.     

Form selection of concomitant polymorphs: A case study informed by crystallization kinetics modeling

Molecular mechanisms and process kinetics of crystallizing concomitant polymorphs remain poorly understood. Solvent-mediated phase transformation and concomitant crystallization are difficult to be distinguished in practice, as multiple forms can be detected at the same time. Herein, we developed a population balance model to simulate a concomitant crystallization process of two polymorphs of tolfenamic acid. Our kinetic modeling aims to understand concomitant crystallization and help guide form selection of such a molecular system. Crystallization kinetics of ethanolic solutions were uncovered from induction time measurements, as well as seeded and unseeded crystallization experiments. Experimental and simulation results demonstrate that the stable form I crystallizes concomitantly with the metastable form II. The faster growing form II results in an intermediate decline in the composition of form I in crystallized samples, a characteristic feature of the concomitantly crystallized system. A four-quadrant scheme of attainable polymorph outcome was simulated under various crystallization conditions.     

Evidence of fluid evolution of Baoshan Cu−Pb−Zn polymetallic deposit: Constraints from in-situ sulfur isotope and trace element compositions of pyrite

In-situ LA-ICP-MS and S isotopes of pyrite from the Baoshan Cu polymetallic deposit were conducted to investigate the ore-forming process and the enrichment mechanism of elements. Three generations of pyrite (Py I, Py II, and Py III) in the skarn-type ores and pyrite in the carbonate-hosted sulfide ores from central, western, and northern (C_Py, W_Py, and N_Py) mining districts are selected for comparison. Compared with Py I and Py III, the contents of most elements in Py II are apparently higher. The As and Se contents are high within a wide range and are decoupled in the growth band of the C_Py. The highest As, Se, and Pb contents were found in W_Py and N_Py. These results indicate the drastic changes in the temperature and fluid mixing during the mineralization. The occurrence of fluctuation and change in temperature and f(O₂) was triggered by intermittent pulses of magmatic-hydrothermal fluids, mixing with meteoric water, and water−rock interactions. The sulfur isotopes of all species of pyrite indicated the magmatic source. The change in the f(O₂) conditions caused slight differences in the sulfur isotope compositions. Consequently, a metallogenic model was proposed to explain the ore-forming processes.