Skip to content

Mass Transfer

Mass transfer by diffusion

Mass transfer inside a bioreactor happens through diffusion.

Diffusion

Diffusion is the thermal motion of fluid particles at any temperature above absolute zero — 0 K, or −273 °C.

Fick's law of diffusion

Fick's law of diffusion states that molecules move from areas of higher concentration to areas of lower concentration.

Dense cluster of component A particles at higher concentration on the left, with arrows passing through a circular cross-sectional area toward a sparser cluster of the same particles at lower concentration on the right

Mass transfer of component A across a cross-sectional area perpendicular to the direction of diffusion can be explained by this law.

JA=NAa=−DABdCAdyJ_A = \frac{N_A}{a} = -D_{AB}\frac{dC_A}{dy}
  • JAJ_A — mass flux of component A (rate of diffusive mass transport per unit cross-sectional area, driven by the concentration gradient), in kilograms per square metre-second (kg/(m2⋅s)\text{kg/(m}^2 \cdot \text{s)})
  • NAN_A — rate of mass transfer of component A, in kilograms per second (kg/s\text{kg/s})
  • aa — cross-sectional area, perpendicular to the direction of diffusion, in square metres (m2\text{m}^2) — not to be confused with the subscript AA, which denotes the diffusing component
  • DABD_{AB} — diffusion coefficient of A through B, describing the ease of diffusion, in square metres per second (m2/s\text{m}^2/\text{s})
  • CAC_A — mass concentration of component A, in kilograms per cubic metre (kg/m3\text{kg/m}^3)
  • yy — distance in the direction of diffusion, in metres (m\text{m})

The negative sign in the equation is because diffusion occurs from higher concentration to lower concentration.

Film theory

Film theory states that at the interface between two phases, there exists a thin, stagnant film or layer of fluid. Within this film, mass transfer occurs only via molecular diffusion.

For example, when a bubble enters a bioreactor, a thin film of the same bulk fluid surrounds the bubble. This film makes it difficult for gas inside the bubble to diffuse out of the bubble and into the fluid.

Concentration profile across the film

Concentration profile of solute A across two phases in contact — organic phase bulk fluid on the left, aqueous phase bulk fluid on the right, separated by a phase boundary flanked by an organic film and an aqueous film; a curve shows concentration falling from CA1 in the aqueous bulk fluid, through CA1i and CA2i at the phase boundary, down to CA2 in the organic bulk fluid

A solute moves from the bulk fluid on one side, across the phase boundary, into the bulk fluid on the other side:

  1. In the bulk fluid of the first phase (e.g. the aqueous phase), the solute is at concentration CA1C_{A1}.
  2. Moving through that phase's film toward the phase boundary, the concentration drops to CA1iC_{A1i} (the interface concentration on this side) — the film poses resistance to transfer, since diffusion is the only way for the solute to move through it.
  3. At the phase boundary itself, there is no transfer resistance — instantaneous thermodynamic equilibrium is established. The concentration can rise or drop here, from CA1iC_{A1i} to CA2iC_{A2i}, because the solubility of the solute differs between the two phases — the partition coefficient decides the concentration on each side.
  4. Moving across the second phase's film (e.g. the organic phase), the concentration continues to drop, from CA2iC_{A2i} down to CA2C_{A2} in that phase's bulk fluid.

Overall: CA1>CA1i>CA2i>CA2C_{A1} > C_{A1i} > C_{A2i} > C_{A2}.

Within each film, the concentration profile is linear. Film theory assumes steady-state molecular diffusion through the film with no accumulation, and integrating Fick's law under that assumption gives a straight-line concentration gradient — the same reasoning that gives a linear temperature profile across a wall in steady-state heat conduction.

Rate of mass transfer (film theory)

Based on film theory, the mass transfer rate is:

NA=kL a (C∗−C)N_A = k_L\, a\, (C^* - C)
  • NAN_A — volumetric rate of mass transfer (how fast solute crosses the interface into the bulk liquid, driven by the gap between saturation and bulk concentration), in gram-mole per cubic metre per second (gmol/(m3⋅s)\text{gmol}/(\text{m}^3 \cdot \text{s}))
  • kLk_L — (liquid) film mass-transfer coefficient (how fast a solute molecule crosses the stagnant film — sets the pace of interfacial transfer), in metres per second (m/s\text{m/s})
  • aa — interfacial area per unit volume, in square metres per cubic metre (m2/m3\text{m}^2/\text{m}^3)
  • C∗C^* — saturation concentration: the maximum concentration of the solute that can dissolve in the liquid, in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)
  • CC — bulk concentration: the current concentration of the solute in the liquid, in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)

kLk_L is how fast the molecule crosses the film: the thinner the film, the higher kLk_L, and the faster the transfer.

C∗−CC^* - C is the driving force for mass transfer: the bigger the difference between the saturation concentration and the bulk concentration, the more strongly the system pulls the solute into the liquid.

kLk_L under film theory

Film theory treats the film as a fixed-thickness layer crossed by steady-state molecular diffusion, which gives kLk_L itself as

kL=DABδk_L = \frac{D_{AB}}{\delta}
  • kLk_L — (liquid) film mass-transfer coefficient (how fast a solute molecule crosses the stagnant film), in metres per second (m/s\text{m/s})
  • DABD_{AB} — diffusion coefficient of the solute through the liquid, in square metres per second (m2/s\text{m}^2/\text{s})
  • δ\delta — thickness of the stagnant film, in metres (m\text{m})

kLk_L is directly (linearly) proportional to DABD_{AB} — doubling the diffusivity, with film thickness δ\delta unchanged, doubles kLk_L.

Worked example: effect of doubling diffusivity on kLk_L

Question. If the molecular diffusivity of a solute doubles while keeping all other physical parameters and film thickness constant, according to Whitman's film theory, the mass transfer coefficient kLk_L will:

  • (a) Increase by a factor of 1.411.41 (2\sqrt{2})
  • (b) Increase by a factor of 2.02.0
  • (c) Increase by a factor of 4.04.0
  • (d) Remain unchanged

Solution.

  1. Film theory gives kL=DAB/δk_L = D_{AB}/\delta — a linear relationship between kLk_L and DABD_{AB}, with δ\delta held fixed.
  2. Doubling DABD_{AB} (with δ\delta constant) doubles kLk_L directly: kL,new=(2DAB)/δ=2 (DAB/δ)=2 kL,oldk_{L,\text{new}} = (2D_{AB})/\delta = 2\,(D_{AB}/\delta) = 2\,k_{L,\text{old}}.
  3. A factor of 2\sqrt{2} (option a) is what penetration or surface renewal theory would predict instead, since those give kL∝DAB0.5k_L \propto D_{AB}^{0.5} — not film theory.

Answer: (b) — kLk_L increases by a factor of 2.02.0.

Interfacial area per unit volume (aa)

