## Abstract

Sedimentation processes are fundamental to solids/liquid separation in water and wastewater treatment, and therefore a robust understanding of the settlement characteristics of mass fractal aggregates (flocs) formed in the flocculation stage is fundamental to optimized settlement tank design and operation. However, the use of settling as a technique to determine aggregates' traits is limited by current understanding of permeability. In this paper, we combine experimental and numerical approaches to assess settling velocities of fractal aggregates. Using a non-intrusive *in situ* digital image-based method, three- and two-dimensional fractal dimensions were calculated for kaolin-based flocs. By considering shape and fractal dimension, the porosity, density and settling velocities of the flocs were calculated individually, and settling velocities compared with those of spheres of the same density using Stokes' law. Shape analysis shows that the settling velocities for fractal aggregates may be greater or less than those for perfect spheres. For example, fractal aggregates with floc fractal dimension, *D*_{f} = 2.61, floc size, *d*_{f} > 320 μm and *d*_{p} = 7.5 μm settle with lower velocities than those predicted by Stokes' law; whilst, for *D*_{f} = 2.33, all aggregates of *d*_{f} > 70 μm and *d*_{p} = 7.5 μm settled below the velocity calculated by Stokes' law for spheres. Conversely, fractal settling velocities were higher than spheres for all the range of sizes, when *D*_{f} of 2.83 was simulated. The ratio of fractal aggregate to sphere settling velocity (the former being obtained from fractal porosity and density considerations), varied from 0.16 to 4.11 for aggregates in the range of 10 and 1,000 μm, primary particle size of 7.5 μm and a three-dimensional fractal dimension between 2.33 and 2.83. However, the ratio decreases to the range of 0.04–2.92 when primary particle size changes to 1.0 μm for the same fractal dimensions. Using the floc analysis technique developed here, the results demonstrate the difference in settlement behaviour between the approach developed here and the traditional Stokes' law approach using solid spheres. The technique and results demonstrate the improvements in understanding, and hence value to be derived, from an analysis based on fractal, rather than Euclidean, geometry when considering flocculation and subsequent clarification performance.

## INTRODUCTION

Coagulation and flocculation are critical steps in drinking water treatment. During coagulation, a coagulant (e.g. aluminium sulphate) is added to the water to neutralize the negatively charged suspended particles, following which shear-induced interactions generated by slow mixing in a flocculator cause their aggregation into mass fractal aggregates (flocs). Clarification processes in water treatment are often based on sedimentation of flocs. In order to optimize sedimentation, it is important to gain a fundamental understanding of floc settling behaviour. Traditional approaches have used Stokes' law with its inherent assumptions of solid, spherical particles. Flocs are clearly far more complex in structure and behaviour and, therefore, it is important that the impacts of these features of flocs are understood in order that their settling behaviour can be assessed and so predicted with accuracy.

Floc aggregation is a dynamic process where mass, surface area, number, and morphology change as functions of the shear stress and time. Furthermore, whilst aggregates may have the same size, they may exhibit different structures due to different particle arrangements during aggregation. The variations in mass, surface area and concentration substantially affect floc behaviour, particularly with regard to collision and to aggregation efficiency (Vahedi & Gorczyca 2012).

Therefore, the irregular shape and porous structure of fractal aggregates may oppose the predicted settling rate using solid spheres. These traits may affect aggregates' density and porosity and hence the drag forces on the floc surface, making velocities (and hydraulic loading rates) slower than the ones predicted by Stokes' law (Vahedi & Gorczyca 2012; Chakraborti & Kaur 2014). This may potentially have influences on the predicted clarification efficiency of settling units, resulting in more solids being dragged out from the tank into the subsequent filters, which may not be designed for such extra loads. Hence, the correct prediction of aggregates' settling rates is crucial for the accurate design of sedimentation tanks and subsequent filters in a full-scale water treatment plant.

Fractal geometry has been extensively used to characterize non-uniform objects and its application to sedimentation will be described below, after a brief description of its main characteristic.

