Non-hydrostatic modeling of exchange flows across complex geometries

Mehmet Ilicak, Tamay M Ozgokmen, Emin Özsoy, Paul F. Fischer

Research output: Contribution to journalArticle

26 Citations (Scopus)

Abstract

Flows between ocean and marginal sea basins are often connected by narrow channels and shallow sills. In this study, we address the effect of lateral and vertical geometric constrictions on mixing and hydraulic control. We conduct a set of numerical experiments of the lock-exchange problem in the presence of lateral and vertical contractions. Large eddy simulations (LES) are carried out using the high-order non-hydrostatic spectral element model Nek5000. A spanwise-averaged 2D non-hydrostatic model, denoted SAM, is also employed. Comparison between 2D and 3D models are conducted on the basis of the shape of the density interface and the time evolution of the background potential energy that quantifies the cumulative effects of the stratified mixing in the system. LES indicates that an S-shaped channel induces the highest amount of mixing possibly due to the additional shears caused by the centrifugal force. Vertical and horizontal constrictions enhance mixing by the trapping and breaking of the internal waves in between the obstacles. On the other hand, vertical and horizontal constrictions overlapping with each other restrain the rate of the exchange flow, reduce vertical shears and the mixing. The lowest mixing is encountered in a {n-ary logical or}-shaped channel. It is also found that SAM appears to be accurate in identifying the hydraulic control points due to both horizontal and vertical constrictions. A good agreement is found between SAM and inviscid two-layer theory regarding the steady-state location of the density interface. However SAM overestimates mixing with respect to LES since overturning eddies tend to merge in 2D, while they break down into smaller scales in 3D. SAM is then further tested a realistic application to model the flow in the Bosphorus Strait. Here, the main challenges of using SAM revolved around a non-trivial reduction of 3D geometry to a 2D mapping function, and excessive diffusion with simple closures. The realism of SAM improves significantly using a comprehensive turbulence closure of Very Large Eddy Simulation, VLES. In conclusion, exchange flows in narrow straits pose significant computational challenges due to the details of domain geometry and their impact on mixing. SAM with the VLES turbulence closure appears to be an attractive modeling tool for a first-order assessment of dynamical problems involving mixing and hydraulic effects.

Original languageEnglish (US)
Pages (from-to)159-175
Number of pages17
JournalOcean Modelling
Volume29
Issue number3
DOIs
StatePublished - 2009

Fingerprint

geometry
Geometry
large eddy simulation
Large eddy simulation
modeling
Hydraulics
hydraulics
strait
Turbulence
turbulence
marginal sea
internal wave
Potential energy
potential energy
sill
contraction
trapping
eddy
ocean
basin

Keywords

  • Bosphorus strait flow
  • Exchange flows
  • Hydraulic control
  • VLES

ASJC Scopus subject areas

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

Cite this

Non-hydrostatic modeling of exchange flows across complex geometries. / Ilicak, Mehmet; Ozgokmen, Tamay M; Özsoy, Emin; Fischer, Paul F.

In: Ocean Modelling, Vol. 29, No. 3, 2009, p. 159-175.

Research output: Contribution to journalArticle

Ilicak, Mehmet ; Ozgokmen, Tamay M ; Özsoy, Emin ; Fischer, Paul F. / Non-hydrostatic modeling of exchange flows across complex geometries. In: Ocean Modelling. 2009 ; Vol. 29, No. 3. pp. 159-175.
@article{1f743f102ce54fcea5b3f458d0f419dc,
title = "Non-hydrostatic modeling of exchange flows across complex geometries",
abstract = "Flows between ocean and marginal sea basins are often connected by narrow channels and shallow sills. In this study, we address the effect of lateral and vertical geometric constrictions on mixing and hydraulic control. We conduct a set of numerical experiments of the lock-exchange problem in the presence of lateral and vertical contractions. Large eddy simulations (LES) are carried out using the high-order non-hydrostatic spectral element model Nek5000. A spanwise-averaged 2D non-hydrostatic model, denoted SAM, is also employed. Comparison between 2D and 3D models are conducted on the basis of the shape of the density interface and the time evolution of the background potential energy that quantifies the cumulative effects of the stratified mixing in the system. LES indicates that an S-shaped channel induces the highest amount of mixing possibly due to the additional shears caused by the centrifugal force. Vertical and horizontal constrictions enhance mixing by the trapping and breaking of the internal waves in between the obstacles. On the other hand, vertical and horizontal constrictions overlapping with each other restrain the rate of the exchange flow, reduce vertical shears and the mixing. The lowest mixing is encountered in a {n-ary logical or}-shaped channel. It is also found that SAM appears to be accurate in identifying the hydraulic control points due to both horizontal and vertical constrictions. A good agreement is found between SAM and inviscid two-layer theory regarding the steady-state location of the density interface. However SAM overestimates mixing with respect to LES since overturning eddies tend to merge in 2D, while they break down into smaller scales in 3D. SAM is then further tested a realistic application to model the flow in the Bosphorus Strait. Here, the main challenges of using SAM revolved around a non-trivial reduction of 3D geometry to a 2D mapping function, and excessive diffusion with simple closures. The realism of SAM improves significantly using a comprehensive turbulence closure of Very Large Eddy Simulation, VLES. In conclusion, exchange flows in narrow straits pose significant computational challenges due to the details of domain geometry and their impact on mixing. SAM with the VLES turbulence closure appears to be an attractive modeling tool for a first-order assessment of dynamical problems involving mixing and hydraulic effects.",
keywords = "Bosphorus strait flow, Exchange flows, Hydraulic control, VLES",
author = "Mehmet Ilicak and Ozgokmen, {Tamay M} and Emin {\"O}zsoy and Fischer, {Paul F.}",
year = "2009",
doi = "10.1016/j.ocemod.2009.04.002",
language = "English (US)",
volume = "29",
pages = "159--175",
journal = "Ocean Modelling",
issn = "1463-5003",
publisher = "Elsevier BV",
number = "3",

}

