Numerical modelling of a melting-solidification cycle of a phase-change material with complete or partial melting

A high accuracy numerical model is used to simulate an alternate melting and solidification cycle of a phase change material (PCM). We use a second order (in time and space) finite-element method with mesh adaptivity to solve a single-domain model based on the Navier-Stokes-Boussinesq equations. An enthalpy method is applied to the energy equation. A Carman-Kozeny type penalty term is introduced in the momentum equation to bring the velocity to zero inside the solid region. The mesh is dynamically adapted at each time step to accurately capture the interface between solid and liquid phases, the boundary-layer structure at the walls and the multi-cellular unsteady convection in the liquid. We consider the basic configuration of a differentially heated square cavity filled with an octadecane paraffin and use experimental and numerical results from the literature to validate our numerical system. The first study case considers the complete melting of the PCM (liquid fraction of 95%), followed by a complete solidification. For the second case, the solidification is triggered after a partial melting (liquid fraction of 50%). Both cases are analysed in detail by providing temporal evolution of the solid-liquid interface, liquid fraction, Nusselt number and accumulated heat input. Different regimes are identified during the melting-solidification process and explained using scaling correlation analysis. Practical consequences of these two operating modes are finally discussed.     

Lattice Boltzmann investigation for enhancing the thermal conductivity of ice using Al2O3 porous matrix

An enthalpy-based Lattice Boltzmann method (LBM) with double-distribution function (DDF) model is used to investigate numerically the effects of inserting a porous matrix on the heat transfer performance of the phase change material (PCM). Simulations are carried out for melting of ice in saturated Al₂O₃ porous matrix encapsulated in a concentric annulus. The process is considered as a conduction/convection controlled phase change problem at a representative elementary volume (REV) scale. The present results are validated by previous published numerical simulations of melting with and without porous media. In this research paper, the effects of decreasing the porosity on the temperature contours, flow patterns within the melt zone, complete melting time of the PCM and average Nusselt number are discussed qualitatively and quantitatively.     

Thermophysical and heat transfer properties of phase change material candidate for waste heat transportation system

A waste heat transportation system––trans-heat (TH) system––is quite attractive that uses the latent heat of a phase change material (PCM). The purpose of this paper is to study the thermophysical properties of various sugars and sodium acetate trihydrate (SAT) as PCMs for a practical TH system and the heat transfer property between PCM selected and heat transfer oil, by using differential scanning calorimetry (DSC), thermogravimetry-differential thermal analysis (TG-DTA) and a heat storage tube. As a result, erythritol, with a large latent heat of 344 kJ/kg at melting point of 117°C, high decomposition point of 160°C and excellent chemical stability under repeated phase change cycles was found to be the best PCM among them for the practical TH system. In the heat release experiments between liquid erythritol and flowing cold oil, we observed foaming phenomena of encapsulated oil, in which oil droplet was coated by solidification of PCM.     

Effects of the heat transfer fluid velocity on the storage characteristics of a cylindrical latent heat energy storage system: a numerical study

A numerical study of the effects of the thermal fluid velocity on the storage characteristics of a cylindrical latent heat energy storage system (LHESS) was conducted. Due to the low thermal conductivity of phase change materials (PCMs) used in LHESS, fins were added to the system to increase the rate of heat transfer and charging. Finite elements were used to implement the developed numerical method needed to study and solve for the phase change heat transfer (melting of PCM) encountered in a LHESS during charging. The effective heat capacity method was applied in order to account for the large amount of latent energy stored during melting of the PCM and the moving interface between the solid and liquid phases. The effects of the heat transfer fluid (HTF) velocity on the melting rate of the PCM were studied for configurations having between 0 and 18 fins. Results show that the overall heat transfer rate to the PCM increases with an increase in the HTF velocity. However, the effect of the HTF velocity was observed to be small in configurations having very few fins, owing to the large residual thermal resistance offered by the PCM. However, the effect of the HTF velocity becomes more pronounced with addition of fins; since the thermal resistance on the PCM side of the LHESS is significantly reduce by the large number of fins in the system.     

Effects of the number and distribution of fins on the storage characteristics of a cylindrical latent heat energy storage system: a numerical study

