AUTHOREA
Log in
Sign Up
Browse Preprints
LOG IN
SIGN UP
Essential Site Maintenance
: Authorea-powered sites will be updated circa 15:00-17:00 Eastern on Tuesday 5 November.
There should be no interruption to normal services, but please contact us at
[email protected]
in case you face any issues.
Simone Silvestri
Public Documents
2
A GPU-based ocean dynamical core for routine mesoscale-resolving climate simulations
Simone Silvestri
and 9 more
May 22, 2024
A document by Simone Silvestri. Click on the document to view its contents.
A new WENO-based momentum advection scheme for simulations of ocean mesoscale turbule...
Simone Silvestri
and 6 more
February 26, 2024
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: (i) two-dimensional turbulence and (ii) 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.