Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations

Tamay M Ozgokmen, Traian Iliescu, Paul F. Fischer

Research output: Contribution to journalArticle

26 Citations (Scopus)

Abstract

Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching Re = O (106), or even higher. In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be O (1000) times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES. This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3 × 103 ≤ Re ≤ 3 × 104, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re = 105 and 106, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67% when Re is increased from Re = 103 to Re = 106, within the computation period, with the most significant increase taking place from Re = 3 × 103 to Re = 105.

Original languageEnglish (US)
Pages (from-to)190-206
Number of pages17
JournalOcean Modelling
Volume30
Issue number2-3
DOIs
StatePublished - 2009

Fingerprint

large eddy simulation
Large eddy simulation
Reynolds number
Direct numerical simulation
Potential energy
potential energy
Turbulent flow
simulation
turbulent mixing
Deconvolution
deconvolution
turbulent flow
water mass
eddy
boundary condition
viscosity
Boundary conditions
Viscosity
salinity

Keywords

  • DNS
  • LES
  • Mixing
  • Reynolds number

ASJC Scopus subject areas

  • Atmospheric Science
  • Oceanography
  • Geotechnical Engineering and Engineering Geology
  • Computer Science (miscellaneous)

Cite this

Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations. / Ozgokmen, Tamay M; Iliescu, Traian; Fischer, Paul F.

In: Ocean Modelling, Vol. 30, No. 2-3, 2009, p. 190-206.

Research output: Contribution to journalArticle

@article{33037eccdfd741fcae6f7b91cb7cddc9,
title = "Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations",
abstract = "Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching Re = O (106), or even higher. In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be O (1000) times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES. This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3 × 103 ≤ Re ≤ 3 × 104, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re = 105 and 106, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67{\%} when Re is increased from Re = 103 to Re = 106, within the computation period, with the most significant increase taking place from Re = 3 × 103 to Re = 105.",
keywords = "DNS, LES, Mixing, Reynolds number",
author = "Ozgokmen, {Tamay M} and Traian Iliescu and Fischer, {Paul F.}",
year = "2009",
doi = "10.1016/j.ocemod.2009.06.013",
language = "English (US)",
volume = "30",
pages = "190--206",
journal = "Ocean Modelling",
issn = "1463-5003",
publisher = "Elsevier BV",
number = "2-3",

}

TY - JOUR

T1 - Reynolds number dependence of mixing in a lock-exchange system from direct numerical and large eddy simulations

AU - Ozgokmen, Tamay M

AU - Iliescu, Traian

AU - Fischer, Paul F.

PY - 2009

Y1 - 2009

N2 - Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching Re = O (106), or even higher. In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be O (1000) times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES. This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3 × 103 ≤ Re ≤ 3 × 104, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re = 105 and 106, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67% when Re is increased from Re = 103 to Re = 106, within the computation period, with the most significant increase taking place from Re = 3 × 103 to Re = 105.

AB - Turbulent mixing of water masses of different temperatures and salinities is an important process for both coastal and large-scale ocean circulation. It is, however, difficult to capture computationally. One of the reasons is that mixing in the ocean occurs at a wide range of complexity, with the Reynolds number reaching Re = O (106), or even higher. In this study, we continue to investigate whether large eddy simulation (LES) can be a reliable computational tool for stratified mixing in turbulent oceanic flows. LES is attractive because it can be O (1000) times faster than a direct numerical simulation (DNS) of stratified mixing in turbulent flows. Before using the LES methodology to compute mixing in realistic oceanic flows, however, a careful assessment of the LES sensitivity with respect to Re needs to be performed first. The main objectives of this study are: (i) to investigate the performance of different LES models at high Re, such as those encountered in oceanic flows; and (ii) to study how mixing varies as a function of Re. To this end, as a benchmark we use the lock-exchange problem, which is described by unambigous and simple initial and boundary conditions. The background potential energy, which accurately quantifies irreversible mixing in an enclosed system, is used as the main criterion in a posteriori testing of LES. This study has two main achievements. The first is that we investigate the accuracy of six combinations of two different classes of LES models, namely eddy-viscosity and approximate deconvolution types, for 3 × 103 ≤ Re ≤ 3 × 104, for which DNS data is computed. We find that all LES models almost always provide significantly more accurate results than cases without LES models. Nevertheless, no single LES model that is persistently superior to others over this Re range could be identified. Then, an ensemble of the four best performing LES models is selected in order to estimate mixing taking place in this system at Re = 105 and 106, for which DNS is presently not feasible. Thus the second achievement of this study is to quantify mixing taking place in this system over an Re range that changes by three orders of magnitude. We find that the background potential energy increases by about 67% when Re is increased from Re = 103 to Re = 106, within the computation period, with the most significant increase taking place from Re = 3 × 103 to Re = 105.

KW - DNS

KW - LES

KW - Mixing

KW - Reynolds number

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

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

U2 - 10.1016/j.ocemod.2009.06.013

DO - 10.1016/j.ocemod.2009.06.013

M3 - Article

AN - SCOPUS:68749099676

VL - 30

SP - 190

EP - 206

JO - Ocean Modelling

JF - Ocean Modelling

SN - 1463-5003

IS - 2-3

ER -