Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Freezing Time

Download as pdf or txt
Download as pdf or txt
You are on page 1of 14

Freezing Loads and Freezing

6 Time Calculation
Gauri S. Mittal
University of Guelph, Guelph, Ontario, Canada

CONTENTS
I. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
II. Freezing Load . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
A. Calculation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
B. Freezing Rate and Thermal Center . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
III. Freezing Time or Rate Prediction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
A. Planks Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
B. Nagaoka et al. Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
C. Levy Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
D. Cleland and Earle Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
E. Cleland et al. Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
F. Pham Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
G. Modified Pham Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
IV. Thawing Time Prediction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
IV. Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Nomenclature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139

I. INTRODUCTION
The purposes of food freezing are: (i) preservation of food; (ii) reducing the activity of enzymes and
microorganisms; (iii) reducing the amount of liquid water for microbial growth; and (iv) reducing
water activity (aw) of foods. Many types of freezers are used for this purpose. Some of these are: (i)
air blast freezers, batch, or continuous; (ii) still air freezers; (iii) belt freezers; (iv) spiral belt free-
zers; (v) fluidized bed freezers; (vi) plate freezers a series of flat plates kept cool by circulating a
coolant; (vii) liquid immersion freezers chilled brine or glycol is used, can also be sprayed; and
(viii) cryogenic freezers liquid N2 or liquid CO2 is used (the boiling point for N2 is 21968C and
for CO2 is 2798C).
Product quality is influenced by ice-crystal size and configuration during the freezing operation.
The advantages of fast freezing operation can be lost during the storage because of the formation of
large ice crystals by joining small crystals. Hence, complete product freezing in the freezer is more
important [1]. The process of ice-crystal formation is based on two operations: (1) nucleation or
crystal formation it influences the type of crystal structure formed in a food product and ice-
crystal nucleation is created by supercooling below initial freezing point, similar to crystallization
process; (2) rate of crystal growth is also supercooling-driven, which depends on (i) diffusion
rate of water molecules from the unfrozen solution to the crystal surface, (ii) the rate at which
heat is removed, and (iii) temperature of the solution.

127
2006 by Taylor & Francis Group, LLC
128 Freezing Loads and Freezing Time Calculation

II. FREEZING LOAD


A. CALCULATION
Freezing load or enthalpy change (DH) to reduce the product temperature (Ti) from some level
above the freezing point (TF) to some desired final temperature (T) is given by

DH sensible heat removed from the product solids (DHs )


sensible heat removed from unfrozen water (DHu )
enthalpy change due to latent heat (DHL )
sensible heat removed from the frozen water (DHI )
T
DHs Ms Cps (Ti  TF ) Ms Cps dT (6:1a)
TF
0
DHu Mu Cpu (Ti  TF ) Mu Cpu (TF  T) (6:1b)
DHL MI LV (6:1c)
T
DHI MI CpI (TF  T) or MI CpI dT (6:1d)
TF

where Ms is the mass of solids, Cps the specific heat of solids, Mu the mass of unfrozen water, Cpu
the specific heat of unfrozen water, MI the mass of ice or unfrozen water, Lv the latent heat of freez-
0
ing, CpI the specific heat of ice, and Cpu the specific heat of unfrozen water below TF. Enthalpy com-
position charts for different food materials using experimental data were provided [2,3]. One
example is given in Figure 6.1.

B. FREEZING RATE AND THERMAL CENTER


The absence of a consistent definition for the freezing time is one of the problems associated with
the published literature on the freezing of foods. This problem arises mainly because foods do not
freeze at a distinct temperature, but rather the phase change takes place over a range of tempera-
tures. A definition of the freezing time requires a definition of the freezing point. A variable temp-
erature distribution exists within the food product during the freezing process, giving different
freezing times depending on the point within the product where the temperature is monitored.
The thermal center, defined as the location in the material which cools most slowly, is generally
used as the reference location. The effective freezing time, defined by the International Institute of

FIGURE 6.1 Riedel plot for grape juice. (From L Riedel, Kaltetechnik 9:38 40, 1957.)

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 129

Refrigeration [4], is the total time required to lower the temperature of a food material at its thermal
center to a desired temperature below the initial freezing point. Other definitions are:

