Submitted to Experiments in Fluids Preprint Experimental investigation of the turbulent Schmidt number in supersonic film cooling with shock interaction Pascal Marquardt Michael Klaas Wolfgang Schröder · · 05.11.2019 Abstract The interaction of an impinging shock and a supersonic helium cooling film is investigated experimentally by high-speed particle-image velocimetry. A lam- inar helium jet is tangentially injected into a turbulent air freestream at a freestream Mach number Ma¥ = 2:45. The helium cooling film is injected at a Mach number Mai = 1:30 at a total temperature ratio T0;i=T0;¥ = 0:75. A deflection b = 8◦ gen- erates a shock that impinges upon the cooling film. A shock interaction case and a reference case without shock interaction are considered. The helium mass frac- tion fluctuations are measured and the turbulent mass flux as well as the turbulent Schmidt number are determined qualitatively. For comparison, large-eddy simulation (LES) results of a comparable flow configuration are used. The streamwise and wall- normal turbulent mass fluxes are in qualitative agreement with the LES data. The turbulent Schmidt number differs significantly from unity. Without shock interaction, the turbulent Schmidt number is in the range 0:5 Sc 1:5 which is in agreement ≤ t ≤ with the literature. With shock interaction, the turbulent Schmidt number varies dras- Pascal Marquardt Chair of Fluid Mechanics and Institute of Aerodynamics Aachen Wüllnerstr. 5a 52062 Aachen Germany Tel.: +49-241-8090420 Fax: +49-241-8092257 E-mail: [email protected] Michael Klaas Chair of Fluid Mechanics and Institute of Aerodynamics Aachen Wüllnerstr. 5a 52062 Aachen Germany Wolfgang Schröder Chair of Fluid Mechanics and Institute of Aerodynamics Aachen Wüllnerstr. 5a 52062 Aachen Germany arXiv:1911.01780v1 [physics.flu-dyn] 5 Nov 2019 2 Pascal Marquardt et al. tically in the vicinity of the shock interaction. Thus, the experimental results confirm the numerical data showing a massively varying turbulent Schmidt number in super- sonic film cooling flows, i.e., the standard assumption of a constant turbulent Schmidt number is valid neither without nor with shock interaction. Keywords particle-image velocimetry supersonic film cooling helium injection · · · turbulent Schmidt number 1 Introduction In supersonic applications with high thermal loads, e.g., scramjet combustors, where the gas temperature exceeds the thermal limit of the surface material, so-called film cooling, i.e., the injection of cold gas along the surface, is often used as a cooling method. As shown in the review of Goldstein (1971), a promising cooling concept for supersonic flows is a cooling configuration with tangential injection. A schematic drawing of this cooling concept beneath a turbulent boundary layer is shown in fig- ure 1. According to Seban and Back (1962) and Juhany and Hunt (1994), the flow can be divided into three main regions. The first region is the potential core right down- stream of the injection, which is bounded by the mixing layer that emanates from the lip and the slot flow boundary layer. The potential core ends where the mixing layer and the slot-flow boundary layer merge. At this location, the wall-jet region, which is characterized by intense mixing, starts. Further downstream, the flow relaxes to an undisturbed turbulent boundary layer which is denoted as boundary-layer region. Shock waves that might be present in the flow field can cause shock induced separa- tion of the cooling film. If the onboard fuel, i.e., hydrogen, is considered as cooling fluid, the flow field exhibits a significant density gradient in the mixing layer. The combination of a free shear layer, a mixing layer, a wall bounded jet, a boundary layer, strong density and species gradients, and shock induced separation results in a complex flow field. To design efficient scramjet engines that consume as little coolant as possible, a precise prediction of the complex flow and temperature field and the distribution of coolant species is essential. This requires a profound understanding of the turbulent mixing mechanisms within the cooling-film flow, especially when flow separation occurs. Turbulent boundary layer Mixing layer Film-cooling flow Potential core Wall-jet region Boundary-layer region Fig. 1 Flow schematic with velocity profiles indicating three distinct flow regions (Seban and Back (1962); Juhany and Hunt (1994)) in a tangential film-cooling configuration (Konopka et al. (2012)). Exp. invest. of the turb. Schmidt number in supersonic film cooling with shock interaction 3 Cost efficient numerical methods to predict the flow in scramjet combustors, i.e., Reynolds-averaged Navier-Stokes (RANS) simulations, face the problem of model- ing the turbulent transport in the flow. Common models, e.g., the k-e model (Jones and Launder (1972)) or the k-w model (Wilcox (1988)), assume a constant relation- ship between the turbulent diffusion of momentum and the turbulent diffusion of heat and mass. The ratio between the momentum eddy diffusivity and the heat eddy dif- fusivity is denoted as turbulent Prandtl number Prt and the respective ratio for the mass eddy diffusivity is called the turbulent Schmidt number Sct. When simulations consider the temperature as a passive scalar, the turbulent Prandtl and the turbulent Schmidt number are identical. The present study is focused on the turbulent mass transport. Hence, the following review is limited to studies concerning the turbulent transport of mass or of a passive scalar. In general, a constant turbulent Schmidt number, which is at least problem dependent, is assumed in RANS simulations. In reviewing the optimal turbulent Schmidt numbers for engineering flow fields relevant to atmospheric dispersion, e.g., jet-in-cross flows and plume dispersion in boundary layers, Tominaga and Stathopoulos (2007) found a range for the optimal turbulent Schmidt number, i.e., the turbulent Schmidt number that showed the best agreement with experimental results, in the range 0:2 Sc 1:3. Therefore, depending on the ≤ t ≤ flow problem, a wide range of turbulent Schmidt numbers must be considered in RANS simulations. The strong effect of the turbulent Schmidt number on RANS of scramjet related flows was shown by Eklund et al. (2001), who performed RANS sim- ulations of an ethylene-fueled scramjet combustor with turbulent Schmidt numbers in the range 0:2 Sc 1:0. The turbulent Schmidt number strongly affects the location ≤ t ≤ of the pre-combustion shock train, the peak pressure, the ignition of the primary fuel, and the amount of heat release within the combustor. The turbulent Schmidt number not only depends on the flow problem, but also varies locally within the flow field. He et al. (1999) performed RANS simulations of a jet-in-cross flow at various turbulent Schmidt numbers and compared the results to experimental results of Crabb et al. (1981) and Kamotani and Greber (1974). Even though the authors found the best agreement between the numerical and the exper- imental results at a turbulent Schmidt number of Sct = 0:2, they concluded that a variable Schmidt number may be needed especially for low momentum flux ratio jet- in-cross flows. Later, Jiang and Campbell (2009) analyzed the temperature distribu- tion within a general combustor experimentally and numerically. They found an opti- mal agreement at a turbulent Schmidt number of Sct = 0:5. Nevertheless, the authors pointed out, that the concept of calculating turbulent scalar transfers based on the modeled momentum transfer and a turbulent Schmidt number is limited and should be improved and new approaches should be developed. Brinckman et al. (2007) pro- posed a scalar-variance model which predicts a variable turbulent Schmidt number. They found that the scalar-variance model showed better agreement with experimen- tal data over a range of reacting and non-reacting flows compared to computations for fixed Sct. Another eddy-viscosity-based turbulence model with variable Sct was proposed by Goldberg et al. (2010). For various benchmark flows including a scram- jet combustor flow the new turbulence model resulted in a better agreement with experimental data compared to simulations with fixed Sct. 4 Pascal Marquardt et al. In the past decade, numerical results which resolve all relevant flow scales, i.e., direct numerical simulations (DNS), became available for reasonably high Reynolds numbers. Li et al. (2009), Wu and Moin (2010), Araya and Castillo (2012), and Li et al. (2016) performed DNS of spatially developing incompressible boundary layers including a passive scalar up to Reynolds numbers Req = 2300. The simulations show variations of the turbulent Schmidt number within the boundary layer. The peak values close to the wall were in the range 1:1 Sc 1:9. Further off the wall, the ≤ t ≤ turbulent Schmidt number drops below unity. The experimental determination of the turbulent Schmidt number and the tur- bulent scalar fluxes is an extremely challenging task, especially in high-speed flows. Different measurement techniques have been used in low-speed water flow. Kooches- fahani et al. (2000) performed molecular tagging velocimetry (MTV) and laser in- duced fluorescence (LIF) measurements of a mixing layer between two water streams. They used a pulsed laser grid to excite a phosphorescent tracer compound and recorded the displacement of the excited molecules within a given time interval. Additionally, the laser excites fluorescence of fluorescein tracers in one of the streams to capture the tracer concentration. They demonstrated the ability to determine the