Pasteurization Water Bath Residence Time

Reference ID: MET-F715 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

The pasteurization water‑bath residence‑time calculation determines the minimum time a product‑filled container must remain in a hot water bath so that the cumulative lethal effect (the “P‑value”) reaches a prescribed target. This is a core design tool in process‑engineering for dairy, juice, and other liquid‑food pasteurization lines, where the cold spot inside the container governs safety. The result is used to size conveyor speed, bath length, or to verify that an existing line meets regulatory lethality requirements.

Methodology & Formulas

The algorithm follows the physics implemented in the reference Python code, expressed with standard engineering symbols.

1. Define constants and geometry

\[ T_{\text{water}} = 85^{\circ}\text{C},\qquad T_{\text{init}} = 20^{\circ}\text{C},\qquad k = 0.6\,\text{W/m·K},\qquad \rho = 1000\,\text{kg/m}^{3},\qquad c_{p} = 4180\,\text{J/kg·K} \]

Thermal diffusivity: \[ \alpha = \frac{k}{\rho\,c_{p}} \]

Container dimensions (short cylinder): radius \(R\) and half‑height \(H\). Characteristic lengths are

\[ L_{c,\text{cyl}} = R,\qquad L_{c,\text{slab}} = H \]

2. Biot numbers (surface‑to‑volume heat‑transfer ratio)

\[ \text{Bi}_{\text{cyl}} = \frac{h\,R}{k},\qquad \text{Bi}_{\text{slab}} = \frac{h\,H}{k} \]

3. Fourier numbers (dimensionless time)

For a given elapsed time \(t\)

\[ \text{Fo}_{\text{cyl}} = \frac{\alpha\,t}{R^{2}},\qquad \text{Fo}_{\text{slab}} = \frac{\alpha\,t}{H^{2}} \]

4. Dimensionless center‑temperature solutions

Domain of validity: The formulas below assume the constant‑surface‑temperature limit (Bi → ∞). For practical engineering this is accurate when Bi > 40. If 0.1 ≤ Bi ≤ 40, the general Heisler series with Bi‑dependent eigenvalues must be used. If Bi < 0.1, a lumped‑capacitance model is sufficient.

Both slab and cylinder contributions are evaluated with a truncated Heisler series (first \(N\) terms). For the slab (infinite plane wall) the eigenvalues and coefficients are

\[ \lambda_{n}^{\text{slab}} = \frac{(2n-1)\pi}{2},\qquad A_{n}^{\text{slab}} = \frac{4(-1)^{\,n+1}}{(2n-1)\pi} \]

Resulting slab dimensionless temperature:

\[ \theta_{\text{slab}}(t) = \sum_{n=1}^{N} A_{n}^{\text{slab}}\, \exp\!\big[-(\lambda_{n}^{\text{slab}})^{2}\,\text{Fo}_{\text{slab}}\big] \]

For the cylinder (infinite solid cylinder) the eigenvalues are the positive zeros of the Bessel function \(J_{0}\) (denoted \(\lambda_{n}^{\text{cyl}}\)) and the coefficients are

\[ B_{n}^{\text{cyl}} = \frac{2}{\lambda_{n}^{\text{cyl}}\,J_{1}(\lambda_{n}^{\text{cyl}})} \]

Resulting cylinder dimensionless temperature:

\[ \theta_{\text{cyl}}(t) = \sum_{n=1}^{N} B_{n}^{\text{cyl}}\, \exp\!\big[-(\lambda_{n}^{\text{cyl}})^{2}\,\text{Fo}_{\text{cyl}}\big] \]

The overall dimensionless temperature at the geometric centre (product of the two solutions) is

\[ \theta_{\text{total}}(t) = \theta_{\text{cyl}}(t)\,\theta_{\text{slab}}(t) \]

5. Physical centre temperature

\[ T_{\text{center}}(t) = T_{\text{water}} + \bigl(T_{\text{init}}-T_{\text{water}}\bigr)\,\theta_{\text{total}}(t) \]

6. Lethal rate and cumulative P‑value

At each time increment \(\Delta t\) the instantaneous lethal rate is

\[ L(t) = 10^{\displaystyle\frac{T_{\text{center}}(t)-T_{\text{ref}}}{z}} \]

where \(T_{\text{ref}} = 85^{\circ}\text{C}\) (reference temperature) and \(z\) is the organism‑specific temperature increase that yields a ten‑fold change in lethality.

The cumulative P‑value after \(n\) steps is

\[ P(t) = \sum_{i=1}^{n} L(t_{i})\,\frac{\Delta t}{60} \]

(the division by 60 converts seconds to minutes). The calculation stops when

\[ P(t) \ge P_{\text{target}} \]

and the corresponding elapsed time \(t\) (converted to minutes) is the required residence time \(t_{\text{res}}\).

7. Translation to process parameters

