Dry coating in a high shear mixer: Comparison of experimental results with dem analysis of particle motions

Similar documents
Experiment and simulation of the dry particle coating

Chemical Engineering Science

DEM modeling of penetration test in static and dynamic conditions

Methylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma.

A new simple recursive algorithm for finding prime numbers using Rosser s theorem

Smart Bolometer: Toward Monolithic Bolometer with Smart Functions

Modeling the mean interaction forces between powder particles. Application to silica gelmagnesium. mixtures.

Full-order observers for linear systems with unknown inputs

Case report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122,

The FLRW cosmological model revisited: relation of the local time with th e local curvature and consequences on the Heisenberg uncertainty principle

Towards an active anechoic room

A numerical study of the penetration test at constant rod velocity

Easter bracelets for years

b-chromatic number of cacti

Solubility prediction of weak electrolyte mixtures

Quasistatic behavior and force transmission in packing of irregular polyhedral particles

Detection and Survey of Interface Defects Within a Pavement Structure with Ultrasonic Pulse Echo

Thomas Lugand. To cite this version: HAL Id: tel

Stator/Rotor Interface Analysis for Piezoelectric Motors

Can we reduce health inequalities? An analysis of the English strategy ( )

Evolution of the cooperation and consequences of a decrease in plant diversity on the root symbiont diversity

Water Vapour Effects in Mass Measurement

Dispersion relation results for VCS at JLab

Direct strain and slope measurement using 3D DSPSI

Passerelle entre les arts : la sculpture sonore

Analysis of Boyer and Moore s MJRTY algorithm

A new approach of the concept of prime number

Vibro-acoustic simulation of a car window

Natural convection of magnetic fluid inside a cubical enclosure under magnetic gravity compensation

Trench IGBT failure mechanisms evolution with temperature and gate resistance under various short-circuit conditions

Solving the neutron slowing down equation

Comment on: Sadi Carnot on Carnot s theorem.

Hook lengths and shifted parts of partitions

RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING

Rolling and sliding dynamics in centrifugal spreading

On size, radius and minimum degree

From Unstructured 3D Point Clouds to Structured Knowledge - A Semantics Approach

STATISTICAL ENERGY ANALYSIS: CORRELATION BETWEEN DIFFUSE FIELD AND ENERGY EQUIPARTITION

Question order experimental constraints on quantum-like models of judgement

Completeness of the Tree System for Propositional Classical Logic

On the longest path in a recursively partitionable graph

On Newton-Raphson iteration for multiplicative inverses modulo prime powers

Entropies and fractal dimensions

On measurement of mechanical properties of sound absorbing materials

L institution sportive : rêve et illusion

A Simple Model for Cavitation with Non-condensable Gases

Numerical modeling of diffusion within composite media

Simultaneous Induction Heating and Electromagnetic Stirring of a Molten Glass Bath

On path partitions of the divisor graph

Soundness of the System of Semantic Trees for Classical Logic based on Fitting and Smullyan

A micromechanical approach to describe internal erosion effects in soils

Ultra low frequency pressure transducer calibration

A numerical analysis of chaos in the double pendulum

On Poincare-Wirtinger inequalities in spaces of functions of bounded variation

New Basis Points of Geodetic Stations for Landslide Monitoring

Numerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral. Equation Approach

Influence of a Rough Thin Layer on the Potential

Mirage detection for electrochromic materials characterization. Application to iridium oxide films

Impulse response measurement of ultrasonic transducers

MODal ENergy Analysis

Theoretical calculation of the power of wind turbine or tidal turbine

On sl3 KZ equations and W3 null-vector equations

Prompt Photon Production in p-a Collisions at LHC and the Extraction of Gluon Shadowing

AC Transport Losses Calculation in a Bi-2223 Current Lead Using Thermal Coupling With an Analytical Formula

A proximal approach to the inversion of ill-conditioned matrices

Cutwidth and degeneracy of graphs

Optically Selected GRB Afterglows, a Real Time Analysis System at the CFHT

Quantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry

Nonlocal computational methods applied to composites structures

Lorentz force velocimetry using small-size permanent magnet systems and a multi-degree-of-freedom force/torque sensor

