Macroscopic Entropy of Non-Equilibrium Systems and Postulates of Extended Thermodynamics: Application to Transport Phenomena and Chemical Reactions in Nanoparticles

In this work, we consider extended irreversible thermodynamics in assuming that the entropy density is a function of both common thermodynamic variables and their higher-order time derivatives. An expression for entropy production, and the linear phenomenological equations describing diffusion and chemical reactions, are found in the context of this approach. Solutions of the sets of linear equations with respect to fluxes and their higher-order time derivatives allow the coefficients of diffusion and reaction rate constants to be established as functions of size of the nanosystems in which these reactions occur. The Maxwell-Cattaneo and Jeffreys constitutive equations, as well as the higher-order constitutive equations, which describe the processes in reaction-diffusion systems, are obtained.


Introduction
Actually, irreversible processes which occur in micro and nanostructures attract more attention of researchers. Extensive development of nanotechologies and creation and construction of devices at micro-and nanoscales has promoted numerous theoretical studies in this field of science [1,2]. One of the most promising routes to study properties of nanosystems and multilayer systems is non-equilibrium thermodynamics. Classical irreversible thermodynamics is based on the local-equilibrium hypothesis, which implies equilibrium in each small volume of the system, although the non-uniform system as a whole remains in non-equilibrium [3][4][5]. Due to the fact that in a locally equilibrium system it is possible to have gradients of state variables, the entropy density in each small volume depends on the internal energy density per unit volume and on the partial densities of the medium components. Such an approach enables to pass from thermodynamic description of equilibrium systems to non-equilibrium ones and to formulate the main theses of irreversible thermodynamics, such as the theorem of minimum entropy production and the general evolutionary criterion [6,7].
If we meanwhile consider the very fast or very abrupt processes, then many properties of non-equilibrium systems cannot be described in the context of local-equilibrium approximation. The validity of local equilibrium is determined by criteria associated to internal microscales of the medium and to external macroscales [8][9][10]. Microscales of a medium are the characteristic scales of the medium microstructure and characteristic relaxation times of dissipative fluxes, which are determined In Sections 7 and 8, the multi-dimensional case of phenomenological equations and the appropriate higher-order constitutive equations, which describe the reaction-diffusion processes in nanosystems, will be considered.

