# Analysis of Heat Transfer Mechanism for Shelf Vacuum Freeze-Drying Equipment

Authors: Hong-Ping Cheng; Shian-Min Tsai; Chin-Chi Cheng
Journal: Advances in Materials Science and Engineering
 (2014)

---

## Abstract

Vacuum freeze-drying technology is applicable to the process of high heat-sensitive products. Due to the long drying period and extremely low processing temperature and pressure, the uniform and efficiency of heat transfer fluid temperature in shelf are critical for product quality. Therefore, in this study, the commercial computer fluid dynamics (CFD) software, FLUENT, was utilized for three-dimension numerical simulation of the shelf vacuum freeze-drying process. The influences of different inlet and outlet positions for shelves on the uniformity of the flow rate and temperature were discussed. Moreover, it explored the impacts on the temperature gradient of shelves after heat exchange of different flow rates and low temperature materials. In order to reduce the developing time and optimize the design, the various secondary refrigerants in different plies of shelves were investigated. According to the effect of heat exchange between different flow rates and low temperature layer material shelves on the temperature gradient of shelves surface, the minimum temperature gradient was 20 L/min, and the maximum was 2.5 L/min.


---

## Body

## 1. Introduction

Vacuum freeze-drying technology is widely used for manufacturing pharmaceuticals. The requirements of this technology are process optimization and long-term storage stability. There are three stages in the complete vacuum freeze-drying process, such as freezing, primary drying, and secondary drying. Most of the water is frozen to the ice during the freezing stage, and then the ice is removed by sublimation during the primary drying stage. During the secondary drying stage, the temperature is usually elevated to remove the residual water of the products [1, 2]. Due to the extremely rigor conditions of vacuum freeze-drying process, such as low temperature, high vacuum, and long processing period, the process stability and efficiency would be crucial for product quality.In the vacuum freeze-drying process, before all the ice is removed from the product, increasing the shelf temperature into the secondary drying would likely cause collapse or eutectic melt of the products. Therefore, the controllability of the temperature of the heat-transfer fluid in the shelf is very important for product quality during the freeze-drying process. Many studies would like to find out how the temperature of the heat-transfer fluid in the shelf affected the freeze-drying process. Liapis and Bruttini [3] constructed the operational control policies by the presented exergy expressions to minimize the irreversibilities occurring in a freeze-drying system and enhance the efficiency of energy utilization in lyophilization. Patel et al. [4] measured the temperature of 5% mannitol and 5% sucrose by thermocouple to determine the setting shelf temperature and understood the time of primary drying end by product temperature. Gieseler et al. [5] investigated the better drying parameters by measuring the total amount of removed water from the product of pure water and a 5% (w/w) mannitol at different pressure and shelf temperature settings, performed in both laboratory and pilot scale freeze dryers. For the pharmaceutical and protein formulations, the collapse temperature and the thermal transition were utilized to find out the optimal heat-transfer fluid temperature of shelf and chamber pressure settings required to achieve the target product temperature [6, 7]. For food products, Sablani and Rahman [8] indicated that the heat-transfer temperature of shelf maintained in −45~15°C during primary drying process would result in the best quality of drying abalone, potato, and brown dates.From the above literatures, one can see that the heat-transfer fluid temperature of shelf is one of the most important factors in drying process. Therefore, some technologies, such as manometric temperature measurement (MTM), had been used to determine the optimal shelf temperature, and the expert system was utilized to estimate the required shelf temperature for maintaining product temperature [9, 10]. In addition, some studies employed the finite element analysis for 2D simulation of sublimation and heat transfer [11, 12], discussing the contact heat transfer between shelf and glass bottle [13, 14] and investigating the temperature distribution between glass bottle and product during the drying cycle [15, 16]. However, the rare literature discussed the flow pattern of fluid inside shelves, design of shelves, and effect of fluid flow and temperature.The objective of this paper is to investigate the heat transfer rate of the fluid in shelf by using CFD method. Several factors, such as the location of inlet and outlet of the fluid in shelf and mass flow rate, are also discussed in order to improve the efficiency of shelf for the laboratory scale vacuum freeze-drying system.


## 2. Operating Conditions