The magnetic field diffusion equation including dynamic, hysteresis: A linear formulation of the problem

Exact Comparison of Quadratic Irrationals

Electromagnetic characterization of magnetic steel alloys with respect to the temperature

On Symmetric Norm Inequalities And Hermitian Block-Matrices

Evaluation of transverse elastic properties of fibers used in composite materials by laser resonant ultrasound spectroscopy

A Study of the Regular Pentagon with a Classic Geometric Approach

A note on the acyclic 3-choosability of some planar graphs

Interactions of an eddy current sensor and a multilayered structure

Widely Linear Estimation with Complex Data

Reduced Models (and control) of in-situ decontamination of large water resources

Some approaches to modeling of the effective properties for thermoelastic composites

PRESSURE FLUCTUATION NUMERICAL SIMULATION IN A CENTRIFUGAL PUMP VOLUTE CASING

Predicting the risk of non-compliance to EMC requirements during the life-cycle

Characterization of Equilibrium Paths in a Two-Sector Economy with CES Production Functions and Sector-Specific Externality

Diurnal variation of tropospheric temperature at a tropical station

A NON - CONVENTIONAL TYPE OF PERMANENT MAGNET BEARING

Exogenous input estimation in Electronic Power Steering (EPS) systems

Characterization of the local Electrical Properties of Electrical Machine Parts with non-trivial Geometry

The influence of the global atmospheric properties on the detection of UHECR by EUSO on board of the ISS

A Simple Proof of P versus NP

Computation and Experimental Measurements of the Magnetic Fields between Filamentary Circular Coils

Eddy-Current Effects in Circuit Breakers During Arc Displacement Phase

Performance analysis of clouds with phase-type arrivals

Climbing discrepancy search for flowshop and jobshop scheduling with time-lags

Multiple sensor fault detection in heat exchanger system

Basic concepts and models in continuum damage mechanics

Using multitable techniques for assessing Phytoplankton Structure and Succession in the Reservoir Marne (Seine Catchment Area, France)

Quasi-periodic solutions of the 2D Euler equation

The sound power output of a monopole source in a cylindrical pipe containing area discontinuities

Sound intensity as a function of sound insulation partition

Transcription:

Dry coating in a high shear mixer: Comparison of experimental results with dem analysis of particle motions Eric Serris, Akira Sato, Alain Chamayou, Laurence Galet, Michel Baron, Philippe Grosseau, Gérard Thomas To cite this version: Eric Serris, Akira Sato, Alain Chamayou, Laurence Galet, Michel Baron, et al.. Dry coating in a high shear mixer: Comparison of experimental results with dem analysis of particle motions. Aibing Yu, Kejun Dong, Runyu Yang, Stefan Luding. Powders and Grains 2013, Jul 2013, Sydney, Australia. American Institut of Physics Publishing LLC, 1542, pp.779, 2013, AIP Conference Proceedings. <10.1063/1.4812047>. <hal-00846905> HAL Id: hal-00846905 https://hal.archives-ouvertes.fr/hal-00846905 Submitted on 25 Jul 2013 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Dry Coating In A High Shear Mixer: Comparison Of Experimental Results With DEM Analysis Of Particle Motions E.Serris (2), A. Sato (1), (2), A. Chamayou (1), L. Galet (1), M. Baron (1), P. Grosseau (2), G. Thomas (2) (1) Université de Toulouse, Mines-Albi, CNRS-UMR 5302, Rapsodee, Campus Jarlard, 81013 Albi Cedex 09, France (2) Ecole Nationale Supérieure des mines, SPIN-EMSE, CNRS UMR 5307 LGF, F-42023 SAINT-ETIENNE, France Abstract. Experimental dry coating of guest particles on the surface of host particles is performed by mechanical forces in a high shear mixer called Cyclomix. The studied system (a mixture of particles of sugar, Suglet as host particles and magnesium stearate as guest particles) was chosen as a model one to achieve better understandings of the phenomena during mixing. To simulate the flow of host/guest particles in the mixer, the Discrete Element Method (DEM) was applied. Experimental results such as flowability and wettability can be explained by particles flows evolutions with different rotational speed or duration treatment inside the Cyclomix. Keywords: Dry coating, Mixer, DEM, particle motion PACS: 81.20.Ev: INTRODUCTION Dry particle coating is a process allowing designated modifications of host particles properties. A coating of fine guest particles onto the surface of host particles is achieved using mechanical forces such as shear and impact forces, generally occurring during mixing. This process has received a widespread attention by many industrial fields because it consists in an environmental friendly and low cost operation avoiding organic solvents. The optimization of a dry coating process requires normally the experimental efforts and empirical techniques [1-3]. Several theoretical approaches have been reported [4,5]. Meanwhile, the Discrete Element Method (DEM), proposed by Cundall and Strack [6], let us trace precisely the motion of each particle. The DEM is one of the most popular numerical methods for simulating and analyzing the solid particle behavior and has been successfully applied in many fields such as soil and rock mechanics [7], agricultural product handling [8], storage and flow of granular materials [9] and so on. DEM has already been applied to several studies involving other dry coating apparatus [10,11]. CFD modeling of particles motion inside the Cyclomix has been applied for granulation experiments [12, 13]. This work contributes to apply DEM to the simulation of the particle motion in a high shear mixer, and its analysis regarding the evolutions of physical properties obtained by experimental work. EXPERIMENTAL Sample Powders The sample powders chosen are sugar particles (Suglets ) and magnesium stearate (MgSt) for host and invited particles respectively. Suglets, products of Colorcon INC., are spherical cores mainly composed of sucrose and maize starch with a highly hydrophilic character. MgSt supplied by Chimiray chosen as invited particles is a fine, white, cohesive and hydrophobic powder. The size distribution of Suglets seems to have a one shape population with median diameter (D50) of about 250 µm, and indicates that this powder has a very homogeneous size, in comparison with that of MgSt having a wide population ranging from 0.1 µm to 50 µm. D50 of MgSt is about 5 µm.

