#### Abstract

The floc size distribution of activated sludge was simulated successfully by population balance model in the previous study (Population Balance Model and Calibration Method for Simulating the Time Evolution of Floc Size Distribution of Activated Sludge Flocculation. Desalination and Water Treatment, 67, 41-50). However, nonignorable errors exist in the simulation for the volume percentage of large flocs. This paper describes the application of a modified population balance model in the simulation of the time evolution of floc size distribution in activated sludge flocculation process under shear-induced conditions. It was found that the application of modified size dependent collision efficiency, modified breakage rate expression by assuming a maximum value, and binominal daughter-particles distribution function could improve the population balance model for activated sludge flocculation and successfully predict the dynamic changes in volume percentage distribution and mean floc size of activated sludge under different shear conditions. The results demonstrate that the maximum breakage rate was independent on the velocity gradient, and both the collision efficiency and breakage rate coefficient show a power-law relationship with the average velocity gradient; the former decreases while the latter increases with the rise of the average velocity gradient. These findings would help to understand the dynamics of activated sludge flocculation.

#### 1. Introduction

Activated sludge process is the most widely used technology for biologically treating wastewater, in which the separation of activated sludge from the effluent is an essential step. The formation of activated sludge flocs through flocculation plays a determinant role in achieving high separation efficiency. If flocculation fails to occur, active biomass would be lost from the system, resulting in not only excess solids being discharged into the environment, but also reduced process efficiency. Therefore, the study on the flocculation dynamics of activated sludge is significant to gain further understanding of the formation of bioflocs.

Flocculation can be regarded as a competition between aggregation and breakage. The mathematical modeling of flocculation usually makes use of the population balance model (PBM). In the PBM the flocculation dynamics are simulated by the change in floc size distribution (FSD) that is induced by simultaneous aggregation and breakage [1, 2]. The modelling of shear-induced flocculation of inorganic particles, such as polystyrene particles [3–5], latex microsphere [2, 6, 7], precipitated calcium carbonate (PCC) [8], and suspension of kaolinite or bentonite [9–12] has been widely reported. These studies focused on the prediction of the change in mean floc size and the steady-state FSD where aggregation and breakage were counterbalanced. In addition, the PBM was also used in modeling the flocculation of natural flocs or mud existing in marine or estuarine sediment ecosystems [13–15].

Different from the inorganic shear-induced flocculation processes, the flocculation of activated sludge is affected by various factors, including the floc’s physical and chemical characteristics and the external environmental conditions, such as dissolved oxygen, pH, etc. [16–18]. Consequently, discrepancies between the predicted and measured results of activated sludge flocculation have always been observed [19–22]. To achieve better modeling, the appropriate models for the aggregation and breakage of activated sludge are the key “issues.” By comparing 16 models with different size-dependent aggregation and breakage expressions, Ding et al. [21] suggested the model for the best description of the dynamic changes in the mass distribution of activated sludge under different shear. However, the relationship between the kinetic parameters and the shear rate might not follow the expected physical relationships. Nopens et al. [23] and Torfs et al. [24, 25] proposed solving the inverse problem as an alternative way to get a grip on the different expressions for aggregation of activated sludge. Although these simulation studies provided valuable findings on the failure of the existing expressions for aggregation to capture the true dynamics, selecting (or developing) appropriate kinetic expressions for aggregation or breakage is still the challenge when the PBM is used to model the flocculation of activated sludge [26].

In our previous study [22], error was found in simulating the volume percentage of large flocs (about larger than 200 *μ*m) of activated sludge flocculation by PBM, which results in smaller simulated mean sizes than their measurements. This work aims to improve the PBM to accurately simulate the time evolution of FSD in the flocculation process of activated sludge. For this purpose, the size-dependent expressions for collision efficiency and breakage rate were modified, the fractal dimension of activated sludge was introduced in modeling collision frequency and breakage rate, and a binominal function was proposed to describe the daughter-particles distribution after floc breakage.

#### 2. Flocculation Experiments

