Accessibility navigation


Adaptively implicit MPDATA advection for arbitrary Courant numbers and Meshes

Weller, H. ORCID: https://orcid.org/0000-0003-4553-7082, Woodfield, J., Kühnlein, C. and Smolarkiewicz, P. K. (2023) Adaptively implicit MPDATA advection for arbitrary Courant numbers and Meshes. Quarterly Journal of the Royal Meteorological Society, 149 (751). pp. 369-388. ISSN 1477-870X

[img] Text (Open Access) - Published Version
· Restricted to Repository staff only
· The Copyright of this document has not been checked yet. This may affect its availability.
· Available under License Creative Commons Attribution.

9MB
[img] Text - Accepted Version
· Restricted to Repository staff only
· The Copyright of this document has not been checked yet. This may affect its availability.

21MB

It is advisable to refer to the publisher's version if you intend to cite from this work. See Guidance on citing.

To link to this item DOI: 10.1002/qj.4411

Abstract/Summary

Advection schemes with time step restrictions are widely used in weather and climate models. This can lead to instability in the presence of high flow speeds (relative to mesh spacing) such as occurs in convective updrafts, regions of mesh convergence or where the winds are unusually high. An adaptively implicit advection scheme is proposed which treats advection implicitly only where the Courant number is high. Flux-correction to ensure monotonicity is adapted to work with implicit time stepping. A version of MPDATA (Multidimensional Positive Definite Advection Transport Algorithm) is derived with an anti-diffusive flux compensating truncation errors of off-centred implicit time stepping. The anti-diffusive flux is gradually reduced as Courant numbers increase above 2 in order to maintain stability at the expense of second-order accuracy at high Courant numbers. Results of two-dimensional advection by deformational flow are presented on various meshes of the sphere. Stability and second-order accuracy are maintained when the Courant number is over 100 in a small region, when strong wind crosses the poles of a rotated latitude-longitude mesh. Good solutions are also obtained on a skipped latitude-longitude mesh, a cubed sphere and hexagonal meshes. Accuracy reverts to first-order when Courant numbers are large over a large fraction of the domain.

Item Type:Article
Refereed:Yes
Divisions:Science > School of Mathematical, Physical and Computational Sciences > Department of Meteorology
ID Code:108775
Uncontrolled Keywords:Advection, Transport, Courant number, implicit, monotonicity, stable
Publisher:Royal Meteorological Society

University Staff: Request a correction | Centaur Editors: Update this record

Page navigation