A new perspective on the modelling of settling behaviour in water resource recovery facilities is introduced. The ultimate goal is to describe in a unified way the processes taking place both in primary settling tanks (PSTs) and secondary settling tanks (SSTs) for a more detailed operation and control. First, experimental evidence is provided, pointing out distributed particle properties (such as size, shape, density, porosity, and flocculation state) as an important common source of distributed settling behaviour in different settling unit processes and throughout different settling regimes (discrete, hindered and compression settling). Subsequently, a unified model framework that considers several particle classes is proposed in order to describe distributions in settling behaviour as well as the effect of variations in particle properties on the settling process. The result is a set of partial differential equations (PDEs) that are valid from dilute concentrations, where they correspond to discrete settling, to concentrated suspensions, where they correspond to compression settling. Consequently, these PDEs model both PSTs and SSTs.
In the conventional treatment of wastewater, both primary settling tanks (PSTs) and secondary settling tanks (SSTs) aim at separating suspended particles from the liquid phase through sedimentation. Although both unit processes are based on the same principle, i.e. settling due to gravity, the modelling approaches have been markedly different. SSTs have been modelled by a single concentration variable for the particle phase (Takács et al. 1991; Plósz et al. 2007; Bürger et al. 2013), which means that this phase is considered as a continuum, as is the liquid phase. Hence, all particles are indirectly assumed to be identical. We call this a concentration-driven model. In contrast, discrete settling in PSTs and the clarification zone of SSTs is known to be governed by distributed settling dynamics driven by individual particle properties (e.g. size, density, shape). However, as the segregation between particles differing in size, density or shape is challenging to model, the clarification process in SSTs is generally lumped into concentration-driven models (by introducing an additional term in the hindered settling velocity function to describe low concentrations) (Takács et al. 1991) and PSTs have mostly been described by simplified models (using linear regression to relate the removal efficiency to certain characteristics of the incoming wastewater) (Amerlinck 2015).
With current focus shifting towards recovery of energy and resources, new challenges arise for operation and control of water resource recovery facilities (WRRFs). This also stimulates a renewed interest in settling as separation process as the separation of interesting fractions is an obvious first step towards their recovery. For PSTs, recent advances have led to a more detailed model based on a number of particle classes with different settling velocities (Bachis et al. 2015). This model is called the particle settling velocity distribution (PSVD) model and allows the simulation of distributed settling of the incoming wastewater particles. On the other hand, more advanced concentration-driven models for SSTs have gained increased attention (Kinnear 2002; Plósz et al. 2007; De Clercq et al. 2008; Bürger et al. 2013; Li & Stenstrom 2014). The focus here has been on the incorporation of compression (i.e. the resistance to hindered settling by the network of flocculated particles that arises at high concentrations) as several studies have shown that hindered settling alone does not capture the complex settling behaviour of activated sludge (De Clercq et al. 2005; Ramin et al. 2014; Torfs et al. 2015). Although compression settling is known to depend on the concentration gradient, which introduces a second-order term in the governing partial differential equation (PDE), the exact relation describing this compressive behaviour is still unknown. A number of studies proposed expressions to describe compression from experimental data (De Clercq et al. 2008; Ramin et al. 2014; Diehl 2015). However, these analyses evidenced that identification of a compression function in a purely concentration-driven model was not possible as the physics of compression cannot be modelled only in dependence of the concentration X (and its gradient) with a constant parameter set. Some approaches introduce an empirical variation in the transition concentration between hindered and compression settling (the so-called critical concentration Xcrit) (De Clercq et al. 2008; Ramin et al. 2014; Locatelli 2015) but a physical explanation for this variability has not yet been provided.
The first part of this contribution provides new experimental evidence that the unexplained variability in thickening and compression behaviour in SSTs has a similar origin as the variability observed in PSTs and the clarification zone of SSTs. Variations in distributed particle properties such as size, shape, porosity, density and flocculation state leading to a distributed settling velocity are shown to have an important influence on the settling behaviour in different settling unit processes (and at different concentrations throughout these unit processes). Hence, although each settling process is still governed by its own specific dynamics, these dynamics can be attributed to a common source, thus calling for a unified framework to describe the different settling unit processes. The experimental evidence presented in this first part provides the background and motivation for the second part of this contribution which addresses the development of a unified settling framework. The conceptual steps to extend current models into a framework that can be applied to different settling unit processes are presented. Finally, in a third part, the potential of the new proposed framework is illustrated with simulation examples.
EXPERIMENTAL EVIDENCE IN SUPPORT OF A UNIFIED FRAMEWORK
Discrete settling in PSTs and the clarification zone of SSTs
Hindered and compression settling in the thickening zone of SSTs
Figure 3 (right) shows how the settling behaviour of a sludge sample changes as it is subject to different amounts of shear stress prior to settling. The settling velocity at the top of the sludge blanket changes noticeably between the different experiments. Since each of these tests was performed at the same initial concentration, the differences in settling behaviour cannot be attributed to variations in concentration but only to variations in particle properties. By stirring the sample, the flocculation state of the sludge is changing as loosely bound flocs are broken up into more stable aggregates. This process may be further facilitated by the release of extracellular polymeric substances acting like a polymer to increase flocculation (Laurent et al. 2009). The resulting aggregates are characterised by better settling properties leading to a faster decrease in the sludge blanket height. Stirring the sample at high shear rates also decreased the sludge's ability to act as a filter for colloids, causing an increase in supernatant turbidity.
Moreover, Figure 3 (right) illustrates the impact of variations in particle properties with respect to the onset of compression settling. For suspensions with sufficiently low critical concentration Xcrit (such as activated sludge), the bend in the batch settling curve corresponds to the point where the sludge blanket enters the compression zone. (Note that for suspension of, for example, hard spherical particles, that undergo hindered settling only, such a bend may also occur, but this does not apply here.) Hence, for the curves in Figure 3 (right) the onset of compression is indicated by the grey dotted lines. At these points, the concentration at the top of the sludge blanket should equal the critical concentration. When shear is applied prior to settling, the sludge water interface reaches the compression zone earlier and, more importantly, at a much lower sludge blanket height (and thus in a more concentrated state). A more concentrated sludge blanket at the onset of the compression zone signifies a higher critical concentration. This indicates that differences in distributed particle properties potentially account for variations in the critical concentration. Due to the applied shear, larger and less stable flocs will be reduced to smaller and denser flocs with different packing properties resulting in a sample that can reach higher concentrations before the particles are in permanent contact.
The combined results of Figure 3 indicate that similar to settling at low concentrations, changes in size, shape and density should also be considered as important factors to describe the variability in hindered and compression settling. This finding has important implications for the modelling of settling processes as it emphasizes the need to step away from purely concentration-driven SST models and opens up new perspectives for a unified settling framework which can be applied to all settling unit processes.
New applications of a unified model framework
A unified framework which can capture the true variability in the settling behaviour would allow expanding the use of settler models beyond their current applications. For conventional SSTs this would not only allow improved predictions of effluent suspended solids but also to model the impact of chemical addition and increased hydraulic loadings. For PSTs, this would allow to model compression of the particles accumulated at the bottom to better predict the underflow concentration (and consequently the feed concentration to an anaerobic digester).
DEVELOPMENT OF A UNIFIED SETTLING MODEL FRAMEWORK
Current state of the art settling models
From the examples above, distributed properties such as particle size, shape and density were shown to be key factors in describing the settling processes in WRRFs. Existing secondary settling models such as the traditional Takács model (Takács et al. 1991) and the more recent Bürger-Diehl framework (Bürger et al. 2013) do not account for distributed behaviour of particles; they simply describe the solids by means of a lumped concentration variable. Moreover, neither of these settling models consider changes in particle properties due to aggregation and break-up processes. Empirical relations describing the effect of polymer addition on hindered settling have been described in literature (Vanderhasselt et al. 1999) but these do not account for the underlying changes in particle property distributions. Recently, a new model for PSTs, called the PSVD model (Bachis et al. 2015), was presented. The PSVD model does include distributed settling behaviour and has been extended to include the effect of chemically enhanced primary treatment (CEPT) by alum dosage (Tik et al. 2016). However, this model only covers the discrete settling regime and can therefore not be applied to SSTs.
In order to achieve a unified description of the processes taking place both in PSTs and SSTs, either existing PST models such as the PSVD have to be extended with hindered and compression dynamics or existing SST models need to be extended to include particle classes and their associated distributed behaviour. In this contribution, it was chosen to start from the existing Bürger-Diehl framework for SSTs (Bürger et al. 2013) as this framework already includes hindered and compression settling in its underlying PDE and most importantly an appropriate numerical scheme to solve this PDE. The Bürger-Diehl framework can be further extended by combining knowledge from several existing approaches. These include the PSVD model in PSTs but can also be found outside WWT modelling such as models for polydisperse sedimentation (Berres et al. 2003) and population balance models (PBMs) describing the mechanisms of flocculation and breakage (Nopens et al. 2015). The remainder of this contribution provides the conceptual steps to extend the Bürger-Diehl framework into a unified framework for the settling dynamics in different settling unit processes.
The Bürger-Diehl secondary settler framework
A number of classes needs to be introduced to represent the distributed properties of the sludge.
These different classes allow to extend the existing framework with distributed settling behaviour such as discrete settling and to account for the effect of distributed particle properties on hindered and compression dynamics.
The model equations can also be augmented with reaction terms to describe changes in particle property distributions due to operational and loading conditions. These different steps are introduced below.
Extension to several particle property classes
A well-flocculated sludge will have a larger concentration of particles in the large/dense/fast class and a lower concentration of particles in the class of small/open/slow particles.
Instead of a single nonlinear convection-diffusion PDE, this will result in a system of PDEs, one for each particle property class. Such a set of PDEs has been presented in literature (Berres et al. 2003) for classes of particles having different sizes and densities. In its most general treatment, each particle class has its own settling velocity function which depends on the concentrations of all classes and their spatial derivatives, making it quite complex. Therefore, in this work, a simpler approach is followed, where only a number of specific dependencies required to capture the dynamics of settling in WRRFs are added.
Specification of distributed settling behaviour
At low concentrations such as those occurring in the PST or in the clarification zone of the SST, sludge particles typically undergo discrete settling. As the discrete settling velocity is considered independent of concentration, each particle class will settle at its own characteristic velocity (v0,i).
The parameter represents the transition concentration between discrete and hindered settling.
By substituting Equation (10) in Equation (7), the corresponding set of PDEs is valid from dilute concentrations , for which they coincide with discrete settling models, to concentrated suspensions in which the sediment is described as a permanently networked, compressible porous layer. It can thus be used to model both PSTs and SSTs.
The distributed properties of the sludge will not only cause distributed settling behaviour but will also influence the transition between the different settling regimes (characterised by Xtrans and Xcrit). For example, aerobic granular sludge is known to have a low tendency to coagulate under reduced hydrodynamic shear (de Kreuk & van Loosdrecht 2004). This feature causes granular sludge to undergo discrete settling at concentrations where conventional activated sludge experiences hindered or compression settling. The presented framework can be applied to granular sludge by simply setting an appropriately high transition concentration Xtrans.
In this function Xcrit,i corresponds to the critical concentration that would arise in a suspension of mono-sized particles of class i. The total critical concentration of a mixture of different particle classes is calculated (at each height in the settling column and each simulation time) as a weighted sum of Xcrit,i according to the concentration of each particle class present. The reasoning behind this definition stems from the idea that faster settling particles (for example the HiCS sludge in Figure 4) would have a less porous structure thus allowing for a denser packing before being subjected to the force of compression. As such, the proposed modelling framework is able to describe the unexplained variability in the critical concentration (reported by De Clercq et al. (2008) and Ramin et al. (2014)) through changes in distributed particle properties. A similar approach can be applied to define the transition between discrete and hindered settling (Xtrans) as a weighted function of the concentrations in the different classes.
Specification of flocculation processes
The reaction terms ri(Qf, Xi, X,Cchem) can be based on flocculation jar experiments (Gong et al. 2011) or can be derived from PBMs which allow detailed modelling of the dynamics of distributions (Nopens et al. 2015).
An alternative and simple approach would be to assume that flocculation/break-up is mostly occurring in the flocculation well. In this case, the flocculation well can be described as a continuous stirred tank reactor (CSTR) prior to the actual settling tank where flow dependent flocculation and break-up processes give rise to a certain particle class distribution that can subsequently be fed as input to the one-dimensional (1D) settler model. This approach would remove the need to add reaction terms to each PDE. However, it is only valid for well-designed clarifiers where flocculation is indeed (approximately) limited to the flocculation well.
This section presents a simulation example to illustrate the specific features of the proposed unified model framework emphasizing its ability to describe variability in settling behaviour that cannot be attributed to changes in concentration only. The model equations were implemented for batch sedimentation only (no bulk flow nor incoming feed flow), and flocculation between classes is not considered thus far (ri= 0). The example includes all three settling regimes (discrete, hindered and compression settling) with discrete and hindered settling velocity functions given by Equation (8).
Introducing different particle classes in a settler model gives rise to a number of coupled second-order PDEs whose discretization might pose significant restrictions on the time step of the solution. Moreover, due to the nature of settling as a separation process, simulation results will give rise to sharp discontinuities in the concentration profiles of all classes. To handle these issues, an IMEX scheme (Boscarino et al. 2015) was used for the numerical implementation of this system. In this scheme, first-order convection terms are discretised in an explicit way and second order diffusive terms are handled in a linearly implicit way ensuring an efficient solution of the proposed model framework. All details concerning the application of an IMEX scheme to the presented model framework can be found in Bürger et al. (in preparation).
Experimental evidence of sludge settling at different concentrations and under different conditions showed that variations in distributed particle properties (such as size, shape, porosity and density leading to a distributed settling velocity) are an important factor influencing the settling behaviour in all settling unit processes in WRRFs. Hence, current modelling practice where the settling behaviour only depends on a lumped variable such as concentration is insufficient to describe the true settling dynamics in WRRF processes.
Therefore, an extension of existing modelling frameworks is proposed, including different particle classes to represent the distribution in particle properties. The result is a unified framework which allows description of the true distributed settling behaviour over the entire concentration range from dilute suspensions, where discrete settling occurs (e.g. for PSTs, granular sludge and clarification in SSTs), to concentrated suspensions in which the sediment is described as a permanently networked, compressible porous layer (e.g. at high concentrations in the sludge blanket of SSTs or at the bottom of PSTs). Moreover, the proposed framework can be further augmented with reaction terms to describe the migration of particles between different classes in order to capture the effect of changes in operational and loading conditions on the sludge composition and the associated settling behaviour. Ultimately, this would allow introducing more rigour into the way settling tanks are modelled, potentially leading to improved operation and control in WRRFs.
Raimund Bürger is supported by Fondef project ID15I10291; CRHIAM, Proyecto Conicyt Fondap 15130015; Fondecyt project 1130154; and BASAL project CMM, Universidad de Chile and CI2MA, Universidad de Concepción. M. Carmen Martí is supported by CONICYT through Postdoctoral 2015 Fondecyt project 3150140. Sophie Balemans is supported by a PhD scholarship from the Research Foundation Flanders (FWO Vlaanderen, number SB-151737). Peter Vanrolleghem is supported by the Natural Sciences and Engineering Research Council of Canada (NSERC/CRSNG), John Meunier and Québec City and holds the Canada Research Chair on Water Quality Modelling. This paper was presented at WWTmod2016 and the fruitful discussions are kindly acknowledged. The authors also acknowledge the support of the Strängs Donationsfond.