METHOD AND APPARATUS FOR FREEZING CONSUMABLE PRODUCTS

Abstract

An apparatus for preserving consumable products comprising an inner housing arranged within an outer insulated housing, wherein walls of the inner housing define a compartment for receiving consumable products, said walls comprising an inlet wall for inflow of a heat exchange fluid into the compartment, an opposed outlet wall for outflow of a heat exchange fluid out of the compartment, side walls and a base, the side walls and base adjoining the inlet wall to the outlet wall, wherein the inlet wall and outlet wall each include a series of apertures to accommodate a continuous heat exchange fluid flow through the apparatus such that, in operation, consumable products received in the compartment of the inner housing are immersed in the heat exchange fluid to exchange heat with the heat exchange fluid.

Claims

1. An apparatus for preserving consumable products comprising an inner housing arranged within an outer insulated housing, wherein walls of the inner housing define a compartment for receiving consumable products, said walls comprising an inlet wall for inflow of a heat exchange fluid into the compartment, an opposed outlet wall for outflow of a heat exchange fluid out of the compartment, side walls and a base, the side walls and base adjoining the inlet wall to the outlet wall, wherein the inlet wall and outlet wall each include a series of apertures to accommodate a continuous heat exchange fluid flow through the apparatus such that, in operation, consumable products received in the compartment of the inner housing are immersed in the heat exchange fluid to exchange heat with the heat exchange fluid.

2. The apparatus of claim 1, wherein the base includes a series of apertures.

3. The apparatus of claim 1, the apparatus including a structure receivable in the compartment for holding the consumable products, wherein the structure is one or more of a tray, a rack and a basket.

4. The apparatus of claim 3, wherein the structure is suspended from a lid of the apparatus.

5. The apparatus of claim 1, wherein the outer housing comprises: an inlet side corresponding to the inlet wall of the inner housing and defining an inlet space between the inlet side and the inlet wall; and an outlet side corresponding to the outlet wall of the inner housing and defining an outlet space between the outlet side and the outlet wall, wherein the inlet side includes at least one inlet communicating from an outside of the outer housing into the inlet space and the outlet side includes at least one outlet communicating from the outlet space to an outside of the outer housing, and wherein, in operation, said heat exchange fluid is introduced into the apparatus via said at least one inlet and removed from the apparatus via said at least one outlet.

6. The apparatus of claim 5, wherein the inlet space and the outlet space are fluidly connected.

7. The apparatus of claim 5, wherein the at least one inlet and the at least one outlet are 80 mm in diameter.

8. The apparatus of claim 1, wherein, in use, the apparatus is connected to an external refrigeration system whereby the heat exchange fluid exchanges heat with a refrigerant.

9. A method of preserving consumable products in a preservation apparatus using a heat exchange fluid, the method comprising: a. determining the total surface area of an approximated geometry of the consumable products to be frozen; b. performing computational fluid dynamics analysis on the consumable products within the apparatus based on flow constraints due both to the tank geometry and predetermined arrangement of food products in the tank and a predetermined increase in temperature of the heat exchange fluid; c. determining heat transfer coefficients of the consumable products at predetermined product surface temperatures; d. selecting a heat transfer coefficient by rounding down the lowest heat transfer coefficient determined at step (c) to the nearest ten; e. dividing the approximated geometry of the consumable products into a predetermined number of equally distributed volume increments; f. estimating thermal properties of the consumable products; g. calculating, based on conservation of energy analysis using the heat transfer coefficient determined at step (d), the time required for each increment to reach a predetermined final temperature from a predetermined initial temperature for a predetermined heat exchange fluid temperature and a predetermined flow rate of the heat exchange fluid; h. determining an achievable average rate of cooling based on the time calculated at step (f); i. if the determined average rate of cooling that can be achieved is 0.5° C. per minute or more, selecting said predetermined heat exchange fluid temperature, and if the determined average rate of cooling that can be achieved is less than 0.5° C. per minute, selecting a heat exchange fluid temperature sufficiently less than the predetermined heat exchange fluid temperature such that the rate of cooling achievable is at least 0.5° C. per minute; j. arranging the consumable products in the apparatus for preservation; and k. subjecting the consumable products to cooling using the temperature determined at step (i) and cooling the consumable products at a rate of at least 0.5° C. per minute.

10. The method of claim 9, wherein the predetermined flow rate is about 4 kg/s.

11. The method of claim 9, wherein the predetermined final temperature is one of: a. −10° C.; and b. −30° C.

12. (canceled)

13. The method of claim 9, wherein the properties of the consumable products include at least one of: a. specific enthalpy, b. water content; c. fat content; and d. protein content.

14.-16. (canceled)

17. The method of claim 9, wherein the desired amount of time is less than 30 minutes.

