Two-equation model computations of high-speed (ma=2.25, 7.2), turbulent boundary layers

Abstract

The objective of this research is to assess the performance of two popularReynolds-averaged Navier-Stokes (RANS) models, standard k-E and k-w, andto suggest modifications to improve model predictions for high-speed flows. Numerical simulations of turbulent ow past a at plate are performed at M1 = 2:25; 7:2.The results from these two Mach number cases are compared with Direct NumericalSimulation (DNS) results from Pirozzoli et al. (2004) and experimental results fromHorstman & Owen (1975). The effect of the Boussinesq coefficient (Cu) and turbulenttransport coefficients (sigmak; sigmaE; sigma; sigma*) on the boundary layer ow is examined. Further,the performance of a new model with realizability-based correction to Cu and corresponding modifications to sigma; sigma* is examined. The modification to Cu is based oncontrolling the ratio of production to dissipation of kinetic energy (P/E=1). The firstchoice of P/E = 1 ensures that there is no accumulation of kinetic energy in stagnation or free-stream regions of the ow. The second choice of P/E= 1:6 holds underthe assumption of a homogeneous shear ow. It is observed that the new model'sperformance is similar to that of the existing RANS models, which is expected for asimple ow over a at plate. Finally, the role of turbulent Prandtl number (Prt) intemperature and density predictions is established. The results indicate that the k-wmodel's performance is better compared to that of the standard k-E model for highMach number flows. A modification to Cu must be accompanied with correspondingchanges to sigmak; sigmaE; sigma; sigma* for an accurate log-layer prediction. The results also indicate that a Prt variation is required across the boundary layer for improved temperatureand density predictions in high-speed flows.