I. Introduction
Beer is one of the most consumed beverages in the world and the most popular among alcoholic ones [1, 2, 3], it has been present in human life for thousands of years [1, 3, 4, 5, 6], and its popularity has grown throughout history. Because of this, the focus on improvements in its production becomes essential.
In the production of beer, the aim is to improve the methods to obtain a quality product, from the preparation of an established formula, to the treatment of the drink in each of the stages. Rodman [7] recognizes fermentation as the most important process when brewing beer, therefore, it is sought to improve the production times of this stage and generate a better yield [8], [9]. Specifically, through improvements in cooling, after the wort boil.
The cooling of the hot wort must be carried out as quickly as possible to avoid the formation of microorganisms or bacteria [1, 3, 5, 10-11], regarding the range of temperatures, Agulló [12] clarifies that it goes from 15 °C to 25 °C, when the yeast is alive. This process requires a lot of time and energy, due to this, shortening its duration and cost implies enormous potential savings [8]; generally, this energy saving is intended through improvements in heat exchangers.
In the brewing industry the cooling process is commonly carried out using plate heat exchangers [13], due to the fact that they are easy to assemble and clean; also, since they are easily opened, they allow the inspection of the entire heat transfer surface [14]. During the heat exchange, incrustations or fouling are generated [15], which result from the caramelized layer of sugars on the surface [5] that reduces the efficiency of the heat exchanger [16].
Qiu and Qiu [17] made a machine for the cooling of craft beer through tubes assembled in a helical way, where a fresh and refrigerated product can be obtained directly from the tap. In turn, they affirmed that conventional machines have the problem of being slow or distributing fluids with low speeds, which is solved by using a system of tubes assembled in a helical way in different distributions, the first creating an internal coil (beer flow) and an external coil (cooling flow). In the other case, both coils were set collinearly, one on top of the other, where the second case obtained the best efficiency.
On the other hand, Liu et al. [18] clarify that all the favorable parameters and, therefore, the efficiency, improve when reducing the diameter of the spiral and increase the pitch of the helical tubes. This by comparing the hydrodynamic performance of a multi-tube helical heat exchanger (MTTHC), between experimental results and CFD simulation. The validation of its data is obtained through an error of 9.81% in pressure and 7.6% in temperature, with respect to the experimental investigation.
Wang et al. [19] carried out a comparison of the obtained experimental results using CFD software for a helically assembled trilobed tube heat exchanger (HCTTT), the tests were carried out with Reynolds number ranges of 2000-40000. Errors were obtained of less than 5% in heat parameters (Q), convection heat transfer coefficient (h), Nusselt number (Nu), and friction factor (f), which makes the research credible. The experiment includes an analysis of HCPT, HCET, and HCTT exchangers with the same diameter of 6 mm, and it was obtained that the thermal performance is higher in the HCTTT up to a maximum of 1.22 times, which increases proportionally to the Reynolds.
El-said et al. [20] compares the results of a helical plate heat exchanger (HPHE) with a Reynolds range of 2000-7000 in hot fluid and 1400-7000 in cold fluid, as well as a passage radius of 0.24-1.31, obtaining that the efficiency varies from 0.1758 to 0.6756 in the established intervals. The error maintained is 13% compared to the experimental data, in addition, it affirms that all the parameters improve, except for the pressure, where there is a slight drop.
Hatef and Ranjbar [21] observed, in their study of the effect of geometric parameters on heat transfer and hydrodynamics characteristics of the helical exchanger, that the heat transfer coefficient increases in relation to the internal diameter of the tube and bending ratio, and due to the decrease in the pitch of the coil. This study is developed by validating the numerical data obtained with the experimental model, using a double helical tube exchanger for their study. As for the Nusselt number, it is affected by the increase in the radius of the exchanger coil [22], resulting in a decrease in the global heat transfer coefficient.
The present research aims to establish the appropriate parameters to reduce the temperature of the beer wort in the shortest possible time, through a CFD simulation, considering water as a cooling fluid, in a helical tube and shell heat exchanger.
II. Methods and Materials
The brewing process is composed of different stages, such as preheating, mashing, hops, cooling, heat exchange, fermentation and, finally, gasification. The heat exchanger is the central part of this analysis and where the fluid is subjected to different changes. The process begins with the preheating of the water, so that it meets the input conditions prior to the maceration; at this point, the grains are added and are taken to a cooking process that allows to obtain the wort, which represents the object of study. This wort will be cooled through the heat exchanger, ensuring that the cooling takes place in the shortest possible time and that the gasification process can continue, where the beer preparation ends.
The fluid with the highest temperature is the wort, which enters the heat exchanger at a temperature of 70 °C in the helical geometry tube, while the water used for cooling enters the shell at 15 °C. This process is shown in Figure 1, the outlet temperature of the wort exchanger should be 25 °C for high fermentation ale of beer type.
A. System Description
The research has foreseen the design of a shell and helical tube heat exchanger. The material used in the design of the shell is commercial steel, while DHP type L copper is considered for the helical tube.
The shape of the helical tube is described with two geometric parameters, the pitch p and the radius of the helix Rh, the selection of the shell exchanger is due to the fact that the hot fluid wort does not contact the cooling fluid that circulates through the geometry of the shell.
One of the most important factors in the fluid cooling process is the correct selection of the geometry of the heat exchanger and especially of the helical tube. The geometric dimensions selected for the system are specified in Table 1.
B. Description of Working Fluids
To generate a system that can be applied, the study parameters are established through the type L copper and steel commercial materials. The thermodynamic properties of the materials are shown in Table 2, for calculation purposes the thermodynamic properties are considered as pure copper, since DHP copper has a purity greater than 99 %.
On the other hand, Table 3 denotes the thermodynamic properties of hot and cold working fluids, such as beer wort and water, respectively.
Through this analysis, it is intended to obtain parameters of temperature and pressure at the outlet of the wort in the heat exchanger, the mesh shown in Figure 2 was obtained after establishing its application with several simulations as valid, and a mesh quality referring to the obliquity and orthogonality.
The convergence of the skew mesh is shown in Figure 3, with an average of 0.20981. On the other hand, the mesh convergence in orthogonality presents an average of 0.82443, with a minimum of 0.11565, as shown in Figure 4. The statistics obtained from the meshing of the heat exchanger correspond to 1696103 nodes with 3227323 elements.
C. Simulation Model
In helical tube modeling, the governing equations that correspond to the ANSYS Fluent 2019 software are considered, within which the k-epsilon RNG (k-ε) turbulence model can be seen.
The fundamental continuity equations are based on Navier Stokes' concept of conservation of mass, momentum, and energy, obtained through numerical methods and expressed in equations 1, 2, 3, respectively.
To calculate the equilibrium of moments, it is necessary to obtain the shear stress generated on the fluids, through equation 4.
Due to the existence of curves in the studied model, the epsilon RNG kappa model is used. In equation 5 the turbulence model for kappa kinetic energy is presented.
On the other hand, equation 6 represents the epsilon diffusivity in the turbulence model.
Within the calculation of helical tubes, centrifugal forces that influence the calculation factors are generated [24], specifically the Reynolds number, therefore, a relationship must be found to correct the deviation that exists due to the aforementioned forces. This relationship is called Dean's number and is obtained through equation 7.
III. Results and Discussion
The simulation starts under the k-epsilon RNG turbulence model, with mass flows of 0.7 kg·s-1and 7.5 kg·s-1 for the wort and water, respectively, and with inlet conditions of 101325 Pa pressure and a temperature of 70 ºC. Since it is a shell exchanger, 18 ºC is set as the outside temperature, corresponding to the water flowing through the shell. Figure 5 shows the behavior of the helical tube with the indicated parameters.
After performing the first simulation, an outlet temperature of 33.43 °C is obtained, which does not correspond to the temperature expected to start the fermentation process.
In the search for the correct parameters for the heat exchanger established with the initial dimensions, the mass flow is reduced to 0.5 kg·s-1, maintaining the temperature and inlet flow for the water. The behavior of this simulation is shown in Figure 6.
The results obtained after the second simulation show a significant improvement, generating a wort outlet temperature of 28.54 °C, satisfactorily approaching the expected result.
Regarding the parameters for the third simulation, an inlet mass flow was established for wort and water, of 0.2 and 7 kg·s-1, respectively; also, the water temperature is kept at 18 °C. The generated behavior is shown in Figure 7.
The results obtained from the third simulation show that the outlet temperature of the wort dropped dramatically to 18.42 °C, exceeding the expected limit of 25 °C.
The input parameters for simulation 4 are 2.5 kg·s-1 and 8 kg s-1for the input mass flow of wort and water, respectively, while the water temperature does not change compared to the previous simulations. The behavior of the wort in the helical tube for the fourth simulation is shown in Figure 8. From the fourth simulation, a wort outlet temperature of 28.44 °C is obtained, which rises compared to simulation 3. However, it falls away from the desired result again.
From the fifth simulation, the inlet temperature of the cooling water is considered an important variable, which is reduced to 15 °C, and the mass flows are considered just as in the first simulation: 0.7 kg·s-1 for wort and 7.5 kg·s-1 for water. The behavior of the wort in the helical tube for these parameters is shown in Figure 9.
With the new established parameters, a wort outlet temperature of 30.72 °C is obtained.
Continuing with the process of emulating the parameters of the initial simulations with the new water inlet temperature, the behavior of the tube is set at 0.5 kg·s-1 for wort, and 7.5 kg·s-1 for water. this behavior can be seen in Figure 10.
Simulation 6 shows a wort outlet temperature of 27.84 °C. The outlet temperature obtained with the simulation meets the requirement for the ale brewing process, a product established for the simulation. Although the change in the geometry of the helical tube [21] allows to obtain good results in terms of heat transfer, it is also necessary to consider the variations of the input mass flow.
IV. Conclusions
Considering the various parameters used for the simulations, better results for reducing the outlet temperature are obtained by changing the mass flow of the wort, where there is a variation of up to 15.16 °C when varying it from 0.7 to 0.2 kg·s-1. This by varying the water inlet temperature, where there is a maximum variation of 2.71 °C and a minimum of 0.01 °C.
The mass flows of 0.2 and 7 kg·s-1 of the wort and water, respectively, are the ideal flows for the operation of the heat exchanger. The error obtained with respect to the mathematical calculations corresponds to 2.71 %, between an estimated temperature of 25 °C and one obtained of 25.67 °C.
Future investigations of helical tube and shell heat exchangers should focus on the behavior and results when varying parameters such as the diameter and pitch of the propeller, since this study did not make considerable modifications to the initial dimensioning.