SIAM Journal on Scientific and Statistical Computing - Special issue on iterative methods in numerical linear algebra
Parallel dynamic graph partitioning for adaptive unstructured meshes
Journal of Parallel and Distributed Computing - Special issue on dynamic load balancing
A Combined Evolutionary Search and Multilevel Optimisation Approach to Graph-Partitioning
Journal of Global Optimization
Developing object-oriented parallel iterative methods
International Journal of High Performance Computing and Networking
Solving coupled geoscience problems on high performance computing platforms
ICCS'05 Proceedings of the 5th international conference on Computational Science - Volume Part II
A portable parallel finite element simulation system
Proceedings of the 7th ACM international conference on Computing frontiers
Are upwind techniques in multi-phase flow models necessary?
Journal of Computational Physics
Sparse matrix and solver objects for parallel finite element simulation of multi-field problems
HPCA'09 Proceedings of the Second international conference on High Performance Computing and Applications
| Hi-index | 0.01 |
Many applied problems in geoscience require knowledge about complex interactions between multiple physical and chemical processes in the sub-surface. As a direct experimental investigation is often not possible, numerical simulation is a common approach. The numerical analysis of coupled thermo-hydro-mechanical (THM) problems is computationally very expensive, and therefore the applicability of existing codes is still limited to simplified problems. In this paper we present a novel implementation of a parallel finite element method (FEM) for the numerical analysis of coupled THM problems in porous media. The computational task of the FEM is partitioned into sub-tasks by a priori domain decomposition. The sub-tasks are assigned to the CPU nodes concurrently. Parallelization is achieved by simultaneously establishing the sub-domain mesh topology, synchronously assembling linear equation systems in sub-domains and obtaining the overall solution with a sub-domain linear solver (parallel BiCGStab method with Jacobi pre-conditioner). The present parallelization method is implemented in an object-oriented way using MPI for inter-processor communication. The parallel code was successfully tested with a 2-D example from the international DECOVALEX benchmarking project. The achieved speed-up for a 3-D extension of the test example on different computers demonstrates the advantage of the present parallel scheme.