1. Introduction
Breaking up solid particles is a widely used industrial process, especially in the paint and coating industries [
1,
2]. Other industries where this process is used are the food, medicine, and chemical industries. Typical applications of the process of breaking up particles are presented in [
3]. The final particle-size distribution has a significant influence on the product properties; therefore, planning and supervising process efficiency and performance are highly important. The product performance depends not only on its ingredients but also on its preparation method [
4]. Titanium dioxide is commonly applied as a paint and coating base in the aqueous suspension of solid particles. Further, owing to its photocatalytic properties and high refractive index, titanium dioxide has applications in the optical industry [
5]. With regard to the technical quality of the substance, titanium dioxide has a white-powder form and consists of clusters with a complex structure. The creation of agglomerates involves several stages, including the nucleation of primary particles, the sticking into primary aggregates, and the formation of clusters called agglomerates [
6]. Primary particles in aggregates are connected to each other via strong crystal bounds, which cannot be broken by the flow [
7]. Connections in agglomerates are established by adhesive forces, which are weaker than chemical bonds; therefore, they can be easily destroyed by shear forces [
8]. Before titanium-dioxide powder can be used as a paint, its structure has to be broken into smaller parts. This process is performed in special devices such as ball mills [
9,
10] and high shear mixers [
11,
12]. Ultrasonication [
5,
13] and high-pressure nozzles [
13,
14,
15,
16] are often applied. The performance features of those devices can be found in [
17]. According to this work, an agitated vessel equipped with a Rushton turbine can generate energy dissipation in a range from 0.1 to 100 W/kg, valve homogenizer at about 10
8 W/kg, and ultrasonication at about 10
9 W/kg. The phenomenon of breaking up solid-particle agglomerates in suspensions occurs when shear stresses generated in the liquid are stronger than the cluster tensile strength [
18]. The size of a particle, created during breakage, depends on the breakage mechanism. Three main possibilities of fragmentation are considered: shattering, rupture, and erosion [
8]. In deagglomeration, all these mechanisms occur simultaneously. An investigation of the changes in particle-size distribution over time enables us to identify the dominant process.
In this study, experiments and modelling of the deagglomeration of titanium dioxide powder conducted in a stirring tank equipped with high-shear impellers were investigated. Based on the division proposed in [
19], particles larger than 1 µm are named suspended particles and those smaller are called colloidal particles. In this study, suspended particles were investigated. The industrial application typically consists of three stages. The first is powder wetting, which removes air from the pores and allows the creation of a suspension. The second is deagglomeration, which imparts the desired properties to the suspension. The third is particle stabilization, which prevents the reagglomeration of particles created during the breakage [
20]. The deagglomeration process was performed in a tank mixer equipped with three types of high-shear impellers for two rotational speeds. CFD simulations were performed for three impeller geometries and three rotational speeds. Based on them, it was possible to divide the fluid zone into some specific regions and check their turbulence properties. Those values were applied in population-balance calculations. Furthermore, simulations allowed us to investigate the liquid path in the tank and its velocity, localize death zones, analyze the surface curvature, and calculate the mixing power and liquid mass flow rates near the impeller. Those values allowed us to characterize the process and estimate the efficiency of each impeller type. As mentioned earlier, turbulence properties were applied in population-balance calculations. For solving the population-balance equations, the direct quadrature method of moment (DQMOM) was used. This method was developed by [
21]. It is based on the quadrature method of moment (QMOM) described in [
22]. The application of this method for solving aggregation and breakage problems is detailed in [
23].
2. Experiments
Experiments were carried out in an unbaffled tank mixer. Its diameter was
mm. The device is presented in
Figure 1. This tank was equipped with a cooling jacket, which allows the prevention of overheating the suspension, and clean water was used as a coolant. In all experiments, an impeller was placed on the same level above the tank bottom
. The suspension height was the same in every experiment, equal to
. Three types of the impeller with diameters equal to
were used. All impellers are presented in
Figure 2. The geometry of impellers was new or hybrid, based on publications of other authors [
24,
25]. In case 1, the impeller had the shape of a saw blade. It was a flat disc with sawtooth cuts on its edge. The number of teeth was the same for every case, equal to 18. The impeller in case 2 was a modification of the previous case. In every tooth, triangular cuts were made. Five cylindrical holes of diameter
were also made. These holes were made on a circumference of diameter
in equal distance. In the last case, described as case 3, the impeller was made similar to that in case 2. Triangular cuts were also made. Instead of cylindrical holes on the circumference of diameter
, gaps were made. Gaps were not drilled perpendicular to the impeller’s face but inclined at
°. The thickness of discs was the same for each impeller, equal to
mm.
The direction of rotation was clockwise, relative to that presented in
Figure 2. For every impeller, experiments were carried out for two values of rotational speed.
The titanium dioxide used in experiments had a rutile structure. The concentration of the solid phase was 40% of the mass mixture. As a continuous phase, ultra-pure demineralized water was used. Its conductance was no higher than 0.05 µS. Samples were analyzed on the particle size analyzer LS 13,320 Beckman Coulter. This device integrates two measurement techniques: laser diffraction (LD) and polarization intensity differential scattering (PIDS). The total size range measured by the apparatus is from 0.04 to 2000 µm.
The deagglomeration process lasted 60 min. Samples were taken during the process without interrupting the working of the impeller. After that, ultra-pure water was added to samples in order to reduce the agglomeration of particles. For every measurement set, at least eight samples were taken. Due to the fast change in size at the beginning of the process, more samples were taken in the first ten minutes of the working of the impeller.
4. Mechanistic Model
To predict the particle-size changes over time in the tank, population-balance methods were employed. Methods of describing the agglomeration and deagglomeration processes are presented in [
23]. To solve the population-balance equations in this case, the direct quadrature method of moments was applied. The details of this method are presented in [
21]. In this paper, only the deagglomeration process was considered and agglomeration was neglected. Due to low values of energy dissipation rate for 500 rpm, only higher values of rotational speed were involved in population balance computations.
Particle clusters break as a result of hydrodynamic stresses generated in the liquid by the power input of an impeller. The values of those stresses can be evaluated via dimensionless analysis. This expression includes the turbulence energy dissipation rate, which can be calculated from the CFD simulation. To simplify this case, a mechanistic model was applied. In the tank domain, three characteristic zones were selected and numbered I, II, and III. The division criterion was the energy dissipation rate, [m2·s−3]. In zone I, the energy dissipation rate was higher than the average over the whole tank (i.e., ). In zone II, was lower than the average but higher than the critical value, above which the breakage process begins (i.e., ). In zone III, the energy dissipation rate was lower than the critical rate (i.e., ); therefore, deagglomeration did not occur in this region. The Fluent software allowed us to calculate the mass flow rate of the liquid near the impeller. The flow rate and zone volumes allowed us to calculate the residence time of the fluid element there. The energy dissipation rate was the average in each zone.
The following assumptions were made in the mechanistic model:
Every fluid element is tracked in the Lagrangian way;
Every fluid element has the same residence time in the selected zone;
Fluid elements are completely segregated and do not mix with each other;
Energy dissipation rates in the zones are constant;
Particles do not aggregate or agglomerate;
In
Figure 5, the liquid circulation through the dissipation zones is presented. The red, orange, and green arrows represent the flows through zones I, II, and III, respectively. The residence time and average energy dissipation rate are presented in
Table 1.
The application of DQMOM methods allows us to solve the population-balance equations with simple linear algebra; source-term specifications for equations are required. Owing to the assumption of negligible agglomeration effects, only breakage terms need to be defined. The population-balance equation for agglomerate breakage is given by the following formula [
18]:
where
is the number density function,
is the particle velocity in the
direction,
is the position vector,
is the time coordinate,
is the particle linear growth rate, and
is the characteristic length of a particle.
and
are birth and death rate functions, respectively, defined by
where
is a fragment distribution function that describes the number and size of particles
. created from the breakage of agglomerates of size
. As mentioned earlier, the deagglomeration mechanism that dominates the process is shattering. For this mechanism, the fragment-distribution function is defined by the following equation:
where
is the Dirac’s delta with a parameter
size, which represents the aggregate size.
is the number of aggregates in the agglomerate of size
. This value can be evaluated by [
14]:
where
kv and
kva are the shape factors of agglomerates and aggregates, respectively.
Df is the fractal dimension of the agglomerate.
Function Γ is a breakage kernel; it describes the frequency of deagglomeration events. It depends on the fluid properties, flow pattern, and particle properties, as reported in [
29]. Breakage occurs when hydrodynamic stresses are greater than the agglomerate tensile strength. This function is defined in two ways. The division criteria are based on the Kolmogorov length microscale [
30], given by
where
[m
2·s
−1] is the kinematic viscosity of the fluid and
[m
2·s
−3] is the turbulence energy dissipation rate. For particles larger than the Kolmogorov microscale,
, deagglomeration occurs in the inertial subrange of turbulence spectrum [
13]. The kernel is given by
where
is the fluid density.
describes the hydrodynamic stresses in the inertial subrange of the turbulence spectrum.
represents the frequency of turbulence events. For particles smaller than or equal to the Kolmogorov microscale,
, breakage occurs in the dissipation subrange of the turbulence spectrum. In this case, the breakage kernel, based on the inversion of the Kolmogorov time scale, is given by
In Equations (8) and (9),
is a constant. This phenomenon is illustrated in
Figure 6.
In
Figure 6,
is the kinetic energy of the turbulence,
is the wave number,
is the characteristic length of the largest eddies in the inertial subrange,
is the Kolmogorov length scale of the eddies, and
and
are the dynamic and kinematic viscosities of the continuous phase, respectively. The agglomerate tensile strength is a result of interactions between its components; it is a result of the occurrence of two forces: attraction and repulsion. The first one is caused by van der Waals forces. The second one is caused by the electrostatic double layer near the particle surface [
31]. The tensile strength in fractal clusters of particles was investigated and reported in [
32]. For agglomerates, this value is given by the following formula:
where
is the cluster porosity,
is the characteristic length of aggregates, and
is the bounding force between agglomerate components. Its value is the sum of attractive and repulsive forces. In this case, only attractive effects are included. The attractive forces are given by
where
is the Hamaker constant,
is the characteristic size of primary particles, and
is the surface-to-surface distance.
The agglomerate porosity,
, is given by the following formula presented in [
13]:
To discover the morphology of the substance, scanning electron microscope images were obtained on two devices. Images of particles before deagglomeration were obtained on the Hitachi SU8230. The accelerating voltage of this device was 10 kV. These samples were not coated before the examination. Images of particles after deagglomeration were obtained on the Phenom G1 from PhoenomWorld company. The accelerating voltage of the device was 5 kV. Before the examination, samples were prepared to prevent charging and were coated with a mixture of gold and palladium (80% of gold and 20% of palladium atomic rate). The coating-layer thickness was approximately 8 nm. This coating was developed by the Q150 T Sputter Turbomolecular pumped coater from the Quorum company. The samples in
Figure 7a,b present the raw material before stirring. The samples presented in
Figure 7c,d were deagglomerated as a water suspension by a high-shear mixer for 60 min.
The analysis results for the images, presented in
Figure 7, suggest that the aggregates and agglomerates had a fractal structure. Agglomerates of titanium dioxide are presented in images (a) and (b). The aggregates presented in images (c) and (d) consisted of smaller particles of submicron sizes. These clusters also had a porous structure. Image (c) shows suspension compounds at 5000× magnification. It shows the aggregates, which were torn from the larger structure. The second image shows the titanium-dioxide aggregate at 20,000× magnification. Based on the SEM images for further calculations, the characteristic size of the aggregates was equal to
= 4 µm.
To define the agglomerate behavior during the breakage, the structures of those clusters were described by a model similar to that presented in [
18]. The parameters required by the model were taken from the mentioned work and were equal;
= π/6,
= π/6,
= 2.85,
= 998.2 kg·m
−3,
= 1.003 × 10
−3 Pa·s,
= 5 × 10
−21 J,
= 1 nm, and
= 20 nm. In the breakage, the kernel constant is equal (
= 1 × 10
−6). This value is lower than 1. This means that not every turbulent event can break up a solid particle. In reference [
18], the authors suggested that this is a result of the intermittency phenomenon. The critical value of energy dissipation rate was found to be
= 0.178 m
2·s
−3 in this case.
5. Results
As a result of the experimental study, particle size distribution changes in time have been obtained. The results for impeller 1 for a rotation speed of 1200 rpm are presented in
Figure 8. The analysis shows that this product consisted of two populations of particles. During the deagglomeration process, the population of larger particles shrank and the population of smaller particles grew. Based on the division proposed in [
33], in this system, the shattering deagglomeration mechanism dominated.
The results of the mechanistic-model application are presented in
Figure 9. The residence time of the liquid element in each zone was found from the zone volume and flow rate near the impeller. Those values were calculated from the CFD results. The mean size of the particles in the system
was calculated based on 4th- and 3rd-order statistical moments
and
, respectively.
The experimental study showed that impeller 3 had the highest ability of deagglomeration. For rotational speeds of 1200 and 3000 rpm, this impeller allowed the smallest final size of the particles to be obtained. This impeller characterized the fastest particle size drop. It shows that the application of this impeller allowed the required particle size to be obtained faster than for other geometries. It shortened the time of the deagglomeration process. The application of this allowed a smaller final size of the particles to be obtained than for other studied geometries. This shows that impeller 3 for the rotational speed 3000 rpm allowed the obtaining of a mean size of 2 µm smaller particles than those of unmodified geometry 1. In addition, for the rotational speed of 1200 rpm, the application of impeller 3 allowed a particle size of 7.589 µm to be obtained. Using impeller 1 for a rotational speed of 3000 rpm allowed a particle size of 7.871 µm to be obtained after 60 min of deagglomeration.
An analysis of
Figure 9 showed good agreement of the experimental data with the population balance methods combined with a mechanistic model of suspension flow. This shows that this method can be applied to predict the results of deagglomeration in the design process. In the first stage of deagglomeration, the particle sizes obtained by experiments were smaller than those predicted by the model. This suggests that in the deagglomeration process for larger agglomerates, a different mechanism of deagglomeration dominated than for smaller one. A potential method for model improvement is to apply a rapid mechanism of deagglomeration for larger particles and a second, slower mechanism for small agglomerates. In
Figure 9f, model prediction ended after 40 min. This was caused by numerical issues. In long or intensive deagglomeration, Dirac deltas used in the quadrature approximation approached the size of the aggregate. It caused the numeric matrix solved in the DQMOM method to be near to the singular matrix, and further calculation was not necessary, because the model reached the final obtainable size possible.
6. Discussion
This experimental study showed that the application of impeller 3 allowed smaller particles than those of other studied geometries to be obtained. In addition, this impeller characterized the fastest particle size drop in time. This suggests that the application of this type of impeller allowed the desirable size of the product to be obtained faster than for other geometries and, consequently, with less energy consumption.
The application of CFD methods and the volume of the fluid multiphase model allowed us to successfully reconstruct the conditions in a tank during particle deagglomeration. It also allowed us to investigate the mixing power that enabled the inspection of the power consumption of the process. CFD simulations indicated the parameters that were useful in the population-balance model; they allowed us to find the properties of zones separated in the mechanistic model. An investigation of the liquid path in the tank, based on the simulation results, indicated the dead zones, which can cause inhomogeneities in the product. The application of population-balance equations in the mechanistic model of suspension flow in the system provided information about the tendencies of breakage ability. Combining computational fluid dynamics and population balance gave a useful tool in the designing process. It allowed the prediction of particle sizes after deagglomeration with good approximation. This approach required less time and hardware resources than the implementation of the population balance equation into the CFD code and its solution in every cell in the tank domain. The application of a mechanistic model of suspension flow simplified the population balance equation, reducing it to a one-dimensional problem—solving it only for the time variable. The computational time needed to obtain the change in particle size with time for 3600 circulations of fluid elements for impeller 3 for 1200 rpm was 174.8 ± 1.1 s, averaged over five measurements. The calculation was made on an AMD Ryzen 7 3700X CPU, on a single core with a frequency up to 4.4 GHz using the Matlab® platform. Differential equations were solved using the implemented algorithm ode45, with a fixed time step. Using an algorithm with an adaptive time step, for example, ode15s, could speed up calculations. The CFD computation required to obtain the results took 48 h, including post-processing, the use of a computing cluster, and parallel computation on a 10-core AMD Opteron 6284 SE CPU.
The proposed method is fast and accurate and can be successfully applied to the selection of mixer geometries in industrial applications. The limitation of this method is connected with the breakage kernel constant Cb that has to be found based on the experiment. The value of this constant will be different for each type of substance. It implicates the necessity of experimental study on each type of substance or literature search. Furthermore, we observed that in the first stage of deagglomeration, the particle sizes obtained by experiments were smaller than those predicted by the model. This suggests that in the deagglomeration process for larger agglomerates, a different mechanism of deagglomeration dominates than for the smaller one. The identification of the breaking mechanisms requires more experimental studies and can be considered as a further development of our research.