Hostname: page-component-89b8bd64d-x2lbr Total loading time: 0 Render date: 2026-05-08T18:52:28.899Z Has data issue: false hasContentIssue false

Phase-field model for the Rayleigh–Taylor instability of immiscible fluids

Published online by Cambridge University Press:  10 March 2009

ANTONIO CELANI
Affiliation:
Institut Pasteur, CNRS, URA 2171, 25 Rue du docteur Roux, 75015 Paris, France
ANDREA MAZZINO
Affiliation:
Department of Physics – University of Genova, and CNISM & INFN – Genova Section, via Dodecaneso 33, 16146 Genova, Italy
PAOLO MURATORE-GINANNESCHI
Affiliation:
Department of Mathematics and Statistics – University of Helsinki, PO Box 4, 00014 Helsinki, Finland
LARA VOZELLA*
Affiliation:
Department of Physics – University of Genova, and CNISM & INFN – Genova Section, via Dodecaneso 33, 16146 Genova, Italy Department of Mathematics and Statistics – University of Helsinki, PO Box 4, 00014 Helsinki, Finland
*
Email address for correspondence: vozella@.sica.unige.it

Abstract

The Rayleigh–Taylor instability of two immiscible fluids in the limit of small Atwood numbers is studied by means of a phase-field description. In this method, the sharp fluid interface is replaced by a thin, yet finite, transition layer where the interfacial forces vary smoothly. This is achieved by introducing an order parameter (the phase-field) continuously varying across the interfacial layers and uniform in the bulk region. The phase-field model obeys a Cahn–Hilliard equation and is two-way coupled to the standard Navier–Stokes equations. Starting from this system of equations we have first performed a linear analysis from which we have analytically rederived the known gravity–capillary dispersion relation in the limit of vanishing mixing energy density and capillary width. We have performed numerical simulations and identified a region of parameters in which the known properties of the linear phase (both stable and unstable) are reproduced in a very accurate way. This has been done both in the case of negligible viscosity and in the case of non-zero viscosity. In the latter situation, only upper and lower bounds for the perturbation growth rate are known. Finally, we have also investigated the weakly nonlinear stage of the perturbation evolution and identified a regime characterized by a constant terminal velocity of bubbles/spikes. The measured value of the terminal velocity is in agreement with available theoretical prediction. The phase-field approach thus appears to be a valuable technique for the dynamical description of the stages where hydrodynamic turbulence and wave-turbulence come into play.

Information

Type
Papers
Copyright
Copyright © Cambridge University Press 2009

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