Coating Device Velocity and motion analysis The coating device used in the experiment is a high shear mixer (Cyclomix, Hosokawa micron B.V. Japan), which is generally used for granulation. This device has been successfully used for dry coating [14]. The schematic diagram of the experimental equipment is shown in Fig. 1. This device presents basically a stationary conical chamber and a vertical rotor that can rotate clockwise direction from 200 rpm to 3000 rpm, with four pairs of flat-bladed impellers from bottom to top. Free volume of the chamber is 1000 cm3. The particle motion, as well as the velocity of each particle has been analyzed from the snapshots of the simulations with different operating conditions. The color shows the normalized dimensionless velocity (red represents 1 and blue represents 0). FIGURE 1 : Schematic diagram of the high shear mixer. RESULT AND DISCUSSION Numerical analysis The particle motion in the Cyclomix has been simulated by the discrete element method (DEM). The particles with the properties of the host particles given in Table 1 were placed in the vessel of the Cyclomix with an hexagonal-close packed structure. TABLE 1. Simulation parameters. Young s E [GPa] modulus Poisson s ratio [-] Density [kg/m3] Coefficient of [-] restitutuion Time step t [µs] T [s] Total simulation time 4.5 0.30 1600 0.16 FIGURE 2. Normalized velocity of particle in color at each rotational speed: (a) 100 rpm, (b) 500 rpm, (c) 1000 rpm and (d) 3000 rpm, particle size 2 mm, filling ratio J = 60 % in cross section. Fig.2 shows the normalized velocity of a particle in color at each rotational speed: (a) 100 rpm, (b) 500 rpm, (c) 1000 rpm and (d) 3000 rpm in cross section views. In the cross section, the particle motion is revealed clearer, at low rotational speed (100 rpm), the particle bed stays below the top impeller, and particles around paddle have relatively high normalized velocity. At 500 rpm, the powder bed reaches the top of the chamber, and at 1000 rpm, most of the particles go to the upper side. Finally at 3000 rpm, the particles stay either on the upper side or bottom, and there are only few particles around the middle part of the chamber. Simulation of host and guest particles without adhesion 2.0 1.5 Gravity effects have been introduced as an external condition imposed to the assembly yet submitted to the forces developed by the rotating impellers. The current computer power is not enough to simulate the actual number of particles. Then a smaller number of particles (about 40 000), but with larger size (2.0 mm) has been chosen. The simulations have been performed in order to know the effect of rotational speed on the mixing of host and guest particles. The rotational speed varies from 250 to 1500 rpm. The filling ratio of samples is fixed at 40 % and 20 % for host and guest particles, respectively. Simulation time is also fixed at 3.0 s. Particle size is also fixed at 3 mm and 1 mm for host and guest particles respectively. The physical parameter used for the host and guest particles are summarized in Table 2. Initially the invited particles are located over the powder bed of host particles.

