Next Article in Journal
Artificial Noise Injection and Its Power Loading Methods for Secure Space-Time Line Coded Systems
Previous Article in Journal
Mimicking Anti-Viruses with Machine Learning and Entropy Profiles
Previous Article in Special Issue
Heat Transfer and Entropy Generation Abilities of MWCNTs/GNPs Hybrid Nanofluids in Microtubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Generation of Forced Convection during Melting of Ice Slurry

by
Beata Niezgoda-Żelasko
Institute for Thermal and Process Engineering Cracow University of Technology, al. Jana Pawła II 37, 31-864 Kraków, Poland
Entropy 2019, 21(5), 514; https://doi.org/10.3390/e21050514
Submission received: 31 March 2019 / Revised: 17 May 2019 / Accepted: 19 May 2019 / Published: 21 May 2019
(This article belongs to the Special Issue Entropy Production in Turbulent Flow)

Abstract

:
This paper looks at entropy generation during ice slurry flow in straight pipes and typical heat exchanger structures used in refrigeration and air-conditioning technology. A dimensionless relationship was proposed to determine the interdependency between flow velocity and the volume fraction of ice, for which the entropy generation rates were at the minimum level in the case of non-adiabatic ice slurry flow. For pipe flow, the correlation between the minimum entropy generation rate and the overall enhancement efficiency was analyzed. As regards heat exchange processes in heat exchangers, the authors analyzed the relationship between the minimum entropy generation rate and the heat exchange surface area and exchanger efficiency.

1. Introduction

Ice slurry-based refrigeration is among the most state-of-the-art refrigeration technologies in indirect installations with and without cold storage. It is also used as a process supporting cold distribution in ice harvester systems. Ice slurry is a mixture of water ice crystals and water or water and a freezing-point depressant (salt, glycol, alcohol, etc.). The ice particle diameters most often amount to 10−4–0.5 mm. The mass fraction of ice crystals in the ice slurry, which may be transported in practice, do not exceed 30%. As an environmentally neutral and natural heat carrier, ice slurry demonstrates a refrigeration potential commensurate to that of refrigerants. Ice slurry flow may be accompanied by phase segregation, which leads to a change in the mean density values and in the value of the dynamic coefficient of fluid viscosity, consequently leading to higher flow resistance values [1]. The heat transfer process in the ice slurry is accompanied by micro-convection of solid particles, resulting in an increase in thermal conduction coefficients and an enhancement of the heat transfer process. Ice slurries differ significantly from single-phase heat carriers in terms of their rheological properties. Ice slurries are non-Newtonian fluids [2,3,4,5,6,7,8,9,10]. The phenomena associated with flow resistance and heat transfer processes involving ice slurries have been extensively discussed in articles, such as [2,3,4,5,6,7,8,9,11,12]. A significant operational problem of ice slurry-fed systems is the elimination of phase segregation and ensuring a homogenous flow. According to [13,14], regardless of the ice content, homogenous flow of ethanol-based ice slurry is possible at a minimum flow velocity of 0.54 ms−1 [13] and 0.75 ms−1 [14] for 0.02 and 0.024 m diameter pipes, respectively. Flow velocity cut-off values for homogenous and heterogeneous flow in d = 0.02 m pipes for 10.5% ethanol-based ice slurry, calculated in [15] on the basis of a solid particle distribution model according to Kitanowski [4,16,17,18], for a 5% and 25% volume fraction of ice amount to 0.68 ms−1 and 0.15 ms−1, respectively.
While dynamic effective viscosity coefficients of slurries are higher than dynamic viscosity coefficients of carrier fluids (and they increase as the mass fraction of ice grows), the slurry flow resistance values are not always higher than those of carrier fluids at the same flow velocities. For an ice content of xs ≥ 10% and the same flow velocity, ice slurry flow resistances may be smaller than the flow resistances of the carrier fluid, provided that the type of flow of the carrier fluid and of the ice slurry is different [3,19,20,21]. In the case of homogeneous flow, these phenomena can be explained by the effect of the absorption of turbulence kinetic energy by ice particles (flow laminarization) and the shifting of the moment of transition to the turbulent flow area by slurries with high mass fractions of ice and higher Reynolds number than in the case of the carrier liquid. For heterogeneous flow, which is favoured by low ice mass fractions and larger dimensions of solid particles, it is difficult to interact with the carrier fluid and solid particles that accumulate in the top part of the tube. This effect favors an increase in turbulence of the carrier liquid and an increase in flow resistance. The character of motion changes at lower speeds and lower Reynolds numbers.
Heat exchange processes involving ice slurries are determined by several phenomena. Solid particle micro-convection associated with moving ice slurries significantly improves the thermal conduction coefficient values with respect to those determined for non-moving slurries, for example, using the Maxwell–Tareff equation [22,23]. An analysis of the Charunyakorn [24] equation implies that higher fluid shear rates correspond to higher thermal conduction coefficient values. The ice melting process is characterized by significantly higher specific thermal efficiency as compared to heat exchange during single-phase flow. These phenomena make ice slurry heat transfer coefficients higher than those of single-phase refrigerants. For example, for ice slurries based on an aquatic solution of ethanol, in the laminar flow area, the presence of solid particles leads to a 2–5.7-fold increase in the heat transfer coefficient with respect to the heat transfer coefficients of the carrier fluid (ethanol xai = 10.6% at −4.5 °C) [2,12,25]. On the other hand, there are some velocities at which for ice slurries with an ice content of xs > 10%–20%, the heat transfer coefficients are lower than those of the carrier fluid. In the turbulent flow area, the ice content is not as significant in terms of its impact on the heat transfer coefficient value, although the values of these coefficients are about 20%–30% higher than the heat transfer coefficients of the carrier fluids [2,12,25].
Considering the maximization of the heat transfer coefficients and the minimization of flow resistance values, for ice slurry-fed installations, the di, w, xs parameters should be selected so that the ice slurry flow is homogenous and laminar, with heat transfer coefficients exceeding those of the carrier fluid. Laminar flow favors lower flow resistances, which are always higher for ice slurries in the turbulent flow area than the flow resistances of the carrier fluid with the same flow velocity.
The large mass fraction of ice makes it possible to limit the mass flux of the refrigerant and the heat exchange surface area of ice slurry-fed heat exchangers. Without a doubt, ice particles enhance the heat transfer process. The efficiency of ice slurry use as a refrigerant can be determined on the basis of the criteria commonly used in heat transfer. The obvious indicators of the efficiency of various methods used to improve heat exchange conditions include the heat transfer coefficient enhancement ratio (εh = h/href) and the flow resistance increase ratio ( ε Δ p = Δ p / Δ p r e f ). The mutual relationship between εh and εΔp is expressed as their quotient (εh−Δp = εh/εΔp) [26], or the so-called overall enhancement efficiency, defined as εNu−f = (Nu/Nuref)/(f/fref)1/3 [27,28]. The use of the second principle of thermodynamics in design is based on an analysis of the entropy generation rate, which enables a quantitative determination of the irreversibility of the analyzed thermal and flow process. It is common practice to use the entropy generation rate for the thermo-dynamic assessment of thermal and flow processes and as an optimization criterion, for example, in the design of heat exchangers. Also used are complex criteria for the assessment of efficiency—a combination of the performance evaluation criteria and the entropy generation minimization (EGM) criterion [29]. In [30], the authors propose a method for the analysis of the entropy generation rate in the design of single-phase thermal processes, heat storage systems or various thermal cycles. The method proposed by Bejan [30] has been successfully used by other authors to perform a thermodynamic assessment of thermal processes during single-phase flows [31,32], multi-phase flows [33,34,35,36,37,38], heat exchange involving non-Newtonian fluids [32,39,40] and design of thermal systems and exchangers [40,41,42,43,44,45,46,47,48].
In the boiling process, an analysis of the entropy generation rate was used, for example, to determine the preferred scope of use of plain and micro-finned tubes [36] or to optimize the geometrical parameters of the pipes of an R22 evaporator [37]. An analysis of the entropy generation rate makes it possible to determine the upper limit of the mass flux density used in micro-finned pipes. The authors of [36] demonstrated that the entropy generation minimization criterion cannot be the only optimization criterion in the design of heat exchangers and that multi-criterion optimization of heat exchanger geometry should be preferred. In this context, entropy generation and material cost minimization was used in [42] to design crossflow plate heat exchangers. The minimization of the dimensionless number of entropy generation made it possible to determine the optimum geometrical parameters for finned heat exchanger plates. Heat flux maximization and entropy generation minimization in the design of fluid coolers and air-cooled condensers was the topic discussed in [41]. The authors stated that while in single-phase heat exchangers entropy generation minimization leads to the maximization of thermal efficiency, this condition is not met in condensers.
Paper [49] is a more general publication about the usefulness of determining entropy generation rates in thermal and flow-related numerical calculations. Analyzing, among other things, the melting process, the authors identified a correlation between the correct determination of the entropy generation rate and numerical errors, convergence criteria and the selection of the time-step. In another study [45], the authors presented the use of the VEG (virtual entropy generation) method in measurement technique to reduce measurement errors in the studies of heat exchangers.
The issues concerning entropy generation in thermal and flow processes with ice slurry were discussed in the studies [38,39,40,46].
The entropy generation rate in the heat exchange process involving ice slurries was discussed in [38], where two types of aquatic slurries with PCM elements and multi-walled carbon nanotubes (MCNT) were compared. An analysis of the entropy generation rates found that for small pipe diameters and low heat flux densities, the use of PCM was preferred. Greater pipe diameters and higher heat flux densities, especially in the laminar flow area, preferred the use of aquatic slurries with the addition of MCNT. In the turbulent flow area, it is preferable to use water as the heat carrier, as it generates the lowest entropy generation rates. In another study [39], the analysis of entropy flux generation was limited to the case of laminar, isothermal flow of ice slurry which was an Ostwald–deWaele fluid. The authors of the paper [46] used the second principle of thermodynamics for the analysis of heat balance and entropy in a heat exchanger supplied with ice slurry of an aqueous solution of ethanol. In this paper, the effect of flow processes and heat transfer on the coolant and ambient side is not taken into account. The assumptions adopted by the authors implied the use of slurries with low mass fractions of ice. A comprehensive approach to the analysis of entropy flux in thermal processes with ice slurry has been presented in a previous publication [40]. The authors analyzed the entropy flux generated in a spiral heat exchanger supplied with ice slurry. It was indicated that it is possible to determine optimal Dean number values, for which the entropy flux is minimal, and the maximum thermal efficiency is achieved for specific mass fractions of ice.
An analysis of the above papers suggests a very broad use of entropy generation rates for the assessment of diverse thermal [32,35,38,43], and flow processes, in which flow parameters [31,33,34,36] and geometrical parameters may be optimized [37,41,42]. The determined entropy generation rates refer to single processes [36,38], or to devices where processes involving various agents take place [41,42,43,44]. Local or global values of the entropy generation rate may be determined analytically [32,33,34,35,38,39,42,44,46], or numerically [31,36,37,40,49,50].
The use of ice slurry in refrigeration processes is a heat transfer enhancement method applied in refrigeration and air-conditioning technology. Ice slurry fed into a heat exchanger is characterized by both mass flux and mass fraction of ice, which determine the heat transfer process intensity and momentum exchange, inducing various entropy generation rates. This article attempts to identify the range of velocity and mass fraction of ice values which ensure the minimum entropy generation rate during the flow of ethanol-based ice slurry, being a non-Newtonian Bingham’s fluid, through straight 0.01–0.02 m pipes. In addition, the study contains aconducted an analysis of the total entropy generation rate in lamelled air coolers and plate fluid coolers in the context of the optimization of the mass fractions of ice and mass fluxes of the ice slurry. In this case, the analysis focused on the relationships between ice slurry’s flow parameters, the heat exchange surface area, the heat exchanger’s efficiency and the total entropy generation rate on the side of the refrigerant and the cooled medium. New aspects presented in the paper refer to a comprehensive approach to thermodynamic analysis of thermal-flow processes with ice slurry, both for flow in straight tubes as well as in heat exchangers. The analysis of minimum total entropy fluxes in heat exchangers allowed to verify the usefulness of the results of the minimization of the entropy flux of ice slurry in the conditions of forced convection of ice slurry in straight tubes. Additionally, the study compares the results of determination of optimal parameters of flow velocity and mass fractions of ice according to the criterion of minimum entropy and the condition of εNu−f > 1. The innovative element presented in the paper was proposing a dimensionless relationship to identify the interdependency between the flow velocity and the mass fraction of ice for which the entropy generation rate was at its minimum level.

