0
Research Papers: Flows in Complex Systems

# Combined Effect of Viscosity, Surface Tension and Compressibility on Rayleigh-Taylor Bubble Growth Between Two Fluids

[+] Author and Article Information
Sourav Roy

Department of Instrumentation Science,
Kolkata-700032, India
e-mail: phy.sou82@gmail.com

L. K. Mandal

Department of Physics,
Budge Budge Institute of Technology,
Kolkata-137, India
e-mail: labakanta@gmail.com

Manoranjan Khan

Department of Instrumentation Science,
Kolkata-700032, India
e-mail: mkhan_ju@yahoo.com

M. R. Gupta

Department of Instrumentation Science,
Kolkata-700032, India

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received January 21, 2013; final manuscript received May 9, 2014; published online July 9, 2014. Assoc. Editor: Stuart Dalziel.

J. Fluids Eng 136(9), 091101 (Jul 09, 2014) (6 pages) Paper No: FE-13-1033; doi: 10.1115/1.4027655 History: Received January 21, 2013; Revised May 09, 2014

## Abstract

The combined effect of viscosity, surface tension, and the compressibility on the nonlinear growth rate of Rayleigh-Taylor (RT) instability has been investigated. For the incompressible case, it is seen that both viscosity and surface tension have a retarding effect on RT bubble growth for the interface perturbation wave number having a value less than three times of a critical value ($kc=(ρh-ρl)g/T$, $T$ is the surface tension). For the value of wave number greater than three times of the critical value, the RT induced unstable interface is stabilized through damped nonlinear oscillation. In the absence of surface tension and viscosity, the compressibility has both a stabilizing and destabilizing effect on RTI bubble growth. The presence of surface tension and viscosity reduces the growth rate. Above a certain wave number, the perturbed interface exhibits damped oscillation. The damping factor increases with increasing kinematic viscosity of the heavier fluid and the saturation value of the damped oscillation depends on the surface tension of the perturbed fluid interface and interface perturbation wave number. An approximate expression for asymptotic bubble velocity considering only the lighter fluid as a compressible one is presented here. The numerical results describing the dynamics of the bubble are represented in diagrams.

<>

## References

Rayleigh, L., 1883, “Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density,” Proc. London Math. Soc., 14, pp.170–177.
Richtmyer, R. D., 1960, “Taylor Instability in a Shock Acceleration of Compressible Fluids,” Commun. Pure Appl. Math., 13, pp. 297–319.
Bernstein, I. B., and Book, D. L., 1983, “Effect of Compressibility on Rayleigh-Taylor Instability,” Phys. Fluids, 26, p. 453.
Yang, Y., and Zhang, Q., 1993, “General Properties of a Multilayer Fluid System,” Phys. Fluids A, 5, p. 1167.
Sharp, D. H., 1984, “An Overview of Rayleigh-Taylor Instability,” Physica D, 12, pp. 3–10.
Plesset, M. S., and Hsieh, D., 1964, “General Analysis of the Stability of Superposed Fluids,” Phys. Fluids, 7, p. 1099.
Baker, L., 1983, “Compressible Rayleigh-Taylor Instability,” Phys. Fluids, 26, p. 950.
Chandrasekhar, S., 1981, Hydrodynamic and Hydromagnetic Stability, Dover Publications, New York, pp. 428–452.
Chhajlani, R. K., and Vaghela, D. S., 1989, “Rayleigh-Taylor Instability of Ionized Viscous Fluids With FLR-Corrections and Surface-Tension,” Astrophys. Space Sci., 155, pp. 257–269.
Mikaelian, K. O., 1993, “Effect of Viscosity on Rayleigh-Taylor and Richtmyer-Meshkov Instability,” Phys. Rev. E, 47, p. 375.
Bhatia, P. K., 1974, “Rayleigh-Taylor Instability of a Viscous Compressible Plasma of Variable Density,” Astrophys. Space Sci., 26, pp. 319–325.
Carlès, P., and Popinet, S., 2002, “The Effect of Viscosity, Surface Tension and Non-Linearity on Richtmyer Meshkov Instability,” Eur. J. Mech. B/Fluids, 21, pp. 511–526.
Sohn, S. I., 2009, “Effects of Surface Tension and Viscosity on the Growth Rates of Rayleigh-Taylor and Richtmyer-Meshkov Instabilities,” Phys. Rev. E, 80, p. 055302.
Roy, S., Gupta, M. R., Khan, M., Pant, H. C., and Srivastava, M. K., 2010, “Effect of Surface Tension on the Rayleigh-Taylor and Richtmyer-Meshkov Instability Induced Nonlinear Structure at Two Fluid Interface and Their Stabilization,” J. Phys. Conf. Ser., 208, p. 012083.
Layzer, D., 1955, “On the Instability of Superposed Fluids in a Gravitational Field,” Astrophys. J., 122, pp. 1–12.
Goncharav, V. N., 2002, “Analytical Model of Nonlinear, Single-Mode, Classical Rayleigh-Taylor Instability at Arbitrary Atwood Numbers,” Phys. Rev. Lett., 88, p. 134502. [PubMed]
Hecht, J., Alon, U., and Shvarts, D., 1994, “Potential Flow Models of Rayleigh-Taylor and Richtmyer-Meshkov Bubble Fronts,” Phys. Fluids, 6, p. 4019.
Drake, R. P., 2006, High Energy Density Physics, Springer, Berlin, p. 175.
Gupta, M. R., Roy, S., Khan, M., Pant, H. C., Sarkar, S., and Srivastava, M. K., 2009, “Effect of Compressibility on the Rayleigh-Taylor and Richtmyer-Meshkov Instability Induced Nonlinear Structure at Two Fluid Interface,” Phys. Plasmas, 16, p. 032303.
Zhang, Q., 1998, “Analytical Solutions of Layzer-Type Approach to Unstable Interfacial Fluid Mixing,” Phys. Rev. Lett., 81, pp. 3391–3394.

