We present a Fourier pseudo-spectral method for solving the resistive ă magnetohydrodynamic equations of incompressible flow in confined domains ă A volume penalisation method allows to take into account boundary ă conditions and the geometry of the domain A code validation is presented ă for the z-pinch test case Numerical simulations of decaying MHD ă turbulence in two space dimensions show spontaneous spin-up of the flow ă in non-axisymmetric geometries which is reflected by the generation of ă angular momentum First results of decaying MHD turbulence in a cylinder ă illustrate the potential of the new method for three-dimensional ă simulations (C) 2010 Elsevier BV All rights reserved