A numerically efficient configuration to simulate turbulent flows is to use triply periodic domains, with numerical forcing techniques to sustain turbulence. Previous homogeneous shear turbulence simulations considered only idealized homogeneous shear flows and not the statistically stationary shear turbulence observed in practical free shear flows. In contrast, the current study mathematically derives the complete forcing technique from the large scales of the turbulent free shear flows. Different statistically stationary free shear flows are considered in this study, namely, a nearly homogeneous shear turbulent flow, turbulent mixing layer, a turbulent planar jet, and a turbulent round jet. The simulations are performed on triply periodic, statistically homogeneous cubic domains in the vicinity of the shear layer in the self-similar region. An a priori analysis is performed to calculate the effects of the different forcing terms and to predict the expected turbulence quantities. The forcing technique is then used to perform direct numerical simulations at different Reynolds numbers. Numerical results for the different cases are discussed and compared with results from experiments and other simulations of free shear turbulent flows. Anisotropy is observed both in the components of velocity and vorticity, with stronger Reynolds number dependence in the anisotropy of vorticity. Energy spectra obtained from the present homogeneous shear turbulence agree well with the spectra from temporally evolving shear layers. The results also highlight the effects of the additional forcing terms that were neglected in previous studies and the role of shear convection and the associated splitting errors in the unbounded evolution of previous numerical simulations.