2. Entropy Generation Rate during Flow in a Straight Pipe

2.1. Entropy Generation during the Melting of Ice Slurry; Flow with Phase Separation

Assuming a homogenous single-dimensional flow of the ice slurry mass flux m ˙ I S in a horizontal pipe (Figure 1), the first and second principle of thermodynamics for the dl test section may be described using Equations (1) and (2), respectively:
d i I S = T I S d s I S + ν I S d p I S
d S ˙ l d l = d ( m ˙ S s s + m ˙ F s F ) d Q ˙ T w
In the test section, the ice slurry entropy changes as per Equation (3):
m ˙ I S d s I S = d ( m ˙ S s s + m ˙ F s F )
Using Equations (1) and (3), Equation (2) can be transformed into (4):
d S ˙ l = d Q ˙ d l ( 1 T w 1 T I S ) + m ˙ I S T I S ν I S ( d p I S d l )
Writing the Fourier’s equation for the ice slurry heat transfer process (5) and using the relationship d Q ˙ = q ˙ d A = q ˙ P d l , Equation (4) can be transformed into (6):
Δ T = T w T I S = q ˙ α I S
d S ˙ l = q ˙ 2 P α I S T I S 2 + q ˙ T I S + m ˙ I S T I S ν I S ( d p I S d l ) = d S ˙ l H + d S ˙ l Δ p
Equation (6) makes it possible to determine the entropy generation rate for a unit length of the pipe. Equation (6) is true for water-based ice slurry with a constant melting point. For ice slurries based, for example, on a 10.6% aquatic solution of ethanol, the temperature glide with complete melting of ice in a 30% ice slurry amounts to 2.2 K. Therefore, it was assumed that in (6), TIS corresponded to the mean temperature of the ice slurry, T I S = T ¯ I S .

2.2. Ice Slurry Flow and Heat Transfer in Pipes

The determination of the entropy generation rate from Equation (6) for non-adiabatic flow of ethanol-based ice slurry requires the identification of both the pressure drop and the heat transfer coefficients for specific heat transfer conditions. The ice slurry under consideration is a Bingham fluid [3,19] and its flow can be treated as a generalized flow of a non-Newtonian fluid. Equation (7) introduces the generalized form of the Reynolds number:
Re K = ρ I S w 2 n * d h n * 8 n * 1 K *
By using Equation (7), Fanning’s factors for technically plain pipes can be found for the laminar and turbulent flow area on the basis of Equations (8) and (9), respectively:
c f = 16 Re K
c f = 0.079 Re K 0.25
Ultimately, ice slurry flow resistance values are calculated using Equation (10).
Δ p = c f 2 l ρ I S w 2 d h
Equation (7) makes it necessary to determine a characteristic flow-behavior index n* and consistency index K*. For Bingham fluid and an arbitrary geometry of cross-section, these parameters are found from equations included in Table 1 [19].
For the generalized non-adiabatic flow of ice slurry, the heat transfer coefficient may be calculated for laminar and flow areas from Equations (11) and (12), respectively [25]:
Nu L = B 1 ( Gz K ) B 2 ( Δ x s K F 100 ) B 3 ( d s d h ) B 4 ( K * ( T I S ) K * ( T w ) ) B 5
Nu T = B 1   ( Pe K ) B 2 ( d s d h ) B 4
The constants Bi in Equations (11) and (12) are shown in Table 2 [25]. Equations (7)–(12) are the result of own experimental studies carried out for ethanol ice slurry, for which measurement stations, the course of experiments, measurement accuracy and results were discussed in detail in studies [2,12,19,25,55]. The experimental scope of the applicability of Equations (7)–(12) is provided in Table 2.

2.3. Results of Calculations