Interfacial area per unit volume (aa) is the total interfacial area between the two phases, divided by the total volume of the liquid (the bulk liquid) in the vessel — not the total vessel volume including the gas phase. In practice, the gas phase (bubbles) occupies only a small fraction of the vessel, so the liquid volume is used as the "total volume" in this ratio.

For a sphere, a=S/V=6/da = S/V = 6/d — so smaller bubbles mean higher aa. Splitting a fixed total gas volume into nn equal bubbles multiplies aa (and the total interfacial area) by n1/3n^{1/3}. This is why sparging aims to break gas into fine bubbles rather than large ones.

Worked example: interfacial area from bubble size

Question. A sparger releases a fixed volume of gas as a single bubble of diameter 10 mm10\ \text{mm}. What is the interfacial area per unit volume (aa) for (a) that single bubble, and (b) the same total gas volume split instead into 100 equal smaller bubbles?

Solution.

  1. For a sphere, a=6/da = 6/d.
  2. Single bubble: d=10 mm=0.01 md = 10\ \text{mm} = 0.01\ \text{m}, so a=6/0.01=600 m2/m3a = 6/0.01 = 600\ \text{m}^2/\text{m}^3.
  3. Splitting a fixed total gas volume into nn equal bubbles scales aa by n1/3n^{1/3} (smaller diameter, more total surface for the same volume).
  4. For n=100n = 100: a=600×1001/3≈600×4.64≈2790 m2/m3a = 600 \times 100^{1/3} \approx 600 \times 4.64 \approx 2790\ \text{m}^2/\text{m}^3.

Answer: a≈600 m2/m3a \approx 600\ \text{m}^2/\text{m}^3 for the single bubble, rising to ≈2790 m2/m3\approx 2790\ \text{m}^2/\text{m}^3 once split into 100 smaller bubbles of the same total volume.

Gas hold-up (εg\varepsilon_g)

Gas hold-up (εg\varepsilon_g) is the fraction of the gas-liquid dispersion's volume occupied by gas bubbles at any instant — not the sparging flow rate, but how much gas is actually held up inside the vessel while it rises through the liquid.

