Thermal conduction in the Fermi-Pasta-Ulam model

Download files
Access & Terms of Use
open access
Copyright: Tempatarachoke, Pisut
Altmetric
Abstract
We conduct a comprehensive and systematic study of the Fermi-Pasta-Ulam (FPU) model using both equilibrium and non-equilibrium molecular dynamics simulations, with the aim being to explain the cause of the anomalous energy-transport behaviour in the model. In the equilibrium scenario, our motivation stems from the lack of a complete understanding of the effects of initial conditions on the energy dissipation among Fourier modes. We also critically reconsider the ¡¥probes' that had been widely used to quantitatively describe the types of energy sharing in a system, and then decide on a preferred choice to be used in our equilibrium study. We establish, from strong numerical evidence, that there exists a critical energy density of approximately 0:1, above which the energy dissipation among the modes becomes independent of initial conditions and system parameters, and that the full equipartition of mode energy is never attained in the FPU model. We report, for the first time, the violation of particle positions in the FPU model at high energies, where the particles are found to pass through one another. In the non-equilibrium scenario, we critically review the Nos„Se-Hoover algorithm thermostatting method largely used by other works, and identify its weaknesses. We also review some other alternative methods and decide on the most appropriate one to be implemented throughout our work. We confirm the divergence of the thermal conductivity of the FPU model as the chain length increases, and that kfpu [symbol] No.41, in agreement with other works. Our study further shows that there exists an upper limit of the anharmonicity in the FPU model, and that any attempt to increase the strength of this anharmonicity will not succeed. We also introduce elastic collisions into the original FPU model and find that the Modified model (FPUC) still exhibits anomalous thermal conductivity. We conclude that a one-dimensional FPU-type model with ¡¥only' nearest-neighbour interaction, regardless of being soft or hard, does not exhibit a finite thermal conductivity as the system size increases, due to the non-chaotic nature of its microscopic dynamics, the origin of which we are unable to account for. Finally, we briefly outline possible research directions.
Persistent link to this record
Link to Publisher Version
Link to Open Access Version
Additional Link
Author(s)
Tempatarachoke, Pisut
Supervisor(s)
Creator(s)
Editor(s)
Translator(s)
Curator(s)
Designer(s)
Arranger(s)
Composer(s)
Recordist(s)
Conference Proceedings Editor(s)
Other Contributor(s)
Corporate/Industry Contributor(s)
Publication Year
2005
Resource Type
Thesis
Degree Type
PhD Doctorate
UNSW Faculty
Files
download whole.pdf 41 MB Adobe Portable Document Format
Related dataset(s)