Impact of Drying Conditions on Food Structure and Composition
Reference ID: MET-3550 | Process Engineering Reference Sheets Calculation Guide
Introduction & Context
In food process engineering, the drying of biological materials is a critical unit operation that dictates final product quality, shelf stability, and structural integrity; understanding the fundamentals and objectives of food drying is essential for optimizing moisture removal, modeling drying kinetics, and preventing thermal degradation of heat‑sensitive components.
Methodology & Formulas
The analysis treats the food product surface as a flat plate subjected to external forced convection. The process begins by determining the flow regime using the Reynolds number, which relates inertial forces to viscous forces:
The Nusselt number, representing the ratio of convective to conductive heat transfer, is determined based on the flow regime identified by the Reynolds number:
\[ \mathrm{Nu}_{L} = C \cdot \mathrm{Re}_{L}^{\,m} \cdot \mathrm{Pr}^{1/3} \]
Finally, the convective heat transfer coefficient is derived from the Nusselt number:
\[ h = \frac{\mathrm{Nu}_{L} \cdot k_{\text{air}}}{L} \]
Flow Regime
Condition
Correlation Constant (C)
Exponent (m)
Laminar
\( \mathrm{Re}_{L} \leq 5 \cdot 10^{5} \)
0.664
0.5
Turbulent
\( \mathrm{Re}_{L} > 5 \cdot 10^{5} \)
0.037
0.8
The structural integrity and porosity of food materials are highly sensitive to the kinetics of moisture removal. Process engineers should consider the following factors:
High drying temperatures often lead to rapid surface crust formation, which can trap internal moisture and result in a porous, honeycomb-like internal structure.
Increased airflow velocity enhances the external mass transfer coefficient, which may prevent case hardening but can lead to structural collapse if the matrix is not sufficiently rigid.
The glass transition temperature of the food matrix must be monitored to prevent structural shrinkage during the transition from a rubbery to a glassy state.
To preserve the nutritional profile and bioactive compounds, engineers should optimize the drying environment by:
Implementing vacuum drying to lower the boiling point of water, thereby reducing the thermal load on the product.
Utilizing multi-stage drying profiles where high temperatures are used during the constant rate period and lower temperatures are applied during the falling rate period.
Reducing oxygen exposure through inert gas flushing to prevent oxidative degradation of vitamins and lipids.
The effective moisture diffusivity is not a constant value and typically decreases significantly as the moisture content drops. Engineers should account for this by:
Adjusting the drying model to incorporate concentration-dependent diffusivity variables.
Recognizing that as the material dries, the internal resistance to moisture movement becomes the rate-limiting step.
Accounting for the shrinkage of the material, which reduces the diffusion path length but also decreases the surface area available for evaporation.
The glass transition temperature is a critical parameter for process control in spray drying and fluidized bed operations. Key considerations include:
If the product temperature exceeds the glass transition temperature, the material enters a rubbery state, leading to increased surface stickiness.
Stickiness can cause fouling on chamber walls and agglomeration of particles, which negatively impacts flowability.
Engineers can mitigate these issues by adding drying aids like maltodextrin to increase the glass transition temperature of the feed mixture.
Worked Example: Forced Convection Drying of a Food Product Bed
A flat bed of food particles is being dried using a stream of hot air. The air flows parallel to the bed surface at a free-stream velocity of \( U_{\infty} = 10.0 \, \text{m/s} \). The bed has a characteristic length along the flow direction of \( L = 0.5 \, \text{m} \). Use the given air properties to determine the flow regime, the Nusselt number, and the convective heat transfer coefficient.
Knowns:
Air density: \( \rho_{\text{air}} = 0.995 \, \text{kg/m}^3 \)
Calculate the Reynolds number to characterize the flow regime over the flat plate:
\[
\mathrm{Re}_L = \frac{\rho_{\text{air}} \cdot U_{\infty} \cdot L}{\mu_{\text{air}}}
\]
Using the values:
\[
\mathrm{Re}_L = \frac{0.995 \cdot 10.0 \cdot 0.5}{2.10 \times 10^{-5}} = 236\,905
\]
The Reynolds number is \( \mathrm{Re}_L \approx 2.37 \times 10^5 \).
Determine the flow regime. Since the Reynolds number (\( 2.37 \times 10^5 \)) is less than \( 5 \times 10^5 \), the flow is classified as laminar.
Evaluate the Prandtl number from the air properties (provided as a single value):
\[
\mathrm{Pr} = 0.711
\]
The Prandtl number is within the valid range for flat plate correlations (0.6 ≤ Pr ≤ 60).
Compute the Nusselt number using the laminar correlation for external flow over a flat plate:
\[
\mathrm{Nu}_L = 0.664 \cdot (\mathrm{Re}_L)^{0.5} \cdot (\mathrm{Pr})^{1/3}
\]
Substituting the known values:
\[
\mathrm{Nu}_L = 0.664 \cdot (236\,905)^{0.5} \cdot (0.711)^{1/3} \approx 288
\]
Calculate the convective heat transfer coefficient using the Nusselt number:
\[
h = \frac{\mathrm{Nu}_L \cdot k_{\text{air}}}{L}
\]
Plugging in the numbers:
\[
h = \frac{288 \cdot 0.0297}{0.5} \approx 17.1 \, \text{W/m}^2\text{·K}
\]
Final Answer:
Flow Regime: laminar
Reynolds Number (\(\mathrm{Re}_L\)): \(2.37 \times 10^5\)
Prandtl Number (\(\mathrm{Pr}\)): 0.711
Nusselt Number (\(\mathrm{Nu}_L\)): 288
Convective Heat Transfer Coefficient (\(h\)): 17.1 W/m²·K
"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
