Current state-of-the-art computer simulations allow us to build 3D dynamical and radiative models of stars from physical first principles with a high degree of realism. The radiative 3D dynamical stellar models obtained with the StellarBox code take into account the effects of turbulence, stellar abundances, a realistic equation of state, and radiative energy transport. In this talk, I will discuss the effects of rotation on the turbulent dynamics and surface structure for a 1.47 Msun star for rotational periods of 1 and 14 days. The simulations are performed with the computational domain at various latitudes. The models reproduce stellar granulation, the subsurface shear layer, structural changes in convection, and the tachocline, which is the interface between the inner radiative zone and the outer convection zone and plays a crucial role in stellar variability. In particular, the model results reveal the formation of differential rotation and meridional circulation.

The work is supported by NASA Astrophysics Theory Program.