In this paper, a commercial CFD software, FLUENT, was adopted to simulate all heat transferring processes in the heat exchanger, in order to optimize the design. Shelf is one of the most important devices in the vacuum freeze-drying system, because of its capability to heat and cool the products. Therefore, several designs and their operating procedures will be investigated in the study. The efficient shelf is able to reduce the required energy of vacuum freeze-drying equipment during drying process and uplifting the quality of dried product. The 50% ethylene glycol solution flowing into the shelf served as the secondary refrigerant. Figure1 shows the schematic drawing of the freeze-drying system, including (a) vacuum chamber, (b) shelf, (c) cold trap, (d) vacuum pump, (e) heat exchanger, (f) condensers, (g) expansion valve, and (h) compressors.Figure 1
Schematic drawing of the freeze-drying system.This study would focus on numerical simulation of the flow field and heat-transfer fluid properties conducted for internal flow channel of shelves, including different inlet positions and inlet velocities. The simulation included five cases, as shown in Figure2 for the detail schematic drawings. Case 1 is the geometric shape of the shelves of vacuum freeze-drying machine. The inlet position is at the top, and the outlet position is at the bottom, as shown in Figure 2(a). For Case 2, the inlet position is at the bottom, and the outlet position is at the top, as shown in Figure 2(b). For Case 3, both the inlet and outlet are at the bottom, as shown in Figure 2(c). For Case 4, both the inlet and outlet are at the top, as shown in Figure 2(d). For Case 5, the increased part is at the bottom, and inlet and outlet are both at the bottom, as shown in Figure 2(e).Geometric shapes of (a) Case 1, (b) Case 2, (c) Case 3, (d) Case 4, and (e) Case 5.
(a)
(b)
(c)
(d)
(e)


### 2.1. Boundary Conditions

A CFD software, FLUENT, was adopted to simulate the heat exchanging process. Several assumptions during the simulation were made as follows:(1) steady-state flow field; (2) the viscous model uses Laminar; (3) the pressure-velocity coupling uses Simple; (4) the convection-diffusion term uses First Order Upwind; (5) the material is maintained at a constant temperature of 233 K in the computational process.The numbers of grids in Case 1 to Case 4 were 1304471 and that in Case 5 were 1308602. The computing range under hypothetical simulation conditions covered one inlet, one outlet, and liquid entering into shelves, while others were solid boundaries. The pipes through the shelves were set as square pipes after hydraulic diameter calculation, so that grids in GAMBIT were used as a hex map. The grids were 1 mm per unit inX-direction after grid testing, 2 mm per unit in Y-direction, and 2 mm per unit in Z-direction, as shown in Figure 3. The three-dimensional graphics of shelves and shelf flow channel as shown in Figure 3(a) and the grids exterior were shown in Figure 3(b) to Figure 3(d).Three-dimensional graphics of shelves and shelf flow channel: (a) overall dimensions of shelves and flow channel structure; (b) local grid structure distribution; (c) side view of overall grid structure distribution; (d) top view of overall grid structure distribution.
(a)
(b)
(c)
(d)The boundary conditions of Cases 1 to 5 are as follows.(1) Inlet and outlet boundaries: the boundary condition of inlet is set as velocity inlet; the inlet velocities are 2.5 L/min, 5 L/min, 10 L/min, and 20 L/min, respectively; the inlet temperature was fixed at 253 K; the outlet boundary was set as outflow; (2) heat transfer fluid entering shelves: the fluid boundary condition settings adopt a 50% ethylene glycol solution; the fluid properties had density of 1085 kg/m3, specific heat of 3490 J/kg-K, heat transfer coefficient of 0.42 w/m-k, and viscosity coefficient of 0.029 kg/m-s; (3) solid boundary: this solid boundary was set as a wall; the material was steel, and the wall surface was set as heat flux; it was assumed that the top contacts low temperature material was at a preset temperature of 233 K.


### 2.2. Method of Analysis

