Observations suggest that the large-scale convective velocities obtained by solar convection simulations might be over-estimated (convective conundrum). One plausible solution to this could be the small-scale dynamo which cannot be fully resolved by global simulations. The small-scale Lorentz force suppresses the convective motions and also the turbulent mixing of entropy between upflows and downflows, leading to a large effective Prandtl number ($\Pr$). We explore this idea in three-dimensional global rotating convection simulations at different thermal conductivity ($\kappa$), i.e., at different $\Pr$. In agreement with previous non-rotating simulations, the convective velocity is reduced with the increase of $\Pr$ as long as the thermal conductive flux is negligible. A subadiabatic layer is formed near the base of the convection zone due to continuous deposition of low entropy plumes in low-$\kappa$ simulations. The most interesting result of our low-$\kappa$ simulations is that the convective motions are accompanied by a change in the convection structure that is increasingly influenced by small-scale plumes. These plumes tend to transport angular momentum radially inward and thus establish an anti-solar differential rotation, in striking contrast to the solar rotation profile. If such low diffusive plumes, driven by the radiative-surface cooling, are present in the Sun, then our results cast doubt on the idea that a high effective $\Pr$ may be a viable solution to the solar convective conundrum. Our study also emphasizes that any resolution of the conundrum that relies on the downward plumes must take into account angular momentum transport as well as heat transport.

It has recently been recognized that the convective velocities achieved in the current solar convection simulations might be over-estimated. The newly-revealed effects of the prevailing small-scale magnetic field within the convection zone may offer possible solutions to this problem. The small-scale magnetic fields can reduce the convective amplitude of small-scale motions through the Lorentz-force feedback, which concurrently inhibits the turbulent mixing of entropy between upflows and downflows. As a result, the effective Prandtl number may exceed unity inside the solar convection zone. In this paper, we propose and numerically confirm a possible suppression mechanism of convective velocity in the effectively high-Prandtl number regime. If the effective horizontal thermal diffusivity decreases (the Prandtl number accordingly increases), the subadiabatic layer which is formed near the base of the convection zone by continuous depositions of low entropy transported by adiabatically downflowing plumes is enhanced and extended. The global convective amplitude in the high-Prandtl thermal convection is thus reduced especially in the lower part of the convection zone via the change in the mean entropy profile which becomes more subadiabatic near the base and less superadiabatic in the bulk.

The main object of the paper is to present the condition of the non-diffusive part of the Reynolds stress for driving the double-cell structure of the solar meridional circulation, which has been revealed by recent helioseismic observations. By conducting a set of mean-field hydrodynamic simulations, we confirm for the first time that the double-cell meridional circulation can be achieved along with the solar-like differential rotation when the Reynolds stress transports the angular momentum upward in the lower part and downward in the upper part of the convection zone. It is concluded that, in a stationary state, the accumulated angular momentum via the Reynolds stress in the middle layer is advected to both the upper and lower parts of the convection zone by each of the two meridional circulation cells, respectively.