Abstract
The objective of this study is to extend the attention of the incompressible smoothed particle hydrodynamics method (ISPH) in the heat transfer field. The ISPH method for the natural convection heat transfer under the Boussinesq approximation in various environments: pure-fluid, nanofluid, and non-Darcy porous medium is introduced. We adopted the improved analytical method for calculating the kernel renormalization factor and its gradient based on a quintic kernel function for the wall boundary treatment in the ISPH method. The proposed method requires no dummy particle layer to meet the impermeability condition and makes the heat flux over the wall boundary easy to implement. We performed four different numerical simulations of natural convection in cavities with increasing complexity in modeling and implementation: the natural convection in a square cavity with constant differentially heated wall temperature, natural convection with the heat flux from the bottom wall for a wide range of Rayleigh numbers, natural convection in a non-Darcy porous cavity fully filled with nanofluid in different flow regimes, and natural convection in a partially layered porous cavity. The results showed excellent agreement with results from literatures and the in-house P1-P1 finite element method code.