A high-intensity swirling flow in a model combustor subjected to large density variations has been examined computationally. The focus is on the Favre-averaged Navier-Stokes computations of the momentum and scalar transport employing turbulence models based on the differential secondmoment closure (SMC) strategy. An updated version of the basic, high-Reynolds number SMC model accounting for a quadratic expansion of both the pressure-strain and dissipation tensors and a nearwall SMC model were used for predicting the mean velocity and turbulence fields. The accompanied mixing between the annular swirling air flow and the central non-swirling helium jet was studied by applying three scalar flux models differing mainly in the model formulation for the pressure-scalar gradient correlation. The computed axial and circumferential velocities agree fairly well with the reference experiment [So et al., NASA Contractor Report 3832, 1984; Ahmed and So, Exp. Fluids 4 (1986) 107], reproducing important features of such a weakly supercritical flow configuration (tendency of the flow core to separate). Although the length at which the mixing was completed was reproduced in reasonable agreement with the experimental results, the mixing activity in terms of the spreading rate of the shear/mixing layer, that is its thickness, was somewhat more intensive. Prior to these investigations, the model applied was validated by computing the transport of the passive scalar in the non-swirling (Johnson and Bennet, Report NASA CR-165574, UTRC Report R81-915540-9, 1981) and swirling (Roback and Johnson, NASA Contractor Report 168252, 1983) flow in a model combustor.