Theoretical Analysis of the Velocity Gradient
When an incompressible fluid flows through a cylindrical tube its velocity relative to the walls changes as a function of the tube radius. In general, this velocity distribution is parabolic: the greatest velocities are achieved at the center of the tube (where R=0) eventually tapering off to 0 at the walls. The parabolic nature of the distribution arises from cylindrical symmetry as well as the fact that the fluid does not move at the walls (the "no-slip" condition).
This gradient in the velocity profile contributes to the force that a floc experiencing roll-up feels. Flocs actually begin to roll up when the velocity at their edge exposed to the flow exceeds some critical value, which is highly dependent on the floc's diameter, its density, and the capture velocity of the system, among other things.
Unable to render embedded object: File (Floc Roll up picture.png) not found.
Applying a force balance to this (where FB = force of buoyancy, FD = force of drag, FG = force of gravity, and theta = 60 degrees):
Unknown macro: {latex}
\large
$$
\sum
Unknown macro: {F = 0}
$$
Unknown macro: {latex}
\large
$$
F_b + F_d - mg\sin (\theta ) = 0
$$
Unknown macro: {latex}
\large
$$
- mg\sin \theta + \rho w gV
Unknown macro: {proj}{{V^2 } \over 2} = 0
\sin \theta + \rho w C_D A
$$
Solving this equation for the velocity:
Unknown macro: {latex}
\large
$$
{{V^2 } \over 2}\rho w C_D A
Unknown macro: {proj}
= \rho p V
g\sin \theta - \rho w gV
Unknown macro: {proj}
\sin \theta
$$
Unknown macro: {latex}
\large
$$
V^2 = {{2g\sin \theta V_
Unknown macro: {proj}
} \over {\rho w C_D A
}}(\rho _p - \rho _w )
$$
In this approximation a sphericity factor of 45/24 is used to take into account the fact that flocs are not spherical.
Unknown macro: {latex}
\large
$$
{{V_
Unknown macro: {proj}
} \over {A_
}} = d_p
$$
This equation and all subsequent equations are incorrect. Check the equations for volume and projected area of a sphere. Change the symbol for the volume of the particle throughout the document. One commonly used symbol is a V with a line through it to distinguish it from velocity.
Assuming that the flow through the tube settlers is creeping flow (with a Reynolds Number much less than one), we can approximate the drag coefficient as 24/Re; when sphericity is taken into account, this becomes 45/Re, giving:
Unknown macro: {latex}
\large
$$
V = \sqrt {{2 \over {45}}g\sin (\theta )d_p {\mathop
Unknown macro: {rm Re}
\nolimits} {{(\rho _p - \rho _w )} \over
Unknown macro: {rho _w }
}}
$$
Replace Re with its definition and then solve the expression for velocity.
This is the velocity that a particle would roll down the tube if there were no flow up the tube.
At this velocity a floc that is sitting on the bottom wall of the tube settler should not move. We should call this something like the critical or failure velocity, since velocity values greater than this quantity (at the edge of the floc exposed to the flow) will cause roll-up.
VELOCITY GRADIENT ANALYSIS
Unknown macro: {latex}
\large
Navier-Stokes for a cylindrical system
$$
Unknown macro: {1 over r}
{\partial \over {\partial r}}(r{{\partial v_z } \over {\partial r}}) =
Unknown macro: {1 over mu }
(\partial P} \over {\partial z)r
$$
Unknown macro: {latex}
\large
Integrate to obtain:
$$
r{{\partial v_z } \over {\partial r}} = {1 \over {2\mu }}(\partial P} \over {\partial z)r^2 + c_1
$$
This constant of integration should be zero, since no velocity gradient exists at the point of maximum velocity, achieved at the center of the tube where R=0.
Unknown macro: {latex}
\large
$$
\partial v_z } \over {\partial r = {1 \over {2\mu }}(\partial P} \over {\partial z)r
$$
Unknown macro: {latex}
\large
$$
\partial P} \over {\partial z= {{\Delta P} \over L}
$$
Unknown macro: {latex}
\large
$$
V_
Unknown macro: {av}
= R^2 \Delta P} \over {8\mu L = {Q \over {\pi R^2 }}
$$
Unknown macro: {latex}
\large
$$
{{\Delta P} \over L} = {{Q8\mu } \over {\pi R^4 }}
$$
Unknown macro: {latex}
\large
$$
\partial v_z } \over {\partial r = {1 \over {2\mu }}{{Q8\mu } \over {\pi R^4 }}r
$$
Unknown macro: {latex}
\large
$$
\partial v_z } \over {\partial r = - {{4Q} \over {\pi R^4 }}r
$$
Unknown macro: {latex}
\large
$$
V_z (r) = {{2Q} \over {\pi R^2 }}[1 - (
Unknown macro: {r over R}
)^2 ] = {{2Q} \over {A_
Unknown macro: {tube}
}}[1 - (
)^2 ]
$$
Unknown macro: {latex}
\large
Finally:,
$$
V_z (r) = V_
Unknown macro: {max }
[1 - (
Unknown macro: {r over R}
)^2 ]
$$
Unknown macro: {latex}
\large
We can now evaluate this expression to determine the velocity at the failure point, which for this model is equal to (tube radius - floc diameter) :
$$
V_z ({\rm{R}}_{{\rm
Unknown macro: {tube}
}} {\rm{ - d}}{\rm{p}} ) = V
Unknown macro: {max }
[(1 - {\rm{R_{{\rm
}} {\rm{ - d}}_{\rm{p}} } \over R})^2 ]
$$
Unknown macro: {latex}
\large
$$
V_z (P_
Unknown macro: {failure}
) = V_
Unknown macro: {max }
[(1 - {\rm{R_{{\rm
Unknown macro: {tube}
}} - d_p } \over R})^2 ] = V_
Unknown macro: {critical}
$$
Unknown macro: {latex}
\large
$$
V_
Unknown macro: {critical}
= \sqrt {{2 \over {45}}g\sin (\theta )d_p {\mathop
Unknown macro: {rm Re}
\nolimits} {{(\rho _p - \rho _w )} \over
Unknown macro: {rho _w }
}}
$$
Unknown macro: {latex}
\large
When the velocity at the critical point exceeds the critical velocity, roll-up occurs. This suggests an inequality of the form:
$$
V_
Unknown macro: {max }
[1 - ({\rm{R_{{\rm
Unknown macro: {tube}
}} - d_p } \over R})^2 ] > \sqrt {{2 \over {45}}g\sin (\theta )d_p {\mathop
Unknown macro: {rm Re}
\nolimits} {{(\rho _p - \rho _w )} \over
Unknown macro: {rho _w }
}}
$$
Unknown macro: {latex}
\large
Solving this inequality for the diameter of the particle allows us to generate a function that will help us to analyze roll-up in terms of a critical floc diameter.
Fix the equations above. Add the floc density from the fractal floc model. Then assuming you hold the capture velocity constant, it should be possible to find the diameter of the tube that causes roll up for a particular floc diameter.
This equation actually happens to be cubic in the diameter of the floc and has the following form:
$$
{\rm{F(}}d_p ) = c_1 d_p ^3 - c_2 d_p ^2 + c_3 d_p - \omega
$$
Unknown macro: {latex}
\large
$$
{\rm{F(}}d_p ) = ({1 \over {R^4 }})d_p ^3 - ({2 \over {R^3 }})d_p ^2 + ({1 \over {R^2 }})d_p - g\sin (\theta ){\mathop
Unknown macro: {rm Re}
\nolimits} {{(\rho p - \rho _w )} \over {\rho _w }}{1 \over {V
Unknown macro: {max }
^2 }}
$$
Unknown macro: {latex}
\large
The coefficients in this function are only functions of the radius of the tube, and an analysis of how they change with the radius could possibly give some physical insight into the phenomenon of floc roll-up as it relates to the geometry of the tube settlers.
$$
{\rm{F(}}d_p ) < 0
$$
No roll-up.
Unknown macro: {latex}
\large
$$
{\rm{F(}}d_p ) = 0
$$
This is the critical diameter at which a floc will not move.
Unknown macro: {latex}
\large
$$
{\rm{F(}}d_p ) > 0
$$
Roll-up occurs.
During Spring 2009 the Plate Settler Spacing team first identified the roll-up phenomenon, and their subsequent analysis is centered around a critical velocity gradient occurring at the walls of the tube settler. While this analysis provides reasonable data, the fundamental force balance has steered the Summer 2009 team in a direction that is more focused on the physical properties of the floc as they relate to roll-up. This is beneficial because a model based entirely on the velocity gradient may lend itself to difficulty in calculation, since this gradient could very well depend on the influent water turbidity and chemistry. The Spring 2009 team's research can be found [here].
It might be useful to analyze the final cubic function's behavior around its inflection points - these may indicate critical values or "states" the system moves through as the capture velocity changes.