This part was presented in detail in the previous work by Li et al. [22]. Flocculation experiments were conducted in a mixing tank using aerobic activated sludge collected from a municipal wastewater treatment plant. Average velocity gradient was set in the range of 28.2~149.8 s^{−1} by adjusting the stirring speed (80~280 r·min^{−1}). FSDs of activated sludge during flocculation were measured using a laser particle size analyzer (S3500, Microtrac). Two-dimensional fractal dimension* D*_{2} of activated sludge was determined based on microscopy and image analysis. The relationship between maximum diameter (i.e., circumscribed circle diameter) (i) and equivalent diameter can be written aswhere is the primary particle size, is the area of floc in class* i*, and is the area of the primary particle.

#### 3. Model Description

The discretized population balance model developed by Kumar & Ramkrishna [28] is used to describe the variation rate in the number of flocs with a given size during activated sludge flocculationwhere is the number concentration of flocs in class* i*, is the volume of flocs in class* i*, is Dirac delta function, *α* is the collision efficiency, *β*(*j, k*) is the collision frequency of flocs in class* j *and* k*,* S*(*i*) is the breakage rate of flocs in class* i*, is the breakage distribution function defining the fraction of daughter-particle in class size* i *breaking from a floc in class size* j*, and* η*_{i} is a proportional coefficient assigning the fraction of the floc from the aggregate ():

##### 3.1. Collision Efficiency

The collision efficiency represents the probability that two flocs stay together after colliding. Generally, the collision efficiency is considered to be constant with a value of 0~1 or a function of the sizes of the two colliding particles [29]. The function proposed by Pruppacher and Klett [30], in which the sticking probability between any two particles is related to the difference of their motion speed which depends on their sizes, has been successfully used to calculate the collision efficiency of natural flocs [31]. In this study, the function was modified as follows:where is a calibration parameter.

##### 3.2. Collision Frequency

The collision frequency due to shear rate has usually been considered as a function of particle size [2, 6]. Considering that the maximum size determines the interactions between flocs more directly, the collision frequency can be calculated fromwhere* G* is the average velocity gradient,where* ε* represents the homogeneous turbulent energy dissipation rate of the mixed tank and

*is the kinematic viscosity of the suspending fluid.*

*ν*##### 3.3. Breakage Rate

The breakage of aggregates is mainly caused by hydrodynamic stress acting on the aggregates, and the breakage rate is always assumed to be relating to particle size [3, 6]. Considering that the maximum size determines the floc breakage behaviors more directly, thenwhere* E* is the breakage rate coefficient relating to the shear rate.

Eq. (8), namely, the traditional expression for breakage rate, demonstrates that the breakage rate is simply proportional to floc size. However, the Kusters [27] expression shows a different feature of the breakage rate that increases with the increasing floc size towards a final maximum [4]:where is the critical turbulent energy dissipation rate at which flocs start the breakage and is related to floc sizewhere* B* is a fitting parameter.

Taking these two expressions together into account, the traditional expression can be modified by setting a final maximum value for the breakage rate:where is the maximum breakage rate.

##### 3.4. Breakage Distribution Function

In terms of daughter-particles distribution, the simple binary breakage function and the normal breakage function are mostly used [2, 3, 5]. As a uniform computational grid, i.e., = ( is the volume of the primary particle), was used to solve the population balance model, the normal breakage function could be approximated to a binomial breakage function. Therefore, can be calculated from Eq. (12) [6, 11]:

According to Eq. (8) the higher probability of daughter-particle distribution would occur when the volume of daughter-particle is close to half of the original-particle. The normal breakage function is superior to the binary breakage function, as the former indicates a probability of daughter-particles distribution. However, the calculation load of a uniform computational grid is comparatively large while the floc size range is wide (e.g., 5~1000*μ*m). For a geometric grid with factor* k* ( = , 1 <* k *≤ 2) adopted to solve the population balance model, a binomial breakage function in another form can be used for the daughter-particle distribution (see Appendix A of Li et al. [32]):

The shape of a binomial distribution is determined by class number* j *and parameter* p*, and the probability mass function reaches its peak at the location* j*×*p*. Defining the parameter to determine the location at which the probability mass function reaches its peak, then

Assuming that the daughter-particle whose volume is the half of the original floc volume has the highest probability, can be calculated from

##### 3.5. Modelling Method

