Enthalpy-porosity Method Research Articles

Mathematical modeling of conductive–convective melting with a two-phase region by the enthalpy-porosity method

Mathematical modeling of conductive–convective melting with a two-phase region by the enthalpy-porosity method

Fins-nanoparticle combination for phase change material enhancement in a triplex tube heat exchanger: A numerical approach to thermal sustainability

Many phase change materials exhibit low thermal conductivity, leading to incomplete melting and solidification processes. To address this issue, researchers have numerically explored the integration of alumina nanoparticles (Al2O3) with paraffin (RT82), a phase change material with a solidification temperature of 65 °C, in a triplex-tube heat exchanger. The phase change material model with internal longitudinal fins employed the both-sides freezing technique and was conducted using Ansys Fluent software, employing the enthalpy-porosity method within a finite-volume framework to model the phase change material behavior during both melting and solidification phases. This methodology aims to significantly improve the thermal performance of thermal energy storage systems by enhancing heat transfer efficiency within the phase change material (PCM), thereby ensuring more effective utilization of stored thermal energy. The numerical findings show that the pure PCM completely solidified in 780 min. By dispersing 1 %, 4 %, 7 %, and 10 % of Al2O3 in the PCM, thermal conductivity improved by 3 %, 12.5 %, 22.5 %, and 32.5 %, respectively. Additionally, the inclusion of nano-PCM reduced the solidification time. This research also compares the overall energy release in two different situations: PCM with and without nanoparticles. The computer-generated simulation results closely correlated with the practical experimentation results.

Numerical analysis of the combined influence of fin shape and location on constrained melting of phase change materials in a spherical capsule with double fins

AbstractThis study presents a novel approach to investigating the combined influence of fin position and shape on the constrained melting behavior of phase change material (PCM) within a spherical capsule (S.C.) through numerical analysis. Unlike previous research, which predominantly focused on single fin shapes or positions, this work uniquely explores the impact of double, simple, and easily manufacturable fin shapes. A two‐dimensional computational model employing the enthalpy–porosity method assesses melting behavior, temperature distribution, and PCM flow. Numerous fin shapes, namely rectangular, trapezoidal converging, trapezoidal diverging stepped, inverse stepped, and triangular, are considered in the analysis. The study reports the influence of the location of two identically shaped fins on the thermal performance. The fins' cross‐sectional area and base thickness are kept equal in all cases. The thermal performance of an S.C.‐integrated fin system is evaluated by analyzing various attributes such as total saving in the duration of melting, enhancement ratio, and Nusselt number. The results indicate that the position of the fins has a more significant impact on melting performance than the fin shape. The best performance is achieved when fins are placed in the lower half of the capsule, followed by the center and upper halves, regardless of fin shape. For rectangular fins, shifting the position of the fin from the bottom half to the center increases the melting time by 24.7% and the top half by 68.3%. The shortest melting time of 93 min is observed for lower‐half rectangular fins, followed by center‐placed triangular fins (94 min). This study offers a theoretical foundation for optimizing the performance of different technologies using latent heat thermal energy storage systems such as packed‐bed, cascaded thermal energy storage systems.

Impact of key parameters on the numerical modelling of Phase-Changing phenomena in horizontal latent heat thermal energy storage

Numerous studies have focused on numerical investigations of the performance of the shell-and-tube Latent Heat Thermal Energy Storage (LHTES). However, many studies have often overlooked rigorous assessment of key parameters’ impact on the numerical simulation accuracy relative to experimental results. Therefore, this paper aims to develop a numerical model to track the solid–liquid interface during phase-changing phenomena in a horizontal LHTES system and investigate the effects of key parameters on the simulation accuracy. For this purpose, an experimental facility was built, and the numerical model was developed in ANSYS using the enthalpy-porosity method. During the simulation of the computational model, consideration was given to thermophysical properties (constant versus variable) of phase change material, contact resistance (zero versus nonzero), and the mushy zone constant (105, 5 × 105, and 106). The simulation results were compared with experimental data obtained from the in-house studies. The results showed that, compared to using constant thermophysical properties, employing variable thermophysical properties more accurately tracked the solid–liquid interface during melting. It substantially improved the simulation accuracy with a Root Mean Square (RMS) error less than 5 % for most of the temperature-measuring locations. Conversely, variable properties did not make a significant difference in the solidification process. The contact resistance due to the air gap was determined to be 0.035 m2·K/W for the current study. Furthermore, incorporating this resistance resulted in a more accurate approximation of the solid–liquid interface during melting, with RMS errors below 5 % across all locations. The study comparing different mushy zone constants recommended utilizing a slightly higher value of 106 for melting, whereas its impact on solidification was found to be insignificant.

