Two-dimensional Kelvin-Helmholtz instabilities of multi-component fluids

Hyun Geun Lee, Junseok Kim

Research output: Contribution to journalArticle

21 Citations (Scopus)

Abstract

The Kelvin-Helmholtz instability of multi-component (more than two) incompressible and immiscible fluids is studied numerically using a phase-field model. The instability is governed by the modified Navier-Stokes equations and the multi-component convective Cahn-Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ the Chorin's projection method for the modified Navier-Stokes equations and the recently developed practically unconditionally stable method for the multi-component Cahn-Hilliard equations. Through our model and numerical solution, we investigate the effects of surface tension, density ratio, magnitude of velocity difference, and forcing on the Kelvin-Helmholtz instability of multi-component fluids. It is shown that increasing the surface tension or the density ratio reduces the growth of the Kelvin-Helmholtz instability. And it is also observed that as the initial horizontal velocity difference gets larger, the interface rolls up more. We also found that the billow height reaches its maximum more slowly as the initial wave amplitude gets smaller. And, for the linear growth rate for the Kelvin-Helmholtz instability of two-component fluids, the simulation results agree well with the analytical results. From comparison between the numerical growth rate of two- and three-component fluids, we observe that the inclusion of extra layers can alter the growth rate for the Kelvin-Helmholtz instability. Finally, we simulate the billowing cloud formation which is a classic example of the Kelvin-Helmholtz instability and cannot be seen in binary fluids. With our multi-component method, the details of the real flow (e.g., the asymmetry in the roll-up and the self-interaction of the shear layer) are well captured.

Original languageEnglish
Pages (from-to)77-88
Number of pages12
JournalEuropean Journal of Mechanics, B/Fluids
Volume49
Issue numberPA
DOIs
Publication statusPublished - 2015 Jan 1

Fingerprint

Kelvin-Helmholtz Instability
Kelvin-Helmholtz instability
Fluid
fluids
Cahn-Hilliard Equation
Modified Equations
Surface Tension
Navier-Stokes Equations
Navier-Stokes equation
Immiscible Fluids
interfacial tension
Phase Field Model
Unconditionally Stable
Projection Method
Incompressible Fluid
binary fluids
Forcing
Difference Method
Asymmetry
incompressible fluids

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Mathematical Physics

Cite this

Two-dimensional Kelvin-Helmholtz instabilities of multi-component fluids. / Lee, Hyun Geun; Kim, Junseok.

In: European Journal of Mechanics, B/Fluids, Vol. 49, No. PA, 01.01.2015, p. 77-88.

Research output: Contribution to journalArticle

@article{ddffdc7531fd43d4a8658a1fdc7a47d4,
title = "Two-dimensional Kelvin-Helmholtz instabilities of multi-component fluids",
abstract = "The Kelvin-Helmholtz instability of multi-component (more than two) incompressible and immiscible fluids is studied numerically using a phase-field model. The instability is governed by the modified Navier-Stokes equations and the multi-component convective Cahn-Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ the Chorin's projection method for the modified Navier-Stokes equations and the recently developed practically unconditionally stable method for the multi-component Cahn-Hilliard equations. Through our model and numerical solution, we investigate the effects of surface tension, density ratio, magnitude of velocity difference, and forcing on the Kelvin-Helmholtz instability of multi-component fluids. It is shown that increasing the surface tension or the density ratio reduces the growth of the Kelvin-Helmholtz instability. And it is also observed that as the initial horizontal velocity difference gets larger, the interface rolls up more. We also found that the billow height reaches its maximum more slowly as the initial wave amplitude gets smaller. And, for the linear growth rate for the Kelvin-Helmholtz instability of two-component fluids, the simulation results agree well with the analytical results. From comparison between the numerical growth rate of two- and three-component fluids, we observe that the inclusion of extra layers can alter the growth rate for the Kelvin-Helmholtz instability. Finally, we simulate the billowing cloud formation which is a classic example of the Kelvin-Helmholtz instability and cannot be seen in binary fluids. With our multi-component method, the details of the real flow (e.g., the asymmetry in the roll-up and the self-interaction of the shear layer) are well captured.",
keywords = "Billowing cloud formation, Kelvin-Helmholtz instability, Linear growth rate, Multi-component fluid flows, Phase-field model",
author = "Lee, {Hyun Geun} and Junseok Kim",
year = "2015",
month = "1",
day = "1",
doi = "10.1016/j.euromechflu.2014.08.001",
language = "English",
volume = "49",
pages = "77--88",
journal = "European Journal of Mechanics, B/Fluids",
issn = "0997-7546",
publisher = "Elsevier BV",
number = "PA",

}

