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.     

An analytical model for gas leakage through contact interface in proton exchange membrane fuel cells

Sealing performance between two contacting surfaces is of significant importance to stable operation of proton exchange membrane (PEM) fuel cells. In this work, an analytical micro-scale approach is first established to predict the gas leakage in fuel cells. Gas pressure and uneven pressure distribution at the interface are also included in the model. At first, the micro tortuous leakage path at the interface is constructed by introducing contact modelling and fractal porous structure theory. In order to obtain the leakage at the entire surface, contact pressure distribution is predicted based on bonded elastic layer model. The gas leakage through the discontinuous interface can be obtained with consideration of convection and diffusion. Then, experiments are conducted to validate the numerical model, and good agreement is obtained between them. Finally, influences of surface topology, gasket compression and gasket width on leakage are studied based on the model. The results show that gas leakage would be greatly amplified when the asperity standard deviation of surface roughness exceeds 1.0 μm. Gaskets with larger width and smaller thickness are beneficial to sealing performance. The model is helpful to understand the gas leakage behavior at the interface and guide the gasket design of fuel cells.     

超临界机组发生FAC的因素及相互关系分析

为了减轻因流动加速腐蚀(FAC)引起的锅炉结垢加速、汽水系统管道厚度减小甚至爆裂现象，对超临界机组发生流动加速腐蚀的机理及其主要影响因素进行了研究，并讨论了管壁内表面粗糙度、蒸汽含汽率、pH值、溶氧量对FAC的影响，以及温度与pH值、温度与流速、pH值与溶解氧量、溶解氧量与氢电导率等影响因素之间的相互作用关系，最后结合实际电厂的运行数据验证了分析结果。研究表明：减小工质流速、管壁粗糙度和氢电导率，增大给水的pH值和溶解氧含量可以使FAC的腐蚀速率减小，超临界加氧处理时pH值应在8.9～9.2之间，溶解氧量范围为45～100μg/L,氢电导率的期望值在0.1μS/cm以下。由于各影响因素之间的作用十分复杂，本文只给出了大致范围和趋势，并未给出准确数据。     

Hall current effects on MHD flow of chemically reacting fluid through a porous medium with heat source

We considered the magnetohydrodynamic (MHD) free convective flow of an incompressible electrically conducting viscous fluid past an infinite vertical permeable porous plate with a uniform transverse magnetic field, heat source and chemical reaction in a rotating frame taking Hall current effects into account. The momentum equations for the fluid flow during absorbent medium are controlled by the Brinkman model. Through the undisturbed state, both the plate and fluid are in a rigid body rotation by the uniform angular velocity perpendicular to an infinite vertical plate. The perpendicular surface is subject to the homogeneous invariable suction at a right angle to it and the heat on the surface varies about a non-zero unvarying average whereas the warmth of complimentary flow is invariable. The systematic solutions of the velocity, temperature, and concentration distributions are acquired systematically by utilizing the perturbation method. The velocity expressions consist of steady-state and fluctuating situations. It is revealed that the steady part of the velocity field has a three-layer characteristic while the oscillatory part of the fluid field exhibits a multi-layer characteristic. The influence of various governing flow parameters on the velocity, temperature, and concentration are analyzed graphically. We also discuss computational results for the skin friction, Nusselt number, and Sherwood number in the tabular forms.     

南水北调中线总干渠藻类的生态调度

南水北调中线总干渠无在线调蓄水库,对藻类生态调度过程中出现的问题开展生态调度实现策略和实施方式研究。主要实现策略包括:划定自身的调蓄区,隔离生态调度对下游的影响;采用高效的渠池运行方式,减少生态调度时蓄量的反复调整;综合考虑安全、快速、平稳等需求,设定生态调度实施进程和方式。具体实施方式包括:将总干渠划分为流速调控区、调蓄区和正常运行区,分别实施等体积、控制蓄量和闸前常水位方式运行;将生态调度过程划分为充水阶段和泄水阶段,基于流速调控目标值、持续时长和水位降幅约束条件,确定各阶段时长和各分区的闸门群调控方案等。基于2018年3月输水工况,采用明渠一维非恒定流模型,仿真总干渠上游15个渠池的藻类生态调度过程。结果表明,生态调度可在3.5 d内完成,各渠池的平均流速由0.48 m/s增至0.93 m/s,持续时间超过2 h。在整个生态调度过程中,水位变化平稳,水位变幅符合安全阈值要求,下游渠道的正常运行未受生态调度明显影响。     

基于移动窗口法的城市遥感生态质量评价研究

