# Extended irreversible thermodynamics of heat transport: a brief introduction/Soojusulekande poordumatute protsesside laiendatud termodunaamika. Luhike sissejuhatus.

1. INTRODUCTION. PHYSICAL MOTIVATIONSCurrent frontiers in technology, e.g. in materials sciences and nanotechnology, require generalized transport equations beyond the classical theory. For instance, transport equations for miniaturized systems whose size is comparable to the mean free path have become an important topic because of the surge of nanotechnology. Analogously, the behaviour of systems submitted to high-frequency perturbations, comparable to the reciprocal of internal relaxation times, is necessary to optimize the operation of high-frequency devices. Equations for heat, mass, charge, and momentum transport have been actively explored in several situations: in miniatarized electronic devices, in nanotubes and nanowires, in theoretical models of energy transport in one-dimensional chains, in rarefied gases, etc. How far thermodynamic formalisms are helpful or necessary in this endeavour is an open question, because the mentioned situations clearly exceed the limits of validity of the classical local-equilibrium thermodynamics. Therefore, their study is a challenge for nonequilibrium thermodynamics to better understand its basic concepts, its limits of application, and its frontiers.

Among the several thermodynamic theories going beyond the local-equilibrium approach, Extended Irreversible Thermodynamics (EIT) provides generalized transport equations that incorporate memory and nonlocal effects, that reduce to the classical transport equations for low frequency and short mean-free paths, and that are able to incorporate the high-frequency behaviour of dissipative signals as, for instance, a finite speed of propagation of heat pulses, and general features of ballistic behaviour for long mean-free paths. Thus, they provide not only small corrections to the usual transport equations, but they tarn out to be useful even under extreme non-equilibrium situations where the classical equations completely fail.

In this paper, we introduce the basic concepts of EIT, namely, the idea of the fluxes as independent variables of the entropy and entropy flux, and the procedure to obtain their evolution equations, which provide generalized transport equations obeying the second law of thermodynamics. The reader interested in this topic and in its relation with other current frontiers of non-equilibrium thermodynamics may go to wider presentations such as Jou et al. (2001a, b) and Lebon et al. (2008).

Finite speed of propagation for dissipative signals was one of the starting stimuli for developing extended thermodynamics, but it appeared initially as an academic motivation because it was thought that the measurement of the relaxation times for the fluxes would be in general rather difficult. However, the development of light-scattering and neutron-scattering techniques showed not only that such relaxation times could be measured with accurateness, but that they were dependent on the wavevector in such a way that nonlocal effects should be taken into consideration.

Accordingly, we will incorporate memory and second-order non-local effects, and will show that this describes a rich phenomenology in heat transport, including thermal waves and phonon hydrodynamics. Though such transport equations had been proposed earlier on different bases, it was generally not realized that they were violating the classical formulation of the second law of thermodynamics, as they led in some cases to a negative entropy production. Extended Irreversible Thermodynamics was able to circumvent this problem by introducing a generalized entropy and entropy flux depending not only on the classical variables but also on the fluxes. Extended Irreversible Thermodynamics provides relevant models to explore the validity of thermodynamic ideas beyond their usual range of application.

However, the inclusion of memory and second-order non-local effects, interesting as it is, cannot cope with the ballistic regime, where most of the particles cross all the system without experiencing practically any collision, except with the walls. To deal with this situation, it is necessary to include an infinite number of higher-order fluxes. Continued-fraction expansions for the effective thermal conductivity are established, leading to the expected limiting behaviours in collisional and ballistic regimes, and providing an interpolation in the intermediate regimes. Systems whose dimensions are comparable to or smaller than the mean free path of the carriers are not necessarily small. For instance, multilayer materials and superlattices, constituted by the superposition of thin layers, may exhibit non-Fourier heat transport when the width of each layer is smaller than the mean free path of the carriers. Such a configuration is characterized by a combination of ballistic regime and a set of barriers separating the successive layers. In ultrafast solidification, a transition from diffusive to ballistic behaviour is also observed (Herlach et al. 2007).

The behaviour of small systems is strongly influenced by non-linear effects, because even a small difference in temperature over a small distance generates very high gradients. As a consequence, the transport equations must explicitly incorporate nonlinear effects. In the present analysis, the heat flux will always have a finite, limited value, which will be of the order of the energy density times the maximum speed of signal propagation. The non-linear transport equations studied here are characterized by saturation in the fluxes for high enough values of the temperature gradient.

The generalized non-equilibrium entropy introduced by EIT, dependent on the fluxes, yields some insights into how to formulate a thermodynamic theory beyond local-equilibrium, an aim shared with other approaches (Eu 1992, 1998; Maugin and Muschik 1994a, b; Verhas 1997; Maugin 1999; Luzzi et al. 2001, 2002; Ottinger 2005). This implies a generalization of equations of state containing non-equilibrium contributions, which are useful to understand some open problems concerning the definition and meaning of temperature.

