A new model for turbulent flows in porous media is developed. The spatial- and time fluctuations in this new model are tied together and treated as a single quantity. This novel treatment of the fluctuations leads to a natural construction of the kappa and epsilon type equations for rigid and isotropic porous media in which all the kinetic energy filtered in the averaging process is modeled. The same terms as those found in the corresponding equations for clear flow, plus additional terms resulting from the interaction between solid walls in the porous media and the fluid characterize the model. These extra terms arise in a boundary integral form, facilitating their modeling. The model is closed by assuming the eddy viscosity approximation to be valid, and using simple models to represent the interaction between the walls in the porous media and the fluid. (C) 2009 Elsevier Ltd. All rights reserved.