Filtration characteristics of a binary multi-layer granular bed filter based on CFD-DEM coupling simulation

The coupled CFD-DEM method with the JKR (Johnson-Kendall-Roberts) model for describing the contact adhesion of dust to filter particles (FPs) is used to simulate the distribution pattern of dust particle deposition in the granular bed filter (GBF) with multi-layer media. The minimum inlet flow velocity must meet the requirement that the contact probability between dust and FPs is in the high contact probability region. The air flow forms vortices on the leeward side of the FPs and changes abruptly at the intersection of different particle size FPs layers. Dust particles form large deposits at the intersection of the first and second layers and the different particle size filter layers. Dual element multilayer GBF can further optimize the bed structure by interlacing filter layers with different particle sizes. Compared with single particle size multi-layer GBF, the bed pressure drop is reduced by 40.24%–50.65% and the dust removal efficiency is increased by 21.93%–55.09%.     

A two‐grid fictitious domain method for direct simulation of flows involving non‐interacting particles of a very small size

The full resolution of flows involving particles whose scale is hundreds or thousands of times smaller than the size of the flow domain is a challenging problem. A naive approach would require a tremendous number of degrees of freedom in order to bridge the gap between the two spatial scales involved. The approach used in the present study employs two grids whose grid size fits the two different scales involved, one of them (the micro‐scale grid) being embedded into the other (the macro‐scale grid). Then resolving first the larger scale on the macro‐scale grid, we transfer the so obtained data to the boundary of the micro‐scale grid and solve the smaller size problem. Since the particle is moving throughout the macro‐scale domain, the micro‐scale grid is fixed at the centroid of the moving particle and therefore moves with it. In this study we combine such an approach with a fictitious domain formulation of the problem resulting in a very efficient algorithm that is also easy to implement in an existing CFD code. We validate the method against existing experimental data for a sedimenting sphere, as well as analytical results for motion of an inertia‐less ellipsoid in a shear flow. Finally, we apply the method to the flow of a high aspect ratio ellipsoid in a model of a human lung airway bifurcation. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental study on filtration performance of the moving bed granular filter with axial flow

The filtration performance of the moving bed granular filter with axial flow (MBGF-AF) is investigated through a large cold experiment. The effect of different operation parameters on the filtration performance (collection efficiency, pressure drop) of the axial-flow moving bed filter is investigated in combination with the dust deposition effect and the mechanism of trapping dust by the capturing particles. The results show that the collection efficiency of MBGF-AF is enhanced by decreasing the superficial gas velocity, increasing the inlet dust concentration properly, or decreasing the moving velocity of the capturing particles. A model covering the above operation parameters is established to calculate the collection efficiency of the moving bed granular filter. It is used in a wide range of operating parameters for the MBGFs.     

Construction of irregular particles with superquadric equation in DEM

Non-spherical particles are widely present in industrial production, and significantly affect the macro and micro characteristics of granular materials. Although the superquadric equation can be used to construct non-spherical particles, its disadvantage is that the particle shape is geometrically symmetric and strictly convex. In this study, two composed approaches are used to describe geometrically asymmetric and concave particle shapes, including a multi-superquadric model and a poly-superquadric model. The multi-superquadric model is a combination of several superquadric elements, and can construct concave and geometrically asymmetric particle shapes.The poly-superquadric model is a combination of eight one-eighth superquadric elements, and can construct convex and geometrically asymmetric particle shapes. Both composed models are based on superquadric equations, and Newton's iterative method is used to calculate the contact force between the elements. Furthermore, superquadric elements, multi-superquadric elements,and poly-superquadric elements are applied for the formation of complex granular beds, and the influences of particle shape on the packing fraction can be successfully captured by the proposed models.     

基于细观方向平均模型的颗粒材料宏观Cosserat连续体本构关系

提出了基于细观微-方向模型（Micro—Directional Model）的宏观Cosserat连续体本构关系。在细观尺度上考虑颗粒旋转自由度及接触力矩,微结构的影响通过接触分布函数体现。给出均质各向同性Cosserat连续体模型弹性常数的细观参数表达式,并建议了二维情况下内尺度参数的细观力学表达式。对颗粒材料宏观行...     

