Heat transfer in a 1D harmonic crystal
A.M. Krivtsov (analytical silution, simulation algorithms), D.V. Tsvetkov (programming, calculation algorithms).
Heat transfer in the simplest discrete systems doesn’t obey common macroscopic laws. Recently experimentalists have observed the similar behavior at nanolevel, in molecular and atomic systems. The simulation below demonstrates heat transfer process in a 1D harmonic crystal. Two graphs are shown: results of molecular dynamics simulation and corresponding continuum solution. You can also compare the results with predictions of other continuum models. The analysis of the system and derivation of the continuum solution are presented in paper: ArXiv:1509.02506 (abstract, pdf). See also information below the simulation.
Use the Restart button to see the process from the beginning.
Discrete model (microlevel)
We consider a one-dimensional crystal, described by the following equations of motion:
whereis the displacement of the th particle, is the particle mass, is the stiffness of the interparticle bond. The crystal is infinite: the index is an arbitrary integer. The initial conditions are
whereare independent random values with zero expectation and unit variance; is variance of the initial velocities, which is a slowly varying function of the spatial coordinate , where is the lattice constant. These initial conditions correspond to an instantaneous temperature perturbation, which can be induced in crystals, for example, by an ultrashort laser pulse. Periodic conditions are used at the boundaries.
The kinetic temperatureis defined as
whereis the Boltzmann constant, , angle brackets stand for mathematical expectation.
Continuum description (macrolevel)
— Reversible heat wave equation: — this equation is derived as a direct consequence of the discrete microscopic equations 
Notations:is time (variable), is the sound speed.
Classic continuum equations
— Heat (Fourier): 
— Heat wave (MCV):
— Wave (d’Alembert): 
Notations:is the relaxation time (constant), is the thermal diffusivity, is the thermal conductivity, is the density, MCV stands for Maxwell-Cattaneo-Vernotte.
- A.M. Krivtsov. On unsteady heat conduction in a harmonic crystal. 2015, ArXiv:1509.02506 (abstract, pdf)
- A.M. Krivtsov. Energy oscillations in a one-dimensional crystal. Doklady Physics, 2014, Vol. 59, No. 9, pp. 427–430 (pdf: 162 Kb)
- A.M. Krivtsov. One-dimensional crystals and heat superconductivity. International Summer School-Conference “Advanced Problems in Mechanics”, 2015, St. Petersburg, Russia. Plenary lecture: pdf 1913 Kb.