A mathematical model capable of predicting the structure of turbulent, underexpanded jets is described. The model is based on solutions of the fluid flow equations obtained using a second-order accurate, finite-volume integration scheme coupled to an adaptive grid algorithm. Turbulence within these jets is modelled using a k-epsilon approach coupled to the compressible dissipation rate model of Sarkar et al. (1991a). Comparison of model predictions and experimental data, reported in the literature, on a number of moderately underexpanded jets demonstrate significant improvements over results derived using the standard k-epsilon approach, and the adequacy of the compressibility corrected turbulence model for predicting such jets.