A numerical study of the effects of the number and distribution of fins on the storage characteristics of a cylindrical latent heat energy storage system (LHESS) was conducted. Due to the low thermal conductivity of phase change materials (PCMs) used in LHESS, fins were added to the system to increase the rate of heat transfer and charging. Finite elements were used to implement the developed numerical method needed to study and solve for the phase change heat transfer (melting of PCM) encountered in a LHESS during charging. The effective heat capacity method was applied in order to account for the large amount of latent energy stored during melting of the PCM and the moving interface between the solid and liquid phases. The effects of increasing the number and distribution of fins on the melting rate of the PCM were studied for configurations having between 0 and 27 fins for heat transfer fluid (HTF) velocities of 0.05 and 0.5?m/s. Results show that the overall heat transfer rate to the PCM increases with an increase in the number of fins irrespective of the HTF velocity. It was also observed that the total amount of energy stored after 12?h increases nearly linearly with the addition of fins up to 12 fins; further addition of fins increasing the total energy stored by ever smaller amounts.     

Conjugate Phase Change Heat Transfer in an Inclined Compound Cavity Partially Filled with a Porous Medium: A Deformed Mesh Approach

In this paper, the melting process of a PCM inside an inclined compound enclosure partially filled with a porous medium is theoretically addressed using a novel deformed mesh method. The sub-domain area of the compound enclosure is made of a porous layer and clear region. The right wall of the enclosure is adjacent to the clear region and is subject to a constant temperature of T_c. The left wall, which is connected to the porous layer, is thick and thermally conductive. The thick wall is partially subject to the hot temperature of T_h. The remaining borders of the enclosure are well insulated. The governing equations for flow and heat transfer, including the phase change effects and conjugate heat transfer at the thick wall, are introduced and transformed into a non-dimensional form. A deformed grid method is utilized to track the phase change front in the solid and liquid regions. The melting front movement is controlled by the Stefan condition. The finite element method, along with Arbitrary Eulerian–Lagrangian (ALE) moving grid technique, is employed to solve the non-dimensional governing equations. The modeling approach and the accuracy of the utilized numerical approach are verified by comparison of the results with several experimental and numerical studies, available in the literature. The effect of conjugate wall thickness, inclination angle, and the porous layer thickness on the phase change heat transfer of PCM is investigated. The outcomes show that the rates of melting and heat transfer are enhanced as the thickness of the porous layer increases. The melting rate is the highest when the inclination angle of the enclosure is 45°. An increase in the wall thickness improves the melting rate.

     

Analytical approximation for solidification processes in PCM storage with internal fins: imposed heat flux

Uses of thermal energy storage systems have expanded notably in recent decades. In thermal energy systems, internal heat transfer enhancement techniques such as fins are often used because of the low thermal conductivity of the phase change materials (PCMs). In this paper, solidification of a PCM is studied in a rectangular storage with horizontal internal plate fins and an imposed constant heat flux on the vertical walls. A simplified analytical solution is presented and its results are compared to those for a numerical approach based on an enthalpy method. The fraction of solidified PCM in storage is calculated with the derived analytical model which determines how much of the storage is solidified after a certain time. The results show that the analytical model satisfactorily estimates the solid–liquid interface and the temperature distribution for the fin, which are useful in the design of PCM-based thermal energy storage or cooling systems.     

Heat transfer enhancement of PCM melting in 2D horizontal elliptical tube using metallic porous matrix

In this study, the melting process of ice as a phase-change material (PCM) saturated with a nickel–steel porous matrix inside a horizontal elliptical tube is investigated. Due to the low thermal conductivity of the PCM, it is motivated to augment the heat transfer performance of the system simultaneously by finding an optimum value of the aspect ratio and impregnating a metallic porous matrix into the base PCM. The lattice Boltzmann method with a double distribution function formulated based on the enthalpy method, is applied at the representative elementary volume scale under the local thermal equilibrium assumption between the PCM and porous matrix in the composite. While reducing or increasing the aspect ratio of the circular tubes leads to the expedited melting, the 90\(^{\circ }\) inclination of each elliptical tube in the case of the pure PCM melting does not affect the melting rate. With the reduction in the porosity, the effective thermal conductivity and melting rate in all tubes promoted. Although the natural convection is fully suppressed due to the significant flow blockage in the porous structure, the melting rates are generally increased in all cases.     

Numerical simulation and experimental analysis of heat transfer through the neck tube into vertical cryogenic insulated cylinders

