The gas kinetic temperature in the centres of starless, high-density cores is predicted to fall as low as 5-6 K. The aim of this study was to determine the kinetic temperature distribution in the low-mass prestellar core Oph D where previous observations suggest a very low central temperature. The densest part of the Oph D core was mapped in the NH3(1,1) and (2,2) inversion lines using the Very Large Array (VLA). The physical quantities were derived from the observed spectra by fitting the hyperfine structure of the lines, and subsequently the temperature distribution of Oph D was calculated using the standard rotational temperature techniques. A physical model of the cores was constructed, and the simulated spectra after radiative transfer calculations with a 3D Monte Carlo code were compared with the observed spectra. Temperature, density, and ammonia abundance of the core model were tuned until a satisfactory match with the observation was obtained. The high resolution of the interferometric data reveals that the southern part of Oph D comprises of two small cores. The gas kinetic temperatures, as derived from ammonia towards the centres of the southern and northern core are 7.4 and 8.9 K, respectively. The observed masses of the cores are only 0.2 M_Sun. Their potential collapse could lead to formation of brown dwarfs or low-mass stars.

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.

We present local numerical models of accretion disk turbulence driven by the magnetorotational instability with varying shear rate. The resulting turbulent stresses are compared with predictions of a closure model in which triple correlations are modelled in terms of quadratic correlations. This local model uses five nondimensional parameters to describe the properties of the flow. We attempt to determine these closure parameters for our simulations and find that the model does produce qualitatively correct behaviour. In addition, we present results concerning the shear rate dependency of the magnetic to kinetic energy ratio. We find both the turbulent stress ratio and the total stress to be strongly dependent on the shear rate.