The highly turbulent nature of solar convection and anisotropies induced by rotation complicate the dynamics of modeling. Molecular dissipation scales on the Sun are at least six orders of magnitude smaller than the depth of the convective envelope. A direct numerical simulation of solar convection would have to resolve this entire range of scales in each of the three dimensions. Consequently, all solar differential rotation models must involve some approximations regarding momentum and heat transport by turbulent motions that are not explicitly computed. Thus, modeling approaches can be classified as either mean-field models or large-eddy simulations according to the approximations.