JFM Papers
Transition in a separation bubble under tonal and broadband acoustic excitation
- John William Kurelek, Marios Kotsonis, Serhiy Yarusevych
-
- Published online by Cambridge University Press:
- 16 August 2018, pp. 1-36
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
Transition and flow development in a separation bubble formed on an airfoil are studied experimentally. The effects of tonal and broadband acoustic excitation are considered since such acoustic emissions commonly result from airfoil self-noise and can influence flow development via a feedback loop. This interdependence is decoupled, and the problem is studied in a controlled manner through the use of an external acoustic source. The flow field is assessed using time-resolved, two-component particle image velocimetry, the results of which show that, for equivalent energy input levels, tonal and broadband excitation can produce equivalent changes in the mean separation bubble topology. These changes in topology result from the influence of excitation on transition and the subsequent development of coherent structures in the bubble. Both tonal and broadband excitation lead to earlier shear layer roll-up and thus reduce the bubble size and advance mean reattachment upstream, while the shed vortices tend to persist farther downstream of mean reattachment in the case of tonal excitation. Through a cross-examination of linear stability theory (LST) predictions and measured disturbance characteristics, nonlinear disturbance interactions are shown to play a crucial role in the transition process, leading to significantly different disturbance development for the tonal and broadband excited flows. Specifically, tonal excitation results in transition being dominated by the excited mode, which grows in strong accordance with linear theory and damps the growth of all other disturbances. On the other hand, disturbance amplitudes are more moderate for the natural and broadband excited flows, and so all unstable disturbances initially grow in accordance with LST. For all cases, a rapid redistribution of perturbation energy to a broad range of frequencies follows, with the phenomenon occurring earliest for the broadband excitation case. By taking nonlinear effects into consideration, important ramifications are made clear in regards to comparing LST predictions and experimental or numerical results, thus explaining the trends reported in recent investigations. These findings offer new insights into the influence of tonal and broadband noise emissions, resulting from airfoil self-noise or otherwise, on transition and flow development within a separation bubble.
Cones of silence, complex rays and catastrophes: high-frequency flow–acoustic interaction effects
- J. T. Stone, R. H. Self, C. J. Howls
-
- Published online by Cambridge University Press:
- 16 August 2018, pp. 37-71
-
- Article
- Export citation
-
In this paper we develop a novel ray solver for the time-harmonic linearized Euler equations used to predict high-frequency flow–acoustic interaction effects from point sources in subsonic mean jet flows. The solver incorporates solutions to three generic ray problems found in free-space flows: the multiplicity of rays at a receiver point, propagation of complex rays and unphysical divergences at caustics. We show that these respective problems can be overcome by an appropriate boundary value reformulation of the nonlinear ray equations, a bifurcation-theory-inspired complex continuation, and an appeal to the uniform functions of catastrophe theory. The effectiveness of the solver is demonstrated for sources embedded in isothermal parallel and spreading jets, with the fields generated containing a wide variety of caustic structures. Solutions are presented across a large range of receiver angles in the far field, both downstream, where evanescent complex rays generate the cone of silence, and upstream, where multiple real rays are organized about a newly observed cusp caustic. The stability of the caustics is verified for both jets by their persistence under parametric changes of the flow and source. We show the continuation of these caustics as surfaces into the near field is complicated due to a dense caustic network, featuring a chain of locally hyperbolic umbilic caustics, generated by the tangency of rays as they are channelled upstream within the jet.
Insights into the dynamics of conical breakdown modes in coaxial swirling flow field
- Kuppuraj Rajamanickam, Saptarshi Basu
-
- Published online by Cambridge University Press:
- 22 August 2018, pp. 72-110
-
- Article
- Export citation
-
The main idea of this paper is to understand the fundamental vortex breakdown mechanisms in the coaxial swirling flow field. In particular, the interaction dynamics of the flow field is meticulously addressed with the help of high fidelity laser diagnostic tools. Time-resolved particle image velocimetry (PIV) (${\sim}1500~\text{frames}~\text{s}^{-1}$) is employed in $y{-}r$ and multiple $r{-}\unicode[STIX]{x1D703}$ planes to precisely delineate the flow dynamics. Experiments are carried out for three sets of co-annular flow Reynolds number $Re_{a}=4896$, 10 545, 17 546. Furthermore, for each $Re_{a}$ condition, the swirl number ‘$S_{G}$’ is varied independently from $0\leqslant S_{G}\leqslant 3$. The global evolution of flow field across various swirl numbers is presented using the time-averaged PIV data. Three distinct forms of vortex breakdown namely, pre-vortex breakdown (PVB), central toroidal recirculation zone (CTRZ; axisymmetric toroidal bubble type breakdown) and sudden conical breakdown are witnessed. Among these, the conical form of vortex breakdown is less explored in the literature. In this paper, much attention is therefore focused on exploring the governing mechanism of conical breakdown. It is should be interesting to note that, unlike other vortex breakdown modes, conical breakdown persists only for a very short band of $S_{G}$. For any small increase/decrease in $S_{G}$ beyond a certain threshold, the flow spontaneously reverts back to the CTRZ state. Energy ranked and frequency-resolved/ranked robust structure identification methods – proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) respectively – are implemented over instantaneous time-resolved PIV data sets to extract the dynamics of the coherent structures associated with each vortex breakdown mode. The dominant structures obtained from POD analysis suggest the dominance of the Kelvin–Helmholtz (KH) instability (axial $+$ azimuthal; accounts for ${\sim}80\,\%$ of total turbulent kinetic energy, TKE) for both PVB and CTRZ while the remaining energy is contributed by shedding modes. On the other hand, shedding modes contribute the majority of the TKE in conical breakdown. The frequency signatures quantified from POD temporal modes and DMD analysis reveal the occurrence of multiple dominant frequencies in the range of ${\sim}10{-}400~\text{Hz}$ with conical breakdown. This phenomenon may be a manifestation of high energy contribution by shedding eddies in the shear layer. Contrarily, with PVB and CTRZ, the dominant frequencies are observed in the range of ${\sim}20{-}40~\text{Hz}$ only. We have provided a detailed exposition of the mechanism through which conical breakdown occurs. In addition, the current work explores the hysteresis (path dependence) phenomena of conical breakdown as functions of the Reynolds and Rossby numbers. It has been observed that the conical mode is not reversible and highly dependent on the initial conditions.
Mean flow generation by an intermittently unstable boundary layer over a sloping wall
- Abouzar Ghasemi, Marten Klein, Andreas Will, Uwe Harlander
-
- Published online by Cambridge University Press:
- 22 August 2018, pp. 111-149
-
- Article
- Export citation
-
Direct numerical simulations (DNS) of the flow in various rotating annular confinements have been conducted to investigate the effects of wall inclination on secondary fluid motions due to an unstable boundary layer. The inner wall resembles a truncated cone (frustum) whose apex half-angle is varied from $18^{\circ }$ to $0^{\circ }$ (straight cylinder). The large inner radius $r_{1}$, the mean rotation rate $\unicode[STIX]{x1D6FA}_{0}$ and the kinematic viscosity $\unicode[STIX]{x1D708}$ were kept constant resulting in the constant Ekman number $E=\unicode[STIX]{x1D708}/(\unicode[STIX]{x1D6FA}_{0}r_{1}^{2})=4\times 10^{-5}$. Flows were excited by time-harmonic modulation of the inner wall’s rotation rate (so-called longitudinal libration) by prescribing the amplitude $\unicode[STIX]{x1D700}\unicode[STIX]{x1D6FA}_{0}$ and the forcing frequency $\unicode[STIX]{x1D714}=\unicode[STIX]{x1D6FA}_{0}$. By steepening the inner wall and hence reducing the effect of the local Coriolis force in the boundary layer three different flow regimes can be realized: a rotation-dominated, a libration-dominated and an intermediate regime. In the present study we focus on the libration-dominated regime. For small libration amplitudes (here $\unicode[STIX]{x1D700}=0.2$), a laminar Ekman–Stokes boundary layer (ESBL) is realized at the librating wall. With the aid of laminar boundary layer theory and DNS we show that the ESBL exhibits an oscillatory mass flux along the librating wall (Ekman property) and an oscillatory azimuthal velocity, which resembles a radially damped wave (Stokes property). For large libration amplitudes (here $\unicode[STIX]{x1D700}=0.8$), the DNS results exhibit an intermittently unstable ESBL, which turns centrifugally unstable during the prograde (faster) part of a libration period. This instability is due to the Stokes property and gives rise to Görtler vortices, which are found to be tilted with respect to the azimuth when the librating wall is at a finite angle relative to the axis of rotation. We show that this tilt is related to the Ekman property of the ESBL. This suggests that linear and nonlinear dynamics are equally important for this intermittent instability. Our DNS results indicate further that the Görtler vortices propagate into the fluid bulk where they generate an azimuthal mean flow. This mean flow is notably different from the mean flow driven in the case of the stable ESBL. A diagnostic analysis of the Reynolds-averaged Navier–Stokes (RANS) equations in the unstable flow regime hints at a competition between the radial and axial turbulent transport terms which act as generating and destructing agents for the azimuthal mean flow, respectively. We show that the balance of both terms depends on the wall inclination, that is, on the wall-tangential component of the Coriolis force.
On the scaling of propagation of periodically generated vortex rings
- H. Asadi, H. Asgharzadeh, I. Borazjani
-
- Published online by Cambridge University Press:
- 22 August 2018, pp. 150-170
-
- Article
- Export citation
-
The propagation of periodically generated vortex rings (period $T$) is numerically investigated by imposing pulsed jets of velocity $U_{jet}$ and duration $T_{s}$ (no flow between pulses) at the inlet of a cylinder of diameter $D$ exiting into a tank. Because of the step-like nature of pulsed jet waveforms, the average jet velocity during a cycle is $U_{ave}=U_{jet}T_{s}/T$. By using $U_{ave}$ in the definition of the Reynolds number ($Re=U_{ave}D/\unicode[STIX]{x1D708}$, $\unicode[STIX]{x1D708}$: kinematic viscosity of fluid) and non-dimensional period ($T^{\ast }=TU_{ave}/D=T_{s}U_{jet}/D$, i.e. equivalent to formation time), then based on the results, the vortex ring velocity $U_{v}/U_{jet}$ becomes approximately independent of the stroke ratio $T_{s}/T$. The results also show that $U_{v}/U_{jet}$ increases by reducing $Re$ or increasing $T^{\ast }$ (more sensitive to $T^{\ast }$) according to a power law of the form $U_{v}/U_{jet}=0.27T^{\ast 1.31Re^{-0.2}}$. An empirical relation, therefore, for the location of vortex ring core centres ($S$) over time ($t$) is proposed ($S/D=0.27T^{\ast 1+1.31Re^{-0.2}}t/T_{s}$), which collapses (scales) not only our results but also the results of experiments for non-periodic rings. This might be due to the fact that the quasi-steady vortex ring velocity was found to have a maximum of 15 % difference with the initial (isolated) one. Visualizing the rings during the periodic state shows that at low $T^{\ast }\leqslant 2$ and high $Re\geqslant 1400$ here, the stopping vortices become unstable and form hairpin vortices around the leading ones. However, by increasing $T^{\ast }$ or decreasing $Re$ the stopping vortices remain circular. Furthermore, rings with short $T^{\ast }=1$ show vortex pairing after approximately one period in the downstream, but higher $T^{\ast }\geqslant 2$ generates a train of vortices in the quasi-steady state.
Confinement effects on regular–irregular transition in shock-wave–boundary-layer interactions
- Ilan J. Grossman, Paul J. K. Bruce
-
- Published online by Cambridge University Press:
- 22 August 2018, pp. 171-204
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
An oblique shock wave is generated in a Mach 2 flow at a flow deflection angle of $12^{\circ }$. The resulting shock-wave–boundary-layer interaction (SWBLI) at the tunnel wall is observed. A novel traversable shock generator allows the position of the SWBLI to be varied relative to a downstream expansion fan. The relationship between the SWBLI, the expansion fan and the wind tunnel arrangement is studied. Schlieren photography, surface oil flow visualisation, particle image velocimetry and high-spatial-resolution wall pressure measurements are used to investigate the flow. It is observed that stream-normal movement of the shock generator downwards (towards the floor and hence the point of shock reflection) is accompanied by (1) growth in the streamwise extent of the shock-induced boundary layer separation, (2) upstream movement of the shock-induced separation point while the reattachment point remains nearly fixed, (3) an increase in separation shock strength and (4) transition between regular and irregular (Mach) reflection without an increase in incident shock strength. The role of free interaction theory in defining the separation shock angle is considered and shown to be consistent with the present measurements over a short streamwise extent. An SWBLI representation is proposed and reasoned which explains the apparent increase in separation shock strength that occurs without an increase in incident shock strength.
The linear instability of the stratified plane Couette flow
- Giulio Facchini, Benjamin Favier, Patrice Le Gal, Meng Wang, Michael Le Bars
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 205-234
-
- Article
- Export citation
-
We present the stability analysis of a plane Couette flow which is stably stratified in the vertical direction orthogonal to the horizontal shear. Interest in such a flow comes from geophysical and astrophysical applications where background shear and vertical stable stratification commonly coexist. We perform the linear stability analysis of the flow in a domain which is periodic in the streamwise and vertical directions and confined in the cross-stream direction. The stability diagram is constructed as a function of the Reynolds number $Re$ and the Froude number $Fr$, which compares the importance of shear and stratification. We find that the flow becomes unstable when shear and stratification are of the same order (i.e. $Fr\sim 1$) and above a moderate value of the Reynolds number $Re\gtrsim 700$. The instability results from a wave resonance mechanism already known in the context of channel flows – for instance, unstratified plane Couette flow in the shallow-water approximation. The result is confirmed by fully nonlinear direct numerical simulations and, to the best of our knowledge, constitutes the first evidence of linear instability in a vertically stratified plane Couette flow. We also report the study of a laboratory flow generated by a transparent belt entrained by two vertical cylinders and immersed in a tank filled with salty water, linearly stratified in density. We observe the emergence of a robust spatio-temporal pattern close to the threshold values of $Fr$ and $Re$ indicated by linear analysis, and explore the accessible part of the stability diagram. With the support of numerical simulations we conclude that the observed pattern is a signature of the same instability predicted by the linear theory, although slightly modified due to streamwise confinement.
Spatio-temporal intermittency of the turbulent energy cascade
- T. Yasuda, J. C. Vassilicos
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 235-252
-
- Article
- Export citation
-
In incompressible and periodic statistically stationary turbulence, exchanges of turbulent energy across scales and space are characterised by very intense and intermittent spatio-temporal fluctuations around zero of the time-derivative term, the spatial turbulent transport of fluctuating energy and the pressure–velocity term. These fluctuations are correlated with each other and with the intense intermittent fluctuations of the interscale energy transfer rate. These correlations are caused by the sweeping effect, the link between nonlinearity and non-locality, and also relate to geometrical alignments between the two-point fluctuating pressure force difference and the two-point fluctuating velocity difference in the case of the correlation between the interscale transfer rate and the pressure–velocity term. All these processes are absent from the spatio-temporal-average picture of the turbulence cascade in statistically stationary and homogeneous turbulence.
Approximate solutions to droplet dynamics in Hele-Shaw flows
- Yoav Green
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 253-270
-
- Article
- Export citation
-
For the past decade, the interaction force between droplets flowing in a Hele-Shaw cell has been modelled as a dipole. In this work, we use the recently derived analytical solution of Sarig et al. (J. Fluid Mech., vol. 800, 2016, pp. 264–277) of a two-droplet system, which satisfies the no-flux condition at both droplet interfaces, and compare it to results of the dipole model, which does not satisfy the no-flux condition. Unfortunately, the recently derived solution is given in terms of infinite Fourier series, making any additional straightforward analysis difficult. We derive simple approximations for these Fourier series. We show that at large spacing the approximations for the interactions reduce to the expected dipole-like solution. We also provide a new lower limit for the velocity for the case of almost touching droplets. For the case of large spacing, the derivation is extended to arbitrary droplet numbers – including an infinite lattice. We present a new correction for the dispersion relation for the perturbations. We investigate the effect of the number of droplets in a lattice, $N$, on the resulting dynamics.
Clarifying the relationship between efficiency and resonance for flexible inertial swimmers
- Daniel Floryan, Clarence W. Rowley
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 271-300
-
- Article
- Export citation
-
We study a linear inviscid model of a passively flexible swimmer, calculating its propulsive performance, eigenvalues and eigenfunctions with an eye towards clarifying the relationship between efficiency and resonance. The frequencies of actuation and stiffness ratios we consider span a large range, while the mass ratio is mostly fixed to a low value representative of swimmers. We present results showing how the trailing edge deflection, thrust coefficient, power coefficient and efficiency vary in the stiffness–frequency plane. The trailing edge deflection, thrust coefficient and power coefficient show sharp ridges of resonant behaviour for mid-to-high frequencies and stiffnesses, whereas the efficiency does not show resonant behaviour anywhere. For low frequencies and stiffnesses, the resonant peaks smear together and the efficiency is high. In this region, flutter modes emerge, inducing travelling wave kinematics which make the swimmer more efficient. We also consider the effects of a finite Reynolds number in the form of streamwise drag. The drag adds an offset to the net thrust produced by the swimmer, causing resonant peaks to appear in the efficiency (as observed in experiments in the literature).
Vibrations of a square cylinder submerged in a wake
- Rajesh Bhatt, Md. Mahbub Alam
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 301-332
-
- Article
- Export citation
-
A numerical investigation is conducted on the flow around and vibration response of an elastic square cylinder (side width $D$) in the wake of a stationary cylinder at Reynolds numbers of $Re=100$ and 200 based on $D$ and the free-stream velocity. The downstream cylinder, referred to as the wake cylinder, is allowed to vibrate in the transverse direction only. The reduced velocity $U_{r}$ is varied from 1 to 30. Cylinder centre-to-centre spacing ratios of $L^{\ast }(=L/D)=2$ and 6 are considered. Simulations are also conducted for a single isolated cylinder, and the results are compared with those for the wake cylinder. The focus is given to vibration response, frequency response, fluctuating lift force, phase relationship between the lift and displacement, work done and the flow structure modification during the cylinder vibration. The results reveal that the dependence of the Strouhal number $St$ on $U_{r}$ can distinguish different branches more appropriately than that of the vibration amplitude on $U_{r}$. The vibration response of the single cylinder at $Re=100$ is characterized by the initial, lower and desynchronization branches. On the other hand, that at $Re=200$ undergoes initial, lower and galloping branches. The galloping involves the characteristics of both the initial and the lower branches or the initial and the desynchronization branches depending on $U_{r}$. For the wake cylinder, the gap flow has a significant impact on the vibration response, leading to (i) the absence of galloping at either $Re$ and $L^{\ast }$, (ii) the presence of an upper branch at $Re=200$, $L^{\ast }=6$ and (iii) an initial branch of different characteristics at $Re=100$, $L^{\ast }=6$. The different facets are discussed in terms of wake structures, work done and phase lag between lift and displacement.
Jet–flap interaction tones
- Peter Jordan, Vincent Jaunet, Aaron Towne, André V. G. Cavalieri, Tim Colonius, Oliver Schmidt, Anurag Agarwal
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 333-358
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
Motivated by the problem of jet–flap interaction noise, we study the tonal dynamics that occurs when an isothermal turbulent jet grazes a sharp edge. We perform hydrodynamic and acoustic pressure measurements to characterise the tones as a function of Mach number and streamwise edge position. The observed distribution of spectral peaks cannot be explained using the usual edge-tone model, in which resonance is underpinned by coupling between downstream-travelling Kelvin–Helmholtz wavepackets and upstream-travelling sound waves. We show, rather, that the strongest tones are due to coupling between Kelvin–Helmholtz wavepackets and a family of trapped, upstream-travelling acoustic modes in the potential core, recently studied by Towne et al. (J. Fluid Mech. vol. 825, 2017) and Schmidt et al. (J. Fluid Mech. vol. 825, 2017). We also study the band-limited nature of the resonance, showing the high-frequency cutoff to be due to the frequency dependence of the upstream-travelling waves. Specifically, at high Mach number, these modes become evanescent above a certain frequency, whereas at low Mach number they become progressively trapped with increasing frequency, which inhibits their reflection in the nozzle plane.
Optimal mixing in two-dimensional stratified plane Poiseuille flow at finite Péclet and Richardson numbers
- F. Marcotte, C. P. Caulfield
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 359-385
-
- Article
- Export citation
-
We consider the nonlinear optimisation of irreversible mixing induced by an initial finite amplitude perturbation of a statically stable density-stratified fluid with kinematic viscosity $\unicode[STIX]{x1D708}$ and density diffusivity $\unicode[STIX]{x1D705}$. The initial diffusive error function density distribution varies continuously so that $\unicode[STIX]{x1D70C}\in [\bar{\unicode[STIX]{x1D70C}}-\unicode[STIX]{x1D70C}_{0}/2,\bar{\unicode[STIX]{x1D70C}}+\unicode[STIX]{x1D70C}_{0}/2]$. A constant pressure gradient is imposed in a plane two-dimensional channel of depth $2h$. We consider flows with a finite Péclet number $Pe=U_{m}h/\unicode[STIX]{x1D705}=500$ and Prandtl number $Pr=\unicode[STIX]{x1D708}/\unicode[STIX]{x1D705}=1$, and a range of bulk Richardson numbers $Ri_{b}=g\unicode[STIX]{x1D70C}_{0}h/(\bar{\unicode[STIX]{x1D70C}}U^{2})\in [0,1]$ where $U_{m}$ is the maximum flow speed of the laminar parallel flow, and $g$ is the gravitational acceleration. We use the constrained variational direct-adjoint-looping (DAL) method to solve two optimisation problems, extending the optimal mixing results of Foures et al. (J. Fluid Mech., vol. 748, 2014, pp. 241–277) to stratified flows, where the irreversible mixing of the active scalar density leads to a conversion of kinetic energy into potential energy. We identify initial perturbations of fixed finite kinetic energy which maximise the time-averaged perturbation kinetic energy developed over a finite time interval, and initial perturbations that minimise the value (at a target time, chosen to be $T=10$) of a ‘mix-norm’ as first introduced by Mathew et al. (Physica D, vol. 211, 2005, pp. 23–46), further discussed by Thiffeault (Nonlinearity, vol. 25, 2012, pp. 1–44) and shown by Foures et al. (2014) to be a computationally efficient and robust proxy for identifying perturbations that minimise the long-time variance of a scalar distribution. We demonstrate, for all bulk Richardson numbers considered, that the time-averaged kinetic-energy maximising perturbations are significantly suboptimal at mixing compared to the mix-norm minimising perturbations, and also that minimising the mix-norm remains (for density-stratified flows) a good proxy for identifying perturbations which minimise the variance at long times. Although increasing stratification reduces the mixing in general, mix-norm minimising optimal perturbations can still trigger substantial mixing for $Ri_{b}\lesssim 0.3$. By considering the time evolution of the kinetic energy and potential energy reservoirs, we find that such perturbations lead to a flow which, through Taylor dispersion, very effectively converts perturbation kinetic energy into ‘available potential energy’, which in turn leads rapidly and irreversibly to thorough and efficient mixing, with little energy returned to the kinetic-energy reservoirs.
Saturation of a turbulent mixing layer over a cavity: response to harmonic forcing around mean flows
- E. Boujo, M. Bauerheim, N. Noiray
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 386-418
-
- Article
- Export citation
-
Turbulent mixing layers over cavities can couple with acoustic waves and lead to undesired oscillations. To understand the nonlinear aspects of this phenomenon, a turbulent mixing layer over a deep cavity is considered and its response to harmonic forcing is analysed with large-eddy simulations (LES) and linearised Navier–Stokes equations (LNSE). The Reynolds number is $Re=150\,000$. As a model of incoming acoustic perturbations, spatially uniform time-harmonic velocity forcing is applied at the cavity end, with amplitudes spanning the wide range 0.045–8.9 % of the main channel bulk velocity. Compressible LES provide reference nonlinear responses of the shear layer, and the associated mean flows. Linear responses are calculated with the incompressible LNSE around the LES mean flows; they predict well the amplification (both measured with kinetic energy and with a proxy for vortex sound production in the mixing layer) and capture the nonlinear saturation observed as the forcing amplitude increases and the mixing layer thickens. Perhaps surprisingly, LNSE calculations based on a monochromatic (single-frequency) assumption yield a good agreement even though higher harmonics and their nonlinear interaction (Reynolds stresses) are not negligible. However, it is found that the leading Reynolds stresses do not force the mixing layer efficiently, as shown by a comparison with the optimal volume forcing obtained from a resolvent analysis. Therefore they cannot fully benefit from the potential for amplification available in the flow. Finally, the sensitivity of the optimal harmonic forcing at the cavity end is computed with an adjoint method. The sensitivities to mean flow modification and to a localised feedback (structural sensitivity) both identify the upstream cavity corner as the region where a small-amplitude modification has the strongest effect. This can guide in a systematic way the design of strategies aiming at controlling the amplification and saturation mechanisms.
On the over-production of turbulence beneath surface waves in Reynolds-averaged Navier–Stokes models
- Bjarke Eltard Larsen, David R. Fuhrman
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 419-460
-
- Article
- Export citation
-
In previous computational fluid dynamics studies of breaking waves, there has been a marked tendency to severely over-estimate turbulence levels, both pre- and post-breaking. This problem is most likely related to the previously described (though not sufficiently well recognized) conditional instability of widely used turbulence models when used to close Reynolds-averaged Navier–Stokes (RANS) equations in regions of nearly potential flow with finite strain, resulting in exponential growth of the turbulent kinetic energy and eddy viscosity. While this problem has been known for nearly 20 years, a suitable and fundamentally sound solution has yet to be developed. In this work it is demonstrated that virtually all commonly used two-equation turbulence closure models are unconditionally, rather than conditionally, unstable in such regions. A new formulation of the $k$–$\unicode[STIX]{x1D714}$ closure is developed which elegantly stabilizes the model in nearly potential flow regions, with modifications remaining passive in sheared flow regions, thus solving this long-standing problem. Computed results involving non-breaking waves demonstrate that the new stabilized closure enables nearly constant form wave propagation over long durations, avoiding the exponential growth of the eddy viscosity and inevitable wave decay exhibited by standard closures. Additional applications on breaking waves demonstrate that the new stabilized model avoids the unphysical generation of pre-breaking turbulence which widely plagues existing closures. The new model is demonstrated to be capable of predicting accurate pre- and post-breaking surface elevations, as well as turbulence and undertow velocity profiles, especially during transition from pre-breaking to the outer surf zone. Results in the inner surf zone are similar to standard closures. Similar methods for formally stabilizing other widely used closure models ($k$–$\unicode[STIX]{x1D714}$ and $k$–$\unicode[STIX]{x1D700}$ variants) are likewise developed, and it is recommended that these be utilized in future RANS simulations of surface waves. (In the above $k$ is the turbulent kinetic energy density, $\unicode[STIX]{x1D714}$ is the specific dissipation rate, and $\unicode[STIX]{x1D700}$ is the dissipation.)
A volume integral implementation of the Goldstein generalised acoustic analogy for unsteady flow simulations
- Vasily A. Semiletov, Sergey A. Karabasov
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 461-487
-
- Article
- Export citation
-
A new volume integral method based on the Goldstein generalised acoustic analogy is developed and directly applied with large-eddy simulation (LES). In comparison with the existing Goldstein generalised acoustic analogy implementations, the current method does not require the computation and recording of the expensive fluctuating stress autocovariance function in the seven-dimensional space–time. Until now, the multidimensional complexity of the generalised acoustic analogy source term has been the main barrier to using it in routine engineering calculations. The new method only requires local pointwise stresses as an input that can be routinely computed during the flow simulation. On the other hand, the new method is mathematically equivalent to the original Goldstein acoustic analogy formulation, and, thus, allows for a direct correspondence between different effective noise sources in the jet and the far-field noise spectra. The implementation is performed for conditions of a high-speed subsonic isothermal jet corresponding to the Rolls-Royce SILOET experiment and uses the LES solution based on the CABARET solver. The flow and noise solutions are validated by comparison with experiment. The accuracy and robustness of the integral volume implementation of the generalised acoustic analogy are compared with those based on the standard Ffowcs Williams–Hawkings surface integral method and the conventional Lighthill acoustic analogy. As a demonstration of its capabilities to investigate jet noise mechanisms, the new integral volume method is applied to analyse the relative importance of various noise generation and propagation components within the Goldstein generalised acoustic analogy model.
Turbulent channel flow of an elastoviscoplastic fluid
- Marco E. Rosti, Daulet Izbassarov, Outi Tammisola, Sarah Hormozi, Luca Brandt
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 488-514
-
- Article
- Export citation
-
We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier–Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarizes. These different behaviours are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centreline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high-speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction.
An analytic solution for the noise generated by gust–aerofoil interaction for plates with serrated leading edges
- Lorna J. Ayton, Jae Wook Kim
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 515-536
-
- Article
- Export citation
-
This paper presents an analytic solution for the sound generated by an unsteady gust interacting with a semi-infinite flat plate with a serrated leading edge in a background steady uniform flow. Viscous and nonlinear effects are neglected. The Wiener–Hopf method is used in conjunction with a non-orthogonal coordinate transformation and separation of variables to permit analytical progress. The solution is obtained in terms of a modal expansion in the spanwise coordinate; however, for low- and mid-range incident frequencies only the zeroth-order mode is seen to contribute to the far-field acoustics, therefore the far-field noise can be quickly evaluated. The solution gives insight into the potential mechanisms behind the reduction of noise for plates with serrated leading edges compared to those with straight edges, and predicts a logarithmic dependence between the tip-to-root serration height and the decrease of far-field noise. The two mechanisms behind the noise reduction are proposed to be an increased destructive interference in the far field, and a redistribution of acoustic energy from low cut-on modes to higher cut-off modes as the tip-to-root serration height is increased. The analytic results show good agreement in comparison with experimental measurements. The results are also compared against nonlinear numerical predictions where good agreement is also seen between the two results as frequency and tip-to-root ratio are varied.
Large-eddy simulation of flow over an axisymmetric body of revolution
- Praveen Kumar, Krishnan Mahesh
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 537-563
-
- Article
- Export citation
-
Wall-resolved large-eddy simulation (LES) is used to simulate flow over an axisymmetric body of revolution at a Reynolds number, $Re=1.1\times 10^{6}$, based on the free-stream velocity and the length of the body. The geometry used in the present work is an idealized submarine hull (DARPA SUBOFF without appendages) at zero angle of pitch and yaw. The computational domain is chosen to avoid confinement effects and capture the wake up to fifteen diameters downstream of the body. The unstructured computational grid is designed to capture the fine near-wall flow structures as well as the wake evolution. LES results show good agreement with the available experimental data. The axisymmetric turbulent boundary layer has higher skin friction and higher radial decay of turbulence away from the wall, compared to a planar turbulent boundary layer under similar conditions. The mean streamwise velocity exhibits self-similarity, but the turbulent intensities are not self-similar over the length of the simulated wake, consistent with previous studies reported in the literature. The axisymmetric wake shifts from high-$Re$ to low-$Re$ equilibrium self-similar solutions, which were only observed for axisymmetric wakes of bluff bodies in the past.
High frequency resonant response of a monopile in irregular deep water waves
- Bjørn Hervold Riise, John Grue, Atle Jensen, Thomas B. Johannessen
-
- Published online by Cambridge University Press:
- 23 August 2018, pp. 564-586
-
- Article
- Export citation
-
Experiments with a weakly damped monopile, either fixed or free to oscillate, exposed to irregular waves in deep water, obtain the wave-exciting moment and motion response. The nonlinearity and peak wavenumber cover the ranges: $\unicode[STIX]{x1D716}_{P}\sim 0.10{-}0.14$ and $k_{P}r\sim 0.09{-}0.14$ where $\unicode[STIX]{x1D716}_{P}=0.5H_{S}k_{P}$ is an estimate of the spectral wave slope, $H_{S}$ the significant wave height, $k_{P}$ the peak wavenumber and $r$ the cylinder radius. The response and its statistics, expressed in terms of the exceedance probability, are discussed as a function of the resonance frequency, $\unicode[STIX]{x1D714}_{0}$ in the range $\unicode[STIX]{x1D714}_{0}\sim 3{-}5$ times the spectral peak frequency, $\unicode[STIX]{x1D714}_{P}$. For small wave slope, long waves and $\unicode[STIX]{x1D714}_{0}/\unicode[STIX]{x1D714}_{P}=3$, the nonlinear response deviates only very little from its linear counterpart. However, the nonlinearity becomes important for increasing wave slope, wavenumber and resonance frequency ratio. The extreme response events are found in a region where the Keulegan–Carpenter number exceeds $KC>5$, indicating the importance of possible flow separation effects. A similar region is also covered by a Froude number exceeding $Fr>0.4$, pointing to surface gravity wave effects at the scale of the cylinder diameter. Regarding contributions to the higher harmonic forces, different wave load mechanisms are identified, including: (i) wave-exciting inertia forces, a function of the fluid acceleration; (ii) wave slamming due to both non-breaking and breaking wave events; (iii) a secondary load cycle; and (iv) possible drag forces, a function of the fluid velocity. Also, history effects due to the inertia of the moving pile, contribute to the large response events. The ensemble means of the third, fourth and fifth harmonic wave-exciting force components extracted from the irregular wave results are compared to the third harmonic FNV (Faltinsen, Newman and Vinje) theory as well as other available experiments and calculations. The present irregular wave measurements generalize results obtained in deep water regular waves.