1. The time required to reduce the product temperature at the slowest cooling location from
the initial freezing point to some desired and specified temperature below the initial freez-
ing point.
2. International Institute of Refrigeration [4] definition: It is the ratio between the minimum
distance from the surface to the thermal center and the time elapsed between the surface
reaching 08C and the thermal center reaching 58C colder than the temperature of initial ice
formation at the thermal center (cm/h).

III. FREEZING TIME OR RATE PREDICTION


It is important to accurately predict the freezing times of foods to assess the quality, processing
requirements, and economical aspects of food freezing. A number of models have been proposed
in the literature to predict freezing times. However, as the freezing process is a moving boundary
problem, that is, one involving a phase change, most of the single-phase, unsteady-state solutions
are unsuitable.
Foods, undergoing freezing, release latent heat over a range of temperatures. Freezing does not
occur at a unique temperature. In addition, foods do not have constant thermal properties during
freezing [5]. As a result, no exact mathematical model exists for predicting the freezing of
foods. Researchers, who have found a solution, have either used numerical finite difference or
finite element methods. So, models for predicting freezing times range from approximate analytical
solutions to more complex numerical methods.
In the past, an extensive amount of work has been done to develop mathematical models for
the prediction of food freezing times. The accuracy of such models is dependent on how closely
the corresponding assumptions approach reality. Most of these models are usually categorized
into one of two forms, analytical or numerical, with the latter generally considered as more
accurate due to the inclusion of a set of assumptions and boundary conditions, which are of a
more realistic nature than those pertaining to the former. Approximately, 30 different methods
to predict freezing and thawing times were reviewed [6]. Details on these models are given
elsewhere [7,8].
The general approach of researchers in the food-freezing field has been to seek approximate or
empirical relationships, rather than to try to derive exact analytical equations. The method can be
classified into two groups: (1) methods relying on analytical approximations, such as those of Refs.
[9 13] or (2) methods relying on regression of computer results or experimental data, such as those
of Refs. [14 18]. The methods vary considerably in complexity and accuracy, the number of arbi-
trary or empirical parameters used ranging from 0 to more than 50 [19].

A. PLANKS EQUATION
Planks equation was derived based on energy balance principle [9]. Heat condition through frozen
region is written as: (Figure 6.2)
 
TS  TF
q kI A (6:2)
X

Convective heat transfer at the surface is given by:


q hc A (T1  TS ) (6:3)

2006 by Taylor & Francis Group, LLC


130 Freezing Loads and Freezing Time Calculation

FIGURE 6.2 Diagram to derive Planks equation.

Total resistance

X 1
Rt (6:4)
k I A hc A

or
DT T1  TF
q(overall) (6:5)
Rt (X=kI A 1=hc A)

This heat transfer should be equal to the latent heat of freezing or

dX
qA rL V (6:6)
dt
dX
the velocity of the freezing front (6:7)
dt

or

dX (T1  TF )A
A rLV  (negative heat transfer) (6:8)
dt (X=kI 1=hc )

or
tF a=2  
LV r 1 X
dt  dX (6:9)
0 T 1  TF 0 hc k I

or
   
LV r a a2 LV r a a2
tF  (6:10)
T1  TF 2hc 8kI TF  T1 2hc 8kI

General form:
 
rLV Pa Ra2
tF (6:11)
TF  T1 hc kI

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 131

where TF is the initial freezing point of the product, TS the surface temperature, kI the thermal con-
ductivity of frozen food, X the thickness of frozen food, hc the convective heat transfer coefficient, A
the surface area, T1 the ambient temperature, Lv the latent heat of freezing, and r the food density.
P and R values for different shaped foods are:

Infinite slab Infinite cylinder Sphere


P 1/2 1/4 1/6
R 1/8 1/16 1/24

For brick-shaped material, Figure 6.3 provides P and R for different b1 and b2 values [20].

Example 1: Lean beef block with dimensions of 1 m  0:25 m  0:6 m, hc 30 W=(m2 K),
T0 58C, T 108C, T1 308C, r 1050 kg/m3 , LV 333:22 kJ/kg, m.c. 74:5%,
kI 1:108 (W/m K), TF 1:758C. Find freezing time using Planks equation.

10
9 0.15

8 0.10
40
7
0.95
2 6
0.9
R
38
5 0.85

0.8
36
4
0.75
34
3 0.7

