Terminal Velocity Particle
Terminal Velocity Particle
Terminal Velocity Particle
A balance of the gravitational, buoyancy and drag forces on the submerged solid
body determines a settling velocity of the body.
πd 3
(ρs − ρf )g = C D πd 2 v 2tsρf (I.1)
6 8
I.1
I.2 INTERMEZZO I
4 (ρs − ρ f ) gd
v ts = (I.2)
3 ρf CD
The drag coefficient CD is sensitive to a regime of the liquid flow round the settling
solid particle and this can be expressed by a relationship CD = fn(Rep). Rep is the
v d
particle Reynolds number Re p = ts .
νf
In a laminar regime (obeying Stokes' law and occurring for Rep < 0.1, i.e.
sand-density particles of d < 0.05 mm approximately) the relation is hyperbolic, CD =
24/Rep, so that
(ρ − ρf ) gd 2
v ts = s (I.3).
18ρ f νf
In a turbulent regime (obeying Newton's law and occurring for Rep >500, i.e.
sand-density particles of d > 2 mm approximately) the drag coefficient is no longer
dependent on Rep because the process of fast settling of coarse particles is governed
by inertial rather than viscous forces. CD = 0.445 for spherical particles and the
terminal velocity is given by
(ρs − ρf )
v ts = 1.73 gd (I.4).
ρf
These two regimes are connected via a transition regime. In this regime a CD value
might be determined from the empirically obtained curve CD = fn(Rep) (see Figure
1.8) or its mathematical approximation (Turton & Levenspiel, 1986)
24 0.657 + 0.413
CD = 1 + 0.173 Re p (I.5).
Re p 1 + 1.63x104 Re −1.09
p
The determination of vts for the transition regime requires an iteration [vts = fn(CD)
is an implicit equation]. Grace (1986) proposed a method for a determination of vts
without necessity to iterate.
SETTLING VELOCITY OF SOLID PARTICLE IN A LIQUID I.3
I.1.3 The Grace method for the solution of the settling-velocity equation
ρ (ρ − ρ f )g
- the dimensionless particle size: d* = d.3 f s
µ f2
ρ f2
- the dimensionless settling velocity: v *ts = v ts 3
µ f (ρs − ρ f )g
Those are mutually related as shown in Figure I.1. Thus using the curve and
rearranging gives directly the velocity vts as a function of particle diameter d. No
iteration is required.
The curve in Figure I.1. can be also described by analytic expressions (Table I.1.)
appropriate for a computational determination of vts according to Grace method.
Table I.1. Correlations for vts* as a function of d*, after Grace (1986).
Range Correlation
The above equations are suitable for the spherical particles. The particles of dredged
solids are not spherical. The non-spherical shape of a particle reduces its settling
velocity.
vt
ξ= (I.6)
v ts
ξ is a function of the volumetric form factor K (K=0.26 for sand, gravel) and the
ρ (ρ − ρ f )g
dimensionless particle diameter d* = 3 f s .d (see Figure I.2.). The
µ f2
terminal velocity for sand particles is typically 50-60 % of the value for the sphere of
the equivalent diameter (see Table I.3 in paragraph I.4).
Figure I.2. The shape factor, ξ, as a function of the dimensionless particle diameter
d* and the volumetric form factor K.
The curves in Fig. I.2 can be analytically approximated by the following correlation
0.0656
log10 ξ = −0.3073 + (I.7).
cosh 2.55log10 d* − 1.114
I.2.1.2 The estimation of the shape factor value for sand and gravel
A different method suggests the ξ value of about 0.7 for sand and gravel so that vt is
simply
vt = 0.7vts.
I.6 INTERMEZZO I
I.2.2 Empirical equations for settling velocity of sand and gravel particles
The above methods for a determination of the terminal velocity of spherical and
non-spherical particles are supposed to be valid for all sorts of solids settling in an
arbitrary fluid. In a dredging practice the terminal settling velocity is often determined
using simple equations found specifically for sand particles. These equations need no
iteration.
Warning: in the following equations the particle diameter, d, is in mm, and the
settling velocity, vt, in mm/s.
For the laminar region, considered for sand particles smaller than 0.1 mm, the Stokes
equation is written as
v t = 424
(Ss − Sf ) d 2 (I.8).
Sf
For the transition zone, 0.1 mm < d < 1 mm, the Budryck equation is used
8.925 (Ss − Sf ) 3
vt = 1 + 95 d − 1 (I.9).
d Sf
The turbulent regime of settling is estimated to occur for sand particles larger than 1
mm. This is described by the Rittinger equation
(Ss − Sf )
v t = 87 d (I.10).
Sf
500
100
50
20
10
Figure I.3a. Terminal settling velocity of sand & gravel particles using Stokes,
Budryck and Rittinger equations.
I.8 INTERMEZZO I
0,1
valsnelheid in water v (mm/s)
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1,0
klasse 1
2,0
3,0
4,0
5,0
.
8925
6,0 v= ⋅ 1 − 95 ⋅ ( 2.65 − 1) ⋅ d3 − 1 [Budryck]
7,0 d
8,0
9,0
10,0
v=
87
⋅ (
2,6 klasse 2
20,0 5−
1) ⋅
d
[Ri
30,0 ttin
ger
40,0 ]
50,0
60,0
70,0
80,0
90,0
v=
100,0
klasse 3
4 24
⋅ (2
200,0
,6 5
− 1)
300,0
⋅d
400,0
6
500,0
[ St o
600,0
700,0
ke s
800,0
900,0
]
1000,0
0,01
0,02
0,03
0,04
0,05
0,06
0,07
0,08
0,09
0,10
0,20
0,30
0,40
0,50
0,60
0,70
0,80
0,90
1,00
2,00
3,00
4,00
5,00
6,00
7,00
8,00
9,00
10,00
korrelgrootte d (mm)
Figure I.3b. Terminal settling velocity of sand & gravel particles using Stokes,
Budryck and Rittinger equations.
SETTLING VELOCITY OF SOLID PARTICLE IN A LIQUID I.9
v th = v t (1 − C v )m (I.11),
in which vth is the hindered settling velocity of solid particle, vt is the terminal
settling velocity of the solid particle and m is the empirical exponent related to the
particle Reynolds number Rep (see Table I.2).
vd
Re p = t m
νf
Rep ≤ 0.2 4.6
0.2 < Rep < 1 4.4Rep-0.03
1 < Rep < 500 4.4Rep-0.1
500 ≤ Rep 2.4
Theoretically, the validity of the Richardson & Zaki equation is limited by the
maximum solids concentration that permits solids particle settling in a particulate
cloud. This maximum concentration corresponds with the concentration in an
incipient fluidized bed (Cv of about 0.57). Practically, the equation was
experimentally verified for concentrations not far above 0.30.
I.10 INTERMEZZO I
Table I.3. Values of parameters associated with a particle settling for quartz sand
particles (ρs = 2650 kg/m3, shape factor K = 0.26) of various sizes in water at room
temperature.
I.5. REFERENCES
Grace, J.R. (1986). Contacting modes and behaviour classification of gas-solid and
other two-phase suspensions. Can. J. Chem. Eng., 64, 353-63.
Richardson, J.F. & Zaki, W.M. (1954). Sedimentation and fluidisation: Part I,
Trans. Instn Chem. Engrs, 32, 35.
Turton, R. & Levenspiel, O. (1986). A short note on the drag correlation for spheres.
Powder Technol., 47, 83-6.
Wallis, G.B. (1969). One Dimensional Two Phase Flow. McGraw-Hill.
Clift, R., Grace, J.R. and Weber, M.E. (1978). Bubbles, Drops and Particles.
Academic Press.
Coulson, J.M., Richardson, J.F., Backhurst, J.R. & Harker, J.H. (1996).
Chemical Engineering. Vol. 1: Fluid Flow, Heat Transfer and Mass Transfer.
Butterworth-Heinemann.
Wilson, K.C., Addie G.R., Sellgren, A. & Clift, R. (1997). Slurry Transport Using
Centrifugal Pumps. Chapman & Hall, Blackie Academic & Professional.
SETTLING VELOCITY OF SOLID PARTICLE IN A LIQUID I.11
CASE STUDY I
Estimate the terminal settling velocity, vt, and hindered settling velocity, vth, for
grains of fine sand (grain size d = 0.120 mm) and gravel (d = 6.0 mm) in the
quiescent column of water at room temperature (for calculations take water density ρf
= 1000 kg/m3 and kinematic viscosity νf = 10-6 m2/s). For both sand and gravel
grains the following values are typical:
solids density ρs = 2650 kg/m3
volumetric shape factor K = 0.26.
For hindered settling consider a cloud of particles settling in a water column. Settling
particles occupy 27% per cent of the total volume of the column so that the
volumetric concentration of solids Cv = 0.27.
Solution:
The Grace fitting correlations (Table I.1.) give for the terminal settling velocity (vts)
of a sphere of d = 0.120 mm the value vts = 11.21 mm/s while for a sphere of d = 6
mm the value vts = 554.13 mm/s. A very similar result is reached if the Grace graph
(Figure I.1.) is used:
1000x1650x 9.81
d* = 0.0001203 = 3.04 gives v*ts = 0.44 for the 0.120 mm sand
0.000001
and
1000x1650x 9.81
d* = 0.0063 = 152 gives v*ts = 22 for the 6 mm gravel.
0.000001
0.001x1650x 9.81
Then v ts = v*ts 3 that gives
1000 2
vts = 11 mm/s for fine sand and vts = 557 mm/s for medium gravel.
The velocity ratio ξ is determined using the fitting correlation (I.7). This gives ξ =
0.522 for the fine sand and ξ = 0.503 for the studied gravel. Again, very similar
results should be obtained using the ξ graph (Figure I.2.).
According to this method (see Eqs. I.8-I.10) the terminal settling velocity of a solid
grain vt is calculated directly, thus without intermediation of the velocity vts. The
equations
2650 − 1000
v t = 87 6 = 273.74 mm/s for the medium gravel of diameter d = 6 mm
1000
(Rittinger, Eq. I.10).
For gravel the result obtained from the Rittinger equation is very similar to that given
by Grace. For fine sand, however, the Grace method predicts considerably lower
terminal settling velocity than the Budryck equation.
Rep = 120 x 10-6 x 9.47 x 10-3 / 10-6 = 1.14 for our fine sand and
Rep = 6 x 10-3 x 0.274 / 10-6 = 1642 for our gravel.
Wallis correlation (Eq. I.12) provides m = 4.286 (fine sand) and m = 2.832 (gravel).