skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Comparison of Implicit and Explicit Simulations of Vortex Breakdown
Vortex breakdown is an instability which occurs in swirling flows and is seen in both incompressible and compressible flow regimes. The choice of numerical algorithm to com- pute this flow is critical to resolving the relevant physical phenomena. In previous work, we described the procedure of combining a Barely Implicit Correction (BIC) algorithm with a fourth-order Flux-Corrected Transport (FCT) algorithm to eliminate the sound speed limit for explicit calculations of low-speed flows, and showed that BIC-FCT could predict the three major vortex breakdown modes (spiral, bubble and double-helix). Here, we show a direct comparison of the implicit BIC and explicit FCT algorithms by using the same time step within the explicit stability limit. The comparison is conducted on simulations of a three-dimensional, swirling jet flow with vortex breakdown, by examining the flow structures and examining the property fields. Differences are observed in the phase of the downstream spiral and double-helix modes. Based on a pressure wave analysis, we attribute this phase shift to the difference of how the explicit FCT and implicit BIC treat an open outflow boundary.  more » « less
Award ID(s):
1839510
PAR ID:
10104190
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
AIAA Scitech 2019 Forum
Page Range / eLocation ID:
1144
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Global instabilities in swirling flows can significantly alter the flame and flow dynamics of swirl-stabilized flames, such as those in modern gas turbine engines. In this study, we characterize the interaction between the precessing vortex core (PVC), which is the consequence of a global hydrodynamic instability, and thermoacoustic instabilities, which are the result of a coupling between combustor acoustics and the unsteady heat release rate. This study is performed using experimental data obtained from a model gas turbine combustor employing two concentric swirling nozzles of air, separated by a ring of fuel injectors, operating at five bar pressure. The flow split between the two streams is systematically varied to observe the impact of flow structure variation on the system dynamics at both non-reacting and reacting conditions. High-speed stereoscopic particle image velocimetry, OH planar laser-induced fluorescence and acetone planar laser-induced fluorescence are used to obtain information about the velocity fields, flame and fuel flow behaviour, respectively. Spectral proper orthogonal decomposition and a complex network analysis are used to identify and characterize the dominant oscillation mechanisms driving the system. In the non-reacting data, a PVC is present in most cases and the amplitude of the oscillation increases with increasing flow through the centre nozzle. In the reacting data, three dominant modes are seen: two thermoacoustic modes and the PVC. Our results show that in the cases where the frequency of the PVC overlaps with either of the thermoacoustic modes, the thermoacoustic modes are suppressed. The complex network analysis coupled with a weakly nonlinear theoretical analysis suggests the mechanisms by which this coupling and suppression of the thermoacoustic mode occur. 
    more » « less
  2. We investigate the effect of inertial particles dispersed in a circular patch of finite radius on the stability of a two-dimensional Rankine vortex in semi-dilute dusty flows. Unlike the particle-free case where no unstable modes exist, we show that the feedback force from the particles triggers a novel instability. The mechanisms driving the instability are characterized using linear stability analysis for weakly inertial particles and further validated against Eulerian–Lagrangian simulations. We show that the particle-laden vortex is always unstable if the mass loading $M>0$ . Surprisingly, even non-inertial particles destabilize the vortex by a mechanism analogous to the centrifugal Rayleigh–Taylor instability in radially stratified vortex with density jump. We identify a critical mass loading above which an eigenmode $$m$$ becomes unstable. This critical mass loading drops to zero as $$m$$ increases. When particles are inertial, modes that fall below the critical mass loading become unstable, whereas modes above it remain unstable but with lower growth rates compared with the non-inertial case. Comparison with Eulerian–Lagrangian simulations shows that growth rates computed from simulations match well the theoretical predictions. Past the linear stage, we observe the emergence of high-wavenumber modes that turn into spiralling arms of concentrated particles emanating out of the core, while regions of particle-free flow are sucked inward. The vorticity field displays a similar pattern which leads to the breakdown of the initial Rankine structure. This novel instability for a dusty vortex highlights how the feedback force from the disperse phase can induce the breakdown of an otherwise resilient vortical structure. 
    more » « less
  3. This paper presents a new representation of curve dynamics, with applications to vortex filaments in fluid dynamics. Instead of representing these filaments with explicit curve geometry and Lagrangian equations of motion, we represent curves implicitly with a new co-dimensional 2 level set description. Our implicit representation admits several redundant mathematical degrees of freedom in both the configuration and the dynamics of the curves, which can be tailored specifically to improve numerical robustness, in contrast to naive approaches for implicit curve dynamics that suffer from overwhelming numerical stability problems. Furthermore, we note how these hidden degrees of freedom perfectly map to a Clebsch representation in fluid dynamics. Motivated by these observations, we introduce untwisted level set functions and non-swirling dynamics which successfully regularize sources of numerical instability, particularly in the twisting modes around curve filaments. A consequence is a novel simulation method which produces stable dynamics for large numbers of interacting vortex filaments and effortlessly handles topological changes and re-connection events. 
    more » « less
  4. We study the emergence of precessing vortex core (PVC) oscillations in a swirling jet experiment. We vary the swirl intensity while keeping the net mass flow rate fixed using a radial-entry swirler with movable blades upstream of the jet exit. The swirl intensity is quantified in terms of a swirl number $$S$$ . Time-resolved velocity measurements in a radial–axial plane anchored at the jet exit for various $$S$$ values are obtained using stereoscopic particle image velocimetry. Spectral proper orthogonal decomposition and spatial cross-spectral analysis reveal the simultaneous emergence of a bubble-type vortex breakdown and a strong helical limit-cycle oscillation in the flow for $$S>S_{c}$$ where $$S_{c}=0.61$$ . The oscillation frequency, $$f_{PVC}$$ , and the square of the flow oscillation amplitudes vary linearly with $$S-S_{c}$$ . A solution for the coherent unsteady field accurate up to $$O(\unicode[STIX]{x1D716}^{3})$$ ( $$\unicode[STIX]{x1D716}\sim O((S-S_{c})^{1/2})$$ ) is determined from the nonlinear Navier–Stokes equations, using the method of multiple scales. We show that onset of bubble type vortex breakdown at $$S_{c}$$ , results in a marginally stable, helical linear global hydrodynamic mode. This results in the stable limit-cycle precession of the breakdown bubble. The variation of $$f_{LC}$$ with $$S-S_{c}$$ is determined from the Stuart–Landau equation associated with the PVC. Reasonable agreement with the corresponding experimental result is observed, despite the highly turbulent nature of the flow in the present experiment. Further, amplitude saturation results from the time-averaged distortion imposed on the flow by the PVC, suggesting that linear stability analysis may predict PVC characteristics for $$S>S_{c}$$ . 
    more » « less
  5. Rotating and swirling turbulence comprises an important class of flows, not only due to the complex physics that occur, but also due to their relevance to many engineering applications, such as combustion, cyclone separation, mixing, etc. In these types of flows, rotation strongly affects the characteristics and structure of turbulence. However, the underlying turbulent flow phenomena are complex and currently not well understood. The axially rotating pipe is an exemplary prototypical model problem that exhibits these complex turbulent flow physics. By examining the complex interaction of turbulent structures within rotating turbulent pipe flow, insight can be gained into the behavior of rotating flows relevant to engineering applications. Direct numerical simulations are conducted at a bulk Reynolds number up to Re_D = 19,000 with rotation numbers ranging from N = 0 to 3. Coherence analysis, including Proper Orthogonal Decomposition and Dynamic Mode Decomposition, are used to identify the relevant (highest energy) modes of the flow. Studying the influence of these modes on turbulent statistics (i.e. mean statistics, Reynolds stresses, turbulent kinetic energy, and turbulent kinetic energy budgets) allow for a deeper understanding of the effects of coherent turbulent flow structures in rotating flows. 
    more » « less