In order to investigate the applicability of modified expressions for collision efficiency, breakage rate, and binomial distribution, especially compared with that of traditional expressions (i.e., constant collision efficiency, size-dependent breakage rate, and binary distribution), four methods were used to simulate the time evolution of FSD during activated sludge flocculation (see Table 1) and their effectiveness was compared.

The initial number concentrations of the models were generated by transforming the volume distribution data coming from the laser particle size analyzer to a number distribution [22]. A geometric grid with factor 1.6 was used, i.e., ≈ 1., aligning with the Fibonacci series rule + ≈ . The binary distribution (in Table 1) was treated by using a simple technique similar to the fixed pivot technology proposed by Kumar and Ramkrishna [28]. Assuming that the floc with volume breaks up into two equal-size daughter-particles, the newly formed floc is represented by assigning proportional fractions to the flocs with volume and .

The integration of Eq. (1) involved computing the discrete Δ*N*_{i}/Δ*t* of dN_{i}/d*t *values. The iterative calculation of Δ*t* was set as 0.1~1s to maintain calculation stability. The minimum error method was used for the parameter calibration:where and are the measured and simulated volume percentage of flocs in class* i* at time* t*, respectively.* t*_{1} and are the start and end time of the flocculation process, respectively. The error of mean floc size was defined as where and are the measured and simulated volume mean floc size at time* t*, respectively, andwhere is the total volume of flocs in class* i*.

#### 4. Results and Discussion

##### 4.1. Simulation Results of FSDs

During the activated sludge flocculation process, the FSD (in the range of 5.5~1132*μ*m) initially increases and finally reaches the steady-state values after about 30 min. Method 1 and method 2 produce similar steady-state FSDs (Figure 1), and the calibrated values of parameters and* E* in the two models are very similar too (Table 2). This result might be attributed to the fact that the binomial distribution is not significantly different from the binary distribution as the binomial distribution parameter calculated from Eq. (14) results in the highest generation probability of daughter-particle whose volume equaled half of the original floc. However, the binomial distribution generated better fitting of FSD and mean size to the experimental data than the binary distribution except for* G* = 149.8s^{−1} (see* Err*_{1} and* Err*_{2} in Table 2).

**(a) G = 28.2s−1**

**(b) G = 64.7s−1**

**(c) G = 101.7s−1**

**(d) G = 149.8s−1**

In Figure 1, compared to the constant collision efficiency used in method 2, the size-dependent collision efficiency adopted in method 3 improves the simulation results of the volume percentage of small flocs (*L*<100*μ*m). However, no obvious improvement was observed for the volume percentage of large flocs (*L*>200*μ*m). Consequently, the simulated mean floc sizes were also smaller than their measurements. Method 4 with the modified expression for breakage rate gave accurate simulation to the volume percentage distributions of large flocs at steady state, and the simulated mean floc size during the flocculation agrees well with the measured results. The simulated time evolution of FSD during the flocculation by using method 4 is more consistent with the experimental data than that by using method 2, which was used in our previous study (Figure 2).

The time evolution of the mean size of activated sludge flocs at different shear rates was shown in Figure 3. The flocculation of activated sludge under shear-induced conditions is a rapid process which would be close to end after about 10 minutes. The simulated mean sizes using method 4 display a better agreement with the experimental data, whereas those obtained by using the other three methods are smaller than the experimental data. This result is related to the simulations of the volume percentage of flocs, especially for the large flocs (about* L* > 200 *μ*m). The volume percentages of the large flocs, despite being relatively small, make a significant contribution to the value of mean size calculated by Eq. (14).

**(a) G = 28.2s−1**

**(b) G = 64.7s−1**

**(c) G = 101.7s−1**

**(d) G = 149.8s−1**

##### 4.2. Parameters of Flocculation Dynamics

The simulation results demonstrate the superiority of the collision efficiency relating to the size discrepancy (Eq. (4)) over a constant collision efficiency, as well as the binominal distribution over the binary distribution for daughter-particles. The previous studies have proved that the constant values of collision efficiency negatively relate to the average velocity gradient [19]. Furthermore, the constant value of collision efficiency (in methods 1 and 2) and the average value calculated by Eq. (4) (in methods 3 and 4) both display a power-law relationship with average velocity gradient and decrease with the rise of average velocity gradient (Figure 4).

