Rarefied gas flows in divergent nozzles were calculated by the DSMC method based upon three temperature model of translational, vibrational and rotational temperatures. As a molecular model, VSS model was adopted, and for partitioning of molecular collision energy, Larsen-Borgnakke method was used. In the calculation, dissociative and associative reactions of N_2, 0_2, NO, N, and O were included. The effects of nozzle angle and wall temperature were discussed in detail. It was pointed out that some improvements were required, particularly for rate coefficients of chemical reactions and modeling of chemical reactions associated with ions and electrons.