If the bath length \(L_{\text{bath}}\) is known, the conveyor speed is

\[ v = \frac{L_{\text{bath}}}{t_{\text{res}}} \]

Conversely, for a fixed speed \(v\) the needed bath length is \(L_{\text{bath}} = v\,t_{\text{res}}\).

Empirical range checks

Criterion	Acceptable Range	Action if violated
Water temperature \(T_{\text{water}}\)	70 °C ≤ \(T_{\text{water}}\) ≤ 100 °C	Adjust bath set‑point
Biot number \(\text{Bi}\)	0.1 ≤ \(\text{Bi}\) ≤ 40	Use lumped‑capacitance if \(\text{Bi}<0.1\); if \(\text{Bi}>40\) switch to constant‑surface series (this sheet)
Fourier number \(\text{Fo}\)	\(\text{Fo}>0.2\) for first‑term approximation	Increase number of terms \(N\) or employ full series
z‑value \(z\)	5 °C ≤ \(z\) ≤ 10 °C (typical microorganisms)	Confirm organism‑specific data
Heat‑transfer coefficient \(h\)	50 W/m²·K ≤ \(h\) ≤ 500 W/m²·K (agitated bath)	Re‑evaluate agitation or bath geometry
Aspect ratio (height/diameter)	≤ 5	Refine geometric model (e.g., include end effects)

Summary of key equations

\[ \begin{aligned} \alpha &= \frac{k}{\rho\,c_{p}}\\[4pt] \text{Bi}_{\text{cyl}} &= \frac{h\,R}{k},\qquad \text{Bi}_{\text{slab}} = \frac{h\,H}{k}\\[4pt] \text{Fo}_{\text{cyl}} &= \frac{\alpha\,t}{R^{2}},\qquad \text{Fo}_{\text{slab}} = \frac{\alpha\,t}{H^{2}}\\[4pt] \theta_{\text{slab}}(t) &= \sum_{n=1}^{N}\frac{4(-1)^{n+1}}{(2n-1)\pi}\, \exp\!\big[-\big(\tfrac{(2n-1)\pi}{2}\big)^{2}\text{Fo}_{\text{slab}}\big]\\[4pt] \theta_{\text{cyl}}(t) &= \sum_{n=1}^{N}\frac{2}{\lambda_{n}^{\text{cyl}}J_{1}(\lambda_{n}^{\text{cyl}})}\, \exp\!\big[-(\lambda_{n}^{\text{cyl}})^{2}\text{Fo}_{\text{cyl}}\big]\\[4pt] \theta_{\text{total}}(t) &= \theta_{\text{cyl}}(t)\,\theta_{\text{slab}}(t)\\[4pt] T_{\text{center}}(t) &= T_{\text{water}} + (T_{\text{init}}-T_{\text{water}})\,\theta_{\text{total}}(t)\\[4pt] L(t) &= 10^{\displaystyle\frac{T_{\text{center}}(t)-T_{\text{ref}}}{z}}\\[4pt] P(t) &= \sum_{i} L(t_{i})\,\frac{\Delta t}{60}\\[4pt] t_{\text{res}} &: \; P(t_{\text{res}})\ge P_{\text{target}} \end{aligned} \]

The minimum residence time is determined by the thermal death time of the target pathogen, typically Listeria monocytogenes or Salmonella, relative to the cold-spot temperature of the product. Process engineers must follow these steps:

Identify the target log reduction required by regulatory standards.
Determine the D-value and z-value for the specific product matrix.
Calculate the come-up time required for the product core to reach the target pasteurization temperature.
Add the required hold time at the target temperature to the come-up time to establish the total residence time.

Several variables can impact the consistency of residence time, potentially leading to under-processing. Key factors include:

Conveyor speed fluctuations or mechanical slippage.
Water bath turbulence and flow patterns, which affect heat transfer coefficients.
Initial product temperature variations at the point of entry.
Product loading density, which can alter the effective thermal mass of the bath.

Validation requires a combination of physical testing and mathematical modeling to ensure the process is robust. Engineers should perform the following:

Conduct heat penetration studies using calibrated thermocouples placed at the slowest heating point of the product.
Perform microbial challenge studies using surrogate organisms to verify log reduction.
Establish a safety margin by adding a buffer to the calculated minimum residence time to account for process deviations.
Document all validation data to support the HACCP plan and regulatory compliance.

Worked Example: Pasteurization Water Bath Residence Time

