## Abstract

The efficiency and ability to control the energy exchanges in thermal energy storage systems using the sensible and latent heat thermodynamic processes depends on the best configuration in the heat exchanger’s design. In 1996, Adrian Bejan introduced the Constructal Theory, which design tools have since been explored to predict the evolution of the architecture in flow systems. This chapter reviews the fundamental knowledge developed by the application of the constructal principle to the energy flows in the design of heat exchangers of thermal energy storage systems. It introduces the Svelteness and scale analysis, as two constructal tools in the evolutionary design of engineering flow systems. It also includes the analysis on essential scales of several configurations, or energy flow architectures, toward establishing the main guidelines in the design of heat exchangers for storing thermal energy.

### Keywords

- thermal energy storage
- heat exchangers
- constructal theory
- phase-change materials
- flow architecture

## 1. Introduction

Engineering systems capture a fraction of the total amount of thermal energy available from renewable sources, and to increase the energy system reliability, the research and development of the flow architecture in thermal energy storage systems is of paramount importance.

One of the crucial issues is the characteristic fluctuations in the availability of energy from renewable resources and wasted energy in industrial processes. The design of efficient thermal energy storage systems is an essential step toward meeting the consumption demands of electricity and heat [1]. Therefore, there is growing attention in the development of thermal energy storage systems to produce adequate energy savings and utilization, with a relevant impact on numerous and diverse applications [2, 3].

There are two basic approaches to thermal energy storage. One using the sensible heat without phase-change (SHS - Sensible Heat Storage), and another using the sensible heat * and*phase-change (LHS - Latent Heat Storage), as depicted in Figure 1. The thermal balance describing each approach is given by

where

In both approaches, the upper limit for the final temperature is the saturation value associated with the vaporization of the liquid. Because of its high heat capacity, water is the most used fluid for SHS. However, theoretically, Huang et al. [4] showed the energy stored in a water-based system is one order of magnitude lower than the energy stored in a PCM. And, experimentally, Kaygusuz [5] showed evidence of a PCM system able to store up to 60% of its theoretical maximum, which represents almost the double value of the theoretical storage capacity in water-based systems. These results were the motivation for further investment in the development of LHS systems. Compared to single-phase heat storage systems, LHS systems store the same amount of energy using more compact systems, reducing production and maintenance costs.

In LHS systems, Figure 1 on the right represents the three theoretical stages of the storage process. Initially, the Thermal Storage Material (TSM) is in its solid-state, and it stores (charges) energy through:

Step 1) the sensible heat until the melting temperature;

Step 2) the latent heat component until a complete phase-change of the TSM from solid to liquid. In this stage, the time of the melting process depends on the advancement of the solid–liquid interface (melting front), according to the configuration of the thermal fluid circuit of the heat exchanger (HE) immersed in the TSM;

Step 3) after the total liquefaction of the TSM, energy storage continues until the liquid reaches the saturation temperature of the next phase-change without compromising the volume of the Thermal Storage facility.

The storage of energy in these steps depends on the thermal properties of TSM, and the thermodynamics of the storage process. Still, the main challenge is the design of heat exchangers, as the engineering system that enables the flow of energy from the sources (renewable and non-renewable) to the TSM, disregarded in recent comprehensive reviews on thermal energy storage [6, 7]. Namely, this design has a significant impact on the charging and discharging times, if using renewable energy sources, given their limited time-window throughout the day.

The standard approach in the design of heat exchangers is to optimize the thermal and hydrodynamic energy flows. It uses an iterative process based on previous work, and typical working conditions, such as the amount of fouling and pressure drop in the system, testing a significant number of trial-and-error designs until the values for the heat transfer performance, hydrodynamic effects and longevity are within pre-established requirements. A common trait in this standard approach is the lack of evaluation criteria grounded on the underlying physical processes. Constructal design distinguishes from the standard approach in providing the evaluation criteria in such a way. Therefore, the purpose of this chapter is to synthesize and present an evolutionary design approach (not optimization) using tools based on constructal theory. Therefore, the novelty is to include the architecture of thermal energy storage systems at the design stage [8] and investigate the best way to introduce the freedom to morph to overcome the shortcomings on charging and discharging periods due to prescribed, rigid and fixed designs when subjected to daily and seasonal changes.

