Differential rotation is a common feature of main-sequence spectral F-type stars. In seeking to make contact with observations and to provide a self-consistent picture of how differential rotation is achieved in the interiors of these stars, we use the three-dimensional anelastic spherical harmonic (ASH) code to simulate global-scale turbulent flows in 1.2 and 1.3 M F-type stars at varying rotation rates. The simulations are carried out in spherical shells that encompass most of the convection zone and a portion of the stably stratified radiative zone below it, allowing us to explore the effects of overshooting convection. We examine the scaling of the mean flows and thermal state with rotation rate and mass and link these scalings to fundamental parameters of the simulations. Indeed, we find that the differential rotation becomes much stronger with more rapid rotation and larger mass, scaling as ΔΩ ∝ M 3.9 Ω 0.6 0 . Accompanying the growing differential rotation is a significant latitudinal temperature contrast, with amplitudes of 1000 K or higher in the most rapidly rotating cases. This contrast in turn scales with mass and rotation rate as ΔT ∝ M 6.4 Ω 1.6 0 . On the other hand, the meridional circulations become much weaker with more rapid rotation and with higher mass, with their kinetic energy decreasing as KE MC ∝ M −1.2 Ω −0.8 0 . Additionally, three of our simulations exhibit a global-scale shear instability within their stable regions that persists for the duration of the simulations. The flow structures associated with the instabilities have a direct coupling to and impact on the flows within the convection zone.