Next Article in Journal
The Role of Cold-Water Thermal Refuges for Stream Salmonids in a Changing Climate—Experiences from Atlantic Canada
Previous Article in Journal
Isotopic Turnover and Fractionation of δ15N and δ13C in Captive Pseudopleuronectes americanus (Walbaum)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation Analysis and Experimental Verification of Freezing Time of Tuna under Freezing Conditions

1
College of Food Science and Technology, Shanghai Ocean University, Shanghai 201306, China
2
Shanghai Professional Technology Service Platform on Cold Chain Equipment Performance and Energy Saving Evaluation, Shanghai 201306, China
3
National Experimental Teaching Demonstration Center for Food Science and Engineering, Shanghai Ocean University, Shanghai 201306, China
4
Quality Supervision, Inspection and Testing Center for Cold Storage and Refrigeration Equipment, Ministry of Agriculture, Shanghai 201306, China
*
Authors to whom correspondence should be addressed.
Fishes 2023, 8(9), 470; https://doi.org/10.3390/fishes8090470
Submission received: 15 August 2023 / Revised: 13 September 2023 / Accepted: 20 September 2023 / Published: 21 September 2023

Abstract

:
In order to predict the regular temperature change in tuna during the freezing process for cold chain transportation, improve the quality of frozen tuna, and reduce the energy consumption of freezing equipment, a three-dimensional numerical model for freezing tuna of different sizes was established. An unsteady numerical simulation of the air velocity and flow field was combined with an analysis of the freezing process of tuna. This paper also studied the effect of air velocity, temperature, and tuna size on the freezing process. The numerical results show that there was a positive correlation between the cold source environment and the tuna-freezing process. Lower temperatures and higher air increased the velocity at which the tuna moved through the maximum ice crystal formation zone, maintaining a better aquatic product quality. In some cases, however, the smaller tuna models achieved a longer freezing time. Due to the difficulty of obtaining the whole tuna sample, the temperature curve and freezing rate over time obtained during the freezing process were tested using a tuna block of a specific size. The maximum error did not exceed 6.67%, verifying the authenticity and feasibility of the simulation.
Key Contribution: A system has been built to simulate the freezing process of tuna by combining the flow field and the temperature field, and the accuracy has been verified via experiments.

Graphical Abstract

1. Introduction

The quality of frozen food, especially aquatic products, is affected by the whole process of food cold chain transportation, which includes the pretreatment, freezing, storage, and transportation of frozen food [1]. Freezing can effectively inhibit the activities of microorganisms in aquatic products, reduce enzyme activity, and effectively prolong the shelf life of frozen foods [2,3,4]. Taking the freezing process of tuna as an example, due to its perishable characteristics and in order to prolong its shelf life, its minimum temperature needs to be lowered to below −18 °C as soon as possible after fishing, and the fish then needs to be placed into a freezer to further reduce its temperature to below −60 °C [5].
During the freezing process, along with the decrease in the temperature of tuna, the water will change from the liquid phase to the solid phase, forming ice crystals [6]. Ice crystals cause damage to microbial cells and cause the death of harmful microorganisms in tuna, and at the same time, excessive ice crystals also damage the organizational structure of aquatic products [7]. Therefore, it is worth considering how to increase the freezing rate of tuna to quickly pass through the maximum ice crystal formation zone and reduce the number of large ice crystals formed, which would lead to reduced negative impact of ice crystals on aquatic products [8].
The freezing methods of aquatic products can be mainly divided into three categories according to the freezing speed, namely, fast freezing (5–20 cm/h), medium-speed freezing (1–5 cm/h), and slow freezing (0.1–1 cm/h) [9]. There are differences in the freezing air velocity; the temperature; and the size, shape, and distribution of the ice crystals. As one of the important freezing technologies, rapid freezing can quickly reduce the temperature of the food center to below −18 °C to achieve the purpose of quickly freezing food and reducing cell and muscle damage [10,11,12].
In order to explore the effect of freezing temperature on the quality of big-eye tuna, Jiang et al. [13] tested factors such as color difference and oxidation stability of tuna frozen at different temperatures and found that the high freezing rate brought about by −40 °C and lower temperatures could effectively inhibit the discoloration and protein changes in aquatic products, and the factors of energy saving were taken into account. However, the study did not explain the difference in freezing rate and quality caused by the freezing of different sizes of tuna. Zhao et al. [14] tested samples of frozen golden pomfret at different temperatures and freezing rates and confirmed that a faster freezing rate was an effective method to reduce the effect on aquatic product quality, but the experiment failed to provide a numerical method to predict the freezing rate and process of fish.
To predict the freezing process of bayberry and the effect of freezing on the quality of it, Zhao et al. [15] conducted a numerical simulation of the freezing process of bayberry and tested the freezing speed and quality of bayberry at different temperatures, but did not consider the cold air velocity during freezing. Wan et al. [16] conducted a numerical simulation of frozen shrimp in refrigerators and predicted the temperature distribution and freezing time, but their study lacked experimental verification and analysis of the effects of freezing conditions on shrimp quality.
Therefore, it is significant to study the freezing rate of tuna in different types of freezing equipment to explore the effects of freezing temperature and air velocity on the quality of aquatic products and to find more energy-efficient freezing methods [17,18]. In this paper, tuna is used as the research object, and different freezing methods are tested. With reference to the equipment temperature and air velocity of the refrigerated equipment, the freezing rate of tuna is analyzed through numerical simulation, and the rationality of the model is verified via experiments. The purpose of this paper is to provide a theoretical reference for studying temperature distribution during the freezing process and for predicting freezing time, to explore reasonable freezing conditions that enable tuna to pass the maximum ice crystal formation zone more quickly, and as well as to explain how to increase the freezing rate to ensure quality while ensuring low energy consumption.
  • ω—tuna water content
  • λ—interface thermal conductivity
  • λl—interface liquid-phase thermal conductivity
  • λs—interface solid-phase thermal conductivity
  • ρ—tuna density
  • CFD—computational fluid dynamics
  • Cp—tuna specific heat capacity
  • Cpl—specific heat capacity before freezing
  • Cps—specific heat capacity after freezing
  • h—specific enthalpy
  • hs—heat transfer coefficient of the tuna surface
  • k—thermal conductivity of tuna
  • n—outer normal line direction
  • T—temperature
  • T0—initial temperature
  • Td—bottom surface temperature
  • Text—air temperature
  • Tl—phase interface liquid-phase temperature
  • Tp—phase interface phase transition temperature
  • Ts—phase interface solid-phase temperature
  • TSo—solidus temperature
  • Tw—surface temperature
  • Tfp—phase transition temperature
  • t—freezing time
  • Q—phase transition enthalpy

