PULSEDYN - A dynamical simulation tool for studying strongly nonlinear chains
Abstract
We introduce PULSEDYN, a particle dynamics program in C++, to solve many-body nonlinear systems in one dimension. PULSEDYN is designed to make computing accessible to non-specialists in the field of nonlinear dynamics of many-body systems and to ensure transparency and easy benchmarking of numerical results for an integrable model (Toda chain) and three non-integrable models (Fermi-Pasta-Ulam-Tsingou, Morse and Lennard-Jones). To achieve the latter, we have made our code open source and free to distribute. We examine (i) soliton propagation and two-soliton collision in the Toda system, (ii) the recurrence phenomenon in the Fermi-Pasta-Ulam-Tsingou system and the decay of a single localized nonlinear excitation in the same system through quasi-equilibrium to an equipartitioned state, and SW propagation in chains with (iii) Morse and (iv) Lennard-Jones potentials. We recover well known results from theory and other numerical results in the literature. We have obtained these results by setting up a parameter file interface which allows the code to be used as a black box. Therefore, we anticipate that the code would prove useful to students and non-specialists. At the same time, PULSEDYN provides scientifically accurate simulations thus making the study of rich dynamical processes broadly accessible.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.