εg=VgVg+VL\varepsilon_g = \frac{V_g}{V_g + V_L}
  • εg\varepsilon_g — gas hold-up (fraction of the dispersion's volume occupied by gas bubbles at any instant — how much gas is retained in the liquid, not the sparging flow rate), dimensionless (fraction, 00–11)
  • VgV_g — volume of gas (bubbles) in the dispersion, in cubic metres (m3\text{m}^3)
  • VLV_L — volume of liquid in the dispersion, in cubic metres (m3\text{m}^3)

εg\varepsilon_g is what turns single-bubble geometry (6/d6/d for one sphere) into a vessel-level interfacial area. The general form folds it in directly:

a=6εgdba = \frac{6\varepsilon_g}{d_b}
  • aa — interfacial area per unit volume (total gas-liquid contact surface available for mass transfer, per unit liquid volume), in square metres per cubic metre (m2/m3\text{m}^2/\text{m}^3)
  • εg\varepsilon_g — gas hold-up (fraction of dispersion volume occupied by gas), dimensionless
  • dbd_b — Sauter mean bubble diameter, in metres (m\text{m})

The single-bubble worked example above used a=6/da = 6/d as shorthand for a fixed total gas volume — that special case works because εg\varepsilon_g is the same in both scenarios being compared, so it cancels out. For any absolute calculation, use the full a=6εg/dba = 6\varepsilon_g/d_b.

What raises εg\varepsilon_g (and generally aa):

  • Higher superficial gas (sparging) velocity — more gas fed per unit time and cross-section means more gas volume held up in the liquid at any instant.
  • Smaller bubbles — they rise more slowly, so they spend longer in the liquid before escaping, raising hold-up for the same gas flow rate.

What lowers it:

  • Coalescence — bubbles merging into larger ones increases dbd_b directly and makes bubbles rise faster (less residence time), lowering both εg\varepsilon_g and aa. Coalescence always works against interfacial area, never for it.
  • High viscosity — thick, non-Newtonian broths often promote coalescence and slow bubble breakup, generally reducing hold-up and aa.

Electrolytes — two separate effects, easy to conflate. Dissolved salts (e.g. NaCl, Na2_2SO4_4) do two unrelated things in a sparged system:

  • They inhibit bubble coalescence. Ions in solution stabilize the thin liquid film between two approaching bubbles, preventing them from merging. This keeps mean bubble diameter smaller and aa (and kLak_La) higher than in pure water, at the same gas flow rate.
  • They lower O2_2 solubility (C∗C^*) via salting-out (see "Oxygen solubility" below) — a completely separate mechanism (ionic strength suppressing gas solubility, not bubble size).

So electrolytes raise aa/kLak_La (via smaller, more stable bubbles) while simultaneously lowering C∗C^* (via salting-out) — two different variables, two different mechanisms, pulling in directions that don't cancel out predictably. There is no "salting-in" effect that raises O2_2 solubility here — that term describes a different phenomenon (certain organic solutes becoming more soluble with certain salts), not this case.

Worked example: electrolyte effects in a sparged bioreactor

Question. Which of the following effects occur when pure electrolytes (salts such as NaCl or Na2_2SO4_4) are added to pure water in a sparged bioreactor? (a) Electrolytes inhibit bubble coalescence. (b) Mean bubble diameter decreases at high salt concentrations. (c) Interfacial area aa increases. (d) Oxygen solubility C∗C^* increases due to the salting-in effect.

Solution.

  1. Ions stabilize the thin film between approaching bubbles, preventing them from merging — (a) is correct.
  2. Coalescence inhibition means bubbles stay smaller instead of merging into larger ones as salt concentration rises — (b) is correct.
  3. Smaller dbd_b at roughly the same hold-up means a=6εg/dba = 6\varepsilon_g/d_b rises — (c) is correct.
  4. Electrolytes lower O2_2 solubility via salting-out, not raise it via "salting-in" — (d) gets both the direction and the name of the effect wrong.

Answer: (a), (b), and (c) are correct; (d) is false.

Worked example: which statements about interfacial area are correct

Question. Select the CORRECT statements regarding the gas-liquid interfacial area per unit volume aa: (a) aa is related to gas hold-up (εg\varepsilon_g) and bubble diameter (dbd_b) by a=6εg/dba = 6\varepsilon_g/d_b. (b) High shear from Rushton turbines reduces mean bubble diameter, thereby increasing aa. (c) Increasing gas sparging velocity increases gas hold-up, generally increasing aa. (d) Coalescence of air bubbles increases aa.

Solution.

  1. (a) is exactly the general interfacial-area equation above — correct.
  2. Rushton turbines are high-shear radial impellers; that shear breaks bubbles into smaller ones, and since a∝1/dba \propto 1/d_b, a smaller dbd_b means a higher aa — (b) is correct.
  3. Higher sparging (superficial gas) velocity increases εg\varepsilon_g (more gas held up in the dispersion at any instant); with dbd_b roughly unchanged, a=6εg/dba = 6\varepsilon_g/d_b rises too — (c) is correct.
  4. Coalescence merges small bubbles into larger ones, increasing dbd_b. Since a∝1/dba \propto 1/d_b, this decreases aa — (d) is false; it directly contradicts (b) and (c), which both work by keeping bubbles small.

Answer: (a), (b), and (c) are correct; (d) is false.

Volumetric mass transfer coefficient (kLak_L a)

kLk_L and aa are hard to measure separately, so they are lumped into a single volumetric mass transfer coefficient, kLak_L a, in reciprocal seconds (s−1\text{s}^{-1}) — the same film-theory rate equation as above, NA=kLa (C∗−C)N_A = k_L a\,(C^* - C), just written with kLk_L and aa combined into one term instead of kept separate.

kLk_L and aa change together (e.g. faster impeller speed raises both), so they're rarely separated in practice. In aerobic fermentation, AA is oxygen, so kLak_L a represents the bioreactor's oxygen transfer capacity — the standard parameter for bioreactor design and scale-up.

Penetration and surface renewal theories

Penetration theory and surface renewal theory replace film theory's stagnant film with short-lived liquid elements that arrive at the interface, exchange solute by unsteady-state diffusion, then get swept back into the bulk before equilibrium is reached.

Higbie's penetration theory

Every liquid element sits at the interface for the same fixed contact time tct_c before being replaced by fresh bulk fluid.

kL=2DABπtck_L = 2\sqrt{\frac{D_{AB}}{\pi t_c}}
  • kLk_L — liquid film mass-transfer coefficient (how fast a solute molecule crosses into a fresh liquid element during its brief time at the interface), in metres per second (m/s\text{m/s})
  • DABD_{AB} — diffusion coefficient of the solute through the liquid, in square metres per second (m2/s\text{m}^2/\text{s})
  • tct_c — contact time of a liquid element at the interface before it is replaced, in seconds (s\text{s}) — assumed identical for every element

This gives kL∝DAB0.5k_L \propto D_{AB}^{0.5} — a square-root dependence, unlike film theory's linear kL∝DAB1.0k_L \propto D_{AB}^{1.0}.

Danckwerts' surface renewal theory

Danckwerts kept the same short-lived-element picture but dropped the assumption of equal contact time: elements are replaced at random, with a fractional rate of surface renewal ss describing the probability per unit time that any given element gets swapped out.

kL=DAB sk_L = \sqrt{D_{AB}\, s}
  • kLk_L — liquid film mass-transfer coefficient (how fast a solute molecule crosses into a fresh liquid element before it is swept away), in metres per second (m/s\text{m/s})
  • DABD_{AB} — diffusion coefficient of the solute through the liquid, in square metres per second (m2/s\text{m}^2/\text{s})
  • ss — fractional rate of surface renewal (probability per unit time that a given liquid element gets swapped out for fresh bulk fluid), in reciprocal seconds (s−1\text{s}^{-1})

Like Higbie's theory, kL∝DAB0.5k_L \propto D_{AB}^{0.5} — the two theories share the same diffusivity dependence even though they disagree on how contact time is distributed.

Comparing the three theories

TheoryInterface pictureContact time assumptionkLk_L vs. DABD_{AB}
Whitman film theorystagnant film, steady-state diffusionnot applicable (steady state)kL∝DAB1.0k_L \propto D_{AB}^{1.0}
Higbie penetration theoryshort-lived elements, unsteady diffusionsame fixed tct_c for every elementkL∝DAB0.5k_L \propto D_{AB}^{0.5}
Danckwerts surface renewal theoryshort-lived elements, unsteady diffusionrandom distribution of ages, rate sskL∝DAB0.5k_L \propto D_{AB}^{0.5}

The DABD_{AB} exponent identifies the theory family — 1.0 for film theory, 0.5 for either penetration or surface renewal. Distinguishing Higbie from Danckwerts then comes down to the contact-time assumption, not the exponent: fixed and equal (Higbie) vs. randomly distributed (Danckwerts).

Worked example: which theory does each statement describe?

Question. Which of the following statements comparing mass transfer theories are correct?

  • (a) Whitman's film theory predicts kLk_L is directly proportional to DABD_{AB} raised to the power 1.01.0.
  • (b) Higbie's penetration theory predicts kLk_L is proportional to DABD_{AB} raised to the power 0.50.5.
  • (c) Danckwerts' surface renewal theory assumes all liquid elements stay at the interface for the exact same contact time.
  • (d) Both Higbie's and Danckwerts' theories predict kLk_L is proportional to DAB\sqrt{D_{AB}}.

Solution.

  1. Film theory models the film as a fixed-thickness, steady-state diffusion path: kL=DAB/δk_L = D_{AB}/\delta, which is linear in DABD_{AB}. (a) is correct.
  2. Higbie's kL=2DAB/(πtc)k_L = 2\sqrt{D_{AB}/(\pi t_c)} is a square-root dependence. (b) is correct.
  3. "Exact same contact time for every element" is Higbie's assumption, not Danckwerts' — Danckwerts specifically replaced it with a random distribution of surface ages. (c) is incorrect, since it describes the wrong theory.
  4. Danckwerts' kL=DAB sk_L = \sqrt{D_{AB}\, s} is also a square-root dependence, matching Higbie's. (d) is correct.

Answer: (a), (b), and (d) are correct; (c) is incorrect.

Oxygen uptake in a bioreactor

Rate of oxygen transfer

NA=kLa (CAL∗−CAL)N_A = k_L a\, (C^*_{AL} - C_{AL})
  • NAN_A — volumetric rate of oxygen transfer (OTR — how fast oxygen crosses from gas bubbles into the broth; must keep pace with cellular consumption for aerobic growth), in gram-mole per cubic metre per second (gmol/(m3⋅s)\text{gmol}/(\text{m}^3 \cdot \text{s}))
  • kLak_L a — volumetric mass transfer coefficient (lumps film coefficient and interfacial area into one measurable term — the standard parameter for bioreactor aeration design and scale-up), in reciprocal seconds (s−1\text{s}^{-1})
  • CAL∗C^*_{AL} — saturation concentration of oxygen in the liquid (its maximum solubility), in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)
  • CALC_{AL} — instantaneous concentration of oxygen in the liquid, in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)

Henry's law and CAL∗C^*_{AL}

Oxygen is a sparingly soluble gas in aqueous media — at 25 ∘C25\,^\circ\text{C}, water in equilibrium with ambient air holds only about 77–9 mg/L9\ \text{mg/L} dissolved O2_2 (≈0.2\approx 0.2– 0.3 mmol/L0.3\ \text{mmol/L}). This small number is exactly why kLak_L a has to be so large in aerobic fermentation: the absolute driving force (CAL∗−CAL)(C^*_{AL} - C_{AL}) available at any moment is tiny.

