Flocculation Design

Important Note: This chapter introduces uncertainty and empirical design. Some of the parameters used to design AguaClara flocculators are based on what has been shown to work in the field, as opposed to having been derived scientifically. To make sure that the reader is aware of these concepts and parameters that don’t yet have a thorough basis in research, they will be highlighted in red when they appear.

AguaClara’s approach to flocculator design is the same as it is for any other unit process. First, critical design criteria, called inputs, are established. These criteria represent the priorities that the rest of the design will be based around. Once these parameters are established, then the other parameters of the design, which are dependent on the inputs, are calculated based on certain constraints.

Take the CDC as an example of this design process in AguaClara Flow Control and Measurement Technologies; its inputs are \(h_{L_{max}}\), \(\sum K\), \(\Pi_{Error}\), and the discrete dosing tube diameters \(D\) that are available at hardware stores or pipe suppliers. Its dependent variables include the number and length of the dosing tubes and the flow through the CDC system.

The flocculator is more complex to design than the CDC, as it has more details and parameters and the equations for those details and parameters are very interdependent. Therefore, there are many ways to design an AguaClara flocculator, and many different sets of critical design criteria to begin with. Below are design approaches for several of the different potential flocculator geometries.

In all cases we will need equations relating head loss, velocity gradient, viscosity, and collision potential. Flocculators must operate effectively over a range of temperatures and the critical design temperature is the coldest operating temperature. The viscosity shown in subsequent equations is the viscosity at the coldest operating temperature.

The following equations describe the design process assuming that a velocity gradient is specified for the coldest operating temperature. Given \(\tilde{G}\), we can very easily find \(\theta\):

(209)\[\theta = \frac{\tilde{G} \theta}{\tilde{G}}\]

Finally, we take retention time \(\theta\) over plant flow rate \(Q\) to get the required volume of the flocculator:

\[\forall_{floc} = \theta Q\]

Now that we have the basic parameters defined, we can start to design the details of the flocculator, starting from the physical dimensions.

Horizontal-Horizontal flow flocculators

Horizontal-Horizontal flocculators can be easily designed to have an optimal \(\Pi_{H_eS}\). Our current understanding is that the optimal value to achieve uniform velocity gradients and to allow the flow to fully expand is when \(\Pi_{H_eS} = 6\). The flow geometry can be optimized to reduce construction costs by varying the relationship between \(S\) and \(W\). The dimension \(W\) is in the vertical direction and thus increasing W will increase the water depth and reduce the plan view area of the flocculator. If land area is constrained and the cost of roofing is considered it is very likely that the optimal solution will require \(W>S\).

Horizontal-Vertical flow flocculators

  1. Input parameters
    • Specify: - \(\tilde{G}\), velocity gradient at the coldest operating temperature - \(\tilde{G} \theta\), collision potential at the coldest operating temperature - \(Q\), plant flow rate - \(H\), height of water at the end of the flocculator - \(L_{max, \, sed}\), max length of a flocculator channel based on clarifier length - \(W_{min, \, human}\) minimum width of a single channel based on the width of the average human hip (someone’s got to go down there…)

    • Find: - \(\theta\), hydraulic retention time - \(\forall_{floc}\), flocculator volume - \(h_{L_{floc}}\), head loss (independent of operating temperature)

  2. Physical dimensions
    • Calculate: - \(L_{channel}\), actual channel length - \(n_{channels}\), amount of channels - \(W_{channel}\), actual channel width

  3. Hydraulic parameters
    • Calculate: - \(H_e\), distance between baffle/obstacle induced flow expansions - \(n_{obstacles}\), amount of obstacles per baffle space - \(S\), baffle spacing, distance between baffles

physical floccuation parameter

Fig. 129 Flocculator geometry definition including the effect of baffle thickness. Accounting for baffle thickness would be particularly important if ferrocement or wood were used for baffles.

