Air cooling systems are currently the most popular and least expensive solutions to maintain a safe temperature in electronic devices. Heat sinks have been widely used in this area, allowing for an increase in the effective heat transfer surface area. The main objective of this study was to optimise the shape of the heat sink geometric model using the Adjoint Solver technique. The optimised shape in the context of minimal temperature value behind the heat sink is proposed. The effect of radiation and trapezoidal fin shape on the maximum temperature in the cooling system is also investigated. Simulation studies were performed in Ansys Fluent software using the Reynolds—averaged Navier–Stokes technique. As a result of the simulation, it turned out that not taking into account the radiation leads to an overestimation of temperatures in the system—even by 14 ${}^{\circ }$C. It was found that as the angle and height of the fins increases, the temperature value behind the heat sink decreases and the heat source temperature increases. The best design in the context of minimal temperature value behind the heat sink from all analysed cases is obtained for heat sink with deformed fins according to iteration 14. The temperature reduction behind the heat sink by as much as 25 ${}^{\circ }$C, with minor changes in heat source temperature, has been achieved.

Changes in the energy sector, associated with the move away from fossil fuels, pose a challenge for appropriate thermal energy management in residential buildings. The important element to deal with the variability of renewable energy in thermal systems is latent heat thermal energy storage. Due to the low thermal conductivity of phase change materials, a number of techniques are proposed to enhance the heat transfer process. In this research, the global sensitivity of fin geometrical parameters on the melting and solidification times and energy efficiency of these processes was investigated. The computational model of the phase change was developed using the finite volume method with the enthalpy-porosity model and Boussinesq approximation. Numerical simulations were carried out according to the design of experiments technique. The multi-dimensional response surface was developed, and the multi-objective optimisation was done. The research shows that the melting process is most influenced by the position of the top fin ($\alpha$ angle) and the solidification process by the position of the bottom fin ($\gamma$ angle). The angle of the tree fin ($\beta$) has a different effect on both processes, with the energy efficiency decreasing during melting and increasing during solidification. Maximum values for the energy efficiencies of melting (${\eta }_m=0.973$) and solidification (${\eta }_s=0.988$) were obtained for $\alpha =18.2^{°}$, $\beta =89.0^{°}$, $L=10.7$ mm and $\gamma =21.0^{°}$.

With the increasing trend toward the miniaturization of electronic devices, the issue of heat dissipation becomes essential. The use of phase changes in a two-phase closed thermosyphon (TPCT) enables a significant reduction in the heat generated even at high temperatures. In this paper, we propose a modification of the evaporation–condensation model implemented in ANSYS Fluent. The modification was to manipulate the value of the mass transfer time relaxation parameter for evaporation and condensation. The developed model in the form of a UDF script allowed the introduction of additional source equations, and the obtained solution is compared with the results available in the literature. The variable value of the mass transfer time relaxation parameter during condensation r${}_c$ depending on the density of the liquid and vapour phase was taken into account in the calculations. However, compared to previous numerical studies, more accurate modelling of the phase change phenomenon of the medium in the thermosyphon was possible by adopting a mass transfer time relaxation parameter during evaporation r${}_e$ = 1. The assumption of ten-fold higher values resulted in overestimated temperature values in all sections of the thermosyphon. Hence, the coefficient r${}_e$ should be selected individually depending on the case under study. A too large value may cause difficulties in obtaining the convergence of solutions, which, in the case of numerical grids with many elements (especially three-dimensional), significantly increases the computation time.