Bioprocess Engineering Formulas: μmax, Kinetics & Cheat Sheet

By BioProcess Tools Team | March 26, 2026 | 7 min read | Last updated: March 2026

The Monod equation μ = μmax · S / (Ks + S) predicts substrate-limited growth, and OTR = kLa · (C* − CL) sets the aeration limit. This bioprocess engineering formulas cheat sheet collects every essential equation in one reference: μmax and Monod kinetics, mass balance and yield coefficients, fed-batch feeding, OTR/kLa, P/V impeller power, heat transfer, downstream sizing, and process economics. Each formula links to the corresponding free calculator, so you can plug in your own numbers.

Figure 1. Key bioreactor engineering parameters and their core equations, mapped to physical locations on the vessel. Each formula links to a dedicated calculator below.

1. Cell Growth Kinetics Formulas

μ = ln(X2 / X1) / (t2t1)
Specific growth rate (h&supmin;¹)
  X = cell concentration (cells/mL or g/L)
  t = time (h)
td = ln(2) / μ
Doubling time (h)
X(t) = X0 × eμt
Exponential growth model
  X0 = initial cell concentration
μ = μmax × S / (Ks + S)
Monod equation — substrate-limited growth
  μmax = maximum specific growth rate (h&supmin;¹)
  S = substrate concentration (g/L)
  Ks = half-saturation constant (g/L)
IVC = ∫ X(t) dt ≈ Σ ½(Xi + Xi+1) × (ti+1ti)
Integral of viable cells (cells·day/mL) — trapezoidal rule
qP = ΔP / ΔIVC
Cell-specific productivity (pg/cell/day or g/cell/h)
  P = product concentration (g/L)

Track Growth in Real Time

Log cell counts, viability, and metabolites. Get automatic growth rate and IVC calculations.

CellTrack App →

2. Mass Balance & Yield Equations

YX/S = ΔX / ΔS
Biomass yield on substrate (g biomass / g substrate)
YX/O = ΔX / ΔO2
Biomass yield on oxygen (g biomass / mol O2)
Carbon balance:
Substratein = Biomass + CO2 + Product + Residual substrate
RQ = CO2 produced / O2 consumed
Respiratory quotient (dimensionless)
  RQ = 1.0 for glucose fully oxidized
  RQ > 1.0 suggests overflow metabolism (e.g., ethanol, acetate)

3. Fed-Batch Feeding Formulas

Exponential feed:
F(t) = (μ × X0 × V0) / (YX/S × Sf) × eμt
  F(t) = feed rate (L/h)
  X0 = initial biomass concentration (g/L)
  V0 = initial volume (L)
  Sf = feed substrate concentration (g/L)
Linear feed:
F(t) = F0 + k × t
  F0 = initial feed rate (L/h)
  k = ramp rate (L/h²)
Substrate mass balance in fed-batch:
dS/dt = F × Sf / Vμ × X / YX/Sm × X
  m = maintenance coefficient (g substrate / g biomass / h)
📈

Calculate Your Feed Profile

Enter biomass, yield coefficients, and feed concentration to generate exponential and linear feed schedules.

Fed-Batch Calculator →

4. Oxygen Transfer

OTR = kLa × (C*CL)
Oxygen transfer rate (mg/L/h or mmol/L/h)
  kLa = volumetric mass transfer coefficient (h&supmin;¹)
  C* = saturated DO concentration (≈7 mg/L at 37°C in air)
  CL = actual dissolved oxygen (mg/L)
Van't Riet correlation (coalescing media):
kLa = 0.026 × (P/V)0.4 × vs0.5

Van't Riet correlation (non-coalescing media):
kLa = 0.002 × (P/V)0.7 × vs0.2

P/V in W/m³, vs in m/s, kLa in s&supmin;¹
OUR = qO2 × X
Oxygen uptake rate (mmol/L/h)
  qO2 = specific oxygen consumption rate (mmol/g/h)

Estimate kLa for Your Bioreactor

Enter vessel dimensions, impeller specs, and gas flow to calculate kLa and OTR.

OTR & kLa Estimator →

5. Impeller & Mixing

P = Np × ρ × N³ × Di&sup5;
Ungassed impeller power draw (W)
  Np = power number (dimensionless, e.g., 5.0 for Rushton)
  ρ = fluid density (kg/m³)
  N = impeller speed (rps)
  Di = impeller diameter (m)
Tip speed:
vtip = π × N × Di
N in rps, Di in m, vtip in m/s
Reynolds number:
Re = ρ × N × Di² / μ
  μ = dynamic viscosity (Pa·s)
  Re > 10,000 = fully turbulent (typical for aqueous fermentation)
Power per unit volume:
P/V = Ptotal / Vliquid
Typical ranges: 0.5–3 W/L (mammalian), 2–10 W/L (microbial)
Mixing time (approximate):
tmix ∝ (V / P)1/3
📊

Scale-Up Calculator

Compare constant P/V, tip speed, kLa, Re, and mixing time across scales.

Scale-Up Calculator →

6. Heat Transfer Formulas