The entropy generation rate was studied for forced convection during melting a 10.6% ethanol solution-based ice slurry with a mass fraction of ice of 0 <xs < 30%. Calculations were made for ice slurry flow through straight pipes with diameters di = 0.01, 0.016, 0.02 m, mass flux of ice slurry 0.04 m ˙ I S 0.62 kgs−1 and heat flux density 2 q ˙ 10 kWm−2. Figure 2 presents the change in the total entropy generation rate as a function of the flow velocity and the mass share of ice during ice slurry flow through a di = 0.016 m pipe with a constant heat flux density q ˙ = 8 kWm−2.
Figure 2 suggests that in the laminar flow area, the lowest entropy generation rates correspond to the maximum fraction of ice in the ice slurry. However, let it be noted that in the entire flow area, the lowest entropy generation rates correspond to low mass fractions of ice xs < 10% and the minimum value of the entropy generation rate is recorded in the turbulent flow area. This is related to the mutual relationship between entropy generation in thermal and flow processes (Figure 3).
In the studied range of flow velocities 0 ≤ w ≤ 3 ms−1 for an ice content greater than xs > 21.9%, the flow is laminar in nature, and the ice slurry flow resistance values are moderate, while the presence of solid particles favors the highest values of heat transfer coefficients [2,12]. In this case, entropy generation is determined by the irreversibility of the thermal processes. For low mass fractions of ice, obtaining suitably high values of heat transfer coefficients depends not only on the presence of solid particles, but also on the appropriate value of flow velocity. As a consequence, this leads to a change into the turbulent flow area, where, in turn, flow resistance values have a significant impact on the induction of the total entropy generation. It also needs to be noted that large mass fractions of ice (xs > 20%) are typically associated with a relatively large range of velocities, for which the total entropy generation rate is at its minimum. This results from the dominant impact of solid particles and phase change, but not flow velocity, on heat transfer coefficients [2]. The calculations made for various heat flux densities and pipe diameters made it possible to determine for specific mass fraction of ice, the flow velocities for which the total entropy generation rate was at its minimum.
The analysis of Equations (6)–(12) and calculation results presented in Figure 2; Figure 3 indicate that in case of non-adiabatic ice slurry flow, the generated entropy flux is determined by flow velocity, mass fraction of ice, heat flux density and tube diameters. It should also be noted that the mass fraction of ice, flow velocity and diameter determine the nature of the flow of the agent. With the increase in the Reynolds number, the effects of particulate matter on the heat transfer coefficients decreases and a greater influence of flow phenomena on the generated entropy fluxes is observed. Therefore, in the turbulent flow area, the entropy flux depends more on flow velocity and diameter than mass fraction of ice and heat flux density. Small tube diameters imply higher values of heat transfer coefficients and low values of the entropy flux generated by the heat exchange process. In general, lower tube diameters correspond to lower total values of the entropy flux. The results presented in Figure 4a,b indicate that in the laminar flow area, the effect of tube diameter on entropy flux is particularly significant: 1.5 S ˙ d i / S ˙ d i o 6.2 for 1.6 d i / d i o 2 . In the turbulent flow area where the values of the Bejan numbers are high (Be > 0.5) the S ˙ d i / S ˙ d i o ratio stabilizes at S ˙ d i / S ˙ d i o = 1.34–1.4 and S ˙ d i / S ˙ d i o = 2.7–3.7 for d i / d i o = 1.6 and d i / d i o = 2, respectively.
Figure 5 presents the relationship w(xs)Smin for various heat flux densities and the flow of ice slurry through a pipe with a diameter of di = 0.016 m. The values of w, xs, determined analytically, were validated using experimental values.
The experimental values of w and xs were determined by calculating, at each measurement point, the entropy generation rate (Equation (6)) using the measured values of heat transfer coefficients and the flow resistance values [2,25]. Next, for the given values of xs, the measurement points were chosen at which S ˙ = S ˙ min and the non-monotonic course of the relationship w(xs)Smin was associated with an area of transition between laminar and turbulent flow. The plotted theoretical w c a l S ˙ min ( X V ) curves provide a qualitatively correct description of the course of analogous curves obtained directly from the measuring points. The mean relative difference between the experimental and analytical flow velocities w S ˙ min ( X V ) does not exceed 7.6%. The maximum difference between the measured velocities and the calculated velocities concerned the case of heat flux density q ˙ = 2 kW and was lower than 30%.
Figure 5 shows the critical values of xsC, and wC found on the basis of the criterion for motion nature change valid for the ice slurry (Equation (13)) [54]:
Re B C = 10000 He ( d s d i ) 0.25 1.25 He + 334.9
All analyzed parameters (xs, di, w, q ˙ ) are independent, but they affect the heat transfer coefficients and flow resistance, and thus the entropy flux and its minimum value, in various ways. Charts similar to the one presented in Figure 5 allow the selection of parameters xs, (Xv Table 1) w, di, q ˙ in such a way that the entropy flux generated by the ice slurry is minimal. The results presented in Figure 5 and similar results for tube diameters di = 0.01 m and di = 0.02 m can be presented in the form of criteria relationship (14). Correlation (14) defines the relationship between the Reynolds number and dimensionless numbers taking into account the effect of geometric (KqX) and flow parameters (Kq, KqX) on the generated minimal entropy flux.
Re I S S min = K q + K q X
where
K q = C 1 ( q ˙ o q ˙ q ˙ o ) C 2
K q X = exp ( C 3 + C 4 ln ( d i d i o ) + C 5 q ˙ r G ˙ o ( X v 100 ) )
The parameters Ci, dio, q ˙ o , G ˙ o are given in Table 3.
Equation (14) makes it possible to determine the optimum velocity values for the condition S ˙ = S ˙ min , for the assumed mass fractions, heat flux densities and pipe diameters with the maximum relative error of 15% (mean error: 1.2%).
Figure 6 presents the relationship between the Reynolds number for the ice slurry (ReIS calculated as for Bingham’s fluid) and q ˙ X v / r .
Figure 7 presents a comparison of the cut-off values of xs, w both for S ˙ = S ˙ min and εNu−f = (Nu/Nupt)/(f/fpt)1/3 = 1, for ice slurry heat transfer during ice slurry flow through di = 0.01 and 0.02 mpipes with constant heat flux density values of q ˙ = 5, 10 kWm−2. The range of xs and w values, for which εNu−f > 1, means that while the same power is required to transport the refrigerant, the use of ice slurry over the use of the carrier fluid as the refrigerant translates into better thermal efficiency. The results shown in Figure 7 suggest that the minimum entropy generation rate condition generally shifts the scope of ice slurry use towards higher flow velocities. The exception is low heat flux densities for the minimum pipe diameter under consideration, di = 0.01 m, for which entropy generation minimization implies lower values of flow velocity with respect to the condition εNu−f > 1. Figure 2 and Figure 7 show that, regardless of the adopted criterion for the assessment of the heat exchange process involving ice slurry, the laminar flow range is especially preferred as regards the selection of xs and w. In design practice, the recommended flow velocities in pipes amount to w < 1 ms−1. This velocity range for the flow of ice slurry in pipes with di < 0.02 m enables effective enhancement of the heat transfer process with the minimum entropy generation rates.
The calculations performed suggest that:
During pipe flow, the lowest entropy generation rates were characteristic of small mass fractions of ice in the turbulent flow area and for the flow velocity of 1.5 < w < 2 ms−1.
In the laminar flow area, the lowest entropy generation rates corresponded to the highest analyzed mass fraction of ice xs = 30%.
Regardless of the share of solid particles, the minimum entropy generation rate criterion requires the application of high flow velocities, which for heat flux density values of q ˙ 10 kWm−2 are greater than w > 1 ms−1.

3. Entropy Generation in Heat Exchangers Fed with Ice Slurry