## Figures

Fig. 1

Solution of the set of Eqs. (30) neglecting the last term of the third equation for the bubble considering initial values ξ1 = 0.1,ξ2 = -0.1,ξ3 = 0.1 and r = 1.545, s = 1. The solid line represents the solution for χ = 0, k2/kc2 = 0; the dashed line represents that for χ = 0, k2/kc2 = 2, and the dotted line represents that for χ = 0.3, k2/kc2 = 2.

Fig. 2

Solution of the set of Eqs. (30) neglecting the last term of the third equation for bubble considering initial values ξ1 = 0.1,ξ2 = -0.1,ξ3 = 0.1 and r = 1.545, s = 1. The solid line represents the solution for χ = 0, k2/kc2 = 0; the dashed line represents that for χ = 0.3, k2/kc2 = 5; the dotted line represents that for χ = 0.7, k2/kc2 = 5; and the dashed–dotted line represents that for χ = 0.7, k2/kc2 = 5.

Fig. 3

Solution of the set of Eqs. (30) for bubble considering initial values ξ1 = 0.1,ξ2 = -0.1,ξ3 = 0.1,c =  100 and r = 1.545, s = 1. The solid line represents the solution for ρh1/ρh0* = 0, ρl1/ρl0* = 0, χ = 0, k2/kc2 = 0; the dashed line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0, k2/kc2 = 0; the dotted line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0, k2/kc2 = 1; the dashed–dotted line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0.3, k2/kc2 = 0 and the dashed–dotted–dotted line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0.3, k2/kc2 = 1.

Fig. 4

Solution of the set of Eqs. (30) for the bubble considering initial values ξ1 = 0.1,ξ2 = -0.1,ξ3 = 0.1,c = 100 and r = 1.545, s = 1. The solid line represents the solution for ρh1/ρh0* = 0, ρl1/ρl0* = 0, χ = 0, k2/kc2 = 0; the dashed line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0.3, k2/kc2 = 6; the dotted line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0.7, k2/kc2 = 6; and the dashed–dotted line represents that for ρh1/ρh0* = -0.002, ρl1/ρl0* = -0.002, χ = 0.7, k2/kc2 = 7.

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections