Biochemical Engineering

James M. Lee
Department of Chemical Engineering
Washington State University
Pullman, WA 99164-2714
[email protected]

Chapter 8 Sterilization.....................................1
8.1. Sterilization Methods.................................................................... 1
8.2. Thermal Death Kinetics................................................................ 2
8.3. Design Criterion ........................................................................... 3
8.4. Batch Sterilization ........................................................................ 4
8.5. Continuous Sterilization ............................................................... 7
8.6. Air Sterilization .......................................................................... 13
8.7. Nomenclature.............................................................................. 23
8.8. Problems ..................................................................................... 25
8.9. References................................................................................... 27

Last Update: August 10, 2001

8−ii Sterilization

© 2001 by James M. Lee, Department of Chemical Engineering, Washington

State University, Pullman, WA 99164-2710.
This book was originally published by Prentice-Hall Inc. in 1992.
You can download this file and use it for your personal study of the subject.
This book cannot be altered and commercially distributed in any form
without the written permission of the author.
If you want to get a printed version of this text, please contact James Lee.
All rights reserved. No part of this book may be reproduced, in any form or
by any means, without permission in writing from the author.
Chapter 8
Sterilization

Most industrial fermentations are carried out as pure cultures in which only
selected strains are allowed to grow. If foreign microorganisms exist in the
medium or any parts of the equipment, the production organisms have to
compete with the contaminants for the limited nutrients. The foreign
microorganisms can produce harmful products which can limit the growth of
the production organisms. Therefore, before starting fermentation, the
medium and all fermentation equipment have to be free from any living
organisms, in other words, they have to be completely sterilized.
Furthermore, the aseptic condition has to be maintained.
8.1. Sterilization Methods
Sterilization of fermentation media or equipment can be accomplished by
destroying all living organisms by means of heat (moist or dry), chemical
agents, radiation (ultraviolet or X-rays), and mechanical means (sonic or
ultrasonic vibrations). Another approach is to remove the living organisms by
means of filtration or high-speed centrifugation.
Heat is the most widely used means of sterilization, which can be
employed for both liquid medium and heatable solid objects. It can be applied
as dry or moist heat (steam). The moist heat is more effective than the dry
heat, because the intrinsic heat resistance of vegetative bacterial cells is
greatly increased in a completely dry state. As a result the death rate is much
lower for the dry cells than for moist ones. The heat conduction in dry air is
also less rapid than in steam. Therefore, dry heat is used only for the
sterilization of glassware or heatable solid materials. By pressurizing a
vessel, the steam temperature can be increased significantly above the boiling
point of water. Laboratory autoclaves are commonly operated at a steam
pressure of about 30 psia, which corresponds to 121°C. Even bacterial spores
are rapidly killed at 121°C.
Chemical agents can be used to kill microorganisms as the result of their
oxidizing or alkylating abilities. However, they cannot be used for the
sterilization of medium because the residual chemical can inhibit the
8-2 Sterilization

fermentation organisms. Chemical agents are frequently employed for

disinfection that commonly implies the treatment to remove or reduce the risk
from pathogenic organisms. Some of the major antimicrobial chemical agents
are (Pelczar and Reid, 1972): phenol and phenolic compounds (phenol,
cresol, orthophenylphenol), alcohol (ethyl, methyl), halogens (iodine,
hypochlorites, chloramines), detergents, dyes, quaternary ammonium
compounds, acids, alkalies, and gaseous chemosterilizers (ethylene oxide, β -
propiolactone, formaldehyde).
Many cellular materials absorb ultraviolet light, leading to DNA damage
and consequently to cell death. Wavelengths around 265 nm have the highest
bactericidal efficiency. However, ultraviolet rays have very little ability to
penetrate matter. Therefore, their use is limited to the reduction of microbial
population in a room where sterility needs to be maintained, such as hospital
operating rooms or clean chambers in a laboratory. X-rays are lethal to
microorganisms and have penetration ability. However, they are impractical
as sterilization tools due to their expense and safety concerns.
Sonic or ultrasonic waves of sufficient intensity can disrupt and kill cells.
This technique is usually employed in the disruption of cells for the purpose
of extracting intracellular constituents rather than as a sterilization technique.
Filtration is most effectively employed for the removal of microorganisms
from air or other gases. In the case of liquid solutions, it is used with
thermolabile medium or products, that is, those easily destroyed by heat, such
as human and animal serums and enzymes.
Among the techniques discussed, moist heat is the most economical and
efficient for the general sterilization requirements of fermentation. Therefore,
the following four sections describe cell death kinetics and sterilization
operations utilizing moist heat.
8.2. Thermal Death Kinetics
Thermal death of microorganisms at a particular temperature can be
described by first-order kinetics:
dn
= −kdn (8.1)
dt
where kd is specific death rate, the value of which depends not only on the
type of species but also on the physiological form of cells. For example, the
value of kd for bacterial spores at 121°C is of the order of 1 min−1, whereas
 Sterilization 8-3

those for vegetative cells vary from 10 to about 1010 min–1 depending on the
particular organism (Aiba et al., 1973).
Integration of Eq. (8.1) yields
n t
ln = −∫ kddt (8.2)
n0 0

or
( )
t
n = n0 exp −∫ kddt (8.3)
0

which shows the exponential decay of the cell population. The temperature
dependence of the specific death rate kd can be assumed to follow the
Arrhenius equation:
kd = kd0 exp − d
E
RT ( ) (8.4)