## 2. Design tools in constructal theory

In 1996, Adrian Bejan [9], professor at Duke University, proposed a Constructal Theory to explain the evolution of configurations in nature stating,

In practice, when using the constructal theory in engineering, one finds the best direction for the flow structures emerging from of what facilitates movement, designated as * constructal design*. And this design establishes a relation between what flows (energy, fluids, people, etc.) and the geometry of the flow architecture, in such a way that it becomes a global property of the engineering system. This property, the

*, is a tool in constructal design.*Svelteness

In thermal energy storage, there are several length and time scales competing in the unfolding heat transfer processes, characterized by mass, momentum and energy balances of the system. However, not all the terms correspond to the dominant scales setting the overall result of charging and discharging of energy. A * scale analysis*is the second tool in constructal design explored in this section, allowing a proper definition of the relevant scales, and their implication to the heat exchanger design.

### 2.1 Svelteness of flow configuration

The * Svelteness*(Sv) is this global geometric property of the flow space, which guides the engineering practice in the assessment of the flow design performance. This property corresponds to the relation between two length scales of the flow system configuration: an external (

According to Bejan and Lorente [10], the evolutionary direction is that of * vascularization*, implying an increase of the Svelteness. Therefore, as an example applied to TES, what is the best constructal design solution for a simple tube inside a tank? Should the design be a PCM in the inner tube and having the thermal fluid flowing through the tank (see [11]), or the opposite (see [12])?

In this example, the inner tube and tank are cylindrical with a length of corresponded to the case where the energy flows inward to the PCM on the inner tube, and Sv

^{i}

*when it flows outward to the PCM in the tank, the relation between Sveltenesses would be*

^{o}

If a larger Svelteness points in the evolutionary design of the outward energy flow from the thermal fluid inside the inner tube, Sv/

^{o}

Considering * d*, and the thermal fluid circulating in the tank, the storage tank volume should equal the total volume of all the tubes with the “necessary” thermal fluid volume,

Considering the previous constructal analysis using the Svelteness, storing energy with PCM material inside the tubes is only worthy when

On the other hand, recent works as that of Agyenim et al. [12], point toward having the PCM in the tank, instead of inside the inner tubes, and in these cases with a single tube through which circulates the thermal fluid,

### 2.2 Method of scale analysis

* Scale analysis*or

*is a problem solving method useful to obtain essencial and expedite information of several energetic processes [13]. It is not the same as a the dimensional analysis performed in fluid mechanics, but to assess the importance of the order of magnitude of the parameters involved in heat transfer processes, and extract the relevant scales from their governing equations. For more details on the principles of scale analysis, see Bejan [13] (pp. 17–20). Here, one uses an example to illustrate the method.*scaling

Consider the example above of an LHS system with the PCM inside a tube and the thermal fluid circulating around it. If there is a sudden change in the thermal fluid temperature (

Scaling means using the symbol

where

which solved for the time scale results in

Bejan [13] contains the synthesis for all the rules in a scale analysis. However, the example above is enough to explain the procedure applied later in section 3. The following section exemplifies the application of constructal theory as an evolutionary design method to develop heat exchangers in sensible and latent heat storage engineering systems.

## 3. Constructal theory in thermal energy storage heat exchangers

One of the essential elements in a constructal theory analysis is the freedom to morph of flowing configurations. Therefore, once we identify what is the flow under analysis, one can better understand what its freedom to morph means. On the other hand, the heat exchanger in thermal energy storage corresponds to the structure obtained after morphing through which energy flows from a source, usually the thermal fluid, to the storage material (e.g. a solid or a phase-change material, PCM). Depending on the storage material, the heat transfer mechanisms vary, and, accordingly, the energy storage scales. For example, if the material is solid, the heat transfer mechanism is diffusion and the mode is the one on the left of Figure 1. But if one uses a PCM, the energy storage story follows the second mode on the right, involving natural convection, a melting process, a solid–liquid interface moving boundary, and all these elements lead to additional complexity of the heat exchanger, affecting the energy storage scales.

### 3.1 Heat exchangers in sensible heat storage

Consider an underground volume of solid with a network of channels through which a thermal fluid transports energy for storage purposes. The storage mechanism is heat diffusion from the channels outer area to the volume of solid. What should be the structure of the channels network? Combelles et al. [14] explored this TES system with tree-shaped configurations of 2D channels made of parallel plates the length

The first step is to characterize the architecture of each configuration type (2D or 3D) in terms of their Svelteness, as the global property of the system, which relates the external length scale given by the total length of the flow network,

which developed to depend on the length and diameter of the first branch (

The Svelteness for both general configurations, considering the relations in Eqs. (13) and (14), depends on two major features of the configuration: its complexity (

and

If the Svelteness indicates the evolution of the flow configuration, one should connect the flow architecture complexity,

The amount of energy stored depends on the material, but in this TES system, the relevant scale is the storage time and the evolution of the temperature in the conductive solid. Considering the solid is, initially, at

This analysis focuses on the timescales of energy storage and the corresponding effect of the configuration complexity (number of bifurcations,

• the timescale of fluid traveling the

since

• and the timescale of thermal diffusion accross the channel;

Combelles et al. [14] argue that in the case where thermal diffusion in the channel’s boundary layer is a slower process than the fluid traveling through the channels,

The third timescale, and the longer, corresponds to the time it takes to store energy in the solid volume,

From the numerical simulations, the evolution of the dimensionless thermal potential

where

with

with as considered in their simulations.

_{k}

Figure 3 shows the results for the complexity degree scale (

In absolute terms, the constructal design of heat exchangers in diffusive TES systems suggests the choice of a 3D configuration, rather than a 2D, for a faster energy storage, since it leads to

### 3.2 Heat exchangers in latent heat storage

Energy storage systems using the latent heat of a certain phase-change material (PCM) rely on the heat transfer mechanisms of diffusion and natural convection. Initially, the PCM is in its solid state and it stores heat by diffusion close to the channel containing the thermal fluid, or fin, until it reaches the fusion (or melting) temperature (

The invasion period corresponds to the time interval until the solid–liquid interface reaches a distance equivalent to the process characteristic length. The consolidation period corresponds to the remaining time until all the PCM in the LHS system is in its liquid state. These periods do not, necessarily, correspond to the timescales associated to the diffusive and convective heat transfer processes. The charging and discharging times of LHS depend on the heat exchanger design and the dominant heat transfer mechanisms through which energy flows from its source to the PCM.

There are several design configurations investigated with a constructal approach for the heat exchangers using phase-change to store energy in PCM. Figure 4 presents three configurations reported in the literature. The configurations with a vertical pipe [15] is the less prone to morphing. The helical pipe [16] in a cylindrical PCM enclosure is fixed, but the ability to vary the number of turns and the diameter of each turn increases the system’s freedom to morph. Finally, an advancing heat source line invading the PCM material aims at the theoretical design with the greatest freedom to morph [17]. One of the novelties in constructal design of engineering systems is determining the Svelteness as expression of its architecture, and the system’s freedom to morph, considering an evolutionary path toward vascularization, i.e. an increase of its Svelteness. The constructal analysis of all designs also implies the investigation of length and timescales associated with heat transfer mechanisms, and the possible effect of the Svelteness in these scales.

Lorente et al. [15] performed a scale analysis to analyze the latent thermal energy storage where energy flows from the thermal fluid circulating inside a central vertical pipe and the surrounding PCM. The dominant heat transfer mechanism is natural convection. The Svelteness in this case would have the height of the enclosure (

The maximum energy one can store in this first TES configuration is

The relation between the maximum amount of energy stored and this time-scale, including Eq. (23) in (24) leads to

Considering the case of Ogoh and Groulx [19], Figure 5 depicts the adding of disc-shape fins to the original vertical pipe in cylindrical PCM enclosure configuration, where each disc-shape fin corresponds to a * construct*(represented on the right).

The total height of the cylinder (

Considering the Svelteness (Sv) for this construct as the ratio between the external characteristic length based on the upper and bottom areas of the disc-shape fin,

Ogoh and Groulx [19] argue for a neglecting effect of convection between annular fins, thus, the energy stored by phase-changing the PCM from solid to its liquid state occur by conduction. In this sense, the most important scale characterizing the melting front (

Assuming a linear temperature distribution across this layer where the heat transfer occurs by diffusion, one could quantify the heat flux as

with

Therefore, since Sv

Figure 6 shows the results for this limit considering the properties reported in [19].

The constructal design analysis of this heat exchanger indicates diminishing returns of less than 10% for a number of fins above

The helical coil illustrated in Figure 4 is an alternative to the vertical pipe, and a geometry where the freedom to morph is larger because of the ability to change the helix diameter, number of turns and pitch angle. Alailami et al. [16] explored the morphing ability of the system in its design stage to optimize the storage of energy analyzing two scales. The timescale of heat penetrating the storage material from the boundaries of the helical coil to the cylinder diameter,

Afterward, focusing the analysis on

Without using the work of Alailami et al. [16], Joseph et al. [21] performed an experiment of this configuration to store energy in a PCM. The authors used unoptimized values for the helical coil (

The last geometry in Figure 4 is theoretical and corresponds to the greatest freedom to morph through an advancing heated line that can bifurcate at some point. The analysis performed by Bejan et al. [17] focus on the invasion (line advances until the storage boundaries) and consolidation (all PCM melts) stages, and introduces a tree invasion pattern with a complexity level up to

A final comment concerns the possible contribution of constructal design to optimize the total cost of heat exchanger design methods, which is not a direct correlation. As shown by Azad and Amidpour [22], since the constructal design provides an evolutionary perspective on the optimum geometric features of heat exchangers, it can lead to a substantial reduction of the total cost, compared to more standard design approaches. Namely, in the aforementioned work, the authors used constructal theory to optimize shell and tube heat exchangers, and the new approach allowed to reduce this cost by 50%. However, the application of a similar reasoning in the development of heat exchangers for thermal energy storage is in need of more research.

## 4. Conclusions

Thermal energy storage is one of the preeminent options to face the energy challenges of this century, providing a high energy saving potential and effective utilization. However, in these systems, the architecture of the heat exchangers through which energy flows, during charge and discharge, is of paramount importance. While most approaches optimize heat exchanger designs, the one presented in this chapter, based on constructal theory, follows an evolutionary design, meaning that the configuration explored at the design stage is dynamic and free to morph. It is not pre-defined, rigid, or still, but considers how it should evolve toward the greater access of the energy currents that flow through it.

Thermal energy storage systems follow two thermodynamic processes using the sensible heat of the energy storage material, or, besides the sensible heat, also the latent heat, as in Phase-Change Material (PCM). After introducing the general considerations on these systems, this chapter presents two design tools in constructal theory: the Svelteness, as a global property of any flow system, which tends to increase and evolve toward vascularization; and the scale analysis, as an expedite problem solving tool that allows obtaining relevant information of the several energetic processes involved.

Using the design tools presented, this chapter reviews and further explores the constructal theory approach in the development of heat exchangers for sensible and latent thermal energy storage configurations. The analysis evidences the explanatory potential of the constructal approach, increasing the sensibility of the engineer to the advantages of including the freedom to morph at the design stage of heat exchangers in thermal energy storage.

## Acknowledgments

The author would like to acknowledge project UIDB/50022/2020 and UIDP/50022/2020 of ADAI for the financial support for this publication.

## Nomenclatures and Abbreviations

Specific heat [J^{−1} K^{−1}]

Diameter [m]

Height [m]

Latent heat of fusion [J/kg]

Thermal conductivity [W^{−1} K^{−1}]

Length [m]

Mass [kg]

Mass flow rate [kg/s]

Number of tubes or Complexity degree [−]

Energy stored [J]

radial coordinate [m]

Rayleigh number [−]

Svelteness [−]

Temperature [K]

time [s]

Final temperature [K]

Initial temperature [K]

Melting temperature [K]

Volume [m

Greek Symbols

Thermal diffusivity [m^{2}/s]

Length [m]

Temperature difference [K]

Scale factor [−]

Scale factor [−]

Normalized temperature difference [−]

Density [kg/m^{3}]

Timescale [s]

Void fraction

Subscripts

average

cylinder

external

fluid, fin

helical

internal, inner

melting

maximum

outward

solid

stored

tube

thermal fluid

Two-dimensional

Three-dimensional

Heat Exchanger

Latent Heat Storage

Phase-Change Material

Sensible Heat Storage

Thermal Energy Storage

Thermal Storage Material