A three-dimensional spectral element model for the solution of the hydrostatic primitive equations

We present a spectral element model to solve the hydrostatic primitive equations governing large-scale geophysical flows. The highlights of this new model include unstructured grids, dual h–p paths to convergence, and good scalability characteristics on present day parallel computers including Beowulf-class systems. The behavior of the model is assessed on three process-oriented test problems involving wave propagation, gravitational adjustment, and nonlinear flow rectification, respectively. The first of these test problems is a study of the convergence properties of the model when simulating the linear propagation of baroclinic Kelvin waves. The second is an intercomparison of spectral element and finite-difference model solutions to the adjustment of a density front in a straight channel. Finally, the third problem considers the comparison of model results to measurements obtained from a laboratory simulation of flow around a submarine canyon. The aforementioned tests demonstrate the good performance of the model in the idealized/processoriented limits. 2003 Elsevier Science B.V. All rights reserved.

[1]  E. Curchitser,et al.  On the transient adjustment of a mid-latitude abyssal ocean basin with realistic geometry: the constant depth limit , 1999 .

[2]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models, Part 1: Scheme Design and Diagnostic Analysis, Part II: Energy, Momentum, and Bottom Torque Consistency , 1998 .

[3]  M. Taylor.  The Spectral Element Method for the Shallow Water Equations on the Sphere , 1997 .

[4]  Beth A. Wingate,et al.  A generalized diagonal mass matrix spectral element method for non-quadrilateral elements , 2000 .

[5]  Dale B. Haidvogel,et al.  A nonconforming spectral element ocean model , 2000 .

[6]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[7]  Centro internazionale matematico estivo. Session,et al.  Advanced Numerical Approximation of Nonlinear Hyperbolic Equations , 1998 .

[8]  D. Haidvogel,et al.  Laboratory–Numerical Model Comparisons of Flow over a Coastal Canyon , 2001 .

[9]  George Em Karniadakis,et al.  Regular Article: A Discontinuous Galerkin ALE Method for Compressible Viscous Flows in Moving Domains , 1999 .

[10]  Dale B. Haidvogel,et al.  A Spectral Element Solution of the Shallow-Water Equations on Multiprocessor Computers , 1998 .

[11]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[12]  Mark A. Taylor,et al.  Global Modelling of the Ocean and Atmosphere Using the Spectral Element Method , 1997 .

[13]  Hong Ma,et al.  A spectral element basin model for the shallow water equations , 1993 .

[14]  J. W. Eastwood,et al.  Springer series in computational physics Editors: H. Cabannes, M. Holt, H.B. Keller, J. Killeen and S.A. Orszag , 1984 .

[15]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[16]  W. Large,et al.  Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization , 1994 .

[17]  M. Iskandarani,et al.  A Spectral Filtering Procedure for Eddy-Resolving Simulations with a Spectral Element Ocean Model , 1997 .

[18]  Dale B. Haidvogel,et al.  Numerical Ocean Circulation Modeling , 1999 .

[19]  Bernardo Cockburn.  An introduction to the Discontinuous Galerkin method for convection-dominated problems , 1998 .

[20]  Mark A. Taylor,et al.  An Algorithm for Computing Fekete Points in the Triangle , 2000, SIAM J. Numer. Anal..

[21]  C. Wunsch,et al.  Dynamics of the long-period tides , 1997 .

[22]  Mark A. Taylor,et al.  Performance of a spectral element atmospheric model (SEAM) on the HP Exemplar SPP2000 , 1998 .

[23]  G. Karami.  Lecture Notes in Engineering , 1989 .

[24]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part I: Scheme Design and Diagnostic Analysis , 1998 .

[25]  Robert L. Higdon,et al.  Stability Analysis of Operator Splitting for Large-Scale Ocean Modeling , 1996 .

[26]  Transient Adjustment of Circulation in a Midlatitude Abyssal Ocean Basin with Realistic Geometry and Bathymetry , 2001 .

[27]  J. Boyd,et al.  A staggered spectral element model with application to the oceanic shallow , 1995 .

[28]  D. Haidvogel,et al.  A semi-implicit ocean circulation model using a generalized topography-following coordinate system , 1994 .