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.

Numerical simulations of the magnetorotational instability (MRI) with zero initial net flux in a non-stratified isothermal cubic domain are used to demonstrate the importance of magnetic boundary conditions.In fully periodic systems the level of turbulence generated by the MRI strongly decreases as the magnetic Prandtl number (Pm), which is the ratio of kinematic viscosity and magnetic diffusion, is decreased. No MRI or dynamo action below Pm=1 is found, agreeing with earlier investigations. Using vertical field conditions, which allow the generation of a net toroidal flux and magnetic helicity fluxes out of the system, the MRI is found to be excited in the range 0.1 < Pm < 10, and that the saturation level is independent of Pm. In the vertical field runs strong mean-field dynamo develops and helps to sustain the MRI.

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect-conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime increases as Rm^(1/2). Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar dynamos.

Using three-dimensional convection simulations it is shown that a sinusoidal variation of horizontal shear leads to a kinematic \alpha effect with a similar sinusoidal variation. The effect exists even for weak stratification and arises owing to the inhomogeneity of turbulence and the presence of impenetrable vertical boundaries. This system produces large-scale magnetic fields that also show a sinusoidal variation in the streamwise direction. It is argued that earlier investigations overlooked these phenomena partly because of the use of horizontal averaging and also because measurements of \alpha using an imposed field combined with long time averages give erroneous results. It is demonstrated that in such cases the actual horizontally averaged mean field becomes non-uniform. The turbulent magnetic diffusion term resulting from such non-uniform fields can then no longer be neglected and begins to balance the \alpha effect.

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

In this paper we make an effort to understand the interaction of turbulence generated by the magnetorotational instability (MRI) with turbulence from other sources, such as supernova explosions (SNe) in galactic disks. First we perform a linear stability analysis (LSA) of non-ideal MRI to derive the limiting value of Ohmic diffusion that is needed to inhibit the growth of the instability for different types of rotation laws. With the help of a simple analytical expression derived under first-order smoothing approximation (FOSA), an estimate of the limiting turbulence level and hence the turbulent diffusion needed to damp the MRI is derived. Secondly, we perform numerical simulations in local cubes of isothermal nonstratified gas with external forcing of varying strength to see whether the linear result holds for more complex systems. Purely hydrodynamic calculations with forcing, rotation and shear are made for reference purposes, and as expected, non-zero Reynolds stresses are found. In the magnetohydrodynamic calculations, therefore, the total stresses generated are a sum of the forcing and MRI contributions. To separate these contributions, we perform reference runs with MRI-stable shear profiles (angular velocity increasing outwards), which suggest that the MRI-generated stresses indeed become strongly suppressed as function of the forcing. The Maxwell to Reynolds stress ratio is observed to decrease by an order of magnitude as the turbulence level due to external forcing exceeds the predicted limiting value, which we interpret as a sign of MRI suppression. Finally, we apply these results to estimate the limiting radius inside of which the SN activity can suppress the MRI, arriving at a value of 14 kpc.

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.

The active RS CVn star II Peg has been spectroscopically monitored for almost 18 years with the SOFIN spectrograph at NOT, La Palma, Spain. In this paper we present five new surface temperature maps of the object for the years 1999 (two maps), 2001 (one map) and 2002 (two maps).

The sunspots and other solar activity indicators tend to cluster on the surface of the Sun.These clusters very often occur at certain longitudes that persist in time.It is of general interest to find new and simple ways to characterize the observed distributions of different indicators and their behaviour in time. In the present work we use Greenwich sunspot data to evaluate statistical but not totally coherent stability of sunspot distribution along latitudes as well as longitudes. The aim was to obtain information on the longitudinal distribution of the underlying spot-generating mechanism rather than on the distribution and migration of sunspots or sunspot groups on the solar surface. Therefore only sunspot groups were included in the analysis, and only the time of their first appearance was used. We use simple nonparametric approach to reveal sunspot migration patterns and their persistency. Our analysis shows that regions where spots are generated tend to rotate differentially as the spots and spot groups themselves do. The activity areas, however, tend to break down relatively fast, during 7-15 solar rotations. This study provides a challenge for solar dynamo models, as our results are consistent with the presence of a non-axisymmetric spot-generating mechanism experiencing differential rotation (known as phase mixing in dynamo theory). The new nonparametric method introduced here, completely independent of the choice of the longitudinal distribution of sunspots, was found to be a very powerful tool for spatio-temporal analysis of surface features.

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.