2. Materials and Methods

2.1. Sample Preparation

Yellowfin tuna from the Pacific Ocean was slaughtered and cut into pieces (with each piece weighing about 500 g) after being caught. Samples were stored after vacuum-freezing at −55 °C. Tuna samples were purchased from Zhengyangfan (Xiamen, China). After being thawed to normal temperature (5 °C), the tuna samples were cut into pieces that were 50 × 30 × 10 mm in size and divided into 5 groups.

2.2. Water Content

Before the experiment, the tuna samples were dried and chopped to make surimi, which was then weighed in a weighing bottle. The water content of the initial sample was measured using a direct drying method to calculate the thermal property parameters [19]. The temperature of the drying oven was adjusted to 100 °C. The weighing bottle was placed in the drying oven with the bottle cap tilted so that it was to the side of the bottle opening. The sample was heated for 1.5 h and left for half an hour, and then the above operation was continued until the quality difference after two drying sessions was less than 2 mg. The water loss after drying was compared with the total mass to determine the water content of the tuna sample. After the above operation was repeated three times, the average value was taken, and the water content obtained during the simulation calculation was determined to be 73%.

2.3. Temperature Curve

A 2640A network multipoint temperature collector (Fluke, Shanghai, China) was connected to a computer before the freezing experiment, and the TT-T-36 calibrated thermocouple temperature probe (Omega Engineering. Inc., Connecticut, America) was fixed at the central point of the tuna block, with the collection interval set at 10 s [20,21]. After the temperature of the refrigerator was adjusted to −18 °C and kept empty for one day before the experiment, the temperature of each point in the freezer was measured with a thermocouple thermometer to ensure that the temperature remained constant at −18 °C.
Based on the above steps, the tuna blocks with fixed thermocouples were taken out of the constant temperature and humidity device and placed in the middle of the refrigerator freezer, making sure that the refrigerator door was rapidly closed, so as to obtain the freezing temperature curves of the center points of tuna pieces under the freezing environment of a temperature of −18 °C and an air velocity of 0 m/s. Because the size of the tuna block is small, the experimental measuring point is not easily fixed; thus, the above steps needed to be repeated to reduce the error. The results show that the temperature curve of the tuna piece at this point is basically the same under the same conditions.

3. Models and Assumptions

3.1. Physical Models

FLUENT was used to simulate the freezing process of tuna, and the changes in the temperature field and basic parameters of tuna under different air velocities, freezing temperatures, and fish meat sizes were calculated. In order to obtain the model with the highest freezing efficiency, this paper added other sizes on the basis of simulating the freezing process of the whole tuna fish.
The structural parameters are: (1) 0.6 m for a whole fish; (2) 0.3 m for fish segment A (taking the widest part of the tuna); (3) 0.15 m for fish segment B (taking the widest part of the tuna); (4) 100 × 50 × 30 mm for piece A; and (5) 50 × 30 × 10 mm for piece B. The tuna is in direct contact with the surface of the freezing device, as shown in Figure 1. During freezing, cold air flows in from the left side of the model to perform convective heat transfer. After heat exchange, the cold air flows out in other directions under the action of wind pressure.

3.2. Mathematical Models

Since the freezing process of tuna is simulated in this paper, and the process involves a phase transition, the thermophysical property parameters required for the simulation mainly include the specific heat capacity, phase transition enthalpy, phase transition temperature, thermal conductivity, and density of the tuna. For the freezing process of most food products, the phase transition does not occur at a definite temperature, but within a relatively small temperature range. The thermophysical parameters such as the thermal conductivity of tuna will change before and after the phase transition. Therefore, the difference in physical parameters before and after freezing the tuna cannot be ignored. In this simulation, the physical parameters such as specific heat capacity during the freezing process of tuna are calculated using formulas that refer to the water content and phase transition, as provided by Gong et al. [22]. Data, such as enthalpy and phase transition temperature data, were used to determine the physical properties of the tuna, as shown in Table 1.

