Abstract
This paper presents the first scientific application of local
time-stepping (LTS) schemes in the Model for Prediction Across
Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global
ocean model that predicts the storm surge around the eastern coast of
the United States during Hurricane Sandy. The variable-resolution meshes
used are of unprecedentedly high resolution in MPAS-O, containing cells
as small as 125 meters wide in Delaware Bay. It is shown that a
particular, third-order LTS scheme (LTS3) produces sea-surface height
(SSH) solutions that are of comparable quality to solutions produced by
the classical four-stage, fourth-order Runge-Kutta method (RK4) with a
uniform time step on the same meshes. Furthermore, LTS3 is up to 35%
faster in the best cases, showing that LTS schemes are viable for use in
MPAS-O with the added benefit of substantially less computational cost.
The results of these performance experiments inform us of the
requirements for efficient mesh design for LTS schemes. In particular,
we see that for LTS to be efficient on a given mesh, it is important to
have enough cells using the coarse time-step relative to those using the
fine time-step, typically at least 1:5.