A Reynolds stress model of turbulence is extended for simulations of turbulent boundary layers over porous walls at high Reynolds numbers using the concept of volume-averaged Navier-Stokes equations. The interface between bulk flow and the flow in porous materials is treated by employing suitable jump conditions for the fluid stresses. The methodology is validated by using DNS-data and new experimental data for flat plate flow and the flow around an airfoil.