A mathematical model is presented which conceptualises fluid flow and heat transfer in cellular metallic foams completely saturated with a fluid in motion. The model consists of a set of elliptic partial differential governing equations describing, firstly, a momentum balance in the fluid by the spatial distribution of its locally mean velocity, and secondly, an energy balance in the fluid and in the solid matrix of the metallic foam, by the spatial and temporal distribution of their locally mean temperatures. The separate energy balance descriptions for the fluid and the solid matrix extend the application of the model to conditions of thermal equilibrium and thermal non-equilibrium between the fluid and the solid matrix. A computational solution algorithm is presented which allows the universal application of the model to porous domains of arbitrary shape, with spatially and temporally variable heat loads in a variety of forms.