Hypersonic flows over blunt bodies subject to magnetic fields are numerically investigated. The magnetogasdynamic equations in the high magnetic Reynolds number formulation form an eight-equation system, with density, momentum, magnetic field, and total energy as unknowns. In the low magnetic Reynolds number approximation, the magnetic field induction is ignored, which leads to a five-equation system, where the magnetic interaction is represented by source terms in the momentum and energy equations. A four-stage modified Runge-Kutta scheme with the Davis-Yee symmetric total variation diminishing model as a postprocessing stage is used to solve the magnetogasdynamic equations. High-temperature effects are simulated by equilibrium and nonequilibrium chemistry models. The equilibrium model computes thermodynamic properties by interpolation from experimental data. The nonequilibrium model is a 1-temperature, 5-species, 17-reaction model solved by an implicit flux-vector splitting scheme. A loosely coupled approach is implemented to communicate between the magnetogasdynamic equations and the chemistry models