For simulating the freeze-drying process of heat transfer fluid in shelf, determination of flow patterns is the first step. Reynolds number is usually adopted to predict the flow patterns in different flow situations. It can be calculated as(1)Re=ρVdiμ,
where ρ is the density (kg/m3),V is the average velocity (m/s), di is the hydraulic diameter (m), and μ is the viscosity (kg/m-s).In this study, the heat transfer fluid in the shelf utilized ethylene glycol solution, with a concentration ratio of 50%. Its thermodynamic properties areρ=1089 kg/m3, di=7.92 mm, and μ=0.0265 kg/m-s. These thermodynamic properties varied with environmental temperature, according to the source of the thermodynamic properties of NIST (National Institute of Standards and Technology). The Reynolds number of the heat transfer fluid in the shelf can be calculated according to (1). The calculated Reynolds number was 12, which is less than 2300. According to the calculated Reynolds number, the inlet flow is decided as the laminar. Thus, the CFD simulation for all the cases was according to the laminar model.


### 2.3. Heat Transfer Rate

This study would focus on the heat-transfer fluid in shelf and heat exchange with freeze products. Therefore, the heat transfer rate was calculated as(2)Qc=hc¯As(Ts-TA),
where Qc is the heat transfer rate, hc¯ is convective heat transfer coefficient, As is the body surface area, Ts is the body wall surface temperature, and TA is the ambient temperature.For calculatingQc, the Nusselt number was calculated first. It can be calculated as
(3)Pr=Cpμk,
where Pr is the Prandtl number, Cp is the specific heat (j/kg-k), μ is the viscosity (kg/m-s), and k is the thermal conductivity (w/m-k). The Prandtl number was located in the range of 0.48 < Pr < 16700.The laminar Nusselt number equation is used for calculatinghc¯ as follows:
(4)Nu=hc¯Dk=1.86(VDρμ)1/3(Cpμk)1/3(DL)1/3(μμw),
where Nu is the Nusselt number, hc¯ is the convective heat transfer coefficient (W/m2-K), D is the hydraulic diameter (m), k is the thermal conductivity (w/m-k), V is the velocity (m/s), μ is the viscosity (kg/m-s), μw is the wall viscosity (kg/m-s), ρ is the density (kg/m3), Cp is the specific heat (j/kg-k), and L is the characteristic length (m).


## 2.1. Boundary Conditions

A CFD software, FLUENT, was adopted to simulate the heat exchanging process. Several assumptions during the simulation were made as follows:(1) steady-state flow field; (2) the viscous model uses Laminar; (3) the pressure-velocity coupling uses Simple; (4) the convection-diffusion term uses First Order Upwind; (5) the material is maintained at a constant temperature of 233 K in the computational process.The numbers of grids in Case 1 to Case 4 were 1304471 and that in Case 5 were 1308602. The computing range under hypothetical simulation conditions covered one inlet, one outlet, and liquid entering into shelves, while others were solid boundaries. The pipes through the shelves were set as square pipes after hydraulic diameter calculation, so that grids in GAMBIT were used as a hex map. The grids were 1 mm per unit inX-direction after grid testing, 2 mm per unit in Y-direction, and 2 mm per unit in Z-direction, as shown in Figure 3. The three-dimensional graphics of shelves and shelf flow channel as shown in Figure 3(a) and the grids exterior were shown in Figure 3(b) to Figure 3(d).Three-dimensional graphics of shelves and shelf flow channel: (a) overall dimensions of shelves and flow channel structure; (b) local grid structure distribution; (c) side view of overall grid structure distribution; (d) top view of overall grid structure distribution.
(a)
(b)
(c)
(d)The boundary conditions of Cases 1 to 5 are as follows.(1) Inlet and outlet boundaries: the boundary condition of inlet is set as velocity inlet; the inlet velocities are 2.5 L/min, 5 L/min, 10 L/min, and 20 L/min, respectively; the inlet temperature was fixed at 253 K; the outlet boundary was set as outflow; (2) heat transfer fluid entering shelves: the fluid boundary condition settings adopt a 50% ethylene glycol solution; the fluid properties had density of 1085 kg/m3, specific heat of 3490 J/kg-K, heat transfer coefficient of 0.42 w/m-k, and viscosity coefficient of 0.029 kg/m-s; (3) solid boundary: this solid boundary was set as a wall; the material was steel, and the wall surface was set as heat flux; it was assumed that the top contacts low temperature material was at a preset temperature of 233 K.


## 2.2. Method of Analysis

