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.
N1 - Copyright:
Copyright 2008 Elsevier B.V., All rights reserved.
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
AN - SCOPUS:3943055376
VL - 132
SP - 1777
EP - 1791
JO - Monthly Weather Review
JF - Monthly Weather Review
SN - 0027-0644
IS - 7
ER -