Title: A scalable multi-level preconditioner for matrix-free ì-finite element analysis of human bone structures
Authors: Arbenz, P ×
van Lenthe, Harry
Mennel, U
Müller, R
Sala, M #
Issue Date: Feb-2008
Publisher: John Wiley & Sons LTD
Series Title: International Journal for Numerical Methods in Engineering vol:73 issue:7 pages:927-947
Abstract: The recent advances in microarchitectural bone imaging disclose the possibility to assess both the apparent density and the trabecular microstructure of intact bones in a single measurement. Coupling these imaging possibilities with microstructural finite element (mu FE) analysis offers a powerful tool to improve bone stiffness and strength assessment for individual fracture risk prediction.
Many elements are needed to accurately represent the intricate microarchitectural structure of bone; hence, the resulting mu FE models possess a very large number of degrees of freedom. In order to be solved quickly and reliably on state-of-the-art parallel computers, the mu FE analyses require advanced solution techniques. In this paper, we investigate the solution of the resulting systems of linear equations by the conjugate gradient algorithm, preconditioned by aggregation-based multigrid methods. We introduce a variant of the preconditioner that does not need assembling the system matrix but uses element-by-element techniques to build the multilevel hierarchy. The preconditioner exploits the voxel approach that is common in bone structure analysis, and it has modest memory requirements, at the same time robust and scalable.

Using the proposed methods, we have solved in 12 min a model of trabecular bone composed of 247 734 272 elements, yielding a matrix with 1 178 736 360 rows, using 1024 CRAY XT3 processors. The ability to solve, for the first time, large biomedical problems with over I billion degrees of freedom on a routine basis will help us improve our understanding of the influence of densitometric, morphological, and loading factors in the etiology of osteoporotic fractures such as commonly experienced at the hip, spine, and wrist. Copyright (c) 2007 John Wiley & Sons, Ltd
ISSN: 0029-5981
Publication status: published
KU Leuven publication type: IT
Appears in Collections:Biomechanics Section
× corresponding author
# (joint) last author

Files in This Item:

There are no files associated with this item.

Request a copy


All items in Lirias are protected by copyright, with all rights reserved.

© Web of science