32
0.65

30 P
0.6
2
28
0.55
26
0.5 24
22
0.45 20
18
1
1 2 3 4 5 6 7 8 9
1

FIGURE 6.3 P and R values for different b1 and b2. (From AJ Ede. Modern Refrigeration 52:52-55, 1949.)

2006 by Taylor & Francis Group, LLC


132 Freezing Loads and Freezing Time Calculation

Solution:

0:6 1
b1 2:4, b2 4, [ P 0:3, R 0:085
0:25 0:25
 
(1050)(333:22  0:745)(1000 J/kJ) 0:3(0:25) 0:085(0:025)2
t 18:7 h
1:75  (30)3600 s/h 30 1:108

The limitations of Planks equation are as follows:

1. It neglects the time required to remove sensible heat above the initial freezing point.
2. It does not consider the gradual removal of latent heat over a range of temperatures during
the freezing process.
3. Constant thermal conductivity assumed for frozen material.
4. It assumes the product to be completely in liquid phase.

Many modifications were suggested on Planks equation to improve its accuracy. Some of these
are given in the subsequent sections.

B. NAGAOKA ET AL. EQUATION


Nagaoka et al. [21] proposed the modifications of Planks equation as follows:
 
DH 0 r Pa Pa2
tF (6:12)
TF  T1 hc kI
DH 0 (1 0:008Ti )Cpu (Ti  TF ) LV CPI (TF  T) (6:13)

where Ti is the initial food temperature, T the final frozen food temperature, Cpu the specific heat of
unfrozen food, and CPI the specific heat of frozen food.

C. LEVY EQUATION
Levy [22] considered the following definition of enthalpy to modify Planks equation:

DH 0 (1 0:008(Ti  TF ))Cpu (Ti  TF ) LV CPI (TF  T) (6:14)

Example 2: Use modified Planks equation to calculate the freezing time for the lean beef block of
1m  0:6m  0:25m, hc 30W=(m2 K), T0 58C, T 108C, T1 308C, r 1050 kg/m3 ,
TF 1:758C, tF ?

DH 333:22 kJ=kg (0:745 m.c.) 248:25 kJ=kg;

Solution:

0:6 1
b2 2:4, b1 4, [ P 0:3, R 0:085
0:25 0:25
Cpu 3:52 kJ=kg K, CPI 2:05 kJ=kg K, kI 1:108 W=m K

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 133

DH 0 {1 0:0085  (1:75)}3:52(5  (1:75)) 248:25 2:05(1:75  (10))


297:59 kJ=kg

 
(1050)(297:59)(1000) 0:3(0:25) 0:085(0:025)2
tf 22:41 h
1:75  (30)3600 30 1:108

D. CLELAND AND EARLE EQUATION


Cleland and Earle [23] modified Planks equation using the nondimensional numbers as
follows:

CPI (TF  T1 )
NSte Stefan number (6:15)
DHref

   
rDHref Pa Ra2 1:65NSte T  T1
tF 1 ln (6:16)
E(TF  T1 ) hc kI kI Tref  T1

Cpu (Ti  TF )
NPK Planks number (6:17)
DHref

P 0:51:026 0:5808NPK NSte (0:2296NPK 0:105) (6:18)

R 0:1251:202 NSte (3:410NPK 0:7336) (6:19)

where Tref is the reference temperature and E is 1 for an infinite slab, 2 for an infinite cylin-
der, and 3 for a sphere. Tref is taken as 2108C and DHref is enthalpy change from TF to Tref.
0:15  NSte  0:35, 0:2  NBi  20, and 0  NPK  0:55.

Example 3 (Cleland and Earle [23] approach): Lamb steak (slab) 0:025 m thick, Ti 208C, T
108C, T1 308C, r 1050 kg/m3 , TF 2:758C, kI 1:35 W/m K, hc
20 W=(m2 K), E 1 for slab, tF ? Cpu 3 kJ=kg K, CPI 1:75 kJ=kg K, DH 240 kJ=kg.
Solution:

DHref 240 1:75(2:75 10) 252:7 kJ/kg

CPI (TF  T1 ) 1:75(2:75 30)


NSte 0:189
DHref 252:7
Cpu (Ti  TF ) 3(20 2:75)
NPK 0:270
DHref 252:7

P 0:51:026 0:5808(0:270) 0:189(0:2296(0:270) 0:105) 0:607

R 0:1251:202 0:189(3:410(0:270) 0:7336) 0:189

2006 by Taylor & Francis Group, LLC


134 Freezing Loads and Freezing Time Calculation

   
rDHref Pa Ra2 1:65 NSte T  T1
tF 1 ln
TF  T1 hc kI kI Tref  T1
 
252:7(1000)(1050) 0:607(0:025) 0:189(0:025)2
(1) 2:289 h
(2:75 30)3600 20 1:35

E. CLELAND ET AL. METHOD


Cleland et al. [24,25] method is based on Calvelo [26] approach, which is given below:
 
1:3179rCPI a2 0:5 0:125 0:9576 0:0550 0:0017NBi 0:1727NPK
tf NSte 10
kI E NBi NSte NSte
  (6:20)
1:65NSte T  T1
 1 ln
kI Tref  T1

Tref is also 108C. NBi is given by ha=kI .

F. PHAM METHOD
Pham method [13] involves total of precooling, phase change, and tempering times.

1 X3
(1 NBii =ai )
tf DHi a (6:21)
E i1 2DTi hc

where
DH1 Cpu (Ti  Tf, ave ) (6:22)
(Ti  T1 )  (Tf;ave  T1 )
DT1 , a1 6 (6:23)
ln ((Ti  T1 )=(Tf;ave  T1 ))
 
hc a hc a
NBii 0:5 (6:24)
kI ku
DH3 CPI (Tf;ave  Tave ), NBi3 NBi2 (6:25)
(Tf;ave  T1 )  (Tave  T1 )
DT3 (6:26)
ln ((Tf;ave  T1 )=(Tave  T1 ))
hc a
DH2 DHf , DT2 Tf;ave  T1 , NBi2 , a2 4 (6:27)
kI
T  T1
Tave T  , a3 6, Tf;ave TF  1:5 (6:28)
2 4=NBi3

where kW is the thermal conductivity of unfrozen food.

G. MODIFIED PHAM METHOD


This modified method of Pham [19] was given after summing precooling, phase change, and
tempering times. E is given by the literature [27]. This method is to calculate the freezing and
thawing time for finite size objects of any shape by approximating them to be similar to an ellipsoid.
The following assumptions were used in developing this method: (i) uniform initial product

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 135

temperature, Ti; (ii) uniform and constant ambient conditions; (iii) a fixed value of final
product temperature, T; and (iv) convective surface heat transfer is following Newtons law of
cooling.
For infinite slab, the freezing time (tslab) is given by

  
ra DH1 DH2 NBi
tslab 1 (6:29)
2hc DT1 DT2 4

Equation (6.21) is valid for the following ranges: 0:02 , NBi , 11, 0:11 , NSte , 0:36, and
0:03 , NPK , 0:61.
The thawing time is given by for thawing to Tf 08C:

 
1:4921Cpu a2 0:5 0:125 1:0248 0:2712 0:061
tslab NSte NPK (6:30)
ku NBi NSte NSte

Equation (6.30) is valid for the following ranges: 0:3 , NBi , 41, 0:08 , NSte , 0:77, and
0:06 , NPK , 0:27.
where

DH1 Cpu (Ti  T3 ) (6:31)


DH2 DH CPI (T3  T) (6:32)
T i  T3
DT1  T1 , DT2 T3  T1 , DT3 1:8 0:263T 0:105T1 (6:33)
2
NBi hc a=kI (6:34)
CPI (TF  T1 )
NSte (6:35)
DHref
Cpu (Ti  TF )
NPK (6:36)
DHref

where Cpu is the specific heat of unfrozen product (J/(kg K)), CPI the specific heat of the frozen
product (J/(kg K)), hc the convective heat transfer coefficient (W/(m2 K)), T1 the ambient temp-
erature (8C), TF the initial freezing temperature (8C), r the product density (kg/m3), kI the thermal
conductivity of frozen product (W/(m K)), ku the thermal conductivity of unfrozen product
(W/(m K)), DH the enthalpy change due to freezing (moisture content) (333 220) (J/kg), R
the characteristic dimension (m), that is radius of cylinder of sphere or half thickness of slab or
other geometries, NBi the Biot number, NSte the Stefan number, NPK the Plank number, and T
the final product temperature (8C).
For other shapes than infinite slab, the following modification is used:

tslab
tellipsoid (6:37)
E
1 2=NBi 1 2=NBi
E 1 (6:38)
b21 2b1 =NBi b22 2b2 =NBi

where V is the volume (m3) and A the smallest cross-sectional area that incorporate R (m2).

2006 by Taylor & Francis Group, LLC


136 Freezing Loads and Freezing Time Calculation

For an infinite slab, E is 1, for an infinite cylinder E is 2, and for a sphere E is 3.

A V
b1 and b2 (6:39)
p R2 b1 (4=3pR3 )

Notes:

1. For Equations (6.11), (6.12), (6.16), (6.20), (6.21), (6.29), (6.30), and (6.34), a is slab
thickness or diameter of cylinder or sphere, or the smallest dimension of brick-shaped
or dissimilar products.
2. DH (moisture content) (latent heat of fusion); DH 0 for modified Planks equations such
as Levys [22] and Nagaoka et al. [21], and DH10 DH CPI (TF  10).

Example 4: Beef slab of 1 m  0:6 m  0:25 m, R 0:25=2 0:125 m, hc 30W=(m2 K),


r 1050 kg/m3 , Ti 58C, T 158C, T1 308C, CPI 2:5 kJ=(kg K), Cpu 3:52 kJ=(kg K),
moisture content 74.5% wet basis, DH 333:22(0:745) 248:25 kJ/kg, m.c. 74:5%,
kI 1:108(W/mK), TF 1:758C, t ?
Solution:

hc a 30(0:125)
NBi 3:3845
kI 1:108
CPI (TF  T1 ) 2:05(1:75 30)
NSte 0:234
DHref 248:25
Cpu (Ti  TF ) 3:52(5 1:75)
NPK 0:0955
DHref 248:25
T3 1:8 0:263T 0:105T1 1:8 0:263(15) 0:105(30) 5:295

DH1 Cpu (Ti  T3 ) 3520(5 5:295) 36238:4 J/kg

DH2 DH CPI (T3  T) 248250 3520(5:295 15) 282411:6 J/kg

T i  T3 (5 5:295)
DT1  T1  (30) 35:1475
2 2
DT2 T3  T1 5:295 30 24:705

The freezing time (tslab) is given by


     
rR DH1 DH2 NBi 1050(0:125) 36238:4 282411:6 3:3845
tslab 1 1
hc DT1 DT2 2 30 35:1475 24:705 2
146789:45 s 40:775 h

A 0:25(0:6)
b1 3:056
pR2 p(0:125)2
V 0:25(0:6)(1)
b2 4 4 3
6:0
b1 3 pR 3
3 p(3:056)(0:125)

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 137

Actual freezing time is given by

tslab
tellipsoid
E
1 2=NBi 1 2=NBi
E 1
b21 2b1 =NBi b22 2b2 =NBi
1 2=3:3845 1 2=3:3845
1 2 1:4939
3:0562 2(3:056)=3:3845 6 2(6)=3:3845

Therefore, t 40.775/1.4939 27.294 h.

IV. THAWING TIME PREDICTION


Although thawing is the opposite process of freezing, the earlier equations on freezing time predic-
tion cannot be readily applied to thawing process. The thawing time is given for thawing to
TF 08C, and can be calculated by one of the following methods.

1. Power law approach to modifying Planks equation as proposed by Calvelo [26] and
Cleland [28]: This and other methods are valid for the following ranges:
0.6 , NBi , 57.3, 0.08 , NSte , 0.77, and 0.06 , NPK , 0.27.

 
1:4921rCpu a2 0:5 0:125 1:0248 0:2712 0:061
tslab NSte NPK (6:40)
ku NBi NSte NSte

where
hc a
NBi (6:41)
kI
CPI (T1  TF )
NSte (6:42)
DH10
Cpu (TF  Ti )
NPK (6:43)
DH10
Here DH10 is the enthalpy change for the temperature change from 0 to 2108C.

2. Linear correction [23]:

 
rCpu a2 P R
t (6:44)
kI E NBi NSte NSte
P 0:50:7754 2:2828NSte NPK  (6:45)
2
R 0:1250:4271 2:1220NSte  1:4847NSte  (6:46)

3. Three-stage calculation method [13]:

rX 3
(1 ha=4kI )
t DHi a (6:47)
E i1 2DTi hc

2006 by Taylor & Francis Group, LLC


138 Freezing Loads and Freezing Time Calculation

where

DH1 Cpu (Tf;ave  T) (6:48)

(Ti Tf;ave )
DT1 T1  , k1 kI (6:49)
2
DH3 Cpu (Tave  Tf;ave ) (6:50)

(Tave Tf;ave )
DT3 T1  , k3 ku (6:51)
2
DH2 DHf (6:52)

DT2 T1  Tf;ave , k2 0:25kI 0:75ku (6:53)

Tf;ave TF  1:5 (6:54)

(T  T1 )
DTave T  , (6:55)
2 4=NBi

4. Correction of Planks equation [13]:

  
rCpu a2 1 1 0:0244 0:6192NPK
tt 0:8941 
ku E 2NBi NSte 8NSte NSte NBi
 
Cpu (Tave  T)
 1 (6:56)
DH10

IV. CONCLUSIONS
Many equations and models have been suggested to calculate freezing time of foods. Whenever a
freezing time prediction method is used, some inaccuracy will be inevitable. This may arise from
one of the three sources: (a) inaccuracy in thermal data; (b) inaccurate knowledge of freezing con-
ditions, particularly the surface heat transfer coefficient; and (c) inaccuracy arising from assump-
tions made in the derivation of the prediction equation. The best freezing time prediction
method will be the one in which the error arising from the category (c) is the least. The method
should require as few input data as possible, and preferably should avoid lengthy or complex oper-
ations or reference to grasp and table. Three important parameters affecting the freezing time pre-
diction are Lv, hc, and D. The parameter hc is the most difficult one to measure accurately, and
therefore, is a major source of error.

NOMENCLATURE
A smallest cross-sectional area that incorporate R (m2)
A surface area (m2)
CPI specific heat of the frozen product (J/(kg K))
Cps specific heat of solids (J/(kgK))
0
Cpu specific heat of unfrozen water below TF (J/(kg K))
Cpu specific heat of unfrozen product (J/(kg K))

2006 by Taylor & Francis Group, LLC


Handbook of Frozen Food Processing and Packaging 139

E shape factor
hc convective heat transfer coefficient (W/(m2 K))
DH enthalpy change due to freezing (moisture content)
(333220 J/kg)
DH freezing load or enthalpy change (J/kg)
D HI sensible heat removed from the frozen water (J/kg)
DHL enthalpy change due to latent heat (J/kg)
DHs enthalpy change of product solids (J/kg)
DHu sensible heat removed from unfrozen water (J/kg)
DH10 enthalpy change for the temperature change from 0 to
2108C J/kg)
kI thermal conductivity of frozen product (W/(m, K))
ku thermal conductivity of unfrozen product (W/(m, K))
kW thermal conductivity of unfrozen food (W/(m, K))
Lv latent heat of freezing (J/kg)
MI mass of ice or unfrozen water (kg)
MS mass of solids (kg)
Mu mass of unfrozen water (kg)
NBi Biot number
NPK Plank number
NSte Stefan number
R characteristic dimension (m), that is radius of cylinder
of sphere or half thickness, of slab or other geometries
tslab freezing time (s)
T final frozen product temperature (8C)
TF initial freezing temperature (8C)
Ti initial food temperature (8C)
Tref reference temperature (8C)
TS surface temperature (8C)
T1 ambient temperature (8C)
V volume (m3)
X thickness of frozen food (m)
r product density (kg/m3)

REFERENCES
1. DJ Cleland, KJ Valentas. Prediction of freezing time and design of food freezers. In: KJ Valentas,
E Rotstein, RP Singh, Eds. Handbook of Food Engineering Practice. Boca Raton, FL: CRC Press,
1997, pp. 71 124.
2. L Riedel. Calorimetric investigations of the freezing of fish meat. Kaltetechnik 8 (12):374 377, 1956.
3. L Riedel. Calorimetric investigations of the meat freezing process. Kaltetechnik 9:38 40, 1957.
4. IIR, Recommendations for the Processing and Handling of Frozen Foods. International Institute of
Refrigeration, Paris, France, 1972.
5. EJ Rolfe. The chilling and freezing of foodstuffs. In: N Blakebrough, Ed., Biochemical and Biological
Engineering Science, Vol. 2, London, UK: Academic Press, 1968, pp. 137 208.
6. K Hayakawa. Estimation of heat transfer during freezing or defrosting of food. In: Freezing, Frozen
Storage, and Freeze Drying, Bulletin of the International Institute of Refrigeration. 1:293 301, 1977.
7. YC Hung. Prediction of cooling and freezing times. Food Technology 44 (5):137 153, 1990.
8. HS Ramaswamy, MA Tung. Review on predicting freezing times of foods. Journal of Food Process
Engineering 7:169 203, 1984.

2006 by Taylor & Francis Group, LLC


140 Freezing Loads and Freezing Time Calculation

9. R Plank. Beitrage zur Berechrung und Bewertung der Gefrigeschwindikeit von Lebensmittelm
zeitschrift fur die gesamte kalte Industrie. Beih Rcihe 3 (10):1 16, 1941.
10. AK Fleming. Immersion freezing small meat products. In: Proceedings of the 12th International
Congress of Refrigeration. Madrid, 2:683 694, 1967.
11. RH Mascheroni, A Calvelo. A simplified model for freezing time calculation in foods. Journal of Food
Science 47:1201 1207, 1982.
12. QT Pham. An approximate analytical method for predicting freezing times for rectangular blocks of
food stuffs. International Journal of Refrigeration 8:43 47, 1985.
13. QT Pham. Extension to Planks equation for predicting freezing times of foodstuffs of simple shapes.
International Journal of Refrigeration 7:377 383, 1984.
14. AC Cleland, RL Earle. A comparison of analytical and numerical methods of predicting the Freezing
times of foods. Journal of Food Science 42:1390 1395, 1977.
15. AC Cleland, RL Earle. Predicting freezing times of food in rectangular packages, Journal of Food
Science 44:964 970, 1979.
16. YC Hung, DR Thompson. Freezing time prediction for slab shape foodstuffs by an improved
analytical method. Journal of Food Science 48:555 560, 1983.
17. J Succar, K Hayakawa. Parametric analysis for predicting freezing time of infinitely slab shaped food.
Journal of Food Science 49:468 477, 1984.
18. C Lacroix, F Castaigne. Simple method for freezing time calculations for infinite flat slabs, infinite
cylinders and spheres. Canadian Institution of Food Science and Technology Journal 20:251 259,
1987.
19. QT Pham. Simplified equations for predicting the freezing times of foodstuffs. Journal of Food
Technology 21:209 219, 1986.
20. AJ Ede. The calculation of the freezing and thawing of foodstuffs. Modern Refrigeration 52:52 55,
1949.
21. J Nagaoka, S Takagi, S Hotani. Experiments on the freezing of fish by air blast freezer. Journal of
Tokyo University of Fisheries 42 (1):65 73, 1956.
22. FL Levy. Calculating freezing time of fish in air blast freezers. Journal of Refrigeration 1:55 58,
1958.
23. AC Cleland, R L Earle. Freezing time prediction for different final product temperatures. Journal of
Food Science 49:1230 1232, 1984.
24. DJ Cleland, AC Cleland, RL Earle. Prediction of freezing and thawing times for multidimensional
shapes by simple formulae. I. Regular shapes. International Journal of Refrigeration 10:156 164,
1987.
25. DJ Cleland, AC Cleland, RL Earle. Prediction of freezing and thawing times for multidimensional
shapes by simple formulae. II. Irregular shapes. International Journal of Refrigeration
10:234 240, 1987.
26. A Calvelo. Recent studies on meat freezing. In: R Lawrie, Ed., Developments in Meat Science 2.
London: Applied Science, 1981, pp. 125 158.
27. MdM Hossain, DJ Cleland, AC Cleland. Prediction of freezing and thawing times for foods of
3-dimensional irregular shape by using a semi-analytical geometric factor. International Journal
of Refrigeration 15:241 246, 1992.
28. AC Cleland, Food Refrigeration Processes Analysis, Design and Simulation. New York: Elsevier
Applied Science, 1990, pp. 79 152.

2006 by Taylor & Francis Group, LLC

You might also like