Pitching & Inoculum — Mathematical Formulations Science · v4.5
Viability — applied only at the commercial purchase point
Viability is a specification of commercial yeast (ADY or CmY): the fraction of total cells/g that are metabolically active. Dead cells don't multiply during propagation and don't ferment. Viability therefore affects one thing only — how much commercial yeast must be purchased to deliver the required number of viable cells to the stage that receives it.
Fermenter cell concentration is determined by inoculum target Xi and working volume V alone, not by the commercial viability spec.
Direct pitch (Ferm is the purchase point):
mcommercial (g) = (Xi · Vw,mL) / (ncells/g · vviability)
Staged (propagation or HDYC chain):
mFerm pitch, viable (g) = (Xi · Vw,mL) / ncells/g ← no viability
The biomass arriving at Ferm is fresh propagate from the upstream stage — assumed essentially 100% viable. Viability of commercial yeast enters only at the most upstream stage (the starter that IS purchased).
mcommercial starter (g) = mstage initial, viable / vviability
Three kinetic models per stage
Pre-Fermentation, HDYC, and the Fermentation Growth Window each independently select from Exponential, Logistic, or Gompertz. Each has a forward (Xi → Xf) and inverse (Xf → Xi) formulation used by Design back-calc, Target back-calc, and Seed-Train Forward modes.
Exponential (unbounded)
Xf = Xi · exp(µ · t) ⟺ Xi = Xf / exp(µ · t)
Sanity-check for short durations only. Overshoots reality when run through a full-duration growth phase because it has no substrate / ethanol / nutrient ceiling.
Logistic (default)
Sigmoidal growth toward an empirical ceiling Xmax that lumps ethanol inhibition, substrate depletion, oxygen limitation, and nutrient exhaustion:
dX/dt = µ · X · (1 − X / Xmax)
Closed-form forward and inverse solutions:
X(t) = X0 · Xmax / [ X0 + (Xmax − X0) · exp(−µ·t) ]
X0 = Xf / [ exp(µ·t) − (Xf/Xmax)·(exp(µ·t) − 1) ]
Gompertz (Zwietering modified)
Adds an explicit lag phase λ before exponential entry — captures rehydration, anaerobic adaptation, or stress recovery:
X(t) = X0 · (Xmax/X0)G(t), G(t) = exp(−exp(µ·e · (λ − t)/A + 1)), A = ln(Xmax/X0)
The forward function is monotone in X0; the inverse is solved by bisection on X0 ∈ (0, Xf] to machine precision.
Generation count Z
Z = log2(Xf / Xi)
Valid across all three kinetic models. For Exponential reduces to Z = µ·t / ln(2).
Fermentation Growth Window
The active growth window is the first 18–22 h of fermentation when cells divide; after this the culture is stationary-phase and ethanol production continues without further division for another 30–50 h. The calculator forward-models (Diagnostic) or back-calculates (Target) this window only.
Mass ⇄ cell-density conversion
Mkg = Xcells/mL · VL / ncells/g ⟺ Xcells/mL = Mkg · ncells/g / VL
Biomass is tracked in the same format-consistent mass basis (kg ADY or kg CmY) as Xmax and the upstream stages; cellsPerG is format-aware via the Yeast Format selector.
Diagnostic mode (forward)
Given end-of-fill Xi, forward-model to end-of-growth density Xf,growth. Verdict compares Xf,growth against the healthy industrial range 1–3×10⁸ cells/mL.
Target mode (back-calc)
User specifies Target Xf (desired end-of-growth density); inverse kinetic solves for required Xi, which then drives upstream seed-train sizing. Infeasibility is flagged when Target Xf ≥ Xmax,ferm.
Seed-Train Forward chain (HDYC topology)
Given commercial starter mass and HDYC/Pre-Ferm/Ferm volumes + kinetics, forward-compute through the chain:
X0,HDYC = mcomm · vviab
Xend,HDYC = forward(X0,HDYC, µHDYC, tHDYC, Xmax,HDYC, λHDYC)
X0,PF = Xend,HDYC · ηHDYC
Xend,PF = forward(X0,PF, µPF, tPF, Xmax,PF, λPF)
mpitch,ferm = Xend,PF · ηPF
Xi,ferm = mpitch,ferm · ncells/g / Vferm,L
Infeasibility at any stage (X0 ≥ Xmax) halts the chain with a descriptive error. Saturation (Xend > 0.98·Xmax) triggers a soft warning — the stage has no headroom and more time won't yield more biomass.
Transfer efficiency
ηk captures physical losses (residual in lines, transfer pumps, etc.) at stage boundary k. In back-calc direction, upstream viable biomass demand scales by 1/η. In forward direction, transferred mass scales by η directly.
Capacity mode (reverse direction)
Given available yeast mass mavailable and target Xi:
Vw,max (mL) = (mavailable (g) · ncells/g · vviab,eff) / Xi
vviab,eff is commercial viability in direct-pitch mode, 1 in staged modes (where the available mass represents propagated biomass).
Fill-cycle geometry (Semi-Batch only)
Fermenter fills linearly from Vheel to Vworking over tfill hours. Yeast is pulsed in over txfer starting at txferStart. Peak cell density occurs at transfer-complete (all cells in, minimum volume):
Vfill-frac(t) = Vheel + (1 − Vheel) · t / tfill
V@xferEnd = Vworking · Vfill-frac(txferStart + txfer)
Peak Xi = (Xi,target · Vworking) / V@xferEnd
Dilution factor = Vworking / V@xferEnd = Peak Xi / Xi,target
At Batch strategy the fill cycle is instantaneous; this section is hidden and values display em-dashes. Current model treats growth during fill as negligible — realistic for short transfers and dilute early-fill conditions; a µfill overlay could be added in a future build.
Fed-batch feed profile (staged modes)
For biomass balance dX/dt = µ·X at constant µ:
F(t) = µ · V · Xviable(t) / (Yx/s · Sfeed)
Peak feed rate is at t = tend. Total sugar delivered:
Stotal = Xi,viable · V · (eµt − 1) / Yx/s
Process-type classification (Propagation auto-classify)
Auto-classified by µ vs the critical threshold µcrit ≈ 0.27 hr⁻¹ for S. cerevisiae:
µ > 0.27 hr⁻¹ → Batch — substrate-unlimited, respirofermentative (Crabtree-on).
µ ≤ 0.27 hr⁻¹ → Fed-Batch — substrate-limited respiratory (Crabtree suppressed).
Pre-Ferm and HDYC Process Strategy are user-editable (no auto-classify) because they're physically controlled by operator choice of feeding strategy, not just µ.
Batch vs Semi-Batch strategy (Fermentation)
- Semi-Batch — fermenter fills over tfill while yeast is pulsed in. Fill Cycle section active.
- Batch — full charge at t=0, instantaneous pitch. Fill Cycle section hidden; Xi represents pitch density at t=0 (equivalent to end-of-fill density in Semi-Batch since fill is instant).
Growth Window math is identical in both strategies — the distinction affects fill dynamics, not post-fill growth.
Xmax format-awareness (ADY ↔ CmY ≈ 1:5.5)
Commercial yeast format conversion preserves biomass in dry-cell-weight terms. When Yeast Format is toggled and Xmax matches a known default, the value auto-swaps. Defaults (ADY / CmY):
- Pre-Ferm: 40 g/L / 220 g/L
- HDYC: 60 g/L / 330 g/L
- Fermentation: 12 g/L / 66 g/L
Custom Xmax values are preserved across format changes (only known defaults auto-swap).
Sig-fig display
Outputs use 3 significant figures: formatSigFigs(1556.73, 3) = "1560", formatSigFigs(0.00423, 3) = "0.00423". Avoids misleading 2-decimal precision across 6 orders of magnitude.
Tolerant number parser
The inoculum input accepts: 10E6, 1e7, 10M, 10,000,000, 10 000 000, 1.0e+07. Internal canonical: scientific notation with uppercase E.
Sensitivity analysis
Per-stage sensitivity tables vary ±10% and ±20% on one parameter (µ or t) holding all else fixed, reporting the corresponding change in Xi mass in the display convention for that stage (commercial purchase for the upstream-most stage, viable biomass for intermediate). Intended as quick-and-dirty "how sensitive am I to this choice?" not a full Monte-Carlo.
Scope and caveats
- Kinetic models are lumped — Xmax absorbs ethanol inhibition, substrate depletion, O2 limitation, and nutrient exhaustion empirically without modeling them separately. For full mechanistic modeling (Monod, Levenspiel product inhibition, Yx/s-vs-µ Crabtree coupling), use the FermAxiom Ethanol Time-Course Simulator.
- Growth Window assumes cells stop dividing at tgrowth (default 20 h) — consistent with industrial observation that ethanol accumulation halts division by hour 18–22. The rest of the ferm run is stationary-phase and not modeled here.
- Fill cycle treats biomass as purely volumetric (no growth during fill). Realistic for short transfers; a µfill kinetic overlay would refine this if needed.
- Fed-batch feed profile assumes constant µ throughout the stage — in reality, µ may drift.
- Strain comparison is mass/cost only; physiological differences (stress tolerance, flocculation, byproduct profile) are not modeled.
- Fresh propagate is assumed ~100% viable; real propagated biomass typically has 95%+ viability but the commercial-viability specification does not apply to it.