Numerical investigation on the mixing mechanism in a cross-torus paddle mixer using the DEM-CFD method
Graphical abstract
Introduction
Powder and particulate materials have been widely utilized in a variety of engineering fields such as food, pharmaceutical, chemical, mineral, and metallurgical engineering, in which the mixing and blending of powders is one of the most demanding unit operations to ensure the quality and performance of industrial products [[1], [2], [3], [4]]. Various types of industrial mixer are currently available with the purpose of achieving an efficient mixing under a specific industrial condition [5]. For instance, according to the features of mixing-driven structure, the industrial mixers can be broadly classified into mixers with rotating vessels (e.g., tumbling drum, V blender, double-cone blender, batch mixer etc.), and mixers with fixed vessels and rotating components (e.g., ribbon blender, paddle blender, ploughshare mixer etc.) [3]. According to the different predominant roles in mixing mechanisms, mixers can also be regrouped into convective, shear, and diffusive dominant mixers [6,7].
To analytically evaluate the granular dynamics, mixing, and segregation within various mixing unit operations, a series of experimental approaches are proposed such as particle sampling, visual tracking, particle image velocimetry, positron emission particle tracking, and magnetic resonance imaging [2,8]. However, those experimental tools may be expensive and intractable to obtain all the detailed information of particles within the powder mixing process [4,9,10]. Incidentally, even if the experiments can be performed successfully, development of theoretical models is impossible to predict mixing state in a mixer. This is because the theoretical solution does not exist in the powder system. In contrast, due to the improvement in advanced numerical models and computer hardware, numerical simulation provides us with a comparably more feasible way to perform the parametric studies and allow us to attain some insights into the optimization of device design and operating conditions [11]. The discrete element method (DEM), originally proposed by Cundall and Strack [12], has been widely applied to the modeling of granular and multiphase flows within a variety of industrial processes [13], such as pneumatic conveying [14,15], screw conveying [16,17], die filling [[18], [19], [20]], fluidized bed [11,[21], [22], [23], [24], [25], [26], [27]], and blending [[28], [29], [30], [31], [32], [33]]. Moreover, with an aim to deal with the modeling of complex shape boundary involved in a dense granular flow especially for the difficult procedures of particle-to-wall collisions, the DEM simulation coupled with the signed distance function (SDF) is currently developed [17], which has been validated in previous studies with a high adequacy as compared to the experiments [4,9,28]. As far as the modeling of arbitrary-shaped boundaries using the computational fluid dynamics (CFD), the immersed boundary method (IBM) has been recently introduced and developed [[34], [35], [36], [37], [38]], in which an arbitrary-shaped boundary can be easily modeled using the local-volume fraction of a solid object. A combination of the SDF and the IBM, referred to the Advanced DEM-CFD method [18,39], herein provides us with a desirable and adequate approach with a broad applicability to a variety of real industrial systems [40].
Over the past decades, to acquire a comprehensive understanding on the mixing features and mechanisms in industrial mixers [41], a large amount of research works including both experiments and simulations have been carried out. Typically, combined with operational conditions, particle properties, and mixer types, the mixing dynamics can exhibit different predominant modes along with all of three mechanisms (i.e., convective, diffusive, and shear mixing) to a higher or lower extent [42,43], which is accompanied by other abundant behaviors in terms of various flow patterns, segregations, avalanches etc. For instance, in a batch mixer such as in a rotating drum, different flow regimes (i.e., slipping, slumping, rolling, cascading, cataracting, and centrifuging) have been experimentally and numerically identified [44,45]. Moreover, for the mixing mechanism in such a mixer, the radial mixing is much stronger than axial mixing, thus, the predominant mixing mode is diffusion [9,43]. In a ribbon mixer, evidence shows that mixing in the axial direction is not as good as mixing in the perpendicular direction; thus, the predominant mixing mode is convection [4,46]. In a horizontal or vertical shaft mixer, it is verified that a high shearing force is exerted by the shear motion of mechanical blades to promise a high-efficient mixing process; hence, the predominant mode is shear [47]. Unfortunately, it should be pointed out that, despite of some achievements on the mechanism analysis of geometry-related properties as indicated above (e.g., batch mixer, ribbon mixer, and shaft mixer), an unambiguously scientific understanding especially for the effects of operational conditions remains limited. Consequently, the design and selection of an appropriate mixing equipment that is suitable for current industrial environments are still not straightforward [3]. Additionally, the mixing mechanism over a broad scope of operational and geometrical parameters has not yet been completely investigated and understood. Considering the critical importance of mixing mechanisms to determine the appropriate parameters to improve the mixing efficiency, it is necessary to provide a novel insight into the predominant mode related to the mixing performance as the change of various operating parameters (e.g., the effect of liquid viscosity, filling level etc.).
Focusing on the above aspects, in the present study, a series of simulations for powder mixing in an elaborately designed mixer, called a cross-torus paddle mixer, is conducted using the Advanced DEM-CFD method. The mixing vessel composed of a hemisphere and a coaxial cylinder remains stationary during the mixing process, inside which a torus paddle attached to a cross support can rotate along one extended rod of the cross support. In this mixer, optimization design, dominant parameters for mixing and mixing mechanism have not been investigated thus far. Hence, the paddle shape and mixing conditions are empirically decided by engineers. Moreover, in such a cross-torus paddle mixer, the mixing state is empirically known to be influenced by rotational speed, filing level, liquid viscosity, and paddle obliquity. Indeed, the liquid viscosity might influence the particle behavior in liquid, and hence effect of the viscosity on mixing state should be examined carefully. Therefore, in this study, the visualization and quantification of the mixing behavior in the cross-torus paddle mixer are discussed in detail over a broad scope of operational and geometrical parameters, i.e., paddle rotational speed, filling level, liquid viscosity, and paddle obliquity. In this study, probability density is newly proposed to examine the convective and diffusive motion during the mixing process in this mixer.
In this paper, the numerical modeling and boundary treatment of Advanced DEM-CFD method are concisely described in Section 2. Further, the simulation models including the geometry of mixer as well as the input parameters are discussed in detail in Section 3. Afterwards, the visualization and quantification of mixing processes regarding the effects of paddle rotational speed, filling level, liquid viscosity, and paddle obliquity are presented in 4.1 Mixing behavior for the reference case, 4.2 Parametric study of mixing state, respectively. By performing the simulation of mixing behavior in such a mono-dispersed system, the present work provides significant evidence and novel insights to further understand the mixing mechanisms.
Section snippets
Solid phase
The current DEM-CFD method is employed to investigate the powder mixing behavior in a container filled with purified liquid. For the solid phase, the translational and rotational motions of a single spherical particle can be written as the following governing equations:where mp, v, t, Fc, Fdrag, p, Vp, I, ω, and T represent the particle mass, translational velocity of particle, time, contact force, drag force, pressure, particle volume, moment of inertia of
Calculated conditions
In this paper, to investigate the mixing performance over a broad scope of operating conditions, a series of parameters including rotational speed, particle number (also known as the filling level), liquid viscosity, and obliquity of rotating paddle are correspondingly considered and performed. A schematic view of computational geometries used in current simulation is shown in Fig. 1, which is mainly composed of three parts, a container, a torus paddle, and a cross support respectively. The
Mixing behavior for the reference case
The mixture can be initially classified into different binary groups along with the vertical axis to visualize the mixing process. Fig. 2 shows the visualization of one typical vertical classification of Case 1 with first rotation in 1.5 s. It can be seen that as the paddle rotates, the distinct homogeneous distribution of binary colored particles (i.e., the transient frame at 0 s) is gradually disturbed, a considerable amount of particles was circulated and lifted to accumulate in the moving
Concluding remarks
Studies on the powder mixing in the industrial mixers are crucial for the improved criteria of mixing equipment design and optimization under the conditions of various mixing unit operations. In this work, with the aim to provide some evidences and novel insights to further understand mixing mechanisms over a broader scope of operating conditions in a cross-torus paddle mixer, a series of numerical simulations using the Advanced DEM-CFD method were newly conducted in a mono-dispersed system.
Nomenclature
- mp
particle mass, kg
- I
inertial moment of particle, kg.m2
- v
particle velocity, m/s
- Fc
contact force, N
- Fdrag
drag force, N
- g
gravitational acceleration, N/m2
- Vp
particle volume, m3
- p
pressure, N/m2
- T
torque, N.m
- k
stiffness, N/m
- e
restitution coefficient, −
- uf
fluid velocity, m/s
- ds
solid particle diameter, m
- Cd
drag coefficient, −
- Res
particle Reynolds number, −
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgment
One of the authors, Shuo LI, gratefully acknowledges the China Scholarship Council (CSC) for the award of a Chinese Government Scholarship for Postgraduates (No. 201906380026).
References (59)
- et al.
Developments in the tools for the investigation of mixing in particulate systems - a review
Adv. Powder Technol.
(2014) - et al.
An experimentally validated DEM study of powder mixing in a paddle blade mixer
Powder Technol.
(2017) - et al.
Numerical study of the mixing efficiency of a ribbon mixer using the discrete element method
Powder Technol.
(2016) Mixing of powders and granular materials by mechanical means - a perspective
Particuology.
(2012)- et al.
Discrete particle simulation of particulate systems: a review of major applications and findings
Chem. Eng. Sci.
(2008) - et al.
Numerical study of the mixing efficiency of a batch mixer using the discrete element method
Powder Technol.
(2016) - et al.
Numerical investigation of a coarse-grain discrete element method in solid mixing in a spouted bed
Chem. Eng. J.
(2018) - et al.
Large-scale discrete element modeling in pneumatic conveying
Chem. Eng. Sci.
(2009) - et al.
CFD-DEM modelling and simulation of pneumatic conveying: a review
Powder Technol.
(2020) - et al.
DEM study of granular transport in partially filled horizontal screw conveyors
Powder Technol.
(2017)
Arbitrary-shaped wall boundary modeling based on signed distance functions for granular flow simulations
Chem. Eng. J.
Numerical investigation on the influence of air flow in a die filling process
J. Taiwan Inst. Chem. Eng.
3D DEM/CFD analysis of size-induced segregation during die filling
Powder Technol.
Numerical simulation of industrial die filling using the discrete element method
Chem. Eng. Sci.
Discrete particle simulation of two-dimensional fluidized bed
Powder Technol.
Study on a large-scale discrete element model for fine particles in a fluidized bed
Adv. Powder Technol.
Verification and validation of a coarse grain model of the DEM in a bubbling fluidized bed
Chem. Eng. J.
Simulating wet gas-solid fluidized beds using coarse-grid CFD-DEM
Chem. Eng. Sci.
Particle scale study of heat transfer in packed and fluidized beds of ellipsoidal particles
Chem. Eng. Sci.
CFD-DEM study of the effects of food grain properties on drying and shrinkage in a fluidised bed
Powder Technol.
Discrete element simulation for the evaluation of solid mixing in an industrial blender
Chem. Eng. J.
Using discrete element method to analyze the mixing of the solid particles in a slant cone mixer
Chem. Eng. Res. Des.
Analysis of the mixing of solid particles in a plowshare mixer via discrete element method (DEM)
Powder Technol.
Large-scale numerical investigation of solids mixing in a V-blender using the discrete element method
Powder Technol.
A DEM study of the mixing of particles induced by a flat blade
Chem. Eng. Sci.
Analysis of solid particle mixing in inclined fluidized beds using DEM simulation
Chem. Eng. J.
Drag and heat transfer closures for realistic numerically generated random open-cell solid foams using an immersed boundary method
Chem. Eng. Sci.
A coupled volume of fluid and immersed boundary method for simulating 3D multiphase flows with contact line dynamics in complex geometries
Chem. Eng. Sci.
Numerical simulation of two-phase flows in complex geometries by using the volume-of-fluid/immersed-boundary method
Chem. Eng. Sci.
Cited by (25)
Numerical investigation on a bimodal mixing system of solid–liquid mixture in an industrial mixing cooker
2024, Advanced Powder TechnologyNumerical study on the elucidation of powder mixing mechanism in a container blender
2023, Advanced Powder TechnologyAdvances in numerical simulation of unit operations for tablet preparation
2023, International Journal of PharmaceuticsOn POD-based modal analysis in simulations of granular flows
2023, Powder TechnologyCitation Excerpt :Owing to advances in computational hardware and algorithms, numerical modeling has become a versatile and reliable tool to elucidate the phenomena behind granular flows. The discrete element method (DEM), proposed by Cundall and Strack [9], has been extensively employed in numerous industrial processes, such as mixing [10–17], milling [18,19], compaction [20], and conveying [21–23]. Typically, in a mixing system, DEM simulation can be used to evaluate mixing performance and associated mechanisms over a broad scope of operational/geometric parameters [12–15].
Development of a reduced-order model for large-scale Eulerian–Lagrangian simulations
2022, Advanced Powder TechnologyCitation Excerpt :Granular flows are commonly encountered in various industrial processes [1], such as fluidization [2–5], conveying [6–8], mixing [9–13], milling [14,15], and die-filling [16–18].