## FRACTAL GEOMETRY

*V*,

*A*and

*L*refer to fractal volume, area and length respectively, whilst

*D*

_{f}is the floc fractal dimension.

Gregory (2009) described flocs as mass fractal objects and found that their fractal structure has important practical implications for floc density. Indeed, several factors affect flocs' fractal structure, including mixer device (Logan & Kilps 1995) and coagulation (Xu *et al.* 2010, 2011), as floc strength is a function of the formation process (He *et al.* 2012). Gregory (2009) found that aggregates formed by perikinetic flocculation display lower fractal dimension than those formed during orthokinetic flocculation, whilst aggregates formed during sweep coagulation are of larger size and fractal dimension compared with those formed during charge neutralization (Kim *et al.* 2001; Li *et al.* 2006).

*D*

_{f}∼ 2, for two-dimensional objects), whereas smaller fractal dimensions result from highly branched structures. It is believed that these more compact aggregates, rather than those exhibiting large overall size, exhibit enhanced performance during sedimentation given that the floc settling velocity depends on the aggregate size and fractal dimension (Johnson

*et al.*1996; Gregory 1997; Chakraborti

*et al.*2000; Vahedi & Gorczyca 2012), as shown in Equation (4):where:

*v*

_{st}is the settling velocity of an individual floc;

*g*is acceleration due to gravity; μ is dynamic viscosity of water;

*D*

_{f}is the fractal dimension;

*d*is floc size;

*d*

_{50}is the median size of particles within a floc (i.e. primary particles);

*θ*is a dimensionless particle shape factor; and

*ρ*

_{p}and

*ρ*

_{w}are the densities of primary particles and water, respectively.

*ν*is the kinematic viscosity of water.

*et al.*(1999) found that the floc effective density (

*ρ*

_{e}), given by , is related to floc size (

*d*) via Equation (8):where:

*A*is a packing factor and is a function of coagulant dose and coagulation pH; and

*n*is a coefficient that is a function of

*D*

_{f}. For aluminium floc and humic substances, the value of

*n*varies from 1.8 to 2.0.

Considering the importance of structural and morphological characteristics of fractal aggregates, the objective of the work reported in this paper is to investigate the porosity and density of aggregates after flocculation based on their fractal dimension calculated individually, measured by a non-intrusive image-analysis-based method, in order to inform our understanding of floc sedimentation and so improve clarification performance.

## METHODS

### Suspensions

Suspensions were prepared in the laboratory from a kaolinite solution following Yukselen & Gregory (2004) to obtain a turbidity of 25 NTU. The kaolinite was characterized by scanning electron microscopy (SEM) using a JEOL JSM-6010LA microscope coupled to an X-ray spectrometer and via laser granulometry using the Malvern Mastersizer 2000 to obtain composition and mean grain size. The suspension was coagulated with commercial aluminium sulphate, and the pH was adjusted with a 1 M sodium hydroxide solution. The optimum coagulation conditions identified by de Oliveira *et al.* (2015) using the same synthetic water were adopted, i.e., 2 mg Al^{+3}·L^{−1}, pH 7.5. Flocculation was performed with mean velocity gradients () of 20–60 s^{−1}, which were determined by torque measurements, for flocculation times of 15 min to obtain flocs of different sizes and shapes.

### Data acquisition

A total of 1,000 floc images were captured immediately post-flocculation (high-speed Miro EX-4 camera with interchangeable lenses, sampling at 25 Hz for 40 s). Image resolution was 800 × 600 pixels with a visual field of 6 × 8 mm. Shutter speed was set to 800 μs and pixel size was 10 μm. Illumination was via a collimated laser beam modified by a cylindrical lens for plane dispersion with a thickness of 2 mm (following Moruzzi *et al.* 2017) placed perpendicular to the focus direction. Nominal laser power was 2,000 mW, producing light at a wavelength of 532 nm (green). A schematic of the experimental apparatus is shown in Figure 1.

