In order to ensure the overall energy balance and the diffusive heat flux continuity across the solid-fluid interface,
a simplified method to deal the boundary layer treatment, namely the harmonic mean of thermal conductivities have
been employed immediately adjacent to the solid surface (about 1 mm thickness) to evaluate the diffusive heat flux
at the interface [15]. Then the diffusive heat flux across the fluid-solid interface is given by,