3.3. Assumptions

On the basis of ensuring the accuracy of the numerical simulation, our model adopts several assumptions: (1) It is assumed that the initial temperature of the model is uniform, the interior is isotropic, and the boundary conditions are constant. (2) It is assumed that the radiation between different products is not included, and only the changes that occur in a single aquatic product during freezing are predicted. Since the radiation temperature between different products in the model was roughly the same, the radiative heat exchange was negligible compared to that resulting from convection. (3) It is assumed that the heat transfer through the wall is negligible, since the wall is considered to be highly insulating due to having a thick insulating layer. (4) It is assumed that moisture evaporation and condensation are negligible [23]. (5) The gas flow and heat transfer predictions are based on the laws of mass momentum and energy conservation, assuming that the gas inside the model follows the ideal gas law.

4. Boundary Conditions and Calculation Methods

Gambit software was used to couple the contact surface between the fluid and tuna in the model. Then, the model was meshed, and the FLUENT software was used to perform numerical simulation calculations on the meshed model. The boundary conditions of the model simulation process were set according to the actual cold chain transportation conditions.

4.1. Inlet and Outlet

In the simulation, the inlet air velocity of the refrigeration equipment was 0 m/s, 1 m/s, 6 m/s, and 20 m/s, and the temperature was −18 °C, −35 °C, and −60 °C. The turbulence is defined by the strength and characteristic of the inlet; here, the turbulence intensity was 5%. To simplify the calculations, the outlet was set as the outflow boundary.

4.2. Boundary Conditions

Since the assumption ignores the radiative heat exchange between aquatic products, the freezing process of tuna is mainly realized by the convective heat exchange between the cold air and the surface of the tuna [24,25].

4.2.1. Internal Heat Conduction Equation

Based on the above premises, the freezing process of tuna is regarded as a variable physical property and unsteady heat conduction problem without an internal heat source [26,27]:
ρ T C p T 𝜕 T 𝜕 t = k T T + h s T T w T e x t

4.2.2. Bottom Equation

The inner wall of the freezing device is insulated to avoid heat transfer. The bottom of the tuna is in direct contact with the wall of the freezing equipment. The temperature of the wall is assumed to be equal to that of the cold air. The temperature of the bottom portion of the tuna was approximately the same as that of the wall:
T d = T e x t
The initial condition is Td = T0, where T0 = 278.15 K is the initial temperature of the tuna.

Surface Convection Equation

When the air speed at the air inlet is 0 m/s, natural convection occurs between the sides of the tuna and the cold air in the freezing device. When there is air velocity at the air inlet, forced convection occurs between the tuna and the air in the freezing device, which is regarded as the third boundary condition [28]:
k ( T ) 𝜕 T 𝜕 n = h s T T w T e x t

Phase-Surface Equation

For the mobile-phase interface, the model satisfies the law of mass conservation and energy conservation:
T s ( x , y , z , t ) = T l ( x , y , z , t ) = T p
λ s 𝜕 T s 𝜕 x + 𝜕 T s 𝜕 y + 𝜕 T s 𝜕 z λ l 𝜕 T l 𝜕 x + 𝜕 T l 𝜕 y + 𝜕 T l 𝜕 z = h ρ d x d t + d y d t + d z d t

4.3. Calculation Method

The simulations were modeled using SolidWorks 2021 and meshed in ANSYS 2021R1. The element order of the model was set to linear, and the inflation option was adjusted for smooth transitions. In order to prove that the simulation results are independent of the number of grids, the calculation was carried out with grids of different sizes. After testing models with elements of different orders of magnitude, it was found that different mesh divisions caused small differences in the simulation. Therefore, it can be assumed that when each model selects the current mesh division method, the calculation results have nothing to do with the mesh number. The mesh quantity details and independence analysis are shown in Table 2.
The left air outlet was set to be the velocity inlet, and the air supply temperature was set according to three different experimental groups: (1) −18 °C (the temperature for general food cold chain transportation); (2) −35 °C (the temperature of the fast-freezing equipment); and (3) −60 °C (the ultra-low temperature freezing temperature required for transporting tuna) [29,30,31].
The air supply speed was set according to three different types of experimental groups: (1) 0 m/s (the air velocity of general direct-cooling refrigeration storage); (2) 1 m/s (the air velocity of air-cooling refrigeration storage and household air-cooling refrigerators); (3) 6 m/s (the air velocity of spiral-type quick-freezing equipment); and (4) 20 m/s (air velocity of impact-type quick-freezing equipment) [32,33].
The temperature at the bottom of the model and the air remains constant, and thus, the air outlets on the other four sides were set as the outflow boundary. The surface of the tuna is the coupling surface, and the software automatically performs the coupling calculation according to the set physical parameters. For the initial setting, the temperature of the tuna was 278.15 K, and the ambient temperature in the freezing area was consistent with the cold air temperature.
The model was set according to the above parameters, the K-ε model and the SIMPLE solution method were selected, and the pressure-based transient solver was used to solve the model. The initial time step was set to 1 s, and the time step was gradually increased to 5 s after the residual curve converged. The convergence accuracy of the energy equation was set to 10−6, and the convergence accuracy of the remaining equations was 10−3. The simulation ended when the maximum temperature of the tuna was detected to be lower than 255.15 K.

