Using the formalism developed in paper I, we treat the case of shear‐driven flows.First, we derive the dynamic equations for the Reynolds stress. The equations are expressed in both tensorial and scalar forms, that is, as a set of coupled differential equations for the functions that enter the expansion of the Reynolds stress in terms of basic tensors. We specialize the general results to (a) axisymmetric contraction, (b) plane strain, and (c) homogeneous shear, for which there is a wealth of DNS, LES, and laboratory data to test the predictions of our model.Second, for homogeneous shear, in the inertial range, the equations for the Reynolds stress spectral function can be solved analytically,E12(k)=−C&egr;1/3Sk−7/3, which is in excellent agreement with recent data. Since the model has no free parameters, we stress that the model yields a numerical coefficientC, which is also in agreement with the data.Third, we derive the general expressions for the rapid and slow parts of the pressure–strain correlation tensors &Pgr;rijand &Pgr;sij. Within the second‐order closure models, the closure of &Pgr;sij(third‐order moment) in terms of second‐order moments continues to be particularly difficult. The general expression for &Pgr;ijare then specialized to the three flows discussed above. When &Pgr;sijis written in the form first suggested by Rotta, we show that the Rotta constant is a nonconstant tensor.Fourth, we discuss the dissipation tensor &egr;ij. In standard turbulence models, one not only assumes that &egr;ij=2/3&egr;&dgr;ij+f(uiuj), wheref(x) is a empirical function of the one‐point Reynolds stressuiuj, but, in addition, one employs a highly parametrized equation for &egr;. In the present model, neither of the two assumptions is required nor adjustable parameters are needed since &egr;ijis computed directly. The model provides thek‐dependentRij(k) as one of the primary quantities. ©1996 American Institute of Physics.