CAL∗C^*_{AL} isn't a fixed property of oxygen — it's set by Henry's law, which relates the partial pressure of a gas above a liquid to its equilibrium (saturation) concentration dissolved in that liquid:

pO2=H CAL∗p_{O_2} = H\, C^*_{AL}
  • pO2p_{O_2} — partial pressure of oxygen in the gas phase, in pascals (Pa\text{Pa})
  • HH — Henry's law constant for oxygen in the liquid (measures how reluctant the gas is to dissolve — higher HH means lower solubility at the same partial pressure), in pascal-cubic-metre per gram-mole (Pa⋅m3/gmol\text{Pa}\cdot\text{m}^3/\text{gmol}) — depends on the gas, the liquid, and temperature
  • CAL∗C^*_{AL} — saturation concentration of oxygen in the liquid, in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)

For a fixed oxygen partial pressure in the sparged gas, CAL∗=pO2/HC^*_{AL} = p_{O_2}/H — so CAL∗C^*_{AL} moves opposite to HH.

Factors affecting CAL∗C^*_{AL}

Per CAL∗=pO2/HC^*_{AL} = p_{O_2}/H, anything that changes pO2p_{O_2} or HH changes CAL∗C^*_{AL}:

  1. Partial pressure of oxygen in the inlet gas, pO2p_{O_2}. Richer oxygen content in the sparged gas (e.g. oxygen-enriched air) directly raises pO2p_{O_2}, and so CAL∗C^*_{AL}.

  2. Total headspace/system pressure. By Dalton's law, pO2=yO2 Ptotalp_{O_2} = y_{O_2}\, P_{total} — for the same gas composition, raising the total pressure (as in a pressurized fermenter) raises pO2p_{O_2} in proportion, and so CAL∗C^*_{AL}. This is why pressurizing a fermenter is a real technique for boosting oxygen supply.

  3. Temperature. HH increases with temperature — dissolving a gas into a liquid is exothermic, so by Le Chatelier's principle, raising temperature shifts the gas–liquid equilibrium back toward the gas phase, and less gas stays dissolved:

    Temperature↑  ⇒  H↑  ⇒  CAL∗↓\text{Temperature}\uparrow \;\Rightarrow\; H\uparrow \;\Rightarrow\; C^*_{AL}\downarrow

    So a hotter broth has a lower, not higher, oxygen saturation concentration — the opposite of what intuition about "warm liquids dissolve more gas" might suggest. This is the same fact as the temperature factor for kLak_L a above, now traced back to its cause: it's HH rising with temperature that shrinks C∗C^* and squeezes the (C∗−C)(C^* - C) driving force.

  4. Dissolved solutes (salts, sugars). Dissolved electrolytes and sugars in the medium raise the effective HH for oxygen — the salting-out effect — so a real fermentation broth has a lower CAL∗C^*_{AL} than pure water at the same temperature and pO2p_{O_2}.

Worked example: which factors affect C∗C^*?

Question. Which of the following physical factors influence the saturation dissolved oxygen concentration (C∗C^*) in a bioreactor?

  • (a) Partial pressure of oxygen in the inlet gas stream
  • (b) Total bioreactor headspace pressure
  • (c) Temperature of the broth
  • (d) Concentration of dissolved salts and sugars in the medium

Solution. All four are genuine factors, via C∗=pO2/HC^* = p_{O_2}/H:

  1. (a) sets pO2p_{O_2} directly. Correct.
  2. (b) sets pO2p_{O_2} indirectly, through Dalton's law (pO2=yO2Ptotalp_{O_2} = y_{O_2} P_{total}). Correct.
  3. (c) sets HH — higher temperature raises HH and lowers C∗C^*. Correct.
  4. (d) also sets HH, via the salting-out effect — dissolved solutes raise HH and lower C∗C^* relative to pure water. Correct.

Answer: all four, (a)–(d), influence C∗C^*.

Worked example: true/false on oxygen transfer fundamentals

Question. Which of the following statements regarding oxygen transfer in aerobic bioreactors is/are TRUE?

  • (a) The driving force for oxygen transfer across the liquid film is (C∗−CL)(C^* - C_L).
  • (b) The addition of surfactants/antifoam agents generally decreases kLk_L by adding interfacial resistance.
  • (c) Smaller bubble size increases the gas-liquid interfacial area per unit volume aa.
  • (d) Dissolved oxygen saturation solubility C∗C^* increases with increasing temperature.

Solution.

  1. The film-theory rate equation is NA=kLa(C∗−C)N_A = k_L a (C^* - C) — the concentration gap between saturation and bulk is exactly the driving force. (a) is true.
  2. Antifoam agents are surfactants: they add a resistive film at the interface, lowering kLk_L. (b) is true.
  3. Smaller bubbles mean more total surface area for the same gas volume (a=6/da = 6/d for a sphere) — smaller dd gives larger aa. (c) is true.
  4. By Henry's law, C∗=pO2/HC^* = p_{O_2}/H, and HH increases with temperature (gas dissolution is exothermic) — so C∗C^* decreases with temperature, not increases. (d) is false.

Answer: (a), (b), and (c) are true; (d) is false.

Worked example: oxygen solubility claims

Question. Select the CORRECT statements regarding oxygen solubility in aqueous fermentation broth: (a) O2_2 is a sparingly soluble gas, with C∗C^* typically around 77–9 mg/L9\ \text{mg/L} at 25 ∘C25\,^\circ\text{C} under ambient air. (b) Increasing bioreactor temperature increases C∗C^*. (c) Increasing vessel headspace pressure increases C∗C^* according to Henry's law. (d) High concentrations of dissolved electrolytes and sugars decrease C∗C^* via the salting-out effect.

Solution.

  1. (a) is the standard textbook solubility value for O2_2 in water at room temperature under air — correct.
  2. HH increases with temperature (gas dissolution is exothermic), so C∗=pO2/HC^* = p_{O_2}/H falls as temperature rises, not rises — (b) is false.
  3. Henry's law makes C∗C^* directly proportional to pO2p_{O_2}; raising headspace pressure raises pO2p_{O_2} (Dalton's law) and therefore C∗C^* — (c) is correct.
  4. Dissolved electrolytes and sugars raise the broth's effective HH (salting-out), which lowers C∗C^* — (d) is correct.

Answer: (a), (c), and (d) are correct; (b) is false.

Path of an oxygen molecule from a gas bubble to a cell's reaction site, numbered in eight steps: (1) interior of the bubble to the gas-liquid interface, (2) across the gas-liquid interface, (3) through the liquid film around the bubble, (4) across the bulk liquid, (5) through the liquid film around the cell, (6) across the liquid-cell interface, (7) to an individual cell within a clump, (8) through the cytoplasm to the reaction site — with k_L a labeling the overall path

