 |
| Titel |
Full Stokes glacier model on GPU |
| VerfasserIn |
Aleksandar Licul, Frederic Herman, Yuri Podladchikov, Ludovic Räss, Samuel Omlin |
| Konferenz |
EGU General Assembly 2015
|
| Medientyp |
Artikel
|
| Sprache |
Englisch
|
| Digitales Dokument |
PDF |
| Erschienen |
In: GRA - Volume 17 (2015) |
| Datensatznummer |
250105434
|
| Publikation (Nr.) |
 EGU/EGU2015-4958.pdf |
|
|
|
|
|
| Zusammenfassung |
| Two different approaches are commonly used in glacier ice flow modeling: models based
on asymptotic approximations of ice physics and full stokes models. Lower order
models are computationally lighter but reach their limits in regions of complex
flow, while full Stokes models are more exact but computationally expansive. To
overcome this constrain, we investigate the potential of GPU acceleration in glacier
modeling.
The goal of this preliminary research is to develop a three-dimensional full Stokes
numerical model and apply it to the glacier flow. We numerically solve the nonlinear Stokes
momentum balance equations together with the incompressibility equation. Strong
nonlinearities for the ice rheology are also taken into account.
We have developed a fully three-dimensional numerical MATLAB application based on
an iterative finite difference scheme. We have ported it to C-CUDA to run it on GPUs. Our
model is benchmarked against other full Stokes solutions for all diagnostic ISMIP-HOM
experiments (Pattyn et al.,2008). The preliminary results show good agreement
with the other models. The major advantages of our programming approach are
simplicity and order 10-100 times speed-up in comparison to serial CPU version of the
code.
Future work will include some real world applications and we will implement the free
surface evolution capabilities.
References:
[1] F. Pattyn, L. Perichon, A. Aschwanden, B. Breuer, D.B. Smedt, O. Gagliardini,
G.H. Gudmundsson, R.C.A. Hindmarsh, A. Hubbard, J.V. Johnson, T. Kleiner,
Y. Konovalov, C. Martin, A.J. Payne, D. Pollard, S. Price, M. Ruckamp, F. Saito,
S. Sugiyama, S., and T. Zwinger, Benchmark experiments for higher-order and
full-Stokes ice sheet models (ISMIP-HOM), The Cryosphere, 2 (2008), 95-108. |
|
|
|
|
|
|