In summary, EIT is not only a conceptually interesting theory, with contributions to such basic questions as the form of the entropy and temperature in nonequilibrium states and the formulation of the second law in far-from-equilibrium situations, but it is also able to describe sophisticated experiments at micro- and nanoscales. Nowadays, a new area is open for EIT: the area of new technological applications.

2. MEMORY EFFECTS AND EXTENDED ENTROPY

The central relations in transport theories are the transport laws, describing the rate of transport of the basic variables in terms of the inhomogeneities of the system, yielding, for instance, the expressions of the fluxes of heat, mass, momentum, or charge as a function of the gradients of temperature, concentration, velocity, and electrical potential. In the classical transport theory, the fluxes J and the conjugated thermodynamic forces X are related by linear transport laws of the form J(t) = L x X (t) with L being suitable phenomenological coefficients. For instance, one may identify J with the heat flux q and X with the temperature gradient X = [nabla][T.sup.-1], and L as L = [lambda][T.sup.2]I, where I denotes the identity tensor. This leads for isotropic systems to the Fourier law

q = -[lambda][nabla]T. (2.1)

This equation implies an instantaneous response of the system. This is in general a good approximation, but it breaks down when the perturbation is very fast or when the response time is very long. In these situations, a flux J(t) at a given time t is related to its respective conjugate force X (t') at previous times t' through

J(t) = [[integral].sup.0.sub.-[infinity]] K(t - t') x X(t')dt', (2.2)

where K(t-t') is a memory function. When K(t-t') = L[delta](t-t') one recovers equations of the form (2.1); the next simplest possibility is to write

J(t) = [[integral].sup.0.sub.-[infinity]] L/[tau] exp[-(t - t')/[tau]] x X(t')dt'. (2.3)

In this case the memory decays exponentially, but it could also decay in other fast ways according to Coleman and Noll's fading memory principle (Wang 1964). The relation (2.3) is equivalent to assuming that the flux is described by a generalized transport equation of the form

[tau] [partial derivative]J/[partial derivative]t = -(J - L x X). (2.4)

Indeed, integration of (2.4), and assuming that the initial value of the flux is zero, leads to (2.3). In the case of heat transport, (2.4) takes the explicit form

[tau] [partial derivative]q/[partial derivative]t = -(q + [lambda][nabla]T), (2.5)

which is generally referred to as the Cattaneo or Maxwell-Cattaneo equation. The term containing the time [tau] represents the heat flux relaxation. When [tau] is negligible or when the time variation of the heat flux is slow, (2.s) will reduce to Fourier's law.

Analogous generalizations are known for other fluxes, but here we will focus our attention on the heat flux. To examine the physical consequences of (2.s), we first consider the behaviour of thermal waves in rigid solids and incompressible fluids at rest. To do this, we use the balance equation of the specific internal energy u, which can be cast in the form

[rho] [partial derivative]u/[partial derivative]t = -[nabla] x q, (2.6)

where du = [c.sub.v]dT, with [c.sub.v] being the heat capacity per unit mass at constant volume. Introducing (2.5) into (2.6) one obtains

[tau] [[partial derivative].sup.2]T/[partial derivative][t.sup.2] + [partial derivative]T/[partial derivative]t - [chi] [[nabla].sup.2]T = 0, (2.7)

where [chi] = [lambda]/[rho][c.sub.v] designates the heat diffusivity, wherein the coefficients [tau] and [lambda] are assumed to be constant. The solutions of (2.7) propagate with finite speed [([chi]/[tau]).sup.1/2], and the amplitude of the pulse decreases in time as a consequence of the irreversible term. For t<<[tau], the first term of equation (2.7) is dominant, so that it reduces to

[tau] [[partial derivative].sup.2]T/[partial derivative][t.sup.2] = [chi] [[nabla].sup.2]T. (2.8)

This is a wave equation with the same wave speed as (2.7); it describes a reversible process as it is invariant with respect to time inversion. In contrast, for time scales longer than [tau](t >> [tau]), the first term of (2.7) is negligible and one obtains the partial differential equation

[partial derivative]T/[partial derivative]t = [chi][[nabla].sup.2]T, (2.9)

associated with diffusion of heat, which is an irreversible process as (2.9) is not invariant when t is changed into -t. Thus, at short times the transport equation (2.5) is reversible and heat propagates with a finite speed (which may be microscopically interpreted as a ballistic motion of the particles before experiencing a collision), whereas at longer times the process becomes irreversible and heat is diffused throughout the system. Then [tau] can be interpreted as the characteristic time for the crossover between ballistic motion and the onset of diffusion. In the context of chaotic deterministic systems, [tau] may be interpreted as the Lyapunov time beyond which predictivity is lost (Nicolis and Prigogine 1989). Thus, from a conceptual perspective, (2.5) is more suitable than the Fourier equation (corresponding to [tau] = 0), as it describes the transition from reversible to irreversible behaviour.

The dynamical properties of equation (2.7) have been thoroughly analysed. For instance, by assuming a solution of the form T(x, t) = [T.sub.0] exp[i(kx-[omega]t)], where [omega] is a (real) frequency and k a (complex) wave-number, it is found that the solution is characterized by a phase speed [v.sub.p], and an attenuation factor [alpha], respectively given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.10)

In the high-frequency limit ([tau][omega] >> 1), the phase speed and attenuation coefficient following from (2.10) are

[v.sub.p[infinity]] [equivalent to] = [([chi]/[tau]).sup.1/2], [[alpha].sub.[infinity]] = 2[([chi][tau]).sup.1/2]. (2.11)

The value of [v.sub.p], diverges when the relaxation time vanishes, thus leading to an infinite speed of propagation. Therefore, the presence of the relaxation time is much more than a minor quantitative correction, as in the high-frequency regime it leads to a completely different behaviour than that predicted by the classical Fourier's law. The analysis of thermal waves has been the topic of much research (Joseph and Preziosi 1989, 1990; Dreyer and Struchtrup 1993; Tzou 1997, Jou et al. 2001a). After application of a thermal pulse, and according to (2.s), the front propagates at a speed [v.sub.p] = [([chi]/[tau]).sup.1/2]; beyond the wavefront, located at [x.sub.f] = [v.sub.p]t, the temperature perturbation is not felt. The temperature behind the wavefront will be higher than the one corresponding to Fourier's law, because the same amount of energy is confined into a smaller volume. For longer times, the difference between the hyperbolic and the parabolic temperature distributions becomes smaller and smaller. Long relaxation time of the heat flux is also typical of superfluids, another topic of interest for extended thermodynamics (Mongiovi and Jou 2007).

Cattaneo's equation (2.5) is appealing in practice, but it is not compatible with the positive definite character of the local-equilibrium entropy production (Jou et al. 2001a; Muller and Ruggeri 1997). Since (2.5) may be considered as an evolution equation for the heat flux, it becomes an independent quantity. To achieve a consistent thermodynamic framework where Cattaneo's equation becomes compatible with the second law, we consider q as an independent variable and assume that the entropy density s is q-dependent. More explicitly, we write:

ds = [T.sup.-1]du - [T.sup.-1]v[[alpha].sub.1]q x dq, (2.12)

where v is the specific volume (the reciprocal of mass density [rho]) and [[alpha].sub.1] a coefficient allowed to depend on T but not on q. From (2.12) and the internal energy balance law (2.6), one obtains for the time derivative s of the entropy

p[??] = -[T.sup.-1][nabla] x q - [T.sup.1][[alpha].sub.1]q x [??], (2.13)

or, equivalently,

p[??] = -[nabla] x ([T.sup.1]q) + q x ([nabla][T.sup.-1] - [T.sup.-1][[alpha].sub.1]q x [??]). (2.14)

This equation can be cast in the general form of a balance equation, namely

p[??] = -[nabla] x [J.sup.s] + [[sigma].sup.s], (2.15)

with [J.sup.s] the entropy flux and [[sigma].sup.s] the entropy production (a positive definite quantity as required by the second law of thermodynamics) given by

[J.sup.s] = [T.sup.-1]q (2.16)

and

[[sigma].sup.s] = q x ([nabla][T.sup.-1] - [T.sup.-1][[alpha].sub.1][??]) [greater than or equal to] 0. (2.17)

Relation (2.17) is a bilinear form [[sigma].sup.s] = q x X in the flux q and the force X, identified as the quantity within parentheses in (2.17); it differs from the classical thermodynamic force, [nabla][T.sup.-1], by the presence of a term in the time derivative of the heat flux.

The simplest way to obtain an evolution equation for q compatible with the positiveness of [[sigma].sup.s] is to assume that the force X is linear in q, namely

[nabla][T.sup.-1] - [T.sup.-1][[alpha].sub.1][??] = [[mu].sub.1]q, (2.18)

where the phenomenological coefficient [[mu].sub.1] may depend on u but not on q because, as previously, third-order contributions in q are omitted. Introduction of (2.18) into (2.17) results in [[sigma].sup.s] = [[mu].sub.1]q x [??] [greater than or equal to] 0, from which it infers that [[mu].sub.1] > 0.

Expression (2.18) contains two coefficients, [[alpha].sub.1] and [[mu].sub.1], which must be identified on physical grounds. By comparing (2.18) with Cattaneo's equation (2.5), one is led to [[mu].sub.1] = [([lambda][T.sup.2]).sup.-1] and [a.sub.1] = [tau]/[lambda]T. With these identifications, the generalized Gibbs equation (2.12) takes the form (Jou et al. 1988, 1999, 2001a)

ds = 1/T du - [tau]/[rho][lambda][T.sup.2] q x dq, (2.19)

wherein it is important to observe that the coefficient of the new term in dq is identified in terms of physical quantities, namely the relaxation time [tau] and the heat conductivity [lambda]. After integration of (2.19), the explicit expression for the entropy outside (local) equilibrium up to second-order terms in q is

[rho]s(u, q) = [rho][s.sub.eq](u) - 1/2 [tau]/[lambda][T.sup.2] q x q. (2.20)

Note that when the relaxation time tends to zero, Cattaneo's equation (2.7) tends to Fourier's equation and (2.19) tends to the classical Gibbs equation of local-equilibrium theory. An alternative way to achieve thermodynamic consistency of the Maxwell-Cattaneo equation is by introducing a vectorial internal variable and obtaining its evolution equation; when the internal variable is proportional to the heat flux, the evolution equation for the latter is seen to be the Maxwelll-Cattaneo equation (Coleman et al. 1982; Morro and Ruggeri 1987; Cimmelli and Kosinski 1991; Verhas 1997).

3. SECOND-ORDER NON-LOCAL EFFECTS

Non-local effects are especially important to describe the transition from diffusive to ballistic regime (Jou et al. 2005; Alvarez and Jou 2007). Indeed, when the mean-free path between successive collisions becomes long, there will be a direct connection between nonadjacent regions, with very different values of the temperature. The surge of interest in nanotechnology has raised new questions about heat transport. In particular, it may be asked whether laws like those of Fourier and Cattaneo remain applicable for nanosystems, where the ratio between the mean-free path and the characteristic length L of the system, the Knudsen number (Kn [equivalent to] l/L), becomes comparable to or higher than 1 and heat transport is no longer diffusive but ballistic. To formulate the problem in simple terms, consider a one-dimensional system of length L, whose opposite boundaries are at temperature T and T - [DELTA]T. Depending on the values of Kn, the heat flux q takes the following limiting forms

q = [lambda] [DELTA]T/L (Kn << 1, diffusive transport), (3.1)

q = [LAMBDA][DELTA]T (Kn >> 1, ballistic transport). (3.2)

The factor [lambda] denotes the thermal conductivity and [LAMBDA] is a heat conduction transport coefficient. In the diffusive limit, the heat flux is proportional to the temperature gradient, according to Fourier's law; in the ballistic regime, it depends only on the temperature difference, but not on the length L of the system.

A simple phenomenological modelling of the transition between the diffusive and the ballistic regimes can be achieved by introducing a heat conductivity [lambda]L(T, l/L) in such a way that in whole generality

q = [lambda](T, l/L)[DELTA]T/L. (3.3)

The limiting values of this generalized conductivity should be

[lambda](T, l/L) [right arrow] [lambda](T) for l/L [right arrow] 0, (3.4a)

[lambda](T, l/L) [right arrow] [lambda](T)/[alpha] L/l [equivalent to] [LAMBDA](T)L for l/L [right arrow] [infinity], (3.4b)

where a is a constant depending on the system. First of all, we introduce second-order non-local effects, which describe weak non-locality (Jou et al. 2001a; Lebon et al. 2008); we will see the possible advantages and the limitations of this situation, but we will show that it does not describe the transition to the ballistic regime. To deal with this transition, an infinite number of higher-order fluxes must be incorporated.

To describe non-local effects, EIT introduces a new extra variable, the flux of the heat flux, described by a second-order tensor Q, and writes instead of the Cattaneo equation (2.7) the following expression (Jou et al. 2001a):

The tensor Q, assumed to be symmetric, may be split in the usual form Q = QI + [Q.sup.s], the scalar Q being one-third of its trace and [Q.sup.s] the deviatoric part. In a relaxational approach, the evolution equations for Q and [Q.sup.s] may be written as

[[tau].sub.0] [partial derivative]Q/[partial derivative]t = -Q + [beta]'[nabla] x q, (3.6)

[[tau].sub.2] [partial derivative][Q.sup.s]/[partial derivative]t = - [Q.sup.s] + 2[beta]"[([nabla][??]).sup.s]. (3.7)

Assuming that the relaxation times [[tau].sub.0] and [[tau].sub.2] are negligibly small and considering only regular solutions for which the time derivatives appearing at the left hand side of (3.6) and (3.7) do not diverge, and substituting (3.6-3.7) into (3.5), one obtains for the evolution of the heat flux

[[tau].sub.1] [partial derivative]q/[partial derivative]t = -(q + [lambda][nabla]T) + [beta]"[[nabla].sup.2]q + ([beta] + 1/3 [beta]") [nabla] ([nabla] x q). (3.8)

Equation (3.8) is comparable with that obtained by Guyer and Krumhansl (1966) from phonon kinetic theory. The phonons in a solid undergo two types of collisions: resistive R-collisions, with defects of the lattice or the boundaries of the crystal, and Umklapp phonon-phonon processes. The former conserve energy but not momentum and have a characteristic time [[tau].sub.R]. The latter, referred to as normal N processes, conserve energy and momentum and have a characteristic time [[tau].sub.N]. Starting from the Boltzmann equation, Guyer and Krumhansl (1966) derived an equation of the form (3.8) with the following identifications of the coefficients:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.9)

