A direct numerical simulation of an oblique shock wave impinging on a turbulent boundary layer at Mach number 2.28 is carried out at moderate Reynolds number, simulating flow conditions similar to those of the experiment by Dupont et al. (J. Fluid Mech., vol. 559, 2006, pp. 255–277). The low-frequency shock unsteadiness, whose characteristics have been the focus of considerable research efforts, is here investigated via the Morlet wavelet transform. Owing to its compact support in both physical and Fourier spaces, the wavelet transformation makes it possible to track the time evolution of the various scales of the wall-pressure fluctuations. This property also makes it possible to define a local intermittency measure, representing a frequency-dependent flatness factor, to pinpoint the bursts of energy that characterise the shock intermittency scale by scale. As a major result, wavelet decomposition shows that the broadband shock movement is actually the result of a collection of sparse events in time, each characterised by its own temporal scale. This feature is hidden by the classical Fourier analysis, which can only show the time-averaged behaviour. Then, we propose a procedure to process any relevant time series, such as the time history of the wall pressure or that of the separation bubble extent, in which we use a condition based on the local intermittency measure to filter out the turbulent content in the proximity of the shock foot and to isolate only the intermittent component of the signal. In addition, wavelet analysis reveals the intermittent behaviour also of the breathing motion of the recirculation bubble behind the reflected shock, and allows us to detect a direct, partial correspondence between the most significant intermittent events of the separation region and those of the wall pressure at the foot of the shock.
Two-way momentum-coupled direct numerical simulations of a particle-laden turbulent channel flow are addressed to investigate the effect of the particle Stokes number and of the particle-to-fluid density ratio on the turbulence modification. The exact regularised point-particle method is used to model the interphase momentum exchange in presence of solid boundaries, allowing the exploration of an extensive region of the parameter space. Results show that the particles increase the friction drag in the parameter space region considered, namely the Stokes number $St_+ \in [2,80]$ , and the particle-to-fluid density ratio $\rho _p/\rho _f \in [90,5760]$ at a fixed mass loading $\phi =0.4$ . It is noteworthy that the highest drag occurs for small Stokes number particles. A measurable drag increase occurs for all particle-to-fluid density ratios, the effect being reduced significantly only at the highest value of $\rho _p/\rho _f$ . The modified stress budget and turbulent kinetic energy equation provide the rationale behind the observed behaviour. The particles’ extra stress causes an additional momentum flux towards the wall that modifies the structure of the buffer and of the viscous sublayer where the streamwise and wall-normal velocity fluctuations are increased. Indeed, in the viscous sublayer, additional turbulent kinetic energy is produced by the particles’ back-reaction, resulting in a strong augmentation of the spatial energy flux towards the wall where the energy is ultimately dissipated. This behaviour explains the increase of friction drag in particle-laden wall-bounded flows.
16 October 2024 in Journal of Parallel and Distributed Computing