Towards adaptive numerical integration of dynamical contact problems

Add to your list(s) Download to your calendar using vCal

• Peter Deuflhard (Zuse Institute Berlin, Free University, Berlin and MATHEON)
• Thursday 07 May 2009, 15:00-16:00
• MR14, CMS.

If you have a question about this talk, please contact ai10.

The author and his research group in computational medicine have come across dynamical contact problems in the context of a collaboration in orthopaedic surgery. Special first attention has focussed on the motion of the patient-specific knee.

As it turned out, the numerical integration of time dependent contact problems has stayed unsatisfactory for decades. The classical Newmark method, which is quite popular in the engineering world, is a real ’’perpetuum mobile’’ in that in generates energy! A rather recent improvement due the Caltech group around Marsden and Ortiz is energy dissipative, but still unsatisfactory, since it produces artificial oscillations (untolerable in the collaboration with surgeons!). For this reason, the author and his coworkers have suggested a further modification called ’’contact-stabilized Newmark method’’. This scheme now is energy dissipative and avoids artificial oscillations. However, the new scheme escapes the usual domain of consistency theory for numerical integration. After a long investigative period, we meanwhile found the theoretical key to this kind of integrators. The new theoretical characterization requires bounded variation in terms of a physical energy functional that includes kinetic energy, elastic energy, and visco-elastic energy. First numerical findings (a few days old) for an adaptive stepsize control in a Hertzian contact problem will be presented.

This talk is part of the Applied and Computational Analysis series.

This talk is included in these lists:

• All CMS events
• All Talks (aka the CURE list)
• Applied and Computational Analysis
• CMS Events
• DAMTP info aggregator
• Featured lists
• Interested Talks
• MR14, CMS
• My seminars
• Type the title of a new list here
• bld31

Note that ex-directory lists are not shown.