Alya: Computational Solid Mechanics for Supercomputers

被引:23
作者
Casoni, E. [1 ]
Jerusalem, A. [3 ,6 ]
Samaniego, C. [1 ]
Eguzkitza, B. [1 ]
Lafortune, P. [5 ]
Tjahjanto, D. D. [3 ,4 ]
Saez, X. [1 ]
Houzeaux, G. [1 ]
Vazquez, M. [1 ,2 ]
机构
[1] CNS, Barcelona Supercomputing Ctr, Barcelona 08034, Spain
[2] CSIC, Artificial Intelligence Res Inst, IIIA, Bellaterra 08193, Spain
[3] Tecnogetafe, IMDEA Mat Inst, Getafe 28906, Spain
[4] KTH Royal Inst Technol, Dept Solid Mech, S-10044 Stockholm, Sweden
[5] Idra Simulat SCP, Barcelona 08004, Spain
[6] Univ Oxford, Dept Engn Sci, Oxford OX1 3PJ, England
关键词
Computational mechanics; Finite element; method; Parallel computing; Chimera; DOMAIN DECOMPOSITION METHODS; DISCONTINUOUS GALERKIN; ITERATIVE SOLVERS; PARALLEL; PRECONDITIONER; EFFICIENT; FRACTURE; DISPLACEMENT; ALGORITHMS; HYBRID;
D O I
10.1007/s11831-014-9126-8
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
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, combining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multi-physics 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.
引用
收藏
页码:557 / 576
页数:20
相关论文
共 74 条
[1]  
Adams M, 1999, UCBCSD991036 EECS DE
[2]   Algebraic multigrid methods for direct frequency response analyses in solid mechanics [J].
Adams, Mark F. .
COMPUTATIONAL MECHANICS, 2007, 39 (04) :497-507
[3]   Multifrontal parallel distributed symmetric and unsymmetric solvers [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2000, 184 (2-4) :501-520
[4]  
[Anonymous], P C HIGH PERF COMP N
[5]   A scalable multi-level preconditioner for matrix-free μ-finite element analysis of human bone structures [J].
Arbenz, Peter ;
van Lenthe, G. Harry ;
Mennel, Uche ;
Mueller, Ralph ;
Sala, Marzio .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2008, 73 (07) :927-947
[6]   A HIGHLY SCALABLE PARALLEL IMPLEMENTATION OF BALANCING DOMAIN DECOMPOSITION BY CONSTRAINTS [J].
Badia, Santiago ;
Martin, Alberto F. ;
Principe, Javier .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2014, 36 (02) :C190-C218
[7]  
Balay S., 2013, ANL9511
[8]   A full-discontinuous Galerkin formulation of nonlinear Kirchhoff-Love shells: elasto-plastic finite deformations, parallel computation, and fracture applications [J].
Becker, G. ;
Noels, L. .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2013, 93 (01) :80-117
[9]  
Benek JA, 1986, TECHNICAL REPORT
[10]  
Bhardwaj M, 2002, SUPERCOMPUTING, P1