from which follows that [[tau].sub.1] is associated with resistive collisions while the non-local contributions find their origin in the normal collisions. Equation (3.8) contains as particular cases the Fourier and the Cattaneo equations for heat transport, for [[tau].sup.1.sub.R] [right arrow] [infinity] and [[tau].sup.1.sub.R] [right arrow] [infinity], and for 1/[[tau].sub.R] finite and 1/[[tau].sub.N] [right arrow] [infinity], respectively. It becomes also clear that the flux of the heat flux Q accounts for the presence of the momentum-preserving N phonon-phonon collisions. Equation (3.8) is also used to describe the so-called phonon hydrodynamics, which arises when the resistive collision time is very long and the normal relaxation time is the predominant one. In this case, (3.8) is analogous to the linearized Navier-Stokes equation, with the heat flux instead of the velocity, [c.sup.2.sub.0][[tau].sub.N]/5 instead of the kinematic viscosity, and [nabla]T playing a role analogous to the pressure gradient in hydrodynamics.

An equation of the form of (3.8) arises also in some other contexts, as for instance in systems composed of two subsystems, each of them with its own temperature. These subsystems are, for instance, the electrons and the lattice in a metal submitted to a short-pulse laser heating, where the electron temperature is much higher than the lattice temperature. This problem is well described by the following evolution equations for the electron and lattice temperatures [T.sub.e] and [T.sub.l] respectively:

