TY - JOUR
T1 - Large eddy simulation of stratified mixing in two-dimensional dam-break problem in a rectangular enclosed domain
AU - Özgökmen, Tamay M.
AU - Iliescu, Traian
AU - Fischer, Paul F.
AU - Srinivasan, Ashwanth
AU - Duan, Jinqiao
N1 - Funding Information:
We greatly appreciate the support of National Science Foundation via grants OCE 0336799, OCE 0352047, DMS 0209326, DMS 0209309, DMS 0322852, and DMS 0513542, the Air Force Office of Scientific Research grants F49620-03-1-0243 and FA9550-05-1-0449, and the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, US Department of Energy, under Contract W-31-109-ENG-38. We thank Hartmut Peters, Mohamed Iskandarani, Marcello Magaldi and Mehmet Ilicak for valuable comments. We also thank the three anonymous reviewers for their constructive comments, which helped greatly improve the manuscript.
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
SN - 1463-5003
VL - 16
SP - 106
EP - 140
JO - Ocean Modelling
JF - Ocean Modelling
IS - 1-2
ER -