(abridged) Aims: To study turbulent transport coefficients that describe the evolution of large-scale magnetic fields in turbulent convection. Methods: We use the test field method together with 3D numerical simulations of turbulent convection with shear and rotation to compute turbulent transport coefficients describing the evolution of large-scale magnetic fields in mean-field theory in the kinematic regime. 1D mean-field models are used with the derived turbulent transport coefficients to compare with direct simulations. Results: The alpha-effect increases monotonically as rotation increases. Turbulent diffusivity, eta_t, is proportional to the square of the turbulent vertical velocity. Whereas eta_t decreases approximately inversely proportional to the wavenumber of the field, the alpha-effect and turbulent pumping show a more complex behaviour. In the presence of shear and no rotation a small alpha-effect is induced which does not seem to show any consistent trend as a function of shear. If the shear is large enough, this small alpha is able to excite a dynamo in the mean-field model. The coefficient responsible for driving the shear-current effect shows several sign changes as a function of depth but is also able to contribute to dynamo action in the mean-field model. The growth rates in these cases are well below those in direct simulations suggesting that an incoherent alpha-shear dynamo may also act in them. If both rotation and shear are present, the alpha-effect is more pronounced. The combination of the shear-current and Omega x J-effects is also stronger than in the case of shear only, but subdominant to the alpha-shear dynamo. The results of direct simulations are consistent with mean-field models where all of these effects are taken into account without the need to invoke incoherent effects.

Observations show that magnetic fields in the interstellar medium (ISM) often do not respond to increases in gas density as would be naively expected for a frozen-in field. This may suggest that the magnetic field in the diffuse gas becomes detached from dense clouds as they form. We have investigated this possibility using theoretical estimates, a simple magneto-hydrodynamic model of a flow without mass conservation and numerical simulations of a thermally unstable flow. Our results show that significant magnetic flux can be shed from dense clouds as they form in the diffuse ISM, leaving behind a magnetically dominated diffuse gas.

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.

(abridged) Aims: Three-dimensional numerical simulations of penetrative compressible convection with uniform horizontal shear are used to study dynamo action and the generation of large-scale magnetic fields. Methods: We consider cases where the magnetic Reynolds number is either marginal or moderately supercritical with respect to small-scale dynamo action in the absence of shear and rotation. The effects of magnetic helicity fluxes are studied by comparing results for the magnetic field with open and closed boundaries. Results: Without shear no large-scale dynamos are found even if the ingredients necessary for the alpha-effect (rotation and stratification) are present in the system. When uniform horizontal shear is added, a large-scale magnetic field develops, provided the boundaries are open. In this case the mean magnetic field contains a significant fraction of the total field. For those runs where the magnetic Reynolds number is between 60 and 250, small-scale dynamo is expected to be excited, but the field distribution is found to be similar to cases with smaller magnetic Reynolds number where the small-scale dynamo is not excited. In the case of closed (perfectly conducting) boundaries, magnetic helicity fluxes are suppressed and no large-scale fields are found. Similarly, poor large-scale field development is seen when vertical shear is used in combination with periodic boundary conditions in the horizontal directions. If, however, open (normal-field) boundary conditions are used in the x-direction, a large-scale field develops. These results support the notion that shear not only helps to generate the field, but it also plays a crucial role in driving magnetic helicity fluxes out of the system along the isocontours of shear, thereby allowing efficient dynamo action.

