1. Introduction
Extended Thermodynamics takes the first steps from the suggestions of kinetic theory of monatomic gases; here, the state of a gas is described by the phase density,
, such that
are the number density of atoms at the point,
, and at time,
t, that have velocities between
and
. The phase density, also called the distribution function, obeys the Boltzmann equation:
where external forces have been neglected and
q takes into account the effects of collisions between the atoms. After that, the moments are defined by:
with
m the molecular mass. Multiplication of the Boltzmann equation by
and integration over all
gives the balance equations for these moments. Some of these equations, and the corresponding moments, may also be considered only through their suitable traces. A corresponding definition can be formulated in the relativistic framework taking, after that, their non-relativistic limit. In this way, we obtain the following balance equations of the many moments model:
where the infinite hierarchy of moments has been closed by choosing two arbitrary numbers,
N and
M, such that
and
are odd numbers.
Other considerations [
1,
2] lead one to think that
N and
M are two subsequent numbers; this aspect does not affect the present work, so that we prefer to maintain their generality.
The quantity, , is an abbreviation for , that is, the tensor, , with , of which we take the trace -times.
Similarly, and .
The moments, and , can be taken as independent variables, while and are their fluxes and and are their production terms. All these quantities are symmetrical with respect to and ; in the particular case of monatomic gases, and are symmetric with respect to all their indexes, but for dense gases and macromolecular fluids, some evolutive terms appear also in the production terms; by taking them in the lefthand side of Equation (1), they may be included in the terms, and , causing the loss of symmetry of and with respect to all their indexes. We remark that the kinetic theory of monatomic gases has been considered only to draw from it suggestions on the properties of moments and balance equations; it is assumed that these same properties hold also for dense gases and macromolecular fluids, except for the above mentioned symmetries.
Moreover, we have
, and , so that Equation (1)1 for , , and the trace of that for are the conservation laws of mass, momentum and energy, respectively. The other production terms, and , are functions of the independent variables, which are zero at equilibrium; in the context of the present macroscopic approach, we may write the linear expressions of these functions with respect to equilibrium and find the sign of the coefficients by imposing that the density of entropy production is not negative; but we omit here these considerations for the sake of brevity.
The following conditions can be considered:
The first one of these conditions surely holds and means that the flux in the conservation law of mass (
i.e., Equation (1)
1 with
) is the independent variable in the conservation law of momentum (
i.e., Equation (1)
1 for
).
The second condition, (2)2, holds only if there is a conservation of angular momentum.
The balance Equation (1) becomes field equations when the fluxes and the productions are known functions of the independent variables. Restrictions on the generality of these functions are furnished by imposing the entropy principle and that of Galilean relativity.
In the particular case
and
, we obtain the 14-moments model, which was firstly studied by Kremer [
3,
4] by imposing, up to the second order with respect to equilibrium, the conditions that come out from the above principles; subsequently, in [
5,
6], an exact non-approximated solution of these same conditions has been found. This result was achieved by writing firstly a relativistic model, for which it is easy to impose the Einstein relativity principle by using the representation theorems for isotropic functions, which are known with this number of variables [
7,
8,
9,
10,
11], and then taking its non-relativistic limit.
An extension of this method, to the case with many moments, is not possible, because there are no known representation theorems for this general case.
In the present paper, we have found a different method to solve the above conditions. We describe these conditions in
Section 2 and solve them in two different ways, which are reported in
Section 3 and
Section 4, respectively. Finally, in
Section 5, we show that, in the particular case of the 14-moments model, they are equivalent to those of [
5,
6], even if these were found with a completely different method.
Obviously, all these considerations belong to the general context of Extended Thermodynamics, whose original and most important papers, in our opinion, are [
12,
13] and whose most important aspects are described in [
14].
Some authors avoid considering Equation (1)
2 and, sometimes, also some components of Equation (1)
1 (see [
15,
16,
17,
18,
19,
20,
21], for example); but, this is not physically acceptable, because it has been shown in [
22,
23] that the form of Equation (1) is suggested by the non-relativistic limit of the relativistic model and that omitting Equation (1)
2 is equivalent to omitting the conservation law of mass or that of momentum energy in the corresponding relativistic model.
A similar result, restricted to the 14-moments model, was previously found in [
24]. See also [
6,
25] for other details.
We conclude this section, remarking that this model has been obtained through a macroscopic approach,
i.e., based on the entropy principle. Another possible formulation of Extended Thermodynamics for dense gases is based on the kinetic approach; in order to permit a comparison between the two approaches, we cite now some references on the kinetic approach. In [
26], Enskog introduced a kinetic theory for dense gases, which yields a very good approximation of the behavior of gases. Later, hydrodynamic-like equations have been derived from the kinetic equation; see, for example, the Chapman-Enskog method [
27].
In 1988, Kremer and Rosa [
28] obtained hydrodynamic equations from the local equilibrium distribution function, as kernel linearizing the collision integral in Enskog’s equation; in this way, they were able to derive sound dispersion relations for monatomic gases by using normal mode analysis. Based on this last paper, in 1991, Marques and Kremer [
29] obtained linearized hydrodynamic equations involving the second order terms of the collision integral; in this way, they improved the results previously known in the literature and, furthermore, they obtained linearized Burnett equations for monatomic gases.
In [
30], Ugawa and Cordero obtained extended hydrodynamic equations derived from Enskog’s equation by using Grad’s moment expansion method in the bi-dimensional case; among other results, they discussed the nature of a simple one-dimensional heat conduction problem and were able to show that, not too far from equilibrium, the non-equilibrium pressure in this case depends on the density, temperature and heat flux vector.
Finally, another model in this context can be found in [
31], and this will surely be the object of further investigations in the future.
2. The Entropy Principle and the Galilean Relativity Principle
The entropy principle states that the supplementary law:
holds for every solution of Equation (1), where
h is the entropy density and
its flux. For Liu’s Theorem [
32], this is equivalent to assuming the existence of Lagrange Multipliers,
and
, such that:
besides a residual inequality, which we leave out for the sake of brevity and because it does not affect the present results.
It is convenient to write Equation (4) in a more compact form by calling
the variables,
and
; in this way, Equation (4) can be written as:
with the obvious meaning of
and
.
In [
33], the idea has been conceived of to define the four-potentials,
and
, as:
so that Equation (5) becomes:
Another idea exposed in [
33] is to take the Lagrange Multipliers as independent variables and call them “main field"; in this way Equation (7) becomes:
so that everything is determined in terms of
and
; moreover, by substituting Equations (8) in (1), we obtain a symmetric system of evolution equations, which is hyperbolic if
is a convex function of the main field. We observe here that it is not necessary to transform the system Equation (1) in the symmetric form, because if it is equivalent to a symmetric hyperbolic system, it itself is hyperbolic. Consequently, we can consider Equations (5)
1 and (7)
2 to still be maintaining the moments as independent variables. In this way, Equations (5)
1 and (7)
2 become:
which allows one to determine the unknown constitutive functions,
, because the matrix,
, is an invertible matrix; otherwise, we could not invert Equation (9)
1 and take the
as independent variables.
Let us introduce now the Galilean Relativity Principle.
To this end, it is firstly necessary to see how the moments and their fluxes transform under a change of frames moving, one with respect to the other, with a translational rectilinear uniform motion with velocity,
. For the variables,
and
, it can be found in [
34,
35] and reads:
where
and
are the counterparts of
and
in the other reference frame; moreover:
Regarding the other variables, we use the identity holding for
:
with
where
S is the set of the couples,
, of integer numbers,
p and
q, such that
,
and
and
.
We omit here the proof of this identity, for the sake of brevity (see [
36] ). Moreover, we define:
Now, we can extend Equation (10) also for the value
and contract the result with
; by using Equations (12) and (14) we obtain:
where, in the second summation, we have changed the index according to
.
Equations (10) and (15) can be written in a more compact form as:
with the obvious meaning of
. This matrix satisfies the following properties:
where
is the constant matrix defined in the following way:
Another property of the matrix,
, is:
The counterparts of these properties when only Equation (1)
1 is considered and Equation (1)
2 is omitted have been already found in [
14]; now, we have found that they hold also for the complete system Equation (1), but we have reported their proofs in a separate paper [
36], for the sake of brevity.
Now, we are ready to impose the Galilean relativity principle for our system. It can be subdivided in two parts:
Vice versa, if we assume that
h and
are non-convective quantities, then Equation (16)
2 will follow as a consequence and, also, the non-convectivity of
. In fact, from Equation (9)
2, it follows:
where, in the last passage, we have used the inverse of Equation (16)
1. Similarly, always with the use of the inverse of Equation (16)
1, we obtain:
Thanks to this equation and to Equation (24), we obtain:
But, from the first equality in Equation (21) written in the new reference frame, we have:
This result allows one to rewrite Equation (25) as:
Now, the hypothesis that
is a non-convective quantity means that:
(Because
,
). By using this equation, we can write Equation (26) as:
Now,
is a non-convective invertible matrix, so it follows that
is a non-convective quantity, and this proves Equation (16)
2.
After that, we note that Equation (22) was deduced without using the non-convectivity of . In other words, we assumed that is a non-convective quantity and proved with Equation (22) that also has this property; similarly, we assume now that is a non-convective quantity and prove through Equation (22) that also satisfies this property.
We conclude that the Galilean relativity principle amounts simply in the conditions (20)1 and (23)2.
2.1. The Galilean Relativity Principle in Terms of the Main Field
A more simple result holds when we take the Lagrange multipliers as independent variables; in fact, from Equations (5) and (16)
1 it follows:
from which, by using Equation (17), it follows:
that is, the equation corresponding to Equation (16)
1, but for the Lagrange multipliers. From Equation (29), it follows:
where Equations (18) and (19) and, again, Equation (29) have been used.
If we start now from the hypothesis that and are non-convective quantities, then the decomposition Equation (16)1 follows as its consequence and of Equations (24)2 and (8).
Similarly, with the use also of Equation (27), we have:
The sum of these two equations gives Equation (16)
2.
Consequently, if we take the Lagrange multipliers as independent variables, the Galilean relativity principle amounts simply in imposing that and are non-convective quantities, that is. and with deduced from Equation (27).
In other words, it is expressed by:
where we have used
Equation (30).
In the next section, we will find the general solution of Equations (20)
1 and (23)
2 in the independent variables
, while in
Section 4, we will find the general solution of Equation (31) in the independent variables
.
Before ending this section, we note that in the earlier papers on Extended Thermodynamics, the independent variables, and , were used, which are defined by Equation (16)1 with instead of , and for .
Similarly, the dependent variables, , were defined by Equation (16)2 with instead of and and instead of .
If we want to express the present results in terms of these other variables, we note that
Equations (8)1 with are still a system of implicit equations for the determination of the Lagrange multipliers in terms of (velocity) and . After that, Equation (31) means that the composite functions, and , do not depend on .
By substituting
into Equation (8)
2, the resulting function,
, satisfy the relation,
. In fact, we have:
where, in the passage denoted with
, we have used Equation (18), in the passage denoted with
, we have used Equation (8) and in the passage denoted with
, we have used Equation (19) and also Equation (30); this last one can be used, because it was deduced, taking into account only the fact that Equation (28) and Equation (29) are invertible relations between
and
, without worrying if these are independent variables or not.
Finally, in the passage denoted with , we have taken into account the fact that the coefficient of inside the square brackets is Equation (31)1 written with E instead of A and, after that, derivated with respect to ; the remaining part is the derivative with respect to of Equation (31)2, written with E instead of A and, after that, derived with respect to .
So, we have proven that the usual decomposition in terms of , and is a consequence of our Equation (31). In this way, the only unknowns are the functions, , and , depending on the variables, .
To find these unknowns, we can use Equation (8) to obtain
as functions of
, substitute them in
,
and
and, after that, calculate this in
, so obtaining
,
and
. However, this is equivalent to calculate Equation (8) in
, so that they become:
After that, we use Equation (32)
1,2 to obtain
and
in terms of
and substitute them in Equation (32)
3 and in the expressions of
and
, calculated in
. This is the essence of the method explained in [
37,
38] and which was already been used in [
39,
40].
3. The General Solution in Terms of the Moments, , as Independent Variables
We find now the general solution of Equations (20)1 and (23)2 in the independent variables, the moments, . Let us begin with Equation (20)1.
To this end, let us call
the variables,
, different from
F and
. With this notation, Equation (20)
1 becomes:
It is not restrictive to consider
h a composite function of:
By using this expression and Equation (18), the Equation (33) becomes:
However, for Equation (19), the first term in the square brackets of Equation (35) becomes:
This allows one to rewrite Equation (35) as:
By using the expression of
reported after Equation (18), our condition becomes
, that is,
H does not depend on
. In other words, the general solution of Equation (20)
1 is:
Let us find now the general solution of Equation (23)
2.
To this end, it will be useful to know firstly the expression of
deduced from the first equality in Equation (21) and from Equation (37), that is:
where we have used Equation (18). Moreover, from Equation (34)
2, we see that the underlined terms contribute with
. Two other terms in Equation (38) are opposite, and there remains:
In other words, also
has an expression like Equation (37), and this is not strange, because we have already said that Equation (23)
1 holds as a consequence of Equation (20)
1; moreover, Equation (23)
1 is the same Equation (20)
1, but with
instead of
h.
Let us find now the general solution of Equation (23)
2, that is
; It is not restrictive to write
which allows one to rewrite our equation as:
where we have taken into account that
satisfies Equation (23)
1. By using Equation (36), there remains
which, by using the expression of
reported after Equation (18), becomes
, that is,
does not depend on
. In other words, the general solution of Equation (23)
2 is:
We note now that Equation (20) looks like Equation (23), but with
instead of
; this allows one to write the most general solution of Equation (20)
2 like Equation (40), that is,
It is interesting to test our general solution in the particular case,
and
, that is, the 14-moments case, where Equations (20)
1 and (23)
2 are:
In this case, the variable Equation (34)
2 is:
and it is easy to verify that,
in composite with Equation (43), are solutions of Equation (42).
4. The General Solution in Terms of the Lagrange Multipliers as Independent Variables
Let us now describe the general solution of Equation (31), using as independent variables. It holds only in the subdomain with , like Equations (20)1 and (23)2 were solved in the previous section, in the subdomain . However, F is the mass density, so that it is obvious that it is different from zero. The same thing cannot be said for , which is zero, quite in equilibrium!
However, if we assume that
in equilibrium is an infinitesimal of a higher order with respect to
, then our solution will be acceptable also at equilibrium. This assumption has been already adopted in [
6,
25] for the 14-moments model and was based on the requirement of the convexity of entropy, which, in turns, assures the hyperbolicity of the resulting system.
Let us call
the variables
, different from
and
. With this notation, Equation (31) becomes:
where we have taken into account that
when
A is the same index of
and
; this fact is evident from the description of
, reported above, after Equation (18).
It is now not restrictive to consider
a composite function of
In this way, Equation (45)
1, by using Equation (18), becomes:
where we have substituted
with
in order not to lose sight the fact that
i is the same index of
.
Now, from the description of
reported above, after Equation (18), we have:
By using these equations, the relation, Equation (47), becomes:
Now, we see that, from the first and last term, in this equation, there remains:
where we have used Equation (48)
2,3. It is now evident that from Equation (49), it remains
, that is,
does not depend on
. In other words, the general solution of Equation (31)
1 is:
It remains now to find the general solution of Equation (31)
2. By defining
from
, it becomes:
where we have used the fact that
already satisfies Equation (31)
1 and taken into account Equation (48). Now, Equation (51) looks like Equation (31)
1, but with
instead of
; by operating in the same way, we obtain that
is a composite function of
and
. In other words, the general solution of Equation (31)
2 is:
Let us now see how the conditions (2) can be solved in terms of the functions,
and
, appearing in Equations (50) and (52), respectively.
4.1. On the Condition
Let us consider now the condition (2)
1, which means that the flux in the conservation law of mass is the momentum density, that is, the independent variable in the subsequent conservation law of momentum. By using Equation (8), it can be expressed as:
However, we note that
λ appears in
, as expressed by Equation (46)
2, only for
, that is, in
; moreover, we have:
Similarly,
appears in
only in
and
; moreover, we have:
In order to prove these properties, we note that from Equation (16)
1, it follows
; by comparing this with the definition (46)
2 of
, we find that
is the coefficient of
in
, calculated in
and
is the coefficient of
in
, calculated in the previous value of
(in effect, this value of
was introduced above when we considered
a composite function of Equation (46). This was inspired by the fact that, from Equation (28)
2 with
, this value
corresponds to the condition
, exactly as for the corresponding passage used in terms of the moments.)
By using Equations (10)
1 and (15)
1, we find:
It follows that the coefficient of
is:
and the coefficient of
is:
where we have omitted the value for
, because, in the other summation, the index,
s, goes from 0 to
n; consequently, for
, there is not the value
.
By using Equations (11) and (14), the coefficient of
becomes:
while the coefficient of
becomes:
By calculating these expressions in
, we find
and
, respectively, from which Equations (54) and (55) easily follow.
With the same passages, we find the other properties, that is, λ appears in only for and appears in only for or .
More precisely, is the coefficient of in , calculated in ; but, λ appears in Equation (56) only for , so that the other summation may assume only the value, , giving only the coefficient of and not of the others, . Similarly, appears in Equation (56) only for , so that the other summation may assume only the value and , giving only the coefficients of and and not of the others, .
By using Equations (54) and (55), we can rewrite Equation (53) as:
which, by using Equations (50) and (52) becomes:
This equation may be easily integrated; in fact, there exists surely a function,
ψ, such that:
After that, Equation (58) says that:
where
is an arbitrary function that does not depend on
and arises from the integration with respect to
.
4.2. On the Eventual Condition
This symmetry condition (2)2 has to be imposed only if there is a conservation of angular momentum, differently from Equation (2)1, which holds surely, for physical reasons.
By using Equation (8), it becomes:
Now, thanks to Equations (54) and (55), we have:
which, by using Equation (57) becomes:
It follows that:
where in the last passage, Equation (52) has been used.
Consequently, the condition (61) becomes:
This equation, together with Equation (58), is nothing else than the integrability conditions, which allow one to obtain the scalar function,
ψ, such that:
Vice versa, from Equation (63), the conditions (58) and (62) easily follow. We may conclude that the four-potentials,
and
, are determined, except for an arbitrary scalar function,
ψ, depending on
and
.
Let us now test the results of the present paper, by comparing their restriction to the 14-moments model with those found, for this particular case and in a completely different way, in [
5,
6].
5. The Particular Case with 14-Moments
A comparison of the definition (46)
2 of
with Equation (28)
2 shows that
is equal to
calculated in
. Now, in Equation (16)
2 of [
5], we find the expressions of
(which in that paper, were called
) in the 14-moments case; so, they allow one to obtain:
where, in the present notation,
and
.
Now, let us assume that
depends on
as a composite function of
,
,
and
and of Equation (64); if we substitute it in (4)
1 of [
5], we find that it is identically satisfied, as we expected for what was said above in the present paper.
Similarly, if we assume that:
with
,
and
given by Equation (64), we obtain that (4)
2 of [
5] is identically satisfied, as we expected.
Regarding the expression of
, we know from the representation theorems [
7,
8,
10] that
with
,
and
arbitrary scalar functions. Moreover, from the same representation theorems, we know that
and
are arbitrary functions of
,
,
,
,
,
,
and
.
Now we expect that there is a correspondence between the solution found in the present paper and that found in [
5], with a completely different method. This is true; in fact, in [
5] it has been found that
is an arbitrary function of eight scalar variables
-
. By comparing their expressions (see Equation 6 of [
5]), with the solution here found, we see that the correspondence between the scalars is the following one:
In order to verify these relations, it suffices to substitute Equation (64) in the lefthand sides and the expressions of
-
(reported in Equation 7 of [
5]) in the righthand sides; after that, we have simply to note that identities are obtained, thanks also to the Hamilton-Cayley theorem (see also on page 32 of [
7]).
We note also that Equation (67) is invertible. In fact, the righthand side of Equation (67)i is a linear function in the variable ; moreover, the coefficient of and the term of zero degree in are functions, depending only on , and ⋯ and . This fact proves that Equation (67) is easily invertible.
For what concerns
, let us analyze firstly the result of [
5]. Equation 5 of this paper gives
and
in terms of four arbitrary functions,
,
,
and
, depending on the variables,
-
. By deducing from (5)
2 the function
in terms of
,
,
and
and by substituting it in (5)
1, we find:
where we have used Equation (6)
1 of [
5].
After that, by using also Equations 6 and 7 of [
5] and the Equation (64) of the present paper, we find the following identities:
by substituting these values in Equation (68), we find:
with
From the arbitrariness of
,
and
, that of
,
and
follows, and
vice versa.
We note now that Equation (70) is equivalent to Equations (65) and (66) of the present article. So, we have found a complete proof that as the result of the present work, restricted to the 14 moments case, is the same of that known in the literature and found with a completely different method.