5. Results and Discussion

5.1. Accuracy of Numerical Simulations

The freezing curve of the tuna piece at the central point under conditions of −18 °C and a 0 m/s air velocity was measured and compared with the simulated value. From the freezing curve of the tuna, it can be seen that the freezing process is divided into three stages [34,35]. In the first stage, the temperature decreases rapidly to the freezing point, and the freezing rate decreases gradually, which might be due to the decrease in the temperature difference between the tuna surface and the environment after heat exchange occurs. In the second stage, the temperature at the center of the tuna reaches the range of the phase transition, and latent heat is released during the phase transition, which greatly reduces the freezing rate and makes the temperature curve smooth [36,37]. In the final stage, the water inside the tuna freezes, and the cooling rate is lower than that at the initial stage because the specific heat of the ice is lower than that of water [38].
From the comparison between the simulated values and the experimental values in Figure 2, it can be seen that the simulated values at the beginning and end of the freezing process have a higher degree of harmony with the experimental values, and there is a certain difference in their temperature curves as the freezing enters the phase transition period. This may be due to the difference between the phase transition temperature in the simulation and the actual phase transition temperature in the experimental process [39]. After calculation, the correlation coefficient between the simulated freezing time curve and the experimental curve is 0.995, and the error between the simulated and actual curve is 6.67%, which shows that this model can accurately predict the freezing time of tuna.

5.2. Analysis of Velocity and Temperature Field

The air velocity of the refrigeration equipment used in cold chain transportation is not 0 m/s. Therefore, let us take the example of freezing the tuna using cold air at a temperature of −60 °C and an air velocity of 1 m/s. Under this condition, the temperature field inside the tuna and the velocity field of the environment are shown in Figure 3. In the temperature field, the tuna meat near the bottom is cooled faster due to the direct contact between the bottom of the tuna and the cooling device. The temperature difference between different areas due to the above reasons can be improved with a faster air velocity and lower temperature.
When the air hits the 100 × 50 × 30 mm tuna block, the flow velocity in the area blocked by the tuna decreases rapidly, while the part close to the tuna wall gains a higher air velocity due to the increased flow. The cold air flows from left to right along the surface of the tuna, and after the heat exchange, the air flows back to the right, so the air velocity around the right part of the tuna is slightly higher than that around the whole tuna block. At the same time, because the gas density increases in a low-temperature environment, and the bottom of the tuna block clings to the wall of the refrigeration equipment that maintains a low temperature, the air tends to flow downward, and the air velocity on the lower side is higher.

5.3. Analysis of Temperature Curve and Nusselt Number

The freezing process of the whole fish was analyzed under an air velocity of 0 m/s, 1 m/s, 6 m/s, and 15 m/s and a freezing temperature of −18 °C, −35 °C, and −60 °C. Figure 4 compares the temperature curve of the highest point of the tuna temperature and the Nusselt number of the tuna surface under different conditions. According to the comparison chart, it can be found that air velocity and temperature have a significant effect on the freezing time.
When the air velocity during the freezing process is fixed, the lower temperature can help the tuna complete the phase transition more quickly and lower it to the temperature needed for storage. When the freezing temperature decreases, the convective heat transfer coefficient will develop in the opposite direction, and the Nusselt number will increase [40].

5.4. Effect of Air Velocity on Freezing of Tuna of Different Sizes

Using varying air velocities, simulations were carried out on tuna samples of different sizes under different temperature environments. The freezing time required to reduce the maximum temperature of the tuna model to −18 °C under different conditions is shown in Figure 5.
According to the histogram of the freezing time of tuna under different conditions, under the premise of other conditions being equal, the freezing time of tuna without forced convection in the direct-cooling refrigeration equipment is much higher than that when the tuna is frozen in the impact freezing equipment at 20 m/s. The freezing time of natural convection can even reach 20 times that of using an impact freezing device, and this difference is compensated for by the lower-temperature environment. However, the time gap between using 6 m/s air-cooling cold storage and using 20 m/s impact freezing equipment is not extremely different. When the requirements are not strict, after balancing the freezing speed and the energy saving effect, the 6 m/s air velocity can be used to freeze tuna.

5.5. Effect of Temperature on Freezing of Tuna of Different Sizes

Under varying temperatures, simulations were carried out on tuna samples of different sizes under different air velocities. The freezing time required to reduce the maximum temperature of the tuna model to −18 °C under different conditions is shown in Figure 6.
According to the histogram of the freezing time of tuna under different conditions, under the premise of other conditions being the same, when the freezing temperature of tuna is −18 °C, the freezing time when using ordinary cold chain transportation refrigeration equipment is higher than that when using fast-freezing equipment at −35 °C. The freezing time when freezing at −18℃ can be more than three times that of using −35 °C fast-freezing equipment. The time difference between using −35 °C fast-freezing equipment and −60 °C ultralow-temperature freezing equipment is not very different. When the requirements are not strict, after balancing the freezing speed and the energy saving effect, the −35 °C fast-freezing equipment can be used to freeze tuna.