where Ed is activation energy, which can be obtained from the slope of the
ln(kd) versus 1/T plot. For example, the activation energy of E. coli is 127
kcal/gmole and that of Bacillus stearothermophilus Fs 7954 is 68.7 (Aiba et
al., 1973).
8.3. Design Criterion
From Eqs. (8.2) and (8.4), the design criterion for sterilization ∇ can be
defined as (Deindoerfer and Humphrey, 1959)
n E
( )
t t
∇ = ln 0 = ∫ kddt = kd0 ∫ exp − d dt (8.5)
n 0 0 RT
which is also known as the Del factor, a measure of the size of the job to be
accomplished. The Del factor increases as the final number of cells decreases.
For example, the Del factor to reduce the number of cells in a fermenter from
1010 viable organisms to one is
1010
∇ = ln = 23 (8.6)
1
The reduction of the number of cells from 1010 to one seems to be
impressive. However, even if one organism is left alive, the whole fermenter
may be contaminated. Therefore, all viable organisms have to be eliminated.
The Del factor to reduce the number of cells to zero is infinity, which means
that it is theoretically impossible to ensure the total destruction of the viable
cells. Therefore, the final number of cells needs to be expressed as the
fraction of one, which is equal to the probability of contamination. For
8-4 Sterilization

example, n = 0.001 means that the chance for a contaminant surviving the
sterilization is 1 in 1000. The Del factor to reduce the number of cells in a
fermenter from 1010 viable organisms to 0.001 is
1010
∇ = ln = 30 (8.7)
0.001
Based on the sterilization criterion calculated, we can design the sterilization
unit.
8.4. Batch Sterilization
Sterilization of the medium in a fermenter can be carried out in batch mode
by direct steam sparging, by electrical heaters, or by circulating constant
pressure condensing steam through heating coil. The sterilization cycles are
composed of heating, holding, and cooling. Therefore, the total Del factor
required should be equal to the sum of the Del factor for heating, holding and
cooling as
∇ total = ∇ heat + ∇ hold + ∇ cool (8.8)
The values of ∇heat and ∇cool are determined by the methods used for the
heating and cooling. The value of ∇hold is determined by the length of the
controlled holding period. The design procedure for the estimation of the
holding time is as follows:
1. Calculate the total sterilization criterion, ∇total.
2. Measure the temperature versus time profile during the heating,
holding, and cooling cycles of sterilization. If experimental
measurements are not practical, theoretical equations for heating and
cooling can be employed, which are of linear, exponential, or
hyperbolic form depending on the mode of heating and cooling. The
suggested equations for different heating and cooling processes are as
follows (Deindoerfer and Humphrey, 1959):
a. For batch heating by direct steam sparging into the medium, the
hyperbolic form is used:
Hmst
T = T0 + (8.9)
c ( M + mst )
b. For batch heating with a constant rate of heat flow such as electrical
heating, the linear form is used:
 Sterilization 8-5

qTt
T = T0 + (8.10)
cM
c. For batch heating with a isothermal heat source such as steam
circulation through heating coil, the exponential form is used:
T = TH + (T0 − TH ) exp  −
UAt 
 (8.11)
 cM 
d. For batch cooling using a continuous nonisothermal heat sink such
as passing cooling water through cooling coil, the exponential form
is used:
   UA   mct 
T = TC0 + (T0 − TC0 ) exp − 1 − exp  −  M  (8.12)
   m c
c  
3. Plot the values of kd as a function of time.
4. Integrate the areas under the kd-versus-time curve for the heating and
the cooling periods to estimate ∇heat and ∇cool, respectively. If using
theoretical equations, integrate Eq. (8.5) numerically after substituting
in the proper temperature profiles. Then, the holding time can be
calculated from
∇ ∇ + ∇ hold + ∇ cool
thold = total = heat (8.13)
kd kd
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
Example 8.1
A fermenter containing 40 m3 of medium (25°C) is going to be sterilized by
the direct injection of saturated steam. The typical bacterial count of the
medium is about 5×1012 m−3, which needs to be reduced to such an extent that
the chance for a contaminant surviving the sterilization is 1 in 1,000. The
steam (345 kPa, absolute pressure) will be injected with a flow rate of 5,000
kg/hr, which will be stopped when the medium temperature reaches 122°C.
During the holding time, the heat loss through the vessel is assumed to be
negligible. After a proper holding time, the fermenter will be cooled by
passing 100 m3/hr of 20°C water through the cooling coil in the fermenter
until the medium reaches 30°C. The coil has a heat-transfer area of 40 m2 and
for this operation the average overall heat-transfer coefficient (U) for cooling
is 2,500 kJ/hr m2 K. The heat-resistant bacterial spores in the medium can be
characterized by an Arrhenius coefficient kd0 of 5.7×1039 hr–1 and an
8-6 Sterilization

activation energy (Ed) of 2.834×105 kJ/kmol (Deindoerfer and Humphrey,

1959). The heat capacity and density of the medium are 4.187 kJ/kg K and
1,000 kg/m3, respectively. Estimate the required holding time.
Solution:
The design criterion can be calculated from Eq. (8.5) as
n ( )( )
∇ = ln 0 = ln  5 × 10 m −3 40m  = 39.8
12 -3 3

n  1 × 10 
The direct injection of steam into the medium can be assumed to follow
the hyperbolic temperature-time profile of Eq. (8.9), which can be used to
calculate the time required to heat the medium from 25°C to 122°C. From the
steam table (Felder and Rousseau, 1986), the enthalpy of saturated steam at
345 kPa and water at 25°C is 2,731 and 105 kJ/kg, respectively. Therefore,
the enthalpy of the saturated steam at 345 kPa relative to raw medium
temperature (25°C) is
H = 2,731 − 105 = 2,626 kJ/kg
From Eq. (8.9),
(2,626 kJ/kg)(5,000 kg/hr)t 78.4t
T = T0 + 3 3
= T0 +
(4.187 kJ/kg K)[(40 m )(1,000 kg/m ) + (5,000 kg/hr)t ] 1 + 0.125t

The solution of the preceding equation for t when T=395°K and T0 =

298°K by using a numerical technique or a trial and error approach yields
that the time required to reach 122°C is 1.46 hrs.
Substitution of the previous equation into Eq. (8.5),
 −2.834 ×105  −1
 