Effectiveness of a PCM-based heat sink with partially filled metal foam for thermal management of electronics

Electronics often experience overheating, which occurs when insufficient or improper cooling causes temperatures to rise quickly, potentially leading to component failure. Recently, a passive technique utilizing phase change material (PCM)-based heat sinks is suitable for modern electronic cooling. However, due to its low thermal conductivity, Aluminum (Al) foam as a thermal conductivity enhancer (TCE) is used. In this study, a partial filling approach is employed, involving a heat sink designed to cool a protruding electronic component (EC) attached to a metal fin. The enthalpy-porosity method and the thermal equilibrium model are applied. Results show that the efficiency of the heat sink is improved by increasing the filling ratio of the Al foam, thereby reducing the EC temperature by 15.85 °C, shortness the melting time by more than 1000 s, and the overall effective conductivity enhanced by 25 times compared to the case without Al foam. In addition, the optimal partial filling ratio is determined to be 2/3, due to the considerable savings of 33.33 % in material cost and overall heat sink weight. The findings suggest that partial filling is a strategy that can be employed to achieve a more comprehensive performance in thermal management using PCM-based heat sinks.

Heat release efficiency Betterment inside a novel-designed latent heat exchanger featuring arc-shaped fins and a rotational mechanism via numerical model and artificial neural network

Thermal energy storage (TES) units featuring phase change materials (PCMs) have a crucial impact on efficient energy management. Since energy demands keep changing, and the world increasingly relies on sustainable energy sources, innovation is continuously required. New technologies, materials, and methods for energy storage must be developed to meet these shifting needs effectively. In this study, a triplex-tube heat exchanger (TTHE) was employed as a TES medium, utilizing PCM within the middle tube. Because of its high thermal energy storage capacity, RT50 was selected as the PCM. Four innovative arc-shaped fins were strategically integrated within the PCM space to accelerate heat release. The system was subjected to rotational speeds of 0.1, 0.3, 0.5, 1, and 1.5 rpm. The investigation was carried out in two stages: initially, the impact of varying rotational speeds on the discharging process of the PCM in a finless TTHE was explored; subsequently, the influence of the same speeds on the solidification behavior of the finned TTHE was analyzed. Natural convection and PCM’s solidification process were examined through the enthalpy-porosity method. The findings represented that in the absence of fins, the device’s rotation had a more noticeable impact on the solidification behavior; however, incorporating fins was much more influential. Up to 2850 s, all the finned systems had solidified entirely in the presence or absence of the rotational mechanism, while in the finless system with a rotational speed of 1.5, 30.52 % of the material was still unsolidified until this time. Finally, the influence of three key structural parameters of the fins: α (outer fins arc angle), β (inner fins arc angle), and S (space between fins) on the discharging time of the PCM was analyzed using an artificial neural network (ANN) model. The study offered critical insights for optimizing fin configuration by systematically varying these parameters within defined ranges. An ANN predictive model was also proposed to assist TES system manufacturers and developers in future design and optimization efforts. The results demonstrated that the optimal setting of the TTHE (with the parameters of α = 60°, β = 60°, S = 10 mm) achieved the liquid fraction of 0.1, 80.33 % faster than the finless system.

Investigating swastika fins to enhance heat transfer and melting process of nano-enhanced phase change materials units

