A symplectic integration method for elastic filaments

Abstract

A new method is proposed for integrating the equations of motion of an elastic filament. In the standard finite-difference and finite-element formulations the continuum equations of motion are discretized in space and time, but it is then difficult to ensure that the Hamiltonian structure of the exact equations is preserved. Here we discretize the Hamiltonian itself, expressed as a line integral over the contour of the filament. This discrete representation of the continuum filament can then be integrated by one of the explicit symplectic integrators frequently used in molecular dynamics. The model systematically approximates the continuum partial differential equations, but has the same level of computational complexity as molecular dynamics and is constraint free. Numerical tests show that the algorithm is much more stable than a finite-difference formulation and can be used for high aspect ratio filaments, such as actin.

References in corpus (7)

• Non-equilibrium mechanics and dynamics of motor activated gels
• Effective medium theory of semiflexible filamentous networks
• Filamin cross-linked semiflexible networks: Fragility under strain
• Non-equilibrium hydrodynamics of a rotating filament
• Sedimentation of pairs of hydrodynamically interacting semiflexible filaments
• Stretching helical nano-springs at finite temperature
• Efficient algorithms for rigid body integration using optimized splitting methods and exact free rotational motion