A calculation code has been developed in order to simulate a spatially developing compressible round jet by solving the axisym metric and three-dimensional compressible Navier-Stokes equations. A fourth-order non-oscillatory scheme and Roe's approximate Riemann solver are used in space discretization, and the third-order Runge-Kutta scheme is utilized as time integration. This code has an advantage that it can be easily extended to any order of accuracy in space and time. The inflow and outflow boundary conditions are important in simulating spatially developing subsonic flow. Though several boundary conditions based on the characteristic method or the wave equation have been proposed, they are still incomplete to get accurate solutions. Here in this study a non-reflecting condition developed by Thompson is employed because of its simplicity without empirical parameters.
The present code is applied to study early stages of spatially developing round jet for two Mach numbers of 0.5 and 1.5. The characteristics of compressible round jet, including the effect of compressibility, are examined. In order to verify the numerical scheme and boundary conditions, calculations are first performed for axisymmetric flow, and then they are extended to three-dimensional jet. In both axisymmetric and three-dimentional calculations, the axial location at which vortex ring rolls up moves downstream as the Mach number increases. In the axisymmetric flow, two vortex rings generated one after the other by Kelvin-Helmholtz instability merge into a larger vortex ring, which is repeated in the downstream. In the three-dimensional flow, the subsonic jet develops three-dimensional rib structures after axisymmetric roll-up and pairing. On the other hand, in the supersonic jet the pairing is prevented by large axial vorticity, which enhances three-dimensionality of vortex. To the contrary, in both subsonic and supersonic jets, longitudinal vortices in the downstream slow down the growth rates of shear layer.
View full abstract