Lessons in Bioreactor Scal-Up P2 a Refresher on Fluid Flow and Mixing
Lessons in Bioreactor Scal-Up P2 a Refresher on Fluid Flow and Mixing
LY
Muhammad Arshad Chaudhry
N
C
O
ommercial-scale production of components. Such force must be applied Figure 1: Familiar examples of laminar
recombinant therapeutic proteins to produce fluid flow (1). Scientists and turbulent flow
routinely involves suspension categorize fluids according to their
N
Turbulent flow
cultures of mammalian cells in density and viscosity as
bioreactors with up to 10,000 L capacity. • compressible (or incompressible) —
IO
With advances made over the past 30 whether their density changes (or not)
years in cellular engineering, basal and in response to pressure, with gases
feed media development, and bioprocess generally being compressible and
S
engineering, expression titers of ~10 g/L, liquids being incompressible
viable-cell densities of >3 × 107 cells/mL, • viscous (all fluids have a finite
IS
and specific productivities of >20 pg/cell/ viscosity, which is responsible for Laminar flow
day are now common. Such high cell “internal friction” during flow,
M
densities (with the potential to go even otherwise expressed as a fluid’s
higher) increase mixing and aeration resistance to motion)
demands and can subject cells to • Newtonian or non-Newtonian, those layers, and that influences the
R
aggressive environments such that they depending on whether they obey velocity of their flow. Particles moving
experience high hydrodynamic stresses. Newton’s law of viscosity. from a slow-moving layer to a faster
E
This installment of my bioreactor Shear stress is calculated by dividing layer will act to reduce its speed; those
scale-up series should refresh your the shear force applied to a fluid by the coming from a fast layer to a slower one
P
knowledge of fluid flow, mixing, and cross-sectional area to which it is will have an accelerating effect. Thus
mass transfer in bioreactors, including applied. Strictly speaking, shear stress are shear stresses induced by fluid flow
H
how the interplay of those parameters is one of two types of fundamental as a consequence of velocity differences.
creates the environment that cells hydrodynamic forces, the other being Depending on fluid velocity, flow can
IT
experience. That in turn influences their “normal stress.” Both are determined be characterized as laminar or turbulent
growth, metabolism, and protein commonly as viscosity multiplied by the (Figure 1). When the overall fluid motion
production. velocity gradient, with flow direction is slow, adjacent fluid layers move in an
W
being perpendicular for shear stress and orderly (more or less parallel) motion
Fluid Flow in Bioreactors parallel for normal stress. Hence, a called laminar flow. In a fast-moving
Fluids — either liquids or gases — are velocity gradient for shear stress is fluid, particles and parcels of fluid
T
substances that undergo continuous written as dUx/dy and for normal stress frequently cross layers randomly and
deformation when subjected to shear as dUx/dx, where U is the velocity of chaotically in a turbulent flow.
IN
force. Shearing causes deformation, a fluid and x and y represent the direction Reynolds number (Re) is a
change in the relative positions of fluid of the force/stress. dimensionless variable for characterizing
R
and manufacturing solid object (zero velocity). The next and viscosity of the fluid, respectively.
fluid layer inward has a bit higher Because the biopharmaceutical
KEYWORDS: Impellers, Reynolds
velocity because of the overall fluid industry is concerned primarily with
number, power number, shear,
motion, the next layer inward has more, fluid flow in stirred-tank bioreactors,
distribution, dispersion, and diffusion
and so on. Within every fluid is a the relevant Reynolds number is defined
LEVEL: Basic continuous exchange of particles among by Equation 2:
H = tank height
* Maximum values during a process at the current bioreactor size ** Calculated from EDR mean
HL = fluid height
Re i = Ni Di2ρ/μ direction, but superimposed on that are Di = impeller diameter
where Ni is the stirrer speed, Di is the secondary chaotic movements of fluid
impeller diameter, and ρ and μ are particles causing flow deviations and
C = impeller
clearance
density and viscosity of the fluid, generating vortices of varying shape,
Sparger
respectively. size, speed, and rotational direction. disc
Conceptually, the Reynolds number Those irregular secondary motions in
expresses the ratio of inertial to viscous turbulent flow have profound
forces. During laminar flow (low consequences for mixing efficiency,
Reynolds number), viscous forces energy loss, and shear intensity. (NRei > 104), the size of the smallest
dominate, and fluid motion is slow. Turbulent flows develop spinning or eddies is given approximately as the
However, with increasing Reynolds swirling fluid structures called eddies, Kolmogorov scale of mixing (λ) or scale
numbers, inertial forces become more which can stretch, coalesce, and divide. of turbulence, defined by Equation 3:
dominant, and flow transitions from Fluid velocity within an eddy undergoes λ = (ν3/ε)1/4
laminar to turbulent — then becomes great change in magnitude and where λ is the characteristic length
fully turbulent at sufficiently high Re direction within a relatively short of the smallest eddies, ν is the
values. distance and time. kinematic viscosity of the fluid, and ε is
Pipe flow is laminar if Re < 2100. For Eddies of different sizes appear in the local rate of turbulent energy
Re values of 2000–4000, the flow is turbulent flow. The size of the largest dissipation per unit mass of the fluid.
considered to be in a transition state eddies is limited by the boundary of the At a steady state, the average rate of
from laminar to turbulent. It is fully flow system. Hence, the diameter of the turbulent energy dissipation throughout
turbulent once Re > 4000. In stirred largest eddies in a bioreactor is similar a bioreactor is equal to the power input
tanks, the value of Re that marks the to (but somewhat smaller than) the by the impeller to the fluid within the
transition from laminar to turbulent flow diameter of the reactor vessel. Large tank. Hence, the greater the power input
depends on tank and impeller geometry. eddies are unstable and give rise to to the impeller, the smaller the eddies
smaller ones, which in turn produce will be. Because λ depends on viscosity,
Turbulent Flow in Bioreactors even smaller eddies, and so on. In a smaller eddies are produced in low-
Turbulent flow is the flow regime that stirred bioreactor, large eddies derive viscosity media for a given power input.
cells encounter in bioreactors. Turbulence energy from the bulk flow generated by An energy input of 0.1 W/kg, which
is essential for effective mixing (mass the impeller, so they contain most of the currently is considered to be high for
and heat transfer in fluids), so achieving turbulent kinetic energy. The smallest animal-cell cultures, yields an estimate
turbulent-flow conditions in bioreactors is eddies cannot sustain rotational motion for the smallest eddy diameter λ of about
vital to culture success. As noted above, and lose their energy as heat through 0.06 mm (60 μm). That represents the
the Reynolds number indicates when the effects of viscosity and fluid friction. smallest scale of mixing achievable by
turbulent flow is generated. As that Maintaining turbulent flow thus requires dispersion under such conditions.
number increases, inertial forces come to continuous supply of energy to make up The characteristic length of the
dominate viscous forces in the fluid, for what is lost as heat (1). smallest eddies (λ) often is an indicator
overcoming the tendency of viscous Dispersion is the process of breaking of the potential for cell damage caused
effects to dampen flow instability. Thus, up the bulk flow into smaller and smaller by fluid shear. Experiments have shown
turbulence can be regarded as a highly eddies to facilitate rapid transfer of that if biological entities (e.g.,
disordered fluid motion resulting from material throughout a tank. The degree mammalian cells) are smaller than λ in a
instabilities growing in an initially of homogeneity possible as a result of bioreactor, then shear damage to such
laminar flow field (1). dispersion is limited by the size of the entities will not occur. Minow et al.
Turbulent flow is very complex by smallest eddies that form. Assuming that presented a comparative study of mixing
nature, so our understanding of it is far isotropic turbulence is prevailing in a and mass transfer in two different
from complete. The flow has an overall stirred tank and that flow is turbulent 1000‑L single-use bioreactors. The
Θ95% (seconds)
edges of each flat impeller blade (Figure
5). The vortices play a critical role to
determine mixing performance because
most mixing takes place near the
vortices issuing from the impeller blades. H/T = 1.60
These trailing vortices also provide gas D/T = 0.33
Θ95% = 10.52/3 · ԑ̄
ε –1/3 · (D/T )–1/3
dispersion in sparged systems. The
associated steep velocity gradients 10
10–2 10–1 1
provide a major source of turbulence. ε̄95%(W/kg)
One reason for widespread use of
Rushton turbines in microbial fermentors
is their ability to handle high gas-flow pitched-blade impeller is about 1. Thus, and diffusion. Distribution is also called
rates and provide efficient gas dispersal. a Rushton turbine has strong form drag, “macromixing,” diffusion is known as
The power required to achieve a given generates high levels of torque, and “micromixing,” and depending on the
stirrer speed depends on frictional forces transmits more power at the same scale of fluid motion, dispersion can be
and form drag (caused by the pressure operating speed than does a pitched classified as either macro- or
difference between the front and rear of blade impeller (1). micromixing.
the impeller blades) that resists impeller Fluid Flow: Rotating impeller blades Distribution is the process by which
rotation. Friction and form drag create pump fluid in a bioreactor vessel. The fluid is transported to all regions of a
torque on the stirrer shaft. The power volumetric flow rate of fluid leaving a stirred vessel by bulk circulation
input can be determined experimentally blade is a characteristic of the impeller currents. It is an important process in
by measuring the induced torque design and varies with operating mixing and often the slowest step
(Equation 4): parameters such as stirrer speed and because fluid has a great distance to
P = 2πNi M blade size. A dimensionless flow number traverse in a long circulation path.
where P is the power, Ni is the stirrer represents the effectiveness of a Dispersion breaks up the bulk flow into
speed, and M is the torque. pumping process (Equation 6): smaller and smaller eddies, which
The power required to mix nongassed Nflow = Q/Ni Di3 facilitates rapid transfer of material
liquids depends on stirrer speed, impeller where Q is the volumetric flow rate of through a vessel. The degree of
shape and size, tank geometry, and liquid fluid leaving the impeller blade, Ni is homogeneity possible from dispersion is
density and viscosity. The relationships the stirrer speed, and Di is the impeller limited by the size of the smallest
among those variables generally is diameter. eddies formed in a particular fluid. To
represented by the dimensionless power A pitched-blade impeller produces achieve mixing on a scale smaller than
number (Equation 5): nearly three to four times the flow of a the smallest eddies, diffusion becomes
NP = P/ρNi3Di5 Rushton turbine for the same power an important process. Such micromixing
where NP is the power number, P is input. This is why the former design is occurs relatively rapidly and over
the power, Ni is the stirrer speed, and Di recognized for high pumping capacity shorter distances than are involved in
is the impeller diameter. and effectiveness of blending operations, distribution.
Consider the power number as whereas the latter are considered to have Mixing time (tm) is a useful parameter
analogous to a drag coefficient for a low pumping efficiency (1). for assessing the overall speed of mixing
stirrer system. The relationship between For mixing to be effective in stirred in stirred vessels. It is the time required
NP and Re generally is determined vessels, the velocity of fluid leaving an to achieve a given degree of homogeneity
experimentally for a range of impeller impeller must be sufficient to carry after starting with completely segregated
and tank configurations. The power material into the most remote regions of materials. Mixing time can be measured
number for different impeller designs a tank. Circulating fluid also must by injecting tracer material — commonly
becomes constant once flow is turbulent sweep the entire vessel in a reasonable acids, bases, or concentrated salt
(1). Under turbulent flow conditions in a time. Fluid flow must be turbulent for solutions — into a vessel and following its
stirred tank, the power number for a good mixing, which can be described as concentration at a fixed point in the tank.
Rushton turbine generally remains a combination of three physical When a small pulse of tracer is added
constant at 5, whereas that of an axial processes (1): distribution, dispersion, to fluid in a stirred tank that already