Fluid Dynamics & Materials Processing

DOI: 10.32604/fdmp.2022.021839

ARTICLE

Thermal Analysis of Melting Occurring Inside a Finned Rectangular Enclosure Equipped with Discrete Pulsed Protruding Heat Sources

Brahim Amahan1, Hamid El Qarnia2,* and Ali El Afif1

1Faculty of Sciences, Department of Physics, Laboratory of Innovation in Sciences Technology and Modeling, Chouaib Doukkali University, El Jadida, Morocco
2Faculty of Sciences Semlalia, Fluid Mechanics and Energetic Laboratory, Department of Physics, Cadi Ayyad University, Marrakesh, Morocco
*Corresponding Author: Hamid El Qarnia. Email: elqarnia@uca.ac.ma
Received: 08 February 2022; Accepted: 01 March 2022

Nomenclature

1  Introduction

The melting of a solid-liquid phase change material (PCM) in a rectangular enclosure equipped with discrete protruding heat sources, flush-mounted on one of its vertical walls, has been the subject of several analytical, experimental and numerical investigations [1–10]. Chu et al. [1] were the first to investigate the effect of heater size, location, aspect ratio and boundary conditions for a laminar flow and arrive at a relationship between the Rayleigh number and the optimum position of the heat source that maximizes heat transfer. Zhang et al. [2] have studied the melting of n-octadecane employing a constant and uniform heat generation and found out a reduction of 50% in the heat source average temperature compared to that obtained by air-cooled natural convection (for a certain period of time). In the case a cavity with thermally insulated horizontal walls and three sources, Ju et al. [6] showed that this reduction can reach up to 70% using a PCM instead of ethylene glycol. Binet et al. [3] showed that a high aspect ratio better controls the heat sources temperature and offers relatively extended melting durations. The influence of changes in the pulses frequency, heat source location and the aspect ratio on the thermal performance of the PCM unit appear to be significant as examined by Krishnan et al. [7]. In order to mitigate the hot spots, Birinci et al. [8] have determined an optimum ratio using six different heat sources with different electrical powers, where the average and total heat generation rates are kept equal for all the examined cases. The Reynolds number varied from 792 to 3962 and all the three known heat transfer mechanisms have been taken into account with an integrated approach to calculate the dimensionless global conductance of the integrated circuit pack. Three numerical simulations have been performed using ANSYS Fluent and measurements have been collected for the surface temperature. The output parameters of the study are the surface and hot spot temperatures, Nusselt number and dimensionless global conductance change with Reynolds number and heat generation ratio. The results have revealed that the relative decrease in the heat generation rates are convenient for the hot spot mitigation.

The inclusion of fins, as an artificial additional heat exchange surface, within a rectangular enclosure filled with PCM and equipped with discrete protruding heat sources has been considered in the past, however in a few limited studies, where the investigation shave focused on of the effect of the fin dimensions, positions and numbers. Joneidi et al. [9] carried out an experimental study to assess the melting phenomenon occurring in a horizontal heat sink with plate fins by changing their height and number. With regard to their results, increasing the number of fins increases the base plate temperature, decreases the melting rate and leads to a much uniform temperature distribution, which demonstrates the correlation between these quantities and the cavity geometrical parameters. On the other hand, an increase in the fin thickness results in a reduction in the melting time [10]. Abdi et al. [11] studied the effect of vertically oriented fins on the heat transfer and energy density in an enclosure heated from the bottom. Their numerical results showed that the number and length of the fins have a significant effect on the natural convection patterns and hence on the heat transfer rates as well as on the energy storage. Using a large number of longer fins leads to the highest melting rate, while small number of shorter fins results in the lowest heat transfer rate. In addition to that, it appears that elongating the fins is significantly important than increasing their number. Joshi et al. [12] have numerically analyzed the effect of the size and location of the fin on the melting time of PCM. They have examined four different configurations and showed that the fin-to-PCM volume ratio decreases by half while the total volume of the thermal energy storage increases.