The ability of nano-enhanced phase change materials (NePCM) to store thermal energy has garnered the interest of scientists. An approach to enhancing thermal conductivity involves the implementation of novel fins that facilitate a more rapid PCM melting process. The present investigation assessed the efficacy of eight distinct fin configurations in augmenting the thermal properties of NePCM. Four occurrences featured swastika fins, whereas the remaining instances exhibited a 45-degree tilt angle. The enthalpy-porosity method was implemented to simulate heat transfer and melting process. The present research investigates the effects of fin tilt angles and swastika fins on the NePCM-thermal energy storage (TES) unit's mean temperature, liquid percentage, and thermal energy storage. Compared to the straight fins (Case 1), the long-arm swastika fins (Case 4) decreased the melting time by a factor of 2.3. In addition, Case 4 demonstrated exceptional thermal energy storage capabilities. Case 7 demonstrated that the tilting angle significantly affected the melting time, with a reduction of 34.9 %, compared to the non-tilted swastika fins utilized in Case 2. However, the increase in the fins' length was more substantial. According to the findings of this study, Case 4 is the optimal fin configuration for NePCM-TES.

Probing the melting dynamics in a phase change Rayleigh–Bénard system under low gravity conditions

The present work aims to investigate dynamic characteristics of Rayleigh–Bénard convection during the solid–liquid phase change process under the influence of variable low-gravity conditions. Low-gravity conditions are crucial to gain physical insights into the complex convection dynamics relevant in terrestrial and space environments that feature diverse gravitational fields. The melting dynamics of a paraffin-based phase change material with Prandtl number Pr ≈ 71 in a square rigid walled enclosure is analyzed by subjecting it to a fixed temperature difference resulting in a Stefan number Ste ≈ 0.33. The enthalpy porosity method is employed to simulate the solid-liquid phase change process in Rayleigh number range 105<Ra <107 to take wide range of gravitational acceleration (0.0g, 0.1g, 0.2g, 0.4g, 0.6g, 0.8g, and g) into account. The simulation results allow detailed insights on the buoyancy-driven convective flow phenomena in a phase-change Rayleigh–Bénard system. Various characteristic flow regimes, the criteria for the onset of convective instabilities, and scaling analysis are presented under varied gravitational acceleration. The findings reveal that the critical Rayleigh number (Racr) for the onset of convection is found to vary between 1351.5 and 1810.5 under the gravity range investigated in contrast to the Racr ≈ 1707.7 for classical Rayleigh-Bénard system under the standard gravity. Furthermore, the scale analysis based on the numerical results shows that the relationship between effective Nusselt number (Nue) and effective Rayleigh number (Rae) follows a power law behavior Nue = 0.27Rae0.25 under each gravity condition, similar to the classical Rayleigh–Bénard system.

Analysis of charging performance of thermal energy storage system with graded metal foam structure and active flip method

The latent heat thermal energy storage (LHTES) technology based on solid-liquid phase change material (PCM) is characterized by high energy storage density, small volume change, and constant operation temperature, which is widely employed in waste heat recovery, solar thermal utilization, and equipment thermal management. This paper introduces active flipping and graded porous to strengthen natural convection and heat conduction to alleviate the issue of low thermal conductivity and non-uniform phase change of PCM. Based on the fixed grid system, the enthalpy-porosity method combined with time-dependent gravity and non-Darcy porous effect is employed to solve the solid-liquid phase change problem under different porosity gradients, dimensionless flipping times (t*), and modified Ra (Ra*) numbers. Furthermore, the effects of the Ra* on the optimal t* and porosity gradient are also discussed. The results show that the flipping method can significantly enhance flow and heat transfer within the container, improve temperature field uniformity, and increase the proportion of latent heat storage. With the increase of porosity gradient and t*, the melting time initially decreases and then increases and the optimal thermal performance is achieved at a 6 % porosity gradient and t* = 0.4. Compared to the uniform porous structure without flipping, it can shorten melting time by 49.37 % and improve thermal energy storage rate (TESR) by 76.23 %. Additionally, the optimal porosity gradient is 6 % under different Ra*, but the optimal t* decreases with the increase of Ra*. This paper explores the charging performance of the thermal energy storage system with the graded metal foam structure and active flip method, which can contribute to the study of heat transfer enhancement in LHTES technology.