5.6. Effect of Fish Size on Freezing of Tuna

In the case of changing the size of the tuna, simulations were carried out under different air velocity and temperature environments. The freezing time required to reduce the maximum temperature of the tuna model to −18 °C when the tuna is cut into different sizes is shown in Figure 7, and the freezing rate is shown in Table 3. According to the histogram of freezing time for tuna of different sizes, it can be seen that the freezing time of tuna is basically positively correlated with the size of the fish.
But the exception to this is that when the size of the fish block is 0.3 m, the freezing rate is sometimes slightly higher than that when freezing the whole tuna fish of 0.6 m, which may be the combined effect of the heat exchange area and tuna volume. Because the whole tuna fish can be approximated to be a spindle type, after the size of the tuna is changed, the cooling capacity obtained through the heat exchange surface is reduced even though the heat exchange area of the tuna is reduced, and the cooling capacity is reduced to the required temperature despite not being significantly lowered.

6. Conclusions

By taking tuna that is frozen in different freezing equipment as the object, this paper establishes different size models of tuna and discusses the effect of different freezing equipment and conditions on the freezing rate of tuna during cold chain transportation. Through the study of the effects of cold air velocity, freezing temperature, and tuna size on the freezing time using data simulation and experimental verification, the following conclusions have been drawn:
When forced convection refrigeration equipment with an air velocity of 20 m/s is used, the freezing speed is the fastest, and the time required to reduce the maximum temperature of tuna to −18 °C is only one-twentieth of that in an environment without air velocity. Even if an air velocity of only 1 m/s is added in an environment with no air velocity, the freezing efficiency and rate can be significantly improved to reduce the effect of ice crystal formation on the quality of tuna.
The shortest freezing speed can be obtained using ultralow-temperature refrigeration equipment with a temperature of −60 °C, and the time required to reduce the maximum temperature of tuna to −18 °C is only less than a quarter of that of ordinary cold chain transportation equipment at −18 °C. But after evaluating the energy loss and freezing time, the −35 °C flash-freezing equipment can also be used as a good option to reduce the temperature of tuna to −18 °C.
The fastest freezing speed and the shortest freezing time were obtained using smaller tuna pieces. However, the smaller 0.3 m tuna segment has a longer freezing time in some cases than the whole 0.6 m tuna, which is the result of the combined effect of the surface heat transfer area of the tuna and volume of the tuna.

Author Contributions