[c.sub.e][partial derivative][T.sub.e]/[partial derivative]t = [nabla] x ([lambda][nabla][T.sub.e]) - C([T.sub.e] - [T.sub.l), (3.10)

[c.sub.l] [partial derivative][T.sub.l]/[partial derivative]t = C([T.sub.e] - [T.sub.l). (3.11)

The constant C refers to the electron-phonon coupling, which accounts for the energy transfer from the electrons to the lattice, and [c.sub.e] and [c.sub.l] are the specific heats of the electrons and lattice per unit volume, respectively. When the solution of (3.11) is introduced into (3.10), one is led to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.12)

Such an equation can also be obtained by eliminating q between the energy balance equation (2.6) and equation (3.8) under the following identifications: [c.sub.l]/C = ([beta]"/[lambda])[rho][c.sub.v], [c.sub.e][c.sub.l]/C = 9[rho][c.sub.v][[tau].sub.R], and [beta]' = 0.

Another context where an equation of the form of (3.8) may be derived is the so-called dual-phase-lag equation (Tzou 1997)

[[tau].sub.1][??] + q = -[lambda]([nabla]T + [tau]' [partial derivative]T/[partial derivative]t), (3.13)

which takes into account relaxation effects in both the heat flux and the temperatare gradient. However, such a procedure yields a parabolic equation, while by starting with the heat flux and the flux of the heat flux as independent variables one is led to a hyperbolic equation that contains the parabolic equation (3.8) in the limit of a vanishing relaxation time of the flux of the heat flux.

The thermodynamics underlying the transport equations (3.5-3.7) is easily derived by following the same procedure as in Section 2, i.e. by introducing an extended Gibbs equation and entropy flux of the forms

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.14)

and

[J.sup.s] = 1/[theta] q + 1/[lambda][T.sup.2][Q.sup.s] x q + 1/[lambda][T.sup.2] Qq. (3.15)

Note, incidentally, that (3.15) agrees with the general form derived by Verhas (1983), according to which the entropy flux is the sum of all fluxes each of them multipied by the suitable intensive quantities appearing in the entropy expression. Combining (3.8) with the energy balance equation and working in the Fourier space, the dispersion relation between the frequency [omega] and the wavenumber k for thermal waves is

i[omega] = -[chi][k.sup.2]/1 + [omega] + i[omega][[tau].sub.R] + [l.sup.2][k.sup.2]/1 + [omega][[tau].sub.N] (3.16)

with [l.sup.2] = 4/3 [beta]" + [beta]' for longitudinal waves and [l.sup.2] = [beta]" for transverse waves. For [omega][[tau].sub.R] >> 1 >> [omega][[tau].sub.N] and l = 0 equation (3.16) leads to the phase velocity [v.sub.p] = [([chi]/[[tau].sub.R]).sup.1/2]. This is the regime referred to as second sound. At higher frequencies, when [omega][[tau].sub.N] >> 1 and [omega][[tau].sub.R] >> 1, the phase velocity obtained from (3.16) for longitudinal waves is

[v.sub.p] = [([chi]/[[tau].sub.R] + [l.sup.2]/[[tau].sub.R][[tau].sub.N]).sup.1/2]. (3.17)

This regime is typical of ballistic propagation of phonons, i.e. a flow in which phonons travel through the whole crystal without suffering any collision.

Equation (3.16) may be interpreted in terms of an effective thermal conductivity given by

[lambda](T, [omega], k) = [lambda](T)/1 + i[omega][[tau].sub.R] + [l.sup.2][k.sup.2]/1 + i[omega][[tau].sub.N]. (3.18)

To examine the transition to a ballistic regime in a steady regime, we assume a zero frequency, and consider that k is of the order of 1/L. In this case, (3.18) becomes

[lambda](T, l/L) = [lambda](T)/1 + a[(l/L).sup.2]. (3.19)

This reduces to [lambda](T) for small values of l/L, but it behaves as [L.sup.2] for high values of l/L, in contrast with the observed data, showing that [lambda] is proportional to L in the ballistic regime, as indicated by (3.4). Therefore, introduction of second-order non-local terms is not sufficient to describe the transition. We will see that the introduction of an infinite number of fluxes leads to the required behaviour.

4. HIGHER-ORDER FLUXES, CONTINUED-FRACTION EXPANSIONS, AND BALLISTIC TRANSPORT

So far, we have taken as non-equilibrium variables the heat flux q and the flux of the heat flux Q. Now, we will incorporate higher-order variables, each one being defined as the flux of the preceding one, and denoted as [J.sup.(1)], [J.sup.(2)], ..., [J.sup.(n)], where [J.sup.(n)], a tensor of the order n, is the flux of [J.sup.(n-1)]. Indeed, it follows from the kinetic theory that the relaxation times of the higher-order fluxes are not shorter than the collision time. Therefore, using only the first-order fluxes as independent variables is not satisfactory to describe high-frequency processes, because when the frequency becomes comparable to the inverse of the relaxation time of the first-order flux, all the higher-order fluxes will behave also like independent variables and must be incorporated in the formalism.

Up to the nth-order moment, the Gibbs equation takes the form (Jou et al. 2001a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.1)

and the entropy flux can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.2)

The symbol [cross product] denotes the complete contraction of the corresponding tensors. We have limited ourselves to the simplest form of the entropy and the entropy flux and have not taken into account cubic coupled contributions. The corresponding evolution equations compatible with positive entropy production are

[[tau].sub.1][[??].sup.(1) = -([J.sup.(1) + [lambda][nabla]T) + [[beta].sub.1][[tau]].sub.1]/[[alpha].sub.1] [nabla] x [J.sup.(2)], (4.3a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.3b)

Since [J.sup.(n)] is the flux of [J.sup.(n-1)], this implies, by the very definition of a flux requiring that [[rho][[??].sup.(n-1)] = [nabla] x [J.sup.(n)] + ..., that the coefficient of [nabla] x [J.sup.(n)] in the equation for [rho][[??].sup.(n-1)] must be equal to -1. As a consequence, [rho][[beta].sub.n] = -[[alpha].sub.n]. Several mathematical aspects of these equations and related ones have been dealt with by Cimmelli and Van (2005).

In the ([omega], k)-Fourier space the hierarchy (4.3a, b) may be written as a generalized thermal conductivity [lambda]([omega], k) given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.4)

with [l.sup.2.sub.n] = [[omega].sup.2.sub.n][([[mu].sub.n][[mu].sub.n+1]).sup.-1] > 0. The result (4.4) allows us to define a k-dependent thermal conductivity [lambda]](k) in the steady state. An especially compact expression for [lambda](k) may be obtained when [l.sub.n] does not depend on n, as for instance assuming [l.sup.2.sub.n] = (1/4) [l.sup.2], it is found that for [omega] = 0,

[lambda](k) = 2[[lambda].sub.0]/[(lk).sup.2] [[square root of (1 + [(lk).sup.2])] -1]. (4.5)

From (4.5) and in the steady state [omega] = 0, one obtains a k-dependent thermal conductivity [lambda](T, k). Since in the situation investigated here the system is characterized by one single length scale, the length L, it seems natural to identify k with the inverse of the length of the system, 1/L. This yields the following for [lambda](T, l/L) (Jou et al. 2005):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.6)