Classical Irreversible Thermodynamics of Heat Conduction, Diffusion and Chemical Reactions
Classical irreversible thermodynamics (CIT) [3][4][5][6][7] is based on the local-equilibrium hypothesis, which leads to the fundamental equation in the local form: where T is the local-equilibrium temperature, µ i are the chemical potential (i = 1, 2, . . . , K); u is the density of the internal energy per unit volume, s is the density of entropy per unit volume, ρ i is the partial density of the i-th component (i = 1, 2, 3, . . . , K, at that all components are assumed to be independent [28]), t is the time.
First, we will consider the classical heat conduction in rigid bodies, which is the simplest thermodynamics theory. Then Equation (1) can be simplified as: that means the entropy density is a function of the internal energy: s = s(u). The balance equation that connects u with the heat flux q is being written, as follows: The fundamental Equation (2) and the balance Equation (3) lead to the standard entropy balance equation (i.e., in the form ∂ t s = −∇ · J s + σ s ): The expression for the entropy flux J s = T −1 q and the entropy production results from Equation (4): where a sign is defined by the second law of thermodynamics. The linear relationship between the flux q and the thermodynamic force ∇T −1 results in one phenomenological equation: which ensures the positivity of entropy production (5) at R q > 0. The phenomenological Equation (6) allows the heat flux to be expressed as: Relation (7), by comparing it with the Fourier's law: defines the coefficient R q : Thus, introducing the linear relationship between the flux and the force, we have adapted the Fourier's law for heat conduction that guaranties the positiveness of entropy production within the framework of CIT.
Next we will consider only dissipative processes [6] at constant density of internal energy u (∂ t u = 0) in the absence of convective motion (diffusion and chemical reactions at v = 0 and dv/dt = 0, where v is the velocity of the multicomponent medium). The fundamental equation in the local form (1) is therefore reduced to the following: As is seen in Equation (9), the entropy density in the considered case is a function of partial densities of the independent components: Any changes in partial densities over time ∂ρ i /∂t obey the balance equations [4,5] where J i is the diffusion flux of the i-th component, w r is the velocity of a chemical reaction r (r = 1, 2, . . . , Q), ν ir is the stoichiometric coefficient of the i-th component in the reaction r. In the present work we will pay attention to the following chemical reactions: The fundamental Equation (9) and the balance Equation (11) define the formalism of classical non-equilibrium thermodynamics. Using Equation (11), we substitute ∂ρ i /∂t in Equation (9). After some transformations, we have the entropy balance equation: where A r = − ∑ i ν ir µ i is the affinity of chemical reaction in Equation (12) (stoichiometric coefficients for the initial components are taken with a negative sign, and those for products are denoted with a positive sign). The balance Equation (13) is written in a standard form, resulting in the entropy flux J s = ∑ i T −1 µ i J i and entropy production: The inequality sign in Equation (14) (positivity of entropy production) is determined by the second law of thermodynamics. Entropy production has a bilinear form, where the vector fluxes J i are multiplied by the appropriate vector forces ∇ T −1 µ i , and the scalar fluxes w r are multiplied by the scalar forces A r /T. According to the Curie principle, there is no interaction between the diffusion fluxes and with the velocity of chemical reaction. Then each term in (14) is non-negative: σ react According to the Equations (15) and (16), each of fluxes J 1 , J 2 , . . . , J K and w 1 , w 2 , . . . , w Q depends only on the corresponding thermodynamic force ∇ T −1 µ 1 , ∇ T −1 µ 2 , . . . , ∇ T −1 µ K and T −1 A 1 , T −1 A 2 , . . . , T −1 A Q , respectively. Further we will focus on the linear correlation between force and flux: where R (1) , R (2) , . . . , R (K) and λ (1) , λ (2) , . . . , λ (Q) are the phenomenological coefficients. Expressing the flux through the thermodynamic force in Equation (17) results in a well-known Fick's law.
Assuming that the chemical potential µ i depends on the component concentration where M i is the molar mass of the i-th component), we Denoting the coefficient before ∇n i through −D 0i , we have Fick's law: In accordance with Equation (19), the concentration gradient ∇n i induces one diffusion flux J i . We have therefore justified the Fick's law in the context of classical non-equilibrium thermodynamics.
The linear dependence of the velocity of chemical reaction on thermodynamic force is being found from Equation (18): On the other hand, the velocity of chemical reaction w r depends on the forward rate is the constant of the forward rate, k − 0r is the constant of the reverse rate. Using the known expression for the chemical potential µ i = µ i 0 + T ln n i , we obtain [7]: where is the universal gas constant. The relationship between the rates of forward and reverse reactions follows from Equation (21): w − r = w + r exp − A r T and the velocity of chemical reaction in the linear approximation is: Comparing Equations (20) and (22), we have: Expressing w + r through the component concentrations, the velocity of reaction is being established from Equation (22) in the form: The analogous equality can be obtained by expressing the velocity of reaction through w − r . The classical thermodynamic theory will be expanded below for both the heat conduction and diffusion and chemical reaction, using in the context various initial postulates of extended irreversible thermodynamics.

