8 Scaleup
8 Scaleup
8 Scaleup
CH-528
Biochemical Engineering
Scale up of bioreactors
By: Prof. Prasenjit Mondal
Department of Chemical Engineering
Contents
2
8.1 Scale-up and Dimensionless
numbers for Scale-up
For the optimum design of a production-scale fermentation system
(prototype), we must translate the data on a small scale (model) to the large
scale. The fundamental requirement for scale-up is that the model and
prototype should be similar to each other.
Two kinds of conditions must be satisfied to insure similarity between
model and prototype. They are:
1. Geometric similarity of the physical boundaries: The model and the
prototype must be the same shape, and all linear dimensions of the
model must be related to the corresponding dimensions of the prototype
by a constant scale factor.
2. Dynamic similarity of the flow fields: The ratio of flow velocities of
corresponding fluid particles is the same in model and prototype as well
as the ratio of all forces acting on corresponding fluid particles. When
dynamic similarity of two flow fields with geometrically similar
boundaries is achieved, the flow fields exhibit geometrically similar flow
patterns.
3
Continued…
The first requirement is obvious and easy to accomplish, but the second is
difficult to understand and also to accomplish and needs explanation. For
example, if forces that may act on a fluid element in a fermenter during
agitation are the viscosity force , drag force on impeller , and gravity force ,
each can be expressed with characteristic quantities associated with the
agitating system. According to Newton's equation of viscosity, viscosity
force is
𝐹 𝑉=𝜇 ( 𝑑𝑢
𝑑𝑦 ) 𝐴 8.1
where is velocity gradient and A is the area on which the viscosity force
acts. For the agitating system, the fluid dynamics involved are too complex
to calculate a wide range of velocity gradients present. However, it can be
assumed that the average velocity gradient is proportional to agitation
speed N and the area A is to which results.
4
Continued…
2
𝐹 𝑉 ∝𝜇 𝑁 𝐷𝐼 8.2
The drag force can be characterized in an agitating system as
𝑃 𝑚𝑜
𝐹𝐷 ∝ 8.3
𝐷𝐼 𝑁
Since gravity force is equal to mass m times gravity constant g,
2
𝐹 𝐺∝ 𝜌 𝐷𝐼 𝑔 8.4
The summation of all forces is equal to the inertial force as,
∫ 𝐹=𝐹𝑉 +𝐹 𝐷 +𝐹 𝐺=𝐹 𝐼 ∝𝜌 𝐷 𝑁 4
𝐼
2
Then dynamic similarity between a model (m) and a prototype (p) is
8.5
achieved if
( 𝐹 𝑉 )𝑚 ( 𝐹 𝐷 )𝑚 ( 𝐹 𝐺)𝑚 ( 𝐹 𝐼 )𝑚
= = = 8.6
( 𝐹 𝑉 )𝑝 ( 𝐹 𝐷) 𝑝 ( 𝐹 𝐺)𝑝 ( 𝐹 𝐼 )𝑝
5
Continued…
or in dimensionless forms:
( ) ( )( ) ( )( ) ( )
𝐹𝐼
𝐹𝑉 𝑝
=
𝐹𝐼
𝐹𝑉 𝑚
𝐹𝐼
𝐹𝐷 𝑝
=
𝐹𝐼
𝐹𝐷 𝑚
𝐹𝐼
𝐹𝐺 𝑝
=
𝐹𝐼
𝐹𝐺 𝑚
8.7
6
Continued…
( 𝑁 ¿ ¿ 𝐹𝑟 ) 𝑃 =( 𝑁 𝐹𝑟 )𝑚 ¿
7
Example 8.1:
( )
2 2
𝑃 𝜌 𝑁𝐷 𝐼 𝑁 𝐷 𝐼
=𝑓 , 8.12
3
𝜌 𝑁 𝐷𝐼
5
𝜇 𝑔
Can you determine the power consumption and impeller speed of a 10,000-
gallon fermenter based on the findings of the optimum condition from a
geometrically similar one-gallon vessel? If you cannot, can you scale up by
using a different fluid system?
Solution:
Since , the scale ratio is,
( 𝐷 𝐼 )𝑝 1/ 3
=1000 =10 8.13
( 𝐷 𝐼 )𝑚
8
Continued…
( 𝑃 𝑚𝑜
𝜌 𝑁
3
𝐷 ) (
5
𝐼 𝑝
=
𝑃 𝑚𝑜
𝜌 𝑁
3
𝐷 )
5
𝐼 𝑚
8.14
( 𝜌 𝑁𝐷2
𝜇
𝐼
) ( 𝑝
=
𝜌 𝑁𝐷 2
𝜇
𝐼
) 𝑚
8.15
( 𝑁2 𝐷𝐼
𝑔 ) ( 𝑝
=
𝑁 2 𝐷𝐼
𝑔 ) 𝑚
8.16
If you use the same fluid for the model and the prototype, and . Canceling
out the same physical properties and substituting Eq. (8.14) to Eq. (8.16)
yields,
9
Continued…
( )
3
5 𝑁𝑝
( 𝑃 𝑚𝑜 )𝑝 =10 ( 𝑃𝑚𝑜 )𝑚 8.17
𝑁𝑚
The equality of the Reynolds number requires
𝑁 𝑃 =0.01 𝑁 𝑚 8.18
which is conflicting with the previous requirement for the equality of the
Reynolds number. Therefore, it is impossible to satisfy the requirement of
the dynamic similarity unless you use different fluid systems.
10
Continued…
If and , to satisfy Eqs. (8.15) and (8.16), the following relationship must
hold.
( 𝜇
𝜌 ) 𝑚
=
1
31.6 ( 𝜇
𝜌 ) 𝑝
8.20
11
Continued…
( 𝑃 𝑚𝑜
𝜌 𝑁
3
𝐷
5
𝐼
) (
𝑝
=
𝑃 𝑚𝑜
𝜌 𝑁
3
𝐷
5
𝐼
)
𝑚
8.21
If the fluid employed for the prototype and the model remains the same,
the power consumption in the prototype is
)[ ]
5
(
3
𝑁𝑝 ( 𝐷 𝐼 )𝑚
( 𝑃 𝑚𝑜 )𝑝 =( 𝑃 𝑚𝑜)𝑚 8.22
𝑁𝑚 ( 𝐷𝐼)𝑝
12
8.2 Scaling of Mass-Transfer Equipment
As discussed by Oldshue, the various quantities which may influence the
product in an agitated industrial reactor do not scale in the same way
with reactor size or impeller rate.
1. The turbulent Reynolds number , determines and thus bubble mass
transfer coefficients
𝜌 𝑢𝑟𝑚𝑠 𝐷 𝜌 𝐷 𝐷 𝑃
( )
1/ 3
Re 𝑡 = ∝ 8.23
𝜇 𝜇 𝜌 𝑉
2. The impeller tip velocity , determines the maximum shear rate , which
in turn influences both maximum stable bubble or microbial floc size
and damage to viable cells.
3. The power input per unit volumethrough , determines mass-transfer
coefficients and particulate sizes. In laminar and transition regimes of
aerators
8.24
13
Continued…
Figure 8.1: Different relationships among air bubbles (B), n-hexadecane droplets
(O), and yeast cells (C), at different stages of batch culture of petrophilum (1 = log
phase, 2 = first half phase, 3 = second half of exponential phase, 4 = after n-
hexadecane exhaustion).
14
Continued…
[ ]
0.45
3 0.44
𝑃𝑎 𝑃 (𝑁𝑖 𝐷 )
2
′ 𝑖 8.28
≈𝑚
𝑉 𝑉 2.22 𝑁 0.56
𝑎
15
Continued…
16
Continued…
where 1 and 2 denote the values in the small and large scale vessel,
respectively.
As indicated in Fig. 8.2a this gives very similar yields of penicillin for
vessels from 5 liters to 200 gallons. However, from the curves in Fig. 8.2a,
we see that, at different P/V values, there are significant differences in yield
at different scales.
Another frequently applied basis for scale-up is constant volumetric
transfer coefficient . Fig. 8.2b shows vitamin yields from bacterial
fermentations at different scales versus corresponding values of ().
The total pressure is included here to correct for the greater driving force
for oxygen transfer at higher pressures which are encountered in large-scale
bioreactors.
17
Continued…
Figure 8.2: (a) Penicillin yields vs. power input in various sized vessels
(After E Gaden, Super Sanita, vol. 1, p. 61, 1961). (b) Vitamin B12 yield
(ug/g) vs. mass-transfer group (
18
Continued…
19
Continued…
When the bubbles rise in plug flow through the vessel but the impeller
still maintains perfect mixing in the liquid phase, varies with position. Over
a differential reactor height dz, the instantaneous loss of oxygen from the
bubble is
𝑑𝑝 𝑂 1
𝐻 𝐴 𝑑𝑧 2
=
𝑔𝑎𝑠 𝑣𝑜𝑙
(
𝑑𝑖𝑓𝑓𝑟𝑒𝑛𝑡𝑖𝑎𝑙
)(
𝑐𝑜𝑛𝑐 . 𝑟𝑎𝑡𝑒 𝑜𝑓
𝑑𝑡 𝑅𝑇 𝑟𝑒𝑎𝑐𝑡𝑜𝑟 𝑣𝑜𝑙 𝑟𝑒𝑎𝑐𝑡𝑜𝑟 𝑣𝑜𝑙𝑢𝑚𝑒 𝑐h𝑎𝑛𝑔𝑒𝑖𝑛𝑏𝑢𝑏𝑏𝑙𝑒 ) 8.32
which equals
− k l 𝑎(𝑐 ¿ ¿ ∗− 𝑐 𝑏) 𝐴𝑑𝑧 ¿
the mass transfer rate into the liquid. Since (M is Henry’s law constant), we
have
∗
𝐻 𝑑𝑝 𝑂 𝐻𝑀 𝑑 𝑐 𝑙 ∗
= = − 𝑎 k l (𝑐 𝑙 − 𝑐 𝑏)
2
8.33
𝑅𝑇 𝑑𝑡 𝑅𝑇 𝑑𝑡
20
Continued…
(
(𝑐 ∗𝑙 − 𝑐𝑏 ) 𝑧 =( 𝑐∗𝑙 − 𝑐 𝑏)𝑖𝑛𝑙𝑒𝑡 exp
The overall mass-transfer rate in the volume Ah is therefore
)− 𝑎 k l 𝑅𝑇 𝑧
𝐻𝑀 𝑢𝑏
8.35
[ ( )]
h
𝐻𝑀 𝑢 𝑏 ∗ − 𝑎 k l 𝑅𝑇 h
∫ 𝑎 k l ( 𝑐 𝑙 −𝑐 𝑏 ) ( 𝑧 ) 𝐴 𝑑𝑧 =
∗
𝑅𝑇
( 𝑐 𝑙 −𝑐 𝑏 )𝑖𝑛𝑙𝑒𝑡 𝐴 1 − exp
𝐻𝑀 𝑢𝑏
0
8.36
21
8.3 Criteria for Scaling-Up Fermentors
Few steps are available for the so-called “scale-up” of fermentors,
starting from the glass apparatus used in fundamental research
investigations. Usually, chemical engineers are responsible for building and
testing pilot-scale fermentors (usually with capacities of ), and subsequently
for designing industrial-scale fermentors based on data acquired from the
pilot-scale system.
One criterion for scaling-up this type of bioreactor is the power input per
unit liquid volume of geometrically similar vessels, which should be
proportional to for the turbulent range and to for the laminar range, where
is the rotational stirrer speed and L is the representative length of the
vessel.
In order to minimize any physical damage to the cells, the product of the
diameter and the rotational speed of the impeller - which should be
proportional to the tip speed of the impeller and hence to the shear rate at
the impeller tip - becomes an important criterion for scale-up.
22
Aerated stirred tanks, bubble columns, and airlifts are usually used for
aerobic fermentations. One criterion of scaling-up aerated stirred tank
fermentor is , approximate values of which can be estimated. For the
turbulent range, a general correlation for in aerated stirred fermentors is of
the following type:
𝑚 3 2 𝑛
𝑘 𝐿 𝑎=𝑐 𝑈 𝐺 (𝑁 𝐿 ) 8.37
where is the superficial gas velocity over the cross-sectional area of the
tank, is the rotational stirrer speed, and is the representative length of
geometrically similar stirred tanks.
For the turbulent regime () should be proportional to the power requirement
per unit liquid volume.
23
Continued…
24
Continued…
Example 8.2: Two geometrically similar stirred tanks with flat-blade turbine
impellers of the following dimensions are to be operated at 30°C as pilot-
scale and production-scale aerobic fermentors.
25
Continued…
Using as the scale-up criterion, estimate the impeller speed and the air rate
for the production-scale fermentor that will give results comparable with
the pilot-scale data. 𝑚 3 2 𝑛
𝑘 𝐿 𝑎=𝑐 𝑈 𝐺 ( 𝑁 𝐿 )
Solution:
In the pilot-scale fermentor:.
At equal , the air rate to the production-scale fermentor should be
0.0295 × 60 × ( )
𝜋
4
×(2.0)2=5.56 𝑚3 𝑚𝑖𝑛−1
With the pilot fermentor .
With this equal value of and for the production fermentor, the production
fermentor should be operated at
1.215
3
𝑁 = = 0.304 𝑁 =0.672 𝑠 −1
4
26
Continued…
Incidentally, the impeller tip speeds in the pilot and production fementors
are calculated as and , respectively.
27
8.4 Design estimation of various Scale-
up Parameters
Most often, power consumption per unit volume is employed as a criterion
for scale-up. In this case, to satisfy the equality of power numbers of a
model and a prototype,
) ( )[ ]
2
( ) (
3
𝑃 𝑚0 𝑃 𝑚0 𝑁𝑝 ( 𝐷𝐼)𝑝
= 8.39
𝐷
3
𝐼 𝑝 𝐷
3
𝐼 𝑚
𝑁𝑚 ( 𝐷 𝐼 )𝑚
Note that represents the power per volume because the liquid volume is
proportional to for the geometrically similar vessels. For the constant ,
[ ]
2
( )
3
𝑁𝑝 ( 𝐷 𝐼 )𝑚 8.40
=
𝑁𝑚 ( 𝐷 𝐼 )𝑝
28
Continued…
[ ]
2 /3
( 𝐷 𝐼 )𝑚
𝑁 𝑝= 𝑁𝑚
( 𝐷𝐼 )𝑝
𝑁 𝑝 =0.34 𝑁 𝑚 8.41
which shows that the impeller speed in a prototype vessel is about one
third of that in a model.
For constant , the Reynolds number and the impeller tip speed cannot
be the same. For the scale ratio of 5,
( 𝑁 𝑅𝑒 )𝑝 =8.5 ( 𝑁 𝑅𝑒 )𝑚
𝑖 𝑖
8.42
( 𝑁 𝐷 𝐼 )𝑝 =1.7 ( 𝑁 𝐷 𝐼 )𝑚 8.43
29
Continued…
Table 8.1 shows the values of properties for a prototype (2,500 gallon)
when those for a model (20-gallon) are arbitrarily set as 1.0 (Oldshue,
1966).
The parameter values of the prototype depend on the criteria used for
the scale-up. The third column shows the parameter values of the
prototype, when is set constant. The values in the third column seem to be
more reasonable than those in the fourth, fifth, and sixth columns, which
are calculated based on the constant value of , and , respectively.
For example, when the Reynolds number is set constant for the two
scales, the reduces to 0.16 percent of the model and actual power
consumption also reduces to 20 percent of the model, which is totally
unreasonable.
30
Continued…
31
Continued…
33
8.5 Power Estimation for Ungassed Liquids
34
Continued…
Figure 8.6: Correlation between Reynolds number (Re) and Power number ().
Curve (a): six-flat blade turbine, four baffles,
curve (b): two-flat blade paddle, four baffles, =0.1 D,
curve (c): three-blade marine propeller, four baffles, =0.1 D.
35
Continued…
In Figure 8.6, the curves a, b, and c correlate data for three types of
impellers, namely, the six-flat blade turbine, two-flat blade paddle, and
three-blade marine propeller, respectively. It should be noted that, for the
range of (Re) >, is independent of (Re). For this turbulent regime it is clear
from Equation 8.44 that
3 5
𝑃 = 𝑐1 𝜌 𝑁 𝑑 8.46
where , is a constant that varies with the impeller types. Thus, for a given
type of impeller varies in proportion to , and liquid density ,but is
independent of liquid viscosity .
For the ranges of Re below approximately 10, the plots are straight lines
with a slope of ; that is, , is inversely proportional to . Then, for this laminar
regime, we can obtain Equation 7.33 from Equations 8.44 and 8.45:
𝑃 = 𝑐2 𝜇 𝑁
2
𝑑
3 8.47
36
Continued…
Thus, for the laminar regime varies in proportion to liquid viscosity , , and
to a constant , which varies with impeller types, although is independent of
the liquid density .
It is worth remembering that the power requirements of geometrically
similar stirred tanks are proportional to in the turbulent regime and to in
the laminar regime.
In the case of scaling-up of geometrically similar stirred tanks, the equal
power input per unit liquid volume, which should be proportional to for the
turbulent regime and to for the laminar regime, is sometimes used as a
criterion for scale-up, because it leads to the same values of in a larger
stirred tank with that in the smaller stirred tank.
37
8.6 Power Estimation for Gas-Sparged Liquids
( ) ( ) ( ) ( )
0.115 2 1.96 (𝑑/ 𝐷)
( )
𝑃𝐺 𝑑
4.38 2
𝑑 𝑁 𝑑𝑁 𝑄
log =−192 8.48
𝑃0 𝐷 𝑣 𝑔 𝑁𝑑
3
where is the impeller diameter (L), is the tank diameter (L), is the
rotational speed of the impeller (), g is the gravitational constant (), Q is the
gas rate (), and is the kinematic viscosity of the liquid ().
38
Continued…
39
Continued…
40
Continued…
41
Continued…
928 𝑔+ 2 ( 𝑔 1
)
. .20000 𝑙 =20928 𝑔
𝑙 2
If no wall growth had been present, the tank would have yielded . Thus,
wall growth, if present, can seriously alter the productivity of a large-scale
system upon scale-up.
42
Continued…
Solution:
Assume the vessel is a cylinder.
Thus,
𝑉 = ( 𝜋
4 ) 𝐷 𝑡2 𝐻
43
Continued…
but , so
3 3
𝑉 =( 𝜋 / 4 )3 𝐷 𝑡 =10,000 𝑐𝑚
Solving for , , and gives:
The scale-up factor is the cube root of the ratio of tank volumes, or
( )
1/ 3
10000 𝑙
=10
10𝑙
To maintain geometric similarity the larger vessel will have its dimensions
increased by a factor of 10. That is:
44
Continued…
a. For constant to apply must be the same in both vessels. Let subscript
refer to the small vessel and to the large vessel.
( )
2/ 3
( )
2 /3
𝐷𝑡1 1
𝑁 2=𝑁 1 =500 𝑟𝑝𝑚 =107 𝑟𝑝𝑚
𝐷𝑡2 10
b. For constant tip speed to apply must be the same in both vessels.
𝑁 2= 𝑁 1
( )
𝐷𝑡1
𝐷𝑡2
=500 𝑟𝑝𝑚
1
10 ( )
=50 𝑟𝑝𝑚
( )
2
𝐷𝑡1
( ) =5 𝑟𝑝𝑚
2
1
𝑁 2= 𝑁 1 =500 𝑟𝑝𝑚
𝐷𝑡2 10
Scale-up on the basis of constant would be the most likely choice unless
the culture was unusually shear sensitive. Scale-up based on constant Re
would not be used.
45