A state-specific vibrational model for CO2, CO, O2 and O system is devised by taking into account the first few vibrational states of each species. All vibrational states with energies at or below 1 eV are included in the present work. Of the three modes of vibration in CO2, the antisymmetric mode is considered separately from the symmetric stretching mode and the doubly degenerate bending modes. The symmetric and the bending modes are grouped together since the energy transfer rates between the two modes are very large due to Fermi resonance. The symmetric and bending modes are assumed to be in equilibrium with the translational and rotational modes. The kinetic rates for the vibrational-translation energy exchange reactions and the intermolecular and intramolecular vibrational-vibrational energy exchange reactions are based on experimental data to the maximum extent possible. Extrapolation methods are employed when necessary. This vibrational model is then coupled with an axisymmetric computational fluid dynamics code to study the expansion of CO2 in a nozzle.