scholarly journals Parabolized stability analysis of jets from serrated nozzles

2016 ◽  
Vol 789 ◽  
pp. 36-63 ◽  
Author(s):  
Aniruddha Sinha ◽  
Kristján Gudmundsson ◽  
Hao Xia ◽  
Tim Colonius

We study the viscous spatial linear stability characteristics of the time-averaged flow in turbulent subsonic jets issuing from serrated (chevroned) nozzles, and compare them to analogous round jet results. Linear parabolized stability equations (PSE) are used in the calculations to account for the non-parallel base flow. By exploiting the symmetries of the mean flow due to the regular arrangement of serrations, we obtain a series of coupled two-dimensional PSE problems from the original three-dimensional problem. This reduces the solution cost and manifests the symmetries of the stability modes. In the parallel-flow linear stability theory (LST) calculations that are performed near the nozzle to initiate the PSE, we find that the serrated nozzle reduces the growth rates of the most unstable eigenmodes of the jet, but their phase speeds are approximately similar. We obtain encouraging validation of our linear PSE instability wave results vis-à-vis near-field hydrodynamic pressure data acquired on a phased microphone array in experiments, after filtering the latter with proper orthogonal decomposition (POD) to extract the energetically dominant coherent part. Additionally, a large-eddy simulation database of the same serrated jet is investigated, and its POD-filtered pressure field is found to compare favourably with the corresponding PSE solution within the jet plume. We conclude that the coherent hydrodynamic pressure fluctuations of jets from both round and serrated nozzles are reasonably consistent with the linear instability modes of the turbulent mean flow.

2011 ◽  
Vol 689 ◽  
pp. 97-128 ◽  
Author(s):  
K. Gudmundsson ◽  
Tim Colonius

AbstractPrevious work has shown that aspects of the evolution of large-scale structures, particularly in forced and transitional mixing layers and jets, can be described by linear and nonlinear stability theories. However, questions persist as to the choice of the basic (steady) flow field to perturb, and the extent to which disturbances in natural (unforced), initially turbulent jets may be modelled with the theory. For unforced jets, identification is made difficult by the lack of a phase reference that would permit a portion of the signal associated with the instability wave to be isolated from other, uncorrelated fluctuations. In this paper, we investigate the extent to which pressure and velocity fluctuations in subsonic, turbulent round jets can be described aslinearperturbations to the mean flow field. The disturbances are expanded about the experimentally measured jet mean flow field, and evolved using linear parabolized stability equations (PSE) that account, in an approximate way, for the weakly non-parallel jet mean flow field. We utilize data from an extensive microphone array that measures pressure fluctuations just outside the jet shear layer to show that, up to an unknown initial disturbance spectrum, the phase, wavelength, and amplitude envelope of convecting wavepackets agree well with PSE solutions at frequencies and azimuthal wavenumbers that can be accurately measured with the array. We next apply the proper orthogonal decomposition to near-field velocity fluctuations measured with particle image velocimetry, and show that the structure of the most energetic modes is also similar to eigenfunctions from the linear theory. Importantly, the amplitudes of the modes inferred from the velocity fluctuations are in reasonable agreement with those identified from the microphone array. The results therefore suggest that, to predict, with reasonable accuracy, the evolution of the largest-scale structures that comprise the most energetic portion of the turbulent spectrum of natural jets, nonlinear effects need only be indirectly accounted for by considering perturbations to the mean turbulent flow field, while neglecting any non-zero frequency disturbance interactions.


2011 ◽  
Vol 681 ◽  
pp. 261-292 ◽  
Author(s):  
M. KEARNEY-FISCHER ◽  
J.-H. KIM ◽  
M. SAMIMY

Mach wave radiation is one of the better understood sources of jet noise. However, the exact conditions of its onset are difficult to determine and the literature to date typically explores Mach wave radiation well above its onset conditions. In order to determine the conditions for the onset of Mach wave radiation and to explore its behaviour during onset and beyond, three ideally expanded jets with Mach numbers Mj = 0.9, 1.3 and 1.65 and stagnation temperature ratios ranging over To/T∞ = 1.0–2.5 (acoustic Mach number 0.83–2.10) were used. Data are collected using a far-field microphone array, schlieren imaging and streamwise two-component particle image velocimetry. Using arc filament plasma actuators to force the jet provides an unprecedented tool for detailed examination of Mach wave radiation. The response of the jet to various forcing parameters (combinations of one azimuthal mode m = 0, 1 and 3 and one Strouhal number StDF = 0.09–3.0) is explored. Phase-averaged schlieren images clearly show the onset and evolution of Mach wave radiation in response to both changes in the jet operating conditions and forcing parameters. It is observed that Mach wave radiation is initiated as a coalescing of the near-field hydrodynamic pressure fluctuations in the immediate vicinity of the large-scale structures. As the jet exit velocity increases, the hydrodynamic pressure fluctuations coalesce, first into a curved wavefront, then flatten into the conical wavefronts commonly associated with Mach wave radiation. The results show that the largest and most coherent structures (e.g. forcing with m = 0 and StDF ~ 0.3) produce the strongest Mach wave radiation. Conversely, Mach wave radiation is weakest when the structures are the least coherent (e.g. forcing with m = 3 and StDF > 1.5).


2019 ◽  
Vol 865 ◽  
Author(s):  
Cheng-Nian Xiao ◽  
Inanc Senocak

We investigate the stability of the Prandtl model for katabatic slope flows using both linear stability theory and direct numerical simulations. Starting from Prandtl’s analytical solution for uniformly cooled laminar slope flows, we use linear stability theory to identify the onset of instability and features of the most unstable modes. Our results show that the Prandtl model for parallel katabatic slope flows is prone to transverse and longitudinal modes of instability. The transverse mode of instability manifests itself as stationary vortical flow structures aligned in the along-slope direction, whereas the longitudinal mode of instability emerges as waves propagating in the base-flow direction. Beyond the stability limits, these two modes of instability coexist and form a complex flow structure crisscrossing the plane of flow. The emergence of a particular form of these instabilities depends strongly on three dimensionless parameters, which are the slope angle, the Prandtl number and a newly introduced stratification perturbation parameter, which is proportional to the relative importance of the disturbance to the background stratification due to the imposed surface buoyancy flux. We demonstrate that when this parameter is sufficiently large, then the stabilising effect of the background stratification can be overcome. For shallow slopes, the transverse mode of instability emerges despite meeting the Miles–Howard stability criterion of $Ri>0.25$. At steep slope angles, slope flow can remain linearly stable despite attaining Richardson numbers as low as $3\times 10^{-3}$.


2015 ◽  
Vol 778 ◽  
pp. 120-132 ◽  
Author(s):  
Mario Weder ◽  
Michael Gloor ◽  
Leonhard Kleiser

We present a decomposition of the temporal growth rate ${\it\omega}_{i}$ which characterises the evolution of wave-like disturbances in linear stability theory for compressible flows. The decomposition is based on the disturbance energy balance by Chu (Acta Mech., vol. 1 (3), 1965, pp. 215–234) and provides terms for production, dissipation and flux of energy as components of ${\it\omega}_{i}$. The inclusion of flux terms makes our formulation applicable to unconfined flows and flows with permeable or vibrating boundaries. The decomposition sheds light on the fundamental mechanisms determining temporal growth or decay of disturbances. The additional insights gained by the proposed approach are demonstrated by an investigation of two model flows, namely compressible Couette flow and a plane compressible jet.


1996 ◽  
Vol 326 ◽  
pp. 399-415 ◽  
Author(s):  
M. Wanschura ◽  
H. C. Kuhlmann ◽  
H. J. Rath

The stability of steady axisymmetric convection in cylinders heated from below and insulated laterally is investigated numerically using a mixed finite-difference/Chebyshev collocation method to solve the base flow and the linear stability equations. Linear stability boundaries are given for radius to height ratios γ from 0.9 to 1.56 and for Prandtl numbers Pr = 0.02 and Pr = 1. Depending on γ and Pr, the azimuthal wavenumber of the critical mode may be m = 1, 2, 3, or 4. The dependence of the critical Rayleigh number on the aspect ratio and the instability mechanisms are explained by analysing the energy transfer to the critical modes for selected cases. In addition to these results the onset of buoyant convection in liquid bridges with stress-free conditions on the cylindrical surface is considered. For insulating thermal boundary conditions, the onset of convection is never axisymmetric and the critical azimuthal wavenumber increases monotonically with γ. The critical Rayleigh number is less then 1708 for most aspect ratios.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1994 ◽  
Vol 279 ◽  
pp. 351-375 ◽  
Author(s):  
Karim Shariff ◽  
Roberto Verzicco ◽  
Paolo Orlandi

Finite-difference calculations with random and single-mode perturbations are used to study the three-dimensional instability of vortex rings. The basis of current understanding of the subject consists of a heuristic inviscid model (Widnall, Bliss & Tsai 1974) and a rigorous theory which predicts growth rates for thin-core uniform vorticity rings (Widnall & Tsai 1977). At sufficiently high Reynolds numbers the results correspond qualitatively to those predicted by the heuristic model: multiple bands of wavenumbers are amplified, each band having a distinct radial structure. However, a viscous correction factor to the peak inviscid growth rate is found. It is well described by the first term, 1 – α1(β)/Res, for a large range of Res. Here Res is the Reynolds number defined by Saffman (1978), which involves the curvature-induced strain rate. It is found to be the appropriate choice since then α1(β) varies weakly with core thickness β. The three most nonlinearly amplified modes are a mean azimuthal velocity in the form of opposing streams, an n = 1 mode (n is the azimuthal wavenumber) which arises from the interaction of two second-mode bending waves and the harmonic of the primary second mode. When a single wave is excited, higher harmonics begin to grow successively later with nonlinear growth rates proportional to n. The modified mean flow has a doubly peaked azimuthal vorticity. Since the curvature-induced strain is not exactly stagnation-point flow there is a preference for elongation towards the rear of the ring: the outer structure of the instability wave forms a long wake consisting of n hairpin vortices whose waviness is phase shifted π/n relative to the waviness in the core. Whereas the most amplified linear mode has three radial layers of structure, higher radial modes having more layers of radial structure (hairpins piled upon hairpins) are excited when the initial perturbation is large, reminiscent of visualization experiments on the formation of a turbulent ring at the generator.


2007 ◽  
Vol 587 ◽  
pp. 173-215 ◽  
Author(s):  
PRASUN K. RAY ◽  
SANJIVA K. LELE

Broadband shock-associated noise is an important component of the overall noise generated by modern airplanes. In this study, sound generated by the weakly nonlinear interaction between linear instability waves and the shock-cell structure in supersonic jets is investigated numerically in order to gain insight into the broadband shock-noise problem. The model formulation decomposes the overall flow into a mean flow, linear instability waves, the shock-cell structure and shock-noise. The mean flow is obtained by solving RANSequations with a k-ε model. Locally parallel stability equations are solved for the shock structure, and linear parabolized stability equations are solved for the instability waves. Then, source terms representing the instability wave/shock-cell interaction are assembled and the inhomogeneous linearized Euler equations are solved for the shock-noise.Three cases are considered, a cold under-expanded Mj = 1.22 jet, a hot under-expanded Mj = 1.22 jet, and a cold over-expanded Mj = 1.36 jet.Shock-noise computations are used to identify and understand significant trends in peak sound amplitudes and radiation angles. The peak sound radiation angles are explained well with the Mach wave model of Tam & Tanna J. Sound Vib. Vol. 81, 1982, p. 337). The observed reduction of peak sound amplitudes with frequency correlates well with the corresponding reduction of instability wave growth with frequency. However, in order to account for variation of sound amplitude for different azimuthal modes, the radial structure of the instability waves must be considered in additionto streamwise growth. The effect of heating on the Mj = 1.22 jet is shown to enhance the sound radiated due to the axisymmetric instability waves while the other modesare relatively unaffected. Solutions to a Lilley–Goldstein equation show that soundgenerated by ‘thermodynamic’ source terms is small relative to sound from ‘momentum’ sources though heating does increase the relative importance of the thermodynamic source. Furthermore, heating preferentially amplifies sound associated with the axisymmetric modes owing to constructive interference between sound from the momentumand thermodynamic sources. However, higher modes show destructive interference between these two sources and are relatively unaffected by heating.


Sign in / Sign up

Export Citation Format

Share Document