Modernizing Titan2D, a Parallel AMR Geophysical Flow Code to Support Multiple Rheologies and Extendability

In this work, we report on strategies and results of our initial approach for modernization of Titan2D code. Titan2D is a geophysical mass flow simulation code designed for modeling of volcanic flows, debris avalanches and landslides over a realistic terrain model. It solves an underlying hyperbolic...

Full description

Saved in:
Bibliographic Details
Published inarXiv.org
Main Authors Simakov, Nikolay A, Jones-Ivey, Renette L, Akhavan-Safaei, Ali, Aghakhani, Hossein, Jones, Matthew D, Patra, Abani K
Format Paper
LanguageEnglish
Published Ithaca Cornell University Library, arXiv.org 02.10.2019
Subjects
Adaptive systems
Arrays
Avalanches
Computer simulation
Data storage
Finite element method
Flow simulation
Geophysics
Hyperbolic systems
Landslides
Layouts
Mass flow
Modernization
Parallel processing
Partial differential equations
Performance indices
Rheological properties
Search algorithms
Terrain models
Online AccessGet full text

Cover

More Information
Summary:In this work, we report on strategies and results of our initial approach for modernization of Titan2D code. Titan2D is a geophysical mass flow simulation code designed for modeling of volcanic flows, debris avalanches and landslides over a realistic terrain model. It solves an underlying hyperbolic system of partial differential equations using parallel adaptive mesh Godunov scheme. The following work was done during code refactoring and modernization. To facilitate user input two level python interface was developed. Such design permits large changes in C++ and Python low-level while maintaining stable high-level interface exposed to the end user. Multiple diverged forks implementing different material models were merged back together. Data storage layout was changed from a linked list of structures to a structure of arrays representation for better memory access and in preparation for further work on better utilization of vectorized instruction. Existing MPI parallelization was augmented with OpenMP parallelization. The performance of a hash table used to store mesh elements and nodes references was improved by switching from a linked list for overflow entries to dynamic arrays allowing the implementation of the binary search algorithm. The introduction of the new data layout made possible to reduce the number of hash table look-ups by replacing them with direct use of indexes from the storage class. The modifications lead to 8-9 times performance improvement for serial execution.
ISSN:2331-8422