18. The method of claim 9, wherein the heat transfer coefficient is one of: a. 130 W/m.sup.2.Math.K; and b. 100 W/m.sup.2.Math.K.

19. (canceled)

20. The method of claim 9, wherein an inlet temperature of the heat exchange fluid is −50° C. and an outlet temperature of the heat exchange fluid is −47° C.

21. The method of claim 9, wherein the determined average rate of cooling is between 0.5° C. per minute and 1.5° C. per minute, and comprising cooling the products at a rate of between 0.5° C. per minute and 1.5° C. per minute.

22. The method of claim 9, wherein the determined average rate of cooling is approximately 1° C. per minute, and comprising cooling the products at a rate of approximately 1° C. per minute.

23. The method of claim 9, wherein the method produces a preserved consumable product such that a quality of at least one of taste, texture and nutrition is retained.

24. (canceled)

25. A method of preserving consumable products using the apparatus of claim 1, wherein the heat exchange fluid is continuously cooled by heat exchange with a refrigerant in a refrigeration cycle external to the apparatus.

Description

BRIEF DESCRIPTION

[0035] Embodiments of the present invention will now be described, by way of non-limiting example only, with reference to the accompanying drawings, in which:

[0036] FIG. 1 is a lower perspective view of a tank for preservation of consumable products;

[0037] FIG. 2 is an upper perspective view of a tank for preservation of consumable products;

[0038] FIG. 3 is an upper perspective view of a tank containing milk bottles;

[0039] FIG. 4 is an upper perspective view of a tank containing trays for fish or meat;

[0040] FIG. 5 is an upper perspective view of a tank containing trays for fish or meat;

[0041] FIG. 6 is an upper perspective view of a tank containing an empty basket for holding bottles;

[0042] FIG. 7 is an upper perspective view of a basket for bottles;

[0043] FIG. 8 is an upper perspective view of a rack for holding trays of meat or fish;

[0044] FIG. 9 is an individual tray for holding meat or fish;

[0045] FIG. 10 is an image showing computational fluid dynamics analysis of temperature in a tank containing 100 1 L bottles of milk;

[0046] FIG. 11 is an image showing computational fluid dynamics of fluid velocity in a tank;

[0047] FIG. 12 is conductivity top view of the analysis shown in FIG. 10;

[0048] FIG. 13 is a cut plot between rows of bottles showing temperature distribution in the tank;

[0049] FIG. 14 is a cut plot through a row of bottles showing temperature distribution in the tank;

[0050] FIG. 15 is a cut plot through a row of bottles showing velocity distribution in the tank;

[0051] FIG. 16 is a graph showing the relationship of specific enthalpy with temperature for whole milk;

[0052] FIG. 17 is a graph showing conductivity of the milk;

[0053] FIG. 18 is a graph showing temperature reduction over time at 11 equally distributed increments for a first case of preserving milk;

[0054] FIG. 19 is a graph showing energy levels over time for a first case of preserving milk;

[0055] FIG. 20 is a graph showing temperature reduction over time at 11 equally distributed increments for a second case of preserving milk;

[0056] FIG. 21 is a graph showing energy levels over time for a second case of preserving milk;

[0057] FIG. 22 is a graph showing temperature reduction over time at 11 equally distributed increments for a third case of preserving milk;

[0058] FIG. 23 is a graph showing energy levels over time for a third case of preserving milk;

[0059] FIG. 24 is a graph showing temperature reduction over time at 11 equally distributed increments for a fourth case of preserving milk;

[0060] FIG. 25 is a graph showing energy levels over time for a fourth case of preserving milk;

[0061] FIG. 26 is a graph showing the relationship of specific enthalpy with temperature for fish;

[0062] FIG. 27 is a graph showing conductivity of the fish product;

[0063] FIG. 28 is a graph showing temperature reduction over time at 11 equally distributed increments for a first case of preserving fish;

[0064] FIG. 29 is a graph showing energy levels over time for a first case of preserving fish;

[0065] FIG. 30 is a graph showing temperature reduction over time at 11 equally distributed increments for a second case of preserving fish;

[0066] FIG. 31 is a graph showing energy levels over time for a second case of preserving fish;

[0067] FIG. 32 is an image showing computational fluid dynamics analysis of temperature in a tank containing meat;

[0068] FIG. 33 is an image showing computational fluid dynamics analysis of velocity in the tank containing meat;

[0069] FIG. 34 is a cut plot through the tray of meat showing temperature distribution in the tank;

[0070] FIG. 35 is a cut plot through the tray of meat showing velocity distribution through the tank;

[0071] FIG. 36 is a graph showing temperature reduction over time at 11 equally distributed increments for a third case of preserving fish;

[0072] FIG. 37 is a graph showing energy levels over time for a third case of preserving fish;