The use of the traditional expression of breakage rate in methods 1, 2, and 3 generates smaller simulated mean sizes and volume percentages of large flocs (about* L*>200*μ*m) than the experimental measurements. It might be due to the fact that the breakage rate of large flocs is overestimated by the traditional expression (Eq. (7)), in which the breakage rate is simply proportional to floc size. In fact, the Kusters [27] expression shows a different feature of the breakage rate that increases with the increasing floc size towards a final maximum. However, the Kusters expression (Eq. (10)) with the introduction of a coefficient* K* gave a better prediction for the volume percentages of large flocs (about* L*>200*μ*m), but the worse prediction for small flocs (about* L*<100*μ*m) than the traditional expression (Figure 5(a)). Moreover, it was found that, for small flocs (about* L*<100*μ*m), the Kusters expression would give smaller breakage rates than the traditional expression (Figure 5(b)). Taking these two expressions together into account, this paper modified the traditional expression by setting a final maximum value for the breakage rate (in method 4). This modification improved the simulation results of FSD and mean size obviously, and it is suggested that the supposed maximum for the breakage rate might be independent on the velocity gradient.

**(a)**

**(b)**

Generally, the relationship between the breakage rate coefficient and the average velocity gradient has been proposed to follow a power law [3, 19]. As the average velocity gradient increases, the breakage rate coefficient also increases. However, in this work, the breakage rate coefficient* E* calibrated in methods 1, 2, and 3 does not follow the relationships (Figure 6). This result might relate to the fact that the parameters are estimated completely by the best fitting to the experimental data, without the consideration of the expected physical relationships between the dynamics parameters and shear rate. Additionally, the breakage rate coefficient of method 4 exhibits an obviously positive correlation with the average velocity.

#### 5. Conclusions

This paper describes the application of PBM with modified expressions for the aggregation and breakage processes in the simulation of FSD (in the range of 5.5~1132*μ*m) in activated sludge flocculation process under shear-induced conditions (average velocity gradient of 28.2-149.8s^{−1}). The following conclusions were drawn:

More accurate simulation of FSD in activated sludge flocculation could be obtained by using size-dependent collision efficiency and a binominal daughter-particles distribution function than using constant collision efficiency and the binary daughter-particles distribution, respectively.

Compared to the traditional expression, in which the breakage rate is set to be proportional to floc size, the modified expression with a supposed maximum for the breakage rate improved the simulation results of FSD and mean size obviously.

The collision efficiency and breakage rate coefficient show a power-law relationship with the average velocity gradient; the former decreases while the latter increases with the rise of the average velocity gradient, and the maximum breakage rate might be independent with the velocity gradient.

However, it should be noted that the determination of appropriate expressions for breakage and aggregation kernels in complex systems such as activated sludge is indeed a challenge, and the applicability of the suggested collision and breakage functions to other systems remains to be tested.

#### Notation

: | Area of floc in class () |

: | Area of the primary particle () |

Parameter in Eq. (9) () | |

: | Binomial breakage function parameter |

: | Two-dimensional fractal dimension |

Breakage rate coefficient () | |

Average velocity gradient () | |

Floc equivalent diameter (μm) | |

: | Floc circumscribed circle diameter or maximum diameter (μm) |

: | Volume mean size (μm) |

: | Primary particle size (μm) |

Coefficient introduced into Kusters [27] expression for breakage rate | |

Parameter, see Eq. (13) | |

: | Number concentration of flocs in class () |

Breakage rate of flocs in class () | |

: | Floc volume of flocs in class () |

: | Volume of a fundamental particle () |

: | Total volume of flocs in class () |

*Greek Letters*

Collision efficiency | |

Collision frequency for flocs of class and () | |

: | Breakage distribution function |

: | Dirac Delta function |

Turbulent energy dissipation rate () | |

: | Critical turbulent energy dissipation rate () |

: | Coefficient, see Eq. (2) |

Kinematic viscosity (). |

#### Data Availability

The experimental data used in this paper are available on request through first author: Zhenliang Li ([email protected]).

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The work was partially supported by Natural Science Foundation of China (51609028), Natural Science Foundation of Chongqing (CSTC2016jcyjA0506), and Chongqing Postdoctoral Science Foundation (No. Xm2017141).