A symplectic integration method for elastic filaments
arXiv:0811.1761 · doi:10.1063/1.3077863
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 (16)
- Nonlinear Elasticity in Biological Gels
- Non-equilibrium mechanics and dynamics of motor activated gels
- Dynamic Fluctuations of Semiflexible Polymers
- Effective medium theory of semiflexible filamentous networks
- Nonlinear dynamic intertwining of rods with self-contact
- Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
- The mechanical response of semiflexible networks to localized perturbations
- Filamin cross-linked semiflexible networks: Fragility under strain
- On micro-bead mechanics with actin filaments
- Efficient numerical integrators for stochastic models
- Non-equilibrium hydrodynamics of a rotating filament
- Sedimentation of pairs of hydrodynamically interacting semiflexible filaments
- Writhing Geometry of Open DNA
- Unfolding cross-linkers as rheology regulators in F-actin networks
- Stretching helical nano-springs at finite temperature
- Efficient algorithms for rigid body integration using optimized splitting methods and exact free rotational motion