TABLE 2. Physical characteristics of host and guest particles material host guest Young s modulus E [GPa] 4.5 2.2 Poisson s ratio [-] 0.30 0.30 Density [kg/m3] 1600 1120 Coefficient of restitution [-] 0.3 0.0 Number of particles [-] 14 148 190 986 Total simulation time T [s] 3.0 To observe the inside of the mixing chamber, the cross sections of the mixing chamber at each rotational speed are shown in Fig. 3 (Guest in Blue and host in Red). At 250 rpm the guest particles are located around the surface of the powder bed and rotating axis (Fig. 3 (a)). FIGURE 4 : Dependence of the contact angle on the operating time t, and on the different rotational speeds ω for a filling ratio J = 60%. FIGURE 3 : Cross section of the mixing chamber at each rotating speed: (a) 250 rpm, (b) 500 rpm, (c) 1500 rpm, T= 3s (host in red and guest in blue). At higher rotational speeds, 500 rpm, the guest particles are located along the side wall (Fig. 3 (b)). It is suggested that the dry coating occurs in this region. At 1500 rpm the guest particles are evenly located around the side wall and especially there are a lot of guest particles seen at the upper part of the mixer (Fig. 3 (c)). Since the guest particles are smaller and lighter than the host particles, they are likely to go upper side. Discussion The sessile water drop test has been carried out to analyze the wettability of the products. Fig. 4 shows the contact angle measured 30 s after the water drop is placed on the powder surface for different processing times t, different speeds of rotation ω and at a filling ratio of J = 60 %. In Fig. 4, at 250 rpm, the contact angle increases with an increase of processing times. From the assumption above, the surface fraction of the coated material will also increase. In other words, the coating process proceeds as a function of the processing time; at 500 rpm, it increases more rapidly than at 250 rpm. At 1000 rpm, (t) increases with time and is asymptote to the contact angle I of guest particles. At 1500 rpm, it increases until 180 s and then starts to decrease gradually. One explanation of this drop of contact angle is that new surfaces of suglets not covered by MgSt are created. It could be said that once fine fragments of the host particles are generated by the attrition or erosion during the coating operation, after that, those fine fragment of host particles stick to the surface of the host particles again. The contact angle and surface fraction of coated particles could decrease at high speed coating operation by those phenomena. Regarding the evolutions of particles motions found by DEM, we will try to explain the evolution of physical properties of composite powders formed by mixing inside the cyclomix. All experimental results are explained [15,16]. We have divided the coating process inside the Cyclomix into several steps strongly depending on the rotational speed. The first step of discrete coating is especially occurring at low rotational speeds. The evolution of wettability is very slow and the maximum contact angle showed in fig 4 for 250 rpm do not reach the value of 132, so we can tell that the coating is not complete at the surface of suglets. If we look at particle motions in fig 2 and 3 we can see that the particles are not individualized and packed at the lower part of the cyclomix. So, it takes a long time for small particle to stick on the surface of big particles.

