A Statistical Theory of Thermally-Driven Turbulent Shear Flows, with the Derivation of a Subgrid Model