Anaerobic sewage treatment systems, notably the upflow anaerobic sludge blanket reactor, have been applied successfully in Brazil in dozens of large scale units. In this paper a small reactor is proposed for use in single family houses, without a connection to a sewerage network. The proposed reactor is much smaller than the traditional septic tank and yet the treatment efficiency is much higher. It comprises a lower cylindrical digestion zone and an upper settling zone, connected by a transition zone. Several configurations with different digestion and settling zone diameters were tested to establish their influence on reactor performance. A high degree of organic matter removal was obtained with 250 L reactors with retention times of 6 to 12 hours. A large sludge mass was retained in the digestion zone and high treatment capacity was maintained throughout the study. Tests showed that sludge settled better in the reactor than it did after discharge with the effluent. In contrast the, methanogenic activity of the sludge in the reactor was the same as that of the sludge expelled.

In regions without a sewerage system, on-site systems must be applied to treat domestic sewage. The main objectives of these activities are to protect public health and the environment, especially water resources, although water and nutrient reuse may also be important (USEPA 2005; Aiyuk et al. 2006; Jorsaraei et al. 2014; Ladu & Lü 2014). The septic tank (ST) has been the most widely applied small-scale sewage treatment system since the last decade of the nineteenth century, even though performance is poor compared to more modern systems (Coelho et al. 2003). This can be attributed to serious ST design flaws, as shown in the schematic in Figure 1. The most important failure is that the sewage passes horizontally over the sludge layer at the bottom of the unit, so that contact between the bacterial sludge mass and incoming organic material is poor. Only the settlable fraction is transported for digestion into the sludge layer (Metcalf and Eddy 2014). ST treatment efficiency can be increased by improving two basic design aspects: (1) ensuring efficient contact between the biological mass within the treatment system and the incoming organic material and (2) maximizing the sludge mass in the treatment system by some form of sludge retention.
Figure 1

Septic tank (ST).

Figure 1

Septic tank (ST).

Close modal
The upflow anaerobic sludge blanket (UASB) reactor has been applied successfully for sewage treatment in regions with high average temperatures, where mesophilic anaerobic digestion can develop (Lew et al. 2011). Conventional UASB reactors (Figure 2) have three zones: the lowest, which receives the influent at the bottom, is the digestion zone. The topmost part is a tranquil zone where suspended solids can settle. Between the two is a transition zone. A characteristic feature of conventional UASB reactors is the three phase separator, which collects the biogas from the digestion zone, and retains and returns solids escaping into the settling zone, so that an effluent substantially free of suspended solids is discharged at the top. Good contact between the sludge mass in the reactor and incoming organic matter is assured by the upflow of liquid through the sludge in the digestion zone, while sludge retention is guaranteed by the settler in the upper zone (Van Haandel & Lettinga 1994; Abassi & Abassi 2012; Chong et al. 2012).
Figure 2

Conventional UASB reactor.

Figure 2

Conventional UASB reactor.

Close modal
Coelho et al. (2003) proposed an adapted UASB reactor as a single family anaerobic treatment unit (Figure 3). This design is characterized by two two-phase separators rather than the single three-phase separator found in conventional UASBs. The design incorporates an inclined side arm through which the effluent flows to discharge, and solid-liquid separation occurs here. Gas is separated from the liquid and solids in the vertical upper part of the reactor. This design was shown to perform significantly better than conventional UASBs because settlement is superior. Even slow-settling particles are retained, and tend to coagulate into larger flocs that settle more readily and can slide down to become part of the biological sludge again.
Figure 3

UASB reactor with two two phase separators.

Figure 3

UASB reactor with two two phase separators.

Close modal

Van Haandel & Lettinga (1994) showed that other anaerobic systems tend to perform less well than the UASB reactor. Almost all large anaerobic treatment systems constructed in Brazil since 1990 are UASB reactors or variants. However, in small-scale units, UASB reactors are less common, although Coelho et al. (2003) showed that a 360 L reactor (as in Figure 3) is superior to the ST/anaerobic filter combination. In the latter case the ST is used mainly to retain gross solids that would otherwise clog the upflow anaerobic filter, which itself provides efficient removal of organic matter from the pre-treated effluent from the ST. In this paper a UASB reactor configuration is proposed that can substitute the traditional ST and be used as single family alternative for decentralized sewage treatment.

Design of the single family units

Building on the design proposed by Coelho et al. (2003), the aim was to develop a single family unit that would provide good contact between incoming organic matter and sludge, efficient sludge retention, and a small volume, as well as sufficient capacity to absorb the sudden hydraulic loads that characterize sewage production in single family households. A balance was sought between cost considerations, which tend to reduce unit volume, and treatment efficiency, which tends to require relatively larger volumes, for sludge storage. All reactors comprised three zones: lower digestion and upper settling zones, both cylindrical, with a conical transition zone between them (see Figure 4). The cylindrical form was chosen as it is easy and inexpensive to make, with materials like polyvinyl chloride or fiber glass.
Figure 4

Proposed Single Family UASB reactor (UASB – SF).

Figure 4

Proposed Single Family UASB reactor (UASB – SF).

Close modal

