No abstract

No abstract

We analyse the low-frequency dynamics of a high-Reynolds-number impinging shockwave/turbulent boundary-layer interaction (SWBLI) with strong mean-flow separation. The flow configuration for our grid-converged large-eddy simulations (LES) reproduces recent experiments for the interaction of a Mach 3 turbulent boundary layer with an impinging shock that nominally deflects the incoming flow by 19.6○ . The Reynolds number based on the incoming boundary layer thickness of Re δ0 ≈ 203 ⋅ 10 3 is considerably higher than in previous LES studies. The very long integration time of 3805 δ 0 U 0 allows for an accurate analysis of low-frequency unsteady effects. Experimental wallpressure measurements are in good agreement with the LES data. Both datasets exhibit the distinct plateau within the separated-flow region of a strong SWBLI. The filtered three-dimensional flow field shows clear evidence of counter-rotating streamwise vortices originating in the proximity of the bubble apex. Contrary to previous numerical results on compression ramp configurations, these Görtler-like vortices are not fixed at a specific spanwise position, but rather undergo a slow motion coupled to the separationbubble dynamics. Consistent with experimental data, power spectral densities (PSD) of wall-pressure probes exhibit a broadband and very energetic low-frequency component associated with the separation-shock unsteadiness. Sparsity-promoting dynamic mode decompositions (SPDMD) for both spanwise-averaged data and wall-plane snapshots yield a classical and well-known low-frequency breathing mode of the separation bubble, as well as a medium-frequency shedding mode responsible for reflected and reattachment shock corrugation. SPDMD of the two-dimensional skin-friction coefficient further identify streamwise streaks at low frequencies that cause large-scale flapping of the reattachment line. The PSD and SPDMD results of our impinging SWBLI support the theory that an intrinsic mechanism of the interaction zone is responsible for the lowfrequency unsteadiness, in which Görtler-like vortices might be seen as a continuous (coherent) forcing for strong SWBLI.

We present and evaluate a detailed multi-species two-phase thermodynamic equilibrium model for large-eddy simulations (LES) of liquid-fuel injection and mixing at high pressure. The model can represent the coexistence of supercritical states and multicomponent subcritical two-phase states. LES results for the transcritical Spray A of the Engine Combustion Network (ECN) are found to agree very well to available experimental data. We also address well-known numerical challenges of trans-and supercritical fluid mixing and compare a fully conservative formulation to a quasi conservative formulation of the governing equations. Our results prove physical and numerical consistency of both methods on fine grids and demonstrate the effects of energy conservation errors associated with the quasi conservative formulation on typical LES grids.

We investigate the shock-induced turbulent mixing between a light and a heavy gas, where a Richtmyer-Meshkov instability (RMI) is initiated by a shock wave with Mach number Ma = 1.5. The prescribed initial conditions define a deterministic multimode interface perturbation between the gases, which can be imposed exactly for different simulation codes and resolutions to allow for quantitative comparison. Wellresolved large-eddy simulations are performed using two different and independently developed numerical methods with the objective of assessing turbulence structures, prediction uncertainties and convergence behaviour. The two numerical methods differ fundamentally with respect to the employed subgrid-scale regularisation, each representing state-of-the-art approaches to RMI. Unlike previous studies, the focus of the present investigation is to quantify the uncertainties introduced by the numerical method, as there is strong evidence that subgrid-scale regularisation and truncation errors may have a significant effect on the linear and nonlinear stages of the RMI evolution. Fourier diagnostics reveal that the larger energy-containing scales converge rapidly with increasing mesh resolution and thus are in excellent agreement for the two numerical methods. Spectra of gradient-dependent quantities, such as enstrophy and scalar dissipation rate, show stronger dependences on the small-scale flow field structures as a consequence of truncation error effects, which for one numerical method are dominantly dissipative and for the other dominantly dispersive. Additionally, the study reveals details of various stages of RMI, as the flow transitions from large-scale nonlinear entrainment to fully developed turbulent mixing. The growth rates of the mixing zone widths as obtained by the two numerical methods are ∼t 7/12 before re-shock and ∼(t − t 0 ) 2/7 long after re-shock. The decay rate of turbulence kinetic energy is consistently ∼(t − t 0 ) −10/7 at late times, where the molecular mixing fraction approaches an asymptotic limit Θ ≈ 0.85. The anisotropy measure a xyz approaches an asymptotic limit of ≈0.04, implying that no full recovery of isotropy within the mixing zone is obtained, even after re-shock. Spectra of density, turbulence kinetic energy, scalar dissipation rate and enstrophy are presented and show excellent agreement for the resolved scales. The probability density function of the heavy-gas mass fraction and vorticity reveal that the light-heavy gas composition within the mixing zone is accurately predicted, whereas it is more difficult to capture the long-term behaviour of the vorticity.