Thermal performance of cascaded latent heat thermal energy storage units constructed based on solid-liquid interface information

A novel approach to constructing cascaded latent heat thermal energy storage units is proposed to enhance the thermal performance of horizontally placed shell-and-tube phase change heat storage systems. This method is based on the solid-liquid interface variation patterns during the solidification and melting processes. It creates interfaces between different phase change materials within the cascaded latent heat thermal energy storage unit to shorten melting time and increase thermal performance by optimizing the unit's internal structure to balance heat transfer rates across all areas. The enthalpy-porosity method was used to numerically study the phase change phenomenon. The inner tube of the units was set as an isothermal boundary, and the outer shell was insulated. Findings from this cascaded latent heat thermal energy storage system demonstrate that, compared to single-latent heat thermal energy storage units, the average heat transfer rates during solidification and melting increased by 27.3% and 36.7%, respectively. The average heat transfer rate during solidification improved by 11.7% compared to traditional concentric circular cascaded latent heat thermal energy storage units. These results indicate that this method can significantly enhance the thermal performance of latent heat thermal energy storage units, providing valuable insights for future structural optimization.

Structural optimization on heat transfer characteristics in polyurethane foam film melting process by deflagration method

Flexible polyurethane foam (FPUF) is extensively utilized in various industrial sectors, due to its distinctive porosity. The deflagration approach is a prevalent method employed to reticulate polyurethane foam. The open porosity of the polymer, along with its reticular framework, significantly shape the structure of FPUF and consequently impact its mechanical properties. This study investigates the impacts of two pivotal parameters, pore size and film thickness, on the open porosity of the foam film. The heat transfer characteristics are studied by standard k-ε model coupling with the melting-solidification model based on the enthalpy-porosity method analysis model, and nine PUF models are selected as the evaluation objects. The experimental results indicate that the reticulation effect is more favorable when the open porosity ranges from 0.7 to 0.85. Based on these findings, the parameters with large variation degrees have a greater influence on the melting process of FPUF film. Furthermore, FPUF composites characterized by thicker films and smaller pores demonstrate a remarkable reticulation performance. The insights garnered from this study hold promise for advancing the comprehension of the underlying mechanisms governing FPUF deflagration, thereby facilitating the optimization of its structural applications in industrial domains.

Simultaneous close-contact melting on two asymmetric surfaces: Demonstration, modeling and application to thermal storage

The present study deals with melting in a geometry suitable for Latent-Heat Thermal Energy Storage (LHTES) systems, which are of importance for future industrial installations utilizing solar energy or waste heat. The intrinsically high thermal resistance of phase-change materials (PCM) can be remedied by taking advantage of close contact melting (CCM), where the solid phase is separated from a hot surface by only a thin liquid layer. The basic CCM takes place on a horizontal flat surface, but during the last decade its application in various finned LHTES units has been demonstrated experimentally and analyzed numerically.Specifically, the configuration studied in the present work is relevant to a horizontal double-pipe concentric storage unit with a longitudinally finned inner tube. While this rather simple geometry has been extensively studied in the past, the present work is completely novel in its demonstration and analysis of simultaneous close-contact melting on two generally asymmetric surfaces created by the longitudinal fins. A unique experimental apparatus is introduced, based on the so-called 'Mercedes' configuration of the fins. It is designed transparent, allowing for real-time observation of melting sequences. Close-contact melting is achieved by supplying heat to the outer shell of the unit: the solid phase is detached from the shell and moves, by translation and rotation, in the liquid phase. The results from these experiments demonstrate the feasibility of CCM on two asymmetric walls within this system, hinting at the potential for optimizing the melting rate by utilizing a larger portion of the extended surface for CCM.A key component of this study is the development of a reliable numerical approach, which goes beyond the limitations of widely-applied enthalpy-porosity method and combines the general enthalpy formulation, convective heat transfer and general rigid body motion that includes rotation. Therefore, a new numerical model has been devised, using an in-house code realized in MATLAB. A full set of the governing conservation equations is solved using a finite-difference framework, integrated with advanced numerical methods for the fluid-solid interaction and an enthalpy formulation for the phase change process. The model is validated carefully, and then numerical studies are conducted to elucidate the melting process.The present paper presents a further proof of the special role that close-contact melting (CCM) can play in properly designed thermal energy storage units and finned systems in general. It is demonstrated that the fins, when properly designed and oriented, can induce CCM in their vicinity, contributing to a very significant increase in the melting rate, which reflects charging of the unit.

