Numerical Investigation of a Multi-Year Sand-Based Thermal Energy Storage System for Building Space Heating Application

Residential space heating in Northern Europe requires long-duration thermal storage to align summer solar gains with winter heating demand. This study investigates a compact sand-based seasonal thermal energy storage integrated with flat-plate solar collectors for an A+ class single-family house in Kaunas, Lithuania. An iterative co-design couples collector sizing with the seasonal charging target and a 3D COMSOL Multiphysics model of a 300 m3 sand-filled, phenolic foam-insulated system, with a 1D conjugate model of a copper pipe heat-exchanger network. The system was charged from March to September and discharged from October to February under measured-weather boundary conditions across three consecutive annual cycles. During the first year, the storage supplied the entire winter heating demand, though 35.2% of the input energy was lost through conduction, resulting in an end-of-cycle average sand temperature slightly below the initial state. In subsequent years, both the peak sand temperature and the residual end-of-cycle temperature increased by 3.7 °C and 3.2 °C, respectively, by the third year, indicating cumulative thermal recovery and improved retention. Meanwhile, the peak conductive losses rate decreased by 98 W, and cumulative annual losses decreased by 781.4 kWh in the third year, with an average annual reduction of 4.15%. These results highlight the progressive self-conditioning of the surrounding soil and demonstrate that a low-cost, sand-based storage system can sustain a complete seasonal heating supply with declining losses, offering a robust and scalable approach for residential building heating applications.