Discrete heat equation for a periodic layered system with allowance for the interfacial thermal resistance: General...

Discrete heat equations for the multilayered periodic systems with allowance for the thermal resistance between the layers and corresponding dispersion relations in ω-k space have been derived and analyzed. The discrete equations imply a finite velocity of thermal disturbances and guarantee the positiveness of the solutions. Analytical expressions for the attenuation distance, and phase and group velocities have been obtained as functions of frequency and thermal resistance between the discrete layers. These functions demonstrate unusual behavior at high frequency compared to the continuum case. Furthermore, the maximum allowed frequency and wave number for the discrete heat equation are limited, whereas there are no such limits in a continuum. The discrete equation contains an infinite hierarchy of continuous partial differential equations, which starts with the Fourier law, proceeds with the hyperbolic equation, the Guyer-Krumhansl (or Jeffreys type) equation, and then with higher-order equations. The partial differential equations with a finite number of terms are only approximations of the discrete equation, which implies that on the ultrashort space and timescales the discrete approach is preferable. This work provides a relatively simple, easy-to-adopt, conceptual tool, together with analytical expressions allowing one to study ultrafast wavelike heat conduction regimes in periodic multilayered metamaterials.