[0073] FIG. 38 is a graph showing temperature reduction over time at 11 equally distributed increments for a fourth case of preserving fish;

[0074] FIG. 39 is a graph showing energy levels over time for a fourth case of preserving fish;

[0075] FIG. 40 is a graph showing the relationship of specific enthalpy with temperature for meat;

[0076] FIG. 41 is a graph showing conductivity of the meat product;

[0077] FIG. 42 is a graph showing temperature reduction over time at 11 equally distributed increments for a first case of preserving meat;

[0078] FIG. 43 is a graph showing energy levels over time for a first case of preserving meat;

[0079] FIG. 44 is a graph showing temperature reduction over time at 11 equally distributed increments for a second case of preserving meat;

[0080] FIG. 45 is a graph showing energy levels over time for a second case of preserving meat;

[0081] FIG. 46 is a graph showing temperature reduction over time at 11 equally distributed increments for a third case of preserving meat;

[0082] FIG. 47 is a graph showing energy levels over time for a third case of preserving meat;

[0083] FIG. 48 is a graph showing temperature reduction over time at 11 equally distributed increments for a fourth case of preserving meat;

[0084] FIG. 49 is a graph showing energy levels over time for a fourth case of preserving meat;

[0085] and

[0086] FIG. 50 is a piping and instrumentation diagram of the refrigeration system.

DETAILED DESCRIPTION

[0087] Immersion Tank

[0088] FIGS. 1 and 2 show an immersion tank 1 for preservation of a consumable product. In this example, the immersion tank 1 is shown containing bottles and has a capacity to hold up to 100 milk bottles. The tank 1 is constructed of steel to conform with ASTM A240. The tank 1 has two heat exchange fluid inlets 2 and two heat exchange fluid outlets 3, the inlets 2 being situated on an inlet wall 4 and the outlets 3 being situated on an outlet wall 5. FIG. 2 shows an inner housing 10 situated internally of the outer walls of the tank 1. Inner housing 10 has an inlet wall 14, outlet wall 15 and base 16, each including apertures 11 to allow inflow and outflow of heat exchange fluid into and out of the inner housing 10. The apertures 11 are provided in four rows of ten on the inlet wall 14 and outlet wall 15, and ten rows of ten on the base 16. The apertures 11 on inlet wall 14 and base 16 are 10 mm in diameter and the apertures 11 on the outlet wall 15 are 20 mm in diameter. The inlet wall 14 is spaced 100 mm away from an inner face 12 of the inlet wall 2, thus providing void 13. A similar void is provided between outlet wall 15 and an inner face of outlet wall 5. A 100 mm void space is further provided in the base. Inner face 12 is defined by steel sheet formwork arranged 50 mm from the inlet wall 2 and secured by brackets, providing a cavity into which polyurethane foam insulation is pumped during manufacture of the tank 1. Insulation is provided in a similar manner in all four walls of the tank from the top of the tank to approximately 595 mm down the walls of the tank. Rows of holes 22 of 30 mm diameter are provided along strips 23 which sit at an angle of approximately 45° between the base 7 and the walls of the tank along the bottom of each wall. The strips 23 are provided to brace the tank structure and can also be used as guides to prevent the trays or basket resting against the walls or base of the inner housing 10. It will be appreciated that other arrangements are possible which also brace the tank structure and perform a guide function. The holes 22 help to reduce stagnation of the heat exchange fluid that may accumulate in these regions of the tank due to the presence of the strips 23.

[0089] A drain 6 is provided from the base 7 of the tank 1 and is shaped as an elbow pipe directed to extend beyond the outlet wall 5 of the tank 1, below the heat exchange fluid outlets 3. The heat exchange fluid inlets 2 and the heat exchange outlets 3 have a diameter of 80 mm.

[0090] The tank 1 further includes a lid formed of steel sheet (not shown). The base 7 of the tank 1 includes four central leg portions 8 supporting the central weight of the tank 1, as well as feet 16 situated at the corners of the tank 1 and formed at the ends of the tank walls. Cut-out portions 9 are provided on the lower ends of the tank walls to provide access for maintenance of the base 7 of the tank. The tank 1 has a height of about 1.105 m and is arranged in a square configuration having side lengths of 1.705 m. In the case of freezing milk, as shown in FIG. 3, a structure comprising a bottle basket 17 is suspended from the lid (not shown) by a bar 18 having hooks at its opposed ends, the bar 18 being attachable to the lid by attachment portion 19.

[0091] As shown in FIG. 4, a structure comprising 40 kg rack trays 21 can be suspended from the lid for preserving meat and fish.

[0092] FIG. 5 shows an upper perspective view of the meat/fish rack trays 21 in the tank 1. FIG. 6 shows the bottle basket 17 in the tank 1. The bottle basket 17 includes receiving holes 20 for holding the sides of the bottles.