which exhibits the required asymptotic behaviours (3.4). This result is in agreement with the experimental data for nanotubes and in some theoretical unidimensional models of heat transport (Alvarez and Jou 2007). The above considerations confirm that EIT provides a consistent modelling for heat transfer processes taking place not only in short times but also at micro- and nanoscales.

5. NON-LINEAR HIERARCHY OF TRANSPORT EQUATIONS AND FLUX LIMITERS

A direct consequence of the finite speed of thermal signals is that, for a given energy density, the heat flux cannot reach arbitrarily high values but will be bounded by a saturation value of the order of the energy density times the maximum speed of propagation. Typical situations arise, for instance, in radiative heat transfer, where the speed of photons is the speed of light c, and where the maximum heat flux is a[T.sup.4]c, or in plasma physics, where the maximum heat flux is [k.sub.B]T[([k.sub.B]T/m).sup.1/2]. Such bounds on heat transfer imply a drastic reduction with respect to the values predicted by Fourier's law at high temperature gradients and play an important role, for instance, in laser-plasma interaction in laser-induced nuclear fusion, or in the collapse of stars.

This saturation effect cannot be described by the classical Fourier law, but can be interpreted by introducing an effective thermal conductivity [lambda](T, [nabla]T), which depends not only on the temperature but also on the temperature gradient. For instance, a typical dependence giving rise to a limited heat flux is