Supervision, D.Y. and J.X.; writing—review and editing, J.X. and Z.Y.; software, Y.H.; investigation, Y.H.; writing—original draft preparation, Y.H.; methodology, D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Shanghai Professional Technology Service Platform on Cold Chain Equipment Performance and Energy Saving Testing Evaluation] grant number [20DZ2292200] and [Science and Technology Innovation Action Plan of Shanghai Science and Technology Commission] grant number [19DZ1207503].

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meneghetti, A.; Magro, F.D.; Romagnoli, A. Renewable energy penetration in food delivery: Coupling photovoltaics with transport refrigerated units. Energy 2021, 232, 120994. [Google Scholar] [CrossRef]
  2. Nakazawa, N.; Okazaki, E. Recent research on factors influencing the quality of frozen seafood. Fish. Sci. 2020, 86, 231–244. [Google Scholar] [CrossRef]
  3. Moureh, J.; Derens, E. Numerical modelling of the temperature increase in frozen food packaged in pallets in the distribution chain. Int. J. Refrig. 2000, 23, 540–552. [Google Scholar] [CrossRef]
  4. Gandotra, R.; Sharma, S.; Koul, M.; Gupta, S. Effect of chilling and freezing on fish muscle. J. Pharm. Biol. Sci. 2012, 2, 05–09. [Google Scholar] [CrossRef]
  5. Ghoneim, G.A.; Youssef, F.Y.; Ahmed, M.B.; Elkhamisy, B.S.; Anees, F.R. Effects of Frozen Storage on Quality Characteristics of Some fishery Products Processed from Bluefin Tuna and Common Carp. Egypt. J. Food Sci. 2022, 50, 203–221. [Google Scholar] [CrossRef]
  6. Yun, Y.-C.; Kim, H.; Ramachandraiah, K.; Hong, G.-P. Evaluation of the Relationship between Freezing Rate and Quality Characteristics to Establish a New Standard for the Rapid Freezing of Pork. Food Sci. Anim. Resour. 2021, 41, 1012–1021. [Google Scholar] [CrossRef]
  7. Shi, Z.; Zhong, S.; Yan, W.; Liu, M.; Yang, Z.; Qiao, X. The effects of ultrasonic treatment on the freezing rate, physicochemical quality, and microstructure of the back muscle of grass carp (Ctenopharyngodon idella). LWT 2019, 111, 301–308. [Google Scholar] [CrossRef]
  8. Tan, M.; Mei, J.; Xie, J. The formation and control of ice crystal and its impact on the quality of frozen aquatic products: A review. Crystals 2021, 11, 68. [Google Scholar] [CrossRef]
  9. Li, D.; Zhao, H.; Muhammad, A.I.; Song, L.; Guo, M.; Liu, D. The comparison of ultrasound-assisted thawing, air thawing and water immersion thawing on the quality of slow/fast freezing bighead carp (Aristichthys nobilis) fillets. Food Chem. 2020, 320, 126614. [Google Scholar] [CrossRef]
  10. Jia, H.; Roy, K.; Pan, J.; Mraz, J. Icy affairs: Understanding recent advancements in the freezing and frozen storage of fish. Compr. Rev. Food Sci. Food Saf. 2022, 21, 1383–1408. [Google Scholar] [CrossRef]
  11. Zhao, X.; Li, H.; Cui, F.; Wang, J.; Yi, S.; Mi, H.; Lv, Y.; Li, X.; Li, J. Effects of four multi-compound freezing medium on the quality of red drum (Sciaenops ocellatus) during frozen storage. Int. J. Food Sci. Technol. 2022, 57, 4400–4410. [Google Scholar] [CrossRef]
  12. Bian, C.; Yu, H.; Yang, K.; Mei, J.; Xie, J. Effects of single-, dual-, and multi-frequency ultrasound-assisted freezing on the muscle quality and myofibrillar protein structure in large yellow croaker (Larimichthys crocea). Food Chem. X 2022, 15, 100362. [Google Scholar] [CrossRef] [PubMed]
  13. Jiang, Q.; Du, Y.; Nakazawa, N.; Hu, Y.; Shi, W.; Wang, X.; Osako, K.; Okazaki, E. Effects of frozen storage temperature on the quality and oxidative stability of bigeye tuna flesh after light salting. Int. J. Food Sci. Technol. 2022, 57, 3069–3077. [Google Scholar] [CrossRef]
  14. Zhao, Y.; Ji, W.; Chen, L.; Guo, J.; Wang, J. Effect of cryogenic freezing combined with precooling on freezing rates and the quality of golden pomfret (Trachinotus ovatus). J. Food Process. Eng. 2019, 42, e13296. [Google Scholar] [CrossRef]
  15. Zhao, Y.; Ji, W.; Guo, J.; Chen, L.; Tian, C.; Wang, Y.; Wang, J. Numerical and experimental study on the quick freezing process of the bayberry. Food Bioprod. Process. 2020, 119, 98–107. [Google Scholar] [CrossRef]
  16. Wan, T.; Xie, J.; Wang, J.; Zhou, J. Calculation of thermal properties and numerical simulation of freezing time of shrimp. Food Mach. 2018, 34, 106–110. (In Chinese) [Google Scholar]
  17. Wang, J.; Zhang, H.; Xie, J.; Yu, W.; Sun, Y. Effects of Frozen Storage Temperature on Water-Holding Capacity and Physicochemical Properties of Muscles in Different Parts of Bluefin Tuna. Foods 2022, 11, 2315. [Google Scholar] [CrossRef]
  18. Nakazawa, N.; Wada, R.; Fukushima, H.; Tanaka, R.; Kono, S.; Okazaki, E. Effect of long-term storage, ultra-low temperature, and freshness on the quality characteristics of frozen tuna meat. Int. J. Refrig. 2020, 112, 270–280. [Google Scholar] [CrossRef]
  19. Shi, G.; Gao, T.; Qian, X.; Xiong, G.; Shi, L.; Wu, W.; Xin, L.; Yu, Q.; Ding, A.; Li, L.; et al. Effects of different quick-freezing treatments on the quality changes of largemouth bass meat during frozen storage. Meat Res./Roulei Yanjiu 2020, 34, 68–74. [Google Scholar]
  20. Rahman, M.S.; Kasapis, S.; Guizani, N.; Al-Amri, O.S. State diagram of tuna meat: Freezing curve and glass transition. J. Food Eng. 2003, 57, 321–326. [Google Scholar] [CrossRef]
  21. Marra, F.; Zell, M.; Lyng, J.; Morgan, D.; Cronin, D. Analysis of heat transfer during ohmic processing of a solid food. J. Food Eng. 2009, 91, 56–63. [Google Scholar] [CrossRef]
  22. Gong, H.; Tao, R.; Yang, Z.; Liu, M.; Shi, Z.; Wang, S. Thermophysical properties of muscle tissue in different varieties of tuna. Meat Res. 2017, 31, 1–5. (In Chinese) [Google Scholar]
  23. Senguttuvan, S.; Rhee, Y.; Lee, J.; Kim, J.; Kim, S.-M. Enhanced heat transfer in a refrigerated container using an airflow optimized refrigeration unit. Int. J. Refrig. 2021, 131, 723–736. [Google Scholar] [CrossRef]
  24. Agnelli, M.E.; Mascheroni, R.H. Cryomechanical freezing. A model for the heat transfer process. J. Food Eng. 2001, 47, 263–270. [Google Scholar] [CrossRef]
  25. Narsaiah, K.; Bedi, V.; Ghodki, B.M.; Goswami, T.K. Heat transfer modeling of shrimp in tunnel type individual quick freezing system. J. Food Process. Eng. 2021, 44, e13838. [Google Scholar] [CrossRef]
  26. Vabishchevich, P. Numerical solution of the heat conduction problem with memory. Comput. Math. Appl. 2022, 118, 230–236. [Google Scholar] [CrossRef]
  27. Wang, L.R.; Jin, Y.; Wang, J.J. A simple and low-cost experimental method to determine the thermal diffusivity of various types of foods. Am. J. Phys. 2022, 90, 568–572. [Google Scholar] [CrossRef]
  28. Deng, S.; Hu, D.; She, S.; Hong, Z.; Hu, X.; Zhou, F. Freezing Effect of Enhancing Tubes in a Freeze-Sealing Pipe Roof Method Based on the Unsteady-State Conjugate Heat Transfer Model. Buildings 2022, 12, 1373. [Google Scholar] [CrossRef]
  29. Mercier, S.; Villeneuve, S.; Mondor, M.; Uysal, I. Time–temperature management along the food cold chain: A review of recent developments. Compr. Rev. Food Sci. Food Saf. 2017, 16, 647–667. [Google Scholar] [CrossRef]
  30. Kono, S.; Kon, M.; Araki, T.; Sagara, Y. Effects of relationships among freezing rate, ice crystal size and color on surface color of frozen salmon fillet. J. Food Eng. 2017, 214, 158–165. [Google Scholar] [CrossRef]
  31. Rodríguez-Jara, E.; Sánchez-De-La-Flor, F.J.; Expósito-Carrillo, J.A.; Salmerón-Lissén, J.M. Thermodynamic analysis of auto-cascade refrigeration cycles, with and without ejector, for ultra low temperature freezing using a mixture of refrigerants R600a and R1150. Appl. Therm. Eng. 2022, 200, 117598. [Google Scholar] [CrossRef]
  32. Muthukumarappan, K.; Tiwari, B.; Swamy, G.J. Refrigeration and freezing preservation of vegetables. In Handbook of Vegetables and Vegetable Processing; Wiley: Hoboken, NJ, USA, 2018; pp. 341–363. [Google Scholar]
  33. Lv, Y.; Chu, Y.; Zhou, P.; Mei, J.; Xie, J. Effects of different freezing methods on water distribution, microstructure and protein properties of cuttlefish during the frozen storage. Appl. Sci. 2021, 11, 6866. [Google Scholar] [CrossRef]
  34. Leng, D.; Zhang, H.; Tian, C.; Xu, H.; Li, P. The effect of magnetic field on the quality of channel catfish under two different freezing temperatures. Int. J. Refrig. 2022, 140, 49–56. [Google Scholar] [CrossRef]
  35. Cuesta, F.; Sánchez-Alonso, I.; Navas, A.; Careche, M. Calculation of full process freezing time in minced fish muscle. MethodsX 2021, 8, 101292. [Google Scholar] [CrossRef]
  36. Diao, Y.; Cheng, X.; Wang, L.; Xia, W. Effects of immersion freezing methods on water holding capacity, ice crystals and water migration in grass carp during frozen storage. Int. J. Refrig. 2021, 131, 581–591. [Google Scholar] [CrossRef]
  37. Zhao, N.; Yang, X.; Li, Y.; Wu, H.; Chen, Y.; Gao, R.; Xiao, F.; Bai, F.; Wang, J.; Liu, Z.; et al. Effects of protein oxidation, cathepsins, and various freezing temperatures on the quality of superchilled sturgeon fillets. Mar. Life Sci. Technol. 2022, 4, 117–126. [Google Scholar] [CrossRef]
  38. Cartagena, L.; Puértolas, E.; de Marañón, I.M. Impact of different air blast freezing conditions on the physicochemical quality of albacore (Thunnus alalunga) pretreated by high pressure processing. LWT 2021, 145, 111538. [Google Scholar] [CrossRef]
  39. Pattanaik, S.; Jenamani, M. Identifying the cooling heterogeneity and quality decay of Indian mangoes during cold chain export by multiphysics modeling. J. Food Process. Eng. 2023, 46, e14250. [Google Scholar] [CrossRef]
  40. Hu, Y.; Zhang, N.; Wang, H.; Yang, Y.; Tu, Z. Effects of pre-freezing methods and storage temperatures on the qualities of crucian carp (Carassius auratus var. pengze) during frozen storage. J. Food Process. Preserv. 2021, 45, e15139. [Google Scholar] [CrossRef]