An analysis of the entropy generation rate in heat exchangers was performed under conditions which differ from those present during flow through pipes. For pipe flow, a specific heat flux density, flow velocity and mean mass share of ice are assumed for which the entropy generation rates are calculated. For heat exchangers, the exchanger’s thermal efficiency, ice slurry mass flux and the inlet value of the mass fraction of ice are assumed. A change in ice slurry enthalpy, the heat exchange surface (heat flux density) and the end temperature of the cooled medium result from the balance equations and, as a consequence, make it possible to determine the entropy generation rate.
The heat transfer process enhancement, greater specific enthalpy of the ice slurry and the almost constant melting point enable ice slurry-fed heat exchangers to have a smaller heat exchanger surface than heat exchangers which rely on single-phase refrigerants. The heat exchange surface area has an impact, for example, on the values of flow resistance in the heat exchanger. In this part of the article, the authors discussed the impact of the ice slurry mass flux and the mass fraction of ice on the change in the total entropy generation rate in heat exchangers. To achieve this, the heat transfer process was discussed for two types of heat exchangers: A lamelled air cooler and a plate heat exchanger, which was used as a milk cooler.

3.1. Air Cooler

Assuming that the melting process of the ice slurry occurs at a constant temperature, it is possible, on the basis of [56], to express the entropy generation rate in the ice slurry melting and air cooling process using Equation (17):
Δ S ˙ H = Q ˙ m T ¯ I S + m ˙ a c p a ln ( T a o u t T a i n )
Taking into consideration the ice melting process in the ice slurry, Equation (17) may be converted into (18)
Δ S ˙ H = m ˙   Δ i I S T ¯ I S + m ˙ a c p a ln ( T a o u t T a i n )
If the fluid flow within the heat exchanger is treated as quasi-adiabatic flow, we can state that
T d S ˙ Δ p = V ˙ d p
Whence, for air treated as an ideal gas ( p V ˙ = m ˙ R T ), Equation (19) can be transformed into (20)
d S ˙ Δ p a = m ˙ a R d p p
Equation (20), following integration, leads to Equation (21):
Δ S ˙ Δ p a = m ˙ a R ln ( p a o u t p a i n )
Disregarding the change in volume of the ice slurry in the melting process, you can transform Equation (19) for a liquid refrigerant into (22)
Δ S ˙ Δ p I S = V ˙ ¯ I S Δ p T ¯ I S
Hence, the total entropy generation rate in the ice slurry-fed air cooler may be calculated from Equation (23):
Δ S ˙ = Δ S ˙ H + Δ S ˙ Δ p = m ˙ Δ i I S T ¯ I S + m ˙ a c p a ln ( T a o u t T a i n ) m ˙ a R ln ( p a o u t p a i n ) V ˙ ¯ I S Δ p T ¯ I S
When determining the value of the entropy generation rate for an air cooler in (23), free flow of air out into the surrounding environment (pa-out = 1 bar) was assumed. On the other hand, air pressure at the inlet to the exchanger was pa-out = pa-out + Δpa.
In order to determine the entropy generation rate in the lamelled air cooler, a calculation algorithm was prepared based on Peclet’s equation and balance equations for the air and the ice slurry [15]. The air-side heat transfer coefficient was calculated from Schmidt’s equation [57], while air-side flow resistance calculations were based on Idelcik’s formula [58]. For the ice slurry, flow resistance values were calculated using Equations (7)–(10). Heat transfer coefficients, in turn, were calculated using Equations (11) and (12).
The calculations were made for:
constant thermal efficiency values Q ˙ = 20; 40 kW,
constant air flow stream m ˙ a = 2.2 kgs−1,
constant inlet air temperature Ta-in = 20 °C,
various ice slurry mass flux values 0.3 m ˙ 1.25 kgs−1,
various mass fractions of ice 5 ≤ xs ≤ 30%,
constant number of pipes and feeds in the heat exchanger made of di = 0.01 m pipes.
The result of the calculations were different heat exchange surfaces (different pipe lengths), which implied different entropy generation rates caused by the flow resistances of the refrigerant and air. A change in the entropy generation rate also resulted from variable mass fluxes of the ice slurry and changes in the specific enthalpy of the ice slurry. Figure 8a shows how the mass fraction of ice affects the entropy generation rate in the air cooler in the event of full melting of the ice slurry (variable mass flux of the ice slurry). Regardless of the thermal efficiency, the minimum entropy generation rate criterion prefers the use of ice slurries with the maximum mass fractions of ice. Low mass fractions of ice generate large ice slurry mass fluxes and high flow velocities (higher than 1 ms−1). This results in a significant increase in entropy generation associated with ice slurry and air-side flow resistances. The Bejan number for a 40 kW cooler efficiency changed from 0.04 to 0.43 for the mass fractions of ice of xs = 30% and xs = 5%, respectively. Figure 8b shows how the mass share of ice and the mass flux of the ice slurry affect the total entropy generation rates in the air coolers (non-complete ice melting scenario). Regardless of the mass flux, the minimum entropy generation rates were obtained for the maximum mass share of ice xs = 30%. For a constant mass flux of the ice slurry, the dominant entropy generation component is the entropy generated by the air-side heat transfer the air-side heat transfer process. Lower mass fractions correspond to higher mean temperatures of the ice slurry. Therefore, the receipt of the same heat flux is conditional upon a larger change in the specific enthalpy of the refrigerant. The entropy generation rate as a function of the mass flux is not a monotonic function. This has been illustrated by Figure 9a.
Calculations performed for a constant thermal efficiency of Q ˙ = 20 kW and for the mass fraction of ice of xs = 30% indicate that for the discussed structure of the exchanger, function Δ S ˙ = f ( m ˙ , x s = 30 % ) has its minimum value at the point corresponding to the mass flux of m ˙ = 0.83 kgs−1. At the point of the minimum of function Δ S ˙ , the Bejan number was 0.2, varying in the entire discussed mass flux range between 0.07 ≤ Be ≤ 0.28. The results presented in Figure 9a refer to various surface areas and therefore also to various heat flux densities. The flow velocity curve makes it possible to determine the ice slurry flow velocity in heat exchanger pipes for particular mass fluxes.
Table 4 shows a comparison of the optimum flow velocities of a 30% ice slurry, determined for flow in a di = 0.01 m pipe. Just as in the case of flow through a straight pipe, in the discussed heat exchanger higher heat flux densities also corresponded to higher optimum flow velocities. The comparison (Table 4) indicates that when designing a heat exchanger, it is possible to determine the operational parameters of the cooler which meet the S ˙ ~ S ˙ min ( d , q ˙ ) condition by finding the ice slurry flow velocity using Equation (14).

3.2. Fluid Cooler