为进一步提高城市生态监测的精确度，在前人开展的城市生态遥感监测研究的基础上引入“区域尺度”的概念，充分考虑小区域范围内不同地物间的交互作用，并利用移动窗口模型(Moving Window-Remote Sensing Ecology Index, MW-RSEI)对沈阳市浑南区的遥感影像进行逐像元分析。实验结果表明，MW-RSEI模型和生态遥感监测模型（Remote Sensing Ecology Index, RSEI)在整体生态评价趋势上表现出一致性。但MW-RSEI模型对城市生态中的细节区域表征更明显。考虑到建筑及裸地周围植被的影响，其较差生态区域的占比仅为6%，生态评价为优的区域占比为11%。研究区生态评价结果表现得较为连续，并未出现明显的断层现象。该结果具有现实意义且与实际生态分布更为符合。MW-RSEI模型在城市生态监测中更为适用，可为相关部门提供一定的技术参考。     

Shock/expansion fan interaction with real gas effects for Earth and Mars atmospheric conditions

Numerical simulations are performed to investigate the real gas effects on shock/expansion fan interaction. Initial perfect gas simulations at low enthalpy capture the flow structures efficiently and outcomes are found to have excellent agreement with the analytical calculations. Furthermore, the simulations with the real gas solver for different enthalpies showed that the variation in enthalpy significantly changes the flow structures. It is observed that an increase in enthalpy leads to a decrease and increase in the postshock and postexpansion fan Mach numbers, respectively. Another important observation is the decrement in the peak pressure ratio with an increment in the enthalpy. These effects are noted to be more pronounced for Mars's environment due to the higher dependency of specific heat on temperature.     

Nonsimilar solution of a boundary layer flow of a Reiner–Philippoff fluid with nonlinear thermal convection

The thermodynamics modeling of a Reiner–Philippoff-type fluid is essential because it is a complex fluid with three distinct probable modifications. This fluid model can be modified to describe a shear-thinning, Newtonian, or shear-thickening fluid under varied viscoelastic conditions. This study constructs a mathematical model that describes a boundary layer flow of a Reiner–Philippoff fluid with nonlinear radiative heat flux and temperature- and concentration-induced buoyancy force. The dynamical model follows the usual conservation laws and is reduced through a nonsimilar group of transformations. The resulting equations are solved using a spectral-based local linearization method, and the accuracy of the numerical results is validated through the grid dependence and convergence tests. Detailed analyses of the effects of specific thermophysical parameters are presented through tables and graphs. The study reveals, among other results, that the buoyancy force, solute and thermal expansion coefficients, and thermal radiation increase the overall wall drag, heat, and mass fluxes. Furthermore, the study shows that amplifying the space and temperature-dependent heat source parameters allows fluid particles to lose their cohesive force and, consequently, maximize flow and heat transfer.     

Impact of viscous dissipation and Dufour on MHD natural convective flow past an accelerated vertical plate with Hall current

The present research work concentrates on viscous dissipation, Dufour, and heat source on an unsteady magnetohydrodynamics natural convective flow of a viscous, incompressible, and electrically conducting fluid past an exponentially accelerated infinite vertical plate in the existence of a strong magnetic field. The presence of the Hall current induces a secondary flow in the problem. The distinguishing features of viscous dissipation and heat flux produced due to gradient of concentration included in the model along with heat source as they are known to arise in thermal-magnetic polymeric processing. The flow equations are discretized implicitly using the finite difference method and solved using MATLAB fsolve routine. Numerical values of the primary and secondary velocities, temperature, concentration, skin friction, Nusselt number, and Sherwood number are illustrated and presented via graphs and tables for various pertinent parametric values. The Dufour effect was observed to strengthen the velocity and temperature profile in the flow domain. In contrast, due to the impact of viscous dissipation, the local Nusselt number reduces. The study also reveals that the inclusion of the chemical reaction term augments the mass transfer rate and diminishes the heat transfer rate at the plate.     

Investigation of liquid water heterogeneities in large area proton exchange membrane fuel cells using a Darcy two-phase flow model in a multiphysics code

Proper management of the liquid water and heat produced in proton exchange membrane (PEM) fuel cells remains crucial to increase both its performance and durability. In this study, a two-phase flow and multicomponent model, called two-fluid model, is developed in the commercial COMSOL Multiphysics® software to investigate the liquid water heterogeneities in large area PEM fuel cells, considering the real flow fields in the bipolar plate. A macroscopic pseudo-3D multi-layers approach has been chosen and generalized Darcy's relation is used both in the membrane-electrode assembly (MEA) and in the channel. The model considers two-phase flow and gas convection and diffusion coupled with electrochemistry and water transport through the membrane. The numerical results are compared to one-fluid model results and liquid water measurements obtained by neutron imaging for several operating conditions. Finally, according to the good agreement between the two-fluid and experimentation results, the numerical water distribution is examined in each component of the cell, exhibiting very heterogeneous water thickness over the cell surface.