Summary: | A method for shock capturing by adaptive filtering for use with high-resolution, high-order schemes for Large Eddy Simulations (LES) is presented. The LES method used in all the examples here employs the Explicit Filtering approach and the spatial derivatives are obtained with sixth-order, compact, finite differences. The adaptation is to drop the order of the explicit filter to two at gridpoints where a shock is detected, and to then increase the order from 2 to 10 in steps at successive gridpoints away from the shock. The method is found to be effective in a series of tests of common inviscid 1D and 2D problems of shock propagation and propagation of waves through shocks. As a prelude to LES, the 3D Taylor−Green problem for the inviscid and a finite viscosity case were simulated. An assessment of the overall performance of the method for LES was carried out by simulating an underexpanded round jet at a Reynolds number of 6.09 million, based in centerline velocity and diameter at nozzle exit plane. Very close quantitative agreement was found for the development of centerline mean pressure when compared to experiment. Simulations on several increasingly finer grids showed a monotonic extension of the computed part of the inertial range, with little change to low frequency content. Amplitudes and locations of large changes in pressure through several cells were captured accurately. A similar performance was observed for LES of an impinging jet containing normal and curved shocks.
|