For simulating the freeze-drying process of heat transfer fluid in shelf, determination of flow patterns is the first step. Reynolds number is usually adopted to predict the flow patterns in different flow situations. It can be calculated as(1)Re=ρVdiμ,
where ρ is the density (kg/m3),V is the average velocity (m/s), di is the hydraulic diameter (m), and μ is the viscosity (kg/m-s).In this study, the heat transfer fluid in the shelf utilized ethylene glycol solution, with a concentration ratio of 50%. Its thermodynamic properties areρ=1089 kg/m3, di=7.92 mm, and μ=0.0265 kg/m-s. These thermodynamic properties varied with environmental temperature, according to the source of the thermodynamic properties of NIST (National Institute of Standards and Technology). The Reynolds number of the heat transfer fluid in the shelf can be calculated according to (1). The calculated Reynolds number was 12, which is less than 2300. According to the calculated Reynolds number, the inlet flow is decided as the laminar. Thus, the CFD simulation for all the cases was according to the laminar model.


## 2.3. Heat Transfer Rate

This study would focus on the heat-transfer fluid in shelf and heat exchange with freeze products. Therefore, the heat transfer rate was calculated as(2)Qc=hc¯As(Ts-TA),
where Qc is the heat transfer rate, hc¯ is convective heat transfer coefficient, As is the body surface area, Ts is the body wall surface temperature, and TA is the ambient temperature.For calculatingQc, the Nusselt number was calculated first. It can be calculated as
(3)Pr=Cpμk,
where Pr is the Prandtl number, Cp is the specific heat (j/kg-k), μ is the viscosity (kg/m-s), and k is the thermal conductivity (w/m-k). The Prandtl number was located in the range of 0.48 < Pr < 16700.The laminar Nusselt number equation is used for calculatinghc¯ as follows:
(4)Nu=hc¯Dk=1.86(VDρμ)1/3(Cpμk)1/3(DL)1/3(μμw),
where Nu is the Nusselt number, hc¯ is the convective heat transfer coefficient (W/m2-K), D is the hydraulic diameter (m), k is the thermal conductivity (w/m-k), V is the velocity (m/s), μ is the viscosity (kg/m-s), μw is the wall viscosity (kg/m-s), ρ is the density (kg/m3), Cp is the specific heat (j/kg-k), and L is the characteristic length (m).


## 3. Results and Discussion

