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.     

Microbial fuel cell (MFC)-based biosensor for combined heavy metals monitoring and associated bioelectrochemical process

To explore the feasibility and related mechanism of MFC biosensor for wastewater detection under the action of combined heavy metals. Cyclic voltammetry (CV) and scanning electrochemical impedance spectroscopy (EIS) were used to explore the related bioelectrochemical process. The response of the reactor to single/combined heavy metals, low/high heavy metal concentrations, and the differences in ohmic resistance (Rs) and charge transfer resistance (Rct) were investigated using Ni as the core heavy metal and the combined action of Cd, Cu and Zn. The results indicated that there was a linear relationship between the concentration and output voltage of the MFC biosensor under the action of combined heavy metals (R² = 0.8803–0.973). However, the internal resistance (Rint) of the MFC biosensor under the action of single heavy metal was far less than that of the combined heavy metal group, and the power density (19.849 W m^?3) was 4 times that of the combined heavy metal group (3.109–4.589 W m^?3). The Rs of the biosensors in the combined heavy metal group were 0.868Ω and 0.860, which were higher than 0.768Ω of the single heavy metal sensor. With the increase of the concentration of heavy metals in the influent, the increase of Rct was more obvious in the combined group, while the Rs in the single group significantly increased (P < 0.05). The results imply that it is possible for MFC biosensors to be used in the detection of actual water polluted by various heavy metals, but the biosensor performance is mainly limited by Rct, which needs to be further improved.     

Analysis of pressure and heat distribution in a dilating or contracting filter chamber with two outlets using multivariate spectral quasilinearization method

The development of efficient filters is an essential part of industrial machinery design, specifically to increase the lifespan of a machine. In the filter chamber design considered in this study, the magnetic material is placed along the horizontal surface of the filter chamber. The inside of the filter chamber is layered with a porous material to restrict the outflow of unwanted particles. This study aims to investigate the flow, pressure, and heat distribution in a dilating or contracting filter chamber with two outlets driven by injection through a permeable surface. The proposed model of the fluid dynamics within the filter chamber follows the conservation equations in the form of partial differential equations. The model equations are further reduced to a steady case through Lie's symmetry group of transformation. They are then solved using a multivariate spectral-based quasilinearization method on the Chebyshev–Gauss–Lobatto nodes. Insights and analyses of the thermophysical parameters that drive optimal outflow during the filtration process are provided through the graphs of the numerical solutions of the differential equations. We find, among other results, that expansion of the filter chamber leads to an overall decrease in internal pressure and an increase in heat distribution inside the filter chamber. The results also show that shrinking the filter chamber increases the internal momentum inside the filter, which leads to more outflow of filtrates.     

A numerical study on nonlinear DPL model for analyzing heat transfer in tissue during thermal therapy

In the present paper, therapeutic treatment of infected tumorous cells has been studied through mathematical modeling and simulation of heat transfer in tissues by using a nonlinear dual-phase lag bioheat transfer model with Dirichlet boundary condition. The components of volumetric heat source in this model such as blood perfusion and metabolism are assumed experimentally validated temperature-dependent function, which gives more accurate temperature distribution in tissues through this model. We have used the finite difference and RK (4, 5) techniques of numerical methods to solve the proposed problem and obtained the exact solution in a particular case. After comparison, we got a good agreement between them. We have used dimensionless quantities throughout this paper. The effect of relaxation and thermalization time with respect to dimensionless temperature distribution has been analyzed in the treatment process.     

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.     

Double pyrovanadates CaMgV2O7: Formation mechanism,phase structure,and thermodynamic properties

A double pyrovanadate CaMgV₂O₇ sample was synthesized via a facile solid-state route under an air atmosphere. The nonequilibrium formation pathways of the CaMgV₂O₇ were investigated via powder X-ray diffraction. A multistep reactions path (metavanadates–pyrovanadates–double pyrovanadate CaMgV₂O₇) was proposed to describe the formation of the CaMgV₂O₇ considering the thermodynamic and kinetic factors. The cell unit parameters of the CaMgV₂O₇ sample indicated the crystallization according to a monoclinic system with space group P12/c1(14), and the lattice parameters of a = 6.756 Å, b = 14.495 Å, c = 11.253 Å, β = 99.12, and V = 108.806 Å³. X-ray photoelectron spectroscopy also confirmed the +5 oxidation state vanadium in CaMgV₂O₇. The endothermic effects at 1033 and 1143 K were related to the incongruent melting and liquidus temperatures of CaMgV₂O₇, respectively. The comprehensive thermodynamic properties of CaMgV₂O₇ were established in both low- and high-temperature regions, utilizing a physical property measurement system and multi-high-temperature calorimetry (96 lines). The heat capacity (200 J mol K⁻¹) and entropy (198 J mol K⁻¹) at 298.15 K were computed based on the low-temperature heat capacity values, and the enthalpy of formation at 298.15 K was also estimated. The fitted high-temperature capacity can be used to obtain the changes in the enthalpy, entropy, and Gibbs free energy. This study is part of building a reliable thermodynamic database of the CaO–MgO–V₂O₅ system.     

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.