The neck tube is an important support structure in cryogenic insulated cylinders. The heat flux from the outside environment through the neck tube into the cryogenic liquid occupies a great proportion of the total heat leak and can be more than half of the total heat loads. In this paper, conjugate convective-conductive heat transfer model between wall and the cold vapor in conditions of natural discharge is numerically investigated. Also a liquid nitrogen boil-off method was adopted in experiments to validate the result of numerical simulation. Experimental results indicate more favorable agreement with conjugate heat transfer (CHT) model compared with simple solid heat conduction (SSHC) model by ANSYS software. And the convection between the wall and vapor is also calculated. The research and results can provide reference in design for neck tube of the cryogenic cylinder.     

Use of experiment and an inverse method to study interface heat transfer during solidification in the investment casting process

A technique to determine the thermal boundary conditions existing during the solidification of metallic alloys in the investment casting process is presented. Quantitative information about these conditions is needed so that numerical models of heat transfer in this process produce accurate results. In particular, the variation of the boundary conditions both spatially and temporally must be known. The method used involves the application of a new inverse heat conduction method to thermal data recorded during laboratory experiments of aluminium alloy solidification in investment casting shell moulds. The resultant heat transfer coefficient for the alloy/mould interface is calculated. An experimental programme to determine requisite mould thermal properties was also undertaken. It was observed that there is significant variation of the alloy/mould heat transfer coefficient during solidification. It is found to be highly dependent on the alloy type and on the vertical position below the initial free surface of the liquid metal. The aluminium casting alloys used in this study were 413, A356, 319 (Aluminum Association designations), and commercially pure aluminium. These alloys have significantly different freezing ranges. In particular, it was found that alloys with a high freezing range solidify with rates of heat transfer to the mould which are very sensitive to metallostatic head.     

Modeling thermal insulation of firefighting protective clothing embedded with phase change material

Experiments and research on heat transport through firefighting protective clothing when exposed to high temperature or intensive radiation are significant. Phase change material (PCM) takes energy when changes from solid to liquid thus reducing heat transmission. A numerical simulation of heat protection of the firefighting protective clothing embedded with PCM was studied. We focused on the temperature variation by comparing different thicknesses and position conditions of PCM combined in the clothing, as well as the melting state of PCM and human irreversible burns through a simplified one-dimensional model. The results showed it was superior to place PCM between water and proof layer and inner layer, in addition, greater thickness increased protection time while might adding extra burden to the firefighter.     

Experimental study of enhanced heat transfer by addition of CuO nanoparticle

An energy storage system has been designed to study the thermal characteristics of paraffin wax with an embedded nano size copper oxide (CuO) particle. This paper presents studies conducted on phase transition times, heat fraction as well as heat transfer characteristics of paraffin wax as phase change material (PCM) embedded with CuO nanoparticles. 40?nm mean size CuO particles of 2, 5 and 10% by weight were dispersed in PCM for this study. Experiments were performed on a heat exchanger with 1.5–10?l/min of heat transfer fluid (HTF) flow. Time-based variations of the temperature distributions are revealed from the results of observations of melting and solidification curves. The results strongly suggested that the thermal conductivity enhances 6, 6.7 and 7.8% in liquid state and in dynamic viscosity it enhances by 5, 14 and 30% with increasing mass fraction of the CNEPs. The thermal conductivity ratio of the composites can be augmented by a factor up to 1.3. The heat transfer coefficient during solidification increased about 78% for the maximum flow rate. The analysis of experimental results reveals that the addition of copper oxide nanoparticles to the paraffin wax enhances both the conduction and natural convection very effectively in composites and in paraffin wax. The paraffin wax-based composites have great potential for energy storage applications like industrial waste heat recovery, solar thermal applications and solar based dynamic space power generation with optimal fraction of copper oxide nanoparticles.     

Conjugate heat transfer in square enclosures

Building elements represented by square vertical enclosures encircled with finite walls or with centered solid body, could maintain the equivalent fluid volumes through the volume ratio scale. Present work aims to investigate the fluid flow and heat transfer in these two building elements. Complete two-dimensional numerical simulation of the conjugate heat conduction and natural convection occurring in both enclosures is carried out. An analytical expression for the minimum size of the inserted body at which the body begins to suppress the natural convection flow is proposed and validated by the numerical results. The fluid flow and heat transfer characteristics are analyzed through the streamlines, heatlines, and total heat transfer rates across both enclosures. Results reveal that heat transfer rates across both enclosures are complex functions of the volume ratio scale, Rayleigh number, and the relative thermal conductivity.     

On cooling behavior of a vertical plate in a phase change material/water composite enclosure under pulsating heat load