We start by making sure that our flocculator will be able to flocculate effectively by defining \(h_{L_{floc}}\) and \(\tilde{G} \theta\). Fixing these two parameters initially allows us to easily find all other parameters which determine flocculator performance. Here are the current standards (effective 2023) in AguaClara flocculators:

  • \(\tilde{G} = 50 Hz\)

  • \(\tilde{G} \theta = 35,000\)

The plant flow rate \(Q\) is defined by the needs of the community that the plant is being desiged for. Additionally, the height of water at the end of the flocculator, \(H\), the maximum length of the flocculator based on the length of the clarifier length, \(L_{max, \, sed}\), and the minimum width of a flocculator channel required for a human to fit inside, \(W_{min, \, human}\), are also defined initially. Ordinarilly in AguaClara plants, the flocculator occupies the same length dimension as the clarifiers, which is why the length constraint exists. See Fig. 130 for a representation of how the flocculator and clarifiers are placed in a plant.

  • \(H = 2 \, {\rm m}\)

  • \(L_{max, \, sed} = 6 \, {\rm m}\)

  • \(W_{min, \, human} = 45 \, {\rm cm}\)

physical design criteria

Fig. 130 Layout of flocculator and clarifiers that was adopted starting with the 2nd AguaClara plant in Tamara, Honduras in 2008.

Physical Dimensions

Deriving the equations required to find the physical dimensions now and the hydraulic parameters (baffle/obstacle design) in the next section requires many steps. To simplify this design explanation the equation derivations are developed in Fluid Mechanics Derivations. All complex equations which seemingly came out of nowhere will be derived in the derivation sheet.

Length