[0093] FIG. 7 shows the bottle basket 17 in isolation. The basket 17 has a base 24 and opposed sides 30. The sides 30 are attached to opposed ends of bar 18 which are to be attached to the lid via attachment portion 19. Bottle stabiliser portion 20a includes bottle-receiving portions 20 and is attached to the basket 17 by clipping into slots 31 provided in opposed sides 30. Further holes 20b are provided between the bottle-receiving portions 20, facilitating fluid flow between the bottles when the basket 17 is filled with bottles and immersed in heat exchange fluid.

[0094] FIG. 8 shows the meat/fish rack tray 21 in isolation. The rack tray 21 also includes opposed sides 30 into which individual trays 28 can clip at slots 31. The individual trays 28 are each closed with latches 25. The rack tray 21 is hung from opposed ends of bar 18 and attached to an underside of the lid of the tank by attachment portion 19. A user can grasp handle 29 to slide out individual trays 28 from the rack 21.

[0095] FIG. 9 shows an individual tray 28 in isolation. The tray has shaped portions 26 for enclosing individual pieces of meat or fish. The tray 28 is closed about a hinge joint 27 and secured by latch 25.

[0096] The tank 1 is filled with heat exchange fluid which does not freeze above −70° C. The heat exchange fluid is pumped into the tank 1 via the heat exchange fluid inlets 2 into cavity 13 at a volumetric flow rate of 17 cubic metres per hour. Pressure is built up in the cavity 13 as heat exchange fluid is forced through the restricted areas of the apertures 11, thus reducing the volumetric flow rate but increasing velocity of the fluid entering the inner housing 10. Some fluid will also travel below the inner housing 10 and be forced up to the opposing cavity in the outlet wall 5, with some fluid also travelling up through apertures 11 provided in the base 7 of the inner housing 10. The apertures 11 provide improved distribution of cold fluid to all parts of the tank (see FIG. 10, for example, showing temperature distribution throughout the tank) and minimise the occurrence of hot spots which would otherwise be likely to occur away from the inlet area. As the heat transfer fluid flows continuously through the tank 1, heat is removed from the milk and milk bottles, and the heated heat exchange fluid leaving the tank 1 will then be exchanged with a refrigeration system which continuously cools the heat exchange fluid. The heat exchange fluid itself exchanges heat with refrigerant in the refrigeration system.

[0097] Preferably, a low range heat transfer fluid is used as the immersion fluid for the tank which, advantageously, has a relatively low viscosity even at very low temperatures, thus reducing the pump power requirements for the system. The below table (Table 1) specifies some of the thermal properties of the heat transfer fluid.

TABLE-US-00001 TABLE 1 Thermal properties of the heat transfer fluid Temperature Density Specific Heat Conductivity Viscosity [° C.] [kg/m3] [J/kg-K] [W/m-K] [Pa-s] −73.3 827 1630 0.1591 0.154 −59 820 1717 0.1585 0.038 −45.6 815 1760 0.1574 0.016 −17.8 791 1840 0.1539 0.005 23.9 755 2010 0.1504 0.0015

[0098] At each of the above temperatures, the heat exchange fluid has a density that is very low and lesser than that of water. Advantageously, if any breakage or spillage were to occur during operation of the tank, the broken or spilled matter will tend to sink to a lower portion of the tank, facilitating drainage of that matter without substantial loss of heat exchange fluid. It will be appreciated that any suitable heat exchange fluid can be used, provided that it has a low enough viscosity that it will not require excessive pump power at the required low temperatures for preservation. It is also preferable that the heat exchange fluid be food safe.

TABLE-US-00002 TABLE 2 Flow rate of the heat transfer fluid for various temperature differences @ 20 kW (−50° C.) Temperature Mass Volumetric difference flow rate flow rate [° C.] [kg/s] [m.sup.3/hr] 1 11.5 50.9 2 5.8 25.5 3 3.9 17 4 2.9 12.7 5 2.3 10.1 10 1.15 5.1

[0099] Table 2 above provides the temperature difference between the tank inlet and outlet for various flow rates of heat exchange fluid, assuming 20 kW of heat is extracted from the fluid in the tank. From Table 2, it can be seen that a temperature difference of 3° C. between inlet and outlet can be achieved using a mass flow rate of approximately 4 kg/s. This temperature difference was deemed an acceptable temperature rise in terms of evaporator duty required as well as cooling of the product required. The acceptable temperature rise must be balanced against costs associated with the maximum number of product that can be processed at once to make the system commercially viable. It will be appreciated that a higher flow rate may be desirable in increasing the heat transfer between the heat exchange fluid and the consumable product. However, a higher flow rate will also cause higher flow resistance and thus a higher pumping power would be required.

[0100] Computational Fluid Dynamics Analysis of Milk

