No Arabic abstract
A delayed detached eddy simulation (DDES) of an overexpanded nozzle flow with shock-induced separation is carried out at a Reynolds number equal to 1.7 10^7. The flow unsteadiness, characterised by self-sustained shock oscillations, induces local unsteady loads on the nozzle wall as well as global off-axis forces. A clear physical understanding of the driving factors of the unsteadiness is still lacking. Under the current conditions, the nozzle operates in a highly-overexpanded regime and comprises a conical separation shock within the nozzle contour, merging into a Mach disk in the nozzle centre. Our current study focuses on the unsteady pressure signature on the nozzle wall, through the use of Fourier-based spectral analysis performed in time and in the azimuthal wavenumber space. The numerical data well agrees with the experimental measurements in terms of mean and fluctuating wall pressure statistics. The frequency spectra are characterised by the presence of a large bump in the low frequency range associated to a breathing motion of the shock system and a broad and high amplitude peak at high frequencies generated by the turbulent activity of the detached shear layer. Moreover, a distinct peak at an intermediate frequency (of the order 1000 Hz) is observed to persist in the wall-pressure spectra along the nozzle wall. The analysis of the pressure signals in the azimuthal wavenumber space indicates that this peak is clearly associated to the first (non-symmetrical) pressure mode and it is thus connected to the generation of side loads. Furthermore, it is found that the unsteady Mach disk is characterised by an intense vortex shedding activity and the interaction of these vortices with the second shock cell is a key factor in the sustainment of an aeroacoustic feedback loop within the nozzle.
Emulsions are omnipresent in the food industry, health care, and chemical synthesis. In this Letter the dynamics of meta-stable oil-water emulsions in highly turbulent ($10^{11}leqtext{Ta}leq 3times 10^{13}$) Taylor--Couette flow, far from equilibrium, is investigated. By varying the oil-in-water void fraction, catastrophic phase inversion between oil-in-water and water-in-oil emulsions can be triggered, changing the morphology, including droplet sizes, and rheological properties of the mixture, dramatically. The manifestation of these different states is exemplified by combining global torque measurements and local in-situ laser induced fluorescence (LIF) microscopy imaging. Despite the turbulent state of the flow and the dynamic equilibrium of the oil-water mixture, the global torque response of the system is found to be as if the fluid were Newtonian, and the effective viscosity of the mixture was found to be several times bigger or smaller than either of its constituents.
A new approach to turbulence simulation, based on a combination of large-eddy simulation (LES) for the whole flow and an array of non-space-filling quasi-direct numerical simulations (QDNS), which sample the response of near-wall turbulence to large-scale forcing, is proposed and evaluated. The technique overcomes some of the cost limitations of turbulence simulation, since the main flow is treated with a coarse-grid LES, with the equivalent of wall functions supplied by the near-wall sampled QDNS. Two cases are tested, at friction Reynolds number Re$_tau$=4200 and 20,000. The total grid node count for the first case is less than half a million and less than two million for the second case, with the calculations only requiring a desktop computer. A good agreement with published DNS is found at Re$_tau$=4200, both in terms of the mean velocity profile and the streamwise velocity fluctuation statistics, which correctly show a substantial increase in near-wall turbulence levels due to a modulation of near-wall streaks by large-scale structures. The trend continues at Re$_tau$=20,000, in agreement with experiment, which represents one of the major achievements of the new approach. A number of detailed aspects of the model, including numerical resolution, LES-QDNS coupling strategy and sub-grid model are explored. A low level of grid sensitivity is demonstrated for both the QDNS and LES aspects. Since the method does not assume a law of the wall, it can in principle be applied to flows that are out of equilibrium.
We numerically examine the mechanisms that describe the shock-boundary layer interactions in transonic flow past an oscillating wing section. At moderate and high angles of incidence but low amplitudes of oscillation, shock induced flow separation or shock-stall is observed accompanied by shock reversal. Even though the power input to the airfoil by the viscous forces is three orders of magnitude lower than that due to the pressure forces on the airfoil, the boundary layer manipulates the shock location and shock motion and redistributes the power input to the airfoil by the pressure forces. The shock motion is reversed relative to that in an inviscid flow as the boundary layer cannot sustain an adverse pressure gradient posed by the shock, causing the shock to move upstream leading to an early separation. The shock motion shows a phase difference with reference to the airfoil motion and is a function of the frequency of the oscillation. At low angles of incidence, and low amplitudes of oscillation, the boundary layer changes the profile presented to the external flow, leads to a slower expansion of the flow resulting in an early shock, and a diffused shock-foot caused by the boundary layer.
Using exact relations between velocity structure functions (Hill, Hill and Boratav, and Yakhot) and neglecting pressure contributions in a first approximation, we obtain a closed system and derive simple order-dependent rescaling relationships between longitudinal and transverse structure functions. By means of numerical data with turbulent Reynolds numbers ranging from $Re_lambda=320$ to $Re_lambda=730$, we establish a clear correspondence between their respective scaling range, while confirming that their scaling exponents do differ. This difference does not seem to depend on Reynolds number. Making use of the Mellin transform, we further map longitudinal to (rescaled) transverse probability density functions.
An essential ingredient of turbulent flows is the vortex stretching mechanism, which emanates from the non-linear interaction of vorticity and strain-rate tensor and leads to formation of extreme events. We analyze the statistical correlations between vorticity and strain rate by using a massive database generated from very well resolved direct numerical simulations of forced isotropic turbulence in periodic domains. The grid resolution is up to $12288^3$, and the Taylor-scale Reynolds number is in the range $140-1300$. In order to understand the formation and structure of extreme vorticity fluctuations, we obtain statistics conditioned on enstrophy (vorticity-squared). The magnitude of strain, as well as its eigenvalues, is approximately constant when conditioned on weak enstrophy; whereas they grow approximately as power laws for strong enstrophy, which become steeper with increasing $R_lambda$. We find that the well-known preferential alignment between vorticity and the intermediate eigenvector of strain tensor is even stronger for large enstrophy, whereas vorticity shows a tendency to be weakly orthogonal to the most extensive eigenvector (for large enstrophy). Yet the dominant contribution to the production of large enstrophy events arises from the most extensive eigendirection, the more so as $R_lambda$ increases. Nevertheless, the stretching in intense vorticity regions is significantly depleted, consistent with the kinematic properties of weakly-curved tubes in which they are organized. Further analysis reveals that intense enstrophy is primarily depleted via viscous diffusion, though viscous dissipation is also significant. Implications for modeling are nominally addressed as appropriate.