MHD simulation