Numerical Simulation of Gravity Current Descending a Slope into a Linearly Stratified Environment

Yakun Guo, Zhiyong Zhang, Bing Shi

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

The accurate prediction of the dilution and motion of the produced denser water (e.g., discharge of concentrated brine generated during solution mining and desalination) is of importance for environmental protection. Boundary conditions and ambient stratification can significantly affect the dilution and motion of gravity currents. In this study, a multiphase model was applied to simulate the gravity current descending a slope into a linearly stratified ambient. The k-ω turbulence model was used to better simulate the near-bed motion. The mathematical model, the initial and boundary conditions, and the details of the numerical scheme are described here. The time-dependent evolution of the gravity current, the flow thickness, and the velocity and density field were simulated for a range of flow parameters. Simulations show that the Kelvin–Helmholtz (K-H) billows are generated at the top of the trailing fluid by the interfacial velocity shear. The K-H instability becomes weaker with the slope distance from the source because of the decrease in interfacial velocity shear along the slope. The ambient stratification restricts and decreases the current head velocity as it descends the slope, which differs from the situation in the homogenous ambient while the head velocity remains in an approximately steady state. Motion of the descending flow into the stratified ambient has two stages: initial acceleration and deceleration at a later stage based on the balance of inertial, buoyancy, and friction forces. When the descending current approaches the initial neutral position at a later stage, it separates from the slope and spreads horizontally into the environment. The simulated results, such as vertical velocity and density profiles and front positions, agree well with the measurements, indicating that the mathematical model can be successfully applied to simulate the effect of the boundary conditions and ambient stratification on the dilution and propagation of gravity currents.
Original languageEnglish
Article number04014061
JournalJournal of Hydraulic Engineering
Volume140
Issue number12
Early online date3 Sep 2014
DOIs
Publication statusPublished - Dec 2014

Fingerprint

Gravitation
gravity
Computer simulation
Dilution
simulation
dilution
stratification
boundary condition
Boundary conditions
Solution mining
solution mining
Mathematical models
dense water
Deceleration
Desalination
Environmental protection
Buoyancy
Turbulence models
desalination
buoyancy

Keywords

  • gravity current
  • numerical models
  • simulation
  • stratification

Cite this

Numerical Simulation of Gravity Current Descending a Slope into a Linearly Stratified Environment. / Guo, Yakun; Zhang, Zhiyong; Shi, Bing.

In: Journal of Hydraulic Engineering, Vol. 140, No. 12, 04014061, 12.2014.

Research output: Contribution to journalArticle

@article{f7cd237f33c746bd8c89bae98d8d8178,
title = "Numerical Simulation of Gravity Current Descending a Slope into a Linearly Stratified Environment",
abstract = "The accurate prediction of the dilution and motion of the produced denser water (e.g., discharge of concentrated brine generated during solution mining and desalination) is of importance for environmental protection. Boundary conditions and ambient stratification can significantly affect the dilution and motion of gravity currents. In this study, a multiphase model was applied to simulate the gravity current descending a slope into a linearly stratified ambient. The k-ω turbulence model was used to better simulate the near-bed motion. The mathematical model, the initial and boundary conditions, and the details of the numerical scheme are described here. The time-dependent evolution of the gravity current, the flow thickness, and the velocity and density field were simulated for a range of flow parameters. Simulations show that the Kelvin–Helmholtz (K-H) billows are generated at the top of the trailing fluid by the interfacial velocity shear. The K-H instability becomes weaker with the slope distance from the source because of the decrease in interfacial velocity shear along the slope. The ambient stratification restricts and decreases the current head velocity as it descends the slope, which differs from the situation in the homogenous ambient while the head velocity remains in an approximately steady state. Motion of the descending flow into the stratified ambient has two stages: initial acceleration and deceleration at a later stage based on the balance of inertial, buoyancy, and friction forces. When the descending current approaches the initial neutral position at a later stage, it separates from the slope and spreads horizontally into the environment. The simulated results, such as vertical velocity and density profiles and front positions, agree well with the measurements, indicating that the mathematical model can be successfully applied to simulate the effect of the boundary conditions and ambient stratification on the dilution and propagation of gravity currents.",
keywords = "gravity current, numerical models, simulation, stratification",
author = "Yakun Guo and Zhiyong Zhang and Bing Shi",
year = "2014",
month = "12",
doi = "10.1061/(ASCE)HY.1943-7900.0000936",
language = "English",
volume = "140",
journal = "Journal of Hydraulic Engineering",
issn = "0733-9429",
publisher = "American Society of Civil Engineers (ASCE)",
number = "12",

}