Oxygen travels from a gas bubble to a cell's reaction site through eight resistances in series:

  1. Transfer from the interior of the bubble to the gas-liquid interface
  2. Transfer across the gas-liquid interface
  3. Diffusion through the liquid film around the bubble
  4. Transport through the bulk liquid
  5. Diffusion through the liquid film around the cell
  6. Movement across the liquid-cell interface
  7. For cells in a clump, diffusion to the individual cell
  8. Transport through the cytoplasm to the reaction site

Step 3 — diffusion through the liquid film around the bubble — is the largest resistance, which is why kLak_L a (built around that film) is the parameter that governs overall oxygen transfer.

This holds for free unicellular cells in well-mixed suspension, and rests on two other facts about the eight steps:

  • Step 2 (gas-film resistance) is negligible. O2_2 diffuses roughly four orders of magnitude faster in the gas phase than in the liquid phase, so for a sparingly-soluble gas like O2_2 almost all the resistance sits on the liquid side of the interface — this is why only the liquid-film coefficient (kLk_L) appears in kLak_L a, with no matching gas-film term.
  • Step 5 (liquid film around the cell) is negligible for small, free cells. A single bacterium or yeast cell is small enough, and its own mass-transfer coefficient high enough, that this film contributes almost nothing next to step 3 — don't confuse it with the bubble's film.

Exception: mycelial pellets. For filamentous fungi growing as pellets or clumps (step 7 in the list), O2_2 must additionally diffuse through the biomass matrix to reach cells at the core. Past a critical pellet radius, this intra-pellet diffusion resistance dominates instead — large enough pellets can starve their core of oxygen entirely (an anoxic core), even when kLak_L a at the bubble is perfectly adequate. So "the liquid film around the bubble is always the dominant resistance" is only true for free-cell suspensions, not pelleted or flocculated growth.

Worked example: which resistance dominates, and when

Question. Which of the following statements about O2_2 transfer resistances from a gas bubble to a microbial cell in a stirred bioreactor are correct? (a) For unicellular bacteria in well-mixed suspension, the primary resistance is the liquid film around the gas bubble. (b) Gas-film resistance is practically negligible. (c) For fungal mycelial pellets, intra-pellet O2_2 diffusion is a significant resistance. (d) The liquid film around an individual bacterial cell is the dominant resistance for all aerobic bioprocesses.

Solution.

  1. For free unicellular cells, step 3 (liquid film around the bubble) is the largest of the eight resistances — (a) is correct.
  2. O2_2's gas-phase diffusivity is far higher than its liquid-phase diffusivity, so the gas film (step 2) contributes negligible resistance — (b) is correct.
  3. Mycelial pellets add an extra resistance step not present for free cells: diffusion through the pellet's own biomass to reach interior cells, which becomes dominant past a critical pellet size — (c) is correct.
  4. (d) conflates two different films: the film around the cell (step 5, negligible for small free cells) with the film around the bubble (step 3, the actual dominant resistance) — and ignores that pellets shift the dominant resistance elsewhere entirely. "For all aerobic bioprocesses" is false.

Answer: (a), (b), and (c) are correct; (d) is false.

Factors affecting cellular oxygen demand

  1. Species of bacteria consuming the oxygen
  2. Growth phase the cells are in
  3. Nature of the carbon source fed to the cells

Rate of oxygen consumption per cell

QO=qO XQ_O = q_O\, X
  • qOq_O — specific oxygen uptake rate: oxygen required by one cell, in gram-mole per gram (dry cell weight) per second (gmol/(g⋅s)\text{gmol}/(\text{g} \cdot \text{s}))
  • XX — biomass (cell) concentration, in grams (dry cell weight) per cubic metre (g/m3\text{g/m}^3)
  • QOQ_O — oxygen uptake rate of the whole bioreactor (OUR), in gram-mole per cubic metre per second (gmol/(m3⋅s)\text{gmol}/(\text{m}^3 \cdot \text{s}))

OTR equals OUR at steady state

At steady state, oxygen is consumed as fast as it is supplied, so the oxygen transfer rate equals the oxygen uptake rate: NA=QON_A = Q_O, which gives

kLa (CAL∗−CAL)=qO Xk_L a\, (C^*_{AL} - C_{AL}) = q_O\, X

Dependence on dissolved oxygen concentration

Specific oxygen uptake rate plotted against dissolved oxygen concentration C_AL — the curve rises from the origin, climbs steeply, then flattens into a plateau past a critical concentration C_crit

  • If CAL<CcritC_{AL} < C_{crit}: there isn't enough dissolved oxygen around the cell, so it can't take up oxygen any faster than that supply allows.
  • If CAL>CcritC_{AL} > C_{crit}: the cell is already taking up oxygen as fast as it can, so more dissolved oxygen doesn't speed it up further.

Two special conditions

Maximum kLak_L a — set CAL=CcritC_{AL} = C_{crit} in NA=QON_A = Q_O:

kLamax=qO XCAL∗−Ccritk_{L}a_{max} = \frac{q_O\, X}{C^*_{AL} - C_{crit}}

This is the minimum aeration/agitation capability a bioreactor must be designed for to keep the culture above CcritC_{crit} at a given cell load.

Maximum cell concentration — set CAL=0C_{AL} = 0 (no oxygen left in the bioreactor) in NA=QON_A = Q_O:

Xmax=kLa CAL∗qOX_{max} = \frac{k_L a\, C^*_{AL}}{q_O}

This gives the theoretical ceiling on biomass density a given aeration setup can sustain — the limit you'd hit before scaling up oxygen supply.

Factors affecting kLak_L a

1. Bubble size

There's a reference bubble size around 22–3 mm3\ \text{mm}:

  • Bubble size >2> 2–3 mm3\ \text{mm}: interfacial area aa drops, so kLak_L a decreases.
  • Bubble size <2< 2–3 mm3\ \text{mm}: surface tension turns tiny bubbles into rigid spheres, which stops internal circulation and slightly lowers kLk_L — but the resulting jump in interfacial area aa easily outweighs that drop, so kLak_L a increases overall.

Sparger design sets the bubble size at the source. A large single orifice sparger releases coarse bubbles; a sintered (micro)sparger forces gas through many fine pores instead, producing much smaller bubbles — pushing aa, and so kLak_L a, up.

2. Viscosity

  • Viscosity↑ → resistance to OTR increases → kLak_L a decreases.
  • Viscosity↓ → resistance to OTR decreases → kLak_L a increases.

3. Stirrer speed

kLak_L a increases directly with stirrer speed.

4. Temperature

Temperature pulls kLak_L a in two directions at once:

  • Helps the kinetics: temperature↑ → viscosity↓ and diffusivity DABD_{AB}↑ → kLk_L↑, so kLak_L a↑.

  • Hurts the driving force: temperature↑ → gas solubility C∗C^*↓ (dissolution is exothermic), shrinking (C∗−C)(C^* - C) in

    NA=kLa (C∗−C)N_A = k_L a\, (C^* - C)

