https://github.com/alexjarosch/sia-fluxlim
Implementation of a MUSCL-superbee flux limiter to a shallow ice approximation flow code based on finite differences
https://github.com/alexjarosch/sia-fluxlim
Last synced: 2 days ago
JSON representation
Implementation of a MUSCL-superbee flux limiter to a shallow ice approximation flow code based on finite differences
- Host: GitHub
- URL: https://github.com/alexjarosch/sia-fluxlim
- Owner: alexjarosch
- License: gpl-3.0
- Created: 2016-01-12T11:03:12.000Z (over 9 years ago)
- Default Branch: master
- Last Pushed: 2019-11-09T10:48:37.000Z (over 5 years ago)
- Last Synced: 2024-11-17T00:33:36.024Z (6 months ago)
- Language: Python
- Size: 33.2 KB
- Stars: 7
- Watchers: 2
- Forks: 3
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
- cryosphere-links - sia-fluxlim - Implementation of a MUSCL-superbee flux limiter to a shallow ice approximation flow code based on finite differences (Ice Sheets and Glaciers / ISG Software)
README
# sia-fluxlim
A python implementation of a MUSCL flux limiter to a shallow ice approximation flow code based on finite differences.
For an overview of possible flux limiter functions (phi) to the MUSCL scheme, see [Wikipedia - flux limiter](https://en.wikipedia.org/wiki/Flux_limiter).
I have implemented now the Sweby (1984) limiter, which uses a parameter called 'beta'. For 1 <= beta <= 2 the limiter is second order TVD and its end members are equal to the limiter discussed in the paper below. If beta = 1, the Sweby limiter is equal to Roe's minmod limiter and if beta = 2 one gets Roe's superbee limiter.
For more details such as the scientific argument as to why this is a good idea, please consult the accompanying publication by **A. H. Jarosch**, C. G. Schoof, and F. S. Anslow, "[Restoring mass conservation to shallow ice flow models over complex terrain](http://www.the-cryosphere.net/7/229/2013/tc-7-229-2013.html),” *The Cryosphere*, Vol. 7, Iss. 1, pp. 229-240, 2013.