298 + 78.4t
1.46
∇heat = 5.7 ×10 ∫ exp 
39
  dt

0
 8.318  1 + 0.125 t  
Numerical integration of the preceding equation by using Advanced
Continuous Simulation Language (ACSL) or other method yields
∇heat = 14.8
During the cooling process, the change of temperature can be
approximated by Eq. (8.12) as
T = 293 + 102exp(−0.531t )
Solving for t when the final temperature is 303°K yields t = 4.38 hrs.
Substitution of the previous equation into Eq. (8.5) gives
 Sterilization 8-7

4.38  −2.834 ×105 

∫  dt = 17.6
39
∇cool = 5.7 ×10 exp
0  8.318 [293 + 102exp(−0.531 t )] 
 
Therefore, the Del factor for the holding time is
∇ hold = ∇ total − ∇ heat − ∇ cool = 39.8 − 14.8 − 17.6 = 7.4
At 122°C, the thermal death constant (kd) is 197.6 hr–1 from Eq. (8.4).
Therefore, the holding time is
∇ 7.4
thold = hold = = 0.037 hr = 2.25 min
kd 197.6
For this example, most of the sterilization was accomplished during the
heating and cooling periods.
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
8.5. Continuous Sterilization
Sterilization can be carried out in a continuous mode rather than in batches.
Continuous sterilization offers several advantages:
1. It simplifies production planning, thus allowing maximum plant
utilization and minimum delays.
2. It provides reproducible conditions.
3. It can be operated at a high temperature (140°C instead of 121°C in
batch sterilization); therefore, the sterilization time can be shortened
(holding time of 1 to 2 minutes).
4. It requires less steam by recovering heat from the sterilized medium.
As a result, it also requires less cooling water.
5. It is easier to automate the process; thus, it is less labor intensive.
A continuous sterilizer consists of three main sections: heating, holding, and
cooling.
Heating Section: Methods of heating can be categorized into two types:
direct steam injection and indirect heating in shell-and-tube or plate-and-
frame heat exchanger. Direct heating is more effective than indirect heating
because there is no barrier between the medium and the heat source. The
steam injector heats the medium to the peak sterilization temperature quickly.
Therefore, sterilization during the heating period is negligible.
8-8 Sterilization

For indirect heating, the plate-and-frame heat exchanger is generally more

effective than the shell-and-tube type for heat transfer due to its larger heat-
transfer area. However, the former is limited to lower pressures (normally
less than 20 atm) due to its weak structural strength compared with the latter.
The plate-and-frame type is also favorable for the sterilization of a high
viscous system.
The temperature change with respect to residence time (τ hold ) as the
medium passes through an isothermal heat source can be approximated as
(Deindoerfer and Humphrey, 1959b),
 UAτ heat 
TC2 = TH − (TH − TC1 )exp −  (8.14)
 cW 
For heating using a countercurrent heat source of equal flow rate and heat
capacity,
∆TUAτ heat
TC2 = TC1 − (8.15)
cW
Holding Section: The heated medium passes through a holding section,
which is usually composed of a long tube. The holding section is maintained
in adiabatic conditions. If the heat loss in the section is negligible, the
temperature can be assumed to be constant. The average residence time in the
holding section is
L
τ hold = (8.16)
u
from which the Del factor can be estimated as
n  E 
∇hold = ln 0 = kdτ hold = kd0 exp− d τ hold (8.17)
n  RT 

where n0 is the number of cells at the beginning of the holding section.

If the medium in a holding section behaves as ideal plug flow, the
residence time of the medium in the section will be exactly the same for all
the medium. Therefore, the degree of the sterilization will be uniform.
However, the slippage due to the viscous nature of fluid, the friction of the
pipe wall, and turbulent eddies of the flowing fluid causes the deviation of
the ideal plug flow. The resulting velocity profile has its maximum value at
the centerline of a pipe, whereas it has its minimum value in the vicinity of
the pipe wall. For laminar flow of Newtonian fluids through a smooth round
pipe, the ratio of the average fluid velocity to the maximum velocity /line
 Sterilization 8-9

u/um is 0.5. The ratio changes rapidly from 0.5 to about 0.75 when laminar
flow changes to turbulent, and then increases gradually to 0.87 when the
Reynolds number is about 106 (McCabe et al., 1985). As a result, if we use
the mean velocity in calculating the required residence time for sterilization,
some portion of medium will be understerilized, which may cause a serious
contamination problem.
.............
.............
.............
.............
.............
................................................. ............. .......................................................
. d(uCn )
Bulk flow uCn S . . . ..
............. uCn S + dxS
.............
.............
dx
.............
.............

dCn .............
dC dC
Dispersion −DS .................................................... .................... .....................................................
.
...................
−DS n + d −DS n dx
dx ............. dx dx dx
.............
.............

........................................

dx
Figure 8.1 Material balances around an elementary section in a holding
tube.
The deviation from ideal plug flow due to the axial mixing can be
described by the dispersion model (Levenspiel, 1972). Let's look at the
differential element with a thickness dx in a holding tube as shown in Figure
8.1. The basic material balance for the microorganisms suspended in the
medium is
In - Out -Killed by Sterilization = Accumulation (8.18)
At steady state, the accumulation term is equal to zero. The input and output
of the microorganisms into or out of the element have both a bulk flow and
an axial diffusion condition. The number of microorganisms entering minus
those leaving by bulk flow is
 d (uCn ) 
uCn S − uCn S + Sdx  (8.19)
 dx 
 
Analogous to the molecular diffusion, the x-directional flux of
microorganisms suspended in a medium due to the axial mixing can be
represented as
dC
J n = −D n (8.20)
dx
where D is the axial dispersion coefficient, characterized by the degree of
backmixing during flow. The mechanism of axial dispersion may be
molecular or turbulent. If D is zero, the velocity distribution approaches that
of the ideal plug flow. At the other extreme, if D is infinitely large, the fluid
in the pipe is well mixed like a fully mixed vessel. For the turbulent flow, the
8-10 Sterilization