So higher temperature speeds up transfer kinetics but weakens the driving force behind it.

5. Gas sparging rate and coalescence

  • Sparging (superficial gas) velocity↑ → gas hold-up εg\varepsilon_g↑ → aa (and kLak_L a)↑, since a=6εg/dba = 6\varepsilon_g/d_b — see "Gas hold-up" above. Equivalently: more gas fed in means more gas held up and more total interfacial area.
  • Coalescence → bubbles merge, dbd_b↑ → aa↓ → kLak_L a↓. Coalescence always reduces interfacial area — it works against everything sparging and impeller shear are trying to achieve.
  • Sparger design matters too: a sintered microsparger produces much finer bubbles than a large single orifice at the same flow rate, raising aa (and kLak_L a) for the same gas input.

6. Antifoam agents

Chemical antifoams (commonly silicone-based surfactants) control foam but come at a direct cost to kLak_L a, especially at high concentration:

  • They lower surface tension in a way that promotes bubble coalescence into fewer, larger bubbles (bigger dbd_b, lower aa).
  • They form a hydrophobic film at the gas-liquid interface that adds its own resistance to O2_2 diffusion, lowering kLk_L directly.

Excess antifoam is a common, easily overlooked way to cripple oxygen transfer even while foam is well controlled — dose the minimum needed, not the maximum tolerated.

Worked example: which changes raise kLak_L a?

Question. Which of the following operational modifications will result in an increase in kLak_L a in an agitated bioreactor? (a) Increasing the impeller rotational speed (NiN_i) (b) Increasing the sparging gas flow rate (c) Adding a chemical antifoam agent at high concentration (d) Replacing a large single orifice sparger with a sintered microsparger.

Solution.

  1. Higher impeller speed increases turbulence and breaks bubbles smaller — kLak_L a increases. (a) is correct.
  2. Higher gas flow rate means more gas hold-up and more total interfacial area — kLak_L a increases. (b) is correct.
  3. High-concentration antifoam promotes coalescence and adds interfacial resistance, both of which lower kLak_L a rather than raise it. (c) is incorrect.
  4. A sintered microsparger produces much finer bubbles than a large single orifice, raising aa — kLak_L a increases. (d) is correct.

Answer: (a), (b), and (d) increase kLak_L a; (c) decreases it.

Worked example: which factors decrease kLak_L a

Question. Which of the following factors will cause a DECREASE in kLak_L a in an agitated fermenter? (a) Addition of excess chemical antifoam agents (b) Increase in broth apparent viscosity due to mycelial fungal growth (c) Increasing the agitation speed from 200200 rpm to 500500 rpm (d) Coalescence of smaller gas bubbles into larger bubbles.

Solution.

  1. Excess antifoam promotes coalescence and adds its own interfacial resistance — both lower kLak_L a — (a) decreases it.
  2. Higher viscosity increases resistance to oxygen transfer, lowering kLak_L a — (b) decreases it.
  3. kLak_L a increases directly with stirrer speed — raising rpm from 200200 to 500500 increases kLak_L a, it doesn't decrease it — (c) is not a decreasing factor.
  4. Coalescence increases dbd_b, which lowers aa and therefore kLak_L a — (d) decreases it.

Answer: (a), (b), and (d) decrease kLak_L a; (c) increases it instead.

Methods of measuring kLak_L a

1. Oxygen (gas) balance method

Measure the oxygen flow rate into and out of the fermenter, and back out NAN_A from the difference:

NA=1VL[(Fg CAG)i−(Fg CAG)o]N_A = \frac{1}{V_L}\Big[(F_g\, C_{AG})_{i} - (F_g\, C_{AG})_{o}\Big]
  • VLV_L — volume of liquid in the fermenter, in cubic metres (m3\text{m}^3)
  • FgF_g — volumetric gas flow rate, in cubic metres per second (m3/s\text{m}^3/\text{s})
  • CAGC_{AG} — concentration of oxygen in the gas phase, in gram-mole per cubic metre (gmol/m3\text{gmol/m}^3)
  • subscripts ii, oo — inlet and outlet gas streams

This measured NAN_A is then equated to NA=kLa (CAL∗−CAL)N_A = k_L a\, (C^*_{AL} - C_{AL}) to solve for kLak_L a.

Why it's useful: this method is non-invasive — it samples the gas streams, not the broth — and can run continuously throughout an active fermentation without interrupting aeration or agitation, unlike the dynamic gassing-out method below. It needs accurate measurement of both inlet and outlet gas flow rates and compositions (FgF_g and CAGC_{AG} on each side); it does not need any knowledge of liquid film thickness — that's a film-theory construct behind kLk_L, not something this macroscopic gas-phase balance ever measures directly.

