This paper aims to introduce developments within the PHASTA flow solver that are needed
for applications in hypersonic CFD research. The PHASTA formulation is described and
unit tests are performed that serve to validate new features that enable the simulation of
thermochemical nonequilibrium. These features include a multi-species reacting flow module,
a two-temperature formulation that includes a model for vibrational temperature, a model
for mass diffusion, and the Spalart-Allmaras RANS turbulence model. The code accurately
predicts the physics for simple problems but requires further development for more complex
flows that represent real-world engineering problems.