Numerical investigation of an innovative through-tubing solution to thermal plug and abandonment of oil wells with thermite reactions

Plug and Abandonment (P&A) procedures are mandatory in the oil and gas industry. Conventional well-plugging methods typically involve the laborious and expensive process of cementing, requiring the removal of production tubing. In response to this challenge, a novel approach, known as Thermal Plug and Abandonment (TP&A), has been explored. TP&A proposes the introduction of an exothermic chemical reaction through the production tubing, generating substantial heat to melt the tubing intentionally. The passage formed by the melting process facilitates the traditional insertion of cement, eliminating the need for tubing removal. In this study, the TP&A process is investigated through numerical computations. The oil well structure is approximated as a two-dimensional axisymmetric domain with multiple layers representing different wellbore materials. A numerical code, incorporating chemical kinetics, phase change, and conjugate heat transfer models, was developed in the OpenFOAM® software. The thermite reaction is modeled using a zero-order kinetic model, and the phase change model employs the well-established enthalpy-porosity method to track material melting and solidification. The study primarily focuses on evaluating heat diffusion through the oil well structure during the TP&A process, with a central emphasis on investigating the melting of the production tubing. It was observed that compacting the mixture and diluting it with alumina up to a certain threshold enhanced the tubing’s melting. Reducing the initial mixture porosity from 0.55 to 0.4 has increased the tubing’s melting volume, constrained to the thermite height, from 60 to approximately 91%. Moreover, this study examined how diluting the thermite mixture with inert alumina affectsthe heat transfer and, consequently, the tubing’s melting. The findings indicated that a 20% dilution can enhance the tubing’s melting volume by up to 87%.

Improving the performance of a shell and tube latent heat thermal energy storage through modifications of heat transfer pipes: A comprehensive investigation on various configurations

The modification of the geometric configurations of heat transfer pipes in shell and tube Latent Heat Thermal Energy Storage (LHTES) systems not only enhances the melting process of the phase change material (PCM) but also improves the overall performance of these systems. This study aims to investigate ways to enhance the performance of LHTES systems by employing heat transfer pipes with various fin and twisted tape arrangements in a horizontal orientation. The Finite Volume Method and Enthalpy-Porosity method are employed to simulate the melting process. Stearic acid is used as the PCM material, while water serves as the heat transfer fluid. Eight different geometric configurations are modelled in the LHTES system: base case, horizontal fins, vertical fins, helical fins, horizontal tape, vertical tape, twisted tape and helical fins with twisted tape. The results show that within the time range of 0 and 29 min, the combined configuration of helical fins with twisted tape consistently demonstrates the fastest melting process. After 29 min, the configuration with vertical fins exhibits a marginally faster melting process than the combined configuration of helical fins with twisted tape. The configurations involving tapes also contribute to accelerated melting, although to a lesser extent than those with fins. Particularly, twisted tape proves highly effective in facilitating faster melting. The complete melting process times for configurations with vertical fins, helical fins, and combined helical fins with twisted tape are 38.7 %, 23.5 % and 32.7 % faster compared to the base case which is ∼69 min. Among the configurations, using tapes results in higher flow resistance and surface area compared to the base case. The attractive features of these configurations make them ideal for creating efficient and space-saving energy storage systems. This study provides crucial insights into essential heat and mass transfer processes, which can be leveraged to develop advanced LHTES systems for enhanced performance and sustainable energy solutions.