TY - JOUR

T1 - Numerical Simulation of Gravity Current Descending a Slope into a Linearly Stratified Environment

AU - Guo, Yakun

AU - Zhang, Zhiyong

AU - Shi, Bing

PY - 2014/12

Y1 - 2014/12

N2 - The accurate prediction of the dilution and motion of the produced denser water (e.g., discharge of concentrated brine generated during solution mining and desalination) is of importance for environmental protection. Boundary conditions and ambient stratification can significantly affect the dilution and motion of gravity currents. In this study, a multiphase model was applied to simulate the gravity current descending a slope into a linearly stratified ambient. The k-ω turbulence model was used to better simulate the near-bed motion. The mathematical model, the initial and boundary conditions, and the details of the numerical scheme are described here. The time-dependent evolution of the gravity current, the flow thickness, and the velocity and density field were simulated for a range of flow parameters. Simulations show that the Kelvin–Helmholtz (K-H) billows are generated at the top of the trailing fluid by the interfacial velocity shear. The K-H instability becomes weaker with the slope distance from the source because of the decrease in interfacial velocity shear along the slope. The ambient stratification restricts and decreases the current head velocity as it descends the slope, which differs from the situation in the homogenous ambient while the head velocity remains in an approximately steady state. Motion of the descending flow into the stratified ambient has two stages: initial acceleration and deceleration at a later stage based on the balance of inertial, buoyancy, and friction forces. When the descending current approaches the initial neutral position at a later stage, it separates from the slope and spreads horizontally into the environment. The simulated results, such as vertical velocity and density profiles and front positions, agree well with the measurements, indicating that the mathematical model can be successfully applied to simulate the effect of the boundary conditions and ambient stratification on the dilution and propagation of gravity currents.

AB - The accurate prediction of the dilution and motion of the produced denser water (e.g., discharge of concentrated brine generated during solution mining and desalination) is of importance for environmental protection. Boundary conditions and ambient stratification can significantly affect the dilution and motion of gravity currents. In this study, a multiphase model was applied to simulate the gravity current descending a slope into a linearly stratified ambient. The k-ω turbulence model was used to better simulate the near-bed motion. The mathematical model, the initial and boundary conditions, and the details of the numerical scheme are described here. The time-dependent evolution of the gravity current, the flow thickness, and the velocity and density field were simulated for a range of flow parameters. Simulations show that the Kelvin–Helmholtz (K-H) billows are generated at the top of the trailing fluid by the interfacial velocity shear. The K-H instability becomes weaker with the slope distance from the source because of the decrease in interfacial velocity shear along the slope. The ambient stratification restricts and decreases the current head velocity as it descends the slope, which differs from the situation in the homogenous ambient while the head velocity remains in an approximately steady state. Motion of the descending flow into the stratified ambient has two stages: initial acceleration and deceleration at a later stage based on the balance of inertial, buoyancy, and friction forces. When the descending current approaches the initial neutral position at a later stage, it separates from the slope and spreads horizontally into the environment. The simulated results, such as vertical velocity and density profiles and front positions, agree well with the measurements, indicating that the mathematical model can be successfully applied to simulate the effect of the boundary conditions and ambient stratification on the dilution and propagation of gravity currents.

KW - gravity current

KW - numerical models

KW - simulation

KW - stratification

U2 - 10.1061/(ASCE)HY.1943-7900.0000936

DO - 10.1061/(ASCE)HY.1943-7900.0000936

M3 - Article

VL - 140

JO - Journal of Hydraulic Engineering

JF - Journal of Hydraulic Engineering

SN - 0733-9429

IS - 12

M1 - 04014061

ER -