This paper presents a numerical and experimental study concerning cooling characteristics of a pulsating heated vertical plate sandwiched between a substrate of phase change material (PCM) and an enclosure of water, forming a composite vertical rectangular enclosure. The vertical plate is assumed to have a uniform pulsating (on/off) volumetric heat source. The PCM considered in the present study is n-Octadecane. In the finite-difference simulation, the two-dimensional buoyancy-driven fluid flows developed in both the water-filled subenclosure and the molten PCM region of the PCM-filled subenclosure were modeled as laminar Newtonian fluid flow adhering to the Boussinesq approximation. Meanwhile, two-dimensional conduction is accounted for the plate heater as well as the solid PCM zone. Numerical results are presented to unveil the cooling behavior of the pulsating heat-generating plate through the PCM substrate and water-filled enclosure. Results of the parametric simulations reveal that the water layer has the better capability of heat dissipation than the PCM substrate. Heat dissipation from the plate through the PCM substrate is mainly via the latent heat absorption as associated with melting phenomenon. Moreover, numerical results obtained are compared with the corresponding experiments.     

Improving energy storage by PCM using hybrid nanofluid [(SWCNTs-CuO)/H2O] and a helical (spiral) coil: Hybrid passive techniques

The aim of this study is the numerical analysis of the melting process of the phase change material (PCM) in a spiral coil. The space between the inner tube and outer shell is filled with RT-50 as PCM. Moreover, the hybrid nanofluid (with a carbon component) flows through the inner tube. The novelty of this work is to use different configurations of fin and different percentage of hybrid nanoparticles (SWCNTs-CuO) on the PCM melting process. In the numerical model created by ANSYS-Fluent, the effect of various inlet temperatures is investigated. The results indicate that the extended surface created by extra fin has a dominant effect on melting time, so by adding the third fin, the melting time is reduced by 39.24%. The next most influential factor in PCM melting is the inlet temperature of the working fluid, so that 10 °C increment of temperature result in the PCM melting time decreased by 35.41%.     

Evaporation of Falling and Shear-Driven Thin Films on Smooth and Grooved Surfaces

One of the most important tasks in development of modern gas turbine combustors is the reduction of NO_x emissions. An effective way to reduce the NO_x emission is using the lean premixed prevaporization (LPP) concept. An important phenomenon taking place in LPP chambers is the evaporation of thin fuel films. To increase the fuel evaporation rate, the use of microstructured walls has been suggested. The wall microstructures make use of the capillary forces to evenly distribute the liquid fuel over the wall, so that the appearance of uncontrolled dry patches can be avoided. Moreover, the wall structures promote the thin film evaporation characterized by ultra-high evaporation rates. An experimental setup was built for the investigation of thin liquid films falling down on the outer surface of vertical tubes with either a smooth or structured surface. In the first testing phase water is used, fuel like liquids will be used later on. The thin film can be heated from both sides, by hot oil flowing inside the tube, and by hot compressed air flowing in co-current direction to the thin film. The film is partly evaporated along the flow. Results for the wavy film structure at different Reynolds numbers are reported. For theoretical investigations a model describing the hydrodynamics and heat transfer due to evaporation of the gravity- and shear-driven undisturbed liquid film on structured surfaces was developed. For low Reynolds numbers or low liquid mass fluxes the wall surface is only partly covered with liquid and the heat transfer is shown to be governed by the evaporation of the ultra-thin film in the vicinity of the three-phase contact line. A numerical model for the solution of a two-dimensional free-surface flow of a liquid film over a structured wall was also developed. The Navier–Stokes equations are solved using the Volume of Fluid (VOF) technique. The energy equation is included in the model. The model is verified by comparison with data from the literature showing favorable agreement. In particular, the proposed model predicts the formation of capillary waves observed in the experiments. The model is used to investigate the flow of liquid on a structured wall. This calculation is the first step towards the modeling of a three-dimensional wavy flow of a gravity- and shear-driven film along a wall with longitudinal grooves. It is found that due to the Marangoni effect, a circulating flow arises within the cavity, thereby leading to an enhancement in the evaporation rate.     

Thermal protection characteristics of a vertical rectangular cell filled with PCM/air layer

