Three-dimensional numerical simulation of melting characteristics of phase change materials embedded with various TPMS skeletons

Abstract

Phase change energy storage technology has great potential for enhancing the efficient conversion and storage of energy. While triply periodic minimal surface (TPMS) structures have shown promise in improving heat transfer, research on their application in phase change heat transfer remains limited. This paper presents numerical simulations of composite phase change materials (PCMs) featuring TPMS skeletons, specifically gyroid, diamond, primitive, and I-graph and wrapped package-graph (I-WP) utilizing the lattice Boltzmann method (LBM). A comparative analysis of the effects of four TPMS skeletons on enhancing the phase change process reveals that the PCM containing the gyroid skeleton melts the fastest, with a complete melting time of 24.1% shorter than that of the PCM containing the I-WP skeleton. The PCM containing the gyroid skeleton is further simulated to explore the effects of the Rayleigh (Ra) number, Prandtl (Pr) number, and Stefan (Ste) number on the melting characteristics. Notably, the complete melting time is reduced by 60.44% when Ra is increased to 10⁶ compared to the case with Ra at 10⁴. Increasing the Pr number accelerates the migration of the mushy zone, resulting in fast melting. Conversely, the convective heat transfer effect from the heating surface decreases as the Ste number increases. The temperature differences caused by the local thermal non-equilibrium (LTNE) effect over time are significant and complex, with peaks becoming more pronounced nearer the heating surface. This study intends to provide theoretical support for the further development of TPMS skeletons in enhancing the phase change process.