Parallel Computing in Geosciences
UPPMAX project no. P2004009, 2004 - 2006
The project represents multi-disciplinary research involving mathematical modelling of problems arising in geophysics and geo-engineering applications, numerical solution methods with an emphasis on preconditioned iterative techniques and their efficient implementation on parallel computers. The geosciences and geo-engineering nowadays consider a variety of problems, which demand the use of sophisticated large-scale computer modelling. The target problems with this project include modelling of stress field changes due to mining and underground constructions, assessment of underground nuclear waste repository projects as well as elastic and viscoelastic effects of the Earth due to other glacial advance and recession.
The very high complexity of these problems is usually determined by the required models. As an example, the mechanical behaviour of rocks can be described by models ranging from linear elasticity to models involving nonlinear material response, permanent deformation, rheology phenomena or influence of thermal changes and flow in porous media. The large-scale arises when we consider large 3D domains to enable monitoring of the far field effects as well as to be able to specify realistic boundary conditions. Simultaneously, some engineering decisions are based on a detailed analysis of the near field processes.
The solution of the above described problems then requires the use of robust and numerically efficient iterative solution methods and their implementation on modern parallel computer architectures. This project addresses both discretization of the mathematical models by the finite element method and solution of the discretized problems by robust and efficient iterative methods using domain decomposition ideas combined with hierarchical substructuring approaches.
The considered numerical methods include Krylov subspace methods preconditioned by Schwarz-type decomposition methods accelerated by an auxiliary coarse grid problem. Such methods are applied to a hierarchy of problems including elasticity, thermo-elasticity, visco-elasticity and poro-elasticity.
Aims of project
- To implement the developed methods on various types of computer architectures using different programming paradigms (MPI and OpenMP) and to improve the overall parallel performance of the methods to insure short simulation times.
- To compare the scalability and performance of the solvers on various computer platforms.
- development of stresses due to mining.
Large-scale 3D geotechnical model, Dolni Rozinka, the Czech Republic.
- modelling thermo-mechanical phenomena.
Large-scale 3D geoenvironmental model, Aspo Prototype Repository, Sweden.
- Lomond - Sun HPC cluster, EPCC, Edinburgh
- Natan - IBM symmetric multiprocessor, IGAS, Ostrava
- Thea - Beowulf cluster, IGAS, Ostrava
- Termit - Beowulf cluster, VSB-TU, Ostrava
- Simba - Sun symmetric multiprocessor, UPPMAX, Uppsala
- Ra - Sun cluster, UPPMAX, Uppsala
Reports and presentations
- Scientific Computing Group
Department of Information Technology
Uppsala University, Sweden
- Department of Applied Mathematics and Computer Science
Institute of Geonics AS CR
Ostrava, Czech Republic
- Scientific Computing Department
Institute for Parallel Processing BAS
- Research Centre for Advanced Remediation Technologies
Technical University of Liberec, Czech Republic