For the laboratory scale vacuum freeze-drying system, most of the process is carried out on the shelves. How to provide heat equally and efficiently on the shelves during the freeze-drying process would be crucial for product quality. Therefore, in five different types design of shelves, two problems are discussed:(1) the effect of different inlet and outlet positions on the same flow rate and temperature when the fluid flows into different plies of shelves; (2) the effects of different flow rates on the same temperature gradients of shelves after heat exchange between fluid in shelves and low temperature material.In the application of a laboratory scale vacuum freeze-drying machine, only the first and second plies of shelves are used, the third ply is mostly idle. Therefore, only the temperature, velocity, and heat transfer rate distribution of the first and second plies are considered in this study. The results of the third ply are used as a reference.The simulated results were shown in Figure4. In Figure 4, the flow rate of the shelves was 2.5 L/min, In Case 1, the difference of inlet velocity between the first ply and the third ply of the shelf is small. The performances of Qc for the first and second plies were 386.689 W and 388.015 W, respectively. The difference between the first and second plies was 1.326 W. At the flow rate of 2.5 L/min in Case 4, the first and the second plies have the same inlet velocity of 1.154 m/s, and the inlet velocity of third ply was 1.144 m/s. The distribution of heat transfer rate Qc for the first, second, and third plies was 390.728 W, 390.035 W, and 389.098 W, respectively. The difference of heat transfer rate Qc between the first and second plies was 0.693 W, and the difference of the third ply is 1.63 W. The heat transfer rate of the three shelves in Case 4 was more uniform than other cases. Therefore, when the flow rate was 2.5 L/min, it was suggested to use the geometrical model of Case 4 as the reference design.Figure 4
Variation of heat transfer rate with respect to inlet velocity at 2.5 L/min in various cases.In Figure5, the inlet flow rate was set as 5 L/min. The inlet velocities of the first, second, and third plies of shelf in Case 1 were 1.909 m/s, 1.855 m/s, and 1.929 m/s, respectively. In terms of performance in heat transfer rate Qc, the heat transfer rate of the first, second, and third plies was 494.631 W, 491.430 W, and 496.265 W, respectively. The difference of heat transfer rates between the first and second plies was 4.551 W. In Case 3, the heat transfer rates Qc of the first, second, and third plies were 496.524 W, 493.756 W, and 500.159 W, respectively. The difference between the first and second plies was 2.768 W, which was the smallest difference. In addition, the difference of heat transfer rate Qc between the first and second plies of shelves in Case 5 was 2.839 W. Therefore, when the inlet flow rate was 5 L/min, Case 1 could be replaced by the geometric shape of Case 3.Figure 5
Variation of heat transfer rate with respect to inlet velocity at 5 L/min in various cases.In Figure6, the inlet flow rate was set as 10 L/min. The simulating results of Case 1 show that the inlet velocity of the first, second, and third plies was 3.859 m/s, 3.531 m/s, and 3.816 m/s, respectively. The inlet velocity of the second ply was lower than that of the other two plies. In terms of distribution of heat transfer rate Qc, the first, second, and third plies were 634.513 W, 619.403, and 633.628 W, respectively. The difference of Qc between the first ply and the second ply was 14.228 W. In the other cases, the differences of Qc between the first and second plies in Cases 2, 3, and 4 were 11.392 W, 15.535 W, and 13.258 W, respectively. The differences in these cases were above 10 W. However in Case 5, the difference of Qc between the first and second plies was only 5.553 W, which was the smallest among all the cases. Therefore, when the flow rate was 10 L/min, the geometric shape of Case 5 could be used for the design.Figure 6
Variation of heat transfer rate with respect to inlet velocity at 10 L/min in various cases.In Figure7, the inlet flow rate was set as 20 L/min. The shelf inlet velocities of the first and second plies were similar in Case 5. The inlet velocities of the first and second plies were 7.767 m/s and 7.678 m/s, respectively. Their difference was only 0.089 m/s, and the difference in other cases was above 0.357 m/s. For heat transfer rate Qc, the difference was only 4.631 W in Case 5, and the difference in the other cases was 16.641 W~13.435 W. Therefore, when the flow rate was 20 L/min, the geometric shape of Case 5 was suggested.Figure 7
Variation of heat transfer rate with respect to inlet velocity at 20 L/min in various cases.The simulating results have a closer relationship between temperature gradient and heat flow rate in all the cases. Average results of all the cases were shown in Figure8. In the flow rate range from 2.5 L/min to 20 L/min, the temperature gradient decreased from 0.062 to 0.018 K/cm monotonically, and the heat transfer rate increased from 388.7 to 775.5 W monotonically. The flow rate had an inverse linear relationship with the temperature gradient, however, a positive linear relationship with the heat transfer rate.Figure 8
Temperature gradient with respect to heat transfer rate in various inlet velocities.


## 4. Conclusions

Vacuum freeze-drying is a long-term process and has a close relation to temperature. This study discussed the effect of different inlet and outlet positions under the same flow rate and temperature of different plies of shelves and investigated the influences of different flow rates under the same temperature gradient of shelves after heat exchange between fluid and low temperature material in the shelves. The conclusions are as follows.(1)
When the flow rate was 2.5 L/min, the inlet velocities of different plies of shelves are closest in Case 4. The difference of heat transfer rateQc between the first and second plies was only 0.693 W. This was the preferred design for the cases of lower flow rates, and it may result in similar quality of materials even with different plies.(2)
Although Case 1 was an actual case, when the flow rate was 5 L/min, the results show that Case 3 was more appropriate, and the difference between the first and second plies was 2.768 W.(3)
Case 5 has a better heat transfer rate when the flow rates were 10 L/min and 20 L/min. The users can adopt the geometric shape of Case 5 when requiring the higher flow rates. In addition, the heat transfer rate in Case 5 was similar to that of Case 3, when the flow rate was 5 L/min. Therefore, the geometric shape of Case 5 allows a wide range of flow rates.(4)
The temperature gradient and heat transfer rate had a close relationship with the flow rate. The flow rate had an inverse linear relationship with the temperature gradient, however, a positive linear relationship with the heat transfer rate.The above results can serve as references for future studies on the geometric production of shelves and setting inlet flow rates for different flow rates.



---

*Source: 515180-2014-05-14.xml*