Figure 1. Physical models of tuna of different sizes: (a) whole fish; (b) fish segment A; (c) fish segment B.
Figure 1. Physical models of tuna of different sizes: (a) whole fish; (b) fish segment A; (c) fish segment B.
Fishes 08 00470 g001
Figure 2. Comparison of experimental and simulated values of central temperature of tuna.
Figure 2. Comparison of experimental and simulated values of central temperature of tuna.
Fishes 08 00470 g002
Figure 3. Nephogram of the temperature field inside tuna and the velocity field of the environment under specific conditions: (a) nephogram of temperature distribution; (b) nephogram of velocity distribution.
Figure 3. Nephogram of the temperature field inside tuna and the velocity field of the environment under specific conditions: (a) nephogram of temperature distribution; (b) nephogram of velocity distribution.
Fishes 08 00470 g003
Figure 4. Maximum temperature and Nusselt number under different conditions: (a) maximum temperature curve of tuna with 1 m/s air velocity; (b) maximum temperature curve of tuna with −60 °C freezing temperature; (c) Nusselt number with 1 m/s air velocity; (d) Nusselt number with −60 °C freezing temperature.
Figure 4. Maximum temperature and Nusselt number under different conditions: (a) maximum temperature curve of tuna with 1 m/s air velocity; (b) maximum temperature curve of tuna with −60 °C freezing temperature; (c) Nusselt number with 1 m/s air velocity; (d) Nusselt number with −60 °C freezing temperature.
Fishes 08 00470 g004
Figure 5. Effect of air velocity on freezing time of tuna of different sizes at different temperatures: (a) 0.6 m, fish; (b) 0.15 m, segment B; (c) 100 × 50 × 30 mm, piece A; (d) 50 × 30 × 10 mm, piece B.
Figure 5. Effect of air velocity on freezing time of tuna of different sizes at different temperatures: (a) 0.6 m, fish; (b) 0.15 m, segment B; (c) 100 × 50 × 30 mm, piece A; (d) 50 × 30 × 10 mm, piece B.
Fishes 08 00470 g005
Figure 6. The effect of temperature on the freezing time of tuna of different sizes under different air velocities: (a) 0.6 m, whole fish; (b) 0.15 m, fish segment B; (c) 100 × 50 × 30 mm, piece A; (d) 50 × 30 × 10 mm, piece B.
Figure 6. The effect of temperature on the freezing time of tuna of different sizes under different air velocities: (a) 0.6 m, whole fish; (b) 0.15 m, fish segment B; (c) 100 × 50 × 30 mm, piece A; (d) 50 × 30 × 10 mm, piece B.
Fishes 08 00470 g006
Figure 7. The effect of tuna size on freezing time with different air velocities: (a) 0 m/s; (b) 1 m/s; (c) 6 m/s; (d) 20 m/s.
Figure 7. The effect of tuna size on freezing time with different air velocities: (a) 0 m/s; (b) 1 m/s; (c) 6 m/s; (d) 20 m/s.
Fishes 08 00470 g007aFishes 08 00470 g007b
Table 1. Thermophysical parameters of tuna.
Table 1. Thermophysical parameters of tuna.
ParametersBefore FreezingAfter Freezing
ValueCalculated withValueCalculated with
Cp (kJ/(kg·K))3.282Cpl = 0.837 + 3.349ω1.754Cps = 0.837 + 1.256ω
λ (W/(m·K))0.508λl = 0.26 + 0.34ω1.519λs = 2ω + 0.22(1 − ω)
ω0.73
ρ (kg/m3)1050
Tp (K)270.1
TSo (K)264.3
Q (J/kg)176,500
Table 2. Mesh quantity details and independence analysis of different models.
Table 2. Mesh quantity details and independence analysis of different models.
ParametersWhole FishSegment ASegment BPiece APiece B
Mesh quality
Nodes of fish33013990376638282496
Elements of fish14,82018,43517,04210,9241875
Nodes of air12,54015,18013,94914,47315,676
Elements of air68,86379,34774,96972,67781,951
Independent analyses
Order of magnitudesFreezing time at −35 °C and −6 m/s air velocity (s)
10425,28525,48718,7932954574
10526,68527,25420,9053408702
10627,45428,05721,5263479731
Table 3. The freezing rate of different tuna sizes and freezing environments.
Table 3. The freezing rate of different tuna sizes and freezing environments.
Freezing ConditionFreezing Rate
TemperatureAir VelocityWhole FishSegment ASegment BPiece APiece B
−18 °C0 m/s0.0200.0150.0240.5705.255
1 m/s0.1380.1360.1850.9987.901
6 m/s0.3090.3110.4102.29810.805
20 m/s0.4670.4270.5733.84515.894
−35 °C0 m/s0.0710.0550.0901.92118.084
1 m/s0.4490.4390.6143.87425.509
6 m/s0.9240.9051.1807.23835.138
20 m/s1.3331.2011.58111.48450.137
−60 °C0 m/s0.1290.1020.1623.46432.372
1 m/s0.7740.7401.0206.75844.931
6 m/s1.5131.4551.85412.21161.361
20 m/s2.1381.9122.44719.04887.784
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huo, Y.; Yang, D.; Xie, J.; Yang, Z. Simulation Analysis and Experimental Verification of Freezing Time of Tuna under Freezing Conditions. Fishes 2023, 8, 470. https://doi.org/10.3390/fishes8090470

AMA Style

Huo Y, Yang D, Xie J, Yang Z. Simulation Analysis and Experimental Verification of Freezing Time of Tuna under Freezing Conditions. Fishes. 2023; 8(9):470. https://doi.org/10.3390/fishes8090470

Chicago/Turabian Style

Huo, Yilin, Dazhang Yang, Jing Xie, and Zhikang Yang. 2023. "Simulation Analysis and Experimental Verification of Freezing Time of Tuna under Freezing Conditions" Fishes 8, no. 9: 470. https://doi.org/10.3390/fishes8090470

APA Style

Huo, Y., Yang, D., Xie, J., & Yang, Z. (2023). Simulation Analysis and Experimental Verification of Freezing Time of Tuna under Freezing Conditions. Fishes, 8(9), 470. https://doi.org/10.3390/fishes8090470

Article Metrics

Back to TopTop