A pseudo-spectral method with volume penalisation for magnetohydrodynamic turbulence in confined domains