Solid Liquid suspension (slurry) 

The use of liquid solid suspension  also called slurry  is common in process industries. Such suspension can for example be found at the very beginning of a process, after a liquid and a solid have been contacted in a reactor or mixer.
There can be different reasons to their use. It can offer an easier handling of the solid, for example, or be made necessary to obtain certain chemical reaction, or certain physical interaction, like the trapping of some aroma of the solids in the liquid.
Although found quite commonly in process industries, those mixtures are nevertheless posing different challenges, especially in terms of rheology where they have most of the time a behavior that is not newtonian,
This page aims at giving to the reader different tools to help to manage industrial situations with liquid solid slurries.
One of the 1st challenges when handling slurries is to avoid that the particles settle and that 2 phases reconstitute. It is advisable to assess the speed at which the slurry can sediment, it will give an idea of the stability of the suspension and of the efforts that will have to be considered during the design to prevent the settling.
The most direct method is carry out a test at pilot plant scale. However, such test is not always available, thus some calculation can also be done.
Calculations actually come from a force balance on a particle settling in the slurry.
A simple method is given below  valid for a single spherical particle
Terminal settling velocity is :
Equation 1 : Terminal settling velcity of a single particle
With the following nomenclature
s=(ρ_{p}/)  specific density of the particle in the fluid carrier
g=gravitational constant (ms2)
CD=Drag coefficient ()
dP=Diameter of particle (m)
The drag coefficient can be calculated depending on the values of the Reynolds particles.
In laminar regime (Re_{p}<2)
C_{D}=24/Re_{p}
Equation 2 : Drag coefficient of particle in laminar regime
In transition regime (2<Re_{p}<500)C_{D}=18.5/Re_{p}^{0.6}
Equation 3 : Drag coefficient of particle in transition regime
In turbulent regime (500<Re_{p}<20000)C_{D}=0.44
Equation 4 : Drag coefficient of particle in turbulent regime
The Reynolds particle can be calculated directly as a function of the Archimede number following
Equation 5 : Reynolds of particle
The value obtained for a single particle is actually greater than the actual velocity observed in a slurry, the different particles interacting between each other, which is reducing the settling speed. A method is proposed by Oroskar and Turian to account for this effect. The terminal settling velocity can be corrected the following way.
V_{sh}=V_{s}(1c_{v})^{n}
Equation 6 : Terminal settling velocity
With the following notation :
c_{v}=% volumetric concentration of solids in slurry
n is a power number function of the Reynolds particle
Figure 1 : n factor as a function of the Reynolds
Note : other correlations can be found in Mecanique et Rheologies des Fluides, Midoux, Lavoisier Tec et Doc, 1993, page 386
Now, knowing the speed at which the suspension will settle is of help in order to design the process. However, it does not answer the following question : when pumping the suspension into a pipe, will the suspension separate and 2 phases will finally appear. As a rule of thumb, it is found in litterature that problem will appear for suspension with a low level of solid (<25% volume) for which the viscosity and low and therefore the settling quick and / or with large particles (>100 microns).
The critical velocity in pipe under which phase separation can happen can be evaluated thanks to the following equation (Hanks 1986) :
Equation 7 : Critical pipe velocity
With the following parameters :
c_{v}=% volumetric concentration of solids in slurry
s=(ρ_{p}/)  specific density of the particle in the fluid carrier
The rheology of slurry can be quite complex since in most of the cases, the rheology will not be Newtonian.
A Newtonian fluid is charecterized by a viscosity that is independant of the shear rate. It is the case of water for example, and many other pure substances. Slurries, on their side, can present a viscosity that will change with the shear rate.
The different types of behaviors that can be found is represented on the following drawing. The shear stress is shown as a function of the shear rate. As reminder, the dynamic viscosity is defined as the ratio shear stress / shear rate
Figure 2 : Rheological behaviors
The different behavior can be better visualized by plotting the dynamic viscosity as a function of the shear rate.
Figure 3 : Rheological behaviors  Dynamic viscosity = f(shear rate)
A very common situation for slurries is to behave according to the pseudoplastic behavior, which means that the liquid will be shear thinning. It presents some advantages since the pressure drop will be less high in the pipes. It should however guide the choice for the pumping equipments for these fluids, a positive displacement pump is indeed probably a better choice than a centrifugal pump in which the viscosity may decrease a lot, preventing the pump to push the fluid in the pipe downstream.
A last rheological phenomenon is worth mentionning : the thixotropy and antithixotropy. Basically, there is a time hysteresis in the viscosity observed. A typical example of thixotropy are paints : agitating the paint will quickly reduce dramatically the viscosity, however, when stopping to agitate, the viscosity will increase again but slowly, giving time to use a very liquid paint.
Different rheological models are available : Powder law (Ostwald de Waele), Bingham, Casson, depending on the actual viscosity behavior.
In practice, a Power Law can be fitted in different parts of the viscosity plot as a function of the shear rate.
The law is in the following format :
Equation 8 : Viscosity powder law
What is important here, from the experience, if to calculate the exponent n.
The Reynolds number as well as the friction factor can then be calculated.
Laminar flow
Turbulent flow
Creating the suspension is the 1st step in a process handling a slurry. The type of agitator is important as well as the determination of the agitation speed necessary to maintain the particles in suspension.
An agitator performing a high shear is advised for contacting the solid and the liquid (typically a turbine with 45 degrees blades). It helps creating the turbulence necessary to overcome the setlling of particles and can also help to break some agregates that may be created.
If the suspension is very homogeneous, then this agitator can then be relayed by another one rotating slower (anchor), if the suspension is settling very quickly, it is probable that it is necessary to keep the turbine on.
In order to determine the minimal agitation speed for having a suspension, the litterature is often refering to the correlation of Zweitering (1958)
Equation 9 : Minimum agitation speed to put solids in suspension
With the following nomenclature
S=Zweitering coefficient, a function of the system geometry
ν=μ/ρ the slurry kinetic viscosity (m2.s1)
d_{p}=the particle diameter (m)
δρ=the difference between volumic weight of particle and of liquid abs(ρ_{p}ρ_{l}) (kg/m3)
D_{i}=agitator diameter (m)
X=mass solid concentration (% weight  please use percentage in the formula)
The values for S can be found in the litterature from different authors and for different agitators. Below are the values from Armenante et al (1988)
Equation 10 : Zweitering Coefficient
Table 1 : Parameters for Zweitering Coefficient Calculation
Coefficient  Disc Turbine vertical blades, rushton type (TPDD)  Turbine 6 vertcial blades (TPP)  Turbine 6 blades tilted, low pumping (TPIB)  HE3 Chemineer 

A  0.99  1.43  2.28  3.49 
a  1.40  1.20  0.83  0.79 
b  2.18  1.95  0.65  0.66 
H_{A} is the height from the bottom of the tank to the agitator (m)
D is the diameter of the agitator (m)
T is the diameter of the tank (m)
Do you want to know more ? www.powderprocess.net is covering all aspects of solids and powder processing and handling The pages below will probably interest you ! 
MIXING All on mixers and mixing principles 
PNEUMATIC TRANSPORT All on powder pneumatic conveying, the different technologies available and design methodologies 
You can find much more through the menus at the top of the page ! 