Timothy C. Andrews,Jemma Shipton,Beth A. Wingate
Abstract:For simulations of time evolution problems, such as weather and climate models, taking the largest stable timestep is advantageous for reducing wall-clock time. A drawback of doing so is the potential reduction in nonlinear accuracy of the numerical solution - we investigate this for the Rotating Shallow Water Equations (RSWEs) on an f-plane. First, we examine how linear dispersion errors can impact the nonlinear dynamics. By deriving an alternate time evolution equation for the RSWEs, the dynamics can be expressed through interactions of three linear waves in triads. Linear dispersion errors may appear in the numerical representation of the frequency of each triad, which will impact the timestepped nonlinear dynamics. A new triadic error quantifies this by composing three stability polynomials from the oscillatory Dahlquist test equation. Second, we design two new test cases to examine the effect of timestep size in a numerical model. These tests investigate how well a timestepper replicates slow nonlinear dynamics amidst fast linear oscillations. The first test case of a Gaussian height perturbation contains a nonlinear phase shift that can be missed with a large timestep. The second set of triadic test cases excite a few linear waves to instigate specific triadic interactions. Two triadic cases are examined: one with a dominant directly resonant triad and another with near-resonances that redistribute fast mode energy into rings in spectral space. Three numerical models, including LFRic from the Met Office, are examined in these test cases with different timesteppers.
What problem does this paper attempt to address?
The problems that this paper attempts to solve are: when performing numerical simulations with relatively large stable timesteps, how does the linear dispersion error affect the accuracy of nonlinear time - discretization, especially in the Rotating Shallow Water Equations (RSWEs) on the f - plane. Specifically, the author focuses on:
1. **The influence of linear dispersion error on nonlinear dynamics**:
- By deriving the time - evolution equations of RSWEs, the author represents the dynamics as the interactions of three linear waves in triads.
- Linear dispersion errors may occur in the numerical representation of the frequencies of each triad interaction, which will affect the nonlinear dynamics after time - discretization.
2. **Design new test cases to evaluate the influence of timesteps**:
- The author designs two new test cases to study the performance of relatively large timesteps in numerical models.
- The first test case is a Gaussian height perturbation, which contains a slow phase shift due to nonlinearity, and this shift may be ignored by relatively large timesteps.
- The second test case is a series of triad interactions, exciting specific triad combinations, including directly resonant and near - resonant triads.
3. **Introduce triadic error as a quantitative index**:
- The author proposes a new quantitative index - triadic error, which is used to measure the influence of linear dispersion error on the nonlinear RSWEs dynamics.
- Through the oscillatory Dahlquist test equation, the author constructs the stability polynomial and calculates the triadic error.
4. **Evaluate the performance of different time - discretization methods**:
- The author compares the triadic errors of classical time - discretization methods and shows the differences between these errors and linear dispersion analysis.
- Through these analyses, the author hopes to find time - discretization methods that can maintain sufficient accuracy with relatively large timesteps.
### Formula summary
- **Triad frequency**:
\[
\Omega_{\alpha, \alpha_a, \alpha_b}^{k, k_a, k_b} = \omega_{\alpha_a}^{k_a} + \omega_{\alpha_b}^{k_b} - \omega_{\alpha}^k
\]
This formula represents the frequency of triad interactions, which is composed of the frequencies of three linear waves.
- **Triad propagator**:
\[
T(\Omega_{\alpha, \alpha_a, \alpha_b}^{k, k_a, k_b}, t) = e^{i \Omega_{\alpha, \alpha_a, \alpha_b}^{k, k_a, k_b}t}
\]
This formula represents the time - evolution operator of triad interactions.
- **Triad error**:
\[
E(k, k_a, k_b, \alpha, \alpha_a, \alpha_b, \Delta t) = \left| T_{\Delta t}(\Omega_{\alpha, \alpha_a, \alpha_b}^{k, k_a, k_b}, \Delta t)-T_N(W_{\alpha_a}^{k_a}, W_{\alpha_b}^{k_b}, W_{\alpha}^k) \right|
\]
This formula defines the triad error, that is, the difference between the truly discretized triad propagator and the numerical approximation.
Through the exploration of these problems, the author aims to provide more efficient and accurate time - discretization methods for meteorological and climate simulations, especially when dealing with multi - scale dynamics problems.