In this paper, we apply a novel immersed boundary method to simulate pore-scale level fluid flow and convective heat transfer in realistic numerically generated open-cell solid foams in a Cartesian computational domain. Five different periodic foam samples of varying porosities (epsilon = [0.877, 0.948]) are generated by numerically mimicking the actual foam formation process (minimizing surface area). The step-by-step procedure for generating the periodic foam geometries is presented. The specific surface areas of the generated foams of different porosities are compared with real foam geometries showing a reasonable agreement. The Reynolds number (Re) is varied from Re approximate to 0 (creeping flow) to Re approximate to 500, and finally drag and Nusselt correlations have been proposed. A detailed analysis is presented on the local velocity and temperature field for the fluid-solid interaction in a complex cellular porous medium. (C) 2018 The Authors. Published by Elsevier Ltd.