A dual-scale coupled micro/macro segregation model for dendritic alloy solidification

The present article reports on the formulation, numerical implementation, and application of a single-domain coupled micro/macroscopic model for simulation of dendritic alloy solidification. Microscopic solutal non-equilibrium effects have been included in the macroscopic modeling of solidification by using a fixed grid dual scale numerical approach. Salient features of the present approach include a continuum model for conservation of mass, momentum, energy and species on the macroscopic scale, a microscopic solute redistribution model, and the solution procedure and auxiliary equations necessary for coupling between the two models. The coupling between macro and micro scale models is practically made possible by introducing an iterative micro/macro time step scheme. The local solidification rate is calculated by implicit iterations of macroscopic conservation equations and the microscopic solute redistribution model. The present model is capable to simulate eutectic reaction, local re-melting, and account for the inter-linkage between micro and macroscopic solute redistribution (micro and macrosegregation).     

Macro–micro relations in granular mechanics

In granular mechanics, macroscopic approaches treat a granular material as an equivalent continuum at macro-scale, and study its constitutive relationship between macro-quantities, such as stresses and strains. On the other hand, microscopic approaches consider a granular material as an assembly of individual particles interacting with each other at micro-scale (i.e., particle-scale), and the physical quantities under study are forces and displacements. This paper aims at linking up the findings from these two scales and to establish the macro–micro relations in granular mechanics.Three aspects of the macro–micro relations are investigated. They are about the internal structure, the stress tensor and the strain tensor. The internal structure is described with geometrical systems at the particle scale. Micro-structural definitions of the stress and strain tensors are derived, which link the macro-stress tensor with the contact forces, and the macro-strain tensor with the relative displacements at contact. In addition to a brief review of the past research work on these topics, further generalizations are made in this paper. In particular, the two cell systems proposed by Li and Li (2009), namely the solid cell system and the void cell system, are introduced and used for the derivation of the macro-structural expressions. The stress expression is derived based on Newton’s second law of motion. The result is valid for both static and dynamic cases. The strain expression is derived based on the compatibility requirement. And the expression is valid for any tessellation subdividing the granular assembly into polyhedral elements.The homogenization for deriving a macroscopic constitutive relationship from microscopic behaviour is discussed. Attention is placed on the macroscopic quantification of the internal structure in terms of a second rank tensor, known as the fabric tensor. Existing definitions of the fabric tensors have been reviewed. The correlations among different fabric tensors and their relations with the stress–strain behaviours have been investigated.     

Simulated configurational temperature of particles and a model of constitutive relations of rapid-intermediate-dense granular flow based on generalized granular temperature

In gas–solid flat-base spout bed with a jet, the flow of particles must go through an intermediate regime where both kinetic/collisional and frictional contributions play a role. In this paper, the statistical framework is proposed to define the generalized granular temperature which sums up the configurational temperature and translational granular temperature. The configurational temperature, translational and rotational granular temperatures of particles are simulated by means of CFD-DEM (discrete element method) in a 3D flat-base spout bed with a jet. The configurational temperatures of particles are calculated from instantaneous overlaps of particles. The translational and rotational granular temperatures of particles are calculated from instantaneous translational and angular velocities of particles. Roughly, the simulated translational and rotational granular temperatures increase, reach maximum, and then decrease with the increase of solids volume fractions. However, the configurational temperature increases with the increase of solids volume fractions. At high solid volume fraction, the predicted configurational temperatures are larger than the translational and rotational granular temperatures, indicating that the rate of energy dissipation do contributes by contact deformation of elastic particles. The generalized granular temperature is proposed to show the relation between the variance of the fluctuation velocity of deformation and the variance of the translational fluctuation velocity of particles. The constitutive relations of particle pressure, viscosity, granular conductivity of fluctuating energy and energy dissipation in rapid-intermediate-dense granular flows are correlated to the generalized granular temperature. The variations of particle pressure, shear viscosity, energy dissipation and granular conductivity are analyzed on the basis of generalized granular temperature in a flat-base spout bed with a jet. The axial velocities of particles predicted by a gas–solid two-fluid model of rapid-intermediate-dense granular flows agree with experimental results in a spout bed.     