10.0

D
1.0
udt

0.1
1.0E+03 1.0E+04 1.0E+05 1.0E+06

dt u ρ
N Re =
µL

Figure 8.2 Correlation for D/udt as a function of Reynolds number.

(Levenspiel, 1958).

dispersion coefficient is correlated as a function of Reynolds number as

shown in Figure 8.2.
The number of microorganisms entering and leaving by axial dispersion is
dC  dC d 
− DS n −  − DS n +  − DS n  dx 
dC
(8.21)
dx  dx dx  dx  
The number of cells killed by sterilization is kdCnSdx. Therefore, by
substituting Eqs. (8.19) and (8.21) into Eq. (8.18) and simplifying,
d  dCn  d (uCn )
D − − kd Cn = 0 (8.22)
dx  dx  dx
For the constant D and u , Eq. (8.22) can be modified into a dimensionless
form,
d 2Cn′ dCn′ k L
− N − N Pe d Cn′ = 0 (8.23)
dx′ dx′
2 Pe
u
where
 Sterilization 8-11

Cn x uL
Cn′ = x′ = N Pe =
Cn0 L D

NPe is known as Péclet number. When NPe = ∞, it is ideal plug flow. The
boundary conditions for the solution of Eq. (8.23) are
dCn′
+ N Pe (1 − Cn′ ) = 0 at x′ = 0
dx′
(8.24)
dCn′
=0 at x′ = 1
dx′
The solution of Eq. (8.23) (Wehner and Wilhelm, 1956) is
4ϕ exp  Pe 
N
Cn′ x′=1 =  2  (8.25)
(1 + ϕ ) exp  ϕ N Pe  − (1 − ϕ ) exp  − ϕ N Pe 
2 2

 2   2 
where
4k d L / u
ϕ = 1+ (8.26)
N Pe
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
Example 8.2
A continuous sterilizer with a steam injector and a flash cooler will be
employed to sterilize medium continuously with the flow rate of 2 m3/hr. The
time for heating and cooling is negligible with this type of sterilizer. The
typical bacterial count of the medium is about 5×1012 m−3, which needs to be
reduced to such an extent that only one organism can survive during two
months of continuous operation. The heat-resistant bacterial spores in the
medium can be characterized by an Arrhenius coefficient (kd0) of 5.7×1039
hr–1 and an activation energy (Ed) of 2.834×105 kJ/kmol (Deindoerfer and
Humphrey, 1959). The sterilizer will be constructed with the pipe with an
inner diameter of 0.102 m. Steam at 600 kPa (gage pressure) is available to
bring the sterilizer to an operating temperature of 125°C. The physical
properties of this medium at 125°C are c = 4.187 kJ/kg K, ρ = 1000 kg/m3,
and µ = 4 kg/m hr.
a. What length should the pipe be in the sterilizer if you assume ideal
plug flow?
8-12 Sterilization

b. What length should the pipe be in the sterilizer if the effect of axial
dispersion is considered?

Solution:
a. The design criterion can be calculated from Eq. (8.5) as
n
∇ = ln 0 = ln (5 × 1012 m −3 )( 2 m3 /hr )( 24 hr/day )( 60 days ) = 37.2
n
Since the temperature at the holding section is constant, Eq. (8.5) is
simplified to
∇ = k dτ hold
From the given data, kd can be calculated by using Eq. (8.4) to yield
kd = 378.6 hr −1
Therefore,
∇ 37.2
τ hold = = = 0.0983 hr
kd 378.6
The velocity of medium is
2 m3 / hr
u= = 245 m/hr
π 2 2
0.102 m
4
The length of the sterilizer is
L = uτ hold = 24.1 m
b. The Reynolds number for the medium flow is
0.102 (245)(1000)
N Re = = 6.24 × 103
4
From Figure 8.2
D
≈ 0.8 for N Re = 6.24 × 103
udt
Therefore,
D ≈ 0.8udt = 20 m 2 /h
Now, substituting all the values given and calculated in this problem to
Eq. (8.25) will result in an equation with only one unknown, L, which
can be solved by using any non-linear equation solver:
 Sterilization 8-13

L = 26.8
Therefore, the holding section should be 26.8 m, which is 2.7 m longer
than the result from the assumption of ideal plug flow.
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
Cooling Section: For the cooling section, a quench cooler with adequate
heat removal capacity is effective. Another technique is to inject the hot
medium through an expansion valve into a vacuum chamber, which is known
as flash cooling. Both of these take a very short time; therefore, the
sterilization during the cooling period can be assumed to be negligible.
A shell-and-tube or a plate-and-frame heat exchanger can also be
employed for cooling. The temperature versus residence time relationship for
cooling using an isothermal heat sink is
 UAτ cool 
TH 2 = TC − (TC − TH1 )exp − (8.27)
 cW 

For cooling using a countercurrent heat sink of equal flow rate and heat
capacity,
∆TUAτcool
TH 2 = TH1 − (8.28)
cW
8.6. Air Sterilization
For aerobic fermentations, air needs to be supplied continuously. Typical
aeration rates for aerobic fermentation are 0.5−1.0 vvm (air volume per liquid
volume per minute). This requires an enormous amount of air. Therefore, not
only the medium but also the air must be free of microbial contaminants. All
of the sterilization techniques discussed for medium can also be employed for
air. However, sterilization of air by means of heat is economically impractical
and is also ineffective due to the low heat-transfer efficiency of air compared
with those of liquids. The most effective technique for air sterilization is
filtration using fibrous or membrane filters.
The cotton plug, routinely used as a closure for tubes or flasks of sterile
solution, is a good example of removal of microorganisms from air by a
fibrous filter. A simple air filter can be made by packing cotton into a
column. However, with cotton filters the pressure drop is high and wetting
can be a breeding ground for the contamination. Therefore, glass fibers are
favorable as filter medium because they give a lower pressure drop and are
less liable to wetting or combustion. Modern fibrous filter systems are
8-14 Sterilization

