Viscous instability of a compressible round jet

The compressible linear stability equations are derived from the Navier Stokes equations in cylindrical polar coordinates. Numerical solutions for locally parallel flow are found using a direct matrix method. Discretization with compact finite difference is found to have better convergence properties than a Chebyshev spectral method for a round jet test case. The method is validated against previous results and convergence is tested for a range of jet profiles. Finally and en method is used to determine the dominant frequency of a Mach 0.9 jet

Surface correlations of hydrodynamic drag for transitionally rough engineering surfaces

Rough surfaces are usually characterised by a single equivalent sand-grain roughness height scale that typically needs to be determined from laboratory experiments. Recently, this method has been complemented by a direct numerical simulation approach, whereby representative surfaces can be scanned and the roughness effects computed over a range of Reynolds number. This development raises the prospect over the coming years of having enough data for different types of rough surfaces to be able to relate surface characteristics to roughness effects, such as the roughness function that quantifies the downward displacement of the logarithmic law of the wall. In the present contribution, we use simulation data for 17 irregular surfaces at the same friction Reynolds number, for which they are in the transitionally rough regime. All surfaces are scaled to the same physical roughness height. Mean streamwise velocity profiles show a wide range of roughness function values, while the velocity defect profiles show a good collapse. Profile peaks of the turbulent kinetic energy also vary depending on the surface. We then consider which surface properties are important and how new properties can be incorporated into an empirical model, the accuracy of which can then be tested. Optimised models with several roughness parameters are systematically developed for the roughness function and profile peak turbulent kinetic energy. In determining the roughness function, besides the known parameters of solidity (or frontal area ratio) and skewness, it is shown that the streamwise correlation length and the root-mean-square roughness height are also significant. The peak turbulent kinetic energy is determined by the skewness and root-mean-square roughness height, along with the mean forward-facing surface angle and spanwise effective slope. The results suggest feasibility of relating rough-wall flow properties (throughout the range from hydrodynamically smooth to fully rough) to surface parameters

A modification of Amiet's classical trailing edge noise theory for strictly two dimensional flows

The aim of this report is to derive theoretical expressions for the far-field pressure generated by disturbances convecting over a trailing edge. First, a general calculation of thefar-field pressure is discussed. Then the classical theory of Amiet (1976b) is reviewed,listing the most relevant assumptions. Amiet's theory is then revised for two-dimensional flows

Direct numerical simulation of turbulent flow over a rough surface based on a surface scan

Typical engineering rough surfaces show only limited resemblance to the artificially constructed rough surfaces that have been the basis of most previous fundamental research on turbulent flow over rough walls. In this article flow past an irregular rough surface is investigated, based on a scan of a rough graphite surface that serves as a typical example for an irregular rough surface found in engineering applications. The scanned map of surface height versus lateral coordinates is filtered in Fourier space to remove features on very small scales and to create a smoothly varying periodic representation of the surface. The surface is used as a no-slip boundary in direct numerical simulations of turbulent channel flow. For the resolution of the irregular boundary an iterative embedded boundary method is employed. The effects of the surface filtering on the turbulent flow are investigated by studying a series of surfaces with decreasing level of filtering. Mean flow, Reynolds stress and dispersive stress profiles show good agreement once a sufficiently large number of Fourier modes are retained. However, significant differences are observed if only the largest surface features are resolved. Strongly filtered surfaces give rise to a higher mean-flow velocity and to a higher variation of the streamwise velocity in the roughness layer compared with weakly filtered surfaces. In contrast, for the weakly filtered surfaces the mean flow is reversed over most of the lower part of the roughness sublayer and higher levels of dispersive shear stress are found

Sound radiation in turbulent channel flows

