The influence of the geometry of open-cell aluminum foam on the thermohydraulic behavior in channel flow is investigated. The mean cell diameter and the strut cross-sectional surface area are chosen as geometrical parameters, ranging respectively between 1.2 and 5.2mm and between 0.0125 and 0.17mm(2). The flow arrangement and the operating conditions are fixed. A numerical model is implemented in a commercial solver, based on volume averaging theory. The model is validated against experimental data. The porous properties, which take the sub-REV scaled physics into account, are written as a function of both geometrical parameters. The thermohydraulic characteristics of 16 well-chosen foams are used to build a surrogate model. An ordinary Kriging model is used for this, indicating that the root mean square error of interpolated results is lower than 0.6 and 3% for, respectively, heat transfer and total pressure. The resulting heat transfer and total pressure difference are nondimensionalized by dividing them by the results obtained from an empty channel. The relative increment of the pressure drop is an order of magnitude higher than the increment observed for heat transfer. Consequently, the applied performance evaluation criterion (defined as the ratio of dimensionless heat transfer versus total pressure) is mainly influenced by the hydraulic performance. For the given application, a clear optimum is found. The proposed method allows performing the parameter study with acceptable computational cost with a sufficient level of detail from an engineering perspective.