Traditional Version of EIT and Higher-Order Heat Fluxes
Extended non-equilibrium thermodynamics is a generalization of classical one by introducing extra variables [11]. In this section we will consider the approach in the context of EIT, which implies that the extra variables are the higher-order fluxes (fluxes of the fluxes): q (1) , q (2) , . . . , q (N) (25) where q (1) = q, q (2) is the flux of the flux q (1) (second-order tensor), q (3) is the flux of q (2) , and so on. Higher-order fluxes (Equation (25)) satisfy a set of the balance equations: where σ (1) , σ (2) , . . . , σ (N) is the production terms of these fluxes. Then the proper fundamental equation will be the generalization of the classical theory Equation (2) where a sign ⊗ means the tensor contraction. By analogy with classical irreversible thermodynamics, the rate of entropy density change is assumed to be representable in the form of a standard equation.
In that case, expressing a flux as: it is easy to obtain the expression for the entropy production: where σ s ≥ 0 obeys the second law of thermodynamics. The expressions in square brackets in Equation (29) are the thermodynamic forces, and the simplest phenomenological equations that satisfy non-negativity of the entropy production are the following linear equations: The coefficients in the phenomenological equations are related to the relaxations times (τ 1 = α 1 /µ 1 , . . . , τ N = α N /µ N ) and correlation lengths (l 2 ). In addition to this fact, α 1 = −β 1 , . . . , α N = −β N . Unlike the unique classical theory of Equation (6), the EIT results in the infinite hierarchy of phenomenological equations at N → ∞ .
Applying the Fourier's transform to the phenomenological Equation (26), we have the following expression for the heat fluxq(ω, k):q whereT(ω, k) is the Fourier's transform of the local-equilibrium temperature, λ(ω, k) is the generalized thermal conductivity coefficient, which can be presented as the continued fraction expansion: where λ 0 (T) is the bulk thermal conductivity, k is the wave vector (we are considering a one-dimensional case). It is assumed that all relaxation times and correlation lengths are equal (τ 1 = . . . = τ N = τ q , l 1 = . . . = l N = l q ) and Equation (32) depends only on one relaxation time τ q and one correlation length l q . Thus, the asymptotic limit is, as follows: If considering only the steady states ( ω → 0 ) and determining the vector k as k = 2πL, where L is the characteristic length, then Equation (33) is reduced to the form: Equation (34) is the required dependence of the generalized thermal conductivity on the parameter l q L = Kn, which is the well-known Knudsen's number. For a stationary state, truncations of Equation (33) are of interest. For example, the second truncation is: where f L/l q is the dimensionless function which takes values from zero to one. Thus, in the context of the considered theory with increasing L/l q , we have λ L/l q → λ 0 , and at L/l q → 0 , we have the linear dependence λ L/l q = λ 0 L/ πl q . The obtained dependence from Equation (34) describes the experimental data well [16].

Higher-Order Time Derivatives of Usual Variables
In this Section we will consider a thermodynamic theory stating that the entropy density is a function of both usual thermodynamic variables and their higher-order time derivatives. It can be considered as the postulate of continuous thermodynamics, as well as a consequence of the local uniformity hypothesis [26], which is less strict if comparing with the local-equilibrium hypothesis. According to the local uniformity hypothesis [26], a system that is non-equilibrium as a whole is considered as uniform in each point, and that is why we may increase the number of state variables. For diffusion and chemical reactions without convection, instead of Equation (3) we have [21,22]: where ∂ t ρ i is the first-order derivative, ∂ t 2 ρ i is the second-order derivative, . . . , ∂ t N−1 ρ i is the (N − 1)-th-order derivative, and i = 1, 2, . . . , K.
Based on Equation (36), we write the fundamental equation which is the generalization of Equation (2), as follows: where θ is the generalized temperature, µ i are the generalized chemical potential (i = 1, 2, . . . , K), Γ i1 , Γ i2 , . . . , Γ iN are the new intensive quantities which correspond to the extra variables (Γ i1 = µ i ). All intensive quantities (θ, µ i , Γ i2 , Γ i3 , . . . , Γ iN , i = 1, 2, . . . , K) depend on both usual and extra variables. The transition to the generalized fundamental Equation (37) means that the local-equilibrium hypothesis is invalid and the thermodynamic system becomes a non-equilibrium as a whole, and at each small element of volume (locally non-equilibrium system). There is meanwhile the thermodynamic formalism identical to classical with preserved basic statements of non-equilibrium thermodynamics [27].
In addition to Equation (11), we will use the equations obtained by differentiating the left and right sides of Equation (11) with respect to time: Like classical thermodynamics, the generalized fundamental Equation (37) and the set of balance Equation (38) allow one to construct the formalism of extended irreversible thermodynamics, namely, to obtain the entropy balance equation, to introduce the generalized fluxes and thermodynamic forces, and to establish a relationship between them. In the extended case, the flux and the entropy production instead Equations (15) and (16) are, as follows: where A r = − ∑ i ν ir µ i is the generalized affinity, B rj = − ∑ i ν ir Γ ij is the analog of chemical affinity corresponding to variables ∂ t j−1 ρ 1 , ∂ t j−1 ρ 2 , . . . , ∂ t j−1 ρ K (j = 1, 2, . . . , N; at j = 1 we have B 1r = A r ). As the result, we obtain a set of phenomenological equations which connect thermodynamic Nr and w r , ∂ t w r , . . . , ∂ t N−1 w r . In the simplest case these quantities are connected by the linear phenomenological equations, which ensure the positivity of entropy production. For the diffusion process, we have: where R (i) αβ are the phenomenological coefficients, α, β = 1, 2, . . . , N. The similar equations can be obtained for the chemical reactions:

Two-Dimensional Case of Phenomenological Equations
Let us further focus on the case N = 2. Then the set of Equation (40) can be written in a much simpler form: It is well known that in EIT the Onsager-Casimir reciprocal relations [4,5] have not been proven. However, we believe that there is the following basis for proving reciprocal relations out of local equilibrium: (i) the proposed theory is based on the local homogeneity hypothesis, which is a less rigorous analogue of the local equilibrium hypothesis; (ii) the classical Onsager-Machlup theory of fluctuations is applicable for systems with memory; (iii) the additional variables considered in the proposed theory are independent variables.
We further assume that Onsager-Casimir reciprocal relations take place in the extended theory. It is obvious that fluxes J i are odd variables. Hence time derivatives ∂ t J i should be viewed as even variables. Therefore, the matrix of phenomenological coefficients in Equations (42)

.
Let us divide Equation (42) by a factor R (i) 11 and make the following substitution by analogy with classical thermodynamics (Equation (19)): being the Maxwell-Cattaneo Equation [11]. We further assume that θ −1 Γ i2 are the homogeneous quantities; i.e., ∇ θ −1 Γ i2 = 0, and the diffusion of the components is defined by only the gradient ∇ θ −1 µ i . Divide the second equation by the factor R that allows it to be written in the matrix form: where: Entropy 2018, 20, 802 10 of 19 As is seen, Equation (44) is the generalized Fick's law (19) with adding a term τ i ∂ t J i . In this case, the gradient ∇n i induces both the diffusion flux J i and the rate of its change ∂ t J i . Equation (44) leads to a delay in the propagation of diffusing components compared to classical thermodynamics. The greater is the relaxation time τ i , the further is the thermodynamic system from the local-equilibrium state. By solving a system of Equations (45) and (46) relative to J i and ∂ t J i using the Kramers formula, we get the following expressions for the flux and its time derivative: where Again, by comparing the first equality in Equation (49) with Fick's law, we get the diffusion coefficient D 0i multiplied by the function f (τ i , τ 0 ) in the extended theory, which takes values from zero to one. In this connection, one can introduce a generalized diffusion coefficient: The greater is the memory effect in the non-equilibrium system, i.e., the greater is τ i , the lower are the diffusion flux J i in Equation (49) and the generalized diffusion coefficient D i in Equation (51). Therefore, the memory effect reduces the diffusion coefficient and, consequently, the rate of diffusion of components. In some cases diffusion may occur in the form of waves [11].
The similar results are obtained for the rate of chemical reaction. In a two-dimensional case (N = 2) a set of Equation (41) is also reduced to two equations: where λ Equation (54) describes the changes over time in the rate of reaction r taking into account the memory effect. By comparing Equations (22) and (54), it is evident that the linear equality (22) is completed with a term τ r ∂ t w r . Then, as for diffusion, we assume B r2 = 0 and divide Equation (53) by a factor λ (r) 11 . It is evident that the ratio λ 11 . Thus, a set of Equations (52) and (53) can be written as: and represented in a matrix notation: where: Solving the linear system of Equations (52) and (53) with respect to w r , ∂ t w r , we get the expressions identical to those for diffusion: where: Taking into consideration the expression for the rate of direct reaction w + r = k + 0r n 1 ν 1r n 2 ν 2r . . . n L ν Lr , we introduce the generalized rate constant for: where the function f (τ r , τ 0 ) takes values from zero to one. A similar expression is valid for the rate constant for the reverse reaction k − r , as well. Comparing the first expression in Equation (58) with a classical one (Equation (22)), it is obvious that the greater is the relaxation time τ r (i.e., the further is the system from the local-equilibrium state), the more is delayed the chemical reaction and the lower are the rates of direct and reverse reactions. Therefore, processes with delay reduce the diffusion coefficients, as well as the rate constant of chemical reactions.