TY - JOUR

T1 - Non-hydrostatic modeling of exchange flows across complex geometries

AU - Ilicak, Mehmet

AU - Ozgokmen, Tamay M

AU - Özsoy, Emin

AU - Fischer, Paul F.

PY - 2009

Y1 - 2009

N2 - Flows between ocean and marginal sea basins are often connected by narrow channels and shallow sills. In this study, we address the effect of lateral and vertical geometric constrictions on mixing and hydraulic control. We conduct a set of numerical experiments of the lock-exchange problem in the presence of lateral and vertical contractions. Large eddy simulations (LES) are carried out using the high-order non-hydrostatic spectral element model Nek5000. A spanwise-averaged 2D non-hydrostatic model, denoted SAM, is also employed. Comparison between 2D and 3D models are conducted on the basis of the shape of the density interface and the time evolution of the background potential energy that quantifies the cumulative effects of the stratified mixing in the system. LES indicates that an S-shaped channel induces the highest amount of mixing possibly due to the additional shears caused by the centrifugal force. Vertical and horizontal constrictions enhance mixing by the trapping and breaking of the internal waves in between the obstacles. On the other hand, vertical and horizontal constrictions overlapping with each other restrain the rate of the exchange flow, reduce vertical shears and the mixing. The lowest mixing is encountered in a {n-ary logical or}-shaped channel. It is also found that SAM appears to be accurate in identifying the hydraulic control points due to both horizontal and vertical constrictions. A good agreement is found between SAM and inviscid two-layer theory regarding the steady-state location of the density interface. However SAM overestimates mixing with respect to LES since overturning eddies tend to merge in 2D, while they break down into smaller scales in 3D. SAM is then further tested a realistic application to model the flow in the Bosphorus Strait. Here, the main challenges of using SAM revolved around a non-trivial reduction of 3D geometry to a 2D mapping function, and excessive diffusion with simple closures. The realism of SAM improves significantly using a comprehensive turbulence closure of Very Large Eddy Simulation, VLES. In conclusion, exchange flows in narrow straits pose significant computational challenges due to the details of domain geometry and their impact on mixing. SAM with the VLES turbulence closure appears to be an attractive modeling tool for a first-order assessment of dynamical problems involving mixing and hydraulic effects.

AB - Flows between ocean and marginal sea basins are often connected by narrow channels and shallow sills. In this study, we address the effect of lateral and vertical geometric constrictions on mixing and hydraulic control. We conduct a set of numerical experiments of the lock-exchange problem in the presence of lateral and vertical contractions. Large eddy simulations (LES) are carried out using the high-order non-hydrostatic spectral element model Nek5000. A spanwise-averaged 2D non-hydrostatic model, denoted SAM, is also employed. Comparison between 2D and 3D models are conducted on the basis of the shape of the density interface and the time evolution of the background potential energy that quantifies the cumulative effects of the stratified mixing in the system. LES indicates that an S-shaped channel induces the highest amount of mixing possibly due to the additional shears caused by the centrifugal force. Vertical and horizontal constrictions enhance mixing by the trapping and breaking of the internal waves in between the obstacles. On the other hand, vertical and horizontal constrictions overlapping with each other restrain the rate of the exchange flow, reduce vertical shears and the mixing. The lowest mixing is encountered in a {n-ary logical or}-shaped channel. It is also found that SAM appears to be accurate in identifying the hydraulic control points due to both horizontal and vertical constrictions. A good agreement is found between SAM and inviscid two-layer theory regarding the steady-state location of the density interface. However SAM overestimates mixing with respect to LES since overturning eddies tend to merge in 2D, while they break down into smaller scales in 3D. SAM is then further tested a realistic application to model the flow in the Bosphorus Strait. Here, the main challenges of using SAM revolved around a non-trivial reduction of 3D geometry to a 2D mapping function, and excessive diffusion with simple closures. The realism of SAM improves significantly using a comprehensive turbulence closure of Very Large Eddy Simulation, VLES. In conclusion, exchange flows in narrow straits pose significant computational challenges due to the details of domain geometry and their impact on mixing. SAM with the VLES turbulence closure appears to be an attractive modeling tool for a first-order assessment of dynamical problems involving mixing and hydraulic effects.

KW - Bosphorus strait flow

KW - Exchange flows

KW - Hydraulic control

KW - VLES

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

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

U2 - 10.1016/j.ocemod.2009.04.002

DO - 10.1016/j.ocemod.2009.04.002

M3 - Article

AN - SCOPUS:67349262547

VL - 29

SP - 159

EP - 175

JO - Ocean Modelling

JF - Ocean Modelling

SN - 1463-5003

IS - 3

ER -