Ballistic-Diffusive Model for Heat Transport in Superlattices and the Minimum Effective Heat Conductivity

There has been much interest in semiconductor superlattices because of their low thermal conductivities. This makes them especially suitable for applications in a variety of devices for the thermoelectric generation of energy, heat control at the nanometric length scale, etc. Recent experiments have confirmed that the effective thermal conductivity of superlattices at room temperature have a minimum for very short periods (in the order of nanometers) as some kinetic calculations had anticipated previously. This work will show advances on a thermodynamic theory of heat transport in nanometric 1D multilayer systems by considering the separation of ballistic and diffusive heat fluxes, which are both described by Guyer-Krumhansl constitutive equations. The dispersion relations, as derived from the ballistic and diffusive heat transport equations, are used to derive an effective heat conductivity of the superlattice and to explain the minimum of the effective thermal conductivity.


Introduction
This work was motivated by the increasing interest in heat transport in superlattices due to their potential applications in a variety of devices, e.g., for the thermoelectric generation of energy, the design of intelligent coatings for temperature conditioning of enclosures, the construction of devices for processing information, the storage of thermal energy, and design of biomedical devices among others [1]. Although the applications of superlattices started several years ago, the nature of the heat transport regime that takes place in them is still being discussed. Particularly, the description of the effects of wave properties of phonons when the film thickness approaches the mean attenuation length of heat carriers continues receiving theoretical and experimental attention [2][3][4][5][6][7].
The heat transport at the nanoscale deviates from that predicted by the Fourier law, which is able to describe only diffusive transport. When the thickness of a heat conducting layer approaches the nanoscale, other non-Fourier types of transport appear, such as ballistic and wave-like [8][9][10]. The first arises when the heat carriers move without dispersion between the walls of the material and do not interact with the bulk. The second is a result of the coherent coupling of heat carriers. It arises when there are no dispersive processes that eliminate their phase information.
The diffusive transport and wave propagation (ballistic/second sound) had been experimentally detected only in cold crystals, liquid He, etc., where both heat transport regimes coexist [11][12][13][14][15]. However, since several years ago, some experimental groups have measured the effective thermal conductivity of superlattices [2][3][4][5] finding clues to the wave properties of phonons in superlattices at environment temperature. Both transport regimes have also been studied in nanostructured heat sources [6].
Saha and collaborators [2] in epitaxial metal/semiconductor superlattices TiN/(Al, Sc), N with constant thickness increased the density of periods and detected a minimum in the thermal conductivity in thicknesses of the order of the wavelength of the heat carriers. The authors used thermo-reflection in the time domain (TDTR) to measure the thermal conductivity of the samples of a TiN/AlScN superlattice. From the wave point of view, the explanation is simple [2]: "For very short periods, when the wavelengths of the dominant phonons involved in heat transport become comparable with the thickness of the period, the individual layers lose their identity and the whole material behaves as an effective medium where the waveforms of the phonon propagate without dispersion in the interfaces, increasing the thermal conductivity". The results by Saha et al. showed a clear minimum (at 4 nm of thickness of period) in the thermal conductivity. In shorter period thicknesses (less than 4 nm), the thermal conductivity increased when the thickness of the period decreased since the phononic modes that do not disperse in the interface contribute significantly to thermal transport. The findings of Saha's group reinforced those obtained some years before by Ravichandran [4] on the existence of a minimum of the effective heat conductivity in superlattices.
In another experimental approach, in measurements at room temperature of effective thermal conductivity as a function of the number of periods in GaAs/AlAs superlattices, keeping constant the thickness of the period [5], it was also found the wave regime of heat propagation. Each internal interface scatters the phonons in a diffuse way, eliminating their phase information. Then the interface behaves as a thermal resistance in such a way that many equivalent interfaces in the series lead to an effective thermal conductivity, in the direction perpendicular to the interfaces, which is approximately independent of the number of layers. However, if the phononic phase is conserved in the interfaces of the superlattice and if the anharmonic dispersion is minimal, the superposition of the Bloch waves creates stop bands and modifies the phononic band structure. In this regime, the mean free paths of the phonons (MFP) are equal to the length of the sample, which leads to a thermal conductivity that is linearly proportional to the total thickness of the superlattice. In addition, the MFP of the phonons is equal to the length of the sample, which leads to a thermal conductivity linearly proportional to the total thickness of the superlattice. The authors carried out measurements of thermal conductivity by time domain thermo-reflection (TDTR) in samples at temperatures between 30 K and 300 K. Below 150 K, linearity of thermal conductivity versus length (number of periods) revealed that heat conduction was of the wave type.
Some microscopic theories about heat transport, based on the Boltzmann equation, treat heat carriers (phonons and electrons) as particles, but ignore their wave properties. All of them predict that the thermal conductivity perpendicular to the layers in a superlattice decreases as the space between layers is reduced. However, as it may be concluded from the experiments mentioned above, as the thickness of the period becomes smaller than the wavelength of the phonon, the thermal conductivity increases. This disagreement is solved by calculations that include the MFP of the phonons. For thinner layers than the MFP, wave theory is applied. While for layers thicker than the MFP, the theory of particles is applied. The combined theory predicts a minimum of thermal conductivity for a certain value of thickness which depends on the average MFP, and therefore on the temperature [16].
In Simkin's work [17], an imaginary part of the form i/l (l is the phononic MFP) is added to the wave number to obtain a minimum in the effective thermal conductivity accordingly with Saha's experiments. The argument used by Simkin is striking: "Phenomenologically, an imaginary part of the wave number of the form i/l can be introduced for ..." No other justification is given. Macroscopically, Entropy 2020, 22, 167 3 of 15 the effective thermal conductivity seems to be determined solely by the steady state which depends on the boundary conditions, the transport coefficients and the interfacial thermal resistance.
The macroscopic models of heat transfer can predict the general form of heat propagation and also can connect the different levels of the hierarchical material behavior [18]. The hyperbolic models with relaxation type dissipation (Cattaneo typer) predict a wave behavior of heat but they are not able to incorporate emerging processes related to the coherent behavior of phonons and also the quantitative modelling of ballistic propagation can be problematic [19]. These aspects of the dynamics of non-Fourier heat propagation and in particular the form of the evolution equations determines the stationary processes, as well. Guyer-Krumhansl type equations have been derived from kinetic theory and shown their validity at room temperatures [20].
The analysis using Fourier propagation modes can shed light on the interaction of heat waves with the microstructure. It can be shown that the long wavelength modes have an attenuation length that tends to infinity when the wave number tends to zero. The larger the wavelength, the larger the penetration of the wave in the material, in such a way that the long wavelength modes do not interact with the interfaces. If we add to this the fact that long wavelength modes transport most of the heat at the nanoscale it should be possible to explain the behavior of thermal conductivity in superlattices. These issues were addressed in this work.

