Abstract
A non-local and non-linear thermodynamical model of heat transfer at nanoscale beyond the well-known Maxwell–Cattaneo theory is derived. The compatibility of the proposed model with second law of thermodynamics is proved. The model is subsequently used to investigate the propagation of a heat pulse in one-dimensional nanosystems in the linear case. The predicted results are compared with those arising from the Maxwell–Cattaneo theory in order to point out the possible influence both of the non-local effects, and of the relaxation effects of the higher-order fluxes. Some problems related to initial data and boundary conditions are also discussed.
Similar content being viewed by others
1 Introduction
In the general context of continuum mechanics, a very important role is played by the so-called balance laws. From the theoretical point of view, a balance law is just the statement of the causes of the contents time-changes of a certain extensive physical quantity over a given domain. Those causes can be only of the following two different types:
-
(i)
Flux. The physical quantity varies owing to possible exchanges with the surroundings through the boundary of the domain.
-
(ii)
Production. The physical quantity varies since some (thermodynamic) process may occur inside the domain.
From the practical point of view, each balance law is expressed in terms of a partial differential equation (i.e., the so-called balance equation) which is usually written in a local form in such a way that it has to be valid in any point of a given domain. Among the different balance equations, in non-equilibrium thermodynamics a great importance deserves the (local) balance equation of the specific entropy s (per unit volume), namely,
with \(J^{\left( s\right) }_{i}\) being the specific-entropy flux, and \(\sigma ^{\left( s\right) }\) the specific-entropy production (per unit volume). In Eq. (1) the term \(J^{\left( s\right) }_{i,_{i}}\) (which accounts for the rate of local-entropy exchanged with the surroundings) may be zero, positive, or negative; according to the second law of thermodynamics, instead, in Eq. (1) the quantity \(\sigma ^{\left( s\right) }\) (deriving from processes occurring inside the system) has to be always non-negative (namely, it is zero at equilibrium, or for reversible transformations, whereas it is larger than zero for irreversible transformations) in any point and at any time, whatever the thermodynamic process is [1, 2]. Since \(\sigma ^{\left( s\right) }\) points out the sources of irreversibility occurring in any thermodynamic process, it should be (directly, or indirectly) related to all quantities which characterize the irreversible process at hand. In non-equilibrium thermodynamic theories, in fact, the physical constraint
plays a very important role: since all solutions of the balance equations have to comply with the condition in Eq. (2), in those theories that unilateral constraint is mainly used to determine the forms either of the constitutive relations [3], or of the evolution equations of the state-space variables [2]. In pursuing these goals, as it will be briefly shown below in the case of two very well-known theories of heat transfer in a rigid body (just for the sake of simplicity), besides the correct choice of the state-space variables, a crucial point is also the right setting of the relation between the specific-entropy flux \(J^{\left( s\right) }_{i}\) and the local heat flux \(q_i\) [4,5,6].
1.1 The Fourier theory
In the case of the Fourier (F-) theory of the heat transfer in a rigid body, the specific internal energy (per unit volume) u is the only state-space variable: in the F-theory, therefore, only the evolution equation of u is needed. The time variations of u, indeed, are stated by the balance law of energy, the local formulation of which reads
According to the point of view of the balance law, as a consequence, from Eq. (3) one may straightforwardly note that \(q_i\) is just the flux of u: the heat flux, therefore, is the only cause of time variations of the internal energy, in the absence of a source term (here and in what follows neglected just for the sake of simplicity). Since in this theory \(s=s\left( u\right) \), moreover, the use of Eq. (3) yields
once the usual relation
is employed,Footnote 1 with \(\theta \) being the non-equilibrium temperature [2, 7, 8]. If we recall that the classical constitutive assumption for the entropy flux states that the specific-entropy flux only depends on the heat flux and the non-equilibrium temperature [1], i.e.,
then Eq. (4) can be also rewritten as
A direct comparison between Eqs. (1) and (7) finally allows to recognize that
In order to avoid any possible violation of the unilateral constraint (2), in the F-theory it is therefore sufficient to assume
as the constitutive equation for the heat flux (which is worldwide known as the F-law), if it is tacitly understood that the positive-valued constant \(\lambda \) means the thermal conductivity. The coupling of Eqs. (8) and (9), in fact, lets that the form of the specific-entropy production reads
from which it is easy matter to see that \(\sigma ^{\left( s\right) }=\sigma ^{\left( s\right) }\left( u\right) \) is always a positive-definite function.
1.2 The Maxwell–Cattaneo theory
The Maxwell–Cattaneo (MC-) theory has been the first attempt to go beyond the validity’s limits of the F-theory. Those limits are particularly evident when the heat transfer is a very fast and/or steep phenomenon (e.g., in the case of ultrasound propagation, the light scattering in gases, shock waves, etc.) [2, 9, 10], or whenever the relaxation time (i.e., the temporal lagging between the application of a temperature gradient and the appearance of a heat flux) becomes large enough to require a description of the transient regime [11].
According to the point of view of Extended Irreversible Thermodynamics (EIT) [2, 7], in the MC-theory of heat transfer in a rigid body the state space is spanned by u and \(q_i\); from the physical point of view, this is tantamount to suppose that both the independent variations of u and \(q_i\) can, in principle, influence the consequent variations of the other thermodynamic fields, as for example the specific entropy.
Whereas u is still ruled by the balance law in Eq. (3), in this theory \(q_i\) can be no longer given by a constitutive relation (as it has been done in the case of the F-theory), but an evolution equation is needed for it. In order to theoretically derive it, we may observe that in the MC-theory \(s=s\left( u,q_i\right) \) in such a way that in the framework of EIT the generalized Gibbs equation read [2, 7]
with the (positive) material function \(\tau _R\) being the relaxation time of \(q_i\), namely, the response time for the onset of the heat flow after a temperature gradient is suddenly imposed. That quantity is substantially the reciprocal of the frequency of resistive phonon collisions, i.e., the mechanisms of phonon scattering in which the phonons’ energy is conserved, but not their momentum.Footnote 2 From Eq. (11) we have that the time derivative of s reads
if Eq. (3) is employed. The use of the constitutive assumption in Eq. (6) allows then to rewrite Eq. (12) as
A direct comparison of Eqs. (1) and (13) yields
as the constitutive equation for the specific-entropy production. In order to guarantee the (thermodynamic) compatibility of the MC-theory with the unilateral constraint (2), it is still sufficient that \(\sigma ^{\left( s\right) }=\sigma ^{\left( s\right) }\left( u,q_i\right) \) is a positive-definite function. For this reason in Eq. (14) we may set
which is also the required evolution equation of \(q_i\) and it is worldwide known as the MC-equation [12]. When Eq. (15) is introduced into Eq. (14), in fact, the following local form of the specific-entropy production for the heat transfer in a rigid body is recovered
in the case of the MC-theory. It seems worth noticing that Eqs. (15) and (16) reduce, respectively, to Eqs. (9) and (10) whenever \(\tau _R\rightarrow 0\). From the above considerations, it arises that the MC-theory is compatible with the classical constitutive equation for \(J^{\left( s\right) }_{i}\) [1], too.
1.3 Aims and scopes
The increasing interest in nanoscale heat-transport sciences led modern non-equilibrium thermodynamic theories to depart both from the F-theory, and from the MC-theory. Dealing with the heat transfer at nanoscale is not properly an easy task; from the theoretical point of view, to tackle with it several interesting approaches can be found in the literature, as it can be seen, for example, in the outstanding review by Guo and Wang [13]. All those approaches lie either on macroscopic, or mesoscopic, or microscopic methods. Among those methods, one spreading over a large span of time is the phonon hydrodynamics model, the usual starting point of which is the so-called Guyer–Krumhansl equation [2, 14], i.e.,
wherein \(\ell \) is the mean-free path of phonons (i.e., quanta of the vibrational mechanical energy arising from oscillating atoms within a lattice) representing the heat carriers in dielectric solids. Equation (17) represents an interesting step beyond Eq. (15) (the latter recovered from the former whenever \(\ell \rightarrow 0\), for example) since it is able to capture non-local effects in the heat transport, a topic of much interest at nanoscale wherein even small temperature differences may lead to high temperature gradients.
In nanosystems, indeed, also non-linear effects may strongly influence the electronic and optical properties: it could be, therefore, important to deeply examine them by introducing generalized nonlinear and non-local heat-transport equations [15]. Whereas the possibility of dealing with a non-local and non-linear model is clearly an important step toward a detailed and more refined description of the heat transfer at nanoscale, it has to be however noted that the introduction of non-linear terms may especially yield additional difficulties among which here we underline the correct setting-up of initial data (ID) and boundary conditions (BCs), which is of itself already a difficult task in the linear case (i.e., when only non-local effects are taken into account). In any model, in fact, both ID, and BCs should be always assigned in such a way that they have to be physically meaningful, besides correctly depicting the experimental reality, and the consequent initial-boundary value problem (IBVP) is mathematically consistent.
Along with the aforementioned observations, in the present paper:
-
In Sect. 2, we derive a non-linear and non-local model of heat transfer in rigid nanosystems on macroscopic grounds and in agreement with the second law of thermodynamics. That model encompasses Eq. (17) as a special case. The aim of this section is to show how non-local and non-linear effects may naturally arise from Eq. (1).
-
In Sect. 3, the above model, in its linear formulation, is employed to analytically solve an IBVP which sketches the propagation of second sound in one-dimensional (1D) nanosystems perturbed by a thermal shock. Therein a comparison with the analytical results arising from the MC-theory is also made. The aim of this section is twofold:
-
1.
to point out the possible influence of the relaxation times of higher-order fluxes, as well as the role played by non-local effects on second-sound propagation;
-
2.
to point out a plausible mathematical way of assigning ID and BCs which are both physically consistent, and depict the exact experimental situation under consideration.
-
1.
-
In Sect. 4 final comments and future perspectives are given.
2 Non-local and non-linear heat conduction in rigid solids
Whereas in a theoretical model (dealing with the heat-transfer phenomenon at nanoscale) the non-local effects may be introduced by the spatial derivatives of the state-space variables, the non-linear effects may be instead substantially understood in the following two different ways:
-
(i)
by letting the several material functions depend on the state-space variables [16, 17];
-
(ii)
by letting products between some state-space variables (or their spatial derivatives) directly enter the model equations [8, 18].
Here we especially aim at analyzing the consequence of the second type of nonlinear effects; thus, throughout this paper we assume constant values for all material functions.
Besides the specific internal energy u and the heat flux \(q_i\), in our approach we also assume as independent (state-space) variable the flux of the heat flux \(Q_{ij}\). Referring interested readers to Refs. [2, 7, 13, 19, 20] for exhaustive theoretical considerations about the physical interpretation of the second-order tensor \(Q_{ij}\), here we only observe that within the kinetic theory of diluted gases one has
wherein \(f\left( \mathbf {r},\mathbf {C},t\right) \) is the distribution function (\(\mathbf {r}\) being the position of the generic particle and m its mass), and \(\mathbf {C}\) is the relative speed of the particles with respect to the barycentric velocity, with \(\Vert \mathbf {C}\Vert \) being its modulus.Footnote 3 Since \(Q_{ij}\) is a symmetric tensor we may split it into its deviatoric (traceless) part \({\mathop {Q_{ij}}\limits ^{\text {o}}}\) and its bulk part \(Q\delta _{ij}\), namely, in the present theory we postulate that the state space is
The basic idea underlying the above selection is that those fluxes, which are genuinely non-equilibrium quantities, constitute a rather natural way of accessing to a system through its boundary.
In agreement with both EIT [2, 7], and the Non-Equilibrium with Internal Variables [21], each of the above state-space variables has to be characterized by its own evolution equation. Recalling that the time rate of u is given by Eq. (3), in this section we therefore derive the evolution equations for the above (state-space) flux variables. To do this, we start to observe that the very general form of the constitutive equation of the specific entropy is
owing to Eq. (18). For this reason, in agreement with EIT (see Part III—Subsec. 9.3.2 in Ref. [2]), in our approach we generalize the Gibbs equation (11) as
wherein the two (positive-valued) material functions \(\tau _1\) and \(\tau _0\) mean the relaxation times of \({\mathop {Q_{ij}}\limits ^{\text {o}}}\) and Q, respectively, and account for the relaxation effects of the higher-order fluxes [2, 7]. In some situations, as for example in diluted gases, the relaxation times of higher-order fluxes are of the same order as that of the heat flux itself, in such a way that when \(q_i\) is taken as an independent variable, all higher-order fluxes should also be taken into account. Further considerations about \(\tau _1\) and \(\tau _0\) are postponed to Sect. 2.2.
At the present stage, we leave undetermined the other two quantities \(\eta _0\) and \(\eta _1\) entering Eq. (19). They however will be clearly identified in Sect. 2.1. For our scope it is enough to claim that they are positive quantities.
The properties of \({\mathop {Q_{ij}}\limits ^{\text {o}}}\) and Q allow to suppose that the ratios \(\frac{\tau _1}{\eta _1}\) and \(\frac{\tau _0}{\eta _0}\) remain constant [2], as it will be shown in Sect. 2.2. As a consequence of this from Eq. (19) we may formally obtain the following time derivative of s:
If we go beyond the classical constitutive assumption in Eq. (6) and therefore further suppose
according to EIT [2], by direct calculations Eq. (20) can be also rewritten as
By comparing Eqs. (1) and (22) we may then point out the following constitutive equation for the specific-entropy production:
Note that Eq. (23) encompasses Eq. (14) as a special case: in fact, the latter form of \(\sigma ^{\left( s\right) }\) is obtained by the former form of \(\sigma ^{\left( s\right) }\) whenever \({\mathop {Q_{ij}}\limits ^{\text {o}}}\rightarrow 0\) and \(Q\rightarrow 0\).
In order to guarantee that the unilateral constraint (2) is always fulfilled, it is sufficient to make assumptions letting the function \(\sigma ^{\left( s\right) }\) be positive definite. Since Eq. (23) substantially displays the bilinear form
wherein A, B and C are strictly negative coefficients, whereas the following quantities
are the so-called thermodynamic forces [1, 2], a possible way of fulfilling Eq. (2) in any situation is to assume
for example. The assumptions in Eqs. (26) straightforwardly lead to
2.1 The model equations
Provided the physical meanings of \(\eta _0\) and \(\eta _1\) are given, Eqs. (27) are proper the evolution equations of the (state-space) flux variables we’re searching for. In order to make those identifications, indeed, we may start to observe that whenever both \(\partial _t{\mathop {Q_{ij}}\limits ^{\text {o}}}\rightarrow 0\), and \(\partial _tQ\rightarrow 0\), the insertion of Eqs. (27b) and (27c) into Eq. (27a) turns out
wherein the second-order tensor
can be substantially meant as an effective thermal conductivity (which could be related to the compelling problem of flux limiters [22,23,24]). Whenever the non-equilibrium-temperature approximation [2, 16] holds, namely, if one has
and
with \(\theta _{\text {eq}}\) being the (constant value of the) local-equilibrium temperature, instead of Eqs. (19) and (21), respectively, the non-linear terms proportional to \({\mathop {Q_{ij}}\limits ^{\text {o}}}\theta _{,_{j}}\) and \(Q\theta _{,_{i}}\) vanish and the effective thermal conductivity above reduces to the bulk thermal conductivity. In that case, moreover, it is easy matter to observe that Eqs. (17) and (28) coincide provided the following identifications hold:
Equations (30) turn out the expressions of \(\eta _1\) and \(\eta _0\) in terms of clear physical quantities. If the assumptions in Eqs. (30) are made, from the local balance of energy (3) and Eqs. (27) we finally have
which are our model equations.
2.2 Comparison with the 9-moment theory
When the non-equilibrium-temperature approximation holds Eqs. (31) become
once the classical thermodynamic relation
is used, with \(c_v\) being the specific heat at constant volume per unit volume (always larger than zero, as experimental evidences point out).
We may indeed recognize that the model in Eqs. (32) is closely related to the model which characterizes the 9-moment theory [25], i.e., to the model based upon
wherein \(\Vert \mathbf {c}\Vert \) means the modulus of the phonons’ Debye speed (defined as some mean value of longitudinal and transversal speed of sound), and
with \(\hbar \) being the Plank’s constant, \(\mathbf {k}\) being the wave vector (the modulus of which is \(\Vert \mathbf {k}\Vert \)), the function \(f\left( \mathbf {x},t,\mathbf {k}\right) \) means the phase density. The relaxation time \(\tau _{\star }\) appearing in Eqs. (34) is such that
where \(\tau _N\) means the reciprocal of the frequency of normal phonon collisions, i.e., the mechanisms of phonon scattering in which both the phonons’ energy, and their momentum is conserved. By comparing Eqs. (32) and (34), indeed, it is possible to set the following parallelisms
between the field variables involved in those models, the only difference being the form of the isotropic part of the two symmetric second-order tensors \(Q_{ij}\) and \(N_{ij}\), which in the latter case is only related to u, as it follows from Eq. (35b). The parallelism between Eqs. (32) and (34), besides giving a common root (very well defined in the kinetic theory) to our model and the 9-moment theory, allows us to give a deeper physical insight into the relaxation times \(\tau _1\) and \(\tau _0\): both of them, in fact, should be related to the relaxation time \(\tau _{\star }\) defined by Eq. (36) which is tantamount to claim that the relaxation times of higher-order fluxes in our model are related both to the resistive scatterings, and to the normal scatterings.
The above considerations justify the assumption of constant values for the ratios \(\frac{\tau _1}{\eta _1}\) and \(\frac{\tau _0}{\eta _0}\) that we have made in Sect. 2. If Eqs. (30) hold, in fact, the above parallelism allows us to state that both \(\frac{\tau _1}{\eta _1}\), and \(\frac{\tau _0}{\eta _0}\) are proportional to \(\Vert \mathbf {c}\Vert ^2\), which can be surely supposed to be constant.
3 Second-sound propagation
In the present section we investigate the second-sound propagation predicted by the theoretical model in Eqs. (31). To this end, for computational needs, we introduce the following dimensionless variables:
with L being the characteristic length of our system, and \(T_0\) a suitable reference temperature (e.g., the initial temperature of the system). For our goal, we also make use of Eq. (33) such that Eqs. (32a) can be rewritten in terms of the time variation of non-equilibrium temperature \(\theta \). In this way Eq. (31) become
wherein we set \(\tau _0=\alpha \tau _R\), \(\tau _1=\beta \tau _R\), and \({\text {Kn}}=\frac{\ell }{L}\) stands for the Knudsen number.Footnote 4 According to what previously said, therefore, in the present analysis the non-dimensional coefficients \(\alpha \) and \(\beta \) account for the relaxation effects of the higher-order fluxes (related both to the resistive scatterings, and to the normal ones), whereas \({\text {Kn}}\) introduces non-local effects.
In closing this part, we note that in deriving Eqs. (39) we estimated the resistive relaxation time of phonons as
with \(\mathrm{V_{0}}\) standing for the average phonon speed (measured at the reference temperature-value \(T_0\)), and we used the well-known Ziman limit [26] for the thermal conductivity, i.e.,
as it is usually assumed in the hydrodynamic theory of second-sound propagation in insulators. Although the assumptions in Eqs. (40) and (41) may lead to a quantitative inaccuracy, it is clear that the results obtained below are not critically dependent upon those assumptions.
3.1 The mathematical formulation
Here we put ourself in the 1D case for the sake of simplicity. From the physical point of view this means that we are only interested in the thermal transport within 1D (indefinitely long) nanostructures (e.g., nanowires, nanofibers, nanotubes, and nanorods) which represent one of the best-defined and controlled classes of nanoscale building blocks. We refer to Refs. [27, 28] and the references cited therein for some comprehensive reviews about the use of these material systems.
Dealing with Eqs. (39) is not properly an easy task, at least if one is searching for analytical solutions which can depict the experimental evidences. One of the main problems in dealing with Eqs. (39) lies, for example, on the question regarding ID and BCs and the way they have to be correctly assigned, as it has been previously observed in Sect. 1. For this reason, as a first attempt to a more general analysis, to reduce our computations to a simpler level, here we put ourself in the non-equilibrium-temperature approximation and therefore we neglect the terms proportional to \({\mathop {\mathsf {H}}\limits ^{\text {o}}}_{ij}\mathsf {T}_{,_{\mathsf {j}}}\) and \(\mathsf {H}\mathsf {T}_{,_{\mathsf {j}}}\) in Eq. (39b). By indicating with \(\mathsf {x}\) the sole (non-dimensional) Cartesian coordinate, as a consequence of the assumptions above, Eqs. (39) reduce to
The above system of coupled PDEs will be solved by prescribing the following ID
which, from the physical point of view, mean that our rigid conductor is initially at equilibrium and at the reference temperature \(T_0\). Moreover, we also assign the following BCs
wherein \(\mathcal {H}\) is the Heaviside unit step function. This is tantamount to mathematically depicting the practical situation in which at the time \(\mathsf {t}=0^+\) the domain boundary at \(\mathsf {x}=0\) is only affected by a step-jump in temperature. The neglect of non-linear terms has allowed us to not assign any BCs for the other unknown functions \(\mathsf {h}_{\mathsf {x}}\left( \mathsf {x},\mathsf {t}\right) \), \({\mathop {\mathsf {H}}\limits ^{\text {o}}}_{\mathsf {xx}}\left( \mathsf {x},\mathsf {t}\right) \) and \(\mathsf {H}\left( \mathsf {x},\mathsf {t}\right) \): we let, in this way, the model furnish the values of those functions also at the boundary. Otherwise, we would have had to assign those BCs, thus incurring the problem of the physical meaning of the prescribed values for the unknown functions at the boundary, a problem which is still open.Footnote 5
The IBVP in Eqs. (42)–(44), therefore, is able to describe the behavior of an infinite 1D nanosystem which is perturbed from its initial equilibrium state by a thermal pulse (i.e., a thermal shock) applied in its left-hand side end. From the practical point of view, in Eq. (44) \(\mathsf {T}^{\bullet }\in \mathbb {R}^+\) stands for a dimensionless constant which captures the energy amount of the thermal shock perturbing the system at hand. Further comments about the ID in Eqs. (43) and the BCs in Eq. (44) will be made in Sect. 3.3.
3.2 The form of the basic fields
The approximate analytical solution of the IBVP in Eqs. (42)–(44) is the following:
wherein we have set
just for the sake of a compact notation. Interested readers can refer to “Appendix” at the end of this paper for a sketch of the solution method used to obtain the form of the basic fields in Eqs. (45); therein considerations about the degree of approximation of Eqs. (45) are also made. Here, instead, we aim at drawing the attention on some consequences arising from Eqs. (45).
In Fig. 1, we plot the theoretical behavior of the non-dimensional temperature \(\mathsf {T}\) arising from Eq. (45a) when \(\alpha =\beta =1\), \(\mathsf {T}^{\bullet }=1\) and \({\text {Kn}}=1\) (just for the sake of computational needs). That figure clearly displays the propagating temperature wave (so pointing out the hyperbolic behavior of Sys. (39) in the 1D linear case): in fact, only the region before the front wave is perturbed (i.e., therein \(\mathsf {T}\ne 0\)) whereas the region ahead the front wave is still in equilibrium (i.e., therein \(\mathsf {T}=0\)). Owing to dissipative effects, from Fig. 1 it is also possible to see how the energy amount of the thermal wave progressively shrinks for increasing time.
In the case of Eq. (45a), it seems also interesting to observe that the dissipative effects are influenced both by the relaxation times of Q and \({\mathop {Q_{ij}}\limits ^{\text {o}}}\) (accounted by \(\alpha \) and \(\beta \), respectively), and by the phonon mean-free path (accounted by \({\text {Kn}}\)).
In particular, we may infer that the temperature-wave energy amount
-
decreases as steeper as \(\alpha \) and \(\beta \) reduce. This can be seen, for example, in Fig. 2 wherein the temperature \(\mathsf {T}\) behavior, as a function of the distance \(\mathsf {x}\) from the thermal-shock source, is investigated for different values of \(\alpha \) and \(\beta \).
-
decreases as steeper as \({\text {Kn}}\) reduces. This can be seen, for example, in Fig. 3 wherein the temperature \(\mathsf {T}\) behavior, as a function of the distance \(\mathsf {x}\) from the thermal-shock source, is investigated for different values of \({\text {Kn}}\).
3.3 Results arising from the entropy balance
From the mathematical point of view, it is clear that the form of the basic fields in Eqs. (45) strictly depends on the model in Eqs. (42), as well as on the ID in Eqs. (43) and on the BCs in Eq. (44). Here we avoid deepening too much into the mathematical details (Eqs. (45) are however an implicit proof of the existence of the solution of the IBVP in Eqs. (42)–(44)), and only look at the physical meaning of Eqs. (43) and (44).
From the practical point of view, we stress that both ID, and BCs should arise from experimental evidences and have their own physical clear meaning. The determination of the initial conditions is generally of no concern either in steady-state situations, or in transient problems, since it is usually admitted that the system at hand is initially in a thermodynamic equilibrium. This is the main reason why ID have been given by Eqs. (43). The right setting-up of BCs, instead, is in general more delicate and controversial. It has to be also noted that in the case of Eqs. (42) the above problem becomes still more acute since therein the four unknown fields \(\mathsf {T}\left( \mathsf {x},\mathsf {t}\right) \), \(\mathsf {h}_{\mathsf {x}}\left( \mathsf {x},\mathsf {t}\right) \), \({\mathop {\mathsf {H}}\limits ^{\text {o}}}_{\mathsf {xx}}\left( \mathsf {x},\mathsf {t}\right) \) and \(\mathsf {H}\left( \mathsf {x},\mathsf {t}\right) \) are coupled.
As a consequence, the following questions may naturally arise: is the IBVP in Eqs. (42)–(44) physically consistent? Since in Sect. 3.1 only the value of the unknown function \(\mathsf {T}\left( \mathsf {x},\mathsf {t}\right) \) at the boundary has been given, are the BCs in Eq. (44) admissible? Whereas the results contained in Sect. 2 make us sure about the physical validity of the theoretical model in Eqs. (42), to answer the above questions in the present section we analyze the consequences of the unilateral constraint (2). In particular, if we introduce the following further dimensionless quantities
from Eqs. (21) and (23), respectively, we have
By neglecting the non-linear terms proportional to \(\mathsf {H}\mathsf {T}_{,_{\mathsf {j}}}\) and \({\mathop {\mathsf {H}}\limits ^{\text {o}}}_{ij}\mathsf {T}_{,_{\mathsf {j}}}\) in Eq. (48b), in the 1D case we have
In Fig. 4, we plot the theoretical behavior of the non-dimensional specific-entropy production \(\mathsf {\Sigma }^{\left( s\right) }\) arising from Eq. (49b) when \(\alpha =\beta =1\), \(\mathsf {T}^{\bullet }=1\) and \({\text {Kn}}=1\) (just for the sake of computational needs). As it was expected, Fig. 4 shows that only the (out-of-equilibrium) region before the front wave is characterized by \(\mathsf {\Sigma }^{\left( s\right) }\) values different from zero (but strictly positive, according to second law of thermodynamics). In the region ahead the front wave, which is instead still at equilibrium, it obviously holds \(\mathsf {\Sigma }^{\left( s\right) }=0\). The behavior of \(\mathsf {\Sigma }^{\left( s\right) }\) is monotonically decreasing: this is logical since in time the temperature wave progressively loses its energy amount, as observed in Sect. 3.2.
The above considerations are enough to claim the validity of the physical setting of the phenomenon at hand, i.e., Eqs. (43) and (44).
3.4 Comparison with the MC-theory
The limits of F-theory, which are well known since decades, have led to several interesting generalizations of it. Those generalizations are obviously not unique: they generally differ in the underlying physical principles and in modeling capabilities. In Ref. [21], for example, two theoretical models are discussed and compared: one based upon Rational Extended Thermodynamics, and the other one based upon Non-Equilibrium Thermodynamics with Internal Variables. In Ref. [29], instead, the Grad-13 equations [30] have been embedded into the framework of GENERIC [31]. Since the comparison between different theories is surely one of the major stimuli both to enhance previous models, and to find new others, by means of the second-sound propagation in this part we compare the results of Sect. 3 with those arising from the MC-theory. The MC-theory is considered here only for the fact that it can be recovered (as special case) from any theory beyond the F-law. We postpone to future works a more complete and refined comparison between the different models proposed in the literature.
Since in the MC-theory the two basic fields \({\mathop {\mathsf {H}}\limits ^{\text {o}}}_{\mathsf {xx}}\left( \mathsf {x},\mathsf {t}\right) \) and \(\mathsf {H}\left( \mathsf {x},\mathsf {t}\right) \) are not taken into account (according to what previously said in Sect. 1.2), then in that theory the IBVP given by Eqs. (42)–(44) reduces to
in such a way that one would have the following analytical solutions
instead of Eqs. (45a) and (45b). A sketch of the solution method used to obtain the form of the basic fields in Eqs. (51) can be found in “Appendix,” too.
In Fig. 5, we compare the behaviors of the analytical solutions in Eqs. (45a) and (51a) when \(\mathsf {t}=0.5\), \(\alpha =\beta =1\), \(\mathsf {T}^{\bullet }=1\), and \({\text {Kn}}=1\) just for the sake of computation. It clearly points out that the temperature wave predicted by Eq. (45a) (i.e., the dashed red curve therein) is faster than the temperature wave predicted by Eq. (51a) (i.e., the solid blue curve therein): in the former case, in fact, the front wave (characterized by the jump discontinuity) is placed at a distance larger than that at which the front wave in the latter case occurs. This result was indeed expected since the second-sound speed predicted by Eqs. (31) is larger than that predicted by the MC-theory [2, 32], whatever the values of the parameter \(\alpha \), \(\beta \) and \({\text {Kn}}\) are.
In the MC-theory the dimensionless 1D specific-entropy flux and production, respectively, read
In Figs. 6 and 7, we compare, instead, the behaviors of the analytical solutions by Eqs. (49b) and (52b) when \(\mathsf {t}=0.8\), \(\alpha =\beta =1\), \(\mathsf {T}^{\bullet }=1\), and \({\text {Kn}}=1\) just for the sake of computation. From that figures, it follows that both the predicted specific-entropy productions tend to a plateau limit which is attained before the front of the moving temperature wave.
4 Final comments and future perspectives
The raise of nanotechnology has required (and, indeed, is still requiring) several efforts to better understand the thermal-transport properties of nanodevices, the performance and reliability of which could be much influenced by relaxation (i.e., memory), non-local and non-linear effects [2, 7, 33,34,35]. For this reason, in Eqs. (31) of this paper we proposed a non-linear and non-local thermodynamical model of heat transfer at nanoscale in a rigid body beyond the usual classical theories. That theoretical model, which encompasses both the F-theory, and the MC-theory as special cases, has been derived in framework of EIT [2, 7] in Sect. 2 and essentially lies on the concept of flux of the heat flux \(Q_{ij}\) (which has been indeed decoupled into its deviatoric and bulk parts), the physical root of which arises from the kinetic theory [2].
Since in the present paper we derived Eqs. (31) only by means of macroscopic theoretical considerations, further comments about those equations are therefore needed. Besides to the form of the specific-entropy flux in Eq. (21), the theoretical model in Eqs. (31) is indeed also strictly related to Eq. (23) which substantially displays the bilinear form (24) in fluxes and forces. Although experimental observations confirm that thermodynamic forces and fluxes are interwoven [1, 2], the exact relations among them are indeed not known a priori. In our approach we considered linear force-flux relations and, therefore, we finally set Eqs. (25). Non-linear flux-force relations could have been also considered in principle [36,37,38,39], indeed, in such a way that a different non-linear and non-local model of heat transfer was obtained. In the present paper, however, we did not purse the approach of non-linear flux-force relations since it is generally related to theoretical questions concerning the interpretation of the second law (see Part I—Sec. 2.4 in Ref. [2], for example) which are surely conceptually interesting, but not properly in the focus of this paper.
By means of the so-called non-equilibrium-temperature approximation (i.e., whenever the difference \(\theta ^{-1}-\theta _{\text {eq}}^{-1}\rightarrow 0\)), in Sect. 3 the non-linear terms in Eqs. (31) have been neglected and the propagation of thermal pulses in 1D nanosystems has been consequently analyzed. That analysis, although simplified, could be however useful in practical applications. For example, in current days the very high circuit densities yield modern chips run so hot that heat is becoming one of the major problems in the design of integrated circuits. As a consequence, one of the aspects which is much being studied at nanoscale is the isolation of very small areas which could be heated by an external source (either static as a hot spot due to a very miniaturized working device, or dynamic as a fast laser pulse): a continuous control of the temperature’s rise, in fact, could be a crucial point for the correct operation of modern devices. In Sect. 3 we have tried to exactly tackle with that problem in 1D nanosystems. The analysis performed therein points out that one should pay attention to the role played by relaxation and non-local effects: as it follows from Fig. 2 (for the relaxation effects) and Fig. 3 (for the non-local effects) they may strongly influence the dissipative effects (i.e., one of the main—natural—tool for heat isolation of different components on a chip), thus not ensuring the right setting-up of modern devices which currently remains a very compelling challenge. That challenge could be also complicated by non-linear effects which may either enhance, or reduce, the influence of memory and non-local effects [15].
Another step toward a more detailed investigation will be surely the analysis of a pulse propagation once in Eqs. (31) no terms are neglected, or when non-linear flux-force relations are considered. We note that if Eqs. (31) are considered in their complete version, the interesting problem of flux limiters [22,23,24] also comes out: in that case, in fact, the form of the unknown fields will have to however guarantee that the effective thermal conductivity in Eq. (29) is positive defined.
The analysis performed in Sect. 3 could be, indeed, also used to indirectly compute a first rough estimation of the relaxation times of the high-order fluxes: whereas different models have been developed in order to estimate \(\tau _R\), very scant information about \(\tau _1\) and \(\tau _0\) (i.e., information on \(\tau _N\) according with the observations in Sect. 2.2) can be found in the literature. Since the wave front depends on the ratios between those relaxation times and \(\tau _R\), as Fig. 2 clearly shows, by checking the wave front in practical applications one could also estimate whether \(\tau _1\) and \(\tau _0\) are larger (or smaller) than \(\tau _R\), in principle, for example.
At the very end, it seems worth noticing that, in general, it is not so simple to claim which is the best (or, alternatively said, the right) proposal among the different theories beyond the F-law. To this end, practical applications in which different theories clearly predict different results could be suitable benchmarks. This is, for example, just the case of the heat-wave propagation we analyzed in the present paper: in fact, the results contained in Sect. 3 can be used to check, between the MC-theory and the model given by Eqs. (31), which one better fits with experimental evidences.
Notes
Whenever the internal energy is no longer the only one state-space variable (as it will be through this paper), Eq. (5) properly reads \(\frac{\partial s}{\partial u}=\frac{1}{\theta }\).
To be more precise we note that in EIT, when the generalized Gibbs equation is introduced, the function \(\frac{\partial s}{\partial q_i}\) is firstly written in a very general form and then, by sequential comparisons, it is completely identified in terms of quantities which have their own clear physical meanings. For the rigorous derivation of Eq. (11) interested readers are therefore referred to Ref. [2] (see therein Part I—Subsec. 2.1.2). Therein comments about the degree of approximation of Eq. (11) can be also found.
It seems worth noticing that in the present paper we assume that \(Q_{ij}\) has the unit of \({\text {Wm}}^{-1}{\text {s}}^{-1}\) along with Ref. [19] (see therein Sec. 2) and Ref. [2] (see therein Part III—Subsec. 9.3.1), for example. In Ref. [7], instead, \(Q_{ij}\) has the unit of \({\text {Wm}}^{-1}\) (see therein Subsec. 7.1.3).
We note that in the case of nanoscale heat transfer one usually has \({\text {Kn}}\sim \mathcal {O}\left( 1\right) \) (i.e., \(\ell \ge L\)), whereas the case \({\text {Kn}}\rightarrow 0^+\) (i.e., \(\ell \ll L\)) corresponds to the case of macroscale heat transfer.
In Ref. [2] interesting considerations can be found in this regard (see therein Box 9.1 in page 213).
References
De Groot, S.R., Mazur, P.: Non-equilibrium Thermodynamics. North-Holland Publishing Company, Amsterdam (1962)
Jou, D., Casas-Vázquez, J., Lebon, G.: Extended Irreversible Thermodynamics, 4th edn. Springer, Berlin (2010)
Muschik, W., Ehrentraut, H.: An amendment to the Second Law. J. Non-Equilib. Thermodyn. 21, 175–192 (1996)
Jou, D., Lebon, G., Mongioví, M.S., Peruzza, R.A.: Entropy flux in non-equilibrium thermodynamics. Physica A 338, 445–457 (2004)
Sellitto, A., Cimmelli, V.A., Jou, D.: Entropy flux and anomalous axial heat transport at the nanoscale. Phys. Rev. B 87, 054302 (2013)
Chowdhury, S.C., Roy, D., Reddy, J.N.: Relating entropy flux with heat flux in two-temperature thermodynamic model for metal thermoviscoplasticity. J. Appl. Mech. 84, 021007 (2017)
Lebon, G., Jou, D., Casas-Vázquez, J.: Understanding Non-equilibrium Thermodynamics. Springer, Berlin (2008)
Cimmelli, V.A., Sellitto, A., Jou, D.: Nonequilibrium temperatures, heat waves, and nonlinear heat transport equations. Phys. Rev. B 81, 054301 (2010)
Uribe, F.J.: Shock waves: the Maxwell-Cattaneo case. Phys. Rev. E 93, 033110 (2016)
Hennessy, M.G., Calvo-Schwarzwälder, M., Myers, T.G.: Modelling ultra-fast nanoparticle melting with the Maxwell-Cattaneo equation. Appl. Math. Model. 69, 201–222 (2019)
Herrera, L., Pávon, D.: Hyperbolic theories of dissipation: Why and when do we need them? Physica A 307, 121–130 (2002)
Cattaneo, C.: Sulla conduzione del calore. Atti Sem. Mat. Fis. Univ. Modena 3, 83–101 (1948)
Guo, Y., Wang, M.: Phonon hydrodynamics and its applications in nanoscale heat transport. Phys. Rep. 595, 1–44 (2015)
Guyer, R.A., Krumhansl, J.A.: Solution of the linearized phonon Boltzmann equation. Phys. Rev. 148, 766–778 (1966)
Sellitto, A., Cimmelli, V.A., Jou, D.: Analysis of three nonlinear effects in a continuum approach to heat transport in nanosystems. Physica D 241, 1344–1350 (2012)
Di Domenico, M., Jou, D., Sellitto, A.: Nonlocal and nonlinear contributions to the thermal and elastic high-frequency wave propagations at nanoscale. Int. J. Heat Mass Transf. 156, 119888 (2020)
Di Domenico, M., Jou, D., Sellitto, A.: Heat-flux dependence of the speed of nonlinear heat waves: Analogies with the Kerr effect in nonlinear optics. Int. J. Therm. Sci. 161, 106719 (2021)
Dong, Y., Cao, B.-Y., Guo, Z.-Y.: Generalized heat conduction laws based on thermomass theory and phonon hydrodynamics. J. Appl. Phys. 110, 063504 (2011)
Lebon, G., Torrisi, M., Valenti, A.: A non-local thermodynamic analysis of second sound propagation in crystalline dielectrics. J. Phys. Condens. Matter 7, 1461–1474 (1997)
Carlomagno, I., Di Domenico, M., Sellitto, A.: High-order fluxes in heat transfer with phonons and electrons: Application to wave propagation. Proc. R. Soc. A 477, 20210392 (2021)
Kovács, R., Madjarević, D., Simić, S., Ván, P.: Non-equilibrium theories of rarefied gases: internal variables and extended thermodynamics. Contin. Mech. Thermodyn. 33, 307–325 (2021)
Levermore, C.D.: Relating Eddington factors to flux limiters. J. Quant. Spec. Rad. Transf. 31, 149–160 (1984)
Anile, A.M., Romano, V.: Covariant flux-limited diffusion theories. Astrophys. J. 386, 325–329 (1992)
Sellitto, A., Cimmelli, V.A.: Flux limiters in radial heat transport in silicon nanolayers. J. Heat Transf. Trans. ASME 136, 071301 (2014)
Dreyer, W., Struchtrup, H.: Heat pulse experiments revisited. Contin. Mech. Thermodyn. 5, 3–50 (1993)
Ziman, J.M.: Electrons and Phonons. Oxford University Press, Oxford (2001)
Chen, R., Hochbaum, A.I., Murphy, P., Moore, J., Yang, P., Majumdar, A.: Thermal conductance of thin silicon nanowires. Phys. Rev. Lett. 101, 105501 (2008)
Cahill, D.G., Braun, P.V., Chen, G., Clarke, D.R., Fan, S., Goodson, K.E., Keblinski, P., King, W.P., Mahan, G.D., Majumdar, A., Maris, H.J., Phillpot, S.R., Pop, E., Shi, L.: Nanoscale thermal transport II. 2003–2012. Appl. Phys. Rev. 1(1), 011305 (2014)
Struchtrup, H., Öttinger, H.C.: Thermodynamically admissible 13-moment equations. Phys. Fluid 34, 017105 (2022)
Grad, H.: On the kinetic theory of rarefied gases. Commun. Pure Appl. Math. 2, 331–407 (1949)
Grmela, M., Öttinger, H.: Dynamics and thermodynamics of complex fluids. I. Development of a general formalism. Phys. Rev. E 56(6), 6620–6632 (1997)
Sellitto, A., Di Domenico, M.: Nonlocal and nonlinear contributions to the thermal and elastic high-frequency wave propagations at nanoscale. Contin. Mech. Thermodyn. 31, 807–821 (2019)
Tzou, D.Y.: Macro to Micro-scale Heat Transfer. The Lagging Behaviour. Taylor and Francis, New York (1997)
Chen, G.: Nanoscale Energy Transport and Conversion–A Parallel Treatment of Electrons, Molecules, Phonons, and Photons. Oxford University Press, Oxford (2005)
Zhang, Z.M.: Nano/Microscale Heat Transfer. McGraw-Hill, New York (2007)
Hurley, J., Garrod, C.: Generalization of the Onsager reciprocity theorem. Phys. Rev. Lett. 48, 1575–1577 (1982)
Verhás, J.: An extension of the governing principle of dissipative processes to nonlinear constitutive equations. Ann. Phys. 495, 189–193 (1983)
Johannessen, E., Kjelstrup, S.: Nonlinear flux-force relations and equipartition theorems for the state of minimum entropy production. J. Non-Equilib. Thermodyn. 30, 129–136 (2005)
Salazar, D.S.P., Landi, G.T.: Nonlinear Onsager relations for Gaussian quantum maps. Phys. Rev. Res. 2, 033090 (2020)
Keiffer, R.S., Jordan, P.M., Christov, I.C.: Acoustic shock and acceleration waves in selected inhomogeneous fluids. Mec. Res. Commun. 93, 80–88 (2018)
Acknowledgements
The authors acknowledge Professor Francesco Oliveri from the University of Messina for fruitful discussions about this work which has been carried out under the auspices of the GNFM (Italian National Group of Mathematical Physics, INdAM).
Funding
Open access funding provided by Universitá degli Studi di Salerno within the CRUI-CARE Agreement.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Andreas Öchsner.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
At the end of this paper, just for the sake of completeness, we give some details about the technique used to solve the above IBVPs, the corresponding results of which are contained in Sect. 3.2.
Analytical solution
By applying the (temporal) Laplace transform \(\mathcal {L}\left[ \cdot \right] \) and recalling the usual relation
being \(\nu \) the (complex) Laplace transform parameter and with \(F_0\left( \mathsf {x}\right) =\displaystyle \underset{\mathsf {t}\rightarrow 0^+}{\lim }\,\,\,F\left( \mathsf {x},\mathsf {t}\right) \), from Eqs. (42) we have the following system of ODEs
once the ID arising from the Laplace transform of Eqs. (43) have been employed in such a way they are embedded in the procedure. In the Laplace domain, moreover, the following BCs, arising from Eq. (44), have to be appended to Eqs. (53)
Since Eqs. (53) substantially constitute a system of linear equations, by coupling Eq. (53a) with the spatial derivatives of Eqs. (53b)–(53d), we may obtain the following second-order ODE (in the unknown Laplace transform \(\mathcal {L}\left[ \mathsf {T}\right] \))
wherein we set
with
just for the sake of a compact notation. Since the characteristic equation of ODE (55) is
by means of BCs in Eq. (54) we are led to
which reduces to
once infinitesimals of order lower than \(\nu ^{-1}\) have been neglected in computing \(\Omega \) (i.e., when the square root of Eqs. (56) is evaluated), and Eq. (46) hold. Considerations about the degree of the above approximation can be found in the second section of this appendix.
From Eqs. (53) we may straightforwardly have
The inverse (temporal) Laplace transforms of Eqs. (60) and (61) furnish Eqs. (45).
In the MC-theory Eqs. (50) holds; therefore in the Laplace domain one has
which leads to
once the BCs in Eqs. (54) is employed. The inverse (temporal) Laplace transforms of Eqs. (63) give Eqs. (51).
Numerical solution
Although we refer to Eqs. (45) as the analytical solution of the IBVP in Eqs. (42)–(44), it seems worth noticing that they are not exact solutions since we neglected infinitesimals of order lower than \(\nu ^{-1}\) in deriving Eq. (60). The small-time approximation is surely logical in very fast phenomena (as the pulse propagation at hand), but a priori the orders of infinitesimals which may be neglected are not known.
To be sure that the approximation we used does not significantly affect Eqs. (45), we also perform a comparison with a numerical solution. In Fig. 8, we show how effective is the approximation used to obtain the analytical expression for the temperature field in Eq. (45a) (which is fundamentally the real benchmark in the present approach). In particular therein the solid (blue) curve refers to Eq. (45a), whereas the dotted (red) curve has been obtained by inverting numerically the exact Laplace transform of the temperature field
through the modified Tzou series [40]
for \(\mathsf {M}=500\), with \(\Omega \) given by Eq. (56).
The very good agreement between the two solutions in Fig. 8 (obtained also for data captions different from those used in that figure) ensures that the aforementioned approximation plays no relevant role in Eqs. (45).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Di Domenico, M., Sellitto, A. & Zampoli, V. Thermal pulse propagation beyond the Maxwell–Cattaneo theory: Application to one-dimensional nanosystems. Continuum Mech. Thermodyn. 34, 1455–1474 (2022). https://doi.org/10.1007/s00161-022-01134-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00161-022-01134-3