Turbulent flowfields within a spatially periodic array were calculated numerically using a finite difference method with a low Reynolds number, two-equation model of turbulence. Exploiting periodic boundary conditions, only a one-structural unit was taken as a calculation domain to simulate a porous medium of regular arrangement in an infinite space. Extensive numerical calculations were carried out for a wide range of Reynolds numbers, to elucidate hydrodynamic behaviors of turbulent flow (post-Forchheimer flow) in porous media. The microscopic numerical results thus obtained at a pore scale ...