The implementation of renewable energy sources and heat pumps with natural refrigerants in the existing 3rd-generation district heating (DH) systems is a promising technology for the conversion to a 4th-generation DH system. This paper aims to investigate this transition through a case study for the existing DH system in Croatia. The district of Rijeka, which is considered in the case study, has an existing 3rd-generation DH system with a capacity of 9.2 MW, which was originally designed for a temperature regime of 130/70 °C and produces thermal energy from natural gas. In order to use heat pumps efficiently in such a system, the temperature of the distribution system and the energy consumption should be reduced. Trnsys software was used to perform complete dynamic simulations of the DH system. Used baseline models were validated according to the actual energy consumption reported by the heat supplier. The application of heat pumps with natural refrigerants in combination with on-site electricity production (PV system) can almost eliminate system CO₂ emissions, reduce primary energy consumption and lower operation costs. Simulations were performed for 18 different scenarios that included 2 different building envelope variants. The results were compared, then discussed and conclusions were drawn.

This study investigates the thermal performance of a passive vertical aluminum heat sink with plate fins through combined experimental measurements and numerical simulations. Using a custom-made experimental apparatus which used water as the heat source, heat transfer rate was determined, and heat transfer coefficient was compared against established empirical correlations, demonstrating good agreement. A 3D steady-state mathematical model was developed to capture the conjugate heat transfer problem of conduction and natural convection, with buoyancy-driven airflow modeled with the incompressible ideal gas law. The problem was solved numerically using the finite volume method through ANSYS Fluent 18.2 solver and validated against experimental data and analytical correlations, exhibiting good agreement throughout. Parametric analysis followed, investigating the influence of various base (50, 65, 80 °C) and ambient (19, 24, 29 °C) temperatures, resulting in base-to-ambient temperature differences from 21 to 61 °C. Increasing this temperature difference led to a significant increase in heat transfer rate, while heat transfer coefficient increased and overall thermal resistance decreased moderately. Additionally, a Nusselt–Rayleigh (Nu–Ra) number correlation, consistent with ranges reported in the literature, was derived, providing the scaling to predict the thermal performance of similar natural convection-governed heat sinks. The validated computational methodology, combined with obtained experimental and numerical results, presents a foundation for future studies focused on more complex heat sink geometries and physics.

In this paper, heat dissipation through a passive vertical plate fin heat sink via natural convection was numerically investigated. The influence of two nondimensional geometric parameters, fin spacing-to-thickness ratio and fin height-to-spacing ratio, on the heat sink’s thermal performance was evaluated. A mathematical model describing the three-dimensional steady-state problem of buoyancy-driven flow and heat transfer was formulated. The solution was obtained numerically using the finite volume method in Ansys Fluent. The model and numerical procedure were validated by comparing the numerical predictions with measurements acquired on a constructed experimental apparatus. The heat sink thermal performance was assessed based on a series of performance metrics: heat dissipation rate, heat transfer coefficient, overall thermal resistance, and fin efficiency. Fin spacing-to-thickness ratio was varied between 1.86 and 4.8. Heat dissipation rate displayed a clear peak at a value of approximately 2.6, which coincided with a minimum in the overall thermal resistance. The heat transfer coefficient initially grew steadily, but at higher values of fin spacing-to-thickness ratio, it began to stagnate. Fin efficiency consistently decreased across the investigated range. Fin height-to-spacing ratio was varied between 1.11 and 7.78. The heat dissipation rate increased almost linearly across this range, but when the mass-specific heat dissipation rate was considered, it yielded diminishing returns. The heat transfer coefficient likewise exhibited a plateauing trend, while fin efficiency decreased steadily across the investigated range of fin height-to-spacing ratio. The obtained numerical results provide guidelines for geometry selection and can serve as a foundation for further analyses and optimizations of passive heat sinks’ thermal performance.