If the same assumptions are made as in the air cooler scenario, and if we take fluid cooling into account, the total entropy generation rate in the cooler can be described using Equation (24):
Δ S ˙ = Δ S ˙ H + Δ S ˙ Δ p = m ˙ Δ i I S T ¯ I S + m ˙ f c p f ln ( T f o u t T f i n ) V ˙ ¯ f Δ p f T ¯ f V ˙ ¯ I S Δ p I S T ¯ I S
In order to determine the entropy generation rate in a plate heat exchanger/fluid (milk) cooler, an algorithm [59] was developed, where the heat transfer coefficients and flow resistance values on the side of the cooled fluid were calculated on the basis of the equations proposed by Tarasov [60,61]. For the ice slurry, flow resistance values were calculated using Equations (7)–(10). Heat transfer coefficients, in turn, were calculated using Equations (11) and (12), supplemented with coefficients Bi for a rectangular cross-section, included in Table 2.
The calculations were made for:
constant thermal efficiency values Q ˙ = 20; 30, 40 kW,
a change in the temperature of the fluid cooled in the exchanger between Tf−in = 12, 20 °C and Tf−in = 2, 5 °C (Tf = 10, 15 K),
various ice slurry mass flux values 0.3 m ˙ 12 kgs−1,
various mass fractions of ice 5 ≤ xs ≤ 30%,
two geometrical configurations of the exchanger: np = 6, ns = 1; np = 12, ns = 2.
Figure 10a shows how the mass share of ice and the cooled fluid parameters affect the entropy generation rate in a fluid cooler with the efficiency of Q ˙ = 30 kW, in a scenario involving total melting of the ice slurry (variable ice slurry mass flux). In such a scenario, the dominant component of the entropy generation is the one generated by thermal processes. The Bejan number in the discussed scenario varied between Be = 0.002–0.05. Due to the form of the second component of the sum (24), higher entropy generation rates are characteristic of greater changes in the temperature of the cooled medium and its lower mass flux. Regardless of thermal efficiency and the geometrical configuration of the exchanger, the use of ice slurries with the maximum mass fraction of ice is preferred, as when the complete ice melting condition is assumed, they imply small flow rates of the refrigerant. Figure 10b presents the impact of the mass share of ice and the mass flux of the ice slurry on the total entropy generation rate in a fluid cooler ( Q ˙ = 30 kW, np = 6, ns = 1). In this calculation variant, the constant mass flux condition implies variable values of changes in the ice slurry’s specific enthalpy and various degrees of ice melting. Just as in the case of the air cooler, regardless of the mass flux, minimum entropy generation rates were recorded for the maximum mass share of ice xs = 30%. For a constant mass flux of the ice slurry, the dominant component of entropy generation is that generated by the ice slurry-side heat transfer process. In a plate heat exchanger, the entropy generation due to flow resistances is determined by the flow of the ice slurry. In the calculations performed by the authors, the flow was of a laminar nature and the greatest flow resistances were generated for an ice slurry with the mass share of ice of xs = 30%, for which the dynamic plastic viscosity coefficient and limit shear stress were highest. The highest Bejan number of 0.15 corresponded to xs = 30%, m ˙ = 1.87 kgs−1.
Figure 9b presents the results of calculations of the entropy generation rate for various thermal efficiencies of the exchanger Q ˙ = 20, 30, 40 kW (np = 12, ns = 2) and a constant mass share of ice xs = 30%. Function Δ S ˙ = f ( m ˙ , x s = 30 % ) demonstrates a minimum which shifts towards higher mass flux values as thermal efficiency increases. At the point marking the minimum of function Δ S ˙ , the Bejan number was 0.27, varying in the entire discussed range of mass flux changes between 0.02–0.93. Note that in the calculations performed, the minimum entropy generation rate (due to the value of ReK = 380−530) was recorded for laminar flow, but for w = 1.8−2.3 ms−1. These flow velocity values are significantly higher than the velocity values used in typical plate exchangers (w < 0.5 ms−1). The application of a flow velocity of w < 1 ms−1 leads to an increase in the entropy generation rate by at least 14%–67%, depending on thermal efficiency.
The curves presented in Figure 9b correspond to a constant efficiency of the exchanger. Therefore, in the case of a melting process, entropy generation minimization is not equivalent to the maximization of exchanger efficiency. A sample relationship between exchanger efficiency and the entropy generation rate was presented in Figure 11a. Figure 11 presents the entropy generation rate for a plate heat exchanger fed with a 30% ice slurry or 15% ethanol. For ice slurry, the exchanger efficiency ( η = Q ˙ / ( W ˙ min Δ T max ) ) is determined by the thermal capacity of the cooled fluid ( W ˙ min = m ˙ f c p f ) and hence it has a constant value. For single-phase cooling, exchanger efficiency may be determined by the thermal capacity of the refrigerant or the cooled fluid ( W ˙ min = M I N ( m ˙ e t c p e t ; m ˙ f c p f ) ). It is worth noting that the use of the minimum entropy generation rate as the only design criterion may, for example, generate high investment costs. Figure 11b presents sample calculations where for low mass fluxes ( m ˙ 0.53 kgs−1), lower entropy generation rates are recorded in an ethanol-fed exchanger than in an ice slurry-fed exchanger. If ethanol is used, however, the required surface area of the exchanger is at least 60% higher than in the case of an ice slurry fed exchanger.

4. Conclusions

Using the results of experimental studies and formula for the calculation of flow resistance and heat transfer coefficient values during the melting of ice slurry under forced convection conditions, an analysis of the entropy generation rates was performed during flow in straight pipes and heat exchangers.
The study determined the area of optimal parameters for ethanol ice slurry (w, xs) for different tube diameters and heat flux densities for which the entropy flux generated during the melting process in the conditions of forced convection of ice slurry is minimal. The original aspect presented in the paper is the description of the set of optimal parameters for w, xs, the criteria relationship between the Reynolds number and dimensionless numbers taking into account the effect of geometric and flow parameters on the generated entropy flux.
A new aspect discussed in this paper was the comparison of different strategies for the selection of flow parameters when using ice slurry as a coolant. Three criteria for the selection of flow parameters of ice slurry were compared: Minimum criterion of entropy flux generated by the slurry during flow in a straight channel, minimum criterion of entropy flux generated in a heat exchanger supplied with ice slurry and generalized overall enhancement efficiency criterion. Attention was drawn to the fact that the fulfilment of the criterion εNu−f > 1 more significantly limits the velocity of ice slurry flow from above than the minimum entropy flux criterion.
Unlike the criterion of minimum entropy flux of ice slurry in a straight channel, the criterion of minimum entropy flux in a heat exchanger always prefers the use of maximum mass fractions of ice. In the first case, high ice mass fractions minimize the entropy flux only for low flow velocities (w < 1 ms−1).

Funding

The research was carried out as part of research project 0759/T10/2002/23 funded by the Polish National Science Centre.

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

Aheat transfer surface, m2
BeBejan number, Be = Δ S ˙ Δ p Δ S ˙
Bicoefficients in Equations (11) and (12) (Table 2)
cfFanning friction factor,
Cicoefficients in Equations (15) and (16) (Table 3)
cpspecific heat, Jkg−1K−1
ddiameter, m
G ˙ mass flux density, kgm−2s−1
GzGraetz number, Gz = Re K Pr B d h l
HeHedström number, He = ρ I S τ p d h 2 μ p 2
ispecific enthalpy, Jkg−1
K*consistency index
KFphase change number, K F = r c p I S ( T w T I S )
llength, m
m ˙ mass flux, kgs−1
n*flow behavior index
npplate number
nsnumber of feeds
NuNusselt number, Nu = α d h λ
ppressure, Pa
Pcircumference, m
PCMphase change materials
q ˙ heat flux density, Wm−2
Q ˙ heat flux, W
Rindividual gas constant, Jkg−1K−1
ReISReynolds number for Bingham fluid, Re I S = ρ I S w d h μ p
ReKgeneralized Reynolds number according to Kozicki, Re K = ρ I S w 2 n * d h n * 8 n * 1 K *
sspecific entropy, Jkg−1K−1
S ˙ l entropy generation rate referred to a unit of length, WK−1m−1
S ˙ entropy generation rate, WK−1
Ttemperature, K
T ¯ mean temperature, K
V ˙ volumetric flow rate, m3s−1
wmean velocity, ms−1
W ˙ thermal capacity, WK−1
xsmass fraction of ice, %
X v volume fraction of ice, %
Greek symbols
νspecific volume, m3kg−1
αheat transfer coefficient, Wm−2K−1
εenhancement efficiency
ηexchanger efficiency
λthermal conduction coefficient, Wm−1K−1
μpdynamic plastic viscosity, Pas
ρdensity, kgm−3
τpyield shear stress, Pa
τwshear stress at pipe wall, Pa
Subscripts
aair
aicarrier fluid
calccalculated value
crit, Ccritical value
etethanol
Ffluid
hhydraulic
Hheat transfer
iinternal, ordinal number
ininlet
ISice slurry
llength
Llaminar
obaseline value
outoutlet
refreference value
Sice, ice crystal
Tturbulent
wwall
Δpflow resistance

