The eddy heat-flux in rotating turbulent convection

¹ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482, Potsdam, Germany e-mail: gruediger@aip.de
² Institute for Solar-Terrestrial Physics, PO Box 4026, Irkutsk 664033, Russia

Received: 15 July 2004
Accepted: 18 October 2004

Abstract

The three components of the heat-flux vector are numerically computed for a stratified rotating turbulent convection using the NIRVANA code in a flat box. The latitudinal component  proves to be negative (positive) in the northern (southern) hemisphere so that the heat always flows towards the poles. As a surprise, the radial heat-flux F_r peaks at the equator rather than at the poles (Taylor numbers ). The same behavior is observed for the radial turbulence intensity which for free turbulence is also believed to peak at the poles (see Eq. ([see full text]) below). As we can show, however, the consequences of this unexpected result (also obtained by Käpylä et al. [CITE], A&A, 422, 793) for the theory of differential rotation are small as mainly the  is responsible to solve the “Taylor number puzzle”. In all our simulations the azimuthal component  proves to be negative so that the rotating turbulence produces a westwards directed azimuthal heat-flux which should be observable. Fluctuations with higher temperature are expected to be anticorrelated with their own angular velocity fluctuations. We find this rotation-induced result as understandable as the  is closely related to the radial Λ-effect which is known to be also negative in stratified and rapidly rotating convection zones.