A New WENO‐Based Momentum Advection Scheme for Simulations of Ocean Mesoscale Turbulence
Name
J Adv Model Earth Syst - 2024 - Silvestri - A New WENO‐Based Momentum Advection Scheme for Simulations of Ocean Mesoscale.pdf
Description
Published version
Size
9.86 MB
Format
Adobe PDF
Checksum (MD5)
cf953d5e23333259d42a849f30780bc2
Author(s)
Silvestri, Simone
Wagner, Gregory L
Campin, Jean‐Michel
Constantinou, Navid C
Hill, Christopher N
Souza, Andre
Ferrari, Raffaele
Date Issued
July 15, 2024
Journal
Journal of Advances in Modeling Earth Systems
Publisher
American Geophysical Union
Citation
Silvestri, S., Wagner, G. L., Campin, J.-M., Constantinou, N. C., Hill, C. N., Souza, A., & Ferrari, R. (2024). A new WENO-based momentum advection scheme for simulations of ocean mesoscale turbulence. Journal of Advances in Modeling Earth Systems, 16, e2023MS004130.
Version
Final published version
Abstract
Current eddy‐permitting and eddy‐resolving ocean models require dissipation to prevent a spurious accumulation of enstrophy at the grid scale. We introduce a new numerical scheme for momentum advection in large‐scale ocean models that involves upwinding through a weighted essentially non‐oscillatory (WENO) reconstruction. The new scheme provides implicit dissipation and thereby avoids the need for an additional explicit dissipation that may require calibration of unknown parameters. This approach uses the rotational, “vector invariant” formulation of the momentum advection operator that is widely employed by global general circulation models. A novel formulation of the WENO “smoothness indicators” is key for avoiding excessive numerical dissipation of kinetic energy and enstrophy at grid‐resolved scales. We test the new advection scheme against a standard approach that combines explicit dissipation with a dispersive discretization of the rotational advection operator in two scenarios: (a) two‐dimensional turbulence and (b) three‐dimensional baroclinic equilibration. In both cases, the solutions are stable, free from dispersive artifacts, and achieve increased “effective” resolution compared to other approaches commonly used in ocean models.
MIT Department
Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences
Terms of Use
Creative Commons Attribution
Persistent DSpace Link
DOI of Published Version
10.1029/2023ms004130
