Heat conduction in 1D harmonic crystal: Discrete and continuum approaches

• Published in: International journal of heat and mass transfer Vol. 176; p. 121442
• Main Authors: Sokolov, Aleksei A., Müller, Wolfgang H., Porubov, Alexey V., Gavrilov, Serge N.
• Format: Journal Article
• Language: English
• Published: Oxford Elsevier Ltd 01.09.2021; Elsevier BV
• Subjects: Anomalous heat transfer; Asymptotic properties; Conduction heating; Conductive heat transfer; Crystal lattices; Differential equations; Discrete media; Energy transfer; Harmonic crystal; Low-dimensional materials; Perturbation; Smoothness; Thermal energy; Thermal processes; Thermodynamics; Thermal processes; Harmonic crystal; Low-dimensional materials; Anomalous heat transfer; Discrete media
• ISSN: 0017-9310; 1879-2189
• DOI: 10.1016/j.ijheatmasstransfer.2021.121442

•Harmonic lattice models are very attractive for investigation of heat transfer since all properties of process can be derived analytically.•Signals in the model of harmonic 1D crystal propagate with infinite velocity causing unphysical behavior.•In context of heat conduction, this artifact can be neglected on bigger scales, when the crystal can be considered as continuum. In this work the energy transfer in a one-dimensional harmonic crystal is investigated. In particular, a comparison between the discrete approach presented by Klein, Prigogine, and Hemmer with the continuum approach presented by Krivtsov is made. In the pioneering work of Klein and Prigogine the transfer of thermal energy is considered. In particular, an expression is obtained, which allows to calculate the thermal energy of each particle as a function of time. Later, Hemmer derived and used similar expressions to solve several particular problems in context of heat conduction. In the work of Krivtsov—in contrast to the discrete approach—a partial differential continuum equation is derived from the lattice dynamics of a 1D harmonic crystal. This so-called ballistic heat equation describes the propagation of heat at a finite speed in a continuous one-dimensional medium. The current work compares analyses based on the discrete equation of Klein, Prigogine, and Hemmer with those from the continuum-PDE-based one by Krivtsov. There is an important difference between the approaches. The continuum approach is derived from the dynamics of the crystal lattice, in which only kinetic degrees of freedom were excited and then thermal equilibration occurred. In contrast to that we consider in the discrete approach explicitly given equal kinetic and potential initial energies. Several exactly solvable initial problems are studied by using both methods. The problem of point perturbation shows a discrepancy in the results obtained in the framework of the continuous and discrete approaches. It is caused by the fact that the smoothness conditions of the initial perturbation is violated for the continuum approach. For other problems it is shown that at large spatial scales, where the one-dimensional crystal can be considered as a continuous medium, the discrete and the continuum relations converge. The asymptotic behavior of the difference between two aforementioned approaches is analyzed.