Exponential integrators for stiff elastodynamic problems

We investigate the application of exponential integrators to stiff elastodynamic problems governed by second-order differential equations. Classical explicit numerical integration schemes have the shortcoming that the stepsizes are limited by the highest frequency that occurs within the solution spectrum of the governing equations, while implicit methods suffer from an inevitable and mostly uncontrollable artificial viscosity that often leads to a nonphysical behavior. In order to overcome these specific detriments, we devise an appropriate class of exponential integrators that solve the stiff part of the governing equations of motion by employing a closed-form solution. As a consequence, we are able to handle up to three orders of magnitude larger time-steps as with conventional implicit integrators and at the same time achieve a tremendous increase in the overall long-term stability due to a strict energy conservation. The advantageous behavior of our approach is demonstrated on a broad spectrum of complex deformable models like fibers, textiles, and solids, including collision response, friction, and damping.

[1]  Andrew Nealen,et al.  Physically Based Deformable Models in Computer Graphics , 2006, Comput. Graph. Forum.

[2]  P. Deuflhard.  A study of extrapolation methods based on multistep schemes without parasitic solutions , 1979 .

[3]  Ronald Fedkiw,et al.  Energy Conservation for the Simulation of Deformable Bodies , 2012 .

[4]  W. Arnoldi.  The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[5]  Michael Griebel,et al.  Numerical Simulation in Molecular Dynamics: Numerics, Algorithms, Parallelization, Applications , 2007 .

[6]  John E. Pearson,et al.  Solving Stiff Differential Equations with the Method of Patches , 1998 .

[7]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[8]  Wolfgang Straßer,et al.  A fast finite element solution for cloth modelling , 2003, 11th Pacific Conference onComputer Graphics and Applications, 2003. Proceedings..

[9]  John C. Platt,et al.  Elastically deformable models , 1987, SIGGRAPH.

[10]  S. Antman.  Nonlinear problems of elasticity , 1994 .

[11]  Reinhard Klein,et al.  Efficient Non-Linear Editing of Large Point Clouds , 2010 .

[12]  S. P. Nørsett.  An A-stable modification of the Adams-Bashforth methods , 1969 .

[13]  John Hart,et al.  ACM Transactions on Graphics , 2004, SIGGRAPH 2004.

[14]  Z. Popovic,et al.  Crystal structure of a monomeric retroviral protease solved by protein folding game players , 2011, Nature Structural &Molecular Biology.

[15]  Eitan Grinspun,et al.  Implicit-Explicit Variational Integration of Highly Oscillatory Problems , 2008, Multiscale Model. Simul..

[16]  Leonard McMillan,et al.  Stable real-time deformations , 2002, SCA '02.

[17]  C F Curtiss,et al.  Integration of Stiff Equations. , 1952, Proceedings of the National Academy of Sciences of the United States of America.

[18]  I. N. Sneddon,et al.  The Solution of Ordinary Differential Equations , 1987 .

[19]  David A. Pope.  An exponential method of numerical integration of ordinary differential equations , 1963, CACM.

[20]  M. Hochbruck,et al.  Exponential integrators , 2010, Acta Numerica.

[21]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[22]  Eitan Grinspun,et al.  Discrete elastic rods , 2008, ACM Trans. Graph..

[23]  Raanan Fattal,et al.  Efficient simulation of inextensible cloth , 2007, SIGGRAPH 2007.

[24]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[25]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[26]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[27]  Gerald Wempner,et al.  Finite elements, finite rotations and small strains of flexible shells , 1969 .

[28]  Ronald Fedkiw,et al.  Invertible finite elements for robust simulation of large deformation , 2004, SCA '04.

[29]  Marlis Hochbruck,et al.  Error analysis of exponential integrators for oscillatory second-order differential equations , 2006 .

[30]  W. Gautschi.  Numerical integration of ordinary differential equations based on trigonometric polynomials , 1961 .

[31]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[32]  J. D. Lawson.  Generalized Runge-Kutta Processes for Stable Systems with Large Lipschitz Constants , 1967 .

[33]  Ernst Hairer,et al.  Long-Time Energy Conservation of Numerical Methods for Oscillatory Differential Equations , 2000, SIAM J. Numer. Anal..

[34]  Robert D. Skeel,et al.  Long-Time-Step Methods for Oscillatory Differential Equations , 1998, SIAM J. Sci. Comput..

[35]  Bernd Eberhardt,et al.  Implicit-Explicit Schemes for Fast Animation with Particle Systems , 2000, Computer Animation and Simulation.

[36]  T. Belytschko,et al.  Applications of higher order corotational stretch theories to nonlinear finite element analysis , 1979 .

[37]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[38]  Olaf Etzmuß,et al.  A High Performance Solver for the Animation of Deformable Objects using Advanced Numerical Methods , 2001, Comput. Graph. Forum.

[39]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[40]  Joseph Hersch.  Contribution à la méthode des équations aux différences , 1958 .

[41]  Mark C. Surles.  An algorithm with linear complexity for interactive, physically-based modeling of large proteins , 1992, SIGGRAPH.

[42]  Andrew P. Witkin,et al.  Large steps in cloth simulation , 1998, SIGGRAPH.

[43]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[44]  Marlis Hochbruck,et al.  A Gautschi-type method for oscillatory second-order differential equations , 1999, Numerische Mathematik.