  Angular momentum transport in accretion discs is often believed to be due to
magnetohydrodynamic turbulence mediated by the magnetorotational instability.
Despite an abundant literature on the MRI, the parameters governing the
saturation amplitude of the turbulence are poorly understood and the existence
of an asymptotic behavior in the Ohmic diffusion regime is not clearly
established. We investigate the properties of the turbulent state in the small
magnetic Prandtl number limit. Since this is extremely computationally
expensive, we also study the relevance and range of applicability of the most
common subgrid scale models for this problem. Unstratified shearing boxes
simulations are performed both in the compressible and incompressible limits,
with a resolution up to 800 cells per disc scale height. The latter constitutes
the largest resolution ever attained for a simulation of MRI turbulence. In the
presence of a mean magnetic field threading the domain, angular momentum
transport converges to a finite value in the small Pm limit. When the mean
vertical field amplitude is such that {\beta}, the ratio between the thermal
and magnetic pressure, equals 1000, we find {\alpha}~0.032 when Pm approaches
zero. In the case of a mean toroidal field for which {\beta}=100, we find
{\alpha}~0.018 in the same limit. Both implicit LES and Chollet-Lesieur closure
model reproduces these results for the {\alpha} parameter and the power
spectra. A reduction in computational cost of a factor at least 16 (and up to
256) is achieved when using such methods. MRI turbulence operates efficiently
in the small Pm limit provided there is a mean magnetic field. Implicit LES
offers a practical and efficient mean of investigation of this regime but
should be used with care, particularly in the case of a vertical field.
Chollet-Lesieur closure model is perfectly suited for simulations done with a
spectral code.