q = - [[lambda].sub.0](T)/[square root of (1 + a[(l [nabla] ln T).sup.2] [nabla T, (5.1)

with l the mean-free path and a a numerical factor. This law reduces to Fourier's law for small temperature gradients and yields a saturation value in the limit of very high temperature gradients. The form of [lambda](T, [lambda]T) has been the subject of many works in radiation hydrodynamics, and several expressions have been proposed (Levermore 1996; Jou et al. 2001a; Anile et al. 2003).

The idea of continued-fraction expansions presented in sections 3 and 4 may be generalized as follows, by taking into account that, in the presence of a temperature gradient, the fluxes of order m, m + 1, and m - 1 may be coupled not only by the wave-vector k but also by the temperature gradient itself; as a consequence the hierarchy of equations (4.3) may be generalized as follows (Zakari 1997):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.2)

In the Fourier-Laplace transform one may write (5.2) as

[[tau].sub.m][[??].sup.(m)] + [J.sup.(m)] = [B.sub.m] x [J.sup.(m+1)] + [C.sub.m][J.sup.(m-1)], (5.3)

wherein [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. This results in a continued-fraction expansion, analogous to that found in the previous section, of the form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.4)

If all the terms [B.sub.i] and [C.sub.i+1] are equal, and considering vanishing values for the frequency and the wavevector, (5.4) will yield the following asymptotic expression for the constitutive equation for the heat flux:

q = - [lambda][nabla]T/1/2 + [square root of (1/4 + [l.sup.2] [([nabla] ln T).sup.2])]. (5.5)

This equation predicts a finite value for the heat flux even when the temperature gradient diverges. The corresponding finite value is given by [q.sub.max] [approximately equal to] [lambda]/l = [rho]uV, where [rho]u is the energy density per unit volume and V the maximum speed for heat pulses.

6. CONCLUSIONS

The current frontiers of technology ask for transport equations far more general than the classical ones, incorporating memory, non-local, and non-linear terms. On the other hand, physical consistency demands that these equations be compatible with the second law of thermodynamics, but this requires understanding how to formulate the second law in the extreme situations considered in technology. Extended Irreversible Thermodynamics provides a limited but fruitful framework to cope with these topics. In this formalism, the fluxes are upgraded to the status of independent variables, and one obtains for them evolution equations satisfying the second law of thermodynamics, formulated in terms of an extended entropy and an extended entropy flux depending themselves on the fluxes. This hierarchy may be expressed as a continued-fraction expansion, which, in the asymptotic case of an infinite number of fluxes leads to a general description ranging from collision-dominated to ballistic dominated regimes. The theoretical predictions agree fairly well with experimental data on heat transport in nanotubes and thin plates.

The consideration of non-local effects requires paying more attention to the boundary conditions and to the meaning of temperature out from equilibrium. Indeed, when non-local effects are important, the influence of the walls is expected to be felt at a considerable distance, and the fact that the mean-free path is long will imply that at each point there is a superposition of populations having different temperatures, thus breaking equipartition. The generalization of the entropy implies also the generalization of the equations of state for temperature, pressure, and chemical potential.

Deeper conceptual questions are found in the case of temperature, which is seen to be a subtler concept than in the classical theory. Indeed, different definitions of temperature in equilibrium may be extrapolated as definitions of effective non-equilibrium temperatures. In equilibrium, all of these temperatures tend to the same value, but out of equilibrium they are different from each other, and a deeper understanding of their relation should be achieved (Caws-Vazquez and Jou 2003).

ACKNOWLEDGEMENTS

We acknowledge financial support from the Direccion General de Investigacion of the Spanish Ministry of Education and Science under grant FIS2006-12296-C02-01 and of the Direccio General de Recerca of the Generalitat of Catalonia, under grant 2005 SGR-00087. G.L was supported by NATO Collaborative Grant PST.CLG 979012.

Received 5 February 2008, in revised form 3 April 2008

REFERENCES

Alvarez, F. X. and Jou, D. 2007. Memory and the non-local effects in heat transport from diffusive and ballistic regimes. Appl. Phys. Lett., 90, 083109 (3 pages).

Anile, A. M., Allegretto, W. and Ringhofer, C. 2003. Mathematical Problems in Semi-conductor Physics. Lecture Notes in Mathematics, Vol. 1823. Springer, Berlin.

Casas-Vazquez, J. and Jou, D. 2003. Temperature in nonequilibrium states: a review of open problems and current proposals. Rep. Prog. Phys., 66, 1937-2023.

Cimmelli, V. A. and Kosinski, W. 1991. Non-equilibrium semi-empirical temperature in materials with thermal relaxation. Arch. Mech., 43, 753-767.

Cimmelli, V. A. and Van, P. 2005. The effect of non-locality on the evolution of higher order fluxes in nonequilibrium thermodynamics. J. Math. Phys., 46, 112901 (15 pages).

Coleman, B. D., Fabrizio, M. and Owen, D. R. 1982. The thermodynamics of second sound in crystals. Arch. Rat. Mech. Anal., 31, 135-158.

Dreyer, W. and Struchtrup, H. 1993. Heat pulse experiments revisited. Continuum Mech. Thermodyn., 5, 3-50.

Eu, B. C. 1992. Kinetic Theory and Irreversible Thermodynamics. Wiley, New York.

Eu, B. C. 1998. Nonequilibrium Statistical Mechanics: Ensemble Method. Kluwer, Dordrecht.

Guyer, R. A. and Krumhansl, J. A. 1966. Solution of the linearized Boltzmann phonon equation. Phys. Rev., 148, 766-778.

Herlach, D., Galenko, P. and Holland-Moritz, D. 2007. Metastable Solids from Undercooled Melts. Pergamon Materials Series, Elsevier, Amsterdam.

Joseph, D. D. and Preziosi, L. 1989. Heat waves. Rev. Mod. Phys, 61, 41-74.

Joseph, D. D. and Preziosi, L. 1990. Heat waves. Addenda, 62, 375-392.

Jou, D., Casas-Vazquez, J. and Lebon, G. 1988. Extended irreversible thermodynamics. Rep. Prog. Phys, 51, 1105-1179.

Jou, D., Casas-Vazquez, J. and Lebon, G. 1999. Extended irreversible thermodynamics revisited (1988-1998). Rep. Prog. Phys, 62, 1035-1142.

Jou, D., Casas-Vazquez, J. and Lebon, G. 2001a. Extended Irreversible Thermodynamics, 3rd ed. Springer, Berlin.

Jou, D., Casas-Vazquez, J. and Criado-Sancho, M. 2001b. Thermodynamics of Fluids under Flow. Springer, Berlin.

Jou, D., Casas-Vazquez, J., Lebon, G. and Grmela, M. 2005. A phenomenological scaling approach for heat transport in nano-systems. Appl. Math. Lett., 18, 963-967.

Lebon, G., Jou, D. and Casas-Vazquez, J. 2008. Understanding Non-equilibrium Thermodynamics. Foundations, Applications, Frontiers. Springer.

Levermore, C. D. 1996. Moment closure hierarchy for kinetic theories. J Stat. Phys., 83, 1021-1065.

Luzzi, R., Vasconcellos, A. R. and Ramos, J. G. 2001. Statistical Foundations of Irreversible Thermodynamics. Teubner, Leipzig.

Luzzi, R., Vasconcellos, A. R. and Ramos, J. G. 2002. Predictive Statistical Mechanics: A Non-equilibrium Ensemble Formalism. Kluwer, Dordrecht.

Maugin, G. A. 1999. The Thereto-mechanics of Nonlinear Irreversible Behaviours. World Scientific, Singapore.

Maugin, G. A. and Muschik, W. 1994a. Thermodynamics with internal variables. Part I. General concepts. J NonEquilib. Thermodyn., 19, 217-250.

Maugin, G. A. and Muschik, W. 1994b. Thermodynamics with internal variables. Part II. Applications. J. NonEquilib. Thermodyn., 19, 250-289.

Mongiovi, M. S. and Jou, D. 2007. Thermodynamcal derivation of a hydrodynamical model of inhomogeneous superfluid turbulence. Phys. Rev. B, 75, 024507.

Morro, A. and Ruggeri, T. 1987. Second sound and internal energy in solids. Int. J. Non-Linear Mech, 22, 27-36.

Muller, I. and Ruggeri, T. 1997. Rational Extended Thermodynamics. Springer, New York.

Nicolis, G. and Prigogine, I. 1989. Exploring Complexity. Freeman, New York.

Ottinger, H. C. 2005. Beyond Equilibrium Thermodynamics. Wiley Interscience, New York.

Tzou, D. Y. 1997. Macro-to-Microscale Heat Transfer. The Lagging Behaviour. Taylor and Francis, New York.

Verhas, J. 1983. On the entropy current. J. Non-Equilib. Thermodyn., 8, 201-206.

Verhas, J. 1997. Thermodynamics and Rheology. Kluwer, Dordrecht.

Wang, C. C. 1964. The principle of fading memory. Arch. Rat. Mech. Anal., 18, 343-366.

Zakari, M. 1997. A continued fraction expansion for flux limners. Physica A, 240, 676-684.

David Jou (a) *, Jose Casas-Vazquez (a), and Georgy Lebon (b)

(a) Universitat Autonoma de Barcelona, 08193 Bellaterra, Catalonia, Spain

(b) Universite de Liege, Sart Tilman, B-4000 Liege, Belgique

* Corresponding author, David.Jou@uab.cat

Printer friendly Cite/link Email Feedback | |

Title Annotation: | PHYSICS |
---|---|

Author: | Jou, David; Casas-Vazquez, Jose; Lebon, Georgy |

Publication: | Proceedings of the Estonian Academy of Sciences |

Date: | Sep 1, 2008 |

Words: | 6639 |

Previous Article: | Preface. |

Next Article: | Objective time derivatives in nonequilibrium thermodynamics/Materiaalsed ajatuletised poordumatute protsesside termodunaamikas. |