Digital images were processed using Image Pro Plus 7.0 software. The images were transformed into binary matrices by the segmentation process using a degree of 128/256 as the threshold and were then subjected to the particle image velocimeter (PIV) processing of the same package. Area and Diameter (max) of each floc were used to track each aggregate, with cut-off values of 200 pixels respectively (Chakraborti *et al.* 2003). In total, 118 aggregates were selected, for which the attributes of interest (Diameter (max), Diameter (min), Diameter (mean), Area, Y coordinate, and Perimeter) were obtained. Further details on data acquisition and treatment can be found in Moruzzi *et al.* (2017, 2019).

### Determination of fractal characteristics

*D*

_{f}and the shape parameter,

*b*, for each aggregate can be obtained from linear regression of the experimental data in the linearized form of Equation (9) (i.e. Equation (10)). Here, the slope of the straight line corresponds to the exponent (

*D*

_{f}) and the intercept corresponds to the shape parameter

*b*:where:

*D*

_{f}is the three-dimensional fractal dimension;

*d*

_{p}is the primary particle diameter (m);

*d*

_{f}is the floc diameter (m);

*N*is the number of primary particles of diameter

*d*

_{p}per volume of floc of diameter

*d*

_{f}; and

*b*is the shape parameter, defined as:where:

*ζ*is the packing factor;

*ξ*is the shape factor; and

*ξ*

_{0}is the primary particle shape factor.

*N*, was determined by rotation of the ellipsoid about the

*x*-axis of the ellipsoid fitted to the highest and lowest dimensions determined by the image analyses, as illustrated in Figure 2 and Equations (12) and (13), following Chakraborti

*et al.*(2000):where:

*A*,

*B*, and

*C*are the integration limits of the

*a*,

*b*, and

*c*axes of the ellipsoid (

*E*) with volume

*V*(

*E*); and

*S*is the surface of the ellipsoid (

*E*).

*D*

_{fp}) was obtained for the set of aggregates by fitting the volumes to Equation (14) in the linearized form:where:

*V*is the volume of the ellipsoid containing the floc (m

^{3}); and

*d*

_{max}is the largest dimension of the floc (m).

*D*

_{fp}′) for the set of aggregates was determined using Equation (15), and the ratio of

*D*

_{fp}:

*D*

_{fp}′ for the data set was then used to calculate the three-dimensional fractal dimension per aggregate according to Equation (16). The value of

*d*

_{p}in Equation (16) was assumed from the kaolin volume distribution, adopting mono-sized primary particles as simplification for the

*D*

_{f}calculation.where:

*A*is the projected floc area on the image plane; and

*D*

_{fp}′ is the two-dimensional fractal dimension.

With the shape, floc dimension, particle dimension, and three-dimensional fractal dimension parameters calculated per particle, the floc porosity and density were determined with Equations (17) and (18), considering a mass balance between the floc, the particle, and the voids occupied by the liquid.

### Settling velocity

*V*

_{sphere}) and fractals (

*V*

_{fractal}). It was assumed that the floc dimensions did not change during sedimentation. For this purpose, the Newton equation (Equation (19)) in equilibrium (Σ

*F*

_{y}= 0) was used for the particular case where the dimensionless Reynolds number (Re) <1, such that the drag coefficient can be described according to Equation (20). Thus,

*V*

_{sphere}can be described by Equation (21) and

*V*

_{fractal}by Equation (22), and they are distinguished from each other in the geometric and density terms of the equations. The geometric term of Equation (22) will always reduce the contribution of size to velocity for fractal aggregates, unless

*d*

_{f}approaches

*d*

_{p}(for

*d*

_{f}>

*d*

_{p}), resulting in the square relation of Equation (21) (

*d*

_{f}

^{2}). Finally, the ratio of the velocities was determined for each particle according to Equation (23).where:

*V*

_{sphere}is the Newton velocity for the Euclidean sphere (m·s

^{−1});

*g*is the acceleration due to gravity (m·s

^{−2}); and

