Critical Freezing Rate Determination

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

Introduction & Context

The Critical Freezing Rate Determination is a fundamental calculation in food process engineering and cryogenics. It is used to predict the time required for a product, modeled as a one‑dimensional infinite slab, to transition from its initial freezing point to a target core temperature. This process is critical for maintaining product quality, as rapid freezing minimizes ice crystal size, thereby preserving cellular structure and texture in biological materials like fish fillets or meat products, and underscores the importance of selecting proper packaging for frozen foods to protect quality during storage and distribution.

This calculation is typically employed during the design phase of industrial blast freezers and contact freezing systems to ensure that the cooling capacity is sufficient to meet safety and quality standards within a specified time frame.

Methodology & Formulas

The determination of the total freezing time is divided into two distinct physical phases: the latent heat removal phase (front arrival) and the sensible cooling phase of the frozen core.

1. Biot Number

The Biot number determines the ratio of internal conductive resistance to external convective resistance. It dictates whether the system can be treated as a lumped mass or if internal temperature gradients must be considered.

\[ Bi = \frac{h \cdot L}{k_{\text{frozen}}} \]

2. Freezing Front Arrival Time (Plank's Equation)

Based on Plank's equation for an infinite slab, this phase accounts for the time required for the freezing front to propagate from the surface to the thermal center, primarily driven by the removal of latent heat of fusion.

\[ t_{\text{front}} = \frac{\rho \cdot h_{sf}}{T_{f} - T_{\infty}} \cdot \left( \frac{L}{h} + \frac{L^{2}}{2 \cdot k_{\text{frozen}}} \right) \]

3. Sensible Cooling Time (First-Term Series Solution)

Once the freezing front reaches the center, the core undergoes sensible cooling from the freezing point to the target temperature. This is calculated using the first-term series solution for a slab with convective boundary conditions, requiring the determination of the first eigenvalue \(\zeta_{1}\) and the coefficient \(C_{1}\).

\[ t_{\text{sens}} = \frac{L^{2}}{\alpha_{\text{frozen}} \cdot \zeta_{1}^{2}} \cdot \ln\left( \frac{C_{1} \cdot (T_{f} - T_{\infty})}{T_{c} - T_{\infty}} \right) \]

Where the thermal diffusivity of the frozen product is defined as:

\[ \alpha_{\text{frozen}} = \frac{k_{\text{frozen}}}{\rho \cdot c_{p,\text{frozen}}} \]

The eigenvalue \(\zeta_{1}\) is the first positive root of the transcendental equation for a plane wall:

\[ \zeta_{1} \cdot \tan(\zeta_{1}) = Bi \]

The first-term series coefficient \(C_{1}\) is calculated as:

\[ C_{1} = \frac{2 \cdot \sin(\zeta_{1})}{\zeta_{1} + \sin(\zeta_{1}) \cdot \cos(\zeta_{1})} \]

4. Total Zone Passage Time and Cooling Rate

The total time to pass through the critical freezing zone is the sum of the front arrival time and the sensible cooling time.

\[ t_{\text{zone}} = t_{\text{front}} + t_{\text{sens}} \]

The average cooling rate across the critical temperature zone (in K/min) is:

\[ \text{Cooling Rate} = \frac{T_{f} - T_{c}}{t_{\text{zone}} \, / \, 60} \]

Regime Validity and Constraints

Parameter	Condition	Implication
Biot Number (\(Bi\))	\(Bi \lt 0.1\)	Internal gradients negligible; lumped capacitance model is sufficient.
Biot Number (\(Bi\))	\(0.1 \le Bi \le 10\)	Distributed-parameter model required. Plank's equation and the convective series solution are valid.
Biot Number (\(Bi\))	\(Bi \gt 10\)	External resistance negligible; constant surface-temperature (\(T \approx T_\infty\)) assumption applies.
Fourier Number (\(Fo\))	\(Fo \lt 0.2\)	First-term series approximation is insufficient; higher-order terms or numerical methods required.
Fourier Number (\(Fo\))	\(Fo \ge 0.2\)	First-term series approximation is valid (error < 2%).

To determine the critical freezing rate, process engineers follow a structured thermal analysis of the product:

Characterize thermal properties: Measure or obtain the frozen-state density (\(\rho\)), thermal conductivity (\(k_{\text{frozen}}\)), specific heat capacity (\(c_{p,\text{frozen}}\)), and latent heat of fusion (\(h_{sf}\)) for the specific formulation.
Define the geometry: Model the product as an infinite slab and determine the half-thickness (\(L\)), which is the characteristic length for symmetric heat transfer from both sides.
Establish process conditions: Specify the freezer medium temperature (\(T_{\infty}\)) and the convective heat transfer coefficient (\(h\)) for the air-blast or contact freezing system.
Identify the critical temperature zone: The zone of maximum ice crystal formation (typically \(-1^{\circ}\mathrm{C}\) to \(-5^{\circ}\mathrm{C}\) for many food products) defines the freezing point (\(T_{f}\)) and target core temperature (\(T_{c}\)).
Calculate the Biot number: Verify that \(0.1 \le Bi \le 10\) so the distributed-parameter model (Plank's equation + series solution) is applicable.
Compute freezing times: Use Plank's equation for the front arrival time and the first-term series solution for the sensible cooling time, then sum them to obtain the total zone passage time (\(t_{\text{zone}}\)).
Compute the rate: Divide the critical temperature interval (\(T_{f} - T_{c}\)) by the zone passage time (in minutes) to obtain the average cooling rate in K/min.

Suboptimal (excessively slow) freezing rates manifest through several physical and quality defects in the final frozen product:

Large, extracellular ice crystals: Slow freezing allows water to migrate out of cells and form large ice crystals in the extracellular spaces, visible under microscopy. These large crystals physically rupture cell walls and membranes.
Excessive drip loss upon thawing: Ruptured cells cannot reabsorb melted water, leading to high drip (purge) loss—a direct indicator of structural damage during freezing.
Degraded texture: The product becomes mushy, spongy, or fibrous after thawing due to the loss of cellular turgor and structural integrity.
Freezer burn and surface dehydration: Prolonged residence time in the freezer increases moisture loss from the product surface, causing oxidative rancidity and dry, leathery patches.
Non-uniform color and appearance: Slow freezing can cause solute concentration effects that alter the optical properties of the tissue, leading to a dull or darkened appearance.

The cooling rate through the critical freezing zone (\(-1^{\circ}\mathrm{C}\) to \(-5^{\circ}\mathrm{C}\)) directly governs the size, morphology, and distribution of ice crystals, which in turn determines the final product quality:

Fast freezing (e.g., > 1 K/min): Produces numerous small, intracellular ice crystals. This preserves cellular structure, minimizes drip loss on thawing, and maintains a firm, natural texture. It is essential for high-value products such as fish fillets, meat cuts, and individually quick-frozen (IQF) vegetables.
Moderate freezing (e.g., 0.2–1 K/min): Results in a mixture of intra- and extracellular ice. Some cellular damage may occur, with moderate drip loss and acceptable texture for many commodity products.
Slow freezing (e.g., < 0.2 K/min): Leads to large, predominantly extracellular ice crystals that cause extensive cell wall rupture. The result is high drip loss, poor texture, and degraded sensory properties. Slow freezing is generally unacceptable for cellular biological materials.
Critical zone importance: The majority of ice formation occurs between \(-1^{\circ}\mathrm{C}\) and \(-5^{\circ}\mathrm{C}\). The residence time in this zone is the single most important determinant of ice crystal size distribution and, consequently, product quality.
Nutritional and sensory retention: Rapid freezing better preserves volatile aroma compounds, water-soluble vitamins, and natural pigments by minimizing the time available for oxidative and enzymatic degradation.

Worked Example: Critical Freezing Rate for a Fish Fillet

Consider a 20 mm thick fish fillet (half-thickness \(L = 0.01\) m) being frozen in an air-blast freezer at \(T_{\infty} = -40.0^{\circ}\mathrm{C}\) with a convective heat transfer coefficient \(h = 100.0\) W/m²K. The target is to pass the thermal center from the initial freezing point \(T_{f} = -1.0^{\circ}\mathrm{C}\) to \(T_{c} = -5.0^{\circ}\mathrm{C}\) (a \(\Delta T_{\text{sens}} = 4.0^{\circ}\mathrm{C}\) critical zone) in no more than 20 minutes. The frozen-state properties are: density \(\rho = 1000.0\) kg/m³, thermal conductivity \(k_{\text{frozen}} = 1.5\) W/mK, specific heat \(c_{p,\text{frozen}} = 2000.0\) J/kgK, and latent heat of fusion \(h_{sf} = 300,\!000\) J/kg. The thermal diffusivity is \(\alpha_{\text{frozen}} = k_{\text{frozen}} \,/\, (\rho \cdot c_{p,\text{frozen}}) = 7.50 \times 10^{-7}\) m²/s.

Known 1: Half-thickness, \(L = 0.01\) m
Known 2: Freezer temperature, \(T_{\infty} = -40.0^{\circ}\mathrm{C}\)
Known 3: Convective coefficient, \(h = 100.0\) W/m²K
Known 4: Initial freezing point, \(T_{f} = -1.0^{\circ}\mathrm{C}\)
Known 5: Target center temperature, \(T_{c} = -5.0^{\circ}\mathrm{C}\)
Known 6: Density, \(\rho = 1000.0\) kg/m³
Known 7: Frozen thermal conductivity, \(k_{\text{frozen}} = 1.5\) W/mK
Known 8: Frozen specific heat, \(c_{p,\text{frozen}} = 2000.0\) J/kgK
Known 9: Latent heat of fusion, \(h_{sf} = 300,\!000\) J/kg
Known 10: Thermal diffusivity, \(\alpha_{\text{frozen}} = 7.50 \times 10^{-7}\) m²/s

Step 1: Biot Number \[ Bi = \frac{h \cdot L}{k_{\text{frozen}}} = \frac{100.0 \cdot 0.01}{1.5} = 0.667 \] Since \(0.1 \le Bi \le 10\), the lumped capacitance assumption is invalid; the combined front-propagation and series solution methodology is appropriate.
Step 2: Freezing Front Arrival Time (Plank's Equation) \[ t_{\text{front}} = \frac{\rho \cdot h_{sf}}{T_{f} - T_{\infty}} \left( \frac{L}{h} + \frac{L^{2}}{2 \cdot k_{\text{frozen}}} \right) = \frac{1000.0 \times 300,\!000}{39.0} \left( \frac{0.01}{100.0} + \frac{0.01^{2}}{2 \times 1.5} \right) = 1025.6 \ \text{s} \] The temperature driving force is \(T_{f} - T_{\infty} = (-1.0) - (-40.0) = 39.0^{\circ}\mathrm{C}\).
Step 3: First Eigenvalue \(\zeta_{1}\) for the Slab The eigenvalue satisfies the transcendental equation \(\zeta_{1} \cdot \tan(\zeta_{1}) = Bi\). Using Newton's method with an initial guess \(\sqrt{Bi} = 0.8165\): \[ \zeta_{1} = 0.7360 \] (Converged after iteration; the exact computed value is \(\zeta_{1} = 0.736005193\).)
Step 4: Coefficient \(C_{1}\) for the First-Term Series \[ C_{1} = \frac{2 \cdot \sin(\zeta_{1})}{\zeta_{1} + \sin(\zeta_{1}) \cdot \cos(\zeta_{1})} = \frac{2 \times 0.671332}{0.736005 + 0.671332 \times 0.741156} = 1.088 \]
Step 5: Sensible Cooling Time (First-Term Series) First compute the argument of the logarithm: \[ \text{arg}_{\ln} = \frac{C_{1} \cdot (T_{f} - T_{\infty})}{T_{c} - T_{\infty}} = \frac{1.088 \times 39.0}{35.0} = 1.213 \] Then: \[ t_{\text{sens}} = \frac{L^{2}}{\alpha_{\text{frozen}} \cdot \zeta_{1}^{2}} \cdot \ln(\text{arg}_{\ln}) = \frac{0.01^{2}}{7.50 \times 10^{-7} \times (0.7360)^{2}} \cdot \ln(1.213) = 47.50 \ \text{s} \] The Fourier number for the sensible cooling phase: \[ Fo_{\text{sens}} = \frac{\alpha_{\text{frozen}} \cdot t_{\text{sens}}}{L^{2}} = \frac{7.50 \times 10^{-7} \times 47.50}{0.01^{2}} = 0.356 \] Since \(Fo_{\text{sens}} = 0.356 \ge 0.2\), the first-term approximation is valid (error < 2%).
Step 6: Total Zone Passage Time \[ t_{\text{zone, sec}} = t_{\text{front}} + t_{\text{sens}} = 1025.6 + 47.50 = 1073.1 \ \text{s} \] Convert to minutes: \[ t_{\text{zone, min}} = \frac{1073.1}{60} = 17.89 \ \text{min} \] This is below the 20-minute threshold, so the freezer design is acceptable.
Step 7: Average Cooling Rate Through the Critical Zone \[ \text{Cooling Rate} = \frac{T_{f} - T_{c}}{t_{\text{zone, min}}} = \frac{4.0^{\circ}\mathrm{C}}{17.89 \ \text{min}} = 0.224 \ \text{K/min} \]

Final Answer: The total time for the thermal center to pass from \(T_{f} = -1.0^{\circ}\mathrm{C}\) to \(T_{c} = -5.0^{\circ}\mathrm{C}\) is 17.89 minutes, yielding an average cooling rate of 0.224 K/min. The Biot number (\(Bi = 0.667\)) and Fourier number (\(Fo = 0.356\)) both lie within the required validity ranges. The design satisfies the 20-minute criterion.

"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