Scenario:
A juice bottle (short cylinder, radius \(R = 0.035\ \text{m}\), half-height \(H = 0.075\ \text{m}\)) enters a hot water bath at \(T_{\text{water}} = 85.0\ ^\circ\text{C}\) with initial product temperature \(T_{\text{init}} = 20.0\ ^\circ\text{C}\). To illustrate the constant‑surface‑temperature limit, the bath is assigned an extremely high heat transfer coefficient \(h = 10000.0\ \text{W/m}^2\text{K}\) (this is not typical of a simple agitated water bath; it approximates condensing‑steam conditions and validates the Bi → ∞ assumption). The product’s thermal properties are water-like: thermal conductivity \(k = 0.6\ \text{W/mK}\), density \(\rho = 1000.0\ \text{kg/m}^3\), specific heat \(c_p = 4180.0\ \text{J/kgK}\), so thermal diffusivity \(\alpha = k/(\rho c_p) = 1.4354 \times 10^{-7}\ \text{m}^2/\text{s}\). The target pasteurization value is \(P_{\text{target}} = 5.0\ \text{minutes}\) at a reference temperature of \(85\ ^\circ\text{C}\) with \(z = 7.0\ ^\circ\text{C}\) for the target microorganism.

Knowns:

Water temperature, \(T_{\text{water}} = 85.0\ ^\circ\text{C}
Initial product temperature, \(T_{\text{init}} = 20.0\ ^\circ\text{C}
Thermal conductivity, \(k = 0.6\ \text{W/mK}
Density, \(\rho = 1000.0\ \text{kg/m}^3
Specific heat, \(c_p = 4180.0\ \text{J/kgK}
Thermal diffusivity, \(\alpha = 1.4354 \times 10^{-7}\ \text{m}^2/\text{s}
z-value, \(z = 7.0\ ^\circ\text{C}
Target P-value, \(P_{\text{target}} = 5.0\ \text{min}
Heat transfer coefficient, \(h = 10000.0\ \text{W/m}^2\text{K}\) (limiting case for constant surface temperature)
Cylinder radius, \(R = 0.035\ \text{m}
Slab half-thickness, \(H = 0.075\ \text{m}
Time step, \(\Delta t = 10.0\ \text{s}
Number of series terms, \(N = 5

Step-by-Step Calculation:

Calculate Biot numbers to confirm the constant-surface-temperature limit: \(Bi_{\text{cyl}} = \frac{h R}{k} = \frac{10000.0 \times 0.035}{0.6} = 583.333\) \(Bi_{\text{slab}} = \frac{h H}{k} = \frac{10000.0 \times 0.075}{0.6} = 1250.0\) Since both Bi ≫ 40, the problem is in the constant‑surface‑temperature regime; we use the Bi → ∞ Heisler series.
Compute Fourier numbers at the elapsed time \(t = 3600.0\ \text{s}\): \(Fo_{\text{cyl}} = \frac{\alpha t}{R^2} = \frac{1.4354 \times 10^{-7} \times 3600.0}{(0.035)^2} = 0.421834\) \(Fo_{\text{slab}} = \frac{\alpha t}{H^2} = \frac{1.4354 \times 10^{-7} \times 3600.0}{(0.075)^2} = 0.091866\)
Evaluate dimensionless center temperatures using 5-term series for slab and cylinder, then take the product-of-solutions. For the slab: \(\theta_{\text{slab}} = 0.960700\) For the cylinder: \(\theta_{\text{cyl}} = 0.139690\) The combined dimensionless temperature: \(\theta_{\text{total}} = \theta_{\text{cyl}} \times \theta_{\text{slab}} = 0.139690 \times 0.960700 = 0.134201\)
Obtain the centre temperature at \(t = 3600\ \text{s}\): \(T_{\text{center}} = T_{\text{water}} + (T_{\text{init}} - T_{\text{water}}) \, \theta_{\text{total}} = 85.0 + (20.0 - 85.0) \times 0.134201 = 76.277\ ^\circ\text{C}\)
Compute the lethal rate at this temperature: \(L = 10^{(T_{\text{center}} - 85.0) / z} = 10^{(76.277 - 85.0) / 7.0} = 0.056735\)
Calculate the cumulative P-value after the full 3600 s (i.e., after 360 time steps of 10 s). Each step contributes \(L_i \times (\Delta t / 60\ \text{min/s})\). The accumulated value from the simulation is: \(P_{\text{cum}} = 0.368\ \text{minutes}\)
Compare with the target: \(P_{\text{cum}} = 0.368\ \text{min} < P_{\text{target}} = 5.0\ \text{min}\). Therefore, the required residence time is not reached within 3600 s. The condition to exit the loop (\(P_{\text{cum}} \ge 5.0\)) was never satisfied, so the returned residence time is \(0.0\) (meaning the target was not met within the simulated 60 minutes).

Final Answer:
With the given parameters (constant‑surface‑temperature limit, container dimensions, and thermal properties), the centre temperature rises sufficiently to produce a cumulative lethal rate of only \(0.368\ \text{minutes}\) after \(60\ \text{minutes}\). Because this is below the required \(5.0\ \text{minutes}\), the pasteurisation target is not achieved within the first hour. The required residence time exceeds \(60\ \text{minutes}\) (the numerical result reports residence_time as \(0.0\) because the loop never broke within the simulation window).

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle