Developing a nonhydrostatic isopycnalcoordinate

Similar documents
Internal Tide Generation Over a Continental Shelf: Analytical and Numerical Calculations

Internal Wave Modeling in Oceanic General Circulation Models

Internal Wave Driven Mixing and Transport in the Coastal Ocean

Erratic internal waves at SIO Pier. data and wavelet analysis courtesy of E. Terrill, SIO

Recent developments of the NERSC version of MICOM

INTERNAL GRAVITY WAVES

Sediment resuspension from shoaling internal solitary waves. Daniel Bourgault, Institut des sciences de la mer de Rimouski, Rimouski, QC, Canada

Paper C Numerical studies of nonhydrostatic pressure effects on the wind forced boundary layer dynamics

Numerical simulations of the interaction of internal waves with a shelf break

Evolution of Internal Waves in Hydrostatic Models

Limitations of hydrodynamic models in simulating Great Lakes processes of varying scale Introduction: Basin-scale waves and wave-driven currents:

Local generation of internal solitary waves in an oceanic pycnocline

On the Calculation of Available Potential Energy in Internal Wave Fields

Goals of this Chapter

Donald Slinn, Murray D. Levine

Energy flux of nonlinear internal waves in northern South China Sea

The effect of a background shear current on large amplitude internal solitary waves

The Application of POM to the Operational Tidal Forecast for the Sea around Taiwan

Laboratory Modeling of Internal Wave Generation in Straits

Generation and Evolution of Internal Waves in Luzon Strait

Stability and consistency of nonhydrostatic free-surface models using the semi-implicit -method

2017 年環境流體力學短期講座 Short Course on Environmental Flows

Hydrodynamics in Shallow Estuaries with Complex Bathymetry and Large Tidal Ranges

OCEAN MODELING. Joseph K. Ansong. (University of Ghana/Michigan)

The Hybrid-Coordinate Ocean Model (HYCOM) Alex Megann, Southampton Oceanography Centre

Collaborative Proposal to Extend ONR YIP research with BRC Efforts

v t + fu = 1 p y w t = 1 p z g u x + v y + w

IWISE: Observations of Internal Wave and Tide Generation in the Luzon Strait Using Shipboard Surveys and Oceanographic Moorings

Internal solitary wave transformation over the bottom step: loss of energy

Modeling internal tides and mixing over ocean ridges

A Two-Level Time-Stepping Method for Layered Ocean Circulation Models: Further Development and Testing

A process study of tidal mixing over rough topography

Instability of a coastal jet in a two-layer model ; application to the Ushant front

Collaborative Proposal to Extend ONR YIP research with BRC Efforts

Climate sensitivity of coupled models with differing ocean components

Concurrent simulation of the eddying general circulation and tides in a global ocean model

An unambiguous definition of the Froude number for lee waves in the deep ocean

A comparison of two formulations of barotropic-baroclinic splitting for layered models of ocean circulation

The Luzon Strait experiment

INTERNAL SOLITARY WAVES IN THE ATMOSPHERE AND OCEAN

Non-hydrostatic pressure in σ-coordinate ocean models

Dynamics of Strongly Nonlinear Internal Solitary Waves in Shallow Water

Contents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition

SMS 303: Integrative Marine

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

Internal Wave Generation and Scattering from Rough Topography

Modeling of Coastal Ocean Flow Fields

Modeling the Impact of Extreme Events on Margin Sedimentation

A Statistical Investigation of Internal Wave Propagation in the Northern South China Sea

A simple model of topographic Rossby waves in the ocean

Split explicit methods

2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven Ocean Fronts

Data Assimilation and Diagnostics of Inner Shelf Dynamics

Modelling the Channel and the Bay of Biscay using the MARS model. Contributions: Dyneco Physed / Pelagos / Benthos EMH

Internal Waves in the Vicinity of the Kuroshio Path

Progress in embedding global tides into HYCOM simulations. Brian K. Arbic Institute for Geophysics The University of Texas at Austin

Oliver Bühler Waves and Vortices

OCN660 - Ocean Waves. Course Purpose & Outline. Doug Luther. OCN660 - Syllabus. Instructor: x65875

Large-scale, realistic laboratory modeling of M 2 internal tide generation at the Luzon Strait

MARINE RESEARCH. Journal of. Interfacial solitary wave run-up in the St. Lawrence Estuary. Volume 63, Number 6

Internal Wave Generation in Straits

Gravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017

Tide-Topography Interactions: Asymmetries in Internal Wave Generation due to Surface Trapped Currents

Assessment of the numerical efficiency of ocean circulation model : Hycom contribution to the COMODO project.

Open boundary conditions in numerical simulations of unsteady incompressible flow

Internal tides in NEMO in the Indonesian seas. A. Koch-Larrouy, D. Nugroho, F. Lyard, D. Allain, B. Tranchant, P. Gaspar, J.

Generation, Propagation, and Dissipation of Internal Waves in Continental Shelf Regimes

OCEAN HYDRODYNAMIC MODEL

The Lagrangian-Averaged Navier-Stokes alpha (LANS-α) Turbulence Model in Primitive Equation Ocean Modeling

Inversion of tsunami data. A. Sladen CNRS, Géoazur 1/35

3.3 Classification Diagrams Estuarine Zone Coastal Lagoons References Physical Properties and Experiments in

SAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany

Thomas Pierro, Donald Slinn, Kraig Winters

NUMERICAL SOLITARY WAVE INTERACTION: THE ORDER OF THE INELASTIC EFFECT

Exact Solitary Wave Solutions in Shallow Water

Modeling Internal Tides and Mixing Over Ocean Ridges

( ) = 1005 J kg 1 K 1 ;

Energy Budget of Nonlinear Internal Waves near Dongsha

Hepeng Zhang Ben King Bruce Rodenborn

Lecture 23. Vertical coordinates General vertical coordinate

Comparison between vertical shear mixing and surface wave-induced mixing in the global ocean

Lecture-1: FVCOM-An unstructured grid Finite-Volume Community Ocean Model

Internal boundary layers in the ocean circulation

Strongly nonlinear long gravity waves in uniform shear flows

Tilting Shear Layers in Coastal Flows

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity

The Effect of Statistical Abyssal Hill Roughness on the Generation of Internal Waves

Nonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation

Semi-enclosed seas. Estuaries are only a particular type of semi-enclosed seas which are influenced by tides and rivers

Ocean self-attraction and loading (SAL) and internal tides dissipation implementation within an unstructured global tide-surge model

Nonlinear Disintegration of the Internal Tide

OCEAN MODELING II. Parameterizations

MAR513 Lecture 10: Pressure Errors in Terrain-Following Coordinates

Numerical Study of Baroclinic Tides in Luzon Strait

Buoyancy Fluxes in a Stratified Fluid

Continuously stratified nonlinear low-mode internal tides

Lecture 1. Equations of motion - Newton s second law in three dimensions. Pressure gradient + force force

Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway

2 Observing the Ocean Ships Navigation The Electronics Era 16

Transcription:

Developing a nonhydrostatic isopycnalcoordinate ocean model Oliver Fringer Associate Professor The Bob and Norma Street Environmental Fluid Mechanics Laboratory Dept. of Civil and Environmental Engineering Stanford University 11 September 2015 Funding: ONR Grants N00014-10-1-0521 and N00014-08-1-0904

Outline Overview of internal gravity waves Nonhydrostatic (Navier-Stokes) modeling Grid-resolution requirements Nonhydrostatic modeling is expensive! A nonhydrostatic isopycnal-coordinate model The cost can be reduced! Conclusions

Internal gravity waves Internal waves z Interfacial waves z r r Venayagamoorthy & Fringer (2007) Arthur & Fringer (2014) Typical speeds in the ocean: 1-3 m/s Frequencies: Tidal (internal tides) - minutes (internal waves) Wavelengths: 100s of km to 10s of m