References

  1. Kitanovski, A.; Poredos, A. Concentration distribution and viscosity of ice-slurry in heterogeneous flow. Int. J. Refrig. 2002, 25, 827–835. [Google Scholar] [CrossRef]
  2. Niezgoda Żelasko, B. Heat transfer of ice slurry flows in tubes. Int. J. Refrig. 2006, 29, 437–450. [Google Scholar] [CrossRef]
  3. Niezgoda-Żelasko, B.; Zalewski, W. Momentum transfer of ice slurries flows in tubes. Experimental investigation. Int. J. Refrig. 2006, 29, 418–428. [Google Scholar] [CrossRef]
  4. Kitanovski, A.; Vuarnoz, A.; Ata-Caesar, D.; Egolf, P.W.; Hansen, T.M.; Doetsch, C. The fluid dynamics of ice slurry. Int. J. Refrig. 2005, 28, 37–55. [Google Scholar] [CrossRef]
  5. Gulipart, J. Experimental study and calculation method of transport characteristics of ice slurries. In First Workshop on Ice Slurries of the International; Institute of Refrigeration: Yverdon-Les-Bains, Switzerland, 1999; pp. 74–82. [Google Scholar]
  6. Christensen, K.G.; Kauffeld, M. Heat transfer measurements with ice slurry. In International Conference of the International Institute of Refrigeration; University of Maryland: Washington, DC, USA, 1997; pp. 127–141. [Google Scholar]
  7. Jensen, E.; Christensen, K.; Hansen, T.; Schneider, P.; Kauffled, M. Pressure drop and heat transfer with ice slurry. In Proceedings of the 2000 International Refrigeration Conference at Purdue; Tree, D.R., Ed.; Purdue University: West Lafayette, IN, USA, 2000; pp. 521–529. [Google Scholar]
  8. Illán, F.; Viedma, A. Experimental study on pressure drop and heat transfer in pipe lines for brine based ice slurry Part II: Dimensional analysis and rheological model. Int. J. Refrig. 2009, 32, 1024–1031. [Google Scholar] [CrossRef]
  9. Mellari, S. Experimental investigations of ice slurry flows in horizontal pipe based on monopropyleneglycol. Int. J. Refrig. 2016, 65, 27–41. [Google Scholar] [CrossRef]
  10. Doetsch, C. Experimentelle Untersuchung und Modellierung des Rheologischen Verhaltens von Ice-Slurries. Ph.D. Thesis, Frauenhofer Institut-UMSICHT, Fachbereich Chemie Technikder Universität, Dortmund, Germany, 2001. [Google Scholar]
  11. Egolf, P.W.; Kitanovski, A.; Ata-Caesar, D.; Stamatiou, E.; Kawaji, M.; Bedecarrats, J.P.; Strub, F. Thermodynamics and heat transfer of ice slurries. Int. J. Refrig. 2005, 28, 51–59. [Google Scholar] [CrossRef]
  12. NiezgodaŻelasko, B.; Żelasko, J. Melting of ice slurry under forced convection conditions in tubes. Exp. Therm. Fluid Sci. 2008, 32, 1597–1608. [Google Scholar] [CrossRef]
  13. Sari, O.; Vuarnoz, D.; Meili, F.; Egolf, P.W. Visualization of ice slurries and ice slurry flows. In Second Workshop on Ice-Slurries of the International; Institute of Refrigeration: Paris, France, 2000; pp. 68–80. [Google Scholar]
  14. Dong, W.I. Experimental study on flow and pressure drop of ice slurry for various pipes. In Fifth Workshop on Ice-Slurries of the International; Institute of Refrigeration: Stockholm, Sweden, 2002. [Google Scholar]
  15. Niezgoda-Żelasko, B. Nowoczesne Systemy Chłodzenia Pośredniego; Wydawnictwo Politechniki Krakowskiej: Kraków, Poland, 2017; pp. 163–176. [Google Scholar]
  16. Doron, P.; Granica, D.; Barnea, D. Slurry flow in horizontal pipes-experimental and modeling. Int. J. Multiph. Flow 1987, 19, 535–547. [Google Scholar] [CrossRef]
  17. Orzechowski, Z. Przepływy Dwufazowe; PWN: Warszawa, Poland, 1990; pp. 61–86. [Google Scholar]
  18. Richardson, J.F.; Zaki, W.N. Sedimentation and fluidisation: PART 1. Trans. Inst. Chem. Eng. 1954, 32, 35–53. [Google Scholar] [CrossRef]
  19. Niezgoda-Żelasko, B.; Żelasko, J. Generalized non-Newtonian flow of iceslurry. Chem. Eng. Process. 2007, 46, 895–904. [Google Scholar] [CrossRef]
  20. Grozdek, M.; Rahmatollah Khodabandeh, R.; Per Lundqvist, P. Experimental investigation of ice slurry flow pressure drop in horizontal tubes. Exp. Therm. Fluid Sci. 2009, 33, 357–370. [Google Scholar] [CrossRef]
  21. Onokokoa, L.; Poirierb, M.; Galanisa, N.; Ponceta, S. Experimental and numerical investigation of isothermal ice slurry flow. Int. J. Therm. Sci. 2018, 126, 82–95. [Google Scholar] [CrossRef]
  22. Skelland, A.H.P. Non-Newtonian Flow and Heat Transfer; John Wiley & Sons Inc.: New York, NY, USA, 1967; pp. 5–469. [Google Scholar]
  23. Alsarraf, J.; Moradikazerouni, A.; Shahsavar, A.; Afrand, M.; Salehipour, H.; Tran, M.D. Hydrothermal analysis of turbulent boehmite alumina nanofluid flow with different nanoparticle shapes in a minichannel heat exchanger using two-phase mixture model. Physica A 2019, 520, 275–288. [Google Scholar] [CrossRef]
  24. Charunkyakorn, P.; Sengupta, S.; Roy, S.K. Forced convective heat transfer in microencapsulated phase change material slurries: Flow in circular ducts. Int. J. Heat Mass Transf. 1991, 34, 819–833. [Google Scholar] [CrossRef]
  25. Niezgoda Żelasko, B.; Żelasko, J. Generalized non-newtonian heat exchange. Flow of ice slurry in pipes. Chem. Process Eng. 2009, 30, 453–473. [Google Scholar]
  26. Mogaji, T.S.; Kanizawa, F.T.; Filho, E.P.B.; Ribatski, G. Experimental study of the effect of twisted-tape insertson flow boiling heat transfer enhancement and pressure drop penalty. Int. J. Refrig. 2013, 36, 504–515. [Google Scholar] [CrossRef]
  27. Gunes, S.; Ozceyhan, V.; Buyukalaca, O. The experimental investigation of heat transfer and pressure drop in a tube with coiled wire inserts placed separetely from tube wall. Appl. Therm. Eng. 2010, 30, 1719–1725. [Google Scholar] [CrossRef]
  28. Akhavan-Behabadi, M.A.; Salimpour, M.R.; Pazouki, V.A. Pressure drop increase of forced convective condensation inside horizontal coiled wire inserted tubes. Int. Commun. Heat Mass Transf. 2008, 35, 1220–1226. [Google Scholar] [CrossRef]
  29. Zimparov, V.D.; Penchev, P.J. Performance evaluation of some tube inserts as heat transfer enhancement techniques. Heat Transf. Eng. 2006, 28, 39–46. [Google Scholar] [CrossRef]
  30. Bejan, A. Method of entropy generation minimization, or modeling and optimization based on combined heat transfer and thermodynamics. Rev. Gén. Therm. 1996, 35, 637–646. [Google Scholar] [CrossRef]
  31. Balaji, C.; Hölling, M.; Herwig, H. Entropy generation minimization in turbulent mixed convection flows. Int. Commun. Heat Mass Transf. 2007, 34, 544–552. [Google Scholar] [CrossRef]
  32. Yilbas, B.S.; Yürüsoy, M.; Pakdemirli, M. Entropy Analysis for non-Newtonian fluid flow in annular pipe: Constant viscosity case. Entropy 2004, 6, 304–315. [Google Scholar] [CrossRef]
  33. Chang, T.B.; Wang, F.J. An analytical inbestigation into the Nusselt number and entropy generation rate of film condensation on horizontal plate. J. Mech. Sci. Technol. 2008, 22, 2134–2141. [Google Scholar] [CrossRef]
  34. Chen, B.C.; Wu, J.R.; Yang, S.A. Entropy generation of forced convection film condensation from downward flowing vapors onto a horizontal tube. J. Chin. Inst. Eng. 2008, 31, 349–354. [Google Scholar] [CrossRef]
  35. Eshafani, J.A.; Modirkhazeni, M. Entropy generation of forced convection film condensation on horizontal elliptical tube. C. R. Mec. 2012, 340, 543–551. [Google Scholar]
  36. Revellin, R.; Lips, S.; Khandekar, S.; Bonjour, J. Local entropy generation for saturated two-phase flow. Energy 2009, 34, 1113–1121. [Google Scholar] [CrossRef]
  37. Ezeora, O.S. Entropy generation analysis and optimum tube length of two-phase flow evaporator tube. In Proceedings of the International Refrigeration and Air Conditioning Conference, West Lafayette, IN, USA, 14–17 July 2008; pp. 2320.1–2320.8. [Google Scholar]
  38. Porras, F.; Guevara, P.; Soriano, G. A study for entropy generation of heat transfer fluid containing Multiwalled Carbon Nanotubes and Microencapsulated Phase Change Materials. In Proceedings of the 11 LACCEI Latain American and Caribbean Conference for Engineering and Technology, Cancun, Mexico, 14–16 August 2013. [Google Scholar]
  39. Bouzid, N.; Saouli, S.; Aiboud-Saouli, S. Entropy generation in ice slurry pipe flow. Int. J. Refrig. 2008, 31, 1453–1457. [Google Scholar] [CrossRef]
  40. Kaymar, A.; Aminossadati, S.M.; Leonardi, C.R. Thermo-hydrodymanics of a helical coil heat exchanger operated with a phase-change ice slurry as a refrigerant. Heat Transf. Eng. 2019, 40, 283–294. [Google Scholar] [CrossRef]
  41. Singh, V.; Aute, V.; Radermacher, R. Usefunles of entropy generation minimization through a heat exchanger modelig tool. In Proceedings of the International Refrigeration and Air Conditioning Conference, West Lafayette, IN, USA, 14–17 July 2008; pp. 2199.1–2199.8. [Google Scholar]
  42. Ahmadi, P.; Hajabdollahi, H.; Dincer, I. Cost and entropy generation minimization of a cross-flow plate heat exchanger using multi-objective genetic algorithm. J. Heat Transf. 2011, 133, 021801. [Google Scholar] [CrossRef]
  43. Qun, C.; Jing, W.; MoRan, W.; Ning, P.; Yuan, G.Z. A comparison of optimization theories for energy conservation in heat exchanger groups. Chin. Sci. Bull. 2011, 56, 449–454. [Google Scholar] [Green Version]
  44. Laskowski, R.; Jaworski, M.; Smyk, A. Entropy generation in a condenser and related correlations. Arch. Thermodyn. 2015, 36, 27–48. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, Z.F.; Jiang, C.; Zhang, Y.F.; Zhou, W.X.; Bai, B.F. Virtual entropy generation (VEG) method in experiment reliability control: Implications for heat exchanger measurement. Appl. Therm. Eng. 2017, 110, 1476–1482. [Google Scholar] [CrossRef]
  46. Strub, F.; Castaing-Lasvignottes, J.; Bedecarrats, J.P.; Peuvrel, C. Application of the second law analysis to a heat exchanger working with ethanol/water ice slurry. Int. J. Thermodyn. 2004, 7, 183–188. [Google Scholar]
  47. Torabi, M.; Zhang, K.; Karimi, N.; Peterson, G.P. Entropy generation in thermal systems with solid structures—A concise review. Int. J. Heat Mass Transf. 2016, 97, 917–931. [Google Scholar] [CrossRef]
  48. Biswa, P.; Basak, T. Entropy generation vs energy efficiency for natural convection based energy flow in enclosures and various applications: A review. Renew. Sustain. Energy Rev. 2017, 80, 1412–1457. [Google Scholar] [CrossRef]
  49. Naterer, G.F.; Camberos, J.A. Entropy and second low fluid flow and heat transfer simulation. J. Thermophys. Heat Transf. 2003, 17, 360–371. [Google Scholar] [CrossRef]
  50. Alnaqi, A.A.; Aghakhani, S.; Pordanjani, A.H.; Bakhtiari, R.; Asadi, A.; Tran, M.D. Effects of magnetic field on the convective heat transfer rate and entropy generation of a nanofluid in an inclined square cavity equipped with a conductor fin: Considering the radiation effect. Int. J. Heat Mass Transf. 2019, 133, 256–267. [Google Scholar] [CrossRef]
  51. Bel, O.; Lallemand, A. Etiude d’un fluide frigoporteur diphasique. 2: Analyse expérimentale du comportement thermique et rhéologique. Int. J. Refrig. 1999, 22, 175–187. (In French) [Google Scholar] [CrossRef]
  52. Melinder, A.; Granryd, E. Using property values of aquaous solutions and ice to estimate ice concentrations and enthalpies of ice slurries. Int. J. Refrig. 2005, 28, 13–19. [Google Scholar] [CrossRef]
  53. Melinder, A. Thermophysical properties of liquid secondary refrigerants. Tables and diagrams for the refrigerants industry. In Tables and Diagrams for the Refrigerants Industry; International Institute of Refrigeration: Paris, France, 1997; pp. 34–36 and 97. [Google Scholar]
  54. Kozicki, W.; Chou, C.H.; Tiu, C. Non-Newtonian flow in ducts of arbitrary cross-sectional shape. Chem. Eng. Sci. 1966, 21, 665–679. [Google Scholar] [CrossRef]
  55. Niezgoda-Żelasko, B. Flow behaviour of ice slurry–experimental investigation. Trans. Inst. Fluid-Flow Mach. 2006, 118, 71–89. [Google Scholar]
  56. Bejan, A. Entropy generation minimization the method and its applications. J. Mech. Eng. 2001, 47, 345–355. [Google Scholar]
  57. VDI-Wärmeatlas. Wärmeübergang an Berippten Oberflächen; VDI Verlag: Düsseldorf, Germany, 1991; pp. Mb1–Mb4. (In Germany) [Google Scholar]
  58. Idielcik, I.E. Handbook of Hydraulic Resistance; CRC Press: London, UK, 1994; pp. 714–756. [Google Scholar]
  59. Niezgoda-Żelasko, B.; Zalewski, W. Chłodnicze i Klimatyzacyjne Wymienniki Ciepła. Obliczenia Cieplne; Wydawnictwo Politechniki Krakowskiej: Kraków, Poland, 2013; pp. 90–97. (In Polish) [Google Scholar]
  60. Tarasov, F.M. O dal’nejšem ulučšenii konstrukcij teploobmennych apparatov piščevoj promyšlennosti. Izvuz Piščevaja Technol. 1970, 1, 127–133. (In Czech) [Google Scholar]
  61. Tarasov, F.M. Teploobmen v plastinčatych apparatch. Izvuz Piščevaja Technol. 1969, 3, 143–147. [Google Scholar]