cylinders made from bonded borosilicate microfibers sheathed in reinforcing

polypropylene mesh in which the layers increase in fineness and density from
the center outward. This type of design can deliver over 3 m3/s of sterile air at
0.1 bar of pressure drop (Quesnel, 1987).
With fibrous filters, airborne particles are collected by the mechanisms of
impaction, interception, and diffusion.
Impaction: When an air stream containing particles flows around a
cylindrical collector, the particle will follow the streamlines until they
diverge around the collector. The particles because of their mass will have
sufficient momentum to continue to move toward the cylinder and break
through the streamlines, as shown in Figure 8.3. The collection efficiency by
this inertial impaction mechanism is the function of the Stokes and the
Reynolds number as:
 C f ρ p d p2 v0 Dc v0 ρ 
ηimp = f ( NSt , N Re ) =   (8.29)
 18 µ D c , µ 
where Cf is known as Cunningham correction factor. The value of Cf can be
estimated from the empirical correlation developed by Davies (Strauss,
1975),
2λ   d 
Cf =1+ 1.257 + 0.400exp  −1.10 p   (8.30)
dp   2λ  
where λ is the mean free path of gas molecules based on the Chapman-
Enskog equation,
 µ  πMw
λ =  (8.31)
 0.499 ρ  8RT

Figure 8.3 Flow pattern around cylindrical fiber, showing the path of
particles collected by inertial impaction.
 Sterilization 8-15

The efficiency ηimp is defined as the fraction of particles approaching the

collector which impact. Various correlations are available in the literature.
An empirical correlation for the efficiency developed by Thom is (Strauss,
1975):
3
NSt
ηimp = 3 2
for N Rec = 10 (8.32)
NSt + 0.77 NSt + 0.22
Another correlation proposed by Friedlander (1967) is
ηimp = 0.075NSt
1.2
(8.33)

The efficiency increases with increasing particle diameter or air flow

velocity.
Interception: The inertial impaction model assumed particles had mass,
and hence inertia, but no size. An interception mechanism is considered
where the particle has size, but no mass, and so they can follow the
streamlines of the air around the collector. If a streamline which they are
following passes close enough to the surface of the fiber, the particles will
contact the fiber and be removed (Figure 8.4). The interception efficiency
depends on the ratio of the particle diameter to the cylindrical collector
diameter (κ =dp/Dc):
1  κ (2 + κ ) 
ηint = (1 + κ ) ln (1 + κ ) − (8.34)
2.002 − ln N Rec  2 (1 + κ ) 

which was developed by using Langmuir's viscous flow equation (Strauss,

1975). The ratio κ is known as interception parameter. The collection
efficiency by interception increases with the increase of the particle size.

Figure 8.4 Flow pattern around cylindrical fiber, showing the interception
collection mechanism.
8-16 Sterilization

Diffusion: Particles smaller than about 1 micron in diameter exhibit a

Brownian motion which is sufficiently intense to produce diffusion. If a
streamline containing these particles is sufficiently close to the collector, the
particles may hit the collector and be removed. Contrary to the previous two
mechanisms, the collection efficiency by diffusion increases with decreasing
particle size or air velocity. The typical size of particles collected by this
mechanism is less than about 0.5 µm. The efficiency of collection by
diffusion can be estimated by an equation analogous to Langmuir's equation,
Eq. (8.34), as (Strauss, 1975):
1  Z (2 + Z ) 
ηdif = (1 + Z ) ln (1 + Z ) − (8.35)
2.002 − ln N Rec  2 (1 + Z ) 

where Z is the diffusion parameter defined as

1
 D 
Z =  2.24 ( 2.002 − ln N Rec ) Br 
3
(8.36)
 vDc 
Friedlander (1967) suggested the following correlation
ηdif = 1.3N Pe−(2 / 3) + 0.7κ 2 (8.37)
where NPe is Péclet number, an important dimensionless parameter in the
theory of convective diffusion. It is defined as
vD
N Pe = c c = N Re NSc (8.38)
DBr
where NSc is the Schmidt number which is defined as
µ
NSc = (8.39)
ρ DBr
The diffusivity due to Brownian motion for submicron size particles can be
estimated from
C kT
DBr = f (8.40)
3πµ d p

where k is Boltzmann's constant (1.38054×10−9 J/K).

Combined Mechanisms: The total collection efficiency of a fibrous filter
is obtained from the combined effect of the preceding three mechanisms. One
straightforward way to combine the collection efficiencies of the different
mechanisms is to add them together, but this implies that a particle can be
 Sterilization 8-17

collected more than once, which does not make sense. A better approach is to
use the following correlation,
ηc = 1 − (1 − ηimp )(1 − ηint )(1 − ηdif ) (8.41)
which allows only the particles not collected by one mechanism to be
collected by the others. Substitution of Eqs. (8.32), (8.34), and (8.35) into
Eq. (8.41) will result in the correlation for the collection efficiency by the
combined mechanisms. Pasceri and Friedlander (1960) correlated the
combined collection efficiency as
6
ηc = 2/3 0.5 + 3κ 2 N Re
0.5
(8.42)
NSc N Rec c

As mentioned earlier, with an increase of the superficial air velocity (v0), ηimp
and ηint increase whereas ηdif decreases. Therefore, the combined collection
efficiency normally decreases to reach a minimum point and then increases
with increasing superficial air velocity.
Effect of Multiple Layers and Packing: All correlations for the collec-
tion efficiency discussed so far are based on the ideal case of a single
cylindrical collector. Now, let's examine a filter unit consisting of randomly
oriented multiple layers. Consider an area (A) of filter at a right angle to the
gas flow and with a depth dh. If the packing density α is defined as the
volume of fiber per unit volume of filter bed, the velocity within the filter
void space is equal to
v0
v= (8.43)
(1 − α )

 
d  n 0  
Cv
v 
Cn 0 A (1 − α ) v  1 − α  dh A (1 − α )
1−α Cn 0 + 
 1−α dh 

dh

Figure 8.5 Shell balances around a differential element of a filter.

A mass balance on the particles for the control volume (Figure 8.5) results
in
Input - Output = Collected by the filter (8.44)
8-18 Sterilization

  Cnv0  
v0 d    A(
Cn A (1 − α ) −  v0  1 − α  dh 1−α )
1−α Cn + 
 1−α dh  (8.45)
v
= 0 Cn ( Adh)ηc Dc L
1− a
where L is the length of cylindrical fiber per unit volume of filter bed, which
is related to the packing density α and the average collector diameter Dc as
π Dc2 L
α= (8.46)
4
Simplifying Eq. (8.45) and substituting Eq. (8.46) for L gives
dC 4αηc
− n= dh (8.47)
Cn π Dc (1 − α )
which, on integration, results
Cn 4B  α 
ln =−  ηc (8.48)
Cn0 π Dc  1 − α 
where B is the filter depth. Therefore, the collection efficiency for the filter
bed can be estimated as
C  4B  α  
η f = 1 − n = 1 − exp  −  ηc  (8.49)
Cn0  π Dc  1 − α  

When fibers are packed together in a filter bed, the velocity will be
increased and the flow pattern will be changed, which increases the
collection efficiency from impaction and interception. Chen (1955) has
determined fiber interference effects experimentally and suggests
ηα = η f (1 + 4.5α ) (8.50)

which is applicable for α < 0.1 and ηf < 1/(1+4.5 α).

In summary, in order to estimate the collection efficiency of a filter bed,
you have to calculate: ηc by using either Eq. (8.41) or Eq. (8.42), ηf using
Eq. (8.49), and ηα using Eq. (8.50). However, it should be noted that the
predictions of the collection efficiency from various correlations vary widely
due to empiricism or the oversimplification in developing the models
represent a complex collection mechanism. Furthermore, since the effect of
the gas velocity on the collection efficiency is large, the collection efficiency
can decrease significantly by increasing or decreasing the gas velocity. To
 Sterilization 8-19

insure sterility in a fermenter system, a conservative approach needs to be

taken which considers the minimal efficiency conditions due to possible
velocity fluctuations and prediction error from various correlations. In
assessing the filtration job to be accomplished, Humphrey (1960)
recommended that the design should permit only a one-in-a-thousand chance
of a single contaminant penetrating the filter during the entire course of the
fermentation.
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
Example 8.3
A filter bed of glass fibers (Dc = 15 µm, the bed depth B =10 cm, and
packing density α = 0.03) is being used to sterilize air (20°C, 1 atm) with an
undisturbed upstream velocity, v0, of 10 cm/s. The air stream contains 5,000
bacteria per cubic meter (dp = 1 µm and ρp = 1 g/cm3).
a. Estimate the single fiber collection efficiency by inertial impaction, by
interception, and by diffusion.
b. Estimate the single fiber collection efficiency based on combined
mechanisms by using Eq. (8.41) and Eq. (8. 42) and compare the
results.
c. Estimate the collection efficiency (ηα) of the filter bed.
d. Show how the superficial velocity v0 affects the various single fiber
collection efficiencies.
Solution:
a. The velocity within the filter void space is from Eq. (8.43)
v0 10
v= = = 10.3 cm/s
(1 − α ) 1 − 0.03
To estimate ηimp by using Eq. (8.32) or Eq. (8.33), we need to calculate
the Reynolds and the Stokes number from the given condition and the
physical properties of air at 20°C and 1 atm (ρ = 1.20×10−3 g/cm3, µ =
1.8 × 10−4 g/cm s).
8-20 Sterilization

Dc v ρ (1.5 × 10−3 )(10.3)(1.20 × 10−3 )

N Rec = = = 0.103
µ 1.80 × 10−4
where all units are in the cgs system. The mean free path λ and the
Cunningham correction factor can be estimated from Eq. (8.31) and
Eq. (8.30) as:
 1.80 × 10−4  π (29)
λ=  −3 
= 6.50 × 10−6 cm
 0.499 (1.20 × 10 )  8(8.314 × 10 )(293)
7

2λ   d p 
Cf =1+ 1.257 + 0.400exp  −1.10   = 1.16
d p   2λ  
The Stokes number is
C f ρ p d p2v 1.16(1) (1 × 10−4 ) (10.3)
2

NSt = = = 0.0247
18µ Dc 18 (1.80 × 10−4 )(1.5 × 10−3 )
Therefore, from Eq. (8.33), the single fiber collection efficiency by
inertial impaction is
ηimp = 0.075 NSt1.2 = 0.075 (0.0247 )1.2 = 8.82 × 10−4
The single fiber collection efficiency by interception is from Eq. (8.34)
for κ = d p Dc = 0.0667 and N Rec = 0.103 ,
1  κ (2 + κ ) 
ηint =  (1 + κ ) ln (1 + κ ) −  = 0.96 × 10−4
2.002 − ln N Rec  2 (1 + κ ) 

The diffusivity due to the Brownian motion can be estimated from

Eq. (8.40),
C f kT 1.16 (1.38 × 10−16 ) 293 δ y ∂ 2Ω
DBr = = = 2. 77 × 10 −5 cm 2 / s
3πµ d p 3π (1.80 × 10−4 )(1 × 10−4 ) δ x ∂u 2

Therefore, the Péclet number is

Dcv 1.5 × 10−3 (10.3)
N Pe = = −5
= 5.57 × 104
DBr 2.77 × 10
The single fiber efficiency by diffusion can be estimated from
Eq. (8.37),
 Sterilization 8-21

ηdif = 1.3N Pe−(2 / 3) + 0.7κ 2 = 4.00 × 10−3

Therefore, the single fiber collection efficiency by inertial impaction,
interception, and diffusion is 8.82×10−4, 9.96×10−4, and 4.00×10−3,
respectively.
b. The combined collection efficiency can be estimated from Eq. (8.41)
as:
ηc = 1 − (1 − 8.82 × 10−4 )(1 − 9.96 × 10−4 )(1 − 4.00 × 10−3 ) = 0.0059
Instead, if we use Eq. (8.42),
µ 1.80 × 10−4
NSc = = = 5.41 × 105
ρ DBr (1.20 × 10 )( 2.77 × 10 )
−3 −7

Therefore,
6
ηc = + 3(0.0667 )2 (0.103)0.5 = 0.0071
(5.41× 10 ) (0.103)
5 2/3 0.5

In this example, diffusion is predominant over impaction and

interception. If we use Eq. (8.32) and Eq. (8.35) for the estimation of
ηimp and ηdif, respectively, ηimp = 6.82×10−5 and ηdif = 6.97×10−4, which
are significantly lower than those predicted from Eq. (8.33) and
Eq. (8.37). Strauss (1975) compared the predicted values from various
correlations and found that the diffusion collection efficiency ηdif from
Eq. (8.35) tends to predict the lower value than the experimental
values.
c. If we choose to use ηc = 0.0059 which was predicted from Eqs. (8. 33),
(8.34), (8.37), and (8.41), the collection efficiency for the filter bed can
be estimated from Eq. (8.49) as
 4 (10 )  0.03  
η f = 1 − exp  − −3   0.0059  = 0.79
 π (1.5 × 10 )  1 − 0.03  
and the collection efficiency including the interference effects is from
Eq. (8.50),
ηα = 0.79 1 + 4.5 ( 0.03)
d. Similarly, we can calculate various collection efficiencies at various
values of v0. As shown in Table 8.1, with an increase of v0, ηimp and ηint
increased whereas ηdif decreased. Therefore, the combined collection
8-22 Sterilization

efficiency ηc decreased initially, reached a minimum value, and then

increased with increasing v0. The minimum collection efficiency was
reached when v0 was about 5 cm/s.
Table 8.1
Various Collection Efficiencies for 1 µm Particle
v0 ηimp×103 ηint×103 ηdif×103 ηc×103 ηc×103
cm/s Eq. (8.33) Eq. (8.34) Eq. (8.37) Eq. (8.41) Eq. (8.42)
0.1 0.004 0.479 22.3 22.8 28.6
1 0.056 0.647 7.25 7.941 0.3
5 0.0384 0.857 4.53 5.76 7.01
10 0.883 0.996 4.00 5.87 7.10
15 1.44 1.10 3.79 6.32 7.54
20 2.03 1.19 3.67 6.88 8.05
50 6.09 1.60 3.42 11.1 10.8
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
 Sterilization 8-23

8.7. Nomenclature
A surface area across which heat transfer occurs during sterilization,
m2
B filter bed depth, m
Cf Cunningham correction factor, dimensionless
Cn cell number density, number of cells/m3
c specific heat of medium, J/kg K
Dc collector diameter, m
dp particle diameter, m
dt pipe diameter, m
D axial dispersion coefficient, m2/s
DBr diffusivity due to the Brownian motion, m2/s
Ed activation energy for thermal cell destruction in Arrhenius equation,
J/kmol
H enthalpy of steam relative to raw medium temperature, J/kg
Jn flux of microorganisms due to axial dispersion, m−2s−1
k Boltzmann's constant, 1.3805×10−23 J/K or 1.3805×10−16 erg/K
kd specific death rate, s−1 or kg/m3s
L length of holding section, m
M initial mass of medium in batch sterilizer, kg
Mw molecular weight of gas molecules, kg/kmol
ms steam mass flow rate, kg/s
mc coolant mass flow rate, kg/s
NPe Péclet number (/line uL / D or v0 Dc / DBr ), dimensionless
NRe Reynolds number ( dt u ρ / µ L ), dimensionless
NRec collector Reynolds number ( Dcv0 ρ / µ ), dimensionless
NSc Schmidt number (µ/ρDBr), dimensionless
NSt Stokes number (Cf ρp dp2 v0/18 µDc), dimensionless
n number of cells in a system
q rate of heat transfer, J/s
R gas constant, 8.314×103 J/kmol K or 8.314×107 erg/mol K
S cross-sectional area of a pipe, m−2
8-24 Sterilization

T absolute temperature, °K
T0 initial absolute temperature of medium, °K
TC absolute temperature of heat sink, °K
TC0 initial absolute temperature of heat sink, °K
TH absolute temperature of heat source, °K
t time, s
td doubling time, s
U overall heat-transfer coefficient, J/s m2 K.
u velocity, m/s
u velocity, m/s
um maximum velocity, m/s
v fluid velocity within the filter void space, m/s
v0 undisturbed upstream fluid velocity, m/s
W mass of medium in a sterilizer, kg
x x-directional distance, m
Z the diffusion parameter defined in Eq. (8.36), dimensionless
α packing density defined as the volume of fiber per unit volume of
filter bed, dimensionless
η collection efficiency, dimensionless
κ the ratio of the particle and the collector diameter (dp/Dc), dimen-
sionless
λ the mean free path of gas molecules, m
µ fluid viscosity, kg/m s
µL liquid viscosity, kg/m s
ρ density, kg/m3
ρp density of particles, kg/m3
τ residence time, s
τ average residence time, s
∇ design criterion for sterilization, dimensionless
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
 Sterilization 8-25

