A multidimensional quasineutral plasma simulation model