For these small rotational speeds we cannot achieve complete coating due to insufficient mixing of both powders clearly shown in Fig 3a. For rotational speed of 500 rpm and 1000 rpm, the complete coating is achieved and the properties do not vary at a steady state step. Wettability of powders reached the same value for a 10 minutes mixing time. If we look at the particle motions, the powder bed reaches the top of the Cyclomix Fig 2b) and Fig 2c) and the fig 2b) shows that that the small particle are on the wall of the Cyclomix. So the coating can occur at these rotational speeds and is be located along the walls of the Cyclomix. For higher rotational speed (faster than 1000 rpm), a particle breakage can occur; this is shown with the drop of contact angle from 132 to 125. The particle motions were quite the same for 1000 rpm and 3000 rpm and the coating occurs rapidly (in less than one minute). After that, this higher energy given to the system could generate brakeage and / or attrition of suglets [15,16]. For this host/invited particles system these rotational speed are not good for dry coating. CONCLUSIONS In this paper, the particle motion of particles in a high shear mixer has been simulated by DEM. The validation of the simulation work was discussed in comparison with experimental data (such as wettability). The velocity of the particles simulated by the DEM can explain the behavior of the coating and helped us to better understand several phenomena. The coating at low rotational speed was not enough efficient because the small particles were not really mixed to the larger particles and they still remain stuck together in large agglomerates. For intermediate rotational speed, the coating occurred and seemed to be located along the Cyclomix walls. For higher rotational speeds, the energy given to the particles was high enough to generate breakage. These results confirm the applicability of the proposed method for simulating the particles motions in the dry coating. This methodology should contribute to the evaluation of the dry particle coating efficiency from mechanical studies (i;e. calculate contacts forces of particles / particles / impellers/wall) of characterization of particles contacts in order to extend predictions to other coating systems. REFERENCES of silica particles by dry coating : Characterization and powder ageing, Powder Technology, 190a (2009 ) 200-209 3 S. Otles, O. Lecoq, J. Dodds, Dry particle high coating of biopowders: An energy approach, Powder Technology, 208 (2011) 378-382 4 T. Iwasaki, M. Satoh, T. Ito, Determination of optimum operating conditions based on energy requirements for dry particle coating in a dry process, Powder Technology 123 (2002) 108-113 5 G. Thomas, Y. Ouabbas, P. Grosseau, M. Baron, A. Chamayou, L. Galet, Modeling the mean interaction forces between powder particles : Application to silica gel-magnesium stearate mixtures, Applied Surface Science, 255 (2009)7500-7507 6 P.A. Cundall, O.D.L. Strack, A discrete numerical model for granular assemblies, Geotechnique, 29 (1979) 47-65 7 R.D. Morrisson, F. Shi, R. Whyte, Modelling of incremental rock breakage by impact For use in DEM models, Minerals Engineering, 20 (2007) 303-309 8 M. Wojtkowski, J. Pecen, J. Horabik, M. Molenda, Rapeseed impact against a flat surface: Physical testing and DEM simulation with two contact models, Powder Technology, 198 (2010) 61-68 9 B. ChaudhuriI, F.J. Muzzio, M.S. Tomassone, Experimentally validated computations of heat transfer in granular materials in rotary calciners, Powder Technology, 198 (2010) 6-15 10 R. Dave, W. Chen, A. Mujumdar, W. Wang, R. Pfeffer, Numerical simulation of dry particle coating processes by the discrete element method, Advanced Powder Technology, 14 (2003) 449-470 11 W. Chen, R. Dave, R. Pfeffer, O. Walton, Numerical simulation of mechanofusion system, Powder Technology, 149 (2004) 121-136 12 B.H. Ng,.Y.L. Ding, M. Ghadiri, Modelling of dense and complex granular flow in high shear mixer granulator-a CFD approach, Chemical engineering science, 64 (2009) 3622-3632 13 A. Hassanpour,.S.J. Antony, M. Ghadiri, Modeling of agglomerates behavior under shear deformation : effect of velocity field of a high shear mixer granulation on the structure of agglomerates, Advanced Powder technology18(6) (2007) 803-811 14 Y. Ouabbas, J. Dodds, L. Galet, A. Chamayou, M. Baron, Particle-particle coating in a cyclomix impact mixer Powder Technology, 189b (2009) 245-252 15 A.Sato, E.Serris, P.Grosseau, G.Thomas, A.Chamayou, L. Galet, M.Baron, Effect of operating conditions on dry particle coating in a high shear mixer, Powder Technology 229 (2012) 97 103 16 A.Sato, E.Serris, P.Grosseau, G.Thomas, A.Chamayou, L. Galet, M.Baron, Chemical engineering science, 86 (2013) 164-172 1. M. Alonso, M. Satoh, K. Myanami, Powder coating rotary mixer rocking motion, Powder Technology, 56 (1988) 135-141 2 Y. Ouabbas, A. Chamayou, L. Galet, M. Baron, G. Thomas, P. Grosseau, B. Guilhot, Surface modification