skip to main content

Title: Axisymmetric contour dynamics for buoyant vortex rings
The present work uses a reduced-order model to study the motion of a buoyant vortex ring with non-negligible core size. Buoyancy is considered in both non-Boussinesq and Boussinesq situations using an axisymmetric contour dynamics formulation. The density of the vortex ring differs from that of the ambient fluid, and both densities are constant and conserved. The motion of the ring is calculated by following the boundary of the vortex core, which is also the interface between the two densities. The velocity of the contour comes from a combination of a specific continuous vorticity distribution within its core and a vortex sheet on the core boundary. An evolution equation for the vortex sheet is derived from the Euler equation, which simplifies considerably in the Boussinesq limit. Numerical solutions for the coupled integro-differential equations are obtained. The dynamics of the vortex sheet and the formation of two possible singularities, including singularities in the curvature and the shock-like profile of the vortex sheet strength, are discussed. Three dimensionless groups, the Atwood, Froude and Weber numbers, are introduced to measure the importance of physical effects acting on the motion of a buoyant vortex ring.
Award ID(s):
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate the motion of a thin vortex filament in the presence of buoyancy. The asymptotic model of Moore & Saffman ( Phil. Trans. R. Soc. Lond. A, vol. 272, 1972, pp. 403–429) is extended to take account of buoyancy forces in the force balance on a vortex element. The motion of a buoyant vortex is given by the transverse component of force balance, while the tangential component governs the dynamics of the structure in the core. We show that the local acceleration of axial flow is generated by the external pressure gradient due to gravity. The equations are then solved for vortex rings. An analytic solution for a buoyant vortex ring at a small initial inclination is obtained and asymptotically agrees with the literature.
  2. This paper investigates the steady axisymmetric structure of the cold boundary-layer flow surrounding fire whirls developing over localized fuel sources lying on a horizontal surface. The inviscid swirling motion found outside the boundary layer, driven by the entrainment of the buoyant turbulent plume of hot combustion products that develops above the fire, is described by an irrotational solution, obtained by combining Taylor's self-similar solution for the motion in the axial plane with the azimuthal motion induced by a line vortex of circulation $2 {\rm \pi}\Gamma$ . The development of the boundary layer from a prescribed radial location is determined by numerical integration for different swirl levels, measured by the value of the radial-to-azimuthal velocity ratio $\sigma$ at the initial radial location. As in the case $\sigma =0$ , treated in the seminal boundary-layer analysis of Burggraf et al. ( Phys. Fluids , vol. 14, 1971, pp. 1821–1833), the pressure gradient associated with the centripetal acceleration of the inviscid flow is seen to generate a pronounced radial inflow. Specific attention is given to the terminal shape of the boundary-layer velocity near the axis, which displays a three-layered structure that is described by matched asymptotic expansions. The resulting composite expansion, dependent onmore »the level of ambient swirl through the parameter $\sigma$ , is employed as boundary condition to describe the deflection of the boundary-layer flow near the axis to form a vertical swirl jet. Numerical solutions of the resulting non-slender collision region for different values of $\sigma$ are presented both for inviscid flow and for viscous flow with moderately large values of the controlling Reynolds number $\Gamma /\nu$ . The velocity description provided is useful in mathematical formulations of localized fire-whirl flows, providing consistent boundary conditions accounting for the ambient swirl level.« less
  3. 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 themore »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.« less
  4. We consider the time evolution in two spatial dimensions of a double vorticity layer consisting of two contiguous, infinite material fluid strips, each with uniform but generally differing vorticity, embedded in an otherwise infinite, irrotational, inviscid incompressible fluid. The potential application is to the wake dynamics formed by two boundary layers separating from a splitter plate. A thin-layer approximation is constructed where each layer thickness, measured normal to the common centre curve, is small in comparison with the local radius of curvature of the centre curve. The three-curve equations of contour dynamics that fully describe the double-layer dynamics are expanded in the small thickness parameter. At leading order, closed nonlinear initial-value evolution equations are obtained that describe the motion of the centre curve together with the time and spatial variation of each layer thickness. In the special case where the layer vorticities are equal, these equations reduce to the single-layer equation of Moore ( Stud. Appl. Math. , vol. 58, 1978, pp. 119–140). Analysis of the linear stability of the first-order equations to small-amplitude perturbations shows Kelvin–Helmholtz instability when the far-field fluid velocities on either side of the double layer are unequal. Equal velocities define a circulation-free double vorticity layer,more »for which solution of the initial-value problem using the Laplace transform reveals a double pole in transform space leading to linear algebraic growth in general, but there is a class of interesting initial conditions with no linear growth. This is shown to agree with the long-wavelength limit of the full linearized, three-curve stability equations.« less
  5. Downburst events initialized at various hours during the evening transition (ET) period are simulated to determine the effects of ambient stability on the outflow of downburst winds. The simulations are performed using a pseudo-spectral large eddy simulation model at high resolution to capture both the large-scale flow and turbulence characteristics of downburst winds. First, a simulation of the ET is performed to generate realistic initial and boundary conditions for the subsequent downburst simulations. At each hour in the ET, an ensemble of downburst simulations is initialized separately from the ET simulation in which an elevated cooling source within the model domain generates negatively buoyant air to mimic downburst formation.

    The simulations show that while the stability regime changes, the ensemble mean of the peak wind speed remains fairly constant (between 35 and 38 m s−1) and occurs at the lowest model level for each simulation. However, there is a slight increase in intensity and decrease in the spread of the maximum outflow winds as stability increases as well as an increase in the duration over which these strongest winds persist. This appears to be due to the enhanced maintenance of the ring vortex that results from the low-level temperature inversion, increased ambientmore »shear, and a lack of turbulence within the stable cases. Coherent turbulent kinetic energy and wavelet spectral analysis generally show increased energy in the convective cases and that energy increases across all scales as the downburst passes.« less