Record Details

A Numerical Study On Absolute Instability Of Low Density Jets

Electronic Theses of Indian Institute of Science

View Archive Info
 
 
Field Value
 
Title A Numerical Study On Absolute Instability Of Low Density Jets
 
Creator Chakravorty, Saugata
 
Subject Aerodynamics
Jets - Fluid Dynamics
Absolute Instability
Jets - Numerical Simulations
Jet Instability
Low Density Jets
Turbulence
Density Jet
 
Description A spectacular instability has been observed in low density round jets when the density ratio of jet fluid to ambient fluid falls below a threshold of approximately 0.6. This phenomenon has been observed in non-buoyant jets of helium in air, heated air jets and heated buoyant jets. The oscillation of the flow near the nozzle is extremely regular and periodic and consists of ring vortices. Even the smaller scale structures that appear downstream exhibit similar regularity. A theory for predicting the onset of this oscillation is based on finding regions of absolute instability from linear stability analysis of parallel flow. However, experiments suggest that the theory is at least incomplete and fortuitous as the oscillation is not a linear process. The present work is to observe and understand the process of regeneration of these oscillations by conducting numerical simulations. Here, two-dimensional, plane jets were simulated because they undergo a qualitatively similar process.

A spatial and temporal picture of a heated jet has been obtained numerically. A perturbation expansion was used to obtain a system of conservation laws for compressible flows which is valid for low Mach numbers. The low Mach number approximation removes the high frequency acoustic waves from the flow field. This enables a larger time step to be taken without making the calculation unstable. To ensure that all the scales of motion are properly resolved, calculations were done at a low Reynolds number. The governing equations were discretized in space using second-order finite difference formulas on a staggered grid. Velocity fields were advanced using a second-order Adams-Bashforth explicit scheme and then corrected by solving for pressure such that continuity is satisfied at every time step. The Poisson problem for pressure requires the time derivative of the density which was approximated by a third-order backward difference formula. Gauss-Siedel iteration was used to find the pressure.

Several numerical tests were conducted prior to simulations of variable density jets to check the stability and accuracy of the code. Two dimensional driven cavity flow calculations were done as a first test. Then a calculation of a forced, spatially developing, incompressible, plane mixing layer was done to check the time accuracy of the code. After obtaining satisfactory performance of the code for the different test cases, two-dimensional, variable density jets were simulated.

Since the plane jet extends ad infinitum in the streamwise direction, a sufficiently large domain was used to capture all the relevant physics in the downstream regions of the jet. An advective boundary condition was imposed at the exit plane. Rigid, slipwall conditions were employed to prescribe lateral boundary conditions.

A 2-D, incompressible plane jet was simulated first. The jet profile was approximated by two hyperbolic tangent shear layers. The most unstable mode of the inviscid shear layer for this profile, along with its first and second harmonics, was imposed on the velocity profile at the inlet plane. The amplitude of oscillation of the harmonics was chosen so as to provide sufficient energy in the perturbation to accelerate the growth of the layer. No explicit phase lag was introduced in the perturbation. The flow was allowed to develop long enough to wash out the effect of the initial condition. The results obtained for this case indicate that experimentally realized phenomena such as vortex pairing were captured in this simulation. Furthermore, to check the convective nature of instability of the incompressible jet, the forcing at the inlet plane was turned off. The disturbances were gradually convected downstream, out of the computational domain.

Next, two-dimensional heated, non-buoyant jets were studied numerically. The effects of the ratio of jet density to ambient density S, the velocity ratio R, and jet width W, on the near field behavior of an initial laminar jet and the regeneration mechanism of the self-sustaining vortices were explored. The theory based on domain of absolute/convective instability identifies these three parameters. No initial perturbation was necessary to start roll-up of the shear layer. For certain choices, e.g., S= 0.75, R = 20, W =10.5, self-sustaining oscillations appeared spontaneously, and these cycles repeated for very long simulation intervals. Waviness on the jet shear layers grow and roll-up into vortices as in constant density shear layers. But unlike the incompressible plane jet, these vortices grow much larger and mixes more with the surrounding fluid. As these vortices evolve, packets of fluid break away as trailing legs similar to side jet expulsions observed in round jets and plumes. The growing vortices disturb the upstream shear layer. Consistently with linear theory, which predicts absolute instability for these parameters, these disturbances are able to grow and roll up. If these disturbances travelled faster than the downstream vortices, it would not be possible for the cycle to repeat. With sufficient shear between the co-flowing streams (R not too small), the entire regeneration process was found to begin from roughly the same streamwise location. Furthermore, it is the symmetric, varicose mode which occurs. At a slightly larger density ratio (S = 0.8, R = 10), self-sustaining oscillations appeared, but each new cycle began slightly farther downstream. It seems likely that these values are close to the boundary in parameter space between self-sustained oscillatory and convectively unstable behaviors. Jet width also influences the selection of these two behaviors. When jet width was reduced, W = 6, even for S = 0.75,R = 20, each new cycle began to shift downstream. For larger jet width (W = 12.3), self-sustaining oscillations occur but the response is now as an asymmetric sinuous mode after a short initial varicose mode.

The detailed processes that have now been revealed in plane jets should serve as guidelines for the study of such processes in the technologically more important round jets.
 
Publisher Indian Institute of Science
 
Contributor Mathew, Joseph
 
Date 2006-07-20T05:52:21Z
2006-07-20T05:52:21Z
2006-07-20T05:52:21Z
2000-05
 
Type Electronic Thesis and Dissertation
 
Format 3824906 bytes
application/pdf
 
Identifier http://etd.iisc.ernet.in/handle/2005/227
null
 
Language en
 
Rights I grant Indian Institute of Science the right to archive and to make available my thesis or dissertation in whole or in part in all forms of media, now hereafter known. I retain all proprietary rights, such as patent rights. I also retain the right to use in future works (such as articles or books) all or part of this thesis or dissertation.