A Numerical Study on Influent Flow Rate Variations in a Secondary Settling Tank

: The secondary settling tank is an essential unit for the biochemical treatment of domestic sewage, and its operational effect influences the quality of effluent. Under the influence of the confluence of rainwater and sewage, wastewater use habits, etc., the inflow of the secondary sedimentation tank changes over time. In this paper, OpenFOAM, an open-source computational fluid dynamics package, was used to study the dynamic behaviors of solid–liquid two-phase flow in the tank under influent flow rate variations. A coupled method including a mixture model, drift equation and a dynamic boundary method is proposed. Numerical investigations were carried out for a 2D axisymmetric sedimentation tank using 12 cases. With increasing influent flow rate, sludge accumulates continuously in the bottom left side of the tank, sludge hopper, and inlet; the sludge blanket thickness near the right end of the tank increases continuously; and the sludge concentration in the tank approximately linearly increases with time, with a low slope. The developed framework is generic and is, therefore, expected to be applicable for modelling sludge sedimentation processes.


Introduction
Secondary settling is referred to in the literature as "the most sensitive and complicated process in activated sludge plants" [1]. As an essential biological treatment unit in municipal wastewater treatment plants, the secondary settling tank (hereinafter called "the tank") restricts the wastewater treatment capacity [2,3]. The sludge sedimentation influences the quality of effluent. Therefore, it is important to study the sludge sedimentation law and its influencing factors in order to improve sludge sedimentation efficiency.
Some new methods have been used in studying sewage treatment systems, such as artificial neural networks [4,5], numerical simulation [6][7][8], and smoothed particle hydrodynamics [9]. In particular, numerical simulation technology has gradually become an important method used to study the sludge sedimentation mechanism and solid-liquid dynamical behavior in the tank [10][11][12]. Different numerical simulation techniques including the mixture model [13], the Navier-Stokes coupled sludge transport equation [10,14], and the Euler-Euler two-phase model [7], have been used to study sludge sedimentation. The results revealed the effects of the geometric structure, temperature distribution, settling velocity, and turbulence model on sludge sedimentation; further, the solid-liquid dynamics behavior and heterogeneous distribution of activated sludge concentration in the tank were discussed [6,8,[15][16][17][18]. Gao and Stenstrom [19] adopted different turbulence models (standard k-ε , realizable k-ε , renormalization group k-ε ) to model the dynamic behavior of solidliquid in the tank. They found that different turbulence modes have little influence on the prediction of the suspended solids concentration at the outlet but a significant impact on the sludge blanket height. Xu et al. [7] developed a new force-based mechanical model in the Euler-Euler framework based on Newton's law rather than the flux theory and used the new model to describe the sludge settling process. They found that the predictions of this model were more accurate than those of flux theory-based models. Ramin et al. [13] developed a new settling velocity model, including hindered, transient, and compression settling, and showed that it could be calibrated to data from a simple, novel settling column experimental setup using the Bayesian optimization method differential evolution adaptive Metropolis (DREAM). In addition, correlations between the Herschel-Bulkley rheological model parameters and sludge concentration were identified with data from batch rheological experiments. Byonghi [11] studied the double perforated baffles in a rectangular secondary settling tank by numerical simulation. He found that the double perforated baffles could improve the quality of effluent water; when the flow rate was large, biological flocculation appeared in the front of the baffles, hindering the longitudinal movement of sludge in the tank. Gao and Stenstrom [20] studied the effects of baffling structures on buoyancy-induced turbulence flow in secondary settling tanks. The results showed that the effects of the secondary settling tank (SST) baffling structures on the buoyancy-induced turbulence can be divided into three zones: Zone A, where the baffling structures had negligible damping effects, and only the buoyancy-coupled turbulence model provided accurate predictions; Zone B, where the baffling structures partially dissipated the effect of buoyancy on turbulence, and the buoyancy-coupled model continued to provide accurate predictions but the buoyancy-decoupled turbulence model only provided qualitative but similar predictions; and Zone C, where the baffling structures fully dissipated the effect of buoyancy on turbulence, and both models provided similar, accurate predictions. Valle Medina and Laurent [6] focused on the development of a computational fluid dynamics (CFD) tool including not only the hydrodynamics of the settling tank but also recent advances in both hindered and compression settling modeling. Then, they studied the effects on the sludge blanket height and sludge concentration using the compression function. Gao and Stenstrom [21] investigated the effects of the buoyancy term (Gb) in the turbulent kinetic energy (TKE) equation and two model parameters on the performance of secondary settling tanks. The results showed that the hydrodynamics can only be correctly predicted by a buoyancy-coupled TKE equation.
However, the above studies were usually based on steady conditions and seldom focused on sludge sedimentation behaviors in secondary settling tanks under influent flow rate variations. While affected by wastewater use habits, water consumption, confluence of rainwater and sewage, etc., the tank operates daily under unsteady operation conditions. Therefore, the flow rate of mixed liquid at the inlet of the secondary settling tank varies with time [22]. The dynamic influent flow rate of the mixed liquid at the inlet not only changes the velocity distribution in the tank but also correspondingly adjusts the turbulence energy and turbulent dissipation rate of the mixed liquid. Meanwhile, the convection and diffusion effect of the mixed liquid is changed. As a result, the sludge sedimentation in the tank is affected by the dynamic flow rate. Therefore, steady numerical results cannot possibly reveal the sludge sedimentation mechanism under dynamic conditions. It is very difficult to reproduce and predict the spatial and temporal evolution of the sludge distribution in the tank. Although it will not need more computational time to simulate the sludge settling process in the tank under dynamic influent flow, the boundary conditions need to consider the influence of time. This means that new boundary conditions need to be adopted. This paper focuses on the liquid-solid dynamic behavior in the secondary settling tank under influent flow rate variations. Furthermore, the influence of reflux conditions is taken into account in the numerical simulation. Different from the previous numerical simulation studies on solid-liquid two-phase flow in a secondary settling tank under steady flow, this study reveals the sludge settling mechanism with respect to influent flow rate variations. The modeling result is expected to aid in the optimization of operating parameters of the tank.

Numerical Model
OpenFOAM has been widely used in numerical simulations in multiphase hydrodynamic behavior research [23,24]. It has also been adopted as a platform to study the solid-liquid behavior in secondary settling tanks [6,13,25].

Mixture Model
The mixture model is a commonly used mathematical model for the numerical simulation of secondary settling tanks, and it takes into account the two-phase slip velocities.
The mixture continuity equation is: where ρ is density; v is velocity; and the subscripts m represent mixture. The mixture velocity is obtained as follows (Equation (2)): where α is the volume fraction of the two phases; the subscripts d and w represent solid and liquid, respectively. There is a relationship between the liquid and solid volume fractions as shown in Equation (3); the mixture density is described in Equation (4).
The mixture momentum equation is expressed as: where v s is the sludge settling velocity and can be calculated via Equation (11); τ is the viscous stress tensor and is described by Equation (7); t τ is the turbulent stress tensor; g is gravitational acceleration; and Pm is pressure.
The drift equation for the prediction of the solids distribution is: where Γ is the diffusion coefficient.
The standard k-ε turbulence equation was used to consider turbulent effects [25,26].

Rheological Model
The Herschel-Bulkley model is used to describe the shear stress [13,[27][28][29], as shown below: where 0 τ is the yield stress and can be calculated via Equation (8); K is the viscosity coefficient and is expressed as Equation (9); n is the rheological index and is described by Equation (10); and γ is the shear rate. However, the sludge concentration distribution in the tank is not uniform. The rheological properties of the mixed liquid in the tank are related to the concentration. According to Ramin's study, an activated sludge sample was collected from the combined recycle flow channel of the Lundtofte wastewater treatment plant on the day of the settling experiments. Sludge samples taken were aerated overnight with coarse air bubbles and the sampled effluent SST was stored at 4 ℃ for the rheological experiments conducted on the following day. The three parameters ( 0 τ , K, n) were estimated by fitting the Herschel-Bulkley model to the rheological measurements at the shear rates between 0.01 and 250 s −1 [13]. The parameters of different sludge concentrations (5, 7.1, 8.5, 10, 12.8 kg/m 3 ) were investigated in a simple laboratory set-up, consisting of a glass column with a Solitax ® total suspended solids sensors [13]. In order to further study the real rheological properties of the mixed liquid, according to the reference by Ramin et al. [13], the following formulae were obtained by fitting the shear stress formula with the concentration in this paper: 57; X is the sludge concentration; and w μ is the viscosity coefficient of water.

Settling Velocity
The hindered settling velocity model by Takacs et al. [30] was used in the numerical simulation. It is also commonly adopted as the settling velocity model to investigate sludge sedimentation in the tank [22,31]. Further, a series of batch settling tests were performed using the sludge samples diluted with secondary settling tank effluent over a concentration range of 1.7-5.1 kg/m 3 , which is typical for feed concentration into secondary settling tanks [13]. Then, the measurement parameters for the sludge rheological properties and settling velocity were analyzed. The formula is as follows: where X is the sludge concentration at different locations in the tank; X0 is the initial sludge concentration at the inlet; and rH and rp are 0.42 and 4.2, respectively. Figure 1 shows the geometrical structure and mesh of the 2D axisymmetric circular sedimentation tank which was used to study dynamic behaviors of solid-liquid two-phase flow in references [32,33]. The secondary settling tank studied is located in Lundtofte wastewater treatment plant (Lyngby, Denmark). It is a circular center-feed conical tank [13]. As shown in Figure 1a, the tank has a length of 12.25 m, height of 6.2 m, inlet width of 0.225 m, outlet width of 0.08 m, and reflux outlet width of 0.58 m. In order to express the material balance in the tank more clearly, it is assumed that the flow rate of the mixed liquid at the inlet is Qin; the sludge concentration at the inlet is Xin; the flow rate at the outlet is Qout; the sludge concentration at the outlet is Xout; the flow rate at the reflux outlet is R*Qout, which is related to the flow rate at the outlet; and the sludge concentration at the reflux outlet is Xrec. Further, Qin = Qout + R*Qout. Figure 1b shows the unstructured mesh of the whole computational domain. As shown in Figure  1b, the mesh is denser around the inlet area, sludge hopper, and outlet. The minimum face area is 0.0002597 m 2 , and the maximum face area is 0.0659 m 2 . In this paper, a grid independence study was performed to ensure the results of independent mesh resolution. The total number of cells in the numerical simulation is 7088.

Parameter Setting
The chosen parameter values were as follows: water density, 1000 kg/m 3 ; sludge density, 1600 kg/m 3 ; gravity acceleration, 9.8 m/s 2 ; flow rate at the inlet, 242 m 3 /h; and R, 1.42. The PIMPLE algorithm was adopted to decouple velocity and pressure. In the numerical simulation, at a time step, the correct number was 4. The relaxation factor of each physical variable was set to 0.7.

Case Setting
It was assumed that the increase of the influent flow rate is linear and that the sludge concentration at the inlet does not change with time. All the cases were modeled for 1800 s under steady conditions (with flow rate at the inlet of 242 m 3 /h and flow rate at the outlet of 100 m 3 /h). As a result, we obtained relatively stable solid-liquid dynamic behavior in the tank. Then, five different cases were used to model the sludge sedimentation under different influent flow rates from 1.5 to 4 during the time from 1800 to 5400 s. At each time step, a program was added to alter the influent flow in a mandatory way to the one at the next time point. This ensured that the influent flow rate increased linearly with time. In addition, the dynamic reflux ratio is also discussed for the cases. The relevant settings for the cases are shown in Table 1.

Ideal Sedimentation Tank
The solid-liquid behavior was simulated in an ideal sedimentation tank in order to verify the accuracy of the above numerical method and numerical model. The ideal sedimentation tank has a length of 8.65 m, width of 1 m, inlet width of 1 m, outlet width of 0.1 m, velocity at the inlet of 0.0191 m/s, and solid volume fraction of 0.001. Brennan described the facilities, methods, and steps of the experiment in detail [34,25]. According to Dahl's study [25], measurements of velocity and dispersed phase concentration were taken along the center plane of the test section at distances of 3.0 m, 5.0 m and 7.0 m from the inlet end wall. Integration of the velocity profiles to give the volumetric flux through each station at the experimentally recorded times shows variation from the claimed influent volume. The mixture velocity was measured by an electro-magnetic induction sensor and the suspended sludge concentration was measured by an optical turbidity meter. Figure 2a shows the structured meshing of the ideal sedimentation tank. The flow rate and solid volume fraction at different locations in the ideal sedimentation tank were given. Figure 2b shows a comparison of the mixture (i.e., mixture of sludge and water) velocity at 3 m, 5 m, and 7 m in the x direction. As shown in Figure 2b, it was found that the mixture velocity was high in the middle of the tank and low at both sides of the tank. The simulated mixture velocity value agreed well with the experimental value.

Dynamic Flow Verification
The influent flow rates of Cases 1-7 were investigated and compared with the analytical values of the flow rate designed in different cases at 2400 s, 3000 s, 3600 s, 4200 s, and 4800 s. It can be seen in Figure 3 that the influent flow rate variations with time were identical to the analytical values at different times, which proves that the algorithm of the boundary conditions guarantees dynamic change of the flow rate at the inlet with time.  Figure 4 shows the streamline in different cases in the tank at 5400 s. Generally speaking, the sludge settling effect proceeds well when the length of the recirculation zone in the secondary sedimentation tank is long [35]. The actual flow field in the tank has complicated structures, usually three-layer vortex structures, which is consistent with the findings of Van and Kranenburg [36] during monitoring of the flow field in an actual secondary settling tank. The reason for this may be that the different sludge concentrations show different rheological properties at different depths in the tank. The upper-layer vortex mainly influences the clarification area above the sludge blanket; the middle-layer vortex affects the flow pattern of the mixture in the middle and lower parts of the tank; and the lower-layer vortex affects the flow pattern of the mixture at the bottom of the tank. In the upper and middle parts of the sludge hopper, there is a main control vortex which affects the flow of the mixture in most areas of the sludge hopper. A small vortex distributes at the bottom of the sludge hopper, which restricts the transport of the sludge from the bottom to the middle part of the sludge hopper. As shown in Figure 4a-e, comparing Cases 1-5, both the number and shape of the vortices in the upper part of the tank change significantly, and those in the middle and lower parts of the tank change insignificantly. Moreover, the flow of the mixture is still controlled by a narrow and long vortex in the lower part of the tank. In Figure 4f,g, when the influent flow rate remains unchanged, there is little difference in the vortex structure in the tank, i.e., a certain difference in the vortex-affected area extending only between the upper right side and the middle and lower parts of the tank. As can be seen from Figure 4h Figure 5 shows the sludge concentration distributions in the tank at 5400 s in different cases. The sludge concentration in the upper part of the tank was very low. With increasing depth in the tank, the sludge concentration increased continuously and reached a peak at the bottom of the sludge hopper. As shown in Figure 5a-e, comparing Cases 1-5, as the influent flow rate increased from 1.5 to 4, the distribution area of sludge concentration of >1.12 g/L widened around the inlet (Block 1 area). Kim et al. [37] hypothesized that the density effect was caused by increased sludge concentration inlet loading on the overall flow field. The sludge constantly accumulated, and the sludge concentration constantly increased at the bottom left side of the tank (Block 2 area), while the sludge blanket constantly thickened near the right side of the tank (Block 3 area); in addition, the distribution area of high sludge concentration (>8.87 g/L) increased at the bottom of the sludge hopper (Block 4 area). This was caused by the following factors: (1) when the mixture flowed into the tank, the sludge particles gradually settled down due to gravity and accumulate. Then, the gap between sludge particles gradually decreased. Consequently, the sludge concentration increased in the tank. (2) Increase in the influent flow rate led to an increase in the total sludge quantity in the whole tank; little sludge was discharged through the outlet, most of the sludge settled in the sludge hopper, and part of the sludge was discharged through the reflux outlet. As a result, the sludge concentration increased and the sludge accumulated continuously in the sludge hopper area. (3) Due to the high viscosity of sludge, the mixture had poor fluidity leading to easier sludge settlement at the bottom of the tank. (4) The flow dead-zone affected the mixture flow at the bottom of the sludge hopper and tank. Therefore, the sludge concentration increased in these areas.

Sludge Concentration Distribution
As shown in Figure 5f,g, the steady flow rate at the inlet in Case 6 (322.7 m 3 /h) was higher than that in Case 7 (242 m 3 /h). Additionally, under steady inflow, the results of Case 6 were consistent with the result by Guyonvarch et al. [27]. Consequently, the sludge concentration in the bottom right side of the tank increased significantly. Cases 7 and 8 showed significantly broader sludge concentration distributions at the bottom of the sludge hopper (sludge concentration >8.87 g/L) than did Cases 1-5. The reason for this may be that the vortex is more stable under steady flow. Thus, it was conducive to sludge sedimentation. As shown in Figure 5h,i, the sludge distribution in Cases 8-12 was similar to that in Cases 1-5; Cases 8-12 showed slightly higher sludge concentration around the bottom right side than did Cases 1-5 if the change in the reflux (Cases 8-12) was not considered. This may be due to larger sludge discharge from the reflux outlet in Cases 8-12 than in Cases 1-5 and, thus, lower total sludge quantity in the whole tank. For a quantitative analysis of the average spatial sludge concentration in the tank, see Figure 6.  Figure 6 shows the average spatial sludge concentrations in different parts of the tank. As shown in Figure 6e, the tank can be divided into three parts-the upper part called Zone 1, the bottom called Zone 2, and the sludge hopper area called Zone 3. As shown in Figure 6a, the different cases all had lower average sludge concentration in Zone 1 than in the other zones. Among them, Case 7 had the lowest average sludge concentration in Zone 1; Cases 2-4 and Cases 9-11 had similar sludge concentration values. However, Cases 5 and 12 had much higher sludge concentrations than the other cases. As can be seen from Figure 6b, different cases had higher average sludge concentrations in Zone 2 than the cases for Zone 1; the average sludge concentrations in Cases 1-5 and 8-12 increased with increasing influent flow rate. It was found by comparing Cases 1 and 8, Cases 2 and 9, Cases 3 and 10, Cases 4 and 11, and Cases 5 and 12 that the dynamic change settings of the flow rate at the reflux outlet did not significantly influence the sludge concentration in Zone 2 and only slightly reduced the average sludge concentration. As shown in Figure 6c, in Zone 3, different cases had the biggest average sludge concentrations in the tank. The fluctuation of sludge concentration in the different cases was smaller in Zone 3 than in other zones. As shown in Figure 6d, the average sludge in the whole tank was a bit higher in Cases 1-5 than that in Cases 8-12. The changes in the sludge concentration in Zone 2 and Zone 3 were similar in different cases. However, the sludge concentration in Zone 3 was lower than that in Zone 2 in different cases. In summary, the sludge concentration in Zone 2 greatly influenced the average sludge concentration in the tank, and the dynamic change settings of the flow rate at the reflux outlet influenced the average sludge concentration in the tank very little.  Figure 7 shows the change in the average spatial sludge concentration in the tank with time. As shown in Figure 7a-c, the sludge concentration in the tank in Cases 1-12 approximately linearly increased with time, with low slope. Among them, Case 5 had the highest average spatial sludge concentration over time, followed by Case 12. The average sludge concentration in the tank was mainly determined by the sludge input, not the dynamic reflux. A comparison of Cases 5 and 1 at 5400 s indicated that the average sludge concentration in Case 5 was 1.47 times that in Case 1. However, the influent flow rate in Case 5 was 2.67 times that in Case 1. That means that much more sludge was discharged through the reflux outlet in Case 5.   Figure 8a, the sludge concentrations in Cases 1-5 increased in a fluctuant manner with time. This is mainly because (1) the flow field of the mixture in the tank remained in a non-steady state and (2) the monitoring site had low sludge concentrations, leading to easy convection and diffusion of the sludge. However, the increase of the influent flow rate with time caused the sludge concentration to increase in this area.
In Figure 8b, the sludge concentrations at the monitoring point increased in a fluctuant manner with time in Cases 6 and 7. Besides this, the sludge concentrations in Case 6 were higher than those in Case 7 at 2400 s, 2900 s, 3400 s, 3900 s, 4400 s, 4900 s, and 5400 s. This is because the steady-state flow rate at the inlet in Case 6 was larger than that in Case 7. As shown in Figure 8c, Cases 8-12 were different in terms of their sludge concentrations at varying time from Cases 1-5, which means that the change in the reflux flow rate influenced the instantaneous sludge concentration in the tank.  Figure 9 shows the velocity distributions in the tank in different cases. The fluid flow behaviors were very complicated in the tank. Among the 12 different cases, the upper part of the tank had lower velocity; the velocity in the lower part was higher than that in the middle part, and the high velocity area (>0.0549 m/s) was distributed from the inlet to the reflux outlet. This is mainly because (1) the mixture at the inlet flowed downwards due to the viscous effect of the fluid and gravity induced sedimentation of sludge in the tank, and (2) one stream of the mixed liquid flowed directly to the reflux outlet along the sludge hopper due to the special geometric structure of the area above the sludge hopper; thus, insufficient attenuation of the flow velocity leaded to larger velocity in the area. Due to the larger velocity distributed in the sludge hopper area, the mixture directly flowed out without complete sedimentation. Therefore, it was not conducive to the efficient utilization of the tank. In addition, the velocity between the middle and lower parts of the tank was low in different cases, which was beneficial for sludge sedimentation in this area.

Velocity Distribution
As shown in Figure 9a-e, comparing Cases 1-5, with increasing influent flow rate, the area with velocity of over 0.0166 m/s expanded at the bottom of the tank. In particular, an obvious S-shaped velocity (>0.0166 m/s) distribution formed between the inlet and reflux outlet. In addition, with increasing influent flow rate, the high velocity area (>0.0166 m/s) was distributed in the right wall area of the tank. As a result, the liquid entrained sludge and moved to the upper outlet, which was not beneficial for sludge settlement in the area. As shown in Figure 9f,g, the velocity at the bottom and middle of the tank in Case 6 was larger than that in Case 7. In Figure 9h,i, Cases 8-12 showed similar velocity distributions to Cases 1-5. The dynamic change of the flow rate at the reflux outlet had little impact on the velocity distribution.  Figure 10 shows the average spatial velocities of different cases in different zones. The average spatial velocities in Cases 5 and 12 were similar but higher than those in other cases. The average velocity in Zone 1 was larger than that in other zones. That is because the sludge concentration in Zone 1 is low and the mixture is very flowable. The average velocity in Zone 3 was lower than that in other zones. The reason for this is that high sludge concentration causes a decrease in the fluidity of the mixture. As shown in Figure 10a, with increasing influent flow rate, the average velocity in the different cases increased in Zone 1, especially in Cases 5 and 12. As can be seen from Figure 10b, the average velocity increased with increasing influent flow rate. As shown in Figure 10c, the average velocity in Cases 1-5 increased more than that in Cases 8-12 because the sludge concentration in Cases 8-12 was high. As shown in Figure 10d, the average velocity in the different cases was much lower in the whole area than in Zones 1 and 2 but higher than in Zone 3.  Figure 11 gives the changes in velocity at the monitoring sites (10 m, 3 m) with time. As shown in Figure 11a, the velocity at the monitoring sites was between 0 and 0.012 m/s and fluctuated with time in Cases 1-5. As shown in Figure 11b, although Case 6 had a much bigger influent flow rate than Case 7, the velocities at 2400 s, 2900 s, 4900 s, and 5400 s were relatively low. As shown in Figure 11c, similar to Cases 1-5, the velocity in Cases 8-12 decreased first, then increased, and finally decreased again. Because the sludge constantly settles down, the flow characteristics of the mixture changed accordingly.

Sludge Settling Velocity Distribution
In order to reveal the different sludge concentration distributions in different cases, the sludge settling velocities at different flow rates are shown in Figure 12. Low sludge settling velocity (>−0.0003 m/s) was distributed in the upper part of the tank and close to the bottom of the tank and the bottom of the sludge hopper; the sludge settling velocity was higher in the bottom and middle parts of the tank. Ramin et al. [13] also found the same research phenomenon in their simulations. From the hindered settling velocity model by Takacs et al. [30], we can explain the result: (1) the sludge concentration was very low, and the sludge was suspended in the upper part of the tank; (2) the type of sludge settling was flocculent or free settling in the bottom and middle parts of the tank, so the sludge settling velocity was higher. As shown in Figure 12a-e, in Cases 1-5, low sludge settling velocity (>−0.0003 m/s) was mainly distributed in the upper part and bottom of the tank, as well as in the bottom of the sludge hopper. The sludge settling velocity in the upper part of the tank was lowest because of the minimal sludge concentration and sludge suspension in this area. However, as the sludge settled down in the bottom of the tank and sludge hopper area, the gap between sludge particles continuously decreased; thus, the sludge possibly settled down under gravity at a higher rate. It was found by comparing Cases 1-5 that the settling velocity in the middle and upper parts of the sludge hopper decreased continuously, meaning that incomplete sludge sedimentation in the sludge hopper area led to a higher sludge concentration value at the reflux outlet. As shown in Figure  12f,g, Cases 6 and 7 had a similar settling velocity distribution to Cases 1-5. In Figure 12h,i, with increasing influent flow rate, the high-value area (>0.0015 m/s) of the sludge settling velocity at the bottom of the tank constantly extended in the x direction but significantly decreased in the y direction.

Conclusions
In this paper, we studied the dynamic behaviors of solid-liquid two-phase flow in a tank under influent flow rate variations using a method coupling a mixture model and drift equation while considering the sludge rheological properties due to the sludge concentration. Meanwhile, a dynamic boundary method was developed to describe changes of the flow velocity with time. Moreover, validation of the inflow changes, velocity, and phase volume was given by comparing with experimental or analyzed values.
With increasing influent flow rate, the vortex structure in the upper part of the tank changed dramatically, while the vortex structures in the middle and bottom of the tank were relatively stable. The sludge accumulated continuously in the bottom left side of the tank, sludge hopper, and inlet; the sludge blanket thickness near the right end of the tank increased continuously; and the sludge concentration in the tank approximately linearly increased with time, with low slope. The settings of the dynamic flow rate at the reflux outlet had a limited influence on the sludge concentration, sludge distribution, mixed liquid velocity distribution, and sludge settling velocity distribution in different areas of the tank.
In actual sewage treatment processes, it is necessary to avoid too fast a flow increase with time, as it leads to a rise in the sludge blanket height and higher sludge concentration in the effluent. In addition, in order to reduce the flow velocity at the inlet, a baffle plate is set around the inlet of the tank. Multiple regulating tanks can be set up to keep the inlet flow stable.
Further studies on the modeling of wastewater treatment processes in secondary sedimentation tanks should consider the dynamic inlet sludge concentration. The biological and physical properties of the sludge also affect the settlement characteristics. Therefore, it is necessary to further study the settling velocity and rheological properties of sludge.