Two Temperature Ballistic-Diffusive Model
We propose to use irreversible thermodynamics with two vectorial dissipative fluxes (internal variables) to incorporate features of the wave-like behavior of the heat carriers and to describe transient processes associated with the results reported by the experiments above mentioned on superlattices. The theory of internal variables is easily applicable to the multilayer system and we hope that it may open new possibilities for analysis.
Ballistic-diffusive models for heat transport are based on the assumption that there exist two kinds of heat carriers in the system. Those which propagate diffusively that are scattered in the core of the material, and those which propagate ballistically, interacting only with the boundaries or the interfaces of the system. This assumption has also been made in other systems where there are two or more types of heat carriers. For example, in thermoelectric materials where it is assumed that heat is transported by electrons and phonons. Also, in polar semiconductors in which there is electron and hole transport as well as phonons. From a microscopic point of view, the ballistic-diffusive models have been obtained from the Boltzmann equation [21]. The idea is to divide the distribution function into two parts, namely, one that describes the phonons that are generated at the boundaries of the material and that travel through it without being scattered (ballistic), and another one that describes the phonons being scattered in the bulk (diffusive). Chen showed that, up to first order in the expansion of the diffusive distribution function in terms of spherical harmonics, the associated heat flux can be described by an equation of the Cattaneo type, while the part of the heat flux transported by the ballistic phonons is obtained as a solution of the Boltzmann equation. The Cattaneo equation for the diffusive flow is coupled with the ballistic heat flow, so that the solution procedure involves first solving the Boltzmann equation to obtain the ballistic flux and then introducing it into the Cattaneo equation to obtain the diffusive heat flux.
In thermodynamic theories of heat transport, ballistic-diffusive models are obtained assuming that the heat flow can be separated into the ballistic and diffusive components. It has been then assumed that the diffusive part is described by the Cattaneo equation and the ballistic part by Guyer-Krumhansl equation [22,23]. We made the basic assumption that there exist two types of heat carriers and that the heat flux is separable, obtaining the constituent equations for the ballistic and diffusive heat flux from the second law of thermodynamics. This was done by including in the space of variables the ballistic and the diffusive heat fluxes. The system consists of two particle populations, namely, ballistic and diffusive heat carriers which have inner energies e b and e d , and heat fluxes q b and q d , respectively. In this way, the variable space becomes e b , e d , q b , q d . Then the balances of the component energies are: where use has been made of the caloric equations e d = C d T d and e b = C b T b , being T d and T b interpreted as quasi-temperatures and C d and C b specific heats of the species. We have assumed that the heat exchange between the ballistic and diffusive parts decomposes to a dissipative part Q and a non-dissipative one that is linear in the quasi-temperature. A is a constant quantity Note that Equations (1) and (2) are not counterbalanced which implies that the total internal energy balance equation has a source term. One can see that Equations (1) and (2) are more general than the corresponding Chen's internal energy balance equations [21]. We can obtain Chen's equations from Equations (1) and (2) by assuming that Q = 0 and A = . In one dimension, we obtain while Chen's equations are obtained approximately if T d ≈ T b . These balances are the most important assumption of Chen's theory, because the internal energy exchanges are considered as external sources in the balances, that is, there is not associated entropy production (see, e.g., [24]). Finally, observe that in Chen's approximation the total internal energy is conserved. Following the usual procedure [25,26], the second law of thermodynamics leads to the constitutive equations for the ballistic and diffusive heat fluxes. In the above expression s is the entropy density and J s is the entropy density flux. In order to get the constitutive equations we express s and J s in the following way and with m 1 and m 2 positive constants and b d and b b second order tensors. These tensors are called current multipliers and they are interpreted as constitutive functions. By substituting Equations (4) and (5) in the second law, inequality Equation (3), and introducing the internal energy balances as constrains we arrive to where I is the identity tensor. One can realize that the second part of the heat exchange, which is linear in the temperature does not contribute to the entropy production. This is an exact realization of the heat supply see, e.g., Wang and Hutter [24]. If the material is one dimensional and isotropic, the linear solution of Equation (6) is where ∂s eq ∂e b and ∂s eq ∂e d have been interpreted as the inverse of the quasi-temperatures T b and T d respectively, and the following conditions must be satisfied by the coefficients l 1 , l 2 , k 1 , k 2 , k 12 , k 21 : that is, the symmetric part of the conductivity matrix must be positive semidefinite, in order to fulfill the inequality of the second law, Equation (6). In a three-dimensional treatment, one needs a complete isotropic representation of the corresponding tensors. The above equations are simplified also by neglecting the cross coefficients in the last two equations, because these are not necessary for the reconstruction of the ballistic-diffusive theory. By eliminating the current multipliers from Equation (7), constitutive equations for the ballistic and diffusive heat fluxes are obtained: respectively. In these equations τ d and τ b are relaxation times of the heat fluxes, λ d , λ b are the heat conducting coefficients of the ballistic and diffusive heat carriers and κ 12 and κ 21 are the related Guyer-Krumhansl coefficients, considering direct damping of heat waves. In case of one temperature systems these are well motivated by kinetic theory calculations. l d and l b are interpreted here as the mean attenuation lengths of the heat carriers. d and b are coupling parameters. They are given in terms of the coefficients l 1 , l 2 , k 1 , k 2 , k 12 , k 21 as follows: . The attenuation length interpretation is supported by dimensional considerations and also by the thermodynamic inequalities.
It should be noted that Equations (8) and (9) are coupled by cross-effects (terms in d , b , κ 12 and κ 21 ).
Then assuming that the cross effects associated to the Laplacian of the heat fluxes are negligible, that is κ 12 = κ 21 = 0, the combination of the constitutive equations with those of internal energy balance (with Λ = 0) is straightforward to eliminate the heat fluxes. Then, two coupled transport equations are obtained which, in dimensionless form, read Entropy 2020, 22, 167 6 of 15 where the timescale was fixed by τ b and the length scale by an external length L. Equations (11) and (12) are the main result of this subsection. They are parabolic type equations due to the presence of terms in third order derivatives. The non-dimensional coefficients are given by the expressions It is interesting to note that these phenomenological coefficients contain information, on the one hand, on the thermal properties of materials and on the other, on properties of the interaction of the heat carriers with the lattice (α and diffusive and ballistic Knudsen numbers K nd and K nb , respectively).
We close this subsection by mentioning that if one assumes that A = , as before to obtain Chen's internal energy balances and b = d = 0 (other couplings neglected), the following diffusive-ballistic transport equations are obtained, These equations may be compared with the diffusive-ballistic model by Lebon and Machrafi [23]. The second one coincides with Lebon's ballistic transport up to constants and the first one coincides with Lebon's diffusive transport equation up to constants and a coupling term which has been neglected here. Lebon and Machrafi's model will be compared with that obtained in this paper in a further subsection.

