Theoretical Analysis of the Velocity Gradient
This analysis will be done for both a tube and plate settler. Although tube settlers are used in the lab, plate settlers are used in the Honduras plants. Therefore, a model that takes into account the differences between the two apparatuses must be developed.
Calculation of Ratio of Settling Velocity to Particle Velocity
To begin the calculations for the critical velocity, it is important to note that the maximum velocity at the center of a pipe is two times the average velocity and the maximum velocity between two parallel plates is 1.5 times the average velocity.
In order to determine the critical velocity at which floc particles will begin to roll up the tube and into the effluent, we compare the settling velocity with the particle velocity experienced from the velocity gradient.
The settling velocity of a particle in a tube settler can be expressed as follows:
Unknown macro: {latex}
\large
$$
V_t = {{gd_0 ^{\left( {3 - D_
Unknown macro: {Fractal}
} \right)} d^{\left( {D_
- 1} \right)} } \over {18\Phi \nu }}\left( {{{\rho _
Unknown macro: {Floc_0 }
} \over {\rho _
Unknown macro: {H_2 O}
}} - 1} \right)
$$
Where:
g = Gravity
Unknown macro: {latex}
\large $$d_o $$
= size of the primary particles
Unknown macro: {latex}
\large $$D_
Unknown macro: {fractal}
$$
= fractal dimension of the floc particles
Unknown macro: {latex}
\large $$\Phi $$
= shape factor for drag on flocs which is equal to
Unknown macro: {latex}
\large $$\nu $$
= viscosity
Unknown macro: {latex}
\large $$\rho _
Unknown macro: {floc}
$$
= density of the floc particle
Unknown macro: {latex}
\large $$\rho _
Unknown macro: {H_2 O}
$$
= density of water
The particle velocity expereinced as a result of the velocity gradient can be expressed as follows:
Unknown macro: {latex}
\large
$$
V_
Unknown macro: {particle}
= V_
Unknown macro: {ratio}
V_\alpha \left[ {1 - \left( {{{{{d_
Unknown macro: {tube}
} \over 2} - d_
Unknown macro: {Floc}
} \over {{{d_
Unknown macro: {Tube}
} \over 2}}}} \right)^2 } \right]
$$
Where
Unknown macro: {latex}
\large $$V_\alpha $$
= directional velocity in the tube settler
Unknown macro: {latex}
\large $$d_
Unknown macro: {tube}
$$
= diameter of the tube settler
Unknown macro: {latex}
\large $$ d_
Unknown macro: {floc}
$$
= the diameter of floc particles
Therefore, the ratio between the settling velocity of the particle and the velocity experienced as a result of the velocity gradient can be expressed as (Explain how the angle figured into the terminal velocity equation.)
Unknown macro: {latex}
\large
$$
\Pi _V = {{{{g\sin (\alpha )d_0 ^2 } \over {18\Phi \nu }}{{\rho _
Unknown macro: {Floc_0 }
- \rho _
Unknown macro: {H_2 O}
} \over {\rho _
}}\left( {{{d_
Unknown macro: {Floc}
} \over
Unknown macro: {d_0 }
}} \right)^{D_
Unknown macro: {Fractal}
- 1} } \over {V_
Unknown macro: {ratio}
{{V_
Unknown macro: {Up}
} \over {\sin (\alpha )}}\left[ {1 - \left( {{{{{d_
Unknown macro: {Tube}
} \over 2} - d_
} \over {{{d_
Unknown macro: {Tube}
} \over 2}}}} \right)^2 } \right]}}
$$
This ratio is a function of particle diameter. When this ratio is greater than one (ie the settling velocity is greater than the velocity experienced by the floc particles in the tube), the flocs will fall back into floc blanket and fail to travel to the effluent. When this ratio is equal to one, the particles will remain stationary in the tube settler. And when the ratio is less than one, the velocity of the particles will exceed the settling velocity and the floc particles will roll up into the effluent, creating a highter turbidity.
Figure 1 shows what size floc could be captured by different tube settler diameters. The lines cross the y value of 1 when the sedimentation velocity matches the upflow velocity at the floc diameter, and this is the critical velocity, i.e. the velocity in which floc roll up will begin, for each of the given tube diameters and plate settler spacing.
Figure 1: The ratio of Sedimentation Velocity to Fluid Velocity vs. Floc Diameter
The graph is cut off at a particle size of 1x10^-3 because there are no floc particles that will be smaller than this. Therefore, for plate spacings of 2 and 5 cm and for tube diameters of 5 cm, there will be no floc roll up under any conditions.
Calculation of the Minimum Diameter of the Flocs that Settle from the Sedimentation Velocity Equation
Assuming an upward flow velocity of 100 m/day, which is standard in all AguaClara plants, the diameter of floc that will roll-up was determined by using a root finding algorithm, and the plate settler spacing or tube diameter was plotted versus the minimum floc diameter.
Figure 2: Plate Spacing or Tube Diameter vs. Minimum Floc Diameter
The minimum floc diameter corresponds to the minimum size of particles that will still settle out of the tube and return to the floc blanket instead of going into the effluent. This graph is based on the capture velocity theory, and must be revised to take into consideration floc roll-up that will occur in the plate settlers.
To ensure that flocs in the plate settler don't roll up, we can calculate the minimum floc diameter that settles from the sedimentation velocity equation. Solving for the floc diameter,
Unknown macro: {latex}
\large
$$
d = d_0 \left( {{{18V_t \Phi \nu _
Unknown macro: {H_2 O}
} \over {gd_0^2 }}{{\rho _
} \over {\rho _
Unknown macro: {Floc_0 }
- \rho _
Unknown macro: {H_2 O}
}}} \right)^{{1 \over {D_
Unknown macro: {Fractal}
- 1}}}
$$
This equation is based on a linearization of the velocity profile at the wall. It uses the velocity gradient at the wall to determine the velocity as a function of distance from the wall. For small floc sizes this linearization is valid and produces an analytical solution. This linearized approach was used because the final results are very similar, and this linearized approach makes it much easier to solve for the critical velocity.
Figure 3 shows how the linearized equations provide only tiny divergence for big flocs in small tubes.
Figure 3: Floc Diameter vs. Spacing
Figure 4 graphically displays the linear velocity gradient solutions. The curves are not quite straight on a log log plot. This is due to the quadratic in the velocity profile. We can obtain a very good approximation by using the velocity gradient at the wall and assuming a linear velocity gradient. That assumption makes an analytical solution possible.
Figure 4: Floc Spacing vs. Floc Diameter
The equation below represents the spacing of the plate settlers to capture a floc particle of a particular size. This is based on the critical velocity model.
Unknown macro: {latex}
\large
$$
S = V_
Unknown macro: {up}
{{108\Phi \nu _
Unknown macro: {H_2 O}
d^2 } \over {g\sin ^2 (\alpha )d_0^3 }}\left( {{
Unknown macro: {d_0 }
\over d}} \right)^{D_
Unknown macro: {Fractal}
} {{\rho _
} \over {\rho _
Unknown macro: {Floc_0 }
- \rho _
Unknown macro: {H_2 O}
}}
$$
Where:
Unknown macro: {latex}
\large$$ alpha $$
= The angle of the tube settler (60 degrees)
Figure 5 represents the minimum spacing that will capture a floc particle with a particular settling velocity. The settling velocity can be translated into a floc diameter by the terminal settling velocity equation.
Figure 5: Minimum spacing vs. Floc Sedimentation Velocity
Figure 6 represents the absolute minimum plate settler spacing that will capture floc particles with a settling velocity of 10 m/day in an AguaClara plant. Theoretically, any spacing below the intersection of two lines would produce a worsened effluent turbidity because larger particles are no longer being captured by the plate settler.
Figure 6: Minimum Plate Settler Spacing vs. Capture Velocity
In order to complete the calculations, the Reynold's Number was checked to ensure that the flow in the settler tube is in fact laminar. Next, the entrance region was checked to ensure that the parabolic velocity profile was fully established.
Figure 7 shows that the entrance regions are always less than the length of the plate setltlers. Thus, based on this analysis, the best designs would include a plate settler with a larger spacing. This maximizes the ratio between the settling velocity and the velocity experienced by the particle, and would ideally prevent the maximum amount of particles from escaping into the effluent.
Figure 7: Ratio of Entrance Length to Plate Length vs. Plate Spacing
*More detail on the calculation process outlined above can be found in the Math CAD File.