A spectral finite-volume method for the shallow water equations

Byoung Ju Choi, Mohamed Iskandarani, Julia Levin, Dale B. Haidvogel

Research output: Contribution to journalArticle

22 Citations (Scopus)

Abstract

A spectral finite-volume (SFV) method is proposed for the numerical solution of the shallow water equations. This is the first phase in the development of a layered (isopycnal) ocean model. Its target applications include, in particular, the simulation of the wind-driven oceanic circulation in geometrically complex basins where layer outcropping and/or isopycnal-bathymetry intersection must be handled explicitly. The present formulation is geometrically flexible and can extend accuracy to arbitrary high order with no change to the basic algorithm. A flux-corrected transport (FCT) algorithm ensures the stability of the computations in regions of vanishing layer thickness and in areas where the flow features are underresolved. The spatial discretization is based on a two-level grid: a globally unstructured elemental grid and a locally structured grid consisting of N X N quadrilateral cells within each element. The numerical solution is continuous within each element but discontinuous across elements; the discontinuity is resolved by upwinding along characteristics. The accuracy and convergence rate of the SFV method are verified on two linearized problems amenable to analytical solution; the SFV solution exhibits a convergence order of N + 1 for smooth solutions. The FCT portion of the model is tested by simulating the formation of an oblique hydraulic jump in a supercritical channel flow. The model is then applied to simulate, in reduced-gravity mode, the double-gyre and wind-driven upper-ocean circulations in a square basin. Finally, the previous experiment is repeated in the North Atlantic basin to illustrate the application of the model in a realistic geometry.

Original languageEnglish (US)
Pages (from-to)1777-1791
Number of pages15
JournalMonthly Weather Review
Volume132
Issue number7
DOIs
StatePublished - Jul 2004

Fingerprint

finite volume method
shallow-water equation
basin
wind-driven circulation
supercritical flow
oceanic circulation
channel flow
upper ocean
gyre
bathymetry
discontinuity
gravity
hydraulics
geometry
ocean
simulation
experiment

ASJC Scopus subject areas

  • Atmospheric Science

Cite this

A spectral finite-volume method for the shallow water equations. / Choi, Byoung Ju; Iskandarani, Mohamed; Levin, Julia; Haidvogel, Dale B.

In: Monthly Weather Review, Vol. 132, No. 7, 07.2004, p. 1777-1791.

Research output: Contribution to journalArticle

Choi, Byoung Ju ; Iskandarani, Mohamed ; Levin, Julia ; Haidvogel, Dale B. / A spectral finite-volume method for the shallow water equations. In: Monthly Weather Review. 2004 ; Vol. 132, No. 7. pp. 1777-1791.
@article{e2e76cd7a46a4ae69da229a236f9e31c,
title = "A spectral finite-volume method for the shallow water equations",
abstract = "A spectral finite-volume (SFV) method is proposed for the numerical solution of the shallow water equations. This is the first phase in the development of a layered (isopycnal) ocean model. Its target applications include, in particular, the simulation of the wind-driven oceanic circulation in geometrically complex basins where layer outcropping and/or isopycnal-bathymetry intersection must be handled explicitly. The present formulation is geometrically flexible and can extend accuracy to arbitrary high order with no change to the basic algorithm. A flux-corrected transport (FCT) algorithm ensures the stability of the computations in regions of vanishing layer thickness and in areas where the flow features are underresolved. The spatial discretization is based on a two-level grid: a globally unstructured elemental grid and a locally structured grid consisting of N X N quadrilateral cells within each element. The numerical solution is continuous within each element but discontinuous across elements; the discontinuity is resolved by upwinding along characteristics. The accuracy and convergence rate of the SFV method are verified on two linearized problems amenable to analytical solution; the SFV solution exhibits a convergence order of N + 1 for smooth solutions. The FCT portion of the model is tested by simulating the formation of an oblique hydraulic jump in a supercritical channel flow. The model is then applied to simulate, in reduced-gravity mode, the double-gyre and wind-driven upper-ocean circulations in a square basin. Finally, the previous experiment is repeated in the North Atlantic basin to illustrate the application of the model in a realistic geometry.",
author = "Choi, {Byoung Ju} and Mohamed Iskandarani and Julia Levin and Haidvogel, {Dale B.}",
year = "2004",
month = "7",
doi = "10.1175/1520-0493(2004)132<1777:ASFMFT>2.0.CO;2",
language = "English (US)",
volume = "132",
pages = "1777--1791",
journal = "Monthly Weather Review",
issn = "0027-0644",
publisher = "American Meteorological Society",
number = "7",

}