The present paper presents a numerical analysis concerning thermal protection characteristics of a vertical rectangular composite cell filled with a solid-liquid phase change material (PCM) and air layer. Inside the composite cell the PCM layer is separated from air layer by a solid partition of negligible thickness. The buoyancy-induced flows developed in both the air-filled layer and the molten PCM zone inside the PCM layer were modeled as two-dimensional laminar Newtonian fluid flow adhering to the Boussinesq approximation. Meanwhile, two-dimensional conduction heat transfer was accounted for the unmelted solid PCM region. Delineation is made via a parametric simulation of the effects of the pertinent parameters:Ste (Stefan number),Sc (subcooling factor),Ra (Rayleigh number), aspect ratio of composite cell,A, and relative thickness ratioA _p /A _a , on the transient thermal protection performance of the composite cell. Results demonstrate that by means of the latent-heat absorption inside the PCM layer, heat penetration across the composite cell can be greatly retarded over an effective duration until a critical instant until the melting front of PCM reaches the partition wall. Such an effective thermal protection duration is found to be a strong function ofRa, Ste, A _p /A _a , andA. In addition, the results of the transient heat transfer rate penetrating through the composite cell are examined as a function of the pertinent parameters of the problem.     

Effects of anisotropy in permeability on the two-phase flow and heat transfer in a porous cavity

This paper reports on the results of a numerical study of convection flow and heat transfer in a rectangular porous cavity filled with a phase change material under steady state conditions. The two vertical walls of the cavity are subject respectively to temperatures below and above the melting point of the PCM while adiabatic conditions are imposed on the horizontal walls. The porous medium is characterized by an anisotropic permeability tensor with the principal axes arbitrarily oriented with respect to the gravity vector. The problem is governed by the aspect ratioA, the Rayleigh numberRa, the anisotropy ratioR and the orientation angle θ of the permeability tensor. Attention is focused on these two latter parameters in order to investigate the effects of the anisotropic permeability on the fluid flow and heat transfer of the liquid/solid phase change process. The method of solution is based on the control volume approach in conjunction with the Landau-transformation to map the irregular flow domain into a rectangular one. The results are obtained for the flow field, temperature distribution, interface position and heat transfer rate forA=2.5,Ra=40, 0≤θ≤π, 0.25≤R≤4. It was found that the equilibrium state of the solid/liquid phase change process may be strongly influenced by the anisotropy ratioR as well as by the orientation angle θ of the permeability tensor. First, for a given set of parametersA,Ra andR, there exists an optimum orientation θ_max for which the flow strength, the liquid volume and the heat transfer rate are maximum. There also exists an orientation θ_min=θ_max+π/2 for which these quantities are minimum. Second, when an anisotropic medium is oriented along the optimum direction θ_max, an increase of the permeability component along that direction will increase the flow and heat transfer rate in a same order while an increase of the other permeability component only has a negligible effect. For the parameter ranges considered in the present study, it was found that the optimum direction is lying between the gravity vector and the dominant flow direction.     

Network simulation method for solving phase-change heat transfer problems with variable thermal properties

 The transient heat conduction equation in a finite slab undergoing phase change (two-phase problem of melting and solidification), with isothermal, adiabatic or convective boundary conduction is studied by the network simulation method; solid phase conductivity and specific heat are assumed to be dependent on temperature. Ablation, as a particular case, is also analysed. A network model is established for a cell and boundary conditions are added to complete the whole network model. No restrictions exist, as to the kinds of linear and non-linear boundary conditions, Stefan number values or the initial conditions (when hypotheses concern of the Stefan problem, numerical and exact solutions are compared for a large interval of Stefan numbers; simulation values show good agreement). Movement of the solid–liquid boundary and thermal fields are determined in all cases. Received on 10 May 2000 / Published online: 29 November 2001     

Thin-film phase change driven by disjoining pressure difference

The thin liquid film at the contact line is gaining increasing attention due to its importance in phase-change heat transfer and wettability control. The intermolecular effects within the thin film are usually represented by of disjoining pressure. In the present work, molecular dynamics were employed to investigate the influence of disjoining pressure on thin–film evaporation and condensation in a nanoscale triple-phase system. The simulation domain was a cuboid with argon gas sandwiched between the two solid platinum walls. The two solid walls were fixed at the same temperature while the liquid films had different initial thicknesses thus corresponding to different disjoining pressures. Spontaneous evaporation and condensation were observed at the thicker and the thinner film, respectively. Disjoining pressure together with thermodynamics theories were employed to qualitatively and quantitatively explain the phenomenon. The evaporation fluxes were measured and compared to the Hertz-Knudsen-Schrage model which is based on the kinetic theory of gases. The resulting non-evaporating thickness was measured and compared to the models based on disjoining theory.