Flocculator length, \(L_{channel}`\) must meet two constraints: it must be less than or equal to the length of the clarifiers, as the flocculator is adjacent to the clarifier. This constraint is \(L_{max, \, sed}\). Next, the flocculator must be short enough to make sure the target volume of the flocculator is met, while still allowing for a human to fit inside \(L_{max, \, \forall}\). The constraint that wins out is the one that results in the *smaller* length value.

\[L_{max, \, sed} = 6 \, {\rm m} L_{max, \, \forall} = \frac{\forall}{n_{min, \, channels} W_{min, \, human} H}\]
Such that:
\(n_{min, \, channels} = 2\)

The reason why \(W_{min, \, human}\) is used is because it represents the absolute minimum of flocculator channel width. If the width ends up being larger, the length will decrease. \(n_{min, \, channels} = 2\) to make sure that the flow ends up on the correct side of the clarifier, as the image below shows. Note that there can only be an even number of flocculator channels, as explained in the image’s caption.

The equation for actual flocculator length is therefore:

\[L_{channel} = {\rm min}(L_{max, \, sed}, \, L_{max, \, \forall})\]
floc channels

Fig. 131 There are an even amount of flocculator channels to keep the AguaClara plant layout consistent for flows greater than 12 L/s. This ensures that the entrance tank, filter box, and filters can be kept in the same places across plants.

Width and Number of Channels

The width of a single flocculator channel must meet the following conditions: - Maintain \(\tilde{G}\) at the value found in the inputs section - Allow for \(3 < \frac{H_e}{S} < 6\). Recall that \(\frac{H_e}{S} = \Pi_{H_eS}\) - Allow for a human to be able to fit into a flocculator channel

The first two conditions are wrapped up into the following equation, which is derived in.

(210)\[W_{Min_{H_eS}} = \frac{\Pi_{H_eS_{min}}Q}{H_e}\left( \frac{K}{2 H_e \nu \tilde{G}^2} \right)^\frac{1}{3}\]

Given that the baffle minor loss coefficient is now known to be a function of \(\Pi_{H_eS}\) we substitute (196).

(211)\[W_{Min_{H_eS}} = \frac{\Pi_{H_eS_{min}}Q}{H_e}\left( \frac{ \left(\frac{\left(1 - \Pi_{vc}^{baffle}\right) ^ 2}{ \Pi_{vc}^{baffle} \Pi_{PlaneJet_{exp}} \Pi_{H_eS_{min}}}\right) ^ 2}{2 H_e \nu \tilde{G}^2} \right)^\frac{1}{3}\]

Group the parameters so that the dimensions inside the exponents are a simple as possible.

(212)\[W_{Min_{H_eS}} = \frac{Q}{\left(\nu \tilde{G}^2 H_e^4\right)^\frac{1}{3}}\left( \frac{ \left(1 - \Pi_{vc}^{baffle}\right) ^ 4 \Pi_{H_eS_{min}}}{2 \left( \Pi_{vc}^{baffle} \Pi_{PlaneJet_{exp}} \right)^2} \right)^\frac{1}{3}\]

This equation represents the absolute smallest width of a flocculator channel if we consider the lowest value of \(\Pi_{H_eS}\) and the highest possible value of \(H_e\):

\(H_e = H_{e_{max}} = H = 2 \, {\rm m}\), this implies that there are no obstacles between baffles \(\Pi_{H_eS} = \Pi_{ {HS}_{min} } = 3\)

Recall our other width constraint, \(W_{min, \, human} = 45 \, {\rm cm}\), which is based on our desire to have a human be able to fit into the channels. The governing constraint is the larger value of \(W_{min}\):

\[W_{min} = {\rm max}(W_{min, \, \Pi_{H_eS}}, \, W_{min, \, human})\]

We can find the number of channels, \(n_{channels}\) and their actual width in one last step, by finding the total flocculator width if there were no channels and dividing that by the minimum flocculator width, \(W_{min}\), found above. The equation for total flocculator width is based on our target volume:

\[W_{total} = \frac{\forall}{H L_{channel}}\]

Finally:

\[n_{channels} = \frac{W_{total}}{W_{min}}\]
Such that:
\(n_{channels}\) is an even number and is not 0. Usually, \(n_{channels}\) is either 2 or 4.

Now that we know \(n_{channels}\), we can find the actual width of a channel, \(W_{channel}\).

\[W_{channel} = \frac{W_{total}}{n_{channels}}\]

Hydraulic Parameters

Now that the physical dimensions of the flocculator have been defined, the baffle module needs to be designed. The parameter on which most others are based is the distance between flow expansions, \(H_e\). Recall that \(H_e = H\) when there are no obstacles in between baffles.

Height Between Expansions

We have a range of possible \(H_e\) values based on our window of \(3 < \frac{H_e}{S} < 6\). However, we have a limitation and a preference which shape how we design \(H_e\). Our limitation is that there can only be an integer number of obstacles. Our preference is to have as few obstacles as possible to make the baffle module as easy to fabricate as possible. Therefore, we want \(\frac{H_e}{S}\) to be closer to \(6\) than it is to \(3\); we are looking for \(H_{e_{max}}\).

We calculate \(H_{e_{max}}\) based on the physical flocculator dimensions. The equation for \(H_e\) is obtained by rearranging Equation (212). Because we have already designed the channel width, we substitute \(W_{channel}\) for \(W_{Min_{H_eS}}\). Since we are looking for \(H_{e_{max}}\), we also substitute \(\Pi_{{H_eS}_{max}}\) for \(\Pi_{H_eS}\). The result is:

(213)\[H_{e_{max}}^\frac{4}{3}= \frac{Q}{W_{Min_{H_eS}} \left(\nu \tilde{G}^2 \right)^\frac{1}{3}}\left( \frac{ \left(1 - \Pi_{vc}^{baffle}\right) ^ 4 \Pi_{{H_eS}_{max}}}{2 \left( \Pi_{vc}^{baffle} \Pi_{PlaneJet_{exp}} \right)^2} \right)^\frac{1}{3}\]

raise to 3/4

(214)\[H_{e_{max}}= \left[ \left(\frac{Q^3}{W_{Min_{H_eS}}^3 \nu \tilde{G}^2}\right)\frac{ \Pi_{{H_eS}_{max}}}{2 \left( \Pi_{vc}^{baffle} \Pi_{PlaneJet_{exp}} \right)^2} \right]^\frac{1}{4} \left(1 - \Pi_{vc}^{baffle}\right)\]

For \(\Pi_{H_eS}\) large enough such that the flow has fully expanded the following equation applies.

(215)\[H_{e_{max}} = \left[ \frac{K}{2 \nu \tilde{G}^2} \left( \frac{Q \Pi_{{HS}_{max}}}{W_{channel}} \right)^3 \right]^\frac{1}{4}\]

Note that this is the maximum distance between flow expansions, and does not account for the limitation that there must be an integer number of obstacles per baffle space. Thus, we need to find the actual distance between flow expansions. To do this, we determine and round up the number of expansions per baffle space using the ceiling function:

\[n_{expansions} = {\rm ceil}\left( \frac{H}{H_{e_{max}}} \right)\]

If we had used the floor() function instead, we would find that \(H_e\) would be larger than our upper bound, \(H_{e_{max}}\). From here, we can easily get to the actual number of flow expansions per baffle spacing:

\[H_e = \frac{H}{n_{expansions}}\]

Finally, we can obtain the number of obstacles per baffle space. The \(- 1\) in the equation is because the baffles themselves provide one flow expansion per baffle space.

\[n_{obstacles} = \frac{H}{H_e} - 1\]

Baffle Spacing

Finally, we can find the space between baffles, \(S\).

Checking the Flocculator Design

We then compare \(n_{spaces, \, required}\) to \(n_{spaces, \, actual}\) to make sure that they are equal.

Average Velocity

As water flows through the flocculators, the flocs will get larger and larger. As a result, their terminal velocity will increase. This is what we want. However, we need to make sure that the flocs don’t settle in the flocculator; that they instead all settle in the clarifier. To make sure of this, we need to make sure that the velocity of water in the flocculator is high enough to scour any flocs that fall to the bottom of the flocculator. The velocity required to scour flocs from the bottom and avoid floc accumulation is around \(v_{scour} = 15 \, {\rm \frac{cm}{s}}\). We can check our average velocity \(\bar v\) against this value. Unfortunately, the only method to increase the velocity through the flocculator is to increase the space between baffles, S, and decrease the channel width, W. Generally those changes aren’t practical and thus there isn’t any way to increase the velocity.

\[\bar v = \frac{Q}{W_{channel} S}\]

Head loss

Head loss is calculated for each expansion using the minor loss coefficient associated with that expansion in Equation (9). The total head loss is simple a sum of the head loss for all of the expansions.

Residence Time

It is now time to make our final check. We need to make sure that our actual residence time is at least as much as we designed for. Fortunately, in our design we did not account for the change in water level throughout the flocculator due to head loss. Therefore, the actual volume of water in the flocculator is actually greater than \(\forall_{floc}\). See Fig. 132 for clarification.

Flocculator_head_loss

Fig. 132 The water level in the flocculator decreases due to head loss. Flocculators may occupy multiple channels, but this extra triangle of water exists in any case.

Thus, the actual average water level in the flocculator is \(H + \frac{h_{L_{floc}}}{2}\). Thus, the actual residence time is:

\[\theta_{actual} = \frac{n_{channels} L_{channel} W_{channel} \left( H + \frac{h_{L_{floc}}}{2} \right)} {Q}\]

Check to see if \(\theta_{actual}\) is greater than \(\theta\).

[Haa98]

Johannes Haarhoff. Design of Around-the-End Hydraulic Flocculators. Journal of Water Supply: Research and Technology - Aqua, 47(3):142–152, June 1998. URL: https://doi.org/10.2166/aqua.1998.20 (visited on 2015-07-03).

[HvdW01]

Johannes Haarhoff and Jeremia J. van der Walt. Towards Optimal Design Parameters for Around-the-End Hydraulic Flocculators. Journal of Water Supply: Research and Technology - Aqua, 50(3):149–159, 2001.

[SO84]

Christopher R. Schulz and Daniel A. Okun. Surface water treatment for communities in developing countries. Wiley, New York, 1984. ISBN 978-0-471-80261-7.