Effective Medium Theory
It is worth mentioning that, in Knudsen's numbers (Equation (13)), information about the microstructure of the superlattice may be introduced through the mean attenuation lengths of the heat carriers, which depends on the period thickness. Clearly, the coefficient α contains information about the microstructure since the relaxation times of the species also depend on the period thickness. The way in which the mean attenuation lengths and relaxation times depend on the number of periods can be obtained from a microscopic theory or from experimental measurements. This issue will be discussed in following sections.
Let us start by obtaining the dispersion relationships of the ballistic and diffusive heat carriers. This step implies a precise definition of the physical nature of the heat carrier. Let us consider that the energy is transported by traveling thermal waves of small amplitude δT 0 , frequency ω and wave number k Substituting this expression in the transport Equations (11) and (12), the dispersion relationships are obtained. These rather long expressions will not be shown at this stage. Let us define now the effective thermal conductivity as This definition is motivated by the dispersion relation of the Fourier theory. It defines the "modes" of heat conduction, like the propagation speeds defines the modes of wave propagation. We assume that expression (17) describes the properties of heat transport in the superlattice. The effective thermal conductivities for the ballistic and diffusive carriers are obtained by using the dispersion relations leaving the result respectively.
In effective medium theories, thermal conductivities may depend not only on the frequency of heat carriers but their wavelength as well. This is suggested several times before us (see Section 4.2). The physical reason is that characteristic length scales related to material heterogeneities or other size related effects are not negligible and influence the heat transfer. With an extended continuum theory, the complete system of equations is considered in calculating the effective properties, and therefore it is possible to calculate explicitly the effective heat conductivity.
If the two heat transport channels (ballistic and diffusive) are considered to act in parallel, the total effective thermal conductivity of the superlattice is given by: This is the main result of this work. Expression (20) with (18) and (19) will be compared with some other formulations. In particular, in Section 4.3, it will be compared with that obtained from Lebon-Machrafi's model [23] and in Section 4.4 with the high order dissipative fluxes formalism [27]. In order to make contact with Lebon-Machrafi's theory we take A = to introduce Chen's balances in Equations (18) and (19). In such cases, the ballistic and diffusive effective conductivities reduce to the following expressions (we also take b = d = 0) : Expression (20) with (21) and (22) will be compared with the mentioned two formalisms in the following sections.