Surface signatures induced by internal gravity waves Straight of Gibraltar South China Sea Surface Convergence (rough) Divergence (smooth) Convergence (rough) SAR Image Courtesy internalwaveatlas.com Propagation direction

Applications of internal gravity waves Breaking of internal tides and waves may provide the necessary mixing to maintain the ocean stratification (Munk and Wunsch 1998). Cold Antarctic water warm equatorial surface waters Cold Arctic water Mixing induced by breaking internal waves prevents the ocean from turning into a "stagnant pool of cold, salty water"... Internal waves are hypothesized to deliver nutrients that sustain thriving coral reef ecosystems (e.g. Florida Shelf, Leichter et al., 2003; Dongsha Atol, Wang et al. 2007) Internal waves influence sediment transport in lakes and oceans and propagation of acoustic signals. Strong internal wave-induced currents can cause oil platform instability and pipeline rupture.

Isopycnal vs z- or sigma-coordinates Advantages of isopycnal coordinates: Reduces the number of vertical grid points from O(100) in traditional coordinates to O(1-10) No spurious vertical (diapycnal) diffusion/mixing Challenges of isopycnal coordinates: Cannot represent unstable stratification Layer outcropping (drying of layers) requires special numerical schemes Hydrostatic

Hydrostatic vs. nonhydrostatic flows Most ocean flows are hydrostatic Long horizontal length scales relative to vertical length scales, i.e. long waves (i.e. Lh >> Lv) Only in small regions is the flow nonhydrostatic Short horizontal length scales relative to vertical scales (i.e. Lh ~ Lv) Can cost 10X more to compute! free surface Long wave (Lh>>Lv) (hydrostatic) Steep bathymetry Lh~Lv (nonhydrostatic) bottom

Nonhydrostatic effects: Overturning Hydrostatic Nonhydrostatic Overturning motions and eddies are not the only nonhydrostatic process

Nonhydrostatic effects: Frequency dispersion of gravity waves Dispersion relation for irrotational surface gravity waves: c 2 g k tanh kd g tanh e, k Deep-water limit: e>>1 (nonhydrostatic) 2 c Shallow-water limit e<<1 (hydrostatic) g k e D L 2 c gd l=2l k=2 /l /L

When is a flow nonhydrostatic? Nonhydrostatic Model Hydrostatic Model Aspect Ratio: e = D/L = 2

Nonhydrostatic Model Hydrostatic Model Aspect Ratio: e = D/L = 1/8 = 0.125 Nonhydrostatic result = Hydrostatic result + e 2

Example 3D nonhydrostatic z-level simulation: Internal gravity waves in the South China Sea 15 o isotherm m From: Zhang and Fringer (2011) China Taiwan Grid Luzon resolution: (Philippines) Horizontal: Dx=1 km Vertical: 100 z-levels (Dz~10 m) Number of 3D cells: 12 million

Generation of weakly nonlinear wavetrains Isotherms: 16, 20, 24, 28 degrees C Long internal tides O(100 km) Short, solitary-like waves O(5 km) How can we determine, apriori, how much horizontal grid resolution is needed to simulate this process? x 100 km

Internal solitary waves Width L Upper layer h 1 density r 1 a density r 2 Lower layer h 2 >> h 1 Speed c Nonlinear effect (steepening): d = a/h 1 Nonhydrostatic effect (frequency dispersion): e = h 1 /L Solitary wave: Balance between nonlinear steepening and nonhydrostatic dispersion. d ~ e 2

The KdV equation When computing solitary waves, the behavior of a 3D, fully nonhydrostatic ocean model can be approximated very well with the KdV (Korteweg and de-vries, 1895) equation: Ocean Model: u t Unsteadiness u u Nonlinear momentum advection p Hydrostatic pressure gradient. q Nonhydrostatic pressure gradient. KdV: t 3 d 2 x x 2 e 6 3 3 x z = +O(e 2 d, e 4 ) The KdV equation gives the well-known solution ( x, t) asech 2 x ct L 0 L 0 4 e 2 3 da