The thermal instability with a piecewise power law cooling function is investigated using one- and three-dimensional simulations with periodic and shearing-periodic boundary conditions in the presence of constant thermal diffusion and kinematic viscosity coefficients. Consistent with earlier findings, the flow behavior depends on the average density, <rho>. When $\bra{\rho}$ is in the range 1-5 x 10^{-24} g cm^{-3} the system is unstable and segregates into cool and warm phases with temperatures of roughly 100 and 10^4 K, respectively. However, in all cases the resulting average pressure <p> is independent of <rho> and just a little above the minimum value. For a constant heating rate of 0.015 erg g^{-1} s^{-1}, the mean pressure is around 24 x 10^{-14} dyn (corresponding to p/k_B ~ K cm^{-3}). Cool patches tend to coalesce into bigger ones. In all cases investigated there is no sustained turbulence, which is in agreement with earlier results. Simulations in which turbulence is driven by a body force show that when rms velocities of between 10 and 30 km/s are obtained, the resulting dissipation rates rates are comparable to the thermal energy input rate. The resulting mean pressures are then about 30 x 10^{-14} dyn, corresponding to p/k_B ~ 2170 K cm^{-3}. This is comparable to the value expected for the Galaxy. Differential rotation tends to make the flow two-dimensional, that is, uniform in the streamwise direction, but this does not lead to instability.

We study convective overshooting by means of local 3D convection calculations. Using a mixing length model of the solar convection zone (CZ) as a guide, we determine the Coriolis number (Co), which is the inverse of the Rossby number, to be of the order of ten or larger at the base of the solar CZ. Therefore we perform convection calculations in the range Co = 0...10 and interpret the value of Co realised in the calculation to represent a depth in the solar CZ. In order to study the dependence on rotation, we compute the mixing length parameters alpha_T and alpha_u relating the temperature and velocity fluctuations, respectively, to the mean thermal stratification. We find that the mixing length parameters for the rapid rotation case, corresponding to the base of the solar CZ, are 3-5 times smaller than in the nonrotating case. Introducing such depth-dependent alpha into a solar structure model employing a non-local mixing length formalism results in overshooting which is approximately proportional to alpha at the base of the CZ. Although overshooting is reduced due to the reduced alpha, a discrepancy with helioseismology remains due to the steep transition to the radiative temperature gradient. In comparison to the mixing length models the transition at the base of the CZ is much gentler in the 3D models. It was suggested recently (Rempel 2004) that this discrepancy is due to the significantly larger (up to seven orders of magnitude) input energy flux in the 3D models in comparison to the Sun and solar models, and that the 3D calculations should be able to approach the mixing length regime if the input energy flux is decreased by a moderate amount. We present results from local convection calculations which support this conjecture.

It has recently been claimed that analysis of Greenwich sunspot data over 120 years reveals that sunspot activity clusters around two longitudes separated by 180 degrees (``active longitudes'') with clearly defined differential rotation during activity cycles.In the present work we extend this critical examination of methodology to the actual Greenwich sunspot data and also consider newly proposed methods of analysis claiming to confirm the original identification of active longitudes. Our analysis revealed that values obtained for the parameters of differential rotation are not stable across different methods of analysis proposed to track persistent active longitudes. Also, despite a very thorough search in parameter space, we were unable to reproduce results claiming to reveal the century-persistent active longitudes. We can therefore say that strong and well substantiated evidence for an essential and century-scale persistent nonaxisymmetry in the sunspot distribution does not exist.

(abridged) Results from kinematic solar dynamo models employing alpha-effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2-3 m/s. The rotation profile of the Sun as obtained from helioseismology is applied. We obtain an estimate of the ratio of the two induction effects, C_alpha/C_Omega \approx 10^-3, which we keep fixed in all models. We also include a one-cell meridional circulation pattern having a magnitude of 10-20 m/s near the surface and 1-2 m/s at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the alpha-effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30 degrees of latitude. Turbulent pumping of the mean fields is predominantly down- and equatorward. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B_r B_phi. Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the alpha-effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5-60 degrees). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar-type solutions. The role of the surface shear layer is found to be important only in the case where the alpha-effect has an appreciable magnitude near the surface.