Heat Carrier Properties
We start this subsection by mentioning that the parameters appearing in Equations (18) and (19) or (20) and (21) of the diffusive and ballistic heat conductivities depend on the period thickness. We now show how expressions for those parameters may be obtained. First, the Knudsen numbers are expressed as where we have taken the superlattice period thickness p as the characteristic length of the system. The diffusive attenuation length l d is considered to depend also on the wave number as a dynamic transition to the coherent transport regime in the diffusive part of the heat flux is involved (this point is discussed below). In principle, see Equation (13), the definition of α is τ d /τ b , the ratio of the mean diffusive and ballistic attenuation times which are given by the phase velocity of diffusive and ballistic heat carriers respectively. In this way, α = l d v b p /l b v d p . From the dispersion relations, derived from Equations (11) and (12), the phase velocities can be obtained as well as the mean attenuation lengths. They also depend on α, since α appears in the dispersion relations. As a consequence, finding α implies solving a transcendent equation. We have made this for a few values of the period thickness and obtained Figure 1.
where we have taken the superlattice period thickness p as the characteristic length of the system.
The diffusive attenuation length d l is considered to depend also on the wave number as a dynamic transition to the coherent transport regime in the diffusive part of the heat flux is involved (this point is discussed below). In principle, see Equation (13), the definition of α is b d τ τ / , the ratio of the mean diffusive and ballistic attenuation times which are given by  (11) and (12), the phase velocities can be obtained as well as the mean attenuation lengths. They also depend on α , since α appears in the dispersion relations. As a consequence, finding α implies solving a transcendent equation. We have made this for a few values of the period thickness and obtained Figure 1. As can be seen, α is a decaying function of the period thickness. We introduce this property in the calculation of the effective heat conductivity by assuming that  (18) and (19) or (21) and (22). They are used to calculate the total effective heat conductivity of the superlattice accordingly with Equation (20). We do this in the following subsection.

