Numerical simulation of axisymmetric jets using differential eddy viscosity models
The work deals with numerical simulation using semi empirical turbulence models which include a differential equation for turbulent viscosity in two-dimensional case. k-ε-µt model [5] and k-ε-µt(Λ) model [6] are used to describe the mean turbulent parameters. The one-parameter models of eddy viscosity (νt -92 [7] and Spalart, Allmaras [8]) are used in one of the jet problems. Averaged flow and turbulence parameters are implemented in Cartesian and cylindrical coordinate systems. Simulation of axisymmetric subsonic jet in quiescent air in initial and transition regions with mean velocity u=87m/s have been conducted using k-ε-μt model. Initial region length has been determined with reasonable accuracy. Difference between calculated and experimental velocities is less than 5%. Verification of turbulence models in the flat plate boundary layer problem with M = 0.5 have been performed. Skin friction coefficient and velocity in logarithmic coordinates are compared with analytical solutions. Underexpanded supersonic turbulent jet simulation at M = 2 in the nozzle exit [13] have been conducted. Calculated static pressure and Mach number distributions along the jet axis are compared with experimental data. The k-ε-μt model reproduced 7 jet barrels. The positions of the maxima and minima are closest to the experiment in comparison with other models. The maxima and minima of the pressure oscillations and Mach in all barrels underestimated. The νt-92 model reproduced 12 jet barrels that is fully consistent with the experiment. The positions of the maxima and minima are distinctly shifted in comparison with experiment. The k-ε-µt(Λ) model reproduced the flow similarly with k-ε-μt model (7 barrels), but predicted flow core length is bigger. Spalart, Allmaras model reproduced only 5 barrels and predicted flow core length is smallest in comparison with other models. Numerical simulation of supersonic axisymmetric jet from converging nozzle in quiescent air [15] has been performed. A comparison with the Pitot pressure and velocity distribution along the axis is shown. The ratio of the chamber total pressure to ambient static pressure equals to NPR=2.5 and 4. In the case of k-ε-µt model at NPR=2.5 the initial region length is in satisfactory agreement with experiment, but the number of oscillations is underestimated. The k-ε-µt(Λ) model overestimates the initial region length. Oscillations of the Pitot pressure and velocity have bigger magnitude for all jet barrels in comparison with k-ε-µt model. In the case of k-ε-μt model at NPR=4 the Pitot pressure in the first barrel is bigger than experimental value. The difference in amplitude decreases downstream. The peaks begin to keep up and become less appreciable. Pitot pressure produced by k-ε-μt(Λ) model is underestimated in the first barrel. The subsequent barrels predicted by the model are much better than the first. Numerical simulation of overexpanded cold jet using the k-ε-μt model in quiescent air at Mach 3.005 in the nozzle exit [16] has been performed. The calculated axis Pitot pressure is underestimated, but the position of the local maxima and minima and local minima values are determined with good accuracy. The biggest difference along the axis between the calculation and experiment is in the two initial jet barrels. The largest deviation from the experiment in all cross-sections is found near the axis.