Characteristic Quantities of Nanosystems
The rates of diffusion and chemical reaction depend not only on τ i , τ r , and τ 0 . It is well known [8][9][10] that the relaxation time of each component τ i is related to the characteristic size of microstructure of the medium l i and to the rate of perturbation propagation in the medium V 0i (the rate of diffusion wave propagation): τ i = l i /V 0i . The same dependence is valid for a chemical reaction: τ r = l r /V 0r , where l r is the characteristic microscale corresponding to the chemical reaction, V 0r is the rate of approach to equilibrium due to the chemical reaction. The characteristic time τ 0 is related to the intrinsic macroscale of the medium L and to the rate of change in macroparameters on account of external effects V: τ 0 = L/V (all these quantities depend neither on the nature of components nor on the chemical reaction). Then: It is evident from Equation (61) that the time ratio τ 0 /τ i linearly depends on two dimensionless magnitudes: the linear parameter ratio L/l i and the rate ratio V 0i /V. In a similar way, τ 0 /τ r depends on L/l r and V 0r /V. Once L >> l i at a fixed rate ratio V 0i /V or V << V 0i at a fixed linear parameter ratio L/l i , the thermodynamic system tends to local equilibrium. Otherwise (L ≤ l i at a fixed rate ratio or V ≥ V 0i at a fixed linear parameter ratio), the latter is non-equilibrium. Similar regularities are observed for parameters of the chemical reaction.
In this paper, the quantity L is a linear size of nanoparticles. We are interested to know whether there is a relationship between the parameters L/l i and L/l r and the rate of diffusion of components in the nanoparticles and the rate of chemical reactions in nanoparticles. In this case V 0i /V = const and V 0r /V = const, and we assume for simplicity that V 0i /V = 1 and V 0r /V = 1 and will consider the dependence of the time ratio τ 0 /τ i and τ 0 /τ r on only the parameters L/l i and L/l r , respectively: Using Equation (62), we present the diffusion coefficient (51) and the rate constant of the direct reaction (Equation (60)) in the form: where f (L/l i ) =

Multi-Dimensional Case
First of all, consider the general case of phenomenological equations at N > 2 in assuming that diffusion and chemical reactions are still described by the Maxwell-Cattaneo Equations (44) and (54). Then the matrix equations, which are a generalization of Equations (47) and (57), are being written as: where the vector L 0 N consists of the components (1, 0, . . . , 0), the vector J N includes the components (J i , ∂ t J i , ∂ t 2 J i , . . . , ∂ t N−1 J i ) or (w r , ∂ t w r , ∂ t 2 w r , . . . , ∂ t N−1 w r ), and the matrix M N depends on the relaxation times τ 0 and τ, where τ denotes τ i or τ r . We have therefore omitted the indices i and r in the matrix M N .
If basing on the Maxwell-Cattaneo constitutive equation, then all coefficients starting from the third in the first row of the matrix M N composed of phenomenological coefficients of Equations (40) and (41) vanish. Thus, the first row of the matrix M N has a form (1, τ, 0, . . . , 0), and the first column in M N is written as (1, −τ, 0, . . . , 0). We further assume that the diagonal coefficients in M N depend on only τ 0 , while the non-diagonal ones are the functions of τ. That's why a three-dimensional case is presented by the following expressions: As is seen from Equation (66), we again assume that the Onsager-Casimir reciprocity relations are valid. A function f 3 (τ, τ 0 ) is the second-order truncation of continuous-fraction expansion, which can be presented in the form: Similarly for all odd N, if a function f N (τ, τ 0 ) is representable as finite truncation of Equation (67), it must be assumed that the matrix M N satisfies the Onsager-Casimir relations. As in traditional version of EIT, at τ 0 /τ → 0 the function f N (τ, τ 0 ) tends to a non-vanishing limit. At N → ∞ for f N (τ, τ 0 ) we have the asymptotic function: In the present work we are interested in a relationship between the nanoparticle dispersity and the rate of irreversible processes which occur in nanoparticles. In this connection, we focus on the case when τ 0 /τ = L/l (Equation (62)). Figure 1 displays a typical function at odd N = 3 and an asymptotical function f N→∞ depending on L/l.  Figure 1 displays the functions and * . As is seen in Figure 1, is closer to the asymptotic limit, rather than * . A similar regularity is observed for all even . Then at even N consider M N which results in a function f N (τ, τ 0 ) as the truncation of continuous fraction expansion. Nevertheless, for this condition at N > 3 we have the matrices, for which the Onsager-Casimir reciprocal relations are invalid. For instance, at N = 4: where the function f 4 (τ, τ 0 ) is representable as truncation of continued-fraction expansion. The matrix M N is easily transformed into the matrix M * N (by replacing a sign before some τ), for which the Onsager-Casimir reciprocal relations hold. However, in this case the function M N cannot be presented as truncation of continued-fraction expansion. For example, at N = 4 the replacement of the sign before τ in the fourth row of the matrix M 4 gives: Figure 1 displays the functions f 4 and f * 4 . As is seen in Figure 1, f 4 is closer to the asymptotic limit, rather than f * 4 . A similar regularity is observed for all even N. Consider now the general case N > 2, when the constitutive equation may contain up to N − 1 time derivatives of fluxes and depends on one relaxation time τ: or:  3 (τ, τ 0 ) will take the form: Similarly, the matrix M 4 and the function f 4 (τ, τ 0 ) can be obtained instead of Equation (73): where we assumed that M 3 and M 4 satisfy the Onsager-Casimir relations. Figure 2 shows the functions f  4 (L/l) that seem to tend to a finite value at L/l → 0 . Meanwhile, if some coefficients in these matrices are assumed to be zero, we obtain the functions f  4 (L/l) which tend to zero at L/l → 0 . This is evident from Figure 2 at: (76) which tend to zero at / → 0. This is evident from Figure 2 at: We can see that Equations (73)-(76) cannot be representable as truncations of continued-fraction expansion.

Constitutive Equations
In this section, we consider some constitutive equations that describe the diffusion process in nanosystems. It should be noted that even in a two-dimensional case the obtained constitutive equations are not limited by the Maxwell-Cattaneo Equation (44) and equations containing the higher-order time derivatives of diffusion fluxes in Equations (71) and (72). Let us consider a case when the generalized temperature and chemical potential ̄ depend on and . Then the thermodynamic force ( ̄ ) is written as: Using Equation (77), we make a replacement in equation: that results in the Jeffreys type constitutive equation (dual-phase-lag transfer equation [22,25]) instead of Equation (44), which describes the diffusion of components: We can see that Equations (73)-(76) cannot be representable as truncations of continued-fraction expansion.