Minimum Effective Heat Conductivity
In this subsection, Equation (20) with Equations (21) and (22) is fitted to the experimental values measured by Saha and collaborators [2,3]. To achieve this, we have four fitting parameters, namely, As can be seen, α is a decaying function of the period thickness. We introduce this property in the calculation of the effective heat conductivity by assuming that α = c 1 /p + c 2 (24) with fitting constants c 1 , c 2 . Mention must be made that the constant c 2 is responsible for the saturation shown by Saha's experimental result. Additional comments on this point are given in the discussion section.
The effective heat conductivities, as functions of the period thickness and wave number, are obtained by substituting the above expressions in Equations (18) and (19) or (21) and (22). They are used to calculate the total effective heat conductivity of the superlattice accordingly with Equation (20). We do this in the following subsection.

Minimum Effective Heat Conductivity
In this subsection, Equation (20) with Equations (21) and (22) is fitted to the experimental values measured by Saha and collaborators [2,3]. To achieve this, we have four fitting parameters, namely, the mean attenuation lengths of the diffusive and ballistic heat carriers and the constants c 1 and c 2 . Figure 1 shows the result when l d = 10 nm, l b = 8 nm, c 1 = 13 nm, c 2 = 0.4.
As may be observed, the effective thermal conductivity given by Equation (20) shows a minimum in the effective heat conductivity around p cr = 4 nm and the corresponding qualitative behaviour for p < p cr and p > p cr , including the saturation for higher values of the period thickness. Figure 2 is the main result of this work.   Figure 2 is the main result of this work.

Preliminar
To summarize, we remark two facts. First, we have obtained a heat transport model to describe effective properties of a superlattice from non-equilibrium thermodynamics with two internal variables, namely, the ballistic and the diffusive heat fluxes. Second, the derived Equation (20), with (21) and (22), reproduces the minimum and the saturation effect of the effective thermal conductivity in the superlattice. Our result shows a reasonable qualitative fitting with the experimental data.

On the Transition to Coherent Heat Transport Regime for Small Period Thickness.
It has been argued in the literature that the minimum of the effective heat conductivity when the period thickness is around 4 nm is due to the transition from coherent to incoherent propagation of the heat carriers. In this short subsection we discuss this point in terms of the power spectrum of the temperature fluctuations as predicted by the diffusive and ballistic transport Equations (14) and (15). We begin by mentioning that fluctuations are introduced through an additive stochastic term in the constitutive equations (Equations (8) and (9)). In the context of linear irreversible thermodynamics, the statistical properties of the stochastic term are determined by the fluctuation-dissipation assuming a Gaussian stochastic term [28,29]. In the case of the diffusive part of the heat flux, the power spectrum of the diffusive temperature shows the transition from diffusive transport to wave propagation of the heat carriers. In Figure 3, the spectrum has been plotted for a Knudsen diffusive number of 1.5 and different values of the wave number, and shows that for the considered period thickness there is a "window" of wave number values where heat is being  (21) and (22). The model predicts the experimentally found minimum around a period thickness of 4 nm as well as the saturation for higher values of the thickness. Black solid line: theoretical values; blue dots: experimental values from Saha et al. [2]. k = 1.2/ nm.

Preliminar
To summarize, we remark two facts. First, we have obtained a heat transport model to describe effective properties of a superlattice from non-equilibrium thermodynamics with two internal variables, namely, the ballistic and the diffusive heat fluxes. Second, the derived Equation (20), with (21) and (22), reproduces the minimum and the saturation effect of the effective thermal conductivity in the superlattice. Our result shows a reasonable qualitative fitting with the experimental data.

