A number of Reynolds-averaged Navier-Stokes-based computational fluid dynamics solutions are presented that detail the effect of various thermochemical models on plume simulations. Four different levels of thermochemical models were considered, each with increasing complexity. In the simplest case, plumes were modeled with a single "perfect air" simulant gas that was calorically perfect. The second model used a single "equivalent" exhaust species that was calorically perfect but with molecular weight and frozen specific heats that mimicked the mass-averaged nozzle exit conditions. The third model used the same equivalent exhaust species but considered caloric imperfections. Finally, a full multispecies model with finite rate chemistry and general thermodynamics was considered. Comparisons are made between the different levels of modeling for power-on simulations involving the Ares I launch vehicle. The results show that, for simulations where plume impingement was very close to the exit of the nozzle or where insignificant plume impingement occurred, the perfect-air model performs well. For simulations in the vicinity of a cavity or where plume impingement occurred moderately downstream of the nozzle, the calorically imperfect equivalent model provides advantages over the other simplified thermodynamic models. Nomenclature A = nozzle cross-sectional area, m 2 C p = specific heat at constant pressure, J=Kg K C v = specific heat at constant volume, J=Kg K D = vehicle reference diameter, m e = internal energy, J=Kg M = Mach number M = molecular weight, Kg=Kg mol p = pressure, Pa R = gas constant, J=Kg K T = temperature, K u = velocity, m=s x = axial distance, m = ratio of specific heats = density, kg=m 3Subscripts eq = property for the equivalent species exit = property at the nozzle exit i = property for species "i" in = property at the nozzle inflow pg = property for a perfect gas rg = property for a real gas = property at the throat