FVCOM validation experiments : comparisons with ROMS for three idealized barotropic test problems

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C07042, doi:10.1029/2007JC004557.The unstructured-grid Finite-Volume Coastal Ocean Model (FVCOM) is evaluated using three idealized benchmark test problems: the Rossby equatorial soliton, the hydraulic jump, and the three-dimensional barotropic wind-driven basin. These test cases examine the properties of numerical dispersion and damping, the performance of the nonlinear advection scheme for supercritical flow conditions, and the accuracy of the implicit vertical viscosity scheme in barotropic settings, respectively. It is demonstrated that FVCOM provides overall a second-order spatial accuracy for the vertically averaged equations (i.e., external mode), and with increasing grid resolution the model-computed solutions show a fast convergence toward the analytic solutions regardless of the particular triangulation method. Examples are provided to illustrate the ability of FVCOM to facilitate local grid refinement and speed up computation. Comparisons are also made between FVCOM and the structured-grid Regional Ocean Modeling System (ROMS) for these test cases. For the linear problem in a simple rectangular domain, i.e., the wind-driven basin case, the performance of the two models is quite similar. For the nonlinear case, such as the Rossby equatorial soliton, the second-order advection scheme used in FVCOM is almost as accurate as the fourth-order advection scheme implemented in ROMS if the horizontal resolution is relatively high. FVCOM has taken advantage of the new development in computational fluid dynamics in resolving flow problems containing discontinuities. One salient feature illustrated by the three-dimensional barotropic wind-driven basin case is that FVCOM and ROMS simulations show different responses to the refinement of grid size in the horizontal and in the vertical.For this work, H. Huang and G. Cowles were supported by the Massachusetts Marine Fisheries Institute (MFI) through NOAA grants DOC/NOAA/NA04NMF4720332 and DOC/ NOAA/NA05NMF472113; C. Chen was supported by NSF grants (OCE0234545, OCE0606928, OCE0712903, OCE0732084, and OCE0726851), NOAA grants (NA160P2323, NA06RG0029, and NA960P0113), MIT Sea grant (2006-RC-103), and Georgia Sea grant (NA26RG0373 and NA66RG0282); C. Winant was supported through NSF grant OCE-0726673; R. Beardsley was supported through NSF OCE—0227679 and the WHOI Smith Chair; K. Hedstrom was supported through NASA grant NAG13– 03021 and the Arctic Region Supercomputing Center; and D. Haidvogel was supported through grants ONR N00014- 03-1-0683 and NSF OCE 043557