Macroscopic properties of inhomogeneous, Navier–Stokes turbulence are calculated using a fundamental theoretical approach. Starting from a set of solenoidal basis vectors suitable for describing a turbulence bounded by infinite, parallel free-slip planes, invoking a random-phase approximation, and using the eddy-damped quasinormal Markovian closure, various quantities of dynamical interest are calculated. Quantities that affect the evolution of both the momentum density and the kinetic energy density of the turbulently evolving fluid are studied. These quantities include the ensemble-averaged values of gradients of the triple-velocity correlation and the pressure–velocity correlation. When the energy spectrum in wave-number space lacks certain reflection symmetries, a nontrivial mean flow velocity will be seen to emerge out of the turbulent eddies. In turn, these mean flow velocities will affect the evolution of the turbulence. Of course, such physically interesting mean-field structures cannot arise in dynamics suffused by the assumption of statistical homogeneity.