The digestion and settling zone diameters must be chosen so that an adequate upward velocity is maintained within those zones. Since this ‘adequate velocity’ was not known at first, six reactors of different diameters were constructed and tested (see Table 1). The important digester parameters were calculated once the basic parameters (HRT, and the digestion, transition and settling zone diameters, as well as the separator's inclination) were defined. Table 2 shows the expressions used to define the various dimensions. The diameter of the conical separator was 0.05 m larger than that of the digestion zone, so that there was an overlap and biogas bubbles could not enter the settling zone.

Table 1

Expressions used to calculate the dimensions of the different parts of the six UASB reactors (Figures 4 and 5)

LineBasic parametersSymbolT1T2T3S1S2S3Unit
Treated flow rate Qa m3.d−1 
Volume Vtot 0.25 0.25 0.25 0.25 0.25 0.25 m3 
Hydraulic Retention Time HRT 
Diam. digestion zone Ddi 0.30 0.30 0.30 0.40 0.40 0.40 
Diam. settling zone Dde 0.60 0.70 0.80 0.65 0.70 0.75 
Inclination of separator α 45 45 45 45 45 45 ° 
Tangent tanα 1.00 1.00 1.00 1.00 1.00 1.00  
Diam. separator Ds 0.35 0.35 0.35 0.45 0.45 0.45 
Free board Hb 0.05 0.05 0.05 0.05 0.05 0.05 
 Derived parameters         
10 Digestion area Adi 0.071 0.071 0.071 0.126 0.126 0.126 m2 
11 Settling area Ade 0.283 0.385 0.502 0.332 0.385 0.442 m2 
12 Velocity in digester zone. vdi 0.59 0.59 0.59 0.33 0.33 0.33 m.h−1 
13 Velocity in settling zone Vs 0.15 0.11 0.08 0.13 0.11 0.09 m.h−1 
14 Overlap 0.025 0.025 0.025 0.025 0.025 0.025 
15 Transition height Htr 0.150 0.200 0.250 0.125 0.150 0.175 
16 Separator height Hse 0.175 0.175 0.175 0.225 0.225 0.225 
17 Settler height Hde 0.125 0.125 0.125 0.175 0.175 0.175 
18 Volume settling zone Vde 0.035 0.048 0.063 0.058 0.067 0.077 m3 
19 Volume transition zone Vtr 0.025 0.041 0.063 0.028 0.036 0.047 m3 
20 Volume digestion zone Vdi 0.190 0.161 0.124 0.164 0.146 0.126 m3 
21 Volumetric proportion in digester zone fdi 0.760 0.643 0.495 0.658 0.585 0.504 – 
22 Volumetric proportion in transition zone ftr 0.099 0.165 0.254 0.110 0.146 0.187 – 
23 Volumetric proportion in settling zone fde 0.141 0.192 0.251 0.232 0.269 0.309 – 
24 Digester height Hdi 2.689 2.274 1.753 1.310 1.165 1.003 
25 Total height Htot 3.014 2.649 2.178 1.659 1.539 1.403 
28 Max. velocity in settling zone vmax 0.223 0.144 0.103 0.241 0.185 0.147 m.h−1 
29 Min. velocity in settling Zone vmin 0.147 0.108 0.083 0.126 0.108 0.094 m.h−1 
LineBasic parametersSymbolT1T2T3S1S2S3Unit
Treated flow rate Qa m3.d−1 
Volume Vtot 0.25 0.25 0.25 0.25 0.25 0.25 m3 
Hydraulic Retention Time HRT 
Diam. digestion zone Ddi 0.30 0.30 0.30 0.40 0.40 0.40 
Diam. settling zone Dde 0.60 0.70 0.80 0.65 0.70 0.75 
Inclination of separator α 45 45 45 45 45 45 ° 
Tangent tanα 1.00 1.00 1.00 1.00 1.00 1.00  
Diam. separator Ds 0.35 0.35 0.35 0.45 0.45 0.45 
Free board Hb 0.05 0.05 0.05 0.05 0.05 0.05 
 Derived parameters         
10 Digestion area Adi 0.071 0.071 0.071 0.126 0.126 0.126 m2 
11 Settling area Ade 0.283 0.385 0.502 0.332 0.385 0.442 m2 
12 Velocity in digester zone. vdi 0.59 0.59 0.59 0.33 0.33 0.33 m.h−1 
13 Velocity in settling zone Vs 0.15 0.11 0.08 0.13 0.11 0.09 m.h−1 
14 Overlap 0.025 0.025 0.025 0.025 0.025 0.025 
15 Transition height Htr 0.150 0.200 0.250 0.125 0.150 0.175 
16 Separator height Hse 0.175 0.175 0.175 0.225 0.225 0.225 
17 Settler height Hde 0.125 0.125 0.125 0.175 0.175 0.175 
18 Volume settling zone Vde 0.035 0.048 0.063 0.058 0.067 0.077 m3 
19 Volume transition zone Vtr 0.025 0.041 0.063 0.028 0.036 0.047 m3 
20 Volume digestion zone Vdi 0.190 0.161 0.124 0.164 0.146 0.126 m3 
21 Volumetric proportion in digester zone fdi 0.760 0.643 0.495 0.658 0.585 0.504 – 
22 Volumetric proportion in transition zone ftr 0.099 0.165 0.254 0.110 0.146 0.187 – 
23 Volumetric proportion in settling zone fde 0.141 0.192 0.251 0.232 0.269 0.309 – 
24 Digester height Hdi 2.689 2.274 1.753 1.310 1.165 1.003 
25 Total height Htot 3.014 2.649 2.178 1.659 1.539 1.403 
28 Max. velocity in settling zone vmax 0.223 0.144 0.103 0.241 0.185 0.147 m.h−1 
29 Min. velocity in settling Zone vmin 0.147 0.108 0.083 0.126 0.108 0.094 m.h−1 

HRT: hydraulic retention time; Diam.: diameter; Vel.: velocity; Vol.: volume; Dig.: digestion; Min.: minimum; Max.: maximum; Set.: settling.

Table 2

Expressions defining aspects of the UASB configurations used in the study

Line n°Derived parametersSymbolFormula
10 Digestion area   
11 Settling area   
12 Velocity in digestion zone vdi =Qa/Adi 
13 Velocity in settling. Zone Vs =Qa/Ade 
14 Overlap =(Ds–Ddi)/2 
15 Transition height Htr =(Dde–Ddi)/2*tanα 
16 Separator height Hse =Ds/2*tanα 
17 Settler height Hde =Hse – Hb 
18 Settler volume   
19 Transition volume   
20 Digestion volume Vdi =Vtot–Vde–Vtr 
21 Volumetric proportion in digestion zone fdi =Vdi/Vtot 
22 Volumetric proportion in transition zone ftr =Vtr/Vtot 
23 Volumetric proportion in settling zone fde =Vde/Vtot 
24 Digester height Hdi =Vdi/Adi 
25 Total height Htot =Hdi + Htr + Hde + Hb 
26 Minimum area in settling zone.   
27 Maximum area in settling zone.   
28 Maximum velocity in settling zone vmax =Qa/Amin 
29 Minimum velocity in settling zone vmin =Qa/Amax 
Line n°Derived parametersSymbolFormula
10 Digestion area   
11 Settling area   
12 Velocity in digestion zone vdi =Qa/Adi 
13 Velocity in settling. Zone Vs =Qa/Ade 
14 Overlap =(Ds–Ddi)/2 
15 Transition height Htr =(Dde–Ddi)/2*tanα 
16 Separator height Hse =Ds/2*tanα 
17 Settler height Hde =Hse – Hb 
18 Settler volume   
19 Transition volume   
20 Digestion volume Vdi =Vtot–Vde–Vtr 
21 Volumetric proportion in digestion zone fdi =Vdi/Vtot 
22 Volumetric proportion in transition zone ftr =Vtr/Vtot 
23 Volumetric proportion in settling zone fde =Vde/Vtot 
24 Digester height Hdi =Vdi/Adi 
25 Total height Htot =Hdi + Htr + Hde + Hb 
26 Minimum area in settling zone.   
27 Maximum area in settling zone.   
28 Maximum velocity in settling zone vmax =Qa/Amin 
29 Minimum velocity in settling zone vmin =Qa/Amax 

The upward velocities in the digestion and settling zones are the basic parameters, when designing the configuration shown in Figure 4, needed to ensure good contact in the digestion zone and good sludge retention in the settling zone. Figure 4 and Table 1 show that, for sewage flow, Qa, the UASB-SF reactor is defined by (1) the reactor volume, Vu, or retention time HRT, (2) the transition zone and separator inclination (α), and (3) the digestion (Dd), settling zone (Ds), and separator (Dp) diameters. These define all other reactor dimensions and were operated at constant flow rates of 0.5 and 1 m3.d−1, corresponding to the minimum and maximum anticipated flow rates for a single family unit. In a second phase, variable flow rates were applied as significant variations in sewage flow rate with respect to time tend to happen in single family units.

The basic parameters in Table 1 influence the values of the derived parameters and must be chosen such that the values of the latter are adequate. The most important variables are:

  • (1) The liquid velocity in the digestion and settling zones. The velocity in the digestion zone must allow retention of a large sludge mass, but be high enough to avoid inorganic solids accumulation there. Velocities of less than 0.3 m.h−1 lead to grit and sand accumulation and are thus too low. The lower the velocity at the top of the settling zone, the more efficient the sludge retention, but an evaluation is needed to determine whether the retention of hard-to-settle particles and their return to the digestion zone is really advantageous, as it might affect sludge settleability. The maximum upward fluid velocity at the bottom of the settling zone is dictated by the difference between the settling zone and phase separator diameters. If it is low, it is easier to return retained particles to the digestion zone.

  • (2) The proportional volume of the different zones. The optimal division of volumes of the settling and digestion zones is not obvious: a large digestion zone allows a large sludge mass, but a large settling zone is required for efficient sludge retention and hence a high sludge concentration in the digestion zone, although it is the settling area, not its volume, that determines particle retention. Nevertheless, the settling volume can be important in promoting flocculation of small particles into better settling flocs. The division of the volume over the settling and the digestion zone was a key point of investigation.

  • (3) Reactor height. A tall reactor has a small footprint, but its implementation may be complicated. In order to avoid pumping, the top of the reactor must be below the level at which the sewage arrives at the treatment unit, but digging deep may be difficult, especially in areas with rocky soil.

It is impossible to determine the optimal values of the basic variables by theoretical analysis. Thus, the performance of a series of reactors with equal volumes but different dimensions was evaluated. Figure 5 shows the experimental set-up for the investigation, with six UASB-SF reactors (see Figure 4), a Y-type reactor (Figure 3) and a conventional unit (UASB-C, Figure 2). All were constructed using fiber glass and had a volume of 250 L, which was assumed to be sufficient for efficient treatment. A smaller volume was considered unlikely to reduce costs significantly but might lead to problems during construction, installation and operation. The inclination of the transition zone and the phase separator was 45° in all cases, as this is usual for UASB reactors, although no specific research into angle optimization could be found.
Figure 5

Experimental pilot plants: three UASB reactors with digestion zone diameters of 0.4 m (numbers 1, 2 and 3) and three UASB reactors with 0.3 m (numbers 4, 5 and 6), a conventional UASB (7) and a UASB-Y reactor (8). Legend: 1- Reactor S1 – Digestion zone diameter 0.4 m; Settler diameter 0.75 m; 2- Reactor S2– Digestion zone diameter 0.4 m; Settler diameter 0.70 m; 3- Reactor S3 – Digestion zone diameter 0.4 m; Settler diameter 0.65 m; 4- Reactor T1 - Digestion zone diameter 0.3 m; Settler diameter 0.8 m; 5 - Reactor T2 - Digestion zone diameter 0.3 m; Settler diameter 0.7 m; 6 - Reactor T3 - Digestion zone diameter 0.3 m; Settler diameter 0.6 m; 7 - UASB conventional reactor, digestion zone diameter 0.3 m; 8- UASB–Y reactor, digestion zone diameter 0.3 m, with two two phase separators; 9 - Distribution tube for influent sewage; 10 - Dosing pumps.

Figure 5

Experimental pilot plants: three UASB reactors with digestion zone diameters of 0.4 m (numbers 1, 2 and 3) and three UASB reactors with 0.3 m (numbers 4, 5 and 6), a conventional UASB (7) and a UASB-Y reactor (8). Legend: 1- Reactor S1 – Digestion zone diameter 0.4 m; Settler diameter 0.75 m; 2- Reactor S2– Digestion zone diameter 0.4 m; Settler diameter 0.70 m; 3- Reactor S3 – Digestion zone diameter 0.4 m; Settler diameter 0.65 m; 4- Reactor T1 - Digestion zone diameter 0.3 m; Settler diameter 0.8 m; 5 - Reactor T2 - Digestion zone diameter 0.3 m; Settler diameter 0.7 m; 6 - Reactor T3 - Digestion zone diameter 0.3 m; Settler diameter 0.6 m; 7 - UASB conventional reactor, digestion zone diameter 0.3 m; 8- UASB–Y reactor, digestion zone diameter 0.3 m, with two two phase separators; 9 - Distribution tube for influent sewage; 10 - Dosing pumps.

Close modal

The biological, chemical and mechanical properties of the sludge that developed in the reactors are also important. The most important biological sludge properties are its activity or methane production potential, and its stability (the fraction of biodegradable organic material it contains). The most relevant chemical property is the fraction of inorganic sludge, while its settleability is the most significant mechanical property. All were evaluated for each reactor as a function of operational conditions.

All reactors had sludge sampling points at several heights. Table 3 shows their levels from bottom to top. The total sludge mass was calculated from the sum of the products of the volume of each reactor segment corresponding to a sampling point and the sludge concentration measured there.

Table 3

Sample point levels from bottom to top for the different reactors (m)

PointsT1T2T3S1S2S3YC
P1 0.20 0.20 0.20 0.10 0.10 0.10 0.20 0.20 
P2 0.40 0.40 0.40 0.25 0.25 0.25 0.40 0.40 
P3 0.60 0.60 0.60 0.40 0.40 0.40 0.60 0.60 
P4 1.00 1.00 1.00 0.70 0.60 0.63 1.00 1.00 
P5 1.40 1.40 1.30 0.97 0.83 – 1.40 1.40 
P6 1.80 1.60 1.60 – – – 1.80 1.80 
P7 2.50 2.10 – – – – – 2.20 
P8 – – – – – – – 2.60 
HT 3.3 2.7 2.0 1.5 1.35 1.15 2.8 3.5 
Dd 0.3 0.3 0.3 0.4 0.4 0.4 0.3 0.3 
Areaa 0.071 0.071 0.071 0.126 0.126 0.126 0.071 0.071 
PointsT1T2T3S1S2S3YC
P1 0.20 0.20 0.20 0.10 0.10 0.10 0.20 0.20 
P2 0.40 0.40 0.40 0.25 0.25 0.25 0.40 0.40 
P3 0.60 0.60 0.60 0.40 0.40 0.40 0.60 0.60 
P4 1.00 1.00 1.00 0.70 0.60 0.63 1.00 1.00 
P5 1.40 1.40 1.30 0.97 0.83 – 1.40 1.40 
P6 1.80 1.60 1.60 – – – 1.80 1.80 
P7 2.50 2.10 – – – – – 2.20 
P8 – – – – – – – 2.60 
HT 3.3 2.7 2.0 1.5 1.35 1.15 2.8 3.5 
Dd 0.3 0.3 0.3 0.4 0.4 0.4 0.3 0.3 
Areaa 0.071 0.071 0.071 0.126 0.126 0.126 0.071 0.071 

HT: total height; Dd: diameter digestion zone.

aArea: Cross section area of the digestion area.

Wastewater and operating conditions

Raw sewage from Campina Grande, Paraíba (Brazil), was used in the study. The reactors were initially operated at constant flow and load conditions for retention times of 12 and 6 hours. The flows were 500 and 1,000 L.d−1, respectively, which were taken as representing the sewage from one and two, four-person, households (125 L.inh−1.d−1). It was expected that these flows could be treated easily, as Coelho et al. (2003) design had been shown to do so efficiently.

During operation of the reactors there was no intentional sludge wastage and the sludge mass built up to its maximum value. Subsequently, sludge was expelled from the reactor at a rate equal to that of its generation. Thus sludge production was determined from the expelled sludge mass in the effluent. The volatile sludge concentration in the effluent was taken to be the difference in VSS concentration in the raw effluent (as it was being released from the reactors) and the settled effluent (after settling the raw effluent for 30 min in an Imhoff cone). Now, knowing the sludge mass in the reactors and the daily volatile sludge expulsion or production (the product of the settleable VSS concentration in the effluent and the flow rate), the sludge age is calculated as:
formula
1
where:
  • Rs = Sludge age

  • MXv = Volatile sludge mass in the reactor (calculated from the concentration profile)

  • Qa = Sewage flow

  • Xve = Volatile sludge concentration in the effluent (the difference of VSS concentrations in the raw and settled effluent).

In Equation (1) the sludge age denotes the value calculated on the basis of the volatile solids. The sludge age on the basis of total sludge can also be calculated but the values of MXv and Xve must be substituted by MXt and Xte, where the index ‘t’ denotes the total solids content of the sludge.

The sludge age in a UASB reactor is known to be its fundamental operating parameter (Santos et al. 2016), determining reactor performance with respect to biogas and sludge production. For any sludge age, the influent chemical oxygen demand (COD) is divided into three fractions: (1) the non-settleable fraction in the effluent, (2) the settleable COD fraction in the effluent, which represents the influent COD fraction that is transformed into sludge, and (3) the COD fraction converted into methane, which is the difference between unity and the sum of the non settleable and settleable COD fractions in the effluent. The objective of anaerobic treatment is that the COD fractions in the effluent and transformed into sludge be as small as possible and consequently the digested COD fraction is maximum. The COD fractions can be determined experimentally from the COD concentrations in the raw and in the settled effluent and the influent as follows:
formula
2a
formula
2b
formula
2c
where:
  • mSe = Non-settleable COD fraction in the effluent

  • mSx = Settleable COD = COD fraction transformed into biological sludge

  • mSd = Digested COD fraction

  • Sse = Non-settleable effluent COD

  • Sbe = Raw effluent COD = total or settleable + non-settleable effluent COD

  • Sta = Influent COD

Analytical procedures

Tests to characterize the performance of the UASBs were carried out over 18 months. Their objective was to determine the division of influent COD into its fractions in the effluent, the sludge and biogas. COD and suspended solids tests of influent and effluent were carried out following procedures described in Standard Methods for the Examination of Water and Wastewater (APHA et al. 2012). COD determinations were done on raw and settled effluent samples, and the COD fraction transformed into sludge was taken as the difference between the two values.

The specific methane production potential (SMPP) is the methane production per unit mass of sludge per unit time unit, under ideal anaerobic digestion conditions. It was determined in accordance with the procedure described by Van Haandel & Lettinga (1994), using 250 mL bottles (effective volume 200 mL) at a constant temperature of 35 °C. Solutions were added to ensure the presence of macro- and micro- nutrients, and that pH = 7 was maintained during the test. A substrate mass of 0.8 g COD was added as acetate solution, (4 g-COD.l−1), capable of producing about 200 mL of methane at 0 °C and 1 bar (standard conditions). The anticipated duration of the test was 3 to 5 days. Cumulative methane production was determined from the volume of solution displaced in a Mariotti flask, and SMPP was determined from the maximum methane production rate and calculated as g-COD.gXv−1.d−1. Substitution of a mixture of acetate and propionate for the acetate feed did not affect the methane production rate, whence it was concluded that biogas production in the treatment units was due to acetotrophic rather than hydrogenotrophic methanogens. SMPP tests were carried out with sludge samples withdrawn from the reactor and samples obtained by settling the raw effluent to obtain the sludge that was expelled from the reactor. Thus it was possible to determine if the SMPP of the sludge in the reactors was the same as that of the sludge expelled from the reactors.

Sludge settleability was evaluated in terms of the equation proposed by Vesilind (1967) to express the sludge settling velocity empirically (Equation (3)):
formula
3
where:
  • v = Sludge settling velocity

  • Xt = Sludge concentration

  • vo, k = Vesilind's settleability constants

The constants k and v0 are linked to the mechanical properties of sludge: v0 reflects the sludge settling velocity of particles in a sludge so dilute that the particles settle individually (unhindered settling). The constant k describes sludge compressibility and the higher its value, the more difficult it is to obtain a high concentration by settling. Vesilind's equation was originally proposed for the determination of activated sludge settling, but tests carried out in this study showed that it can also be used for anaerobic sludge such as generated in the UASB reactors. The Vesilind constants k and vo were determined for all reactors, again, not only with samples taken from the reactors but also with sludge that was expelled from the reactors.

Reactor performance

The performance of all reactors was evaluated in terms of influent organic matter removal efficiency as well as digestion efficiency. In this context, digestion efficiency is given by Equation (2c). The removal efficiency is the fraction of the influent COD that is found in the settled effluent and is expressed by Equation (2a). The digestion efficiency is the fraction of influent COD that is converted to methane and is expressed by Equation (2c). Removal efficiency is given by Equation (2a) plus the fraction transformed into settleable solids (Equation (2b)). Table 4 shows the organic matter removal efficiency and sludge retention data for the reactors operating at 12 and 6 hour HRTs.

Table 4

Characterization of COD removal and sludge retention for the different reactors operated at 12 and 6 h HRTs

VariableSymbUnitT1T2T3S1S2S3RYRC
HRT = 12 h 
 Influent COD Sta mg.l−1 587 587 587 805 805 805 587 587 
 Effluent COD (sol) Sre mg.l−1 191 169 155 244 248 251 185 189 
 COD to sludge Sx mg.l−1 216 116 97 51 48 43 128 141 
 Removal efficiency – 67 71 74 70 69 69 68 69 
 Digestion efficiency – 31 50 57 61 63 61 47 44 
 Volumetric Load rate – g/l.d−1 1.2 1.2 1.2 1.7 1.7 1.7 1.2 1.2 
HRT = 6 h 
 Influent COD Sta mg.l−1 460 460 460 460 460 460 – – 
 Effluent COD (sol) Sre mg.l−1 139 141 166 154 138 159 – – 
 COD to sludge Sx mg.l−1 91 47 73 74 43 101 – – 
 Removal efficiency – 70 69 64 67 70 65 – – 
 Digestion efficiency – 52 59 48 49 61 45 – – 
 Volumetric Load Rate – g/l.d−1 1.9 1.9 1.9 1.9 1.9 1.9 – – 
VariableSymbUnitT1T2T3S1S2S3RYRC
HRT = 12 h 
 Influent COD Sta mg.l−1 587 587 587 805 805 805 587 587 
 Effluent COD (sol) Sre mg.l−1 191 169 155 244 248 251 185 189 
 COD to sludge Sx mg.l−1 216 116 97 51 48 43 128 141 
 Removal efficiency – 67 71 74 70 69 69 68 69 
 Digestion efficiency – 31 50 57 61 63 61 47 44 
 Volumetric Load rate – g/l.d−1 1.2 1.2 1.2 1.7 1.7 1.7 1.2 1.2 
HRT = 6 h 
 Influent COD Sta mg.l−1 460 460 460 460 460 460 – – 
 Effluent COD (sol) Sre mg.l−1 139 141 166 154 138 159 – – 
 COD to sludge Sx mg.l−1 91 47 73 74 43 101 – – 
 Removal efficiency – 70 69 64 67 70 65 – – 
 Digestion efficiency – 52 59 48 49 61 45 – – 
 Volumetric Load Rate – g/l.d−1 1.9 1.9 1.9 1.9 1.9 1.9 – – 

Table 4 reveals some important aspects:

  • (1) The COD removal efficiency (COD transformed into methane or sludge) was influenced only marginally by the configuration in the range used, apart from design T1 whose performance was inferior to that of the others at both 6 and 12 h HRTs. The influence of hydraulic and organic loads on removal efficiency was very limited in the investigated range: reducing HRT from 12 to 6 h led to an insignificant decrease. Similarly the anaerobic digestion efficiency (COD fraction converted to methane) was also affected little by configuration or HRT.

  • (2) The COD fraction converted into settleable sludge was about half of the non-settleable total. Thus it is very important to state whether raw or settled effluent was used to determine effluent COD.

  • (3) The non-biodegradable content of soluble material in the influent is atypically high for municipal sewage, which must be attributed largely to the fact that the sewerage outfall was partially blocked and conveyed only a fraction of the generated sewage. Because of this, extensive anaerobic digestion was taking place in the sewerage network, reflected in a low influent COD concentration, with high non-biodegradable soluble and low non-biodegradable particulate fractions, due to settling in the sewer.

  • (4) The soluble effluent COD consisted largely of non-biodegradable organic matter. This was apparent from the biochemical oxygen demand/COD ratio, which was less than 0.25 in all tests and had an average value of 0.17.

  • (5) The volatile fatty acids (VFA) concentration was consistently low at less than 0.5 mol.l−1.

  • (6) As reactor configuration was shown to have little influence on treatment efficiency, other factors will determine the preferred design. These include the total reactor height (difficulty of installation, especially in rocky terrain) and instability caused by shock loads. Reactors with low upward fluid velocities, which tend to be shallow, are affected by shock loads less than those with relatively higher velocities.

  • (7) All reactors had much smaller footprints than equivalent STs, which is a clear advantage in housing schemes with small plots. The Brazilian norm for the minimum area of ST is 1.5 m2, which is 3 to 5 times larger than the reactors that were operated in this investigation.

Characterization of sludges in reactors

Sludge concentration profile and mass

In Table 5 the total and volatile sludge masses in each reactor are calculated from the concentration profiles. Also included are the averages of daily expelled sludge mass (total and volatile), and the sludge ages calculated from the ratio of the retained and expelled masses, for the different reactors for both total and volatile sludge, by using Equation (1).

Table 5

Sludge masses and ages in and expelled from different reactors, calculated for constant flow rates of 500 and 1.000 l.d−1

VariableT1T2T3S1S2S3YC
Influent flow rate = 500 l.d1 
 Total sludge (g−17,300 9,938 6,465 4,123 6,084 4,933 10,626 10,817 
 Volatile sludge (g.l−13,763 5,120 3,205 2,379 3,512 2,775 4,790 5,218 
 Expelled total sludge (g.l−10.16 0.103 0.082 0.099 0.128 0.11 0.114 0.136 
 Expelled volatile sludge (g.l−10.116 0.074 0.056 0.072 0.102 0.062 0.067 0.076 
 Total sludge age(d−1129 308 380 136 176 216 333 208 
 Volatile sludge age (d−196 221 284 107 127 217 258 182 
Influent flow rate = 1,000 l.d1 
 Total sludge (g−12,578 4,282 5,952 4,789 4,115 4,520 – – 
 Volatile sludge (g.l−11,493 2,338 3,302 2,640 2,221 2,163 – – 
 Expelled total sludge (g.l−10.130 0.110 0.119 0.208 0.110 0.141 – – 
 Expelled volatile sludge (g.l−10.094 0.051 0.091 0.121 0.056 0.090 – – 
VariableT1T2T3S1S2S3YC
Influent flow rate = 500 l.d1 
 Total sludge (g−17,300 9,938 6,465 4,123 6,084 4,933 10,626 10,817 
 Volatile sludge (g.l−13,763 5,120 3,205 2,379 3,512 2,775 4,790 5,218 
 Expelled total sludge (g.l−10.16 0.103 0.082 0.099 0.128 0.11 0.114 0.136 
 Expelled volatile sludge (g.l−10.116 0.074 0.056 0.072 0.102 0.062 0.067 0.076 
 Total sludge age(d−1129 308 380 136 176 216 333 208 
 Volatile sludge age (d−196 221 284 107 127 217 258 182 
Influent flow rate = 1,000 l.d1 
 Total sludge (g−12,578 4,282 5,952 4,789 4,115 4,520 – – 
 Volatile sludge (g.l−11,493 2,338 3,302 2,640 2,221 2,163 – – 
 Expelled total sludge (g.l−10.130 0.110 0.119 0.208 0.110 0.141 – – 
 Expelled volatile sludge (g.l−10.094 0.051 0.091 0.121 0.056 0.090 – – 

The results in Table 5 show that there was a considerable reduction in sludge mass in the system when the influent flow increased from 0.5 to 1 m3.d−1. Reactors T1, T2 and T3 (digestion zone diameter 0.3 m) lost about 70, 40 and 10%, respectively, of the initial mass. This can be attributed to (1) the inoculum used in these reactors, which was taken from an anaerobic pond with high concentrations of inert and inorganic matter, and (2) reactors T1 and T2 had relatively small settlers, so that the sludge could not be retained when the flow rate was increased. Reactors S1, S2 and S3, like UASB-Y, were inoculated with sludge from T1, T2 and T3, and this conventional sludge mass was much less affected by flow increases, because the inoculum was better and they had larger settling areas.

The sludge ages obtained by applying Equation (1) show that the reactors with bigger sedimentation zones had better biomass retention, expelling less solids in the effluent while maintaining a longer sludge age. Other research has shown that there is no advantage in operating a UASB reactor at a sludge age exceeding 100 days, because removal of biodegradable matter was essentially complete (Santos et al. 2016). From the data in Table 5 it is noted that the sludge age determined on the basis of total sludge mass and expulsion rate is 10 to 30% longer than that based on the volatile sludge. This indicates that there was selective retention of inorganic solids in the reactors, especially when the sludge age was long. For shorter sludge ages the difference between the total and volatile sludge ages was much smaller.

The total and volatile solids concentrations in the eight reactors at each sample point (Table 3) are shown in Figure 6(a)6(c) as a function of digestion zone height, by reactor, for upflow rates of 0.5 and 1 m3.d−1. The reactor configurations, and the applied hydraulic and organic loads, affected the sludge masses.
Figure 6

(a) Sludge and total volatile solids concentrations versus depth in digestion zone for UASB – Y and conventional reactor operated with HRT 12 hours. (b) Sludge and total volatile solids concentrations versus depth in digestion zone for reactors T1, T2 and T3 operated at 12 and 6 h. (c) Sludge and total volatile solids concentrations versus depth in digestion zone for reactors S1, S2 and S3 operated at 12 and 6 h.

Figure 6

(a) Sludge and total volatile solids concentrations versus depth in digestion zone for UASB – Y and conventional reactor operated with HRT 12 hours. (b) Sludge and total volatile solids concentrations versus depth in digestion zone for reactors T1, T2 and T3 operated at 12 and 6 h. (c) Sludge and total volatile solids concentrations versus depth in digestion zone for reactors S1, S2 and S3 operated at 12 and 6 h.

Close modal

Figure 6(a)6(c) show the concentration profiles in the reactors. Considerable concentration stratification is evident. This is attributed to changes in sludge settleability with depth, since the upflow velocity of the liquid phase was constant and the mixing energy of the rising bubbles – from biogas generation – (<1 W.m−3) was insufficient to induce significant turbulence.

The experimental data underlying Figure 6(a)6(c) also indicate qualitative stratification, with a greater proportion of inorganic matter in the lower regions of the reactors. The volatile fraction also decreased with the cross-sectional area of the digestion zone, especially when the upflow rate was low. It is possible that the higher inorganic fraction improved sludge settleability. This is not unexpected as selective expulsion of organic sludge could mean that inorganic sludge, which tends to have better settling characteristics, was retained more efficiently and accumulated in the reactor's lower regions. This may have been enhanced by slow sintering, with the inorganic sludge gradually acquiring greater mechanical rigidity. This was observed but not quantified during the study. Selective inorganic sludge retention is also reflected in the calculated sludge age (the ratio of sludge mass retained to daily mass expelled), as shown in Table 5, the total sludge age calculated is systematically higher than that calculated on the bases of volatile solids.

Figure 6(a)6(c) also show that the sludge layer in most reactors did not extend upward through the digestion zone. In other words, it was not possible to build the sludge up over the full height of this zone, even if hydrodynamic conditions were adequate for sludge settling and accumulation. It is believed that sludge retention by gravity settling was counter-balanced by break-up of the sludge flocs, with parts of the flocs with poor settling characteristics being swept out of the reactor. The sludge particles leaving the digestion zone may be discharged with the effluent, but it is also possible that these latter will flocculate in the settling zone and form bigger particles with better settling characteristics that can then return to the digestion zone.

Sludge activity

The biological and mechanical sludge properties were also investigated. The key biological and mechanical properties, respectively, were the SMPP, and the settleability as defined by the Vesilind constants. These are important because they indicate the treatment capacity limits of the reactors. The settleability constants enable estimation of the maximum theoretical sludge mass that could be retained. On the other hand, using the sludge mass in the reactors and the SMPP, the maximum theoretical extent of methanogenesis can be calculated. These theoretical values can then be compared to the actual experimental observations.

Table 6 presents the SMPPs for the 6 and 12 hour HRTs evaluated. Tests were carried out for both HRTs using samples withdrawn from the reactors as well as expelled sludge samples. The data indicate that the SMPP was approximately the same in all reactors, with no significant differences between the collected and expelled sludges. Reactor load seems to affect SMPP value: longer retention times and lower loads tend to have larger SMPPs. This is as expected, the higher the load, the greater the volatile fraction – the latter comprising flocculated particulate influent organic matter, reducing the fraction of methanogens in the volatile sludge (Van Haandel and Lettinga 1994).

Table 6

Average values of SMPP, expressed as g-COD.gXv.d−1 for sludge from the reactors

 T1T2T3S1S2S3YC
Sewage flow = 500 L.d1 
 Average in reactor 0.14 0.14 0.17 0.16 0.13 0.13 0.10 0.11 
 Expelled from reactor 0.15 0.13 0.13 0.10 0.09 0.10 0.10 0.11 
Sewage flow = 1,000 L.d1 
 Average in reactor 0.17 0.14 0.18 0.15 0.15 0.13 – – 
 Expelled from reactor 0.13 0.12 0.17 0.11 0.10 0.13 – – 
 T1T2T3S1S2S3YC
Sewage flow = 500 L.d1 
 Average in reactor 0.14 0.14 0.17 0.16 0.13 0.13 0.10 0.11 
 Expelled from reactor 0.15 0.13 0.13 0.10 0.09 0.10 0.10 0.11 
Sewage flow = 1,000 L.d1 
 Average in reactor 0.17 0.14 0.18 0.15 0.15 0.13 – – 
 Expelled from reactor 0.13 0.12 0.17 0.11 0.10 0.13 – – 

Sludge settleability

The settleability tests confirmed that Vesilind's formula (Equation (3)) describes settling perfectly in all cases. Table 7 shows the test results for sludge within and expelled from the reactors. The sludge characteristics indicated that the sludge in all reactors had excellent settleability, all with very low compressibility constant, k. There were no great differences between the constants in different reactors, but there was a clear difference between the sludge in the reactors and that expelled from them. In particular, k, for the expelled sludge is significantly higher on average, indicating that expelled sludge settleability was inferior to that of the sludge in the reactors.

Table 7

Vesilind's settleability constants as functions of sampling point levels

ReactorT1T2T3S1S2S3YC
Compressibility constant k (l.g1) for sewage flow = 500 L.d1 
 Average in reactor 0.10 0.08 0.07 0.14 0.07 0.08 0.10 0.09 
 Expelled from reactor 0.14 0.10 0.13 0.18 0.09 0.11 0.16 0.12 
Compressibility constant k(l.g1) for sewage flow = 1,000 L.d1 
 Average in reactor 0.11 0.08 0.08 0.09 0.08 0.10 – – 
 Expelled from reactor 0.12 0.11 0.10 0.12 0.11 0.12 – – 
Unhindered settling velocity v0(m.h1) for sewage flow = 500 L.d1 
 Average in reactor 10 15 10 17 10 18 14 
 Expelled from reactor 18 10 16 17 12 15 15 13 
Unhindered settling velocity v0(m.h1)for sewage flow = 1,000 L/d1 
 Average in reactor 10 11 15 10 – – 
 Expelled from reactor 16 10 17 14 13 – – 
ReactorT1T2T3S1S2S3YC
Compressibility constant k (l.g1) for sewage flow = 500 L.d1 
 Average in reactor 0.10 0.08 0.07 0.14 0.07 0.08 0.10 0.09 
 Expelled from reactor 0.14 0.10 0.13 0.18 0.09 0.11 0.16 0.12 
Compressibility constant k(l.g1) for sewage flow = 1,000 L.d1 
 Average in reactor 0.11 0.08 0.08 0.09 0.08 0.10 – – 
 Expelled from reactor 0.12 0.11 0.10 0.12 0.11 0.12 – – 
Unhindered settling velocity v0(m.h1) for sewage flow = 500 L.d1 
 Average in reactor 10 15 10 17 10 18 14 
 Expelled from reactor 18 10 16 17 12 15 15 13 
Unhindered settling velocity v0(m.h1)for sewage flow = 1,000 L/d1 
 Average in reactor 10 11 15 10 – – 
 Expelled from reactor 16 10 17 14 13 – – 

Having determined the settleability constants, the theoretical sludge concentration that can be maintained in the reactor can be calculated, using Equation (3a), for any upward liquid velocity, vup, in a UASB:
formula
3a

If the expelled sludge's Vesilind constants are used in Equation (3a) it can be shown that the theoretical sludge concentration that could be retained greatly exceeds, in all cases, the concentrations achieved (see Figure 6). In other words, sludge particles with poor settling characteristics are expelled from the digestion zone because they cannot be retained. They tend to flocculate in the settling zone, however, and acquire better settleability, so that they can settle again and become part of the sludge in the digestion zone. It is possible, therefore, that, if the settled solids from a settler were returned to the UASB reactor, much of the material would be retained. If this were so, a larger sludge mass could be maintained and the reactor operated at longer sludge age, and thus at higher efficiency and with lower sludge production.

Table 8 shows the theoretical sludge concentrations that could be retained, calculated using Equation (3a), for the constant values for the different reactors and the average sludge concentration as determined experimentally. The theoretical sludge mass values tend to be larger than the experimental equivalents for most reactors. As can be seen in Figure 6, however, the actual sludge concentration in the bottoms of the reactors is much higher than the theoretical maximum calculated for expelled sludge. This is because some sludge had better settleability than the average value used to calculate the theoretical concentration. Equally, in the digestion zone upper sections, the sludge concentration is much lower than the theoretical maximum, confirming that mechanisms other than settling are relevant to sludge retention in the digestion zone, otherwise, the digestion zone would fill entirely with sludge at a concentration equal to the theoretical maximum. Thus, the low sludge concentrations in the upper sections of the reactors is attributed to floc break-up, with small particles carried upwards by the liquid flow. Some of these can flocculate in the transition or settling zones, and return to the sludge bed. Others will be discharged with the effluent.

Table 8

Comparison of the theoretical and experimental, treatment and sludge retention capacities for UASBs studied

VariablesT1T2T3S1S2S3YC
Sewage flow = 500 l.d1 
 Theoretical sludge concentration (g.l−135 49 50 33 57 52 41 43 
 Experimental sludge concentration (g.l−129 39 29 40 36 30 49 46 
 Theoretical treatment capacity (g-COD.d−1739 1,108 1,469 960 739 960 960 818 
 Experimental treatment capacity (g-COD.d−190 151 167 255 254 256 137 129 
Sewage flow = 1,000 l.d1 
 Theoretical sludge concentration (g.l−126 37 40 34 42 34 – – 
 Experimental sludge concentration (g.l−111 18 25 39 37 28 – – 
 Theoretical treatment capacity (g-COD.d−1367 504 1,425 1,548 1,728 1,279 – – 
 Experimental treatment capacity (g-COD.d−1294 272 231 200 278 233 – – 
VariablesT1T2T3S1S2S3YC
Sewage flow = 500 l.d1 
 Theoretical sludge concentration (g.l−135 49 50 33 57 52 41 43 
 Experimental sludge concentration (g.l−129 39 29 40 36 30 49 46 
 Theoretical treatment capacity (g-COD.d−1739 1,108 1,469 960 739 960 960 818 
 Experimental treatment capacity (g-COD.d−190 151 167 255 254 256 137 129 
Sewage flow = 1,000 l.d1 
 Theoretical sludge concentration (g.l−126 37 40 34 42 34 – – 
 Experimental sludge concentration (g.l−111 18 25 39 37 28 – – 
 Theoretical treatment capacity (g-COD.d−1367 504 1,425 1,548 1,728 1,279 – – 
 Experimental treatment capacity (g-COD.d−1294 272 231 200 278 233 – – 

It is also possible to compare the theoretical and actual methane production capacities. The theoretical production level is the product of the volatile sludge mass and SMPP for each reactor. The true production can be calculated as the product of the digested COD fraction and the applied load. Both theoretical and actual are shown in Table 8. In all cases the theoretical treatment capacity or methane production is much greater than that achieved, i.e., methane production was not the limiting process for anaerobic treatment. This is confirmed because the VFA concentration (including acetate, the main methanogenic substrate) was always low.

  • (1) A simple, 250 L, UASB-based reactor can treat sewage from a single family efficiently. Its organic matter removal efficiency is much greater than can be achieved in a traditional ST, or combination ST and anaerobic filter.

  • (2) A set of 8 × 250 L reactors with differing geometries could all digest raw sewage efficiently at 12 and 6 h HRTs. The experiments showed clearly that, within the size range investigated, geometry had no significant influence on performance. Thus, configuration preferences can be based on considerations of ease of manipulation and installation, and a shallow unit (depth 1.5 m) with a relatively large footprint (diameter 0.7 m) can be chosen.

  • (3) Sludge retention in the settling zone, above the digestion zone, contributes to sludge mass maintenance, but settling is not the sludge retention determining mechanism in UASBs.

  • (4) Sludge age was confirmed as the fundamental operating parameter determining UASB performance and, particularly, division of the influent COD into three fractions: (i) discharged in the effluent, (ii) converted to sludge, and (iii) converted to methane. If the sludge age is the same, UASBs treating the same sewage at different HRTs will tend to have the same effluent quality and sludge production rate, and therefore the same digestion efficiency. The lower the sludge age, the higher the proportions of influent COD going to the effluent or the excess sludge.

  • (5) The SMPP of all sludge investigated in the study was within the range commonly found for sludges from treating raw sewage, i.e. 0.10 to 0.15 g-COD.g−1Xv.d−1. The SMPP of the sludge in the reactors did not differ significantly from that of the expelled sludge. In all cases, it was much greater than the actual specific methane production, showing that methanogenesis was not the efficiency limiting process for anaerobic digestion. In fact the incomplete removal of biodegradable organic material arises mainly because the preceding hydrolytic, acidogenic and acetogenic processes could not develop fully.

  • (6) The settleability of the sludge within the reactors was always superior to that of the sludge expelled from them, although sludge retention was not the sludge mass retention determining mechanism. In most reactors the digestion zone was not completely filled with sludge, showing that a break-up mechanism exists, producing small particles that settle poorly and tend to escape from the reactor. Some of these tend to flocculate in the settling zone and may return to the digestion zone as recovered sludge, others are discharged with the effluent.

The experimental investigation was carried out with financial support from agencies of the Brazilian Government (CNPq, ANA and FINEP).

Aiyuk
S.
Forrez
I.
Lieven
K.
Van Haandel
A.
Verstraete
W.
2006
Anaerobic and complementary treatment of domestic sewage in regions with hot climates, a review
.
Bioresource Technology
97
(
17
),
2225
2241
.
Coelho
A. L. S. S.
Nascimento
M. B. H.
Cavalcanti
P. F. F.
Haandel
A. C.
2003
The UASB reactor as an alternative for the septic tank for on-site sewage treatment
.
Water Science Technology
48
(
12
),
659
664
.
Jorsaraei
A.
Gougol
M.
Van Lier
J. B.
2014
A cost effective method for decentralized sewage treatment
.
Process Safety and Environmental Protection
92
,
815
821
.
Ladu
J. L. C.
X.
2014
Effects of hydraulic retention time, temperature, and effluent recycling on efficiency of anaerobic filter in treating rural domestic wastewater
.
Water Science and Engineering
7
(
2
),
168
182
.
Lew
B.
Lusting
I.
Beliavski
M.
Tarre
S.
Green
M.
2011
An integrated UASB-sludge digester system for raw domestic wastewater treatment in temperate climates
.
Bioresource Technology
102
(
7
),
4921
4924
.
Metcalf and Eddy
2014
Wastewater Engineering: Treatment and Resource Recovery
.
McGraw-Hill
,
New York
,
USA
.
Standard Methods for the Examination of Water and Wastewater
2012
22nd edn,
American Public Health Association/American Water Works Association/Water Environment Federation
Washington DC, USA
.
United States Environmental Protection Agency – USEPA
2005
Handbook for Managing Onsite and Clustered (Decentralized) Wastewater Treatment Systems
.
EPA/832-B-05-001
.
EPA: Environmental Protection Agency
,
Washington, DC
.
Van Haandel
A. C.
Lettinga
G.
1994
Anaerobic Sewage Treatment: A Practical Guide for Regions with a Hot Climate
.
John Wiley and Sons
,
Chichester
,
UK
.
Vesilind
P. A.
1967
Theoretical considerations: design of prototype thickeners from batch settling tests
.
Water and Sewage Work
115
(
7
),
302
307
.