TY - JOUR

T1 - A spectral finite-volume method for the shallow water equations

AU - Choi, Byoung Ju

AU - Iskandarani, Mohamed

AU - Levin, Julia

AU - Haidvogel, Dale B.

PY - 2004/7

Y1 - 2004/7

N2 - A spectral finite-volume (SFV) method is proposed for the numerical solution of the shallow water equations. This is the first phase in the development of a layered (isopycnal) ocean model. Its target applications include, in particular, the simulation of the wind-driven oceanic circulation in geometrically complex basins where layer outcropping and/or isopycnal-bathymetry intersection must be handled explicitly. The present formulation is geometrically flexible and can extend accuracy to arbitrary high order with no change to the basic algorithm. A flux-corrected transport (FCT) algorithm ensures the stability of the computations in regions of vanishing layer thickness and in areas where the flow features are underresolved. The spatial discretization is based on a two-level grid: a globally unstructured elemental grid and a locally structured grid consisting of N X N quadrilateral cells within each element. The numerical solution is continuous within each element but discontinuous across elements; the discontinuity is resolved by upwinding along characteristics. The accuracy and convergence rate of the SFV method are verified on two linearized problems amenable to analytical solution; the SFV solution exhibits a convergence order of N + 1 for smooth solutions. The FCT portion of the model is tested by simulating the formation of an oblique hydraulic jump in a supercritical channel flow. The model is then applied to simulate, in reduced-gravity mode, the double-gyre and wind-driven upper-ocean circulations in a square basin. Finally, the previous experiment is repeated in the North Atlantic basin to illustrate the application of the model in a realistic geometry.

AB - A spectral finite-volume (SFV) method is proposed for the numerical solution of the shallow water equations. This is the first phase in the development of a layered (isopycnal) ocean model. Its target applications include, in particular, the simulation of the wind-driven oceanic circulation in geometrically complex basins where layer outcropping and/or isopycnal-bathymetry intersection must be handled explicitly. The present formulation is geometrically flexible and can extend accuracy to arbitrary high order with no change to the basic algorithm. A flux-corrected transport (FCT) algorithm ensures the stability of the computations in regions of vanishing layer thickness and in areas where the flow features are underresolved. The spatial discretization is based on a two-level grid: a globally unstructured elemental grid and a locally structured grid consisting of N X N quadrilateral cells within each element. The numerical solution is continuous within each element but discontinuous across elements; the discontinuity is resolved by upwinding along characteristics. The accuracy and convergence rate of the SFV method are verified on two linearized problems amenable to analytical solution; the SFV solution exhibits a convergence order of N + 1 for smooth solutions. The FCT portion of the model is tested by simulating the formation of an oblique hydraulic jump in a supercritical channel flow. The model is then applied to simulate, in reduced-gravity mode, the double-gyre and wind-driven upper-ocean circulations in a square basin. Finally, the previous experiment is repeated in the North Atlantic basin to illustrate the application of the model in a realistic geometry.

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

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

U2 - 10.1175/1520-0493(2004)132<1777:ASFMFT>2.0.CO;2

DO - 10.1175/1520-0493(2004)132<1777:ASFMFT>2.0.CO;2

M3 - Article

VL - 132

SP - 1777

EP - 1791

JO - Monthly Weather Review

JF - Monthly Weather Review

SN - 0027-0644

IS - 7

ER -