The aim of this dissertation is to construct a virtual element method (VEM) for models in magneto-hydrodynamics (MHD), an area that studies the behavior and properties of electrically conducting fluids such as a plasma. MHD models are a coupling of the Maxwell’s equations for electromagnetics and models for fluid flow....