The effect of solutal undercooling on double‐diffusive convection and macrosegregation during binary alloy solidification: a numerical investigation

In this paper, we present a macroscopic numerical model that is capable of capturing the interaction between the double‐diffusive convective field and a localized fluid flow on account of solutal undercooling during non‐equilibrium solidification of binary alloys. The model is essentially based on a fixed‐grid enthalpy based control volume approach. In the present model, microscopic features pertaining to non‐equilibrium effects on account of solutal undercooling are incorporated through the formulation of a modified partition‐coefficient. The effective partition‐coefficient is numerically modelled by means of a number of macroscopically observable parameters related to the solidifying domain. This feature has made the present treatment different from micro‐macro modelling of alloy solidification, which involves certain parameters that may not be macroscopically resolvable. Numerical simulations are performed for the case of two‐dimensional transient solidification of Pb–Sn alloys (both hypoeutectic and hypereutectic) in a rectangular cavity, employing the present model. The simulation results are also compared with the corresponding experimental results quoted in the literature, and the agreement is excellent. From the results, it can be concluded that non‐equilibrium effects on account of solutal undercooling result in a more enhanced macrosegregation. Copyright © 2002 John Wiley & Sons, Ltd.     

基于交变电场的月表除尘方法研究进展

实施月球表面探测时, 月表尘埃物质容易在探测设备表面沉积并造成严重不利影响. 近年来交变电场 除尘方法被认为是在月表进行尘埃防护的有效手段, 但要使该方法实用化, 仍需要弄清除尘的机理, 确定相关 影响的因素, 以便进行优化设计. 本文从实现交变电场的电帘出发, 分别介绍电帘表面的交变电场分布、颗粒 在电帘表面的受力状态、颗粒群运动过程中的摩擦电荷分布、颗粒运动的模拟方法、颗粒运动过程、除尘效 率和能耗的理论和实验研究, 并对自清洁太阳能电池板系统的设计与制备进行了总结. 对颗粒的起跳过程、静 电的生成和除尘效率影响因素等普遍关注的问题进行了细致的讨论和总结, 指出了颗粒的带电量对起跳过程 的关键作用. 最后, 对该除尘防尘存在的问题进行了归纳, 并展望了该领域的发展趋势.      

Dust distribution of solid and adhesive mixed dust in a granular bed filter

Granular bed filters can effectively filter adhesive dust in high-temperature flue gas. In this study, polyvinyl chloride (PVC) powder was used as adhesive dust, and the mixture of PVC and ash powder was used to simulate solid and adhesive mixed dust. The effects of gas temperature, velocity, and inlet adhesive dust mass content on dust distribution in granular bed (GBF) were discussed. Results show that the mixed dust mainly accumulates on the upper part of the granular bed, and the mass of the collected dust decreases exponentially from the upper layer to the bottom layer in the GBF. The adhesive dust content collected in each layer differs from that of the incoming dust, and their deviation varies approximately linearly along with the depth of the bed. The total dust distribution and adhesive dust content deviation are influenced by gas temperature and inlet adhesive dust content but independent of gas velocity. The correlations of dust distribution of solid and adhesive mixed dust are presented based on the experimental results.     

Effects of grid size on predictions of bed expansion in bubbling fluidized beds of Geldart B particles: A generalized rule for a grid-independent solution of TFM simulations

Numerical simulations of gas–solid fluidized beds based on the kinetic theory of granular flow exhibit a significant dependence on domain discretization. Bubble formation, bubble size and shape all vary greatly with the discretization, and the use of an inappropriate scale resolution leads to inaccurate predictions of fluidization hydrodynamics. In this study, grid-independent solutions of the two fluid model were examined by comparing the bed expansions obtained from numerical simulations with experimental results and empirical predictions, based on bubbling fluidized beds of Geldart B particles. Grid independence was achieved with a grid resolution equal to 18 times the particle diameter. The simulation results were compared with previously published data for verification purposes. The results of this work should provide a guideline for choosing the appropriate grid size and thereby minimize the time and expense associated with large simulations.     

Inertial steady flow of a bed of granular material down a surface with microrelief with allowance for finite time of intergranular contacting