[0101] Computational fluid dynamics analysis is performed on the tank in order to visualise how the heat exchange fluid flows in the tank, to estimate the heat transfer coefficient between the heat exchange fluid and the consumable product, and to determine the pressure loss as the heat exchange fluid flows through the tank.

[0102] FIG. 10 is a side view of a simulated flow analysis of heat exchange fluid in the tank for the case where the bottle surface temperature is specified as −30° C. The figure is shade-coded to represent the various temperatures at different points within the tank, with black representing the warmest zones and white representing the coldest zones. The arrows indicate the flow trajectories of the fluid as it moves through the tank and around the bottles. Cold fluid at −50° C. is able to surround the bottles due to the distribution provided by apertures 11. As can be seen from this figure, the heat exchange fluid enters the tank at approximately −50° C. through the inlet wall at the left side of the tank in this view. As the heat exchange fluid absorbs heat from the consumable product (in this case, milk bottles), slightly warmer heat exchange fluid at approximately −47° C. accumulates at the top of the tank before it exits through the outlet wall. A 3° C. increase in temperature of heat exchange fluid was considered acceptable.

[0103] FIG. 11 is a side view of a simulated flow analysis of heat exchange fluid in the tank for the same case, though the figure is shade-coded to represent the various flow velocities. The velocities in the tank around the milk bottles is very low, i.e. less than 0.02 m/s. However, the velocities are high at the inlet wall (on the left of the image) as the fluid passes through the apertures, at various aperture points in the base and somewhat higher at the aperture of the outlet wall (on the right of the image).

[0104] FIG. 12 is a top view of the tank under the same conditions as those for FIGS. 10 and 11, shade-coded for temperature distribution.

[0105] A similar process of analysis is performed at various other bottle surface temperatures and the below table (Table 3) is produced, showing the average heat transfer coefficients determined between the milk (including bottles) and the heat exchange fluid.

TABLE-US-00003 TABLE 3 Results of computational fluid dynamics analysis of milk bottles Bottle Bottles Rate of heat Heat Surface Inlet total transferred Transfer Temperature Temperature surface area from bottles Coefficient [° C.] [° C.] [m.sup.2] [W] [W/m.sup.2K] 0 −50 7.1655 65138 181 −10 −50 7.1655 38964 135 −20 −50 7.1655 29204 135 −30 −50 7.1655 20307 141 −40 −50 7.1655 11892 165

[0106] As the flow velocity of the heat exchange fluid was set very low, the convective heat transfer from the bottles to the heat exchange fluid is largely driven by natural movement of the fluid caused by the local heating of the fluid and the corresponding density differences, i.e. buoyancy-driven flows. Pressure drop is from the inlet to the outlet is found to be approximately 660 Pascal for a flow rate of 4 kg/s.

[0107] Heat Transfer Analysis on Milk

[0108] Further analysis is performed by investigating the influence of varying input parameters of the refrigeration system. This may include the geometry of the product, the temperature of the product, the characteristics of the packaging and the characteristics of the racking systems utilised. The method involves dividing the consumable product into geometrical increments (e.g. cylindrical shells for bottles). For every one of these increments, a conservation of energy equation is solved, i.e. for a given time-step, a certain amount of energy is removed from a shell, resulting in a decrease in temperature of that shell. The amount of energy removed is a function of the temperatures of the adjacent shells, as well as the resistance to heat flow between the shells. This involves taking into account thermal properties of the consumable product as a function of temperature.

[0109] Analysis is performed assuming that the milk can be treated as a solid mass having a starting temperature of 2° C. and the properties given in the ASHRAE Engineering Handbook—Refrigeration, Thermal Properties of Foods (including protein, fat and water contents). FIG. 16 is a plot of specific enthalpy values of whole milk at different temperatures obtained from the Handbook. FIG. 14 is a graph showing the conductivity of whole milk as determined by the Kopelman method.

[0110] A number of cases are investigated to assess different geometries of milk bottle products and different temperatures of the heat exchange fluid. The table below (Table 4) provides a summary of investigated cases.