Simplifying the gas flow rates: in general FgF_g can differ between inlet and outlet, because O2_2 is consumed and CO2_2 is produced at different molar rates. If the respiratory quotient (RQ=RQ = moles CO2_2 produced // moles O2_2 consumed) is close to 11, the loss of O2_2 moles is offset almost exactly by the gain of CO2_2 moles, so the total inlet and outlet volumetric gas flow rates end up nearly equal (Fg,i≈Fg,oF_{g,i} \approx F_{g,o}) — simplifying the balance considerably.

Worked example: what the off-gas balance method does and doesn't need

Question. Select the CORRECT statements regarding the off-gas oxygen balancing method for measuring OTR: (a) It is a non-invasive, continuous method that can be applied during active fermentation. (b) It requires accurate measurement of inlet and exit gas flow rates and gas compositions. (c) It depends heavily on knowing the exact value of the liquid film thickness. (d) If respiratory quotient (RQ) is close to 11, inlet and outlet total volumetric gas flow rates are nearly equal.

Solution.

  1. The method only samples gas streams, never the broth, and runs continuously without interrupting the process — (a) is correct.
  2. The defining equation needs exactly FgF_g and CAGC_{AG} on both the inlet and outlet sides — (b) is correct.
  3. Film thickness is a film-theory construct behind kLk_L, not something this macroscopic gas-phase mass balance ever measures — (c) is false.
  4. RQ≈1RQ \approx 1 means moles of O2_2 consumed ≈\approx moles of CO2_2 produced, so the total gas-phase mole (and volumetric flow) count barely changes between inlet and outlet — (d) is correct.

Answer: (a), (b), and (d) are correct; (c) is false.

2. Dynamic gassing-out method (most important)

The oxygen supply is interrupted at t0t_0, dissolved oxygen falls, supply is restored, and CALC_{AL} climbs back up toward a steady level CˉAL\bar{C}_{AL}. Two points on that recovery curve, (t1,CAL1)(t_1, C_{AL_1}) and (t2,CAL2)(t_2, C_{AL_2}), give:

Dissolved oxygen C_AL plotted against time: a flat line drops sharply once air flow is stopped at t0 as the culture consumes oxygen (OUR), bottoms out, then climbs back up through two marked recovery points (t1, C_AL1) and (t2, C_AL2) toward a steady reoxygenated level C̄_AL

kLa=ln⁡(CˉAL−CAL1CˉAL−CAL2)t2−t1k_L a = \frac{\ln\left(\dfrac{\bar{C}_{AL} - C_{AL_1}}{\bar{C}_{AL} - C_{AL_2}}\right)}{t_2 - t_1}
  • kLak_L a — volumetric mass transfer coefficient (the aeration parameter being measured, read off the slope of the reoxygenation recovery curve), in reciprocal seconds (s−1\text{s}^{-1})
  • t0t_0 — time at which the oxygen supply is cut off
  • t1t_1, t2t_2 — two time points during reoxygenation, at which the instantaneous oxygen concentrations are CAL1C_{AL_1} and CAL2C_{AL_2}
  • CˉAL\bar{C}_{AL} — oxygen level reached at steady state after reoxygenation

(ln⁡\ln is the natural log; if only base-10 log tables are available, ln⁡x=2.303log⁡x\ln x = 2.303 \log x.)

This shared formula covers two distinct named variants, which differ in how de-oxygenation happens and — critically — in what CˉAL\bar{C}_{AL} actually equals.

Variant A: dynamic method with respiring cells (Taguchi–Humphrey)

De-oxygenation happens through the culture's own respiration: cutting off the air supply is what starves the broth of oxygen, since the cells themselves consume the residual dissolved oxygen (this decline phase can even be used separately to read off the OUR). Air is then switched back on and CALC_{AL} climbs back toward CˉAL\bar{C}_{AL}.

CˉAL\bar{C}_{AL} vs. CAL∗C^*_{AL} — don't conflate them. Because the respiring culture keeps consuming oxygen even during recovery, the level CALC_{AL} actually settles at (CˉAL\bar{C}_{AL}) sits below the true no-consumption saturation value (CAL∗C^*_{AL}): CˉAL<CAL∗\bar{C}_{AL} < C^*_{AL}. Using CAL∗C^*_{AL} in the ratio instead of CˉAL\bar{C}_{AL} gives the wrong kLak_L a. This variant specifically requires a living, respiring culture — unlike the variant below, and unlike the sulphite oxidation method, which needs oxygen consumption to be purely chemical, not biological.

Variant B: physical gassing-out method (nitrogen stripping)

The alternative way to de-oxygenate: actively strip dissolved O2_2 by sparging nitrogen gas through the broth, with nothing biological consuming it. Air is then switched back on, and because nothing consumes oxygen during recovery, the system approaches the true saturation value directly — here CˉAL=CAL∗\bar{C}_{AL} = C^*_{AL}, so the formula can safely be written with C∗C^*. This variant is useful when respiration alone is too slow, or when characterizing a broth's kLak_L a (e.g. its rheology effects) independent of any organism — it can be run on cell-free medium or before inoculation.

Dissolved salts don't interfere with either variant — the method works in real fermentation media exactly as it does in plain water.

Probe lag is a real limitation, in either variant. The method assumes the dissolved-oxygen probe responds instantly to changes in CALC_{AL}, but every DO probe has its own response time τp\tau_p (set by its membrane and internal electrode geometry). If τp\tau_p is comparable to or larger than 1/kLa1/k_L a, the probe can't keep up with the true recovery curve and its reading lags behind reality — the raw data must then be corrected for probe dynamics (e.g. fitting a first-order probe-response model) before kLak_L a can be extracted correctly. A fast probe (τp≪1/kLa\tau_p \ll 1/k_L a) avoids the correction entirely.

What 1/kLa1/k_L a physically means. kLak_L a has units of s−1\text{s}^{-1}, so its reciprocal is a time constant — the characteristic time scale of reoxygenation, exactly like τ\tau in any first-order approach to equilibrium. After one time constant (t=1/kLat = 1/k_L a) the system has closed 63%63\% of the gap to CˉAL\bar{C}_{AL}; after three, ∼95%\sim95\%; after five, essentially done. A well-aerated reactor with kLa=0.1 s−1k_L a = 0.1\ \text{s}^{-1} reoxygenates on a 1/kLa=10 s1/k_L a = 10\ \text{s} timescale; a poorly aerated one with kLa=0.001 s−1k_L a = 0.001\ \text{s}^{-1} takes 1/kLa≈17 min1/k_L a \approx 17\ \text{min}. Comparing τp\tau_p to 1/kLa1/k_L a is just comparing the probe's own reaction speed to the reaction speed of the thing it's trying to measure — a probe with τp=8 s\tau_p = 8\ \text{s} is badly lagged against the fast 10 s10\ \text{s} case above but effectively instantaneous against the slow 17 min17\ \text{min} one.

Worked example: which assertions about the respiring-cells variant are correct?

Question. Select the CORRECT assertions regarding the dynamic gassing-out method (Variant A, with a respiring culture) for determining kLak_L a:

  • (a) The broth is first de-oxygenated by sparging nitrogen gas.
  • (b) Dissolved oxygen is monitored as a function of time after air sparging is resumed.
  • (c) The slope of ln⁡[(C∗−C)/C∗]\ln[(C^* - C)/C^*] versus time yields kLak_L a.
  • (d) The method cannot be applied if living, respiring microorganisms are present.

Solution.

  1. In this variant, de-oxygenation comes specifically from the culture's own respiration once air is cut off — nitrogen sparging belongs to Variant B, not this one. (a) is incorrect here.
  2. This is exactly the recovery phase the method measures: air resumes and CALC_{AL} is tracked climbing back up. (b) is correct.
  3. The correct ratio uses CˉAL\bar{C}_{AL} (the level actually reached under respiration), not CAL∗C^*_{AL} (the no-consumption saturation value) — see the distinction above. (c) is incorrect as stated.
  4. The opposite is true: this variant is specifically designed to be used with a living, respiring culture. (d) is incorrect.

Answer: (b) only.

Worked example: correctness checks across both variants

Question. Regarding the dynamic gassing-out method for determining kLak_L a in a fermenter, which of the following statements is/are CORRECT? (a) De-oxygenation of the broth is achieved by sparging nitrogen gas or shutting off aeration while respiring cells consume oxygen. (b) The slope of the line obtained by plotting ln⁡[(C∗−C1)/(C∗−C2)]\ln[(C^*-C_1)/(C^*-C_2)] against (t2−t1)(t_2-t_1) during re-aeration directly yields kLak_L a. (c) If the DO probe response time (τp\tau_p) is comparable to or larger than 1/kLa1/k_L a, probe lag must be corrected. (d) The method cannot be used if the liquid broth contains dissolved salts.

Solution.

  1. Both routes are real — they're just the two different named variants above (respiring-cells vs. nitrogen-stripping), not two options within a single technique. (a) is correct as a statement about the method family as a whole.
  2. This is the two-point slope form of the recovery curve. Note that C∗C^* here is only literally exact for Variant B (nitrogen-stripped); for Variant A the correct baseline is CˉAL\bar{C}_{AL}, not the true C∗C^* (see above). Read at the level this question intends, the underlying "slope of ln⁡(driving force)\ln(\text{driving force}) vs. time" technique is what yields kLak_L a either way. (b) is correct, with that caveat in mind.
  3. A slow DO probe can't track the true recovery curve once its own response time approaches the process time scale 1/kLa1/k_L a, so its lag must be corrected for — (c) is correct.
  4. Dissolved salts are present in essentially every real fermentation broth and don't interfere with this method — (d) is false.

Answer: (a), (b), and (c) are correct; (d) is false — but remember which baseline (C∗C^* vs. CˉAL\bar{C}_{AL}) is actually correct depends on which variant you're running.

3. Sulphite oxidation method

Sodium sulphite reacts with dissolved oxygen as fast as it enters the liquid, catalysed by trace amounts of Cu2+^{2+} or Co2+^{2+} (added as a copper or cobalt sulphate salt):

Na2SO3+0.5 O2→Cu2+/Co2+Na2SO4\text{Na}_2\text{SO}_3 + 0.5\, \text{O}_2 \xrightarrow{\text{Cu}^{2+}/\text{Co}^{2+}} \text{Na}_2\text{SO}_4

This only works because, at moderate-to-high sulphite concentrations (above roughly 0.04 kmol/m30.04\ \text{kmol/m}^3), the reaction is zero-order with respect to sulphite — the rate is independent of how much sulphite is left, so it stays fast enough to consume oxygen the instant it dissolves. That's what makes the method valid: the measured rate reflects pure mass transfer, not reaction kinetics. (At low sulphite concentration this breaks down — the reaction becomes rate-limiting instead, and the method no longer measures kLak_L a correctly.)

Reaction order describes how sensitive a reaction's rate is to a reactant's concentration, from the empirical rate law rate=k[A]n\text{rate} = k[\text{A}]^n. It's the reason zero-order is a requirement here, not just an incidental fact:

Order (nn)Rate lawIf sulphite were consumed during the assayEffect on the method
Zero (n=0n=0, actual)rate =k=kRate stays constant — sulphite level doesn't matterWorks: measured rate = pure O2_2 mass transfer
Half (n=0.5n=0.5)rate =k[SO32−]0.5=k[\text{SO}_3^{2-}]^{0.5}Rate drifts mildly as sulphite depletesPartially confounded — some kinetic drift mixed into the reading
First (n=1n=1)rate =k[SO32−]=k[\text{SO}_3^{2-}]Rate falls proportionally with sulphiteConfounded — can't tell a rate change from less O2_2 arriving vs. less sulphite left
Second (n=2n=2)rate =k[SO32−]2=k[\text{SO}_3^{2-}]^2Rate falls off sharply (quadratically)Worst case — small sulphite depletion swings the rate hard; kinetics dominate over mass transfer

As nn rises, the measured oxidation rate stops being a clean proxy for kLak_L a and starts reflecting how much sulphite is left instead. The method is only valid in the concentration range where the reaction has settled into zero order — the one regime where the chemistry gets out of the way and lets the measurement see the physics (mass transfer) directly.

Because the sulphite binds oxygen immediately, no free oxygen accumulates in solution — CAL=0C_{AL} = 0 — so the driving force collapses to just C∗C^*:

OTR=kLa C∗⇒kLa=OTRC∗\text{OTR} = k_L a\, C^* \qquad\Rightarrow\qquad k_L a = \frac{\text{OTR}}{C^*}

Worked example: catalyst and reaction order in the sulphite method

Question. During kLak_L a measurement by the sulphite oxidation method, sodium sulphite is converted to sodium sulphate in the presence of a catalyst. The catalyst commonly used and the reaction order with respect to sulphite at moderate concentrations are: (a) Cu2+^{2+}, first-order (b) Co2+^{2+}, zero-order (c) Fe2+^{2+}, second-order (d) Zn2+^{2+}, half-order.

Solution.

  1. Co2+^{2+} (cobaltous sulphate) is a standard catalyst for this reaction, alongside Cu2+^{2+} — rules out (c) and (d), which use catalysts not associated with this method.
  2. The method's validity depends on the reaction being fast enough to never be rate-limiting: at moderate-to-high sulphite concentrations, the rate is independent of sulphite concentration — zero-order.
  3. (a) pairs a valid catalyst with the wrong order (first-order would make the reaction rate-limiting, defeating the method's purpose) — so only (b) pairs a valid catalyst with the correct kinetics.

Answer: (b) Co2+^{2+}, zero-order.

Key points

  • Fick's law (JA=−DAB dCA/dyJ_A = -D_{AB}\, dC_A/dy) always runs high concentration to low; film theory assumes a linear concentration profile across each stagnant film — steady-state diffusion, no accumulation, same reasoning as linear conduction through a wall.
  • Film theory gives kL∝DAB1.0k_L \propto D_{AB}^{1.0}; penetration (Higbie) and surface renewal (Danckwerts) theories both give kL∝DAB0.5k_L \propto D_{AB}^{0.5} — the two square-root theories differ only in their contact-time assumption: fixed and equal (Higbie) vs. randomly distributed (Danckwerts).
  • At the phase boundary itself there is no resistance — instantaneous equilibrium — and it's the partition coefficient, not diffusion, that decides the concentration jump between the two phases there.
  • kLk_L and aa move together and can't be measured separately, so they're lumped into kLak_La, the volumetric mass transfer coefficient — in aerobic fermentation this is the reactor's oxygen transfer capacity, the standard scale-up parameter.
  • Of the eight resistances between a gas bubble and a cell's reaction site, diffusion through the liquid film around the bubble is the largest — which is why kLak_La (built around that film) governs overall O2_2 transfer.
  • At steady state, OTR = OUR: kLa(CAL∗−CAL)=qOXk_La(C^*_{AL} - C_{AL}) = q_O X. Setting CAL=CcritC_{AL} = C_{crit} gives the minimum kLak_La a reactor needs to design for; setting CAL=0C_{AL} = 0 gives the maximum biomass a given aeration setup can sustain.
  • Temperature is a trade-off for kLak_La: it speeds up the kinetics (viscosity↓, diffusivity↑ → kLk_L↑) but shrinks the driving force (solubility C∗C^*↓, since dissolution is exothermic).
  • The dynamic gassing-out method is the standard way to measure kLak_La: cut the air supply, let dissolved O2_2 fall, restart it, and read kLak_La off the slope of the recovery curve between two timepoints.