Aims. The alpha- and gamma-effects, which are responsible for the generation and turbulent pumping of large scale magnetic fields, respectively, due to passive advection by convection are determined in the rapid rotation regime corresponding to the deep layers of the solar convection zone. Methods. A 3D rectangular local model is used for solving the full set of MHD equations in order to compute the electromotive force (emf), E = <u x b>, generated by the interaction of imposed weak gradient-free magnetic fields and turbulent convection with varying rotational influence and latitude. By expanding the emf in terms of the mean magnetic field, E_i = a_ij <B_j>, all nine components of a_ij are computed. The diagonal elements of a_ij describe the alpha-effect, whereas the off-diagonals represent magnetic pumping. The latter is essentially the advection of magnetic fields by means other than the underlying large-scale velocity field. Comparisons are made to analytical expressions of the coefficients derived under the first-order smoothing approximation (FOSA). Results. In the rapid rotation regime the latitudinal dependence of the alpha-components responsible for the generation of the azimuthal and radial fields does not exhibit a peak at the poles, as is the case for slow rotation, but at a latitude of about 30 degrees. The magnetic pumping is predominantly radially down- and latitudinally equatorward as in earlier studies. The numerical results compare surprisingly well with analytical expressions derived under first-order smoothing, although the present calculations are expected to lie near the limits of the validity range of FOSA.

Aims: The Strouhal number (St), which is a nondimensional measure of the correlation time, is determined from numerical models of convection. The Strouhal number arises in the mean-field theories of angular momentum transport and magnetic field generation, where its value determines the validity of certain widely used approximations, such as the first order smoothing (hereafter FOSA). More specifically, the relevant transport coefficients can be calculated by means of a cumulative series expansion if St < St_crit \approx 1. Methods: We define the Strouhal number as the ratio of the correlation and turnover times, which we determine separately, the former from the autocorrelation of velocity, and the latter by following test particles embedded in the flow. Results: We find that the Strouhal numbers are, generally, of the order of 0.1 to 0.4 which is close to the critical value above which deviations from FOSA become significant. Increasing the rotational influence tends to shorten both timescales in such a manner that St decreases. However, we do not find a clear trend as function of the Rayleigh number for the parameter range explored in the present study.

We study the mixing length concept in comparison to three-dimensional numerical calculations of convection with rotation. In a limited range, the velocity and temperature fluctuations are linearly proportional to the superadiabaticity, as predicted by the mixing length concept and in accordance with published results. The effects of rotation are investigated by varying the Coriolis number, Co = 2 Omega tau, from zero to roughly ten, and by calculating models at different latitudes. We find that \alpha decreases monotonically as a function of the Coriolis number. This can be explained by the decreased spatial scale of convection and the diminished efficiency of the convective energy transport, the latter of which leads to a large increase of the superadibaticity, \delta = \nabla - \nabla_ad as function of Co. Applying a decreased mixing length parameter in a solar model yields very small differences in comparison to the standard model within the convection zone. The main difference is the reduction of the overshooting depth, and thus the depth of the convection zone, when a non-local version of the mixing length concept is used. Reduction of \alpha by a factor of roughly 2.5 is sufficient to reconcile the difference between the model and helioseismic results. The numerical results indicate reduction of \alpha by this order of magnitude.

We determine the Strouhal number (hereafter St), which is essentially a nondimensional measure of the correlation time, from numerical calculations of convection. We use two independent methods to estimate St. Firstly, we apply the minimal tau-approximation (MTA) on the equation of the time derivative of the Reynolds stress. A relaxation time is obtained from which St can be estimated by normalising with a typical turnover time. Secondly, we calculate the correlation and turnover times separately, the former from the autocorrelation of velocity and the latter by following test particles embedded in the flow. We find that the Strouhal number is in general of the order of 0.1 to 1, i.e. rather large in comparison to the typical assumption in the mean-field theories that St << 1. However, there is a clear decreasing trend as function of the Rayleigh number and increasing rotation. Furthermore, for the present range of parameters the decrease of St does not show signs of saturation, indicating that in stellar convection zones, where the Rayleigh numbers are much larger, the Strouhal number may indeed be significantly smaller.