Numerical discretization of KdV Many ocean models discretize the equations with second-order accuracy in time and space. (e.g. SUNTANS, Fringer et al. 2006; POM, Blumberg and Mellor, 1987; MICOM, Bleck et al., 1992; MOM, Pacanowski and Griffes, 1999). A second-order accurate discretization of the KdV equation using leapfrog (i.e. POM) is given by Use the Taylor series expansion to determine the modified equivalent form of the terms, e.g.

Modified equivalent KdV equation The discrete form of the KdV equation produces a solution to the modified equivalent PDE (Hirt 1968) which introduces new terms due to discretization errors: KdV Modified equivalent KdV: 2 Dx Numerical dispersion K h 1 Physical dispersion 2 Kl l=dx/h 1 grid "lepticity" (Scotti and Mitran, 2008) K=O(1) constant. The numerical discretization of the first-order derivative produces numerical dispersion. Note that the errors in the nonlinear term are smaller by a factor d. For numerical dispersion to be smaller than physical dispersion, l < 1. Vitousek and Fringer (2011)

Hydrostatic vs. nonhydrostatic for l=0.25 (Dx=h 1 /4) Numerical dispersion is 16 times smaller than physical dispersion. Vitousek and Fringer (2011)

Hydrostatic vs. nonhydrostatic for l=8 (Dx=8h 1 ) "Numerical solitary waves!" Numerical dispersion is 64 times larger than physical dispersion. Vitousek and Fringer (2011)

Nonhydrostatic isopycnal model? Zhang et al. simulation: 12 million cells, Dx=1 km = 5 h 1 To begin to resolve nonhydrostatic effects, Dx=200 m = h 1 300 million cells! With Dx=100 m, 1.2 billion! The z-level SUNTANS model requires O(100) z-levels to minimize numerical diffusion of the pycnocline. Solution: Isopycnal model with O(2) layers = 50X reduction in computation time. Nonhydrostatic isopycnal coordinate model. 2-layer hydrostatic result with isopycnal model of Simmons, U. Alaska Fairbanks.

Essential features of the nonhydrostatic isopycnal-coordinate model Staggered C-grid layout Split Montgomery potential into Barotropic (implicit) & Baroclinic (explicit) parts IMplicit-EXplicit (IMEX) multistep method. Durran (2012) MPDATA for upwinding of layer heights Implicit theta method for vertical diffusion Explicit horizontal diffusion Predictor/corrector method for nonhydrostatic pressure Second-order accurate in time and space Vitousek and Fringer (2014)

Nonhydrostatic test cases No stratification Note: density need not change in each layer Two-layer Smooth pycnocline

Hydrostatic internal seiche 2 layers 100 layers

Nonhydrostatic internal seiche 2 layers 100 layers

Dispersion relation speed = function(wavelength) Vitousek and Fringer (2014)

Internal solitary wave formation z-level model leads to numerical diffusion, or thickening of the pycnocline. Isopycnalcoordinate model: eliminates spurious numerical diffusion captures solitary wave behavior at 1/50 cost Vitousek and Fringer (2014)

Internal Solitary Waves in the (idealized) South China Sea (SCS) 10-layer isopycnal model following Buijsman et al. (2010) Vitousek and Fringer (2014)

Internal Wave Runup LES model 2-layer model 2-layer isopycnal model vs. an LES model (Bobby Arthur, 2014) Vitousek and Fringer (2014) Test case similar to: Michallet & Ivey 1999 Bourgault & Kelley 2004

Conclusions Simulation of nonhydrostatic effects in the SCS requires Dx<h 1 O(billion) grid cells in 3D with z-level model. We have developed a nonhydrostatic isopycnalcoordinate model using stable higher-order timestepping. More isopycnal layers are needed: To resolve stratification To resolve nonhydrostatic effects Most oceanic/lake processes are weakly nonhydrostatic and so <10 layers suffice for many applications. The result is a reduced computational cost by O(10). Ongoing work: Development of unstructured-grid model.