An inertial flow of a granular material can be described by the laws of conservation of mass, momentum, and energy of random motion of solid particles by invoking some closing relations. In this work, these closing relations are inferred from the dimensional theory. The system of equations obtained is used to determine characteristics of a steady flow of a bed of a granular material down an inclined surface with a microrelief for various Richardson numbers and finite contact times of the particles during their collisions. Novosibirsk Military Institute, Novosibirsk 630103. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 6, pp. 128–132, November–December, 1999.     

A model for expansion ratio in liquid-solid fluidized beds

Liquid-solid fluidized beds are used in mineral processing industries to separate particles based on parti- cle size, density, and shape. Understanding the expanded fluidized bed is vital for accurately assessing its performance. Expansion characteristics of the fluidized bed were studied by performing several experi- ments with iron ore, chromite, quartz, and coal samples. Using water as liquid medium, experiments were conducted to study the effects of particle size, particle density, and superficial velocity on fluidized bed expansion. The experimental data were utilized to develop an empirical mathematical model based on dimensional analysis to estimate the expansion ratio of the fluidized bed in terms of particle character- istics, operating and design parameters. The predicted expansion ratio obtained from the mathematical model is in good agreement with the experimental data.     

Penetration into a granular bed in the presence of upward gas flows

We present a numerical study on the penetration of spherical projectiles into a granular bed in the presence of upward gas flows. Due to the presence of interstitial fluid, the force chains between particles in the granular bed are weakened significantly, and this distinguishes the penetration behavior from that in the absence of fluid. An interesting phenomenon, namely granular jet, is observed during the penetration, and the mechanism for its formation and growth is attributed to the merging of granular vortices generated by the interaction between the intruder and primary particles. Moreover, both the final penetration depth and the maximum diameter of the crater are found to follow a power-law dependence with the impact velocity, and the maximum height reached by the granular jet tends to increase linearly as the impact velocity increases, agreeing well with the experimental results reported in the literature.     

Chebyshev super spectral viscosity solution of a two‐dimensional fluidized‐bed model

The numerical solution of a model describing a two‐dimensional fluidized bed by a Chebyshev super spectral viscosity (SSV) method is considered. The model is in the form of a hyperbolic system of conservation laws with a source term, coupled with an elliptic equation for determining a stream function. The coupled elliptic equation is solved by a finite‐difference method. The mixed SSV/finite‐difference method produces physically shaped bubbles, on a very coarse grid. Fine scale details, which were not present in previous finite‐difference solutions, are present in the solution. Copyright © 2003 John Wiley & Sons, Ltd.     

Coupled Navier–Stokes—Molecular dynamics simulations using a multi‐physics flow simulation framework

Simulation of nano‐scale channel flows using a coupled Navier–Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi‐physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub‐domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier–Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier–Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro‐scale MD problems are constructed to provide viscous stresses for the Navier–Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite‐volume discretization are evaluated using MD. The second example deals with fully developed three‐dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier–Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite‐volume discretization of the macro‐domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro‐scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd.     

Solving the Savage–Hutter equations for granular avalanche flows with a second‐order Godunov type method on GPU

In this article, we apply Davis's second‐order predictor‐corrector Godunov type method to numerical solution of the Savage–Hutter equations for modeling granular avalanche flows. The method uses monotone upstream‐centered schemes for conservation laws (MUSCL) reconstruction for conservative variables and Harten–Lax–van Leer contact (HLLC) scheme for numerical fluxes. Static resistance conditions and stopping criteria are incorporated into the algorithm. The computation is implemented on graphics processing unit (GPU) by using compute unified device architecture programming model. A practice of allocating memory for two‐dimensional array in GPU is given and computational efficiency of two‐dimensional memory allocation is compared with one‐dimensional memory allocation. The effectiveness of the present simulation model is verified through several typical numerical examples. Numerical tests show that significant speedups of the GPU program over the CPU serial version can be obtained, and Davis's method in conjunction with MUSCL and HLLC schemes is accurate and robust for simulating granular avalanche flows with shock waves. As an application example, a case with a teardrop‐shaped hydraulic jump in Johnson and Gray's granular jet experiment is reproduced by using specific friction coefficients given in the literature. Copyright © 2014 John Wiley & Sons, Ltd.