Geometric numerical integration of inequality constrained, nonsmooth Hamiltonian systems

DOI10.1137/100800105MaRDI QIDQ4903737zbMATH OpenOpenAlexFDO

Authors Danny M. Kaufman, Dinesh K. Pai

Publication date 24 January 2013

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1007.2233

zbMATH Keywords

numerical experimentsHamiltonian systemsgeometric numerical integrationinequality-equality-constrained systems

Mathematics Subject Classification ID

Completely integrable finite-dimensional Hamiltonian systems, integration methods, integrability tests (37J35) Numerical methods for Hamiltonian systems including symplectic integrators (65P10) Discretization methods and integrators (symplectic, variational, geometric, etc.) for dynamical systems (37M15)

Abstract: We consider the geometric numerical integration of Hamiltonian systems subject to both equality and "hard" inequality constraints. As in the standard geometric integration setting, we target long-term structure preservation. We additionally, however, also consider invariant preservation over persistent, simultaneous and/or frequent boundary interactions. Appropriately formulating geometric methods to include such conditions has long-remained challenging due to the inherent nonsmoothness they impose. To resolve these issues we thus focus both on symplectic-momentum preserving behavior and the preservation of additional structures, unique to the inequality constrained setting. Leveraging discrete variational techniques, we construct a family of geometric numerical integration methods that not only obtain the usual desirable properties of momentum preservation, approximate energy conservation and equality constraint preservation, but also enforce multiple simultaneous inequality constraints, obtain smooth unilateral motion along constraint boundaries and allow for both nonsmooth and smooth boundary approach and exit trajectories. Numerical experiments are presented to illustrate the behavior of these methods on difficult test examples where both smooth and nonsmooth active constraint modes persist with high frequency.

Recommendations

• Stabilized Integration of Hamiltonian Systems with Hard-Sphere Inequality Constraints
• Geometric integrators for piecewise smooth Hamiltonian systems
• Symplectic Integration of Constrained Hamiltonian Systems
• Symplectic integrators for index 1 constraints
• Simulating Hamiltonian Dynamics

Cited in

(14)

• Using a half-implicit integration scheme for the SPH-based solution of fluid-solid interaction problems
• Modeling granular material dynamics and its two-way coupling with moving solid bodies using a continuum representation and the SPH method
• Discrete Lagrangian mechanics for nonseparable nonsmooth systems
• Nonholonomic systems with inequality constraints
• Half-explicit timestepping schemes on velocity level based on time-discontinuous Galerkin methods
• Geometric integrators for piecewise smooth Hamiltonian systems
• Chrono: An Open Source Multi-physics Dynamics Engine
• On the use of multibody dynamics techniques to simulate fluid dynamics and fluid-solid interaction problems
• Using Nesterov's method to accelerate multibody dynamics with friction and contact
• Accurate and efficient simulations of Hamiltonian mechanical systems with discontinuous potentials
• A regular fast multipole method for geometric numerical integrations of Hamiltonian systems
• An explicit time-integrator with singular mass for non-smooth dynamics
• Hybrid Monte Carlo methods for sampling probability measures on submanifolds
• Stabilized Integration of Hamiltonian Systems with Hard-Sphere Inequality Constraints

This page was built for publication: Geometric numerical integration of inequality constrained, nonsmooth Hamiltonian systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4903737)