Adaptively Implicit Advection for Atmospheric Flows

Hilary Weller,Christian Kuehnlein,Piotr K. Smolarkiewicz

2024-06-13

Abstract:Implicit time-stepping for advection is applied locally in space and time where Courant numbers are large, but standard explicit time-stepping is used for the remaining solution which is typically the majority. This adaptively implicit advection scheme facilitates efficient and robust integrations with long time-steps while having negligible impact on the overall accuracy, and achieving monotonicity and local conservation on general meshes. A novel and important aspect for the efficiency of the approach is that only one linear solver iteration is needed for each advection solve. The implementation in this paper uses a second-order Runge-Kutta implicit/explicit time-stepping in combination with a second/third-order finite volume spatial discretisation. We demonstrate the adaptively implicit advection in the context of deformational flow advection on the sphere and a fully compressible model for atmospheric flows. Tracers are advected over the poles of latitude-longitude grids with very large Courant numbers and through hexagonal and cubed-sphere meshes with the same algorithm. Buoyant flow simulations with strong local updrafts also benefit from adaptively implicit advection. Stably stratified flow simulations require a stable combination of implicit treatment of gravity and acoustic waves as well as advection in order to achieve long stable time-steps.

Fluid Dynamics,Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve the problem of time - step limitations in atmospheric flow simulations, especially the instability and inefficiency caused by high Courant numbers. Specifically, the paper proposes an adaptive implicit convection scheme that applies the implicit time - stepping method locally in space and time to handle high - Courant - number regions, while other regions use the standard explicit time - stepping method. This method can significantly improve computational efficiency while maintaining stability and accuracy, allowing for simulations with longer time steps.

### Main Contributions
1. **Efficient and Stable Long - Time - Step Simulations**: By locally applying the implicit time - stepping where needed, this method can effectively handle high - Courant - number regions, thereby achieving efficient and stable simulations with long time steps.
2. **Minimizing Additional Costs**: Although implicit calculations are introduced, this method minimizes the additional computational cost by applying implicit steps only when necessary and requiring only one linear solver iteration per convection solution.
3. **Maintaining Monotonicity and Conservation**: This method can achieve monotonicity and local conservation on general grids, which is especially important for atmospheric flow simulations.
4. **Applicable to Multiple Grid Types**: This method is applicable not only to deformed - flow convection on a sphere but also to hexagonal and cubed - sphere grids, as well as to buoyancy - flow simulations that include strong local updrafts.

### Technical Details
- **Time Discretization**: Use the second - order Runge - Kutta implicit/explicit time - stepping method.
- **Space Discretization**: Combine second - order/third - order finite - volume space discretization.
- **Adaptive Strategy**: Dynamically select implicit or explicit time - stepping according to the magnitude of the Courant number.
- **Linear Stability Analysis**: Determine the value ranges of key parameters (such as α, β, γ) through von - Neumann stability analysis to ensure the stability and accuracy of the method.

### Application Examples
- **Deformed - Flow Convection on a Sphere**: Demonstrates the superior performance of this method in handling large Courant numbers on polar - rotating latitude - longitude grids.
- **Buoyancy - Driven Flow**: Verifies the advantages of implicitly handling gravity waves and sound waves.
- **Stably Stratified Flow over Topography**: Demonstrates the stability and accuracy of the method in handling complex topography.

In conclusion, this paper proposes an innovative adaptive implicit convection scheme that solves key problems in atmospheric flow simulations and improves the efficiency and stability of simulations.