We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors, corresponding to locating the computational domain at different latitudes ranging from the pole to the equator of the star. We perform three-dimensional numerical simulations in the same parameter regimes for comparison. The free parameters appearing in the closure models are calibrated by two fitting methods using simulation data. Unique determination of the closure parameters is possible only in the non-rotating case or when the system is placed at the pole. In the other cases the fit procedures yield somewhat differing results. The quality of the closure is tested by substituting the resulting coefficients back into the closure model and comparing with the simulation results. To eliminate the possibilities that the results obtained depend on the aspect ratio of the simulation domain or suffer from too small Rayleigh numbers we performed runs varying these parameters. The simulation data for the Reynolds stress and heat fluxes broadly agree with previous compressible simulations. The closure works fairly well with slow and fast rotation but its quality degrades for intermediate rotation rates. We find that the closure parameters depend not only on rotation rate but also on latitude. The weak dependence on Rayleigh number and the aspect ratio of the domain indicates that our results are generally valid

Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution. We vary the Reynolds (Re) and P\'eclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.

We determine the timescales associated with turbulent diffusion and isotropization in closure models using anisotropically forced and freely decaying turbulence simulations and to study the applicability of these models. We compare the results from anisotropically forced three-dimensional numerical simulations with the predictions of the closure models and obtain the turbulent timescales mentioned above as functions of the Reynolds number. In a second set of simulations, turning the forcing off enables us to study the validity of the closures in freely decaying turbulence. Both types of experiments suggest that the timescale of turbulent diffusion converges to a constant value at higher Reynolds numbers. Furthermore, the relative importance of isotropization is found to be about 2.5 times larger at higher Reynolds numbers than in the more viscous regime.

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

To study the Reynolds stresses which describe turbulent momentum transport from turbulence affected by large-scale shear and rotation. Three-dimensional numerical simulations are used to study turbulent transport under the influences of large-scale shear and rotation in homogeneous, isotropically forced turbulence. We study three cases: one with only shear, and two others where in addition to shear, rotation is present. These cases differ by the angle (0 or 90\degr) the rotation vector makes with respect to the z-direction. Two subsets of runs are performed with both values of \theta where either rotation or shear is kept constant. When only shear is present, the off-diagonal stress can be described by turbulent viscosity whereas if the system also rotates, nondiffusive contributions (\Lambda-effect) to the stress can arise. Comparison of the direct simulations are made with analytical results from a simple closure model. We find that the turbulent viscosity is of the order of the first order smoothing result in the parameter regime studied and that for sufficiently large Reynolds numbers the Strouhal number, describing the ratio of correlation to turnover times, is roughly 1.5. This is consistent with the closure model based on the minimal tau-approximation which produces a reasonable fit to the simulation data for similar Strouhal numbers. In the cases where rotation is present, separating the diffusive and nondiffusive components of the stress turns out to be challenging but taking the results at face value, we can obtain nondiffusive contributions of the order of 0.1 times the turbulent viscosity. We also find that the simple closure model is able to reproduce most of the qualitative features of the numerical results provided that the Strouhal number is of the order of unity.