Our aim in this study is to examine the effect of a horizontal rectangular fin position on the melting process of a phase change material inside a rectangular enclosure heated by multiple pulsed protruding discrete heat sources. As mentioned earlier, the inclusion of fins inside cavities has not received enough attention even though it is of interest to engineering and applied applications. In the next section, we formulate a two-dimensional model capable of describing heat transfer features of the proposed finned rectangular enclosure during the melting process of a PCM. The melting process is initiated by heat extracted from three pulsed protruding heat sources mounted on the left wall of the enclosure. In order to describe the fluid flow with natural convection, mass, momentum and energy conservation equations are derived, appropriately normalized, numerically solved using the finite volume method and finally validated against experimental data selected from the literature.

2  Model Formulation

The system under consideration, shown in Fig. 1, consists of a phase change material (PCM), n-eicosane, filling a rectangular enclosure of height l and width w attached to a substrate of thickness e. Three identical pulsed protruding heat sources (simulating electronic components), each of height lc and thickness ec, are flush-mounted on the left wall. The distance between two consecutive heat sources is Γ, and the distance between the bottom enclosure wall and the lower heat source is δ. The length and thickness of the horizontal fin are H and ef, respectively. The thermal conductivity, kc, the specific heat capacity at constant pressure, cp,c, and mass density, ρc of the heat sources are different from those of the wall (substrate), ks, cp,s and ρs. The enclosure can be used either as a latent heat storage unit to store thermal energyorasa heat sink to cool electronic devices such as the heat sources, the conductive walls and the substrate. The pulsed power generated within the heat sources per unit length, denoted by Qv, is displayed in Fig. 2 and thus Q¯  = Qv, min t1 + Qv, max t2, where t1 and t2 are the durations of minimum (Qv, min) and maximum (Qv, max) pulsed power generated by the heat sources set in this study to 10 and 50 W/m.

Figure 1: The schematic view of the rectangular enclosure

Figure 2: Pulsed power generated by the heat sources during the first two cycles

In all solid media (solid PCM, left wall, heat sources and fin), the heat transfer occurs by pure conduction, while in the liquid PCM it is ensured by advection. We track the moving solid-liquid interface that separates the solid and liquid phases of the PCM by solving the equation for the liquid volume fractionf1. The heat transfer and melting process are governed by mass, momentum and energy conservation equations where a number of assumptions have been used:

•   The liquid PCM is Newtonian and incompressible;

•   The flow is laminar;

•   The Boussinesq approximation is adopted (assuming a linear variation of density with temperature);

•   The thermophysical properties are held constant over the range of temperatures considered;

•   The volume change of the PCM during melting is neglected;

•   The phase change is isothermal;

•   The viscous dissipation is neglected in the liquid PCM.

Under these assumptions, the independent state variables necessary for an adequate description of such a physical problem are the x and y components of the transient velocity field u(x, y, t) and v(x, y, t) and the field of temperature T(x, y, t). We have appropriately scaled the governing equations using the following dimension less space and time, physical parameters and the independent state variables:

X=xl0, Y=yl0, τ=αm,ll02t ,U = u l0αm,l,  V = v l0αm,l, θ=T−TmΔT,P=p lo2ραm,l2 (1)

where l0=lw−3lcec−He is the constant characteristic length, αm,l the thermal diffusivity, Tm the melt temperature, ρ the overall mass density and p is pressure. The mass of PCM remains unchanged in the present study. The normalized governing equations can be written under the following form:

Mass conservation equation

∂U∂X+∂V∂Y=0     (2)

Momentum equations

∂U∂τ+U∂(U)∂X+V∂(U)∂Y=Pr(∂2U∂X2+∂2U∂Y2)−∂P∂X+SU  (3)

∂V∂τ+U∂(V)∂X+V∂(V)∂Y=Pr(∂2V∂X2+∂2V∂Y2)−∂P∂Y+SV (4)

Energy conservation equation

∂θ∂τ+U∂(θ)∂X+V∂(θ)∂Y=α¯(∂2θ∂X2+∂2θ∂Y2)+Sθ   (5)

where the different source terms involved in the model are given by

