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.