Figure 1. Outline of the heat transfer process in the ice slurry during its flow through a pipe.
Figure 1. Outline of the heat transfer process in the ice slurry during its flow through a pipe.
Entropy 21 00514 g001
Figure 2. Total entropy generation rate as a function of the mass fraction of ice and ice slurry flow velocity: di = 0.016 m, q ˙ = 8 kWm−2.
Figure 2. Total entropy generation rate as a function of the mass fraction of ice and ice slurry flow velocity: di = 0.016 m, q ˙ = 8 kWm−2.
Entropy 21 00514 g002
Figure 3. Relationship between entropy generation rates S ˙ H , S ˙ Δ p , S ˙ , for di = 0.016 m, q ˙ = 8 kWm−2: (a) xs = 11.7%; (b) xs = 26.6%.
Figure 3. Relationship between entropy generation rates S ˙ H , S ˙ Δ p , S ˙ , for di = 0.016 m, q ˙ = 8 kWm−2: (a) xs = 11.7%; (b) xs = 26.6%.
Entropy 21 00514 g003
Figure 4. Effect of tube diameter on the generated entropy flux: (a) q ˙ = 2 kW, effect of the mass fraction of ice; (b) xs = 11.7%, effect of heat flux density.
Figure 4. Effect of tube diameter on the generated entropy flux: (a) q ˙ = 2 kW, effect of the mass fraction of ice; (b) xs = 11.7%, effect of heat flux density.
Entropy 21 00514 g004
Figure 5. The calculated and experimentally determined parameters w and xs, for which the ice slurry entropy generation rate is at its minimum, di = 0.016 m, 2 kWm−2 q ˙ 8 kWm−2.
Figure 5. The calculated and experimentally determined parameters w and xs, for which the ice slurry entropy generation rate is at its minimum, di = 0.016 m, 2 kWm−2 q ˙ 8 kWm−2.
Entropy 21 00514 g005
Figure 6. A comparison of the values calculated on the basis of correlation (13) with the determined ReIS and q ˙ X v / r , for which S ˙ = S ˙ min .
Figure 6. A comparison of the values calculated on the basis of correlation (13) with the determined ReIS and q ˙ X v / r , for which S ˙ = S ˙ min .
Entropy 21 00514 g006
Figure 7. Curves w ( x s ) S ˙ min , w ( X v ) ε N u f = 1 : (a) di = 0.01 m; (b) di = 0.02 m.
Figure 7. Curves w ( x s ) S ˙ min , w ( X v ) ε N u f = 1 : (a) di = 0.01 m; (b) di = 0.02 m.
Entropy 21 00514 g007
Figure 8. Entropy generation rate in a lamelled air cooler: (a) Δ S ˙ ( Q ˙ , x s ) , m ˙ a = 2.2; 3.2 kgs−1; (b) Δ S ˙ ( m ˙ , x s ) , m ˙ a = 2.2 kgs−1, Q ˙ = 20 kW.
Figure 8. Entropy generation rate in a lamelled air cooler: (a) Δ S ˙ ( Q ˙ , x s ) , m ˙ a = 2.2; 3.2 kgs−1; (b) Δ S ˙ ( m ˙ , x s ) , m ˙ a = 2.2 kgs−1, Q ˙ = 20 kW.
Entropy 21 00514 g008
Figure 9. Minimum entropy generation rate Δ S ˙ min ( m ˙ ) , xs = 30%: (a) lamelled air cooler Q ˙ = 20 kW; (b) ice water cooler Q ˙ = 20–40 kW.
Figure 9. Minimum entropy generation rate Δ S ˙ min ( m ˙ ) , xs = 30%: (a) lamelled air cooler Q ˙ = 20 kW; (b) ice water cooler Q ˙ = 20–40 kW.
Entropy 21 00514 g009
Figure 10. Entropy generation rate in a plate heat exchanger, milk cooler Q ˙ = 30 kW: (a) Δ S ˙ ( m ˙ f , x s ) , ns = 1, 2, np = 6, 12; (b) Δ S ˙ ( m ˙ , x s ) , ns = 1, np = 6.
Figure 10. Entropy generation rate in a plate heat exchanger, milk cooler Q ˙ = 30 kW: (a) Δ S ˙ ( m ˙ f , x s ) , ns = 1, 2, np = 6, 12; (b) Δ S ˙ ( m ˙ , x s ) , ns = 1, np = 6.
Entropy 21 00514 g010
Figure 11. Entropy generation rate in a plate heat exchanger (milk cooler) Q ˙ = 30 kW fed with a 30% ice slurry or 15% ethanol: (a) a comparison of the efficiency η of plate exchangers (b) a comparison of the surface areas of heat exchangers Aet/AIS.
Figure 11. Entropy generation rate in a plate heat exchanger (milk cooler) Q ˙ = 30 kW fed with a 30% ice slurry or 15% ethanol: (a) a comparison of the efficiency η of plate exchangers (b) a comparison of the surface areas of heat exchangers Aet/AIS.
Entropy 21 00514 g011
Table 1. Equations used to calculate the values of physical properties of the ice slurry.
Table 1. Equations used to calculate the values of physical properties of the ice slurry.
PropertiesFormulaComments
Enthalpy of ice slurry [51,52] i I S = ( x s 100 ) i s + ( 1 x s 100 ) i a -
Enthalpy of carrier liquid [51,52] i a = Δ i m i x + c p a ( T 273.15 ) Tref = 273.15 K, ia = ia(xa,T) [51]
Enthalpy of ice [52] i s = r + c p s ( T 273.15 ) r = 332.4 kJ kg−1
Mean specific heat of ice slurry [51] c p I S = ( x s 100 ) c p s + ( 1 x s 100 ) c p a cpa [53], cps [51]
Heat conductivity of ice slurry λIS,w=0 [22] λ IS , w = 0 = λ a [ 2 λ a + λ s 2 ( X v 100 ) ( λ a λ s ) 2 λ a + λ s + ( X v 100 ) ( λ a λ s ) ] λa [53], λs [51]
Ice slurry density ρ I S = ( X V 100 ) ρ s + ( 1 X V 100 ) ρ a ρa [53], ρs [51]
Yield shear stress of ice slurry for xai = 10.6%; ds = 0.1–0.15 mm [3] τ p = 0.013 1.4284 ( x s 100 ) + 73.453 ( x s 100 ) 2 394.64 ( x s 100 ) 3 + 835.82 ( x s 100 ) 4 -
Plastic viscosity of ice slurry for xai = 10.6%; ds = 0.1–0.15 mm [3] μ p = 0.0035 + 0.0644 ( x s 100 ) 0.7394 ( x s 100 ) 2 + 5.6963 ( x s 100 ) 3 19.759 ( x s 100 ) 4 + 26.732 ( x s 100 ) 5 -
Parameter K* K * = ( c μ p ) n * τ w 1 + d n * / c [ c c + d τ w 1 + d / c c d τ w d / c τ p + c 2 d ( c + d ) τ p 1 + d / c ] n * εB= τpw
c,d [54]
Parameter n* n * = c [ 1 ε B 1 + d / c c + d ε B ( 1 ε B d / c ) d ] 1 ε B d [ 1 ε B 1 + d / c c + d ε B ( 1 ε B d / c ) d ]
X v X v = ( 1 + 1 x s x s ρ s ρ a ) 1 -
Table 2. Values of coefficients in Equations (11) and (12).
Table 2. Values of coefficients in Equations (11) and (12).
Cross SectionType of FlowParameterValueApplicability Range
Pipe
0.01 d i 0.02
laminarB12.523% < xs < 30%
wm > 0.1 ms−1
200 < ReK < 2100
B20.11
B3−0.10
B4−0.35
B50.052
turbulentB10.00963% < xs < 30%
wm <4.5 ms−1
2100 < ReK < 11,000
B20.70
B4−0.10
Rectangular
and slit cross-section
0.0055 d h 0.012
laminarB13.665.6% < xs <3 0%
wm > 0.5 ms−1
30 < ReK < 2300
B20.16
B3−0.28
B4−0.12
B50.16
turbulentB10.00323% <xs < 30%
wm < 3.1 ms−1
1900 < ReK < 6000
B20.86
Table 3. Values of coefficients in Equations (15) and (16).
Table 3. Values of coefficients in Equations (15) and (16).
ParameterValue
C1 C 1 = c 11 ( d i d 0 ) 2 + c 12 ( d i d 0 ) + c 13 ; c11 = −2292.0231, c12 = 10,483.068, c13 = −4879.7335
C2 C 2 = c 21 ( d i d 0 ) 2 + c 22 ( d i d 0 ) + c 33 ; c21 = −0.0708, c22 = 0.2925, c23 = 0.6583
C311.3766
C41.15317
C552.0511
q ˙ o 10,000
G ˙ o 380
dio0.01
Table 4. Optimum velocity values according to Equation (14) and for Δ S ˙ min from Equation (23) di = 0.01, xs = 30%.
Table 4. Optimum velocity values according to Equation (14) and for Δ S ˙ min from Equation (23) di = 0.01, xs = 30%.
qwtubewHEBetubeBeHE
80.930.720.080.08
101.191.230.1250.14

Share and Cite

MDPI and ACS Style

Niezgoda-Żelasko, B. Entropy Generation of Forced Convection during Melting of Ice Slurry. Entropy 2019, 21, 514. https://doi.org/10.3390/e21050514

AMA Style

Niezgoda-Żelasko B. Entropy Generation of Forced Convection during Melting of Ice Slurry. Entropy. 2019; 21(5):514. https://doi.org/10.3390/e21050514

Chicago/Turabian Style

Niezgoda-Żelasko, Beata. 2019. "Entropy Generation of Forced Convection during Melting of Ice Slurry" Entropy 21, no. 5: 514. https://doi.org/10.3390/e21050514

APA Style

Niezgoda-Żelasko, B. (2019). Entropy Generation of Forced Convection during Melting of Ice Slurry. Entropy, 21(5), 514. https://doi.org/10.3390/e21050514

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop