scholarly journals Wavepackets and trapped acoustic modes in a turbulent jet: coherent structure eduction and global stability

2017 ◽  
Vol 825 ◽  
pp. 1153-1181 ◽  
Author(s):  
Oliver T. Schmidt ◽  
Aaron Towne ◽  
Tim Colonius ◽  
André V. G. Cavalieri ◽  
Peter Jordan ◽  
...  

Coherent features of a turbulent Mach 0.9, Reynolds number$10^{6}$jet are educed from a high-fidelity large eddy simulation. Besides the well-known Kelvin–Helmholtz instabilities of the shear layer, a new class of trapped acoustic waves is identified in the potential core. A global linear stability analysis based on the turbulent mean flow is conducted. The trapped acoustic waves form branches of discrete eigenvalues in the global spectrum, and the corresponding global modes accurately match the educed structures. Discrete trapped acoustic modes occur in a hierarchy determined by their radial and axial order. A local dispersion relation is constructed from the global modes and found to agree favourably with an empirical dispersion relation educed from the simulation data. The product between direct and adjoint modes is then used to isolate the trapped waves. Under certain conditions, resonance in the form of a beating occurs between trapped acoustic waves of positive and negative group velocities. This resonance explains why the trapped modes are prominently observed in the simulation and as tones in previous experimental studies. In the past, these tones were attributed to external factors. Here, we show that they are an intrinsic feature of high-subsonic jets that can be unambiguously identified by a global linear stability analysis.

Author(s):  
Thomas L. Kaiser ◽  
Thierry Poinsot ◽  
Kilian Oberleithner

The hydrodynamic instability in an industrial, two-staged, counter-rotative, swirled injector of highly complex geometry is under investigation. Large eddy simulations show that the complicated and strongly nonparallel flow field in the injector is superimposed by a strong precessing vortex core. Mean flow fields of large eddy simulations, validated by experimental particle image velocimetry measurements are used as input for both local and global linear stability analysis. It is shown that the origin of the instability is located at the exit plane of the primary injector. Mode shapes of both global and local linear stability analysis are compared to a dynamic mode decomposition based on large eddy simulation snapshots, showing good agreement. The estimated frequencies for the instability are in good agreement with both the experiment and the simulation. Furthermore, the adjoint mode shapes retrieved by the global approach are used to find the best location for periodic forcing in order to control the precessing vortex core.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 212
Author(s):  
Miles Owen ◽  
Abdelkader Frendi

The results from a temporal linear stability analysis of a subsonic boundary layer over a flat plate with a straight and wavy leading edge are presented in this paper for a swept and un-swept plate. For the wavy leading-edge case, an extensive study on the effects of the amplitude and wavelength of the waviness was performed. Our results show that the wavy leading edge increases the critical Reynolds number for both swept and un-swept plates. For the un-swept plate, increasing the leading-edge amplitude increased the critical Reynolds number, while changing the leading-edge wavelength had no effect on the mean flow and hence the flow stability. For the swept plate, a local analysis at the leading-edge peak showed that increasing the leading-edge amplitude increased the critical Reynolds number asymptotically, while the leading-edge wavelength required optimization. A global analysis was subsequently performed across the span of the swept plate, where smaller leading-edge wavelengths produced relatively constant critical Reynolds number profiles that were larger than those of the straight leading edge, while larger leading-edge wavelengths produced oscillating critical Reynolds number profiles. It was also found that the most amplified wavenumber was not affected by the wavy leading-edge geometry and hence independent of the waviness.


2018 ◽  
Vol 48 (11) ◽  
pp. 2779-2797 ◽  
Author(s):  
Megan A. Stamper ◽  
John R. Taylor ◽  
Baylor Fox-Kemper

AbstractMotivated by recent observations of submesoscales in the Southern Ocean, we use nonlinear numerical simulations and a linear stability analysis to examine the influence of a barotropic jet on submesoscale instabilities at an isolated front. Simulations of the nonhydrostatic Boussinesq equations with a strong barotropic jet (approximately matching the observed conditions) show that submesoscale disturbances and strong vertical velocities are confined to a small region near the initial frontal location. In contrast, without a barotropic jet, submesoscale eddies propagate to the edges of the computational domain and smear the mean frontal structure. Several intermediate jet strengths are also considered. A linear stability analysis reveals that the barotropic jet has a modest influence on the growth rate of linear disturbances to the initial conditions, with at most a ~20% reduction in the growth rate of the most unstable mode. On the other hand, a basic state formed by averaging the flow at the end of the simulation with a strong barotropic jet is linearly stable, suggesting that nonlinear processes modify the mean flow and stabilize the front.


