A priori analysis of subgrid scale pressure and heat flux in high pressure mixing and reacting shear layers

Abstract

Direct Numerical Simulation (DNS) data on high pressure H2/O2 and H2/air flames using the compressible flow formulation, detailed kinetics, a real fluid equation of state, and generalised diffusion are analysed. The DNS is filtered over a range of filter widths to provide exact terms in the Large Eddy Simulation (LES) governing equations, including unclosed terms. The filtered pressure and the filtered heat flux vector are extensively compared with the pressure and the heat flux vector calculated as a function of the filtered primitive variables (i.e. the exact LES term is compared with its form available within an actual LES). The difference between these forms defines the subgrid pressure and the subgrid heat flux vector. The analyses are done both globally across the entire flame, as well as by conditionally averaging over specific regions of the flame; including regions of large subgrid kinetic energy, subgrid scalar dissipation, subgrid temperature variance, flame temperature, etc. In this work, although negligible for purely mixing cases, the gradient of the subgrid pressure is shown to be of the same order as, and larger than, the corresponding divergence of the turbulent subgrid stresses for reacting cases. This is despite the fact that all species behave essentially as ideal gases for this flame and holds true even when the ideal gas law is used to calculate the pressure. The ratio of the subgrid pressure gradient to the subgrid stress tensor divergence is shown to increase with increasing Reynolds number. Both the subgrid heat flux vector and its divergence are found to be substantially larger in reacting flows in comparison with mixing due to the associated larger temperature gradients. However, the divergence of the subgrid heat flux vector tends to be significantly smaller than other unclosed terms in the energy equation with decreasing significance with increasing Reynolds number.