Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain

Tamay M Ozgokmen, Traian Iliescu, Paul F. Fischer, Ashwanth Srinivasan, Jinqiao Duan

Research output: Contribution to journalArticle

53 Citations (Scopus)

Abstract

Mixing in both coastal and deep ocean emerges as one of the important processes that determines the transport of pollutants, sediments and biological species, as well as the details of the global thermohaline circulation. Both the observations, due to their lack in space and time resolution, and most coastal and general circulation models due to inadequate physics, can only provide partial information about oceanic mixing processes. A new class of nonhydrostatic models supplemented with physically based subgrid-scale (SGS) closures, or so-called large eddy simulation (LES), is put forth as another tool of investigation to complement observational and large-scale modeling efforts. However, SGS models have been developed primarily for homogeneous, isotropic flows. Here, four SGS models based on Smagorinsky eddy viscosity and diffusivity are tested for stratified flows in the context of 2D dam-break problem in a rectangular enclosed domain. This idealized testbed leads to a number of simplifications about the initial conditions, boundary conditions and geometry, while exhibiting the dynamically complex characteristics of stratified flows involving the interaction of shear-induced mixing and internal waves. Direct numerical simulations (DNS) at high resolutions are taken as benchmark solutions. Under-resolved simulations without SGS terms (so-called DNS*) are used to quantify the impact of SGS stresses. The performance of LES is assessed by using the time evolution of the volume fraction of intermediate density water masses generated by mixing. The simulations are conducted using a nonhydrostatic high-order spectral element model Nek5000 developed to exhibit minimal numerical dissipation and dispersion errors, which is advantageous to quantify accurately the impact of SGS stresses. It is found that all tested SGS models lead to improved results with respect to those from DNS*. Also, SGS models allow for simulations with coarse resolutions that blow up in DNS* due to lack of adequate dissipation where needed. The SGS model in which the vertical eddy diffusion is modulated via a function that depends on the Richardson number Ri shows the most faithful reproduction of mixed water masses at all resolutions tested. The sensitivity of the results to the tunable parameter of the SGS model, to changes in the Ri-dependent function and resolution of the turbulent overturning scales is shown.

Original languageEnglish (US)
Pages (from-to)106-140
Number of pages35
JournalOcean Modelling
Volume16
Issue number1-2
DOIs
StatePublished - 2007

Fingerprint

large eddy simulation
Large eddy simulation
Dams
dam
Direct numerical simulation
simulation
stratified flow
water mass
dissipation
eddy
thermohaline circulation
Richardson number
internal wave
diffusivity
Testbeds
general circulation model
boundary condition
physics
viscosity
Water

ASJC Scopus subject areas

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

Cite this

Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain. / Ozgokmen, Tamay M; Iliescu, Traian; Fischer, Paul F.; Srinivasan, Ashwanth; Duan, Jinqiao.

In: Ocean Modelling, Vol. 16, No. 1-2, 2007, p. 106-140.

Research output: Contribution to journalArticle

Ozgokmen, Tamay M ; Iliescu, Traian ; Fischer, Paul F. ; Srinivasan, Ashwanth ; Duan, Jinqiao. / Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain. In: Ocean Modelling. 2007 ; Vol. 16, No. 1-2. pp. 106-140.
@article{852dbfce2f004080aab014bbd038f2f6,
title = "Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain",
abstract = "Mixing in both coastal and deep ocean emerges as one of the important processes that determines the transport of pollutants, sediments and biological species, as well as the details of the global thermohaline circulation. Both the observations, due to their lack in space and time resolution, and most coastal and general circulation models due to inadequate physics, can only provide partial information about oceanic mixing processes. A new class of nonhydrostatic models supplemented with physically based subgrid-scale (SGS) closures, or so-called large eddy simulation (LES), is put forth as another tool of investigation to complement observational and large-scale modeling efforts. However, SGS models have been developed primarily for homogeneous, isotropic flows. Here, four SGS models based on Smagorinsky eddy viscosity and diffusivity are tested for stratified flows in the context of 2D dam-break problem in a rectangular enclosed domain. This idealized testbed leads to a number of simplifications about the initial conditions, boundary conditions and geometry, while exhibiting the dynamically complex characteristics of stratified flows involving the interaction of shear-induced mixing and internal waves. Direct numerical simulations (DNS) at high resolutions are taken as benchmark solutions. Under-resolved simulations without SGS terms (so-called DNS*) are used to quantify the impact of SGS stresses. The performance of LES is assessed by using the time evolution of the volume fraction of intermediate density water masses generated by mixing. The simulations are conducted using a nonhydrostatic high-order spectral element model Nek5000 developed to exhibit minimal numerical dissipation and dispersion errors, which is advantageous to quantify accurately the impact of SGS stresses. It is found that all tested SGS models lead to improved results with respect to those from DNS*. Also, SGS models allow for simulations with coarse resolutions that blow up in DNS* due to lack of adequate dissipation where needed. The SGS model in which the vertical eddy diffusion is modulated via a function that depends on the Richardson number Ri shows the most faithful reproduction of mixed water masses at all resolutions tested. The sensitivity of the results to the tunable parameter of the SGS model, to changes in the Ri-dependent function and resolution of the turbulent overturning scales is shown.",
author = "Ozgokmen, {Tamay M} and Traian Iliescu and Fischer, {Paul F.} and Ashwanth Srinivasan and Jinqiao Duan",
year = "2007",
doi = "10.1016/j.ocemod.2006.08.006",
language = "English (US)",
volume = "16",
pages = "106--140",
journal = "Ocean Modelling",
issn = "1463-5003",
publisher = "Elsevier BV",
number = "1-2",

}