SU=−C (1−fl)2fl3+bU,SV=−C (1−fl)2fl3+bV+RaPrθ,  Sθ=δ1((δ2−1)1Ste∂fl∂τ+δ23EcLc) (6)

Written in terms of the Heaviside function

δ1= {1   heat source and PCM0          left wall

δ2={1heatsource0PCM (7)

Two constants appear in Eq. (6), C=1025kgm−3s−1 with a relatively high value to cancel the velocity field in the solid PCM and b = 0.005 to avoid a division by zero. The dimensionless governing equations involve the well-known dimensionless numbers:

Pr = Jm,lαm,l, Ra = gβl03ΔTJm,lαm,l,Ste = cml,pΔTΔH,α¯=ααm,l,f = 1τ1+ τ2 ,ΔT = 3Q¯km,l (8)

The temperature and heat flux density are continuous functions at all interfaces: wall-heat sources, wall-PCM, heat sources-PCM, fin-PCM and fin-wall. The solid surfaces are regarded as impermeable and for which the no-slip condition applies. Seeking solutions of the coupled governing partial differential equations requires the knowledge of initial and boundary conditions. The initial conditions for the system under consideration are:

θ=U=V=0 (9)

The boundary conditions are given in the following way:

Adiabatic walls:

∂θ∂η|wall=0;η⊥wall (10)

Wall-heat source interface:

θs=θcandKc∂θs∂X|X=0=Ks∂θc∂X|X=0 (11)

Wall-PCM interface:

θs=θmandKs∂θs∂X|X=0=Km∂θm∂X|X=0 (12)

Wall-fin interface:

θs=θfandKs∂θs∂X|X=0=Kf∂θf∂X|X=0 (13)

Heat sources–PCM interface:

θc=θmandKc∂θ∂η|c=Km∂θ∂η|m (14)

Heat fin–PCM interface:

θf=θmandKf∂θ∂η|c=Km∂θ∂η|m (15)

Wall:

U=V=0 (16)

Here, the quantity η refers to the coordinate normal to the heat source-PCM interface. The above governing equations are discretized using the finite volume method (FVM) in a staggered mesh, which consists of M and N nodes in X and Y directions, respectively. The power law scheme is used for the evaluation of the total flux combining both advective and conductive terms. The SIMPLE routine is used to couple pressure and velocity equations [13]. The energy equation (Eq. (5)) is solved to determine the temperature field in all parts of the enclosure (PCM, wall, heat sources and fin). Note that the source term Sϴ = 0, except in the PCM, where Sθ=−1Ste∂fl∂τ . The energy equation for PCM is solved using the enthalpy fixed-grid technique developed by Voller et al. [14]. The source term Sϴ, which is the central feature of this technique, keeps track of the latent heat evolution; and its driving element is the local liquid volume fraction fl given by:

{fl =  1               if  θ> 0fl =  0              if  θ< 00 <fl< 1     if  θ =  0 (17)

In the numerical implementation, the value of liquid volume fraction fl is determined iteratively from the solution of the energy equation. A tri-diagonal matrix iterative method is used to solve the algebraic equations for U, V and θ. The iterative procedure is applied until reaching a well appropriate converging solution for the flow and energy fields at each time step, that is, when the criteria relating mass and energy balances are smaller than 10−9 and 10−3, respectively. Numerical calculations were performed to check the grid size and the time step effects on fl and ϴm. The results reveal that the grid size set to 50 × 60 and the time step to 4.509 × 10−4 ensure the best compromise between the execution time and the accuracy of the results.

3  Results and Discussion

The numerical predictions of the model developed in this study are compared to experimental data obtained by Casano et al. [15] who numerically and experimentally studied a one-dimensional phase-change process dominated by heat conduction. In their experimental arrangement sketched in Fig. 3, a plane slab of PCM (n-octadecane) was heated from above by a periodically heat flux density qH( 1 +  cos 2 πpe t ), with qH=8 W/mandpe=4 h, while the bottom cold wall was maintained at Tc = 15°C. Thermocouples placed at eight locations (see Table 1) are used to measure the temperature distribution. The first and the last ones are put in contact with the cold and hot walls, respectively.

Figure 3: Schematic of the plan slab of PCM

Fig. 4 displays our predicted and their measured temperatures corresponding to the eight locations in the PCM and a rather nice agreement is obtained.

Figure 4: Comparison between computed (lines) and measured (in symbols) temperature distributions vs. time for a period of 4 h. Symbols are data from [15]

Numerical simulations are performed to highlight the effect of the fin position on the thermal behavior of the PCM enclosure. The PCM used for simulations is the n-eicosane. The geometrical parameters and the thermophysical properties of the components of the physical model are listed in Tables 2 and 3. The dimensionless frequency of the pulsed power generated in each heat source is f = 35.76 (period Pe = 310 s). We have carried out our simulations for five cases one without a fin and the remaining four cases correspond to four different positions for a fin position:

•   P1: below the lower heat source.

•   P2: between the lower and the central heat sources.

•   P3: between the central and the upper heat sources.

•   P4: above the upper heat source.

Fig. 5 displays the dimensionless maximum temperature of the heat sources against the normalized time for four fin positions. The temperature oscillates due to the periodic nature of the power generated in the heat sources and its maximum is reduced by the presence of the fin. During the first six cycles, the largest maximum temperature is the one that is associated with the PCM enclosure without a fin. The numerical analysis also shows that the time evolution of the maximum temperature is characterized by three different regimes. The first regime corresponds to the first two periods and is characterized by an increase in both the minima and maxima of the maximum temperature for all the examined cases of the fin positions. During this regime, which also corresponds to the inception of the melting process, the heat transfer occurring in the liquid PCM layer is mainly dominated by conduction. The second regime starts when the liquid PCM zone largely expands and the natural convection movement takes place within the melt state. The growth of the liquid layer with time intensifies the natural convection and results in an improvement of the heat transfer extracted from the heat sources. Note that during this regime, the difference between the minima associated with the different fin positions decreases from one cycle to another. The same is also detected in the maxima of the maximum temperature. The analysis of the figure also shows that the quasi-steady state is reached after five cycles. As the time progresses, the PCM continues to melt and stores thermal energy essentially as a latent heat. The liquid layer expands and the liquid volume fraction increases with time as illustrated in Fig. 6. After ten cycles (τ = 0.34), the liquid volume fraction reaches 70% of the total PCM volume, and the melt stores thermal energy, not only in a latent form but also in a sensible one. This results in an increase in the average temperature of the liquid PCM and hence in a reduction in the temperature difference between the liquid PCM and the heat emanating from different surfaces (heat sources and substrate surfaces). Therefore, the extracted heat transfer from heat sources drops and this leads to an augmentation in both the maxima and minima of the maximum temperature of the heat sources.

Figure 5: Effect of fin position on the temporal variation of the dimensionless maximum temperature

Figure 6: Effect of the fin position on the temporal variation of the liquid fraction

Among all the fin positions considered in this study, the position of the fin placed between the central and the upper heat sources (P3) corresponds to the lowest values for both the minima and maxima of the maximum temperature generated by the heat sources that have the ability to evacuate more heat compared to the other position cases. In addition, such a position results in the slowest melting process corresponding tothe relatively smallest melting time. These findings might be of great importance for electronics cooling applications. Moreover, Fig. 6 demonstrates that the melting time of the PCM is almost the same for the three already defined positions P1, P2 and P4 and for the case of an enclosure without a fin.

4  Conclusions

We have investigated, through a modelling effort, numerical analysis and validation, the effect of the fin position on the thermal behavior and the cooling performance of a heat sink consisting of a PCM (n-eicosane) filling a rectangular enclosure equipped with three-pulsed discrete heat sources. The presence of a fin appears to be favorable for cooling electronic modules. The main outcomes of this study areas as follows:

•   The time evolution of the dimensionless maximum temperature of the heat sources exhibits three distinct regimes;

•   The lowest value of the maximum temperature and the largest melting time are obtained for a fin located between the central and the upper heat sources, while the opposite is obtained in the case of an enclosure without a fin.

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.