TABLE-US-00004 TABLE 4 Investigated cases for different milk bottle geometries and inlet temperature of heat exchange fluid Heat Exchange Heat Geometry Fluid Exchange (diameter Outer Temperature Fluid Flow HTC # Product and volume) Layer [° C.] [kg/s] [W/m.sup.2-K] 1 Milk 85 mm 1 mm glass −50 4 130 bottle - 1000 mL 2 Milk 85 mm 1 mm glass −70 4 130 bottle - 1000 mL 3 Milk 69 mm 1.5 mm plastic −50 4 130 bottle - 500 mL 4 Milk 69 mm 1.5 mm plastic −70 4 130 bottle - 500 mL

[0111] In view of the results obtained from the computational fluid dynamics analysis shown in Table 3, a conservative heat transfer coefficient was estimated as 130 W/m.sup.2.Math.K, which was obtained by taking the lowest simulated heat transfer coefficient and rounding it down to the nearest ten. The table below (Table 5) shows the freezing time results for the above investigated cases as predicted by simulation software.

TABLE-US-00005 TABLE 5 Freezing time results for milk Time for Heat product Time for Exchange average product core Fluid temperature temperature Temperature to reach −30° C. to reach −30° C. # Geometry [° C.] [minutes] [minutes] 1 85 mm bottle −50 39.5 43.5 2 85 mm bottle −70 25.5 32 3 69 mm bottle −50 39 41 4 69 mm bottle −70 26 29

[0112] FIG. 18 is a graph showing the temperature-time relationship at the eleven geometrical equally distributed increments from the surface of the bottles to the core for investigated case #1. As can be seen from FIG. 18, a rapid drop in temperature is observed within the first five minutes for the outermost layer of the bottle, whereas a significant drop in temperature at the core is not observed until approximately 41 minutes. Additionally, the core temperature does not reach −30° C. until approximately 43.5 minutes. It is hypothesised that this outcome may minimise ice crystal formation.

[0113] FIG. 19 shows the total enthalpy over time for investigated case #1, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively. These target temperatures of −10° C. and −30° C. were selected as relevant targets based on a number of factors. For example, microbial activity is killed or slowed at different temperatures based on the microbe concerned, as well as at different rates of cooling and duration of time at the target temperature. In particular, storage of foods at −10° C. for a few days is lethal to certain microbes including Toxoplasma gondii, entamoeba and trypanosomes. However, these temperatures must be achieved within a reasonable time frame for commercial viability and, additionally, quick freezing is necessary to minimise ice crystal formation resulting in cell damage. The selected values balance these factors of economic considerations, potential cell damage (resulting in sensory deterioration of food products) and microbe destruction and/or retardation.

[0114] It was decided that the refrigeration system must preserve 100 1 L milk bottles starting at 4° C. to a final temperature of −30° C. within 35 minutes.

[0115] Mass of milk is calculated at 103.8 kg

[0116] Specific enthalpy of milk at 4° C. is 369 kJ/kg

[0117] Specific enthalpy of milk at −30° C. is 18.5 kJ/kg

[0118] The required energy removal from the milk is given by:

m(h.sub.initial−h.sub.final)  (1)

[0119] where

[0120] m=mass of milk

[0121] h.sub.initial=initial specific enthalpy of milk (at 4° C.)

[0122] h.sub.final=final specific enthalpy of milk (at −30° C.)

[0123] Therefore, the heat removed from the milk is:

103.8 (369−18.5)=36 382 (kJ)

[0124] For a freezing period of 35 minutes, the required rate of heat removal is:

36 382/(35 (mins)×60 (secs)=17.3 (kW)

[0125] Assuming a gain of 1 kW heat from atmosphere, the required refrigeration system duty is 18.3 kW. Freezing of milk is seen as the highest load on the refrigeration system thus the system must be specified at least to be able to preserve 100 1 L bottles of milk in 35 minutes, i.e. at least 18.3 kW. Accordingly, an evaporator duty of 20 kW is assumed as an initial design condition for the refrigeration system, including a safety margin of 1.7 kW. The inlet temperature of the tank may be −50° C. or −70° C.

[0126] Returning to case #1, the 85 mm diameter glass bottles of milk are being cooled by a heat exchange fluid having an inlet temperature of −50° C. and a flow rate of 4 kg/s, and a heat transfer coefficient of 130 W/m.sup.2.Math.K. The thickness of the glass outer layer is 1 mm.

[0127] FIG. 20 shows the temperature profiles at equally distributed increments of case #2. In this case, the milk is contained in glass bottles of 85 mm diameter, 1 mm thickness, and a capacity of 1 L. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 130 W/m.sup.2.Math.K.

[0128] FIG. 21 shows the calculated energy over time for investigated case #2, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0129] FIG. 22 shows the temperature profiles at equally distributed increments of case #3. In this case, the milk is contained in plastic bottles of 69 mm diameter, 1.5 mm thickness and 500 mL volume. The heat exchange fluid has an inlet temperature of −50° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 130 W/m.sup.2.Math.K.

[0130] FIG. 23 shows the calculated energy over time for investigated case #3, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0131] FIG. 24 shows the temperature profiles at equally distributed increments of case #4. In this case, the milk is contained in plastic bottles of 69 mm diameter, 1.5 mm plastic and 500 mL volume. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 130 W/m.sup.2.Math.K.

[0132] FIG. 25 shows the calculated energy over time for investigated case #4, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0133] As shown by comparison of FIGS. 18, 20, 22 and 24, the freezing times between the 85 mm diameter bottles and the 69 mm bottles did not significantly vary. The 1.5 mm thickness plastic bottles have higher heat resistance than the 1 mm thickness glass bottles. The surface area of the 500 mL bottles is significantly smaller than that of the 1000 mL bottles.

[0134] Heat Transfer Analysis on Fish

TABLE-US-00006 TABLE 6 Results of computational fluid dynamics on meat and fish pieces Temperature Total Rate of heat Temperature of heat surface transferred Heat of product exchange area of from Transfer surface fluid at inlet products products Coefficient [° C.] [° C.] [m2] [W] [W/m2K] 0 −50 9.98 68839 138 −20 −50 9.98 38215 128 −30 −50 9.98 23673 119 −40 −50 9.98 10147 102

[0135] The starting temperature of fish was also assumed to be 2° C. The type of fish is assumed to be cod, the properties of which were obtained from the ASHRAE Engineering Handbook—Refrigeration, Thermal Properties of Foods. FIG. 26 shows the obtained specific enthalpy values for fish. FIG. 24 shows the conductivity values obtained for fish.

[0136] In view of the results obtained from the computational fluid dynamics analysis shown in Table 6, a conservative heat transfer coefficient was estimated as 100 W/m.sup.2.Math.K, which was obtained by taking the lowest simulated heat transfer coefficient and rounding it down to the nearest ten.

[0137] The cases investigated for fish are given in the table below (Table 7).

TABLE-US-00007 TABLE 7 Investigated cases for different fish geometries and different heat exchange fluid inlet temperatures Heat Exchange Heat Fluid Exchange Temperature Fluid Flow HTC # Product Geometry Outer Layer [° C.] [kg/s] [W/m.sup.2-K] 1 Fish 50mm thick, 0.1 mm plastic −50 4 100 100 × 50 mm 2 Fish 50 mm thick, 0.1 mm plastic −70 4 100 100 × 50 mm 3 Fish 25 mm thick, 0.1 mm plastic −50 4 100 100 × 100 mm 4 Fish 25 mm thick, 0.1 mm plastic −70 4 100 100 × 100 mm

[0138] As shown above, the flow rate of heat exchange fluid and thus the heat transfer coefficient (HTC) remain the same for each of the above cases. The table below (Table 7) shows the freezing time results for the above investigated cases as predicted by simulation software.

TABLE-US-00008 TABLE 8 Freezing time results for fish Time for Heat product Time for Exchange average product core Fluid temperature temperature Temperature to reach −30° C. to reach −30° C. # Geometry [° C.] [minutes] [minutes] 1 50 mm thick, −50 26.5 29.5 100 × 50 mm 2 50 mm thick, −70 17.5 20 100 × 50 mm 3 25 mm thick, −50 23 24 100 × 100 mm 4 25 mm thick, −70 15 16 100 × 100 mm

[0139] The temperature changes over time for eleven equally distributed increments from surface to core of the fish are calculated for each of the above cases. FIG. 28 shows the temperature profiles for case #1. In this case, the fish is 50 mm thick and has a length and width of 100 mm and 50 mm respectively. The fish is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −50° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0140] FIG. 29 shows the calculated energy over time for investigated case #1, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0141] FIG. 30 shows the temperature profiles for case #2. In this case, the fish is 50 mm thick and has a length and width of 100 mm and 50 mm respectively. The fish is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0142] FIG. 31 shows the calculated energy over time for investigated case #2, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0143] FIG. 36 shows the temperature profiles for case #3. In this case, the fish is 25 mm thick and has a length and width of 100 mm and 100 mm respectively. The fish is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −50° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0144] FIG. 37 shows the calculated energy over time for investigated case #3, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0145] FIG. 38 shows the temperature profiles for case #4. In this case, the fish is 25 mm thick and has a length and width of 100 mm and 100 mm respectively. The fish is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0146] FIG. 39 shows the calculated energy over time for investigated case #4, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0147] As can be seen from the above results, the thinner pieces of fish can be frozen more quickly than the thicker pieces of fish.

[0148] Heat Transfer Analysis on Meat

[0149] The starting temperature of meat was also assumed to be 2° C. The type of meat is assumed to be lean sirloin, the properties of which were obtained from the ASHRAE Engineering Handbook—Refrigeration, Thermal Properties of Foods. FIG. 40 shows the obtained specific enthalpy values for fish. FIG. 41 shows the conductivity values obtained for fish.

[0150] The cases investigated for meat are given in the table below (Table 9).

TABLE-US-00009 TABLE 9 Investigated cases for different meat geometries and different heat exchange fluid inlet temperatures Heat Heat Exchange Exchange Fluid Fluid Temperature Flow HTC # Product Geometry Outer Layer [° C.] [kg/s] [W/m.sup.2-K] 1 Meat 50 mm thick, 0.1 mm plastic −50 4 100 100 × 50 mm 2 Meat 50 mm thick, 0.1 mm plastic −70 4 100 100 × 50 mm 3 Meat 25 mm thick, 0.1 mm plastic −50 4 100 100 × 100 mm 4 Meat 25 mm thick, 0.1 mm plastic −70 4 100 100 × 100 mm

[0151] As shown above, the flow rate of heat exchange fluid and thus the heat transfer coefficient (HTC) remain the same for each of the above cases. The table below (Table 10) shows the freezing time results for the above investigated cases as predicted by simulation software.

TABLE-US-00010 TABLE 10 Freezing time results for meat Time for Heat product Time for Exchange average product core Fluid temperature temperature Temperature to reach −30° C. to reach −30° C. # Geometry [° C.] [minutes] [minutes] 1 50 mm thick, −50 25.5 28 100 × 50 mm 2 50 mm thick, −70 16.5 19.5 100 × 50 mm 3 25 mm thick, −50 22 22.5 100 × 100 mm 4 25 mm thick, −70 14.5 15 100 × 100 mm

[0152] The temperature changes over time for eleven equally distributed increments from surface to core of the meat are calculated for each of the above cases. FIG. 42 shows the temperature profiles for case #1. In this case, the meat is 50 mm thick and has a length and width of 100 mm and 50 mm respectively. The meat is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −50° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0153] FIG. 43 shows the calculated energy over time for investigated case #1, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0154] FIG. 44 shows the temperature profiles for case #2. In this case, the meat is 50 mm thick and has a length and width of 100 mm and 50 mm respectively. The meat is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0155] FIG. 45 shows the calculated energy over time for investigated case #2, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0156] FIG. 46 shows the temperature profiles for case #3. In this case, the meat is 25 mm thick and has a length and width of 100 mm and 100 mm respectively. The meat is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −50° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0157] FIG. 47 shows the calculated energy over time for investigated case #3, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0158] FIG. 48 shows the temperature profiles for case #4. In this case, the meat is 25 mm thick and has a length and width of 100 mm and 100 mm respectively. The meat is contained within 0.1 mm plastic. The heat exchange fluid has an inlet temperature of −70° C. and a mass flow rate of 4 kg/s, and the heat transfer coefficient is 100 W/m2.Math.K.

[0159] FIG. 49 shows the calculated energy over time for investigated case #4, as compared against the energy level at which the average temperature of the product is at −10° C. and −30° C. respectively.

[0160] As can be seen from the above results, the thinner pieces of meat can be frozen more quickly than the thicker pieces of meat.

[0161] In the preferred embodiments of the present invention, the consumable products are contained in packaging or within the basket or tray structures used for immersion in the tanks. Advantageously, this prevents sublimation of the product during the preserving process, which is a problem with existing immersion preservation methods that do not use packaging due to issues with cracking.

[0162] Advantageously, using the above analyses, the system can be used to reduce the temperature of the food product to a target temperature (e.g., −50° C.) within a desired amount of time (e.g. 30 minutes) at a desired flow rate (e.g. 4 kg/s). Input parameters of the apparatus are set based on the findings of the analyses. This allows a user to simply select the product (e.g., “fish fillet”) when using the system. Preferably, thin packaging is used with minimal air spaces, such as vacuum sealing or heat shrink wrapping.

[0163] While a theoretical freezing rate of 10° C. per minute or 100° C. per minute may result in even greater elimination results for the microbe levels in the consumable product, as well as less cellular damage, such rates are not viable for a commercial freezing system. Moreover, surprisingly it has been found that significant levels of microbe reduction can be achieved at a slower cooling rate, such as a rate of about 1° C. per minute, while retaining a high degree of sensory qualities and preventing damage to packaging.

[0164] Although heat transfer analysis in the abovementioned examples was performed on 11 volume increments of each consumable food product, it will be appreciated that this number can be varied.

[0165] The table below provides results for milk preserved over a cooling period of 45 minutes to a temperature of −30° C. immersed in heat exchange fluid at −50° C. The results have been described as similar to those of the pasteurisation process.

TABLE-US-00011 TABLE 11 Analysis Results of preserved raw milk Sample Analysis Results Description: Test Result Unit Raw Milk Total Viable Count 380 cfu/g E. coli <10 cfu/g Coliform <10 cfu/g Listeria spp. Not Detected in 25 g Salmonella spp. Not Detected in 25 g Coagulase + ve <100 cfu/g Staphylococcus Yeast <10 cfu/g Mould <10 cfu/g

[0166] Refrigeration System

[0167] FIG. 50 is a piping and instrumentation diagram of the refrigeration system that continuously cools the heat exchange fluid.

[0168] The refrigeration system includes a heat exchanger for exchanging heat between the heat transfer fluid and the refrigerant, which can be R404A.