2015 ◽  
Vol 765 ◽  
pp. 150-166 ◽  
Author(s):  
Xiaoyi Lu ◽  
Carlos Pantano

AbstractOne-dimensional planar premixed flames propagating in a uniform flow are susceptible to hydrodynamic instabilities known (generically) as Darrieus–Landau instabilities. Here, we extend that hydrodynamic linear stability analysis to include a lateral shear. This generalization is a situation of interest for laminar and turbulent flames when they travel into a region of shear (such as a jet or shear layer). It is shown that the problem can be formulated and solved analytically and a dispersion relation can be determined. The solution depends on a shear parameter in addition to the wavenumber, thermal expansion ratio, and Markstein lengths. The study of the dispersion relation shows that perturbations have two types of behaviour as wavenumber increases. First, for small shear, we recover the Darrieus–Landau results except for a region at small wavenumbers, large wavelengths, that is stable. Initially, increasing shear has a stabilizing effect. But, for sufficiently high shear, the flame becomes unstable again and its most unstable wavelength can be much smaller than the Markstein length of the zero-shear flame. Finally, the stabilizing effect of low shear can make flames with negative Markstein numbers stable within a band of wavenumbers.


2015 ◽  
Vol 45 (9) ◽  
pp. 2189-2208 ◽  
Author(s):  
Florian Sévellec ◽  
Thierry Huck

AbstractA weakly damped mode of variability, corresponding to the oceanic signature of the Atlantic multidecadal oscillation (AMO) was found through the linear stability analysis of a realistic ocean general circulation model. A simple two-level model was proposed to rationalize both its period and damping rate. This model is extended here to three levels to investigate how the mode can draw energy from the mean flow, as found in various ocean and coupled models. A linear stability analysis in this three-level model shows that the positive growth rate of the oscillatory mode depends on the zonally averaged isopycnal slope. This mode corresponds to a westward propagation of density anomalies in the pycnocline, typical of large-scale baroclinic Rossby waves. The most unstable mode corresponds to the largest scale one (at least for low isopycnal slope). The mode can be described in four phases composing a full oscillation cycle: 1) basin-scale warming of the North Atlantic (AMO positive phase), 2) decrease in upper-ocean poleward transport [hence a reduction of the Atlantic meridional overturning circulation (AMOC)], 3) basin-scale cooling (negative AMO), and 4) AMOC intensification. A criterion is developed to test, in oceanic datasets or numerical models, whether this multidecadal oscillation is an unstable oceanic internal mode of variability or if it is stable and externally forced. Consistent with the classical theory of baroclinic instability, this criterion depends on the vertical structure of the mode. If the upper pycnocline signature is in advance of the deeper pycnocline signature with respect to the westward propagation, the mode is unstable and could be described as an oceanic internal mode of variability.


Author(s):  
Jens S. Müller ◽  
Finn Lückoff ◽  
Kilian Oberleithner

The fundamental impact of the precessing vortex core (PVC) as a dominant coherent flow structure in the flow field of swirl-stabilized gas turbine combustors has still not been investigated in depth. In order to do so, the PVC needs to be actively controlled to be able to set its parameters independently to any other of the combustion system. In this work, open-loop actuation is applied in the mixing section between the swirler and the generic combustion chamber of a non-reacting swirling jet setup to investigate the receptivity of the PVC with regard to its lock-in behavior at different streamwise positions. The mean flow in the mixing section as well as in the combustion chamber is measured by stereoscopic particle image velocimetry and the PVC is extracted from the snapshots using proper orthogonal decomposition. The lock-in experiments reveal the axial position in the mixing section that is most suitable for actuation. Furthermore, a global linear stability analysis is conducted to determine the adjoint mode of the PVC which reveals the regions of highest receptivity to periodic actuation based on mean flow input only. This theoretical receptivity model is compared with the experimentally obtained receptivity data and the applicability of the adjoint-based model for the prediction of optimal actuator designs is discussed.


2021 ◽  
Vol 927 ◽  
Author(s):  
Romain Canu ◽  
Marie-Charlotte Renoult

We performed a linear stability analysis of a Newtonian ferrofluid cylinder surrounded by a Newtonian non-magnetic fluid in an azimuthal magnetic field. A wire is used at the centre of the ferrofluid cylinder to create this magnetic field. Isothermal conditions are considered and gravity is ignored. An axisymmetric perturbation is imposed at the interface between the two fluids and a dispersion relation is obtained allowing us to predict whether the flow is stable or unstable with respect to this perturbation. This relation is dependent on the Ohnesorge number of the ferrofluid, the dynamic viscosity ratio, the density ratio, the magnetic Bond number, the relative magnetic permeability and the dimensionless wire radius. Solutions to this dispersion relation are compared with experimental data from Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and, more recently, Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502). A better agreement than the inviscid theory and the theory that only takes into account the viscosity of the ferrofluid is shown with the data of Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and those of Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502) for small wavenumbers.


2017 ◽  
Vol 832 ◽  
pp. 97-145 ◽  
Author(s):  
Kuchimanchi K. Bharadwaj ◽  
Debopam Das

The present work investigates the puffing instability of circular buoyant plumes by performing global linear stability analysis and experiments. In the non-dimensional parameter space investigated, plumes exhibit global instability only for axisymmetric perturbations with two unstable modes, which are of oscillatory type. The frequencies of these two unstable global modes agree well with the experiments which suggest that puffing occurs in buoyant plumes as a result of linear global instability. A comprehensive investigation on the effect of various non-dimensional parameters and inlet velocity profiles on frequency and growth rates of the global modes is carried out. The results are used to delineate the stability boundaries for these global modes and to obtain scaling laws for the associated oscillation frequencies. The analysis demonstrates that the two buoyancy parameters, Froude number and source-to-ambient density ratio, play dominant roles in impacting plume transition and oscillation frequencies. Results from global linear stability analysis and earlier experiments have majorly differed in two aspects. The earlier experiments reported a switch in puffing frequency scaling in Richardson number range 100–500, while the instability analysis predicts this switch at around 6000. Also, the instability analysis predicts the occurrence of puffing at density ratios higher than the critical value 0.5–0.6 reported in earlier experiments. To address these differences and validate the results obtained from global linear stability analysis, experiments are performed in a set-up that has been carefully designed to minimize the settling chamber disturbances. The present experiments corroborate the findings of global linear stability analysis. The mechanisms responsible for global instability in plumes have been identified using perturbation vorticity transport equation.


Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner ◽  
Laurent Selle ◽  
Thierry Poinsot

Abstract Shape optimization with respect to the suppression or enhancement of dynamical flow structures is an important topic in combustion research and beyond. In this paper, we investigate the flow in an industrial fuel injection system by experimental means, as well as large eddy simulation (LES) and linear stability analysis (LSA) for two configurations of the swirler. In the first configuration, the reference geometry, a precessing vortex core (PVC) occurs. In the second configuration, a center body is mounted in the interior of the injector. It is shown by both experiments and LES that the PVC is suppressed by the presence of the center body, while the mean flow remains nearly unaffected. The method of LSA is applied in order to explain the effect of the geometry change. The work shows that LSA is capable of explaining the occurrence or disappearance of coherent structures evolving on the turbulent flows if the geometry is changed. This is an important step in using LSA in the context of shape optimization of industrial fuel injectors.


Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner ◽  
Laurent Selle ◽  
Thierry Poinsot

Abstract Shape optimization with respect to the suppression or enhancement of dynamical flow structures is an important topic in combustion research and beyond. In this paper, we investigate the flow in an industrial fuel injection system by experimental means, as well as Large Eddy Simulation (LES) and Linear Stability Analysis (LSA) for two configurations of the swirler. In the first configuration, the reference geometry, a Precessing Vortex Core (PVC) occurs. In the second configuration, a center body is mounted in the interior of the injector. It is shown by both experiments and LES that the PVC is suppressed in the presence of the center body, while the mean flow remains nearly unaffected. The method of LSA is applied in order to explain the effect of the geometry change. The work shows that LSA is capable of explaining the occurrence or disappearance of coherent structures evolving on the turbulent flows if the geometry is changed. This is an important step in using LSA in the context of shape optimization of industrial fuel injectors.


Sign in / Sign up

Export Citation Format

Share Document