TY - JOUR

T1 - Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain

AU - Ozgokmen, Tamay M

AU - Iliescu, Traian

AU - Fischer, Paul F.

AU - Srinivasan, Ashwanth

AU - Duan, Jinqiao

PY - 2007

Y1 - 2007

N2 - Mixing in both coastal and deep ocean emerges as one of the important processes that determines the transport of pollutants, sediments and biological species, as well as the details of the global thermohaline circulation. Both the observations, due to their lack in space and time resolution, and most coastal and general circulation models due to inadequate physics, can only provide partial information about oceanic mixing processes. A new class of nonhydrostatic models supplemented with physically based subgrid-scale (SGS) closures, or so-called large eddy simulation (LES), is put forth as another tool of investigation to complement observational and large-scale modeling efforts. However, SGS models have been developed primarily for homogeneous, isotropic flows. Here, four SGS models based on Smagorinsky eddy viscosity and diffusivity are tested for stratified flows in the context of 2D dam-break problem in a rectangular enclosed domain. This idealized testbed leads to a number of simplifications about the initial conditions, boundary conditions and geometry, while exhibiting the dynamically complex characteristics of stratified flows involving the interaction of shear-induced mixing and internal waves. Direct numerical simulations (DNS) at high resolutions are taken as benchmark solutions. Under-resolved simulations without SGS terms (so-called DNS*) are used to quantify the impact of SGS stresses. The performance of LES is assessed by using the time evolution of the volume fraction of intermediate density water masses generated by mixing. The simulations are conducted using a nonhydrostatic high-order spectral element model Nek5000 developed to exhibit minimal numerical dissipation and dispersion errors, which is advantageous to quantify accurately the impact of SGS stresses. It is found that all tested SGS models lead to improved results with respect to those from DNS*. Also, SGS models allow for simulations with coarse resolutions that blow up in DNS* due to lack of adequate dissipation where needed. The SGS model in which the vertical eddy diffusion is modulated via a function that depends on the Richardson number Ri shows the most faithful reproduction of mixed water masses at all resolutions tested. The sensitivity of the results to the tunable parameter of the SGS model, to changes in the Ri-dependent function and resolution of the turbulent overturning scales is shown.

AB - Mixing in both coastal and deep ocean emerges as one of the important processes that determines the transport of pollutants, sediments and biological species, as well as the details of the global thermohaline circulation. Both the observations, due to their lack in space and time resolution, and most coastal and general circulation models due to inadequate physics, can only provide partial information about oceanic mixing processes. A new class of nonhydrostatic models supplemented with physically based subgrid-scale (SGS) closures, or so-called large eddy simulation (LES), is put forth as another tool of investigation to complement observational and large-scale modeling efforts. However, SGS models have been developed primarily for homogeneous, isotropic flows. Here, four SGS models based on Smagorinsky eddy viscosity and diffusivity are tested for stratified flows in the context of 2D dam-break problem in a rectangular enclosed domain. This idealized testbed leads to a number of simplifications about the initial conditions, boundary conditions and geometry, while exhibiting the dynamically complex characteristics of stratified flows involving the interaction of shear-induced mixing and internal waves. Direct numerical simulations (DNS) at high resolutions are taken as benchmark solutions. Under-resolved simulations without SGS terms (so-called DNS*) are used to quantify the impact of SGS stresses. The performance of LES is assessed by using the time evolution of the volume fraction of intermediate density water masses generated by mixing. The simulations are conducted using a nonhydrostatic high-order spectral element model Nek5000 developed to exhibit minimal numerical dissipation and dispersion errors, which is advantageous to quantify accurately the impact of SGS stresses. It is found that all tested SGS models lead to improved results with respect to those from DNS*. Also, SGS models allow for simulations with coarse resolutions that blow up in DNS* due to lack of adequate dissipation where needed. The SGS model in which the vertical eddy diffusion is modulated via a function that depends on the Richardson number Ri shows the most faithful reproduction of mixed water masses at all resolutions tested. The sensitivity of the results to the tunable parameter of the SGS model, to changes in the Ri-dependent function and resolution of the turbulent overturning scales is shown.

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

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

U2 - 10.1016/j.ocemod.2006.08.006

DO - 10.1016/j.ocemod.2006.08.006

M3 - Article

AN - SCOPUS:36749037433

VL - 16

SP - 106

EP - 140

JO - Ocean Modelling

JF - Ocean Modelling

SN - 1463-5003

IS - 1-2

ER -