Many-body dynamics problems are expected to handle millions of unknowns when, for instance, investigating the three-dimensional flow of granular material. Unfortunately, the size of the problems tractable by existing numerical solution techniques is severely limited on convergence grounds. This is typically the case when the equations of motion embed a differential variational inequality problem that captures contact and possibly frictional interactions between rigid and/or flexible bodies. As the size of the physical system increases, the speed and/or the quality of the numerical solution decreases. This paper describes three methods – the gradient projected minimum residual method, the preconditioned spectral projected gradient with fallback method, and the modified proportioning with reduced gradient projection method – that demonstrate better scalability than the projected Jacobi and Gauss–Seidel methods commonly used to solve contact problems that draw on a differential-variational-inequality-based modeling approach.

Using Krylov subspace and spectral methods for solving complementarity problems in many-body contact dynamics simulation / Toby, Heyn; Mihai, Anitescu; Tasora, Alessandro; Dan, Negrut. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - 95:7(2013), pp. 541-561. [10.1002/nme.4513]

Using Krylov subspace and spectral methods for solving complementarity problems in many-body contact dynamics simulation

TASORA, Alessandro;
2013-01-01

Abstract

Many-body dynamics problems are expected to handle millions of unknowns when, for instance, investigating the three-dimensional flow of granular material. Unfortunately, the size of the problems tractable by existing numerical solution techniques is severely limited on convergence grounds. This is typically the case when the equations of motion embed a differential variational inequality problem that captures contact and possibly frictional interactions between rigid and/or flexible bodies. As the size of the physical system increases, the speed and/or the quality of the numerical solution decreases. This paper describes three methods – the gradient projected minimum residual method, the preconditioned spectral projected gradient with fallback method, and the modified proportioning with reduced gradient projection method – that demonstrate better scalability than the projected Jacobi and Gauss–Seidel methods commonly used to solve contact problems that draw on a differential-variational-inequality-based modeling approach.
2013
Using Krylov subspace and spectral methods for solving complementarity problems in many-body contact dynamics simulation / Toby, Heyn; Mihai, Anitescu; Tasora, Alessandro; Dan, Negrut. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - 95:7(2013), pp. 541-561. [10.1002/nme.4513]
File in questo prodotto:
File Dimensione Formato  
IJNME_2013_krylov.pdf

non disponibili

Tipologia: Documento in Pre-print
Licenza: Creative commons
Dimensione 2.61 MB
Formato Adobe PDF
2.61 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11381/2621247
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 55
  • ???jsp.display-item.citation.isi??? 44
social impact