We study stellar convection using a local three-dimensional MHD model, with which we investigate the influence of rotation and large-scale magnetic fields on the turbulent momentum and heat transport. The former is studied by computing the Reynolds stresses, the latter by calculating the correlation of velocity and temperature fluctuations, both as functions of rotation and latitude. We find that the horisontal correlation, Q_(theta phi), capable of generating horisontal differential rotation, is mostly negative in the southern hemisphere for Coriolis numbers exceeding unity, corresponding to equatorward flux of angular momentum in accordance with solar observations. The radial component Q_(r phi) is negative for slow and intermediate rotation indicating inward transport of angular momentum, while for rapid rotation, the transport occurs outwards. Parametrisation in terms of the mean-field Lambda-effect shows qualitative agreement with the turbulence model of Kichatinov & R\"udiger (1993) for the horisontal part H \propto Q_(theta phi)/cos(theta), whereas for the vertical part, V \propto Q_(r phi)/sin(theta), agreement only for intermediate rotation exists. The Lambda-coefficients become suppressed in the limit of rapid rotation, this rotational quenching being stronger for the V component than for H. We find that the stresses are enhanced by the presence of the magnetic field for field strengths up to and above the equipartition value, without significant quenching. Concerning the turbulent heat transport, our calculations show that the transport in the radial direction is most efficient at the equatorial regions, obtains a minimum at midlatitudes, and shows a slight increase towards the poles. The latitudinal heat transport does not show a systematic trend as function of latitude or rotation.

"The purpose of numerical models is not numbers but insight." (Hamming) In the spirit of this adage, and of Don Cox's approach to scientific speaking, we discuss the questions that the latest generation of numerical models of the interstellar medium raise, at least for us. The energy source for the interstellar turbulence is still under discussion. We review the argument for supernovae dominating in star forming regions. Magnetorotational instability has been suggested as a way of coupling disk shear to the turbulent flow. Models make evident that the unstable wavelengths are very long compared to thermally unstable wavelengths, with implications for star formation in the outer galaxy and low surface brightness disks. The perennial question of the factors determining the hot gas filling factor in a SN-driven medium remains open, in particular because of the unexpectedly strong turbulent mixing at the boundaries of hot cavities seen in the models. The formation of molecular clouds in the turbulent flow is also poorly understood. Dense regions suitable for cloud formation clearly form even in the absence of self-gravity, although their ultimate evolution remains to be computed.

In this letter, we present numerical calculations made to investigate the possible resolution dependence of the Shakura & Sunyaev (1973) viscosity parameter alpha from local magnetohydrodynamic simulations of the magnetorotational instability (MRI). We find that the values of alpha do indeed depend significantly on the numerical resolution but also that when the highest resolutions attainable by the computational resources available are used, the growth of the alpha-parameter seems to saturate. The values of alpha are at most of the order of 10^(-3), which indicates that the sole presence of turbulence due to dynamo generated magnetic field in the disk is not enough to reproduce alphas of the order unity which could explain some observational results (e.g. Cannizzo 1993).

In this paper we describe a numerical method designed for modelling different kinds of astrophysical flows in three dimensions. Our method is a standard explicit finite difference method employing the local shearing-box technique. To model the features of astrophysical systems, which are usually compressible, magnetised and turbulent, it is desirable to have high spatial resolution and large domain size to model as many features as possible, on various scales, within a particular system. In addition, the time-scales involved are usually wide-ranging also requiring significant amounts of CPU time. These two limits (resolution and time-scales) enforce huge limits on computational capabilities. The model we have developed therefore uses parallel algorithms to increase the performance of standard serial methods. The aim of this paper is to report the numerical methods we use and the techniques invoked for parallelising the code. The justification of these methods is given by the extensive tests presented herein.