Application of the entropy generation minimization method to a solar heat exchanger: A pseudo-optimization design process based on the analysis of the local entropy generation maps
Highlights
► An entropy generation minimization method is applied to a solar heat exchanger. ► The approach is heuristic and leads to a pseudo-optimization process with CFD as main tool. ► The process is based on the evaluation of the local entropy generation maps. ► The geometry with pin-fins in general outperforms all other configurations. ► The entropy maps and temperature contours can be used to determine the optimal pin array design parameters.
Introduction
During the last three decades the EGM (Entropy Generation Minimization) method has become a well-established procedure in thermal science and engineering. Its theoretical foundations are to be found in the well-known Guoy–Stodola theorem which states that for whatever open process whose evolution can be approximated as a succession of quasi-equilibrium states and during which the system is in thermal contact with an ambient at T0, the lost available power , i.e., the difference between the ideally produced power and the one really extracted, is proportional to the global rate of entropy generation [1]. The lost power is always positive, regardless of whether the system is a power source (e.g., an expander) or a power sink (e.g., a compressor). Hence, a reduction in the entropy generation rate leads directly to an increase in the efficiency of the system. The EGM method relies on a simultaneous heat transfer and engineering thermodynamics approach that leads to the construction of realistic models for the systems object of study, be they heat transfer devices, power components or complex plants. Once a model has been established, the strong link between the global entropy generation rate and the efficiency is exploited by the designer who can compare different configurations (“design options”) that either produce the same output with less irreversible losses or use the same amount of resource input to generate a larger output; both cases corresponding of course to a higher resource-to-end use efficiency. Naturally, the minimization of the entropy generation is not an easy task in real design applications; the level of complexity of the model clearly increases when complicated boundary conditions apply and/or when the operating point is varying in time. Nevertheless, several applications of the method can be found in the recent technical literature. Such applications may be divided in two categories, those based on deterministic approach and those based on a heuristic approach. Following the work done in Refs. [1], [2], the overwhelming majority of applications for heat transfer problems belong to the first category. This technique is applied to “classic” problems in Refs. [1], [2], [3], [4], [5]; pin-fins geometries are optimized with this method in Refs. [6], [7], while plate-fins heat sinks are optimized in Refs. [8], [9]. The fundamental characteristic of the deterministic approach is the analytical definition of the global rate of entropy generation ( [W/K]) as a function of critical design parameters, like geometry, dimensions and working conditions. is studied numerically and the optimal geometry is defined by the set of parameters which minimize the global rate of entropy generation. Since the function is determined using the correlations for average heat transfer rates and fluid friction available in literature, the deterministic approach ultimately results in a lumped-sum parameter model, that is the thermodynamic behaviour of the entire system (or control volume) is described globally by a single mathematical function. While this analytical approach has clear advantages, it is not applicable to complex problems where no reliable and explicit correlations for the mean heat transfer and fluid friction are available. Typical examples are turbomachinery and (convective) heat exchangers design problems. For such kind of problems, the thermodynamic behaviour of the system is more accurately described by a local, distributed-parameter model, and the minimization problem is best solved with a heuristic approach: the initial configuration is successively improved by introducing design changes based on a critical analysis of the local entropy generation maps obtained by means of CFD (Computational Fluid Dynamics) simulations. While this approach has been already adopted in some turbomachinery problems [10], [11], examples for heat exchangers like the one covered in this work appear quite rarely in the archival literature: one example can be found in Ref. [12]. The heuristic approach consists in focussing the attention first on the local entropy generation rates , and in considering the global one only after having carefully studied the implications of the local irreversibility on the overall design. The calculation of the local entropy generation rates requires that the velocity and temperature fields are known; therefore a CFD solver, i.e. a distributed-parameter model, needs to be employed. Once a solution is obtained, through the visual analysis of the local entropy maps the practitioner is able to identify the areas where the entropy is produced at the higher rates, that is the areas which demand design modifications. Even if the heuristic approach strongly emphasises the role of the local entropy rates, it must be noted that is nevertheless the quantity that reflects the thermodynamic performance of the whole system, that is quantity that allows the designer to discern the better performer from a second-law perspective. The aim of the present work is to show an application of the heuristic approach to an industrial heat transfer problem.
This work is organized as follows: more details about the adopted heuristic approach are given in the second paragraph; the third paragraph describes the heat transfer problem the method has been applied to; the fourth paragraph deals with the numerical simulations; the fifth paragraph presents the results and in the last paragraph conclusions are drawn.
Section snippets
Strategy for probing the solution space and description of the heuristic procedure
The entropy generation rate can be shown, for the case in study (i.e. in the absence of phase changes and chemical reactions), to consist of two parts [1]: one, called “viscous” , that depends on the physical viscosity, on the local temperature of the fluid and on the second power of the local velocity gradient, and another, called “thermal” , that depends on the physical conductivity, on the square of the local temperature of the fluid and on the second power of the local temperature
The TAK plant
The object of this study is the solar roof tile, called TAK, which is part of a complex system (the TAK plant) to provide the heating and cooling of a house (Fig. 1). The TAK plant can be divided in four main sub-systems:
- •
TAK sub-system: consisting of south-facing roof-tiles through which a secondary fluid, called XEN (a detailed description of which is given in the next paragraph) flows by forced circulation.
- •
Heat-pump sub-system: a first heat exchanger (EX1) allows connecting the TAK sub-system
Alternative geometries
As stated before, the (“family” of) geometries chosen for further development are geometry A and E. In the next paragraphs the variations on each of these two families are presented. In Section 6, the results and the rationale behind each design change is discussed.
Meshing
Four starting geometries (A.0, B, C, D in Fig. 4(a)) were made available by the GreenMind Company which provided the relative IGES (Initial Graphics Exchange Specification) files. Once imported into GAMBIT, the mesh generator, all of the geometric features of the TAK deemed inessential for the goals of the present simulations were removed: only the fluid channel, the inlet and outlet ducts and the aluminium slab are retained. The obvious symmetry of the design, and with the exception of
Geometry E
The central idea is that we expect low values of the local entropy rates in areas where similar boundary layers merge. The word “similar” means that the two boundary layers must draw origin from the same boundary conditions: that is the same surface geometry and the same ΔT. On the contrary, we expect high values of and inside the boundary layers, since these are zones where velocities and temperatures change rapidly; so, from a thermodynamic perspective, we would like to reduce them.
Conclusions and suggestions
A careful and detailed description of the steps of the heuristic optimization of a particular type of heat exchanger has been presented. The thermodynamic fields, in all simulations, have been evaluated on a satisfactorily refined grid, using the thermal or the total entropy generation rate as an objective function. In this way, a sufficiently large and reliable database of numerical data has been obtained. The pseudo-optimization process described in this work is an effective tool in the hands
References (19)
Entropy generation through heat and fluid flow
(1982)A study of entropy generation in fundamental convective heat transfer
Journal of Heat Transfer
(1979)- Benedetti P, Sciubba E. Numerical calculation of the local rate of entropy generation around a heated finned-tube. The...
- et al.
Minimization of the local rates of entropy production in the design of air-cooled gas turbine blades
Journal of Engineering for Gas Turbines and Power
(1999) - et al.
A CFD-aided design procedure for compact heat exchangers
Numerical Heat Transfer
(2005) - et al.
Second-law analysis on a pin-fin array under cross flow
International Journal of Heat and Mass Transfer
(1997) - et al.
Optimization of pin-fin heat sinks using entropy generation minimization
Components and Packaging Technologies, IEEE Transactions on
(2005) - et al.
Optimization of plate fin heat sinks using entropy generation minimization
Components and Packaging Technologies, IEEE Transactions on
(2001) - et al.
Optimal design methodology of plate-fin heat sinks for electronic cooling using entropy generation strategy
Components and Packaging Technologies, IEEE Transactions on
(2004)
Cited by (56)
Performance evaluation criteria based on the first and second laws of thermodynamics
2024, International Communications in Heat and Mass TransferEnergy and entropy generation analysis in a supercapacitor for different operating conditions
2022, EnergyCitation Excerpt :Ramírez-Minguela et al. [26] performed an entropy generation analysis in the solid oxide fuel cell to evaluate the effects of four different biogas compositions in the thermodynamic irreversibilities. Giangaspero et al. proposed an entropy generation minimization method to optimize the forced convective cooling system of LED-based spotlight [27] and to optimize the solar heat exchanger configurations [28]. However, although some entropy generation studies are implemented in the literature [29–33], the entropy generation analysis related to the supercapacitor cell accounting for the thermodynamic irreversibilities, to the best of our knowledge, has never been reported yet.