On the Transition to Coherent Heat Transport Regime for Small Period Thickness
It has been argued in the literature that the minimum of the effective heat conductivity when the period thickness is around 4 nm is due to the transition from coherent to incoherent propagation of the heat carriers. In this short subsection we discuss this point in terms of the power spectrum of the temperature fluctuations as predicted by the diffusive and ballistic transport Equations (14) and (15). We begin by mentioning that fluctuations are introduced through an additive stochastic term in the constitutive equations (Equations (8) and (9)). In the context of linear irreversible thermodynamics, the statistical properties of the stochastic term are determined by the fluctuation-dissipation assuming a Gaussian stochastic term [28,29]. In the case of the diffusive part of the heat flux, the power spectrum of the diffusive temperature shows the transition from diffusive transport to wave propagation of the heat carriers. In Figure 3, the spectrum has been plotted for a Knudsen diffusive number of 1.5 and different values of the wave number, and shows that for the considered period thickness there is a "window" of wave number values where heat is being transported wavely. In the case of Figure 3, the window goes from k = 0.39/nm (green line) to k =1.2/nm (brown line) with a period 13.91 nm. This is the result of a coherent transport of the heat carriers represented by the heat waves, Equation (16). The existence of two dominant frequencies for wave numbers within the window is an expression of the conservation of the phase in scattering processes of the heat carriers in a specific length scale which depends on the period thickness as it may be inferred from Figure 4. In this figure the amplitude of the wavely transport window is shown as a function of the period thickness. As the wavelength is a wave parameter, Figure 4 shows the length scale where the heat carrier phase is conserved. numbers within the window is an expression of the conservation of the phase in scattering processes of the heat carriers in a specific length scale which depends on the period thickness as it may be inferred from Figure 4. In this figure the amplitude of the wavely transport window is shown as a function of the period thickness. As the wavelength is a wave parameter, Figure 4 shows the length scale where the heat carrier phase is conserved. It must be noted that the window practically does not exist for higher periods as Figure 4 shows, while for the smaller periods it is maximum around the minimum of the heat conductivity (not shown in Figure 4).   It must be noted that the window practically does not exist for higher periods as Figure 4 shows, while for the smaller periods it is maximum around the minimum of the heat conductivity (not shown in Figure 4).  It must be noted that the window practically does not exist for higher periods as Figure 4 shows, while for the smaller periods it is maximum around the minimum of the heat conductivity (not shown in Figure 4).
Finally, the mean attenuation length for the diffusive part of the heat flux as a function of the period thickness for different frequencies can be seen in Figure 5. It was calculated from the imaginary part of the wave number as derived from the dispersion relation obtained from Equation (14). As expected, this length increases as the superlattice period thickness diminishes. For the smallest frequencies, there exists a certain period thickness such that for smaller periods the mean attenuation length is bigger than the corresponding value of the period. These facts have an effect on the effective thermal conductivity through the dependence of the non-dimensional parameter α on the period thickness. All of them reflect the physics, explaining the existence of the minimum in the effective heat conductivity. It is worth mentioning that, as implied by the properties of the power spectrum of the diffusive part of the heat flux, the system is not in a pure heat transport regime. For the higher values of the period thickness considered there exist heat carriers being transported diffusively and ballisticaly in the system, whereas for smaller values of the period there exist, in fact, three transport modes, namely ballistic, diffusive, and wave transport. the period thickness. All of them reflect the physics, explaining the existence of the minimum in the effective heat conductivity. It is worth mentioning that, as implied by the properties of the power spectrum of the diffusive part of the heat flux, the system is not in a pure heat transport regime. For the higher values of the period thickness considered there exist heat carriers being transported diffusively and ballisticaly in the system, whereas for smaller values of the period there exist, in fact, three transport modes, namely ballistic, diffusive, and wave transport.  (14). Equation (24) was used for the parameter α with the same constants 1 c and 2 c used for fitting the total effective thermal conductivity to the experimental data in Figure 2.
We close this subsection by mentioning that the ballistic part of the heat flux does not show a dynamic transition like the one found for the diffusive part. [22,23] In this subsection we display the effective thermal conductivity from ballistic-diffusive model by Lebon et al. [22,23], where two transport equations for the quasi-temperatures of ballistic and diffusive heat carriers were obtained under the assumption that the ballistic heat flux satisfies a Cattaneo constitutive equation and the diffusive heat flux a Guyer-Krumhansl equation. This differs from the scheme derived in this work wherein both of the fluxes satisfy Guyer-Krumhansl type equations. Our results show that this is crucial to understand the appearance of the minimum in the effective heat conductivity of the superlattice. Moreover, their heat transport equations are coupled by a term in the equation for the diffusive quasi-temperature which is linear in the ballistic quasi-temperature. The transport equation for the ballistic quasi-temperature does not contain coupling terms. In our formalism, an additional coupling term proportional to the time derivative of the ballistic quasi-temperature was obtained in the diffusive equation. It is worth mentioning that both terms come from a bit more general representation expressions of vector fields, Equation (7). In Figure 5. Diffusive mean attenuation length as a function of the period thickness for different frequencies (Blue: 100 Hz, Orange: 3 Hz, Green: 1.2 Hz and Red: 1 Hz). The mean attenuation length was calculated from the imaginary part of the wave number as derived from Equation (14). Equation (24) was used for the parameter α with the same constants c 1 and c 2 used for fitting the total effective thermal conductivity to the experimental data in Figure 2.

Comparison with Lebon, Machrafi and Grmela model
We close this subsection by mentioning that the ballistic part of the heat flux does not show a dynamic transition like the one found for the diffusive part. [22,23] In this subsection we display the effective thermal conductivity from ballistic-diffusive model by Lebon et al. [22,23], where two transport equations for the quasi-temperatures of ballistic and diffusive heat carriers were obtained under the assumption that the ballistic heat flux satisfies a Cattaneo constitutive equation and the diffusive heat flux a Guyer-Krumhansl equation. This differs from the scheme derived in this work wherein both of the fluxes satisfy Guyer-Krumhansl type equations. Our results show that this is crucial to understand the appearance of the minimum in the effective heat conductivity of the superlattice. Moreover, their heat transport equations are coupled by a term in the equation for the diffusive quasi-temperature which is linear in the ballistic quasi-temperature. The transport equation for the ballistic quasi-temperature does not contain coupling terms. In our formalism, an additional coupling term proportional to the time derivative of the ballistic quasi-temperature was obtained in the diffusive equation. It is worth mentioning that both terms come from a bit more general representation expressions of vector fields, Equation (7). In obtaining the heat transport equations, we have neglected the source term in the balance Equations (1) and (2). The contribution of this term to the coupling problem deserves to be investigated. Our internal energy balance Equations (1) and (2) contain Chen's theory as an approximation.

Comparison with Lebon, Machrafi and Grmela model
We finish this part by displaying the effective heat conductivity as predicted by Lebon-Machrafi model. We use the same scenario for the movement of the heat carriers described in Section 3.1. The result can be seen in Figure 6. It is observed that, not only does this model not reproduce the minimum of the effective thermal conductivity, but the fitting to the experimental values (blue dots) is not so good either.
(1) and (2). The contribution of this term to the coupling problem deserves to be investigated. Our internal energy balance Equations (1) and (2) contain Chen's theory as an approximation.
We finish this part by displaying the effective heat conductivity as predicted by Lebon-Machrafi model. We use the same scenario for the movement of the heat carriers described in Section 3.1. The result can be seen in Figure 6. It is observed that, not only does this model not reproduce the minimum of the effective thermal conductivity, but the fitting to the experimental values (blue dots) is not so good either.  [27].

Comparison with High Order Dissipative Fluxes Formalism by Álvarez and Jou
In this subsection, we compare the result of the previous section concerning the total effective heat conductivity to those predicted by the high order dissipative fluxes version of extended irreversible thermodynamics for a rigid solid by Álvarez and Jou [27]. In this formalism, the generalized thermal conductivity is obtained as ...
is obtained from the continued fraction expression (24). This spectral expression establishes the contribution of propagating modes of wave number k and frequency ω to the total thermal conductivity.  [27] In this subsection, we compare the result of the previous section concerning the total effective heat conductivity to those predicted by the high order dissipative fluxes version of extended irreversible thermodynamics for a rigid solid by Álvarez and Jou [27]. In this formalism, the generalized thermal conductivity is obtained as

Comparison with High Order Dissipative Fluxes Formalism by Álvarez and Jou
being K b (T) the bulk thermal conductivity. The coefficients τ i and l i have been interpreted as characteristic relaxation times and lengths of heat carriers, respectively. If all the relaxation times and characteristic lengths are considered to be equal, the asymptotic expression is obtained from the continued fraction expression (24). This spectral expression establishes the contribution of propagating modes of wave number k and frequency ω to the total thermal conductivity. From Equation (26), the thermal conductivity can be obtained in terms of the size L of the system by taking k = 2π/L, which in the steady state takes the form: where l is the mean attenuation length of heat carriers in the material. Equation (27) should be understood as the contribution of just the propagating mode with wave number k = 2π/L to the heat conductivity of the system. Nevertheless, it has been useful to describe size effects in nanoscaled systems and reproduces satisfactorily experimental results [27]. Thus, the equation for the heat flux can be written as follows q(x, t) = −K(T, l)∇T This is a Fourier-like heat conduction law, but it differs from the classical Fourier law in a fundamental way. It includes (through the generalized heat conductivity Equation (25)) the dynamic effects of the complete hierarchy of high order dissipative fluxes. Equation (27) corresponds to the stationary state for particular values of the wave number. In this way, size effects have been incorporated to the effective heat conductivity. At this stage it may be assumed that the total heat flux decomposes into a ballistic and a diffusive part [27,30]. Now we proceed further by introducing the expression for the mean attenuation length of the diffusive heat carriers in Equation (27) used in Section 3.1 (in order to fit our model to Saha's experimental results) and take the characteristic length L as the period thickness.
In Figure 7 it can be seen a plot of the effective thermal conductivity vs. period thickness as predicted by Equation (27). In general, the agreement of the calculated with the experimental data is good but it does not show the minimum in the effective thermal conductivity.
heat conductivity of the system. Nevertheless, it has been useful to describe size effects in nanoscaled systems and reproduces satisfactorily experimental results [27]. Thus, the equation for the heat flux can be written as follows This is a Fourier-like heat conduction law, but it differs from the classical Fourier law in a fundamental way. It includes (through the generalized heat conductivity Equation (25)) the dynamic effects of the complete hierarchy of high order dissipative fluxes. Equation (27) corresponds to the stationary state for particular values of the wave number. In this way, size effects have been incorporated to the effective heat conductivity. At this stage it may be assumed that the total heat flux decomposes into a ballistic and a diffusive part [27,30]. Now we proceed further by introducing the expression for the mean attenuation length of the diffusive heat carriers in Equation (27) used in Section 3.1 (in order to fit our model to Saha`s experimental results) and take the characteristic length L as the period thickness.
In Figure 7 it can be seen a plot of the effective thermal conductivity vs. period thickness as predicted by Equation (27). In general, the agreement of the calculated with the experimental data is good but it does not show the minimum in the effective thermal conductivity.

Conclusions
Some conclusions derived from the above study are presented in order.
The ballistic-diffusive model described in this work for the heat transport in superlattices contains the essence of the wave properties of heat carriers which explain their transport properties.

Conclusions
Some conclusions derived from the above study are presented in order.
The ballistic-diffusive model described in this work for the heat transport in superlattices contains the essence of the wave properties of heat carriers which explain their transport properties. The wave-like behavior of the heat carriers is represented by the non-dimensional coefficients in Equations (11) and (12). The so obtained heat transport model reveals that the heat transport regime in the superlattice is composed of a part of heat carriers being diffusively transported and another part being ballistically transported through the system. The two regimes coexist and together constitute two heat transport channels with specific transport properties. The microstructure enters the formalism through the characteristic length of the superlattice, namely, the period thickness, the mean attenuation length and the phase velocity of heat carriers. Remarkably, the transition found from diffusive transport to wave propagation of heat carriers is in line with the physical explanation of the existence of the minimum of the effective heat conductivity of superlattices with small period thickness. This transition introduces a third heat transport regime which influence the dynamics through the parameter α explaining the minimum of the effective heat conductivity.
The developed ballistic-diffusive model requires that both components of the total heat flux satisfy Guyer-Krumhansl type equations which lead to coupled heat transport equations for the ballistic and diffusive quasi-temperatures. These equations are of the parabolic type and this seems to be another determining property to predict the existence of a minimum in the effective thermal conductivity. Accordingly, with the scenario depicted here in Sections 3.1 and 4.2 for the properties of the heat carriers, the minimum is due to the fact that the long wavelength wavely propagating modes have a mean attenuation length larger than the characteristic length of the superlattice and even larger than its total length. This theoretical framework may shed light to the dependence of the mean attenuation length of the heat carriers on the microstructure. Through the proposed procedure to fit the model to the experimental data it can be found the explicit expression for the ballistic and diffusive Knudsen numbers as a function of the mean attenuation lengths. It would be expected that the prediction of the minimum of the effective thermal conductivity remain valid and a better fitting to the experimental data.