Qmet = OUR × V × 460 kJ/mol
Metabolic heat generation (W)
  OUR in mol/L/s, V in L
  460 kJ/mol is the heat of combustion per mol O2
Q = U × A × LMTD
Heat transfer rate (W)
  U = overall heat transfer coefficient (W/m²/K)
  A = heat transfer area (m²)
  LMTD = log-mean temperature difference (K)
LMTD = (ΔT1ΔT2) / ln(ΔT1 / ΔT2)
  ΔT1 = Tbroth − Tcoolant,in
  ΔT2 = Tbroth − Tcoolant,out
Coolant flow rate:
coolant = Q / (Cp × ΔT)
  Cp = specific heat capacity (J/kg/K, water ≈ 4184)
  ΔT = coolant temperature rise (K)

Heat Transfer Calculator

Calculate metabolic heat load and verify your cooling system can handle it.

Heat Transfer Calculator →

7. Downstream Processing

DBC = mass bound / column volume
Dynamic binding capacity (mg/mL resin)
Measured at 10% breakthrough (DBC10%)
LRV = log10(load titer / filtrate titer)
Log reduction value — viral clearance
  LRV ≥ 4 per step is typical regulatory expectation
Yield = (massout / massin) × 100%
Step yield (%)
Pool volume:
Vpool = CV × Vcolumn
  CV = number of column volumes in the elution pool

Related tools for applying these downstream formulas:

8. Centrifugation

RCF = 1.118 × 10−5 × r × N²
Relative centrifugal force (× g)
  r = radius (cm)
  N = rotational speed (RPM)
Sigma factor (tubular bowl, simplified):
Σ = (2π × N² × L × r²) / g
  L = bowl length (m)
  r = bowl radius (m)
  g = 9.81 m/s²
Scale-up by Q/Σ equivalence:
(Q / Σ)lab = (Q / Σ)production
Maintain constant Q/Σ to preserve separation performance

Centrifugation Scale-Up

Calculate sigma factors and scale between lab and production centrifuges.

Centrifugation Calculator →

9. Process Economics Formulas

COGS/g = Total Cost / (Batches × Volume × Titer × Yield)
Cost of goods sold per gram of product
  Total Cost = raw materials + labor + facility + overhead
Productivity = Titer / Culture Duration
Volumetric productivity (g/L/day)
Annual output:
Output = Batches/year × V × Titer × DSP Yield
in kg/year or g/year
💰

Model Your Process Economics

Build a full cost model for your upstream and downstream process.

Fermentation Economics →

10. Frequently Asked Questions

What is μmax in bioprocess?

In bioprocess, μmax (the maximum specific growth rate) is the highest rate at which biomass can grow per unit biomass when no nutrient is limiting, with units of inverse time (h−1). It is the asymptotic ceiling of the Monod model μ = μmax · S / (Ks + S) as substrate concentration S becomes large. Typical μmax values are E. coli 0.8–2.0 h−1 in rich medium (LB, K-12 vs BL21) and 0.7–0.8 h−1 in minimal glucose medium; Pichia pastoris 0.15–0.30 h−1 on glycerol; Saccharomyces cerevisiae 0.45 h−1 on glucose; CHO 0.025–0.04 h−1; HEK293 0.025–0.035 h−1. μmax is measured in batch culture from the slope of ln(X) vs t during early exponential phase, where substrate is in clear excess. It anchors fed-batch design: the set growth rate is chosen below μmax to avoid overflow metabolism (acetate in E. coli, lactate in CHO).

What is bioprocess kinetics?

Bioprocess kinetics is the quantitative description of cell growth, substrate uptake, and product formation rates inside a bioreactor. The three core equations are: (1) growth: dX/dt = μ · X, with the Monod substrate dependence μ = μmax · S / (Ks + S); (2) substrate uptake: dS/dt = −μ/YX/S · X − ms · X, where YX/S is the biomass yield on substrate and ms is the maintenance coefficient; (3) product formation: dP/dt = qP · X, where qP (g product / g cell / h) follows growth-associated (Luedeking–Piret α · μ), non-growth-associated (β), or mixed forms. These rate equations are integrated against the mass balance for batch, fed-batch, or continuous (chemostat) reactors. At steady state in a chemostat, μ equals the dilution rate D by definition.

What is the Monod equation?

The Monod equation describes how specific growth rate depends on the concentration of a single limiting substrate: μ = μmax · S / (Ks + S). Here μmax is the maximum growth rate, S is the limiting substrate concentration (g/L), and Ks is the half-saturation constant (the S at which μ = μmax/2). At S much greater than Ks, μ approaches μmax; at S much less than Ks, the relationship is approximately first-order in S. Typical Ks values are very low for glucose: E. coli 4 mg/L, S. cerevisiae 25 mg/L; Pichia on glycerol 10–30 mg/L; CHO on glucose 50–100 mg/L. The Monod model is the foundation for chemostat steady-state design and fed-batch substrate-limited feeding strategies. See the Monod kinetics deep dive for derivations and worked examples.

What is the formula for kLa in a bioreactor?