*C*

_{d}is the drag coefficient.

## RESULTS AND DISCUSSION

Figure 3 shows the kaolinite (Kaolin, Fluka) particle size distribution: 0.4–100 μm in volume with median 7.5 μm and in number (*N*) with median 1.0 μm, which is in agreement with other studies (Zbik & Smart 1998; Aparicio *et al.* 2004). An example of one scanning electron microscopy image taken from the kaolin dry sample is also presented in detail, showing qualitatively the range of size, shapes and textures of the primary particles. This result supports the definition of the appropriate pixel size for the determination of the smaller cluster of the aggregate, and primary particle size adopted herein as well.

Figure 4 shows a post-flocculation (*G*_{f} = 20 s^{−1}, *T*_{f} = 15 min) aggregate characterization image. It is clear that any assumption of solid sphericity, and hence use of traditional Stokes' law approach, is inappropriate for representing the shape of the aggregates. Further, the existence of voids within the floc alters the porosity and density of the aggregate, both of which affect floc terminal velocity.

Figure 5 shows the relationship between log volume of the encased ellipsoid and the log of the ratio of *d*_{max}:*d*_{p}, from which the values *b**=* 0.78 and *D*_{fp} = 2.35 were determined (Equation (10)). The value of *D*_{fp} represents elongated medium aggregates, far from a perfect sphere, and agrees with other studies (Johnson *et al.* 1996; Chakraborti *et al.* 2003; Li *et al.* 2006). Plotting log aggregate area against log (*d*_{max}/*d*_{p}) yielded a two-dimensional fractal dimension, *D*_{fp}′, value of 1.50 (Figure S1, Supplementary Information).

These results were then applied to Equations (16)–(18) to determine the three-dimensional fractal dimension, porosity and density for the flocs.

Table 1 presents descriptive statistics obtained from all the aggregates in terms of the highest dimension (*D*_{max}), lowest dimension (*D*_{min}), mean dimension (*D*_{mean}), aggregate fractal dimension (*D*_{f}), porosity (*ɛ*), sphericity (*Ψ*) and density (*ρ*). Results show that the mean of the flocs' highest dimension (*D*_{max}) was 362 μm, with a maximum of 816 μm, a standard deviation of 121 μm, and a confidence interval (0.05) of 21 μm. For the flocs' lowest dimension (*D*_{min}), a mean of 138 μm was obtained, with a maximum of 310 μm, a standard deviation of 47 μm, and a confidence interval for the mean (0.05) of 8 μm. The average floc size obtained in the experiments was 231 μm, with a maximum of 451 μm, a standard deviation of 70 μm, and a confidence interval for the mean (0.05) of 12 μm. The three-dimensional fractal dimension (*D*_{f}) calculated per aggregate has a mean of 2.61 for the interval 2.33–2.83, indicating a variety of floc shapes, ranging from the most elongated to those close to spheres, in the limits of the interval. For the experiments, density was obtained for the flocs (*ρ*_{floc}) with a mean of 1,068 kg·m^{−3} within the range of 1,024–1,138 kg·m^{−3}. The mean porosity (*ɛ*) of the aggregates was 0.76 with an interval of 0.53–0.91. Sphericity (*Ψ*) and aspect ratio (*D*_{max}/*D*_{min}) of 0.58 and 2.62 average, respectively, showed that flocs are majority elongated structures, far from spherical.

Attribute . | Average . | Standard deviation . | Max . | Min . | 95%* . |
---|---|---|---|---|---|

D_{max} (μm) | 362 | 121 | 816 | 196 | 21 |

D_{min} (μm) | 138 | 47 | 310 | 73 | 8 |

D_{average} (μm) | 231 | 70 | 451 | 156 | 12 |

D_{f} (−) | 2.61 | 0.09 | 2.83 | 2.33 | 0.02 |

ρ_{f} (kg m^{−3}) | 1,068 | 22 | 1,138 | 1,024 | 4 |

ɛ (−) | 0.76 | 0.07 | 0.91 | 0.53 | 0.01 |

Ψ (−) | 0.58 | 0.14 | 0.97 | 0.14 | 0.02 |

D_{max}/D_{min} | 2.62 | 0.67 | 5.14 | 1.53 | 0.12 |

Attribute . | Average . | Standard deviation . | Max . | Min . | 95%* . |
---|---|---|---|---|---|

D_{max} (μm) | 362 | 121 | 816 | 196 | 21 |

D_{min} (μm) | 138 | 47 | 310 | 73 | 8 |

D_{average} (μm) | 231 | 70 | 451 | 156 | 12 |

D_{f} (−) | 2.61 | 0.09 | 2.83 | 2.33 | 0.02 |

ρ_{f} (kg m^{−3}) | 1,068 | 22 | 1,138 | 1,024 | 4 |

ɛ (−) | 0.76 | 0.07 | 0.91 | 0.53 | 0.01 |

Ψ (−) | 0.58 | 0.14 | 0.97 | 0.14 | 0.02 |

D_{max}/D_{min} | 2.62 | 0.67 | 5.14 | 1.53 | 0.12 |

*95% significance.

Figure 6 shows the relationship between porosity and fractal dimension, with the more open structures (*D*_{f} < 2.5) having a greater porosity (>80%) than more closed structures. To a lesser extent than shape, large size flocs have also shown greater porosity for the same fractal dimension. The average porosity found here (approximately 76%) is in agreement with Gorczyca & Ganczarczyk (1999) and Vahedi & Gorczyca (2012) for large flocs. There is no doubt that shape and porosity are dependant floc traits which may affect settling rate, however, the nonhomogeneous mass distribution and pore population within aggregate structure are also important for terminal velocity, and should be considered for further developments, as suggested by Vahedi & Gorczyca (2014).

The direct consequence of the shape and porosity relation on floc settling rates is the influences on an aggregate's density and drag. Unless water can flow through floc voids, reducing the resistance to settling, an elongated aggregate is more likely to have lower settling rates than those closer to a sphere-shape. The so-called flow-through effect is still not well understood and there is also no consensus on the contribution of permeability to the terminal velocities of porous aggregates. Whilst initially Adler (1987) showed that the effect of permeability on floc settling is minor and that it is unlikely water can flow through pores during sedimentation, nevertheless Johnson *et al.* (1996) showed the more porous the floc is, the more likely fluid can flow through it, increasing permeability and so reducing resistance to settling. In general, the literature has shown that elongated aggregates are likely to settle at lower rates than spheres and, according to Bushell *et al.* (2002), it is difficult to explain the greatly reduced resistance to settling described by Johnson *et al.* (1996).

Figure 7 shows that more spherical, less elongated aggregates (with increased *D*_{f} values) display a greater density, contributing to higher terminal velocity of aggregates. Therefore, the density of the aggregates varies inversely with porosity. If the mass was the predominant effect on floc settling velocity, compact spheres would always exhibit higher settling velocities. In fact, several researchers have shown that fractal aggregates can settle with slower velocities than those predicted by Stoke's law (Tambo & Watanabe 1979; Khelifa & Hill 2006; Jarvis *et al.* 2008; Vahedi & Gorczyca 2012).

Consequently, terminal velocities of aggregates will mostly depend upon the dual effect of porosity and density, determined to a great extent by aggregate shape, i.e fractal dimension.

Common practice for designers of sedimentation tanks is to adopt average floc density for perfect spheres, so that settling velocities of aggregates can be simulated using a characteristic aggregate size and Stokes' law, as Equation (21). However, results shown here suggest that flocs are elongated, porous structures, and far from spherical.

The impact of this oversimplification of floc shape is that the drag force changes as result of the cross-sectional area, thus altering the settling velocity. In general, it is supposed that deviation from sphericity will result in increasing drag, irrespective of floc orientation and permeability (Bushell *et al.* 2002). However, Johnson *et al.* (1996) suggest that the actual drag is lower than that calculated for fractal aggregates, even when permeability is considered. Figure 8 shows calculated settling velocity for spheres and flocs (based on Equations (21) and (22) for *d*_{p} of 7.5 μm, this being median particle size by volume) against floc size for different shapes, i.e. fractal dimensions (*D*_{f} = 3.00 for spheres; 2.33, 2.61 and 2.83 for fractals).

Figure 8 shows that fractal aggregates of size between 100 and 300 μm settle with velocities between 0.3 and 7.9 mm/s, which is in agreement with results presented by Vahedi & Gorczyca (2012), Khelifa & Hill (2006) and Jarvis *et al.* (2008), who performed both experiments and simulations. The authors showed settling velocities varying from 0.1 to 7.1 mm/s for aggregate size between 100 and 300 μm, corroborating that the simulations presented here are within experimental measurements performed by several authors. Nevertheless, results compiled by Khelifa & Hill (2006) reveal there is considerable scatter in settling velocity of fractal aggregates, varying up to 100-fold from each other for the same floc size, possibly due to the nonhomogeneous mass distribution and pore population mentioned by Vahedi & Gorczyca (2014). Here, aggregates with a fractal dimension of 2.83 settled with higher velocities than spheres of the same size, for all ranges of floc sizes, due to the usual assumption that spherical flocs have the same density as fractal aggregates, and so less mass than aggregates of equal size.

Conversely, for aggregates with fractal dimension of 2.33 and 2.61 there is a size threshold above which fractals settle with lower velocities than those calculated using Stokes' law for compact spheres. For small floc sizes, close to the size of the primary particle (*d*_{f} ≈ *d*_{p}), the geometric term of Equation (22) approaches the term *d*_{f}^{2} from Equation (21) and, therefore, the differential density between particle and liquid is the predominant driving effect over settling rate, surpassing the reduction caused by the fractal geometry. However, the differential density is not enough to surpass the reduction in the geometry term of Equation (22), caused by the low fractal dimension of elongated large aggregates. Results presented by Vahedi & Gorczyca (2012) have also shown that multi-fractal aggregates can settle with lower velocities than those predicted by the modified Stokes' law, by introducing a porous effect into the original equation. The authors simulated fractal aggregates with multi-fractal dimension varying from 2.6 to 2.7 and aggregate size less than 320 μm. In contrast, experiments carried out by Johnson *et al.* (1996) showed that settling velocities of fractal aggregates were between four and eight times greater than those predicted by Stokes' law, for aggregate size in the range of 100–1,000 μm. Here, fractal aggregates (*D*_{f} = 2.61, *d*_{f} > 320 μm) settle with lower velocities than those predicted by Stokes' law; whilst, for *D*_{f} = 2.33, all aggregates of *d*_{f} > 70 μm and *d*_{p} = 7.5 μm settled below the velocity calculated by Stokes' law for spheres.

The ratio of the fractal to Euclidean velocities (Γ) is shown in Figure 9(a) for different values of the aggregates' mean equivalent diameter (*d*_{f}) in the situation (*D*_{f} = 3.00 for spheres; 2.33, 2.61 and 2.83 for fractals). The ratio of fractal aggregate to sphere settling velocity (the former being obtained from fractal porosity and density considerations), varied from 0.16 to 4.11 for aggregates in the range of 10 to 1,000 μm, primary particle size of 7.5 μm and three-dimensional fractal dimension between 2.33 and 2.83. This emphasizes that fractal aggregates can behave differently, settling with higher or lower velocities, compared with Stokes' law, once settling velocities of aggregates depend upon the dual effect of porosity and density, determined by the aggregate's shape. Although large elongated flocs contain higher mass, their shape results in a lower contribution of the geometric term of Equation (22) to the settling rate than small flocs. This would only be overcome if water could flow through the flocs' pores whilst settling, as a result of macro-pore distribution within the aggregates, as mentioned by Vahedi & Gorczyca (2012). Again, there is no consensus about the permeability effect on settling rate of fractal aggregates, and despite the fact that permeability was not considered here for modelling, the findings are in agreement with a wide range of experiments, like those performed by Vahedi & Gorczyca (2012) and Johnson *et al.* (1996), for instance.

On the other hand, if primary particle size is changed, results vary for the same fractal geometry. Figure 9(b) shows the effect of primary particle size (*d*_{p}) on Γ. Results were taken following the same procedures used for Figure 9(a), and it is clear that the *d*_{p} of 1.0 μm can change the ratio of fractal aggregate to sphere settling velocity markedly to values between 0.04 and 2.92, i.e. in a lower range when compared with the *d*_{p} of 7.5 μm. For the same fractal dimension, the lower the *d*_{p}, the higher the porosity and, therefore, the lower the density, thus affecting the mass contribution over the terminal velocities of the aggregates. Further, the results showed all velocities simulated for fractal dimension of 2.33 were far lower than those predicted by Stokes' law for spheres (Γ <1), and only aggregates with size lower than 40 μm settled with higher velocities than Stokes' law for spheres when a fractal dimension of 2.61 was simulated, although, Γ was always higher than 1 for a fractal dimension of 2.83.

Whether lower or higher, simulations have shown that settling velocities of fractal aggregates can be very different from those predicted by spheres using Stokes' law. The factors affecting settling velocities of fractal aggregates depend upon size, shape, porosity, permeability, and primary particle size, and are far more complex than assumed by Stokes' law for spheres. Results presented in this paper have shown that settling velocities may vary strongly when floc shape changes from spherical to fractal aggregates, and the accuracy of predictions varies with floc size, fractal dimension and primary particle size. In practice, settling velocities of flocs are a function of floc size and fractal dimension, which are controlled by coagulation flocculation units.

## CONCLUSIONS

In this study, the porosity and density of aggregates formed after the flocculation of water containing kaolin were calculated using fractal dimension, and the terminal velocities of the Euclidean sphere and of the fractal aggregates were determined using images of a series of 118 flocs measured individually.

Settling was not used as a means to characterize aggregates' behaviour, but image analysis was used instead. Therefore, the simulations here performed were derived from aggregates' characteristics taken by image analysis, to input data into both Stokes' equation and the modified equation, based on fractal geometry. Findings are in agreement with a wide range of aggregates' traits and settling rates reported in the literature, confirming that the results are reliable.

A consistent increase of aggregate porosity with decrease of fractal dimension was observed, and the opposite was observed for the density of aggregates. Therefore, more spherical aggregates display a greater density, contributing to a higher terminal velocity of aggregates.

It was found that fractal aggregates can behave differently, settling with higher or lower velocities, compared with Stokes' law, once settling rates of aggregates depend upon the dual effect of porosity and density, determined by the aggregate's size and shape. For small floc sizes, close to the size of the primary particle (*d*_{f} ≈ *d*_{p}), the differential density between particle and liquid is the dominant effect on settling rate, surpassing the geometry reduction yielded by fractal aggregates. The opposite was observed for large fractal aggregates, where the differential density was not big enough to surpass the reduction caused by the low fractal dimension of elongated aggregates.

The results obtained differ from other work in the field by calculating settling velocities from fractal dimension and demonstrated the importance of advancing the analysis of particles considering their sizes and shapes beyond those described by Euclidean geometry. Specifically, applying fractal geometry to determine the porosity and density of the flocs is an important evaluation tool, with far-reaching implications for sedimentation tank design and operation.

## ACKNOWLEDGEMENTS

Rodrigo B. Moruzzi is grateful to São Paulo Research Foundation (Fundação de Amparo à Pesquisa do Estado de São Paulo – FAPESP) Grant 2017/19195-7 for financial support and to CNPq for the fellowship Grant 301210/2018-7.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this paper is available online at https://dx.doi.org/10.2166/wst.2020.171.

## REFERENCES

_{13}polymer