8.8. Problems
8.1 A fermenter containing 10 m3 of medium (25°C) is going to be
sterilized by passing saturated steam (500 kPa, gage pressure)
through the coil in the fermenter. The typical bacterial count of the
medium is about 3×1012 m−3, which needs to be reduced to such an
extent that the chance for a contaminant surviving the sterilization
is 1 in 100. The fermenter will be heated until the medium reaches
115°C. During the holding time, the heat loss through the vessel is
assumed to be negligible. After the proper holding time, the
fermenter will be cooled by passing 20 m3/hr of 25°C water through
the coil in the fermenter until the medium reaches 40°C. The coils
have a heat-transfer area of 40 m2 and for this operation the average
overall heat-transfer coefficient (U) for heating and cooling are
5,500 and 2,500 kJ/hr m2K, respectively. The heat resistant bacterial
spores in the medium can be characterized by an Arrhenius
coefficient (kd0) of 5.7×1039 hr–1 and an activation energy (Ed) of
2.834×105 kJ/kmol (Deindoerfer and Humphrey, 1959). The heat
capacity and density of the medium are 4.187 kJ/kgK and 1,000
kg/m3, respectively. Estimate the required holding time.
8.2 A continuous sterilizer with a steam injector and a flash cooler will
be employed to sterilize medium continuously. The time for heating
and cooling is negligible with this type of sterilizer. The typical
bacterial count of the medium is about 5×1012 m−3, which needs to
be reduced to such an extent that only one organism can survive
during the three months of continuous operation. The heat resistant
bacterial spores in the medium can be characterized by an Arrhenius
coefficient (kd0) of 5.7×1039 hr–1 and an activation energy (Ed) of
2.834×105 kJ/kmol (Deindoerfer and Humphrey, 1959). The
holding section of the sterilizer will be constructed with 20 m-long
pipe with an inner diameter of 0.078 m. Steam at 600 kPa (gage
pressure) is available to bring the sterilizer to an operating
temperature of 125°C. The physical properties of this medium at
125°C are c =4.187 kJ/kgK and ρ =1,000 kg/m3, and µ = 4 kg/m hr.
a. How much medium can be sterilized per hour if you assume
ideal plug flow?
b. How much medium can be sterilized per hour if the effect of
axial dispersion is considered?
8-26 Sterilization

8.3 You need to design a filter for a 10,000-gallon fermenter that will
be aerated at a rate of 535 ft3/min (at 20°C and 1 atm). The bacterial
count in the air is 80 per ft3. Average size of the bacteria is 1 µm
with density of 1.08 g/cm3. You are going to use glass fibers (Dc =
15 µm) with packing density α = 0.03. The cross-sectional area of
the filter will be designed to give a superficial air velocity v0 of 5
ft/s.
a. What depth of the filter would you recommend to prevent
contamination?
b. How is the answer in (a) changed if v0 is decreased to 1 ft/s?
Explain the results.
8.4 Based on the combined mechanisms of impaction, interception, and
diffusion, a minimum efficiency should result for spheres
depositing on cylindrical collectors.
a Develop an expression for the particle size corresponding to this
minimum efficiency. As an approximation, ignore the effect of
particle size on the Cunningham correction for slip.
b. What is the particle size corresponding to the minimum
efficiency for the filter bed of glass fibers described in Example
8.3?
––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––
 Sterilization 8-27

8.9. References
Aiba, S., A. E. Humphrey and N. F. Millis, Biochemical Engineering (2nd
ed.), pp. 242−246. Tokyo, Japan: University of Tokyo Press, 1973.
Chen, C. Y., “Filtration of Aerosols by Fibrous Media,” Chem. Rev. 55
(1955):595−623.
Deindoerfer, F. H. and A. E. Humphrey, “Analytical Method for
Calculating Heat Sterilization Time,” Appl. Micro. 7 (1959a):256−264.
Deindoerfer, F. H. and A. E. Humphrey, “Principles in the Design of
Continuous Sterilizers,” Appl. Micro. 7 (1959b):264−270.
Felder, R. M. and R. W. Rousseau, Elementary Principles of Chemical
Processes (2nd ed.) pp. 630−635. New York, NY: John Wiley & Sons,
1986.
Friedlander, S. K., “Aerosol Filtration by Fibrous Filters,” in Biochemical
and Biological Engineering Science, vol 1., ed. N. Blakebrough.
London, England: Academic Press, Inc., 1967, pp. 49−67.
Humphrey, A. E., “Air Sterilization,” Adv. Appl. Micro. 2
(1960):301−311.
Levenspiel, O., “Longitudinal Mixing of Fluids Flowing in Circular
Pipes,” Ind. Eng. Chem. 50 (1958):343−346.
Levenspiel, O., Chemical Reaction Engineering (2nd ed.), p. 272. New
York, NY: John Wiley & Sons, 1972.
McCabe, W. L., J. C. Smith, and P. Harriott, Unit Operations of Chemical
Engineering (4th ed.), pp. 76−90. New York, NY: McGraw-Hill Book
Co., 1985.
Pasceri, R. E. and S. K. Friedlander, “The Efficiency of Fibrous Aerosol
Filters,” Can. J. Chem. Eng., 38 (1960):212−213.
Pelczar, M. J. and R. D. Reid, Microbiology, pp. 441−461. New York,
NY: McGraw-Hill Book Co., 1972.
Quesnel, L. B., “Sterilization and Sterility,” in Basic Biotechnology, eds.
J. Bu'lock and B. Kristiansen. London, England: Academic Press,
1987, pp. 197−215.
8-28 Sterilization

Strauss, W., Industrial Gas Cleaning (2nd ed.), pp. 182, 278−297. Oxford,
England: Pergamon Press Ltd., 1975.
Wehner, J. F. and R. H. Wilhelm, “Boundary Conditions of Flow
Reactor,” Chem. Eng. Sci. 6 (1956):89−93.