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. This work describes and tests the calibration process of the chemical-diffusive model (CDM) for the simulation of non-premixed diffusion flames. The CDM is an alternative, simplified approach for incorporating the effects of combustion in a fluid simulation, based on the ideas of regulating the rate of energy release such that the properties of combustion waves (e.g. flames and detonations) are reproduced. Past implementations of the CDM have considered single-stoichiometry fuel-air mixtures or mixtures with variable stoichiom- etry but with premixed modes of combustion. In this work, the CDM is tested and shown to work for non-premixed, low-Mach-number flames (i.e., diffusion flames) by incorporat- ing it into a numerical model which solves the reactive and compressible Navier-Stokes equations with the barely implicit correction (BIC) algorithm, which removes the acoustic limit on the integration time-step size. Simulations of one-dimensional premixed laminar flames reproduce the required premixed laminar flame speed, thickness, and temperature. A two-dimensional, steady-state, laminar coflow diffusion flame is computed, and the result demonstrates the ability of the algorithm to compute a non-premixed flame. Lastly, a two- dimensional simulation of two opposing jets of fuel and air show that the CDM approach can compute the structure of a counter-flow diffusion flame. 
    more » « less