• search hit 2 of 2
Back to Result List

An Inexact Interior Point Method for the Large-Scale Simulation of Granular Material

  • Non-smooth contact dynamics provides an increasingly popular simulation framework for granular material. In contrast to classical discrete element methods, this approach is stable for arbitrary time steps and produces visually acceptable results in very short computing time. Yet when it comes to the prediction of draft forces, non-smooth contact dynamics is typically not accurate enough. We therefore propose to combine the method class with an interior point algorithm for higher accuracy. Our specific algorithm is based on so-called Jordan algebras and exploits the relation to symmetric cones in order to tackle the conical constraints that are intrinsic to frictional contact problems. In every interior point iteration a linear system has to be solved. We analyze how the interior point method behaves when it is combined with Krylov subspace solvers and incomplete factorizations. We show that efficient preconditioners and efficient linear solvers are essential for the method to be applicable to large-scale problems. Using BiCGstab as a linear solver and incomplete Cholesky factorizations, we substantially improve the accuracy in comparison to the projected Gauss-Jacobi solver.

Download full text files

Export metadata

Author:Jan Kleinert, Bernd Simeon, Martin Obermayr
URN (permanent link):urn:nbn:de:hbz:386-kluedo-36944
Document Type:Preprint
Language of publication:English
Publication Date:2014/01/28
Year of Publication:2013
Publishing Institute:Technische Universität Kaiserslautern
Date of the Publication (Server):2014/01/29
Number of page:21
Faculties / Organisational entities:Fachbereich Mathematik
DDC-Cassification:5 Naturwissenschaften und Mathematik / 510 Mathematik
MSC-Classification (mathematics):65-XX NUMERICAL ANALYSIS
Licence (German):Standard gemäß KLUEDO-Leitlinien vom 10.09.2012