Lighthill’s acoustic analogy is formulated for turbulent channel flow with pressure as the acoustic variable, and integrated over the channel width to produce a two-dimensional inhomogeneous wave equation. The equivalent sources consist of a dipole distribution related to the sum of the viscous shear stresses on the two walls, together with monopole and quadrupole distributions related to the unsteady turbulent dissipation and Reynolds stresses respectively. Using a rigid-boundary Green function, an expression is found for the power spectrum of the far-field pressure radiated per unit channel area. Direct numerical simulations (DNS) of turbulent plane Poiseuille and Couette flow have been performed in large computational domains in order to obtain good resolution of the low-wavenumber source behaviour. Analysis of the DNS databases for all sound radiation sources shows that their wavenumber–frequency spectra have non-zero limits at low wavenumber. The sound power per unit channel area radiated by the dipole distribution is proportional to Mach number squared, while the monopole and quadrupole contributions are proportional to the fourth power of Mach number. Below a particular Mach number determined by the frequency and radiation direction, the dipole radiation due to the wall shear stress dominates the far field. The quadrupole takes over at Mach numbers above about 0.1, while the monopole is always the smallest term. The resultant acoustic field at any point in the channel consists of a statistically diffuse assembly of plane waves, with spectrum limited by damping to a value that is independent of Mach number in the low-M limit

Surface-sampled simulations of turbulent flow at high Reynolds number

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.Comment: Author accepted version. Accepted for publication in the International Journal for Numerical Methods in Fluids on 26 April 201

Instability and low-frequency unsteadiness in a shock-induced laminar separation bubble

Three-dimensional direct numerical simulations (DNS) of a shock-induced laminar separation bubble are carried out to investigate the flow instability and origin of any low frequency unsteadiness. A laminar boundary-layer interacting with an oblique shock-wave at M = 1:5 is forced at the inlet with a pair of monochromatic oblique unstable modes, selected according to local linear stability theory (LST) performed within the separation bubble. Linear stability analysis is applied to cases with marginal and large separation, and compared to DNS. While the parabolized stability equations approach accurately reproduces the growth of unstable modes, LST performs less well for strong interactions. When the modes predicted by LST are used to force the separated boundary-layer, transition to deterministic turbulence occurs near the reattachment point via an oblique-mode breakdown. Despite the clean upstream condition, broadband low-frequency unsteadiness is found near the separation point with a peak at a Strouhal number of 0:04, based on the separation bubble length. The appearance of the low-frequency unsteadiness is found to be due to the breakdown of the deterministic turbulence, filling up the spectrum and leading to broadband disturbances that travel upstream in the subsonic region of the boundary-layer, with a strong response near the separation point. The existence of the unsteadiness is supported by sensitivity studies on grid resolution and domain size that also identify the region of deterministic breakdown as the source of white noise disturbances. The present contribution confirms the presence of low-frequency response for laminar flows, similarly to that found in fully turbulent interactions

Large Eddy Simulations and modal reconstruction of laminar transonic buffet

Transonic buffet refers to the self-sustained periodic motion of shock waves observed in transonic flows over wings and limits the flight envelope of aircraft. Based on the boundary layer characteristics at the shock foot, buffet has been classified as laminar or turbulent and the mechanisms underlying the two have been proposed to be different (Dandois et al., 2018, J. Fluid Mech., vol. 18, pp. 156-178). The effect of various flow parameters (freestream Mach and Reynolds numbers and sweep and incidence angles) on laminar transonic buffet on an infinite wing (Dassault Aviation's supercritical V2C aerofoil) is reported here by performing Large-Eddy Simulations (LES) for a wide range of parameters. A spectral proper orthogonal decomposition identified the presence of a low-frequency mode associated with buffet and high-frequency wake modes related to vortex shedding. A flow reconstruction based only on the former shows periodic boundary-layer separation and reattachment accompanying shock wave motion. A modal reconstruction based only on the wake mode suggests that the separation bubble breathing phenomenon reported by Dandois et al. is due to this mode. Together, these results indicate that the physical mechanisms governing laminar and turbulent buffet are the same. Buffet was also simulated at zero incidence. Shock waves appear on both aerofoil surfaces and oscillate out of phase with each other indicating the occurrence of a Type I buffet (Giannelis et al., 2018, Aerosp. Sci. Technol., vol. 18, pp. 89-101) on a supercritical aerofoil. These results suggest that the mechanisms underlying different buffet types are the same.Comment: 40 pages, 31 figures, submitted to Journal of Fluid Mechanic