SWIFT: A Monotonic, Flux-Form Semi-Lagrangian Tracer Transport Scheme for Flow with Large Courant Numbers

Thomas M. Bendall,James Kent

2024-10-31

Abstract:Local conservation of mass and entropy are becoming increasingly desirable properties for modern numerical weather and climate models. This work presents a Flux-Form Semi-Lagrangian (FFSL) transport scheme, called SWIFT, that facilitates this conservation for tracer variables, whilst maintaining other vital properties such as preservation of a constant, monotonicity and positivity. Importantly, these properties all hold for large Courant numbers and multi-dimensional flow, making the scheme appropriate for use within a dynamical core which takes large time steps. The SWIFT scheme presented here can be seen as an evolution of the FFSL methods of Leonard et al and Lin and Rood. Two-dimensional and three-dimensional schemes consist of a splitting into a sequence of one-dimensional calculations. The new SWIFT splitting presented here allows monotonic and positivity properties from the one-dimensional calculations to be inherited by the multi-dimensional scheme. These one-dimensional calculations involve separating the mass flux into terms that correspond to integer and fractional parts of the Courant number. Key to achieving conservation is coupling the transport of tracers to the transport of the fluid density, through re-use of the discrete mass flux that was calculated from the fluid density in the transport of the tracers. This work also describes how these properties can still be attained when the tracer is vertically-staggered from the density in a Charney-Phillips grid.

Numerical Analysis,Atmospheric and Oceanic Physics,Fluid Dynamics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a new transport scheme to meet the requirements of local mass conservation and entropy conservation in modern numerical weather and climate models. Specifically, the paper introduces a flux - form semi - Lagrangian (FFSL) transport scheme called SWIFT, aiming to solve the following problems:

1. **Maintaining local mass conservation**: Ensure the mass conservation of the tracer variables under large Courant numbers and multi - dimensional flow conditions.
2. **Maintaining monotonicity and positivity**: Ensure that the transport scheme can still maintain the monotonicity and positivity of the solution under large Courant numbers, preventing the occurrence of unphysical negative concentrations or new extreme points.
3. **Applicability to the atmospheric dynamical core**: Be able to maintain computational efficiency when using large time steps without sacrificing other required properties such as consistency and stability.

### Specific problem description

#### 1. Local mass conservation
Modern numerical weather and climate models increasingly require the property of local mass conservation. This means that the mass budget of each grid cell should be closed, that is, the inflow and outflow of mass should be balanced. This is very important to ensure that the model does not become a net source or sink of the tracer substances.

#### 2. Monotonicity and positivity
- **Monotonicity**: Prevent the introduction of new minimum or maximum values and ensure that the shape of the solution remains unchanged.
- **Positivity**: Ensure that no unphysical negative concentrations are generated, which is crucial for ensuring the physical rationality of the solution.

#### 3. Stability under large Courant numbers
Many operational numerical weather prediction models use relatively large time steps. Therefore, an ideal transport scheme should be able to maintain stability under large Courant numbers without affecting computational efficiency and other required properties.

#### 4. Compatibility with the atmospheric dynamical core
The transport scheme must adapt to the parameter settings of the atmospheric dynamical core, such as using hexahedral cells, C - grid staggered grids in the horizontal direction and Charney - Phillips staggered grids in the vertical direction. These settings are very important for a good representation of gravity wave dispersion.

### Solution overview

The SWIFT scheme solves the above problems in the following ways:
- **One - dimensional to multi - dimensional splitting**: Decompose the multi - dimensional transport problem into a series of one - dimensional calculations to ensure that monotonicity and positivity are inherited from one - dimensional calculations to the multi - dimensional scheme.
- **Separating the integer and fractional parts of the Courant number**: Ensure consistency and stability under large Courant numbers by separating the integer and fractional parts in the mass flux.
- **Re - using the discrete mass flux**: Transport the tracer variables by re - using the discrete mass flux calculated when transporting the density to ensure the consistency of the transport.

In summary, the SWIFT scheme aims to provide an efficient transport method that can maintain mass conservation, monotonicity and positivity under large Courant numbers and multi - dimensional flow conditions, and is suitable for modern numerical weather and climate models.