On the effect of nanoparticles in different PCM configurations in the thermal performance of a packed bed energy storage system

When it comes to storing energy and releasing it when needed, packed bed storage systems are considered to be one of the most efficient systems available. This paper studies the thermal performance of the melting process in a packed bed system using spherical capsules through CFD modeling. The effective packed bed model was used to simulate the fluid flow in the packed bed and enthalpy porosity method was used to solve the phase change material (PCM) domain. The capsules were filled with three various types of phase change materials, while water was utilized as the heat transfer fluid (HTF). Due to the importance of configurations of materials in the packed beds, the effects of various layouts of them on the performance of the system were investigated. Moreover, Al2O3 nanoparticles as a cheap and easily available nanoparticles were added to the materials. Also, the change in inlet flow rate of the HTF was examined. According to the results, positioning the materials in decreasing order of melting temperature in the flow direction is able to decline the charging time of the system by about 1600 s and increase the charging efficiency by 6.8 %. The addition of Al2O3 nanoparticles was able to decrease the charging time by about 1900 s in the best case and increase the system's efficiency by about 7.5 %. Increasing the inlet flow rate of HTF from 2.5 × 10−5 to 5 × 10−5 (1.5–3 Lit/min) reduced the charging time about 1800 s.

Evaluation of wavy wall configurations for accelerated heat recovery in triplex-tube energy storage units for building heating applications

Efficient thermal energy storage (TES) is crucial for integrating intermittent renewable energy sources and managing fluctuations in energy supply and demand. Among TES methods, latent heat TES (LHTES) using phase change materials (PCMs) is highly promising due to its high energy storage density and nearly isothermal operation during phase transitions. However, the inherently low thermal conductivity of PCMs hinders heat transfer rates, negatively impacting charging and discharging performance. This study investigates novel wavy channel geometries for the PCM container in a triplex-tube LHTES unit to enhance heat transfer and improve solidification rates during the discharge process. A comprehensive computational model based on the enthalpy-porosity method is developed to simulate the PCM (paraffin wax RT35) solidification inside wavy channels with sinusoidal, zigzag, and step-function profiles. The effects of channel geometry, heat transfer fluid (HTF) velocity, and inlet temperature on solidification rate and heat recovery are systematically evaluated. Results demonstrate that the step-function geometry significantly accelerates solidification compared to straight channels, reducing discharge time by 65.1 % and increasing heat recovery rate by 147.9 %. Decreasing waviness width from 15 mm to 5 mm further reduces solidification time by 70.1 %, while increasing waviness height from 5 mm to 15 mm leads to a 74.4 % reduction. Moreover, increasing HTF flow from Re = 250 to 1000 enhances heat recovery by 92.3 %, and lowering inlet temperature from 20 °C to 10 °C improves it by 76.9 %. The findings provide valuable insights for designing efficient LHTES units with wavy channel geometries for renewable energy integration and thermal management applications.

Enhancing cold energy harvest with daytime passive radiative cooling and phase change materials integration: A numerical exploration

Daytime passive radiative cooling (DPRC) is an emerging technology to dissipate heat to the ultimate cold sink, the universe, with zero energy input. However, the achieved cold surface temperature of an object by DPRC can readily increase close to the ambient temperature. Sensible heat readily induces changes in the temperature of objects, but it generates a heat leakage path to surroundings with the largely varying temperature gradients. To address the challenge, we designed a DPRC coating integrated with phase change materials (PCM) so as to harvest and store the cold energy, respectively from the universe. To evaluate the feasibility of this concept, a numerical model, which integrates a self-programmed MATLAB code to calculate the net cooling power and a FLUENT-assisted enthalpy-porosity method to simulate the phase change process, was developed and validated. Also, the cold energy harvesting performance of the proposed device was explored under different climate types and conditions. We found that the surface temperature of the PCM with a DPRC coating can be reduced to −15.9 °C and −4.5 °C in a day in the temperate continental climate and the subtropical desert climate, respectively, when the non-radiative heat transfer is completely avoided. The temperature difference in the PCM is less than 1.0 °C when aluminum fins are added to enhance the heat transfer performance of the PCM. This study provides new insights into the continuous steady cooling of DPRC using PCMs, and the harvested cold energy can be used to cool water for the condenser of an air conditioner to save energy consumption.