Constitutive Equations
In this section, we consider some constitutive equations that describe the diffusion process in nanosystems. It should be noted that even in a two-dimensional case the obtained constitutive equations are not limited by the Maxwell-Cattaneo Equation (44) and equations containing the higher-order time derivatives of diffusion fluxes in Equations (71) and (72). Let us consider a case when the generalized temperature θ and chemical potential µ i depend on n i and ∂ t n i . Then the thermodynamic force ∇ θ −1 µ i is written as: Using Equation (77), we make a replacement in equation: that results in the Jeffreys type constitutive equation (dual-phase-lag transfer equation [22,25]) instead of Equation (44), which describes the diffusion of components: where D 0i = 1 Further, using the integral representation of the diffusion flux, we obtain the higher-order diffusion equations. Let us consider the expression [32,33] for the diffusion flux: where κ is a constant. Let us represent the right-hand side of the Expression (80) in the form Using (80) we obtain from (81) the Maxwell-Cattaneo type equation: where τκ = D 0i Similarly, from (81) we can obtain the equality: which can be cast into equation: The comparison of Equations (84) with (78) allows us to write the right-hand side of Equation (84) in the form 1/R As a result, we obtain the corresponding phenomenological equation: where: In this case the simplest matrix M 4 and function f 4 (τ, τ 0 ) have the form: It is obvious that higher-order equations can also be obtained in the same way. At L/l → 0 one observes a decrease in the rates of diffusion and chemical reactions, and the higher-order constitutive equations can be used for description of irreversible processes which occur in nanosystems.
We have, thus, obtained a wider class of constitutive equations compared to the traditional EIT. For the completeness of thermodynamic theory, the following question is very important: whether these constitutive equations can be transformed into the form which is independent of the frame (material-invariant equations)? Several approaches are known for the Maxwell-Cattaneo equation. The most well-known approaches use the Jaumann derivative instead of the usual local time derivative [11]. The modified version of the objective time derivative was used to solve the problem of Marangoni-Bénard instability [34]. In [35,36], the material derivative and the Oldroid's derivative was used to modify the Maxwell-Cattaneo equation. In this case, it was possible to obtain a single equation for the temperature field [36] and to prove its Galilean invariance.
However, the question remains: is it possible to apply these approaches to other equations considered in this paper? There is a problem even for the Jeffreys Equation (79) where the term D 0i η∇∂ t n i is added, and is more complex for higher-order Equations (71), (72) and (84). In addition, there is an opinion that frame-indifference is not satisfied in several disciplines [11]. Therefore, there is doubt for its applicability in the general case.

Discussion and Conclusions
Thermodynamic theory of nanosystems assumes various approaches, which are based on the introduction of additional variables. One of the first theories is inspired by statistical mechanics, where an ensemble of equivalent systems is considered [37]. In this case, the number of systems is an additional variable. Using this approach, we can generalize the equations of classical thermodynamics and investigate systems at the nanoscale. Unfortunately, transfer processes are not considered in this theory.
Another approach is based on the introduction of internal variables of nanoparticles (volume or size of particles as additional variables) [38,39]. It is assumed that the system is in the local equilibrium and the well-known formalism of thermodynamics of irreversible processes is used. In this case, stochastic methods based on the Fokker-Planck equation and the Smoluchowski equation can be applied.
In this paper, we proceeded from classical thermodynamic theory and considered the traditional version of EIT of heat conduction, wherein the extra variables are the higher-order heat fluxes (fluxes of the fluxes). Within this approach the heat transport in the nanosystems and the thermal conductivity as a function of the size of the nanoscale structures were considered [1,2].
At the same time, assuming that the entropy density is a function of usual thermodynamic variables and their higher-order time derivatives, we studied the irreversible processes which occur in reaction-diffusion nanosystems. The entropy balance equations were obtained, and the linear phenomenological equations were constructed based on the expression of the entropy production. The characteristic relaxation times, which define the memory effects of irreversible processes connected with the coefficients of phenomenological equations, were introduced, as well. Thus, the diffusion coefficients and the rate constant of the chemical reaction were established as the functions of size of nanosystems, using in the context the linear theory. The proposed thermodynamic theory can be developed in the framework of the linear EIT and proceeds from the constitutive equations, at that the matrix of phenomenological coefficients satisfies the Onsager-Casimir reciprocal relations in the case of the odd number of phenomenological equations.
Solutions of a linear system of phenomenological equations relative to fluxes and their higher-order time derivatives depend on the matrix of phenomenological coefficients. As in the case of the traditional version of EIT, a series of matrices leads to the solutions which are representable in the form of finite truncations of continued-fraction expansion. Therefore, obtained solutions with increased number of phenomenological equations converge to the finite limit. The established expressions for fluxes and their time derivatives reveal that the rates of diffusion and chemical reactions decreases with reducing dimensions of nanosystems (Equations (63) and (64)). For a constitutive equation with the higher-order time derivatives of flux one can also construct the linear theory, according to which the greater is the dispersion of nanosystems, the lower are the rates of diffusion and chemical reactions.