kLa (the volumetric oxygen mass transfer coefficient, units h−1) is most commonly estimated using the Van’t Riet correlation: kLa = a · (P/V)b · (vs)c, where P/V is gassed power per unit volume (W/m³), vs is superficial gas velocity (m/s), and a, b, c are empirical constants. Typical exponents are b = 0.4–0.7 and c = 0.4–0.5; for coalescing media (water) a is around 0.026 with b = 0.4 and c = 0.5, while for non-coalescing media (salt solutions, broth) a is around 0.002 with b = 0.7 and c = 0.2. kLa can also be measured experimentally with the dynamic gassing-out method: degas the broth with N2, switch to air, fit dC/dt = kLa · (C* − C). The full kLa guide covers derivations, correlations, and three measurement methods.

How is power per volume (P/V) calculated for a bioreactor?

Power per unit volume (P/V, W/m³) is the standard scale-up criterion for geometrically similar bioreactors: P/V = (Np · N³ · Di5 · ρ) / V, where Np is the impeller power number (dimensionless), N is rotational speed (rev/s), Di is impeller diameter (m), ρ is broth density (kg/m³), and V is liquid volume (m³). Power numbers depend on impeller geometry in the turbulent regime: Rushton 6-bladed disc 5.0–5.5, pitched-blade (4-blade, 45°) 1.2–1.7, hydrofoil (elephant-ear / SC-3) 0.3–0.8. P/V scales as N³ · Di² / V at constant geometry — useful for scale-up because doubling the impeller diameter requires reducing speed by a factor of 22/3 to hold P/V constant. Typical industrial values are 1–5 kW/m³ for microbial fermentations and 30–100 W/m³ for mammalian cell culture to control shear.

What is the Reynolds number formula for a bioreactor?

The impeller Reynolds number, used to assess mixing regime in a stirred-tank bioreactor, is Re = (ρ · N · Di²) / μ, where ρ is broth density (kg/m³), N is impeller rotational speed (rev/s), Di is impeller diameter (m), and μ is dynamic viscosity (Pa·s). Regimes for stirred tanks: laminar Re < 10, transitional 10 < Re < 1×104, and turbulent Re > 1×104. Most industrial bioreactors operate firmly in the turbulent regime (Re ≈ 105–106), which is required for the Van’t Riet kLa correlation and Np-table power numbers to apply. Mycelial broths and high-viscosity cultures can drop into the transitional regime — in which case Np rises sharply and mixing time gets dramatically longer.

What is the most-used equation in bioprocess design?

The most-used equation in bioprocess design is the Monod equation, μ = μmax · S / (Ks + S), which links specific growth rate to limiting substrate concentration and underpins batch, fed-batch, and chemostat models. Runners-up by frequency are the volumetric oxygen transfer relation OTR = kLa · (C* − CL) for aeration sizing, the exponential feed rule F(t) = (μ · X0 · V0) / (YX/S · Sf) · eμt for fed-batch dosing, and the impeller power draw P = Np · ρ · N³ · Di5 for scale-up. Together with the yield-coefficient definitions (YX/S, YP/S), these four sit inside almost every upstream process model. Published bioreactor textbooks (Doran 2013; Shuler and Kargi 2017) build their entire upstream chapters around them, and every calculator on this cheat sheet reduces to one of these forms.

How do you calculate kLa in a stirred bioreactor?

In a stirred bioreactor, kLa is calculated either from an empirical correlation or from a direct dynamic experiment. The Van’t Riet correlation kLa = a · (P/V)b · vsc uses gassed power per unit volume P/V (W/m³) and superficial gas velocity vs (m/s); for water-like coalescing broths a is around 0.026 with b = 0.4 and c = 0.5, and for non-coalescing salt-rich media a is around 0.002 with b = 0.7 and c = 0.2. For direct measurement, the dynamic gassing-out method degases the broth with N2, switches to air, and fits dC/dt = kLa · (C* − C) to the recorded DO trace; the negative slope of ln(C* − C) versus time gives kLa in h−1. Typical microbial values reach 200–800 h−1 at industrial scale; mammalian bioreactors sit at 5–50 h−1 because impellers must run at lower tip speed to protect cells. See the OTR/kLa estimator to run both routes side by side.

What is the difference between yield coefficient YX/S and YP/S?

YX/S (yield of biomass on substrate, g biomass / g substrate) and YP/S (yield of product on substrate, g product / g substrate) both quantify how efficiently a limiting substrate is converted, but they track different sinks. YX/S measures how much cell mass forms per gram of substrate consumed: typical values are 0.4–0.5 g/g for E. coli on glucose, 0.3–0.5 g/g for S. cerevisiae, and 0.02–0.06 g dry cell weight per gram glucose for CHO. YP/S measures how much target product (protein, metabolite, or organic acid) forms per gram of substrate and is set by pathway stoichiometry plus flux distribution; ethanol from S. cerevisiae reaches 0.48 g/g (near the 0.51 g/g theoretical maximum), while a recombinant protein might sit at 0.05–0.20 g/g. Overall carbon closure requires YX/S + YCO2/S + YP/S + Yresidual/S = 1.0 on a carbon-mole basis, which is how off-gas RQ measurements validate a fed-batch run.

📚 Resources & Further Reading