The ground source heat pump heating system is considered as one of the best solutions for the transition towards green heating under the continental climate conditions like Kazakhstan. In this paper, experimental and theoretical investigations were carried out to develop a ground source heat pump-based heating system under the weather conditions in Kazakhstan and to evaluate its thermodynamic performance. The water-to-water heat pump heating system, integrated with a ground source heat exchanger and used refrigerant R134a, was designed to provide hot water to meet the requirements for space heating. The predicted values of the coefficient of performance and the experimental results were found to be in good agreement within 6.2%. The thermodynamic performance of the system was also assessed using various environment-friendly refrigerants, such as R152a, R450A, R513A, R1234yf and R1234ze, as potential replacements for R134a. Although R152a is found to be a good alternative for R134a in terms of coefficient of performance and total equivalent warming impact, its flammability hinders its application. The heating system using refrigerants R450A, R513A, R1234yf and R1234ze shows 2–3% lower coefficient of performance than that of R134a. The highest exergy destruction is found to be attributed to the compressor, followed by the expansion valve, evaporator, and condenser. Considering their low flammability and low environmental impact, R450A, R513A, R1234yf and R1234ze are identified as valuable replacements for R134a.

The purpose of the present work was to determine the thermal performance of borehole heat exchangers, considering the influences of their geometric configurations and the thermophysical properties of the soil, grout and pipe wall material. A three-dimensional model was developed for the heat and mass transfer in soil (a porous medium) and grout, together with one-dimensional conductive heat transfer through the pipe walls and one-dimensional convective heat transfer of the heat transfer fluid circulating in the pipes. An algorithm was developed to solve the mathematical equations of the model. The COMSOL Multiphysics software was used to implement the algorithm and perform the numerical simulations. An apparatus was designed, installed and tested to implement the thermal response test (TRT) method. Two wells of depth 50 m were drilled in the Almaty region in Kazakhstan. Gravel and till/loam were mainly found, which are in accordance with the stratigraphic map of the local geological data. In each well, two borehole heat exchangers were installed, which were an integral part of the ground source heat pump. The TRT measurements were conducted using one borehole heat exchanger in one well and the data were obtained. The present TRT data were found to be in good agreement with those available in literature. The numerical results of the model agreed well with the present TRT data, with the root-mean-square-deviation within 0.184 °C. The TRT data, together with the predictions of the line-source analytical model, were utilized to determine the soil thermal conductivity (${\lambda }_g$ = 2.35 W/m K) and the thermal resistance of the borehole heat exchanger from the heat transfer fluid to the soil ($R_b$ = 0.20 m K/W). The model was then used to predict the efficiencies of the borehole heat exchangers with various geometric configurations and dimensions. The simulation results show that the spiral borehole heat exchanger extracts the highest amount of heat, followed by the multi-tube, double U-type parallel, double U-type cross and single U-type. It is also found that the spiral configuration can save 34.6% drilling depth compared with the conventional single U-type one, suggesting that the spiral configuration is the best one in terms of the depth and the maximum heat extracted. The simulation results showed that (i) more heat was extracted with a higher thermal conductivity of grout material, in the range of 0.5–3.3 W/m K; (ii) the extracted heat remained unchanged for a thermal conductivity of pipe material higher than 2.0 W/m K (experiments in the range of 0.24–0.42 W/m K); (iii) the extracted heat remained unchanged for a volumetric flow rate of water higher than 1.0 m³/h (experimental flow rate 0.6 m³/h); and (iv) the heat extracted by the borehole heat exchanger increased with an increase in the thermal conductivity of the soil (experiments in the range of 0.4–6.0 W/m K). The numerical tool developed, the TRT data and simulation results obtained from the present work are of great value for design and optimization of borehole heat exchangers as well as studying other important factors such as the heat transfer performance during charging/discharging, freezing factor and thermal interference.

This paper presents a numerical analysis of two hot water storage tank configurations—one equipped with an external heat exchanger (Tank-1) and the other with an internal heat exchanger (Tank-2). The objective is to evaluate and compare their thermal performance during charging and discharging processes. The numerical model is developed by solving a system of ordinary differential equations using the 4th-order Runge–Kutta method, implemented in the Python programming language. The results indicate that Tank-1 demonstrated a higher charging efficiency of 94.6%, achieving full charge in approximately 2 h and 20 min. In comparison, Tank-2 required 3 h and 47 min to reach full charge, with a charging efficiency of 85.9%. During discharge, both configurations exhibited similar behavior, with an efficiency of 13.63% over approximately 33 min. The analysis showed that the external heat exchanger configuration led to more effective thermal stratification, supported by the Richardson number analysis, which indicated a significant effect of buoyancy during charging. This design advantage makes Tank-1 particularly suitable for applications requiring rapid heating and minimal heat loss, such as in cold climates or intermittent demand systems. The numerical model demonstrated reliable predictive accuracy, achieving an RMSE of 6.1% for the charging process and 6.8% for the discharging process, thereby validating the model’s reliability. These findings highlight the superior performance of the external heat exchanger configuration for fast and efficient energy storage, particularly for applications in cold climates.

Thermal stratification strongly affects the efficiency and operational reliability of sensible thermal energy storage (TES) tanks in energy systems. This study numerically investigates the combined influence of inlet configuration and mass flow rate on the charging performance of a vertical cylindrical TES tank (H = 3 m, D = 1 m) using transient CFD simulations. Five inlet designs—open, orifice, groove, shower, and shower-groove are analyzed at three flow rates: $Q_1$ = 0.0003 m³/s, ${Q_2=Q}_1/2$, and $Q_3=Q_1/3$. System performance is evaluated using key thermal and stratification metrics. Increasing the flow rate from $Q_3$ to $Q_1$ enhances convective heat transfer and energy and exergy efficiencies, but significantly intensifies mixing and degrades thermal stratification. At $Q_1$, the groove inlet achieves the highest capacity ratio and exergy efficiency (0.87), while exhibiting increased mixing. Reducing the flow rate to $Q_2$ and $Q_3$ limits inlet-induced momentum, leading to improved stratification for all configurations. The shower-groove inlet reaches a maximum stratification level (tail factor) of 1.13 at $Q_3$, indicating superior thermal layering, albeit with lower energetic efficiency (≈0.40–0.45). The groove inlet provides the best overall compromise at $Q_2$, combining high efficiency with stable stratification. These results demonstrate a clear efficiency-stratification trade-off and highlight the importance of selecting inlet-flow combinations according to application-specific objectives.

Short-term (daily) and long-term (seasonal) thermal energy storage allows efficient use of renewable thermal energy by replacing fossil fuel systems. In the present research, a three-dimensional numerical simulation of the thermal efficiency of a single-stage and three-stage cascaded shell-and-tube type latent heat thermal energy storage device is carried out using various phase change materials. The mathematical model is based on the fundamental conservation laws of mass, momentum, and energy. Numerical implementation was carried out using built-in solvers of COMSOL Multiphysics v.5.6 software. The numerical solver was verified by comparison with experimental data with acceptable agreement. A comparison of calculations between single and cascade configurations showed that the use of a multi-stage configuration allows for an increase in the temperature range of considered thermal energy storage devices. The use of the internal finned structure of containers almost halves the charging and discharging time. The developed numerical calculation tool can be used in the future to study the thermal efficiency of more complex thermal energy storage device configurations, considering three dimensions and phase transitions.