TY - JOUR

T1 - Two-dimensional Kelvin-Helmholtz instabilities of multi-component fluids

AU - Lee, Hyun Geun

AU - Kim, Junseok

PY - 2015/1/1

Y1 - 2015/1/1

N2 - The Kelvin-Helmholtz instability of multi-component (more than two) incompressible and immiscible fluids is studied numerically using a phase-field model. The instability is governed by the modified Navier-Stokes equations and the multi-component convective Cahn-Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ the Chorin's projection method for the modified Navier-Stokes equations and the recently developed practically unconditionally stable method for the multi-component Cahn-Hilliard equations. Through our model and numerical solution, we investigate the effects of surface tension, density ratio, magnitude of velocity difference, and forcing on the Kelvin-Helmholtz instability of multi-component fluids. It is shown that increasing the surface tension or the density ratio reduces the growth of the Kelvin-Helmholtz instability. And it is also observed that as the initial horizontal velocity difference gets larger, the interface rolls up more. We also found that the billow height reaches its maximum more slowly as the initial wave amplitude gets smaller. And, for the linear growth rate for the Kelvin-Helmholtz instability of two-component fluids, the simulation results agree well with the analytical results. From comparison between the numerical growth rate of two- and three-component fluids, we observe that the inclusion of extra layers can alter the growth rate for the Kelvin-Helmholtz instability. Finally, we simulate the billowing cloud formation which is a classic example of the Kelvin-Helmholtz instability and cannot be seen in binary fluids. With our multi-component method, the details of the real flow (e.g., the asymmetry in the roll-up and the self-interaction of the shear layer) are well captured.

AB - The Kelvin-Helmholtz instability of multi-component (more than two) incompressible and immiscible fluids is studied numerically using a phase-field model. The instability is governed by the modified Navier-Stokes equations and the multi-component convective Cahn-Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ the Chorin's projection method for the modified Navier-Stokes equations and the recently developed practically unconditionally stable method for the multi-component Cahn-Hilliard equations. Through our model and numerical solution, we investigate the effects of surface tension, density ratio, magnitude of velocity difference, and forcing on the Kelvin-Helmholtz instability of multi-component fluids. It is shown that increasing the surface tension or the density ratio reduces the growth of the Kelvin-Helmholtz instability. And it is also observed that as the initial horizontal velocity difference gets larger, the interface rolls up more. We also found that the billow height reaches its maximum more slowly as the initial wave amplitude gets smaller. And, for the linear growth rate for the Kelvin-Helmholtz instability of two-component fluids, the simulation results agree well with the analytical results. From comparison between the numerical growth rate of two- and three-component fluids, we observe that the inclusion of extra layers can alter the growth rate for the Kelvin-Helmholtz instability. Finally, we simulate the billowing cloud formation which is a classic example of the Kelvin-Helmholtz instability and cannot be seen in binary fluids. With our multi-component method, the details of the real flow (e.g., the asymmetry in the roll-up and the self-interaction of the shear layer) are well captured.

KW - Billowing cloud formation

KW - Kelvin-Helmholtz instability

KW - Linear growth rate

KW - Multi-component fluid flows

KW - Phase-field model

UR - http://www.scopus.com/inward/record.url?scp=84907486465&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84907486465&partnerID=8YFLogxK

U2 - 10.1016/j.euromechflu.2014.08.001

DO - 10.1016/j.euromechflu.2014.08.001

M3 - Article

AN - SCOPUS:84907486465

VL - 49

SP - 77

EP - 88

JO - European Journal of Mechanics, B/Fluids

JF - European Journal of Mechanics, B/Fluids

SN - 0997-7546

IS - PA

ER -