Parallel performance evaluation of MIC(0) preconditioning algorithm for voxel µFE simulation

  • Authors:
  • Ivan Lirkov;Yavor Vutov;Marcin Paprzycki;Maria Ganzha

  • Affiliations:
  • Institute for Parallel Processing, Bulgarian Academy of Sciences, Sofia, Bulgaria;Institute for Parallel Processing, Bulgarian Academy of Sciences, Sofia, Bulgaria;Systems Research Institute, Polish Academy of Sciences, Warsaw, Poland;Systems Research Institute, Polish Academy of Sciences, Warsaw, Poland

  • Venue:
  • PPAM'09 Proceedings of the 8th international conference on Parallel processing and applied mathematics: Part II
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

Numerical homogenization is used for up-scaling of a linear elasticity tensor of strongly heterogeneous micro-structures. Utilized approach assumes presence of a periodic micro-structure and thus periodic boundary conditions. Rotated trilinear Rannacher-Turek finite elements are used for the discretization, while a parallel PCG method is used to solve arising large-scale systems with sparse, symmetric, positive semidefinite matrices. Applied preconditioner is based on modified incomplete Cholesky factorization MIC(0). The test problem represents a trabecular bone tissue, and takes into account only the elastic response of the solid phase. The voxel micro-structure of the bone is extracted from a high resolution computer tomography image. Numerical tests performed on parallel computers demonstrate the efficiency of the developed algorithm.