Alya: Computational Solid Mechanics for Supercomputers

Abstract

While solid mechanics codes are now conventional tools both in industry and research, the increasingly more exigent requirements of both sectors are fuelling the need for more computational power and more advanced algorithms. For obvious reasons, commercial codes are lagging behind academic codes often dedicated either to the implementation of one new technique, or the upscaling of current conventional codes to tackle massively large scale computational problems. Only in a few cases, both approaches have been followed simultaneously. In this article, a solid mechanics simulation strategy for parallel supercomputers based on a hybrid approach is presented. Hybrid parallelization exploits the thread-level parallelism of multicore architectures, com- bining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multiphysics code. Hybrid parallelization is specially well suited for the current trend of supercomputers, namely large clusters of multicores. The strategy is assessed through transient non-linear solid mechanics problems, both for explicit and implicit schemes, running on thousands of cores. In order to demonstrate the flexibility of the proposed strategy under advance algorithmic evolution of computational mechanics, a non-local parallel overset meshes method (Chimera-like) is implemented and the conservation of the scalability is demonstrated.

Analysis of hemodynamics and wall mechanics at sites of cerebral aneurysm rupture

“…The maximum eigenvalue of the resulting stress tensor was calculated and averaged along the vessel thickness and used as a measure of the local wall stress. All CSD simulations were carried out in parallel using an in-house finite element structural solver 13. Each CSD run took less than 1 h on 128 processors of a distributed memory system equipped with 3056 nodes composed of two eight-core Intel SandyBridge-EP 2.6 GHz processors.…”

Domain Decomposition Methods for Domain Composition Purpose: Chimera, Overset, Gluing and Sliding Mesh Methods

“…The method was applied to the Navier-Stokes equations as a gluing mesh strategy in [35], as a Chimera method in [52], and to solid mechanics problems in [20]. Discretization error is nodally zero if the exact solution belongs to the finite element space and the L 2 norm of the solution error is consistent with that of the original scheme.…”

Fully coupled fluid‐electro‐mechanical model of the human heart for supercomputers

“…The four problems share some common features, described below. All of them are implemented and solved in Alya, the BSC's in‐house parallel simulation code for coupled multi‐physics problems, which is designed to run efficiently on high performance computers . The four problems are discretized in space using the finite element method (FEM) on nonstructured meshes.…”