Limitations of the enthalpy-porosity method for numerical simulation of close-contact melting on inclined surfaces

One of the main challenges for latent thermal energy storage (LTES) systems is low heat transfer rates due to the low thermal conductivity of most phase change materials (PCM). Close-contact melting (CCM) can accelerate melting times in LTES systems but the current numerical techniques for solid liquid phase change have difficulties with accurately predicting this process. In this study, close-contact melting of PCM on an inclined surface is simulated using the enthalpy-porosity method in ANSYS Fluent. All PCM properties, including density, are temperature-dependent. In this way, phenomena such as natural convection, volume change and buoyancy between the solid and liquid are taken into account. The volume change is compensated by a gaseous expansion volume. Both 2D and 3D simulations are used to show discrepancy between state-of-the-art enthalpy porosity modelling and experimentally observed phenomena in the case of CCM. The mushy zone constant, which is set to 105 to allow motion of the solid bulk, causes the solid phase to deform as a highly viscous fluid instead of moving as a rigid body. The velocity differences inside the solid are more than 50 % of its sinking velocity. As a result, the movement of the solid resembles creep behaviour and the obtained CCM patterns are not physically accurate. Furthermore, the density difference between the solid and liquid phases causes an avalanching effect in the mushy zone, which artificially strengthens convection. In conclusion, the enthalpy porosity method exhibits significant limitations in accurately capturing close-contact melting phenomena.

Numerical modelling of thermal hysteresis in melting and solidification of phase change materials

This research focuses on exploring the impact of thermal hysteresis effects on the performance of Phase Change Materials (PCMs) used in thermal storage or buffering applications. Computational Fluid Dynamics modeling in ANSYS Fluent is employed to investigate thermal hysteresis phenomena during solid-liquid phase change. The traditional enthalpy-porosity method is extended by introducing an alternative relationship between the liquid fraction and PCM temperature, implemented in Fluent through User Defined Functions. The study demonstrates the developed model by simulating a rectangular PCM enclosure with a heated wall featuring a time-dependent temperature boundary condition. Considering convective heat transfer, a comparison is made between the newly developed hysteresis model and the default solidification-melting model of ANSYS Fluent. The new model facilitates numerical analysis of thermal hysteresis during partial and complete melting and solidification processes, enabling the study of hysteresis effects on the part-load operation of PCM systems.

Dynamic Response of Phase Change Heat Exchange Unit with Layered Porous Media for Pulsed Electronic Equipment

Effective heat dissipation challenges transient high-power electronic devices in hypersonic vehicle cabins. This study introduces a Phase Change Heat Exchange Unit with Layered Porous Media (PCHEU–LPM) employing pulsed heat flow at the top and forced convection at the bottom. The primary aim was a comparative parametric study analyzing the thermal response of the heating surface under pulsed heat flow conditions. The geometric model was generated using electron microscopy images of manufactured objects and the numerical model was established based on the enthalpy–porosity method. Numerical simulations explored amplitude and frequency effects on pulsed thermal excitation, evaluating temperature and phase fields. A comprehensive time-frequency transformation assessed the temperature response. The results indicated an initial decrease and subsequent increase in interface temperature fluctuation with pulse heat flux amplitude growth. Temperature field uniformity correlated with natural convection strength in two-phase and liquid-phase regions. At mid and low frequencies, the phase change process increasingly suppressed interface temperature fluctuations. Optimal pulse thermal excitation selection was crucial for minimizing temperature fluctuations while maintaining the interface temperature within the expected phase transition range. In conclusion, a novel design concept is posited herein, aiming to enhance surface temperature uniformity and broaden the applicability of electronic devices through the manipulation of porosity rates.