Direct numerical simulations of planar and cylindrical density currents

Mariano I. Cantero, S. Balachandar, Marcelo H. García, James P. Ferry

Research output: Contribution to journalArticle

Abstract

The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof(Gr) numbers are investigated (10 5 and 1.5×10 6) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial threedimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.

Original languageEnglish (US)
Pages (from-to)923-930
Number of pages8
JournalJournal of Applied Mechanics, Transactions ASME
Volume73
Issue number6
DOIs
StatePublished - Nov 1 2006

Fingerprint

Direct numerical simulation
direct numerical simulation
Current density
Flow structure
current density
slumping
Grashof number
Fluids
Vortex shedding
Navier Stokes equations
Numerical models
Vortex flow
Boussinesq approximation
three dimensional flow
vortex shedding
fluids
configurations
lobes
Navier-Stokes equation
slip

ASJC Scopus subject areas

  • Condensed Matter Physics
  • Mechanics of Materials
  • Mechanical Engineering

Cite this

Direct numerical simulations of planar and cylindrical density currents. / Cantero, Mariano I.; Balachandar, S.; García, Marcelo H.; Ferry, James P.

In: Journal of Applied Mechanics, Transactions ASME, Vol. 73, No. 6, 01.11.2006, p. 923-930.

Research output: Contribution to journalArticle

Cantero, Mariano I. ; Balachandar, S. ; García, Marcelo H. ; Ferry, James P. / Direct numerical simulations of planar and cylindrical density currents. In: Journal of Applied Mechanics, Transactions ASME. 2006 ; Vol. 73, No. 6. pp. 923-930.
@article{4c12421e90a742b7bf55da610ced192f,
title = "Direct numerical simulations of planar and cylindrical density currents",
abstract = "The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof(Gr) numbers are investigated (10 5 and 1.5×10 6) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial threedimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.",
author = "Cantero, {Mariano I.} and S. Balachandar and Garc{\'i}a, {Marcelo H.} and Ferry, {James P.}",
year = "2006",
month = "11",
day = "1",
doi = "10.1115/1.2173671",
language = "English (US)",
volume = "73",
pages = "923--930",
journal = "Journal of Applied Mechanics, Transactions ASME",
issn = "0021-8936",
publisher = "American Society of Mechanical Engineers(ASME)",
number = "6",

}

TY - JOUR

T1 - Direct numerical simulations of planar and cylindrical density currents

AU - Cantero, Mariano I.

AU - Balachandar, S.

AU - García, Marcelo H.

AU - Ferry, James P.

PY - 2006/11/1

Y1 - 2006/11/1

N2 - The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof(Gr) numbers are investigated (10 5 and 1.5×10 6) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial threedimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.

AB - The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof(Gr) numbers are investigated (10 5 and 1.5×10 6) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial threedimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.

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

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

U2 - 10.1115/1.2173671

DO - 10.1115/1.2173671

M3 - Article

AN - SCOPUS:33751113889

VL - 73

SP - 923

EP - 930

JO - Journal of Applied Mechanics, Transactions ASME

JF - Journal of Applied Mechanics, Transactions ASME

SN - 0021-8936

IS - 6

ER -