Hostname: page-component-77f85d65b8-7lfxl Total loading time: 0 Render date: 2026-04-18T03:22:39.856Z Has data issue: false hasContentIssue false

Combined Rayleigh–Taylor and Kelvin–Helmholtz instabilities on an annular liquid sheet

Published online by Cambridge University Press:  22 December 2016

M. Vadivukkarasan
Affiliation:
Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai, 600036, India
Mahesh V. Panchagnula*
Affiliation:
Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai, 600036, India
*
Email address for correspondence: mvp@iitm.ac.in

Abstract

This paper describes the three-dimensional destabilization characteristics of an annular liquid sheet when subjected to the combined action of Rayleigh–Taylor (RT) and Kelvin–Helmholtz (KH) instability mechanisms. The stability characteristics are studied using temporal linear stability analysis and by assuming that the fluids are incompressible, immiscible and inviscid. Surface tension is also taken into account at both the interfaces. Linearized equations governing the growth of instability amplitude have been derived. These equations involve time-varying coefficients and have been analysed using two approaches – direct numerical time integration and frozen-flow approximation. From the direct numerical time integration, we show that the time-varying coefficients evolve on a slow time scale in comparison with the amplitude growth. Therefore, we justify the use of the frozen-flow approximation and derive a closed-form dispersion relation from the appropriate governing equations and boundary conditions. The effect of flow conditions and fluid properties is investigated by introducing dimensionless numbers such as Bond number ( $Bo$ ), inner and outer Weber numbers ( $We_{i}$ , $We_{o}$ ) and inner and outer density ratios ( $Q_{i}$ , $Q_{o}$ ). We show that four instability modes are possible – Taylor, sinuous, flute and helical. It is observed that the choice of instability mode is influenced by a combination of both $Bo$ as well as $We_{i}$ and $We_{o}$ . However, the instability length scale calculated from the most unstable wavenumbers is primarily a function of $Bo$ . We show a regime map in the $Bo,We_{i},We_{o}$ parameter space to identify regions where the system is susceptible to three-dimensional helical modes. Finally, we show an optimal partitioning of a given total energy ( $\unicode[STIX]{x1D701}$ ) into acceleration-induced and shear-induced instability mechanisms in order to achieve a minimum instability length scale ( ${\mathcal{L}}_{m}^{\ast }$ ). We show that it is beneficial to introduce at least 90 % of the total energy into acceleration induced RT instability mechanism. In addition, we show that when the RT mechanism is invoked to destabilize an annular liquid sheet, ${\mathcal{L}}_{m}^{\ast }\sim \unicode[STIX]{x1D701}^{-3/5}$ .

JFM classification

Information

Type
Papers
Copyright
© 2016 Cambridge University Press 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

Article purchase

Temporarily unavailable