A review of what numerical simulations tell us about the internal rotation of the sun

The simulated solar differential rotation from two independent numerical modeling efforts agree with each other and with present solar observations. The models solve the nonlinear, three-dimensional, time-dependent, anelastic equations of motion for thermal convection in a stratified, rotating, spherical shell. The simulated angular velocity in the convection zone is constant on cylinders coaxial with the rotation axis, maximum at the equator and decreasing with depth. The latitudinal variation of this angular velocity at the surface is in agreement with Doppler measurements of the solar surface rotation rate. The radial variation through the convection zone is consistent with the analysis of the rotational frequency splitting of solar oscillations. 15 refs., 5 figs