loading page

Storm Surge Modeling as an Application of Local Time-stepping in MPAS-Ocean
  • +3
  • Jeremy Lilly,
  • Giacomo Capodaglio,
  • Mark R. Petersen,
  • Steven R Brus,
  • Darren Engwirda,
  • Robert Higdon
Jeremy Lilly
Oregon State University

Corresponding Author:[email protected]

Author Profile
Giacomo Capodaglio
Los Alamos National Lab
Author Profile
Mark R. Petersen
Los Alamos National Laboratory (DOE)
Author Profile
Steven R Brus
Los Alamos National Laboratory (DOE)
Author Profile
Darren Engwirda
Columbia University / NASA-GISS
Author Profile
Robert Higdon
Oregon State University
Author Profile

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.