Lifts of Symmetric Tensors: Fluids, Plasma, and Grad Hierarchy

Geometrical and algebraic aspects of the Hamiltonian realizations of the Euler’s fluid and the Vlasov’s plasma are investigated. A purely geometric pathway (involving complete lifts and vertical representatives) is proposed, which establishes a link from particle motion to evolution of the field variables. This pathway is free from Poisson brackets and Hamiltonian functionals. Momentum realizations (sections on T*T*Q) of (both compressible and incompressible) Euler’s fluid and Vlasov’s plasma are derived. Poisson mappings relating the momentum realizations with the usual field equations are constructed as duals of injective Lie algebra homomorphisms. The geometric pathway is then used to construct the evolution equations for 10-moments kinetic theory. This way the entire Grad hierarchy (including entropic fields) can be constructed in a purely geometric way. This geometric way is an alternative to the usual Hamiltonian approach to mechanics based on Poisson brackets.


Introduction
Mechanics can be typically cast into a Hamiltonian form; the state variables can be for instance particle position and momentum, rigid body angular momentum [1], distribution functions in kinetic theory [2,3], hydrodynamic fields [4,5], electromagnetic fields [6], etc. An advantage of the geometric formulation is that it provides additional leads towards proper coupling of the particular theories, e.g., MHD [7], as well as automatic consistency with mechanics. Such properties can be used also in the numerical simulations [8][9][10][11][12]. As a manifestation of the skew symmetry of the Poisson bracket, the Hamiltonian function H is conserved throughout the motion. That means that Hamiltonian systems are energy-preserving, which manifests the reversible character of the Hamiltonian dynamics, see also [13] for discussion of time-reversal symmetry and Onsager-Casimir reciprocal relations. It is possible to generalize the Hamiltonian framework in a way that is proper also for irreversible systems.
{F, H} g * (µ) = ± µ, δF δµ , δH δµ g , where F and H are two function(al)s defined on the dual space g * , and the pairing at the right hand side is the one between g * and g. The bracket inside the pairing in Equation (2) is the Lie algebra bracket on g, and the notation δH/δµ stands for the Fréchet derivative of the functional H. In this case, the dynamics is governed by the Lie-Poisson equationṡ where ad * denotes the coadjoint representation of g on g * , obtained by the dualization of the adjoint action of g on itself. Here, the adjoint action is defined to be the Lie algebra bracket on g. Diffeomorphism Groups. The configuration spaces of some of the physical systems, such as fluids and plasma theories, are diffeomorphism groups [31,[35][36][37]. Diffeomorphism groups are infinite dimensional Lie groups [38][39][40]. In these cases, the Lie-Poisson formulation, presented in the previous paragraph, takes the following particular form. Assume that a continuum rests in a manifold M in R 3 , then diffeomorphism group Di f f (M) acts from left on the particle space M by evaluation, whereas right action commutes with the particle motion and constitutes an infinite dimensional symmetry group of the kinematic description. This is called the particle relabeling symmetry [41]. We assume the Lie algebra of Di f f (M) as the space X (M) of smooth vector fields. Here, the Lie algebra bracket is the minus of the Jacobi-Lie bracket of vector fields, that is, where L X denotes the Lie derivative operator. We define the dual space X * (M) of the Lie algebra as the space of one-form densities Λ 1 (M) ⊗ Den (M) on M. Here, the pairing between a vector field X and a dual element Π ⊗ dµ is defined as where the pairing inside the integral is the canonical one between the covector Π(z) and the vector X(z). Here, dµ is a density, that is, a volume form on M. The adjoint action of the Lie algebra onto itself is defined by the Lie algebra bracket in (4). A simple calculation shows that the coadjoint action of the Lie algebra X (M) on its dual X * (M) is where div dµ X denotes the divergence of the vector field X with respect to the volume form dµ. At this point, without lost of generalization, we fix the volume form dµ, so that we particularly consider a dual element as a one-form Π. For this choice, a Hamiltonain functional H generates the Lie-Poisson equations on the dual space X * (M), given bẏ In the divergence-free case (corresponding to volume-preserving motion), where the second term on the right hand side of Equation (8) identically vanishes, we obtain equationṡ Π = −L X Π. (9) A geometric pathway to the dynamics of the continuum. To describe the motion of a continuum, one may start to write down the whole microscopic data, involving the interactions, which is very difficult. The kinetic theory uses statistical concepts to handle practical problems of microscopic theory. In previous publications [42,43], a geometric pathway has been proposed for the incompressible Lie-Poisson Equations (9). Even though we present this geometry in the upcoming section, let us briefly summarize this geometrization. In order to arrive at the Lie-Poisson dynamics, start with the most basic ingredient of the theory, a vector field X generating the motion of a single particle, then, by applying purely geometric operations such as complete cotangent lifts and vertical representatives, define a generalized vector field VX c * on the space of sections. This results with the Lie-Poisson dynamics presented in (9). Explicitly, it has been proved in the works by the authors of [42,43] that the Lie-Poisson equations can be written in the form oḟ where VX c * is the vertical (evolutionary), representative of the complete cotangent lift X c * of the vector field X generating particle motion. Here,Π v is a vector field defined as the vertical lift of the one-forṁ Π, see, for example, the work by the authors of [44]. Note that, the geometrization procedure in Equation (10) for the Lie-Poisson Equations (9) is free of a Hamiltonian functional, and thus free of a Poisson structure. As we shall argue in the main body of the paper we believe that such a geometrization looks promising for continuum theories. Let us make a mathematical remark here to depict this geometrization in more technical terms. A single particle traces a curve in the configuration space, so that its motion is determined by an ordinary differential equation (ODE) with time as the independent variable. On the other hand, the motion of the whole continuum is determined by a partial differential equation (PDE) governing a (scalar or vectorial) field. That is, mathematically, the single particle motion is a submanifold of a tangent bundle, whereas the motion of the continuum is a submanifold a jet bundle. Therefore, in order to find a link between these two motions, one needs to propose a geometric operator taking a tangent vector to a jet bundle element. Further, the final product of this operator must forget the motion on the base level in order just to concentrate on the field parameters. The geometrization in Equation (10) does both of these two tasks simultaneously, by relating the base motion to the motion of the field Π as well as by removing affects of the dynamics on the base level.
The goal of the present paper. In the literature, the geometrization (10) has only been studied for volume preserving dynamics [42,43], and no relationship with compressible dynamics has been established so far. Our aim in the present work is to extend the application area of the geometric pathway (10) to compressible systems and establish purely geometrical relations between all possible realizations (momentum or the usual field equations). The geometrization (10) depends on the duality between vector and covector fields, our first novel result is to put the averaged 2D-Euler equation into this framework by redefining the duality by a Sobolev norm. For further generalization, we define an injective Lie algebra homomorphism called Generalized Complete Cotangent Lift (GCCL). GCCL is a mapping from the space of symmetric tensors on a manifold to Hamiltonian vector fields on the cotangent bundle of the manifold. The dual of GCCL is a Poisson mapping relating the Vlasov dynamics in the momentum formulation to (both compressible and incompressible) Euler's fluids in the momentum formulation. Then, using GCCL, compressible fluid motion is put into the framework of (10). As a byproduct, a unique decomposition of the space of symmetric tensors is established as a matched pair Lie algebra. We further equip the compressible fluid flow with entropy. Then, we study 10-moment hierarchy, which paves the way towards to whole Grad hierarchy including entropic moments. In summary, we provide an alternative purely geometric construction of mechanics in kinetic theory and Grad hierarchy.
Let us sum up the main novelty of this work. The Lie-Poisson dynamics is determined by the Lie-Poisson bracket and a Hamiltonian function(al), see, for example, the work by the authors of [30]. In the present work, we first establish that (see Section (2.3)) one can write the Lie-Poisson dynamics for incompressible fluid and Vlasov plasma without referring to a Poisson bracket or a Hamiltonian function(al). Inspired by a series of papers on the moments of the Vlasov dynamics, see, for example, the works by the authors of [45][46][47][48][49], we propose a Lie algebra homomorphism GCCL from the symmetric contravariant tensors on a manifold to the Hamiltonian vector fields on its cotangent bundle (see Section 3.2). Dualization of this mapping determines the moments of the Lie-Poisson equations of form (9) (see Section 3.5). Double application of the GCCL operator and geometric pathway (10) lead to the compressible fluid dynamics (see Section 4.1) and the 10-moment approximation (see Section 4.2) free from a Poisson bracket and a Hamiltonian function(al).
The contents. This paper comprises three main sections. In the following section, we present the geometrical and algebraic foundations of the geometrization (10) as well as comment on underlying physical intuitions. Examples of this geometry are incompressible fluid flow and Vlasov's plasma. We close this section by redefining the duality with a Sobolev metric and writing an averaged 2D-Euler's equation in the form of (10). In the third section, Lie algebra of symmetric contravariant tensor fields and a matched pair decomposition of this space is established. An injective Lie algebra homomorphism, which we call Generalized Complete Cotangent Lift (GCCL), is defined. It is shown that the dual of GCCL is a Poisson mapping from the plasma level to the fluid level both in momentum formulations. In the fourth section, we generalize the geometrization (10) to compressible fluid motion including entropy density. Finally, the procedure is generalized to also cover the 10-moment Grad hierarchy in kinetic theory and incorporate a complimentary hierarchy of kinetic moments.
Notation. Let M and Q be finite dimensional manifolds equipped with the local coordinates (x a ) and (q i ), respectively. Cotangent bundles T * M and T * Q carry canonical symplectic two-forms, in Darboux' coordinates given by Ω M = dx a ∧ dΠ a and Ω Q = dq i ∧ dp i , respectively. The Hamilton's equation is generated by a Hamiltonian function H on T * Q, and it is defined as where ι is the contraction operator. Further, it is assumed that all the requirements of functional analytic issues, such as existence, uniqueness, regularity, and convergence, are satisfied.

Complete Tangent and Cotangent Lifts
Let M be an m−dimensional manifold equipped with local coordinates (x a ). A vector field X on M generates a flow on the manifold M, say φ t : M → M, which describes motion on the manifold with respect to a parameter t (representing the time). Locally, this reads that the vector field X(x) = X a (x)∂/∂x a generates infinitesimal transformations:

Complete Tangent Lift
Let us first recall the concept of complete cotangent lift, see, e.g., the work by the authors of [50], including its physical interpretation for clarity. Attaching a tangent plane to each point of the manifold, i.e., constructing the tangent bundle TM with induced local coordinates (x a , v a ), the vector field X tells which vector of the tangent plane at each point describes the motion. In other words, the vector field tells the direction and the velocity with which the motion continues from the point, and it can be interpreted as a mapping from M to TM. To see this, consider the flow φ c t on the tangent bundle TM defined by the following equation.
where τ M is the tangent bundle projection mapping a vector in TM to its initial point in M. Flow φ c t constitutes a one-parameter group of diffeomorphisms on TM called the complete tangent lift of the flow. From the differentiation of Equation (13) with respect to t at t = 0 we define the complete tangent lift X c of X as follows Notice that X c is a vector field on TM. Consider now a vector v ∈ T x M, which by definition corresponds to a curve γ(t) passing through x satisfyingγ = v. The transformation (12) maps the curve to a new curvex(γ(t)). Therefore, the transformation maps the vector v to a new tangent vector with components From the perspective of the tangent bundle TM, a vector field X induces both motion in the ∂/∂x a direction as well as motion in the ∂/∂v b direction, that is,

Complete Cotangent Lift
Consider the cotangent bundle T * M equipped with conjugate coordinates (x a , Π a ). The complete cotangent lift of a flow ϕ t (of a vector field X) on M is a one-parameter group of diffeomorphisms ϕ c * t on T * M satisfying where π M is the natural projection defined on T * M to M. The vector field X c * on T * M, which has the flow ϕ c * t , is called the complete cotangent lift of X [44]. The infinitesimal version of the Equation (17) determines X c * as follows, The complete cotangent lift X c * is dual to the complete tangent lift X c in the following sense. Taking the duality between elements of the tangent planes and elements of the adjoint cotangent spaces, the duality between the cotangent and tangent lifts must be the same, that is, Therefore, we put which means thatΠ a = Π a − tΠ b ∂X b /∂x a . From the perspective of the cotangent bundle T * M, a vector field X also induces motion in the momentum, which is the complete cotangent lift of vector field to X(T * Q), see, e.g., the work by the authors of [50]. The complete cotangent lift of X expresses how both position on the manifold Q and covectors (one-forms) vary along the motion induced by the field; a diagram summarizing the discussions done so far follows.

Lemma 1. Maps given in
are Lie algebra isomorphism intos, that is, for all X, Y ∈ X (M) .

From Jet Bundle to Tangent Bundle
Triple (E , π, M) denotes a smooth bundle with coordinates (x a ) on the base manifold M and x a , u λ on the total manifold E . J 1 π is the first order jet manifold associated with (E , π, M) with induced coordinates x a , u λ , u λ a . There exist fibrations, of J 1 π on E and M, respectively, [52]. Contact forms determine holonomic sections of the jet bundle fibration.

Holonomic Lifts of Vector Fields
Consider a fiber bundle (E , π, M). Let X be a vector field on M and consider a section σ of the fibration π. The lie derivative (directional derivative) of a smooth function, F, defined on the total space, E , with respect to the vector field X can be computed by means of σ as follows, L X (F • σ). This definition reads the definition of the holonomic lift X hol of the vector field X by the identity Note that in terms of the local coordinates, for X = X a ∂/∂x a , the holonomic lift defined in (25) is computed to be Note that X hol is not a classical vector field on E , as its coefficients depend on the first order jet bundle. Such kind of sections are called generalized vector fields [52][53][54][55]. In order to justify the term holonomic, we see that the values of X hol at the contact one forms ϑ a in Equation (24) vanish identically. For a vector field Y = Y a ∂/∂x a + Y α ∂/∂u α , we define the holonomic part HY of Y as

Lie Algebra of Generalized Vector Fields
Assume a fibration (E , π, M). In general, a generalized vector field on E takes the form of The first order prolongation pr 1 ξ of ξ is defined by where D x a is the total derivative operator with respect to x a , and u λ ba is an element of the second order jet bundle. The Lie bracket of two first order generalized vector fields ξ and η is the unique first order generalized vector field If ξ and η are two classical vector fields on E , then [•, •] pro reduces to the Jacobi-Lie bracket of vector fields [54], which results in the following lemma.

Lemma 2.
The mapping H : Y → HY in (27), taking a vector field to its holonomic part, is a Lie algebra isomorphism from the space of projectable vector fields into the space of generalized vector fields of order one.

Vertical Representatives
Note that the holonomic lift of a vector field reads the action of the vector field X on the fiber coordinates. So that, for a projectable vector field Y on E , to read just the vertical motion, that is, the dynamics governing the sections, one needs to subtract holonomic part inside Y, so that the vertical motion is the one obtained by subtracting the holonomic part HY of the vector field from itself [54]. We call this the vertical representative: Note that VY lies in the kernel of Tπ. On the other hand, the generalized bracket of vertical representatives satisfies B is a vertical-vector valued two-form where the brackets are as in (30), see [42]. We require that generalized vector fields are projectable [53].

Lie-Poisson Dynamics of Incompressible Systems
Now, we concentrate on the case of the cotangent lift and its vertical representative. For this, we consider the vector bundle (T * M, π M , M), and let Π be a section of this bundle. Let X be a vector field on X, then its holonomic part is The vertical representative of the cotangent lift X c * , that is, We have then the mapping taking a vector field, X, on Q to a vertical vector field on T * Q, by first taking the cotangent lift, then taking the vertical representative. For the case of complete lifts, the vector-valued two-form B defined in Equation (33) vanishes. This reads as the following lemma.

Lemma 3.
The mapping V c * in (37) is a Lie algebra isomorphism into, that is, Vertical representative is sought for a vector field X c * on cotangent bundle T * T * Q. The cotangent bundle has a natural fiber structure T * T * Q over the base T * Q. The vector field thus has a component in the direction (∂/∂q i , ∂/∂p i ) as well as in (∂ Π i , ∂ Π i ). For a given point (q i , p i ), the component of the vector field in the (∂ Π i , ∂ Π i ) directions describes evolution of the momentum variables (Π i , Π i ), which can be seen as motion along the fiber attached to (q i , p i ). The momentum coordinate can be seen as a function of the point (q i , p i ), and the time-evolution equation for (Π i , Π i ) is the vertical component of the vector field. However, as also the point (q i , p i ) is subject to motion, evolution of (q i , p i ) affects the momentum coordinate (Π i , Π i ). By subtracting this evolution, X hol h , from the vertical vector field, X c * h , we obtain the vertical representative of the total vector field, which determines the evolution of the momentum variableΠ This is a purely geometric way for extending dynamics on a manifold to the cotangent bundle on the manifold.

Vertical Lifts of One-Forms
Consider the cotangent lift T * π M : T * M → T * T * M of the projection π M : T * M → M, and recall the isomorphism Ω T * M : T * T * M → TT * M associated with the symplectic two-form Ω T * M on the cotangent bundle T * M. Euler vector field is defined as which is a vertical-valued vector field, that is, of the one-from α is a vertical-valued vector field on T * M [44]. Taking the coordinates (x a , y b ) on T * M, the Euler vector field is computed as X E = −y a ∂/∂y a and the vertical lift of the one-form

Geometry of Lie-Poisson Equations
Now, we will collect all the geometric structures introduced in this section to arrive at the Lie-Poisson equations in terms of lifts and vertical representatives. Starting with a differential one-form Π = Π a (x)dx a , an element of the space Λ 1 (M) of one-form sections on M. Assuming that Λ 1 (M) is a vector space, its tangent space equals to the product TΛ 1 (M) = Λ 1 (M) × Λ 1 (M). Let us consider a curve Π(t) in Λ 1 (M). The time derivativeΠ is an element of the tangent space. As we have identified T Π Λ 1 (M) = Λ 1 (M), the geometric objectΠ is a differential one-formΠ =Π a (x)dx a which is in Λ 1 (M). Now, we compute the vertical lift of this one-form section as explained in (40), this locally Therefore, a direct observation results in the following proposition.

Proposition 1.
If the motion of a single particle is governed by a volume preserving vector field X on a manifold M, then the Lie-Poisson equation governing the motion of the continuum consisting of such particles can be written asΠ Proof. In order to prove this observation, recall the Lie-Poisson equations given in Equation (9). See that minus of the Lie derivative on the right hand side equals to VX c * that is the vertical representative of the complete cotangent lift of X. By employing Equation (41), one immediately arrives at the required result (42).
We remark that in this geometrization, one of the crucial step is to determine the dual space. Here, we are showing this fact explicitly in the following example.

Example: Incompressible Fluid Flow
For an ideal incompressible fluid in a bounded compact region, Q ⊂ R 3 , the configuration space is the group Di f f vol (Q) of volume preserving diffeomorphisms on Q. The Lie algebra X div (Q) of Di f f vol (Q) is the algebra of divergence-free vector fields parallel to the boundary of Q, and the dual space of one-form modulo exact one-form densities on Q. Here, [Υ] = {Υ + dp :p ∈ F (Q)} ∈ Λ 1 (Q)/dF (Q) denotes the equivalence class containing Υ, and the volume three-form d 3 q is the Euclidean volume on R 3 [5,56]. Let (x a , Υ b ) be induced coordinates and X = X a ∂/∂x a be a divergence-free vector field. Then, according to Equation (36), and the equations of motion for the dynamics generated by VX c * are For a generic element Υ + dp ∈ [Υ] , Equation (45) becomes Euler's equations for ideal fluid, that (45) becomes the Euler's equation in vorticity form ∂ω/∂t + L X ω = 0.

Example: Vlasov's Plasma
Configuration space of collisionless and non-relativistic plasma motion is group Di f f can (T * Q) of canonical diffeomorphisms on the phase space T * Q of configuration manifold Q ⊂ R 3 of individual charged particles [2,57,58].

Lie Algebra of the Canonical Diffeomorphisms
Lie algebra of the group is space X ham (T * Q) of Hamiltonian vector fields on T * Q equipped with minus of the Jacobi-Lie bracket [•, •]. We can identify the space X ham (T * Q) with the space of nonconstant smooth functions (more terminologically nonconstant Hamiltonian functions) on T * Q equipped with the canonical Poisson bracket {•, •} as the Lie algebra bracket, that is, This is a manifestation of the identity

The Dual Space
The identification (47) reads that Lie algebra of the canonical diffeomorphism can also be considered as F (T * Q)/R. In this case, the dual space of the Lie algebra is the space Den(T * Q) of densities on T * Q. By fixing the symplectic volume Ω 3 Q , we can further identify the dual space with the smooth functions F (T * Q). This fits the classical approach, where elements in the function space are accommodated as the plasma density functions. We are now returning to the very first definition of the Lie algebra consisting of the Hamiltonian vector fields, and try to define a dual to that space consisting of the differential one-forms.

Lemma 4.
The following identity holds, where Ω T * Q : Π → Π is induced from the symplectic two-form Ω T * Q .
Proof. With this definition of the dual space the L 2 -pairing of the Lie algebra and its dual becomes nondegenerate provided we take the symplectic volume dµ = Ω 3 where we have applied integration by parts in the second line ( [30], internet supplement). The calculation can be also carried out in the Darboux coordinates using X h = L · dh, where L is the Poisson bivector (inverse of Ω Q ).
Proof. Then, the dual of the Lie algebra isomorphism h → X h is and it is a momentum map. Notice that the operator div Ω Q takes the one-form of a real valued function.
defines the plasma density function. This calculation leads us to add a subscript f to the notation of the one-form section Π f , so that we have div Ω Q Π f = f .

Momentum-Vlasov Equations [43]
We start with the total energy function h = p 2 /2m + eφ Π of a single particle. Here, φ Π is the potential. Locally, the Hamiltonian vector field for the Hamiltonian function h is computed to be In Darboux' coordinates q i , p i ,Π i , Π i on T * T * Q, the complete cotangent lift of X h reads The vertical representative of the cotangent lift X c * h , is a vertical valued generalized vector field of order 1. Vertical lift of the one-from [44]. Hence, m-Vlasov Equations (56) can be recast in the formΠ that is given explicitly byΠ

Vlasov Equation.
Let us recall the identity (42) and apply it in the present case. The momentum-Vlasov equations can be compactly written aṡ Lemma 5. The operator div Ω Q and the Lie derivative L X h commute for Hamiltonian vector fields X h , that is Proof. Recall the calculation in (49) and apply the present case as follows On the other hand, we have an integration that by part reads, Comparing the first and the second calculations in the proof for an arbitrary function k the proof is completed.
When the dual mapping in Equation (51) is employed, the momentum-Vlasov equations then turn to div If, in particular, the dualization is determined by the L 2 pairing on the function space, we arrive at This reads the Eulerian dynamics in density variables, that is, Vlasov's equation, If h is the total energy of a single particle, we have In the work by the authors of [57], the accompanying Poisson equation has been obtained by a momentum mapping coming from the gauge symmetry of the Hamiltonian dynamics.

Example: Averaged 2D-Euler Equation
In the previous example we have employed L 2 pairing of the functions, that is simply multiply-and-integrate. This determines the structure of the Lie-Poisson equation. Let us consider that Q equals to R with coordinates x, so that the cotangent bundle turns out to be T * Q = R 2 with coordinates (x, y). Now, we change the pairing to the Sobolev H 1 -pairing, given by where λ is a real parameter.
Here, ∇ f is the gradient of f . After applying integration by-parts to the second term and omitting the total divergence terms we write H 1 -pairing in terms of the L 2 -pairing as follows In this framework, the momentum map defined in (51) takes the form of To see that, consider the following equalities where we have employed Lemma (4) for the case of T * Q = R 2 and used the Sobelev metric (67) on the function space. Note that in this case we omit the subscript in order not to mix this mapping with the one in (51). Let us apply the momentum mapping (68) to both sides of the equationΠ = VX c * h . From the left hand side, one arrives at for the left hand side, one computes where we have used (36) in the first equality and used the commutation relation (5) in the second equality. Therefore, we have Assume, in particular, that f = Ω is the vorticity of an ideal inviscid incompressible homogeneous fluid and h = Ψ is the stream function determined by the equation Ω = ∆Ψ; then, this system turns out to be an averaged 2D-Euler equation, [59,60] (

Generalized Complete Cotangent Lift
Before starting to elaborate the fluid theories and the kinetic moments of the plasma dynamics, we study some geometrical arguments motivating from [47].

Schouten Concomitant
Direct product TQ = ⊕ ∞ n=0 T n Q of spaces T n Q of symmetric contravariant tensor fields on a tmanifold Q ⊂ R 3 of all orders constitutes an infinite dimensional vector space. In a local coordinate system q i on Q, an element of TQ can be written in the form of where X n ∈ T n Q is a symmetric contravariant tensor field of order n and X i 1 i 2 ...i n are real valued coefficient functions. T 0 Q is the space F (Q) of smooth functions and T 1 Q is the space X (Q) of vector fields on Q. Schouten concomitant is a Lie algebra structure on the space TQ [39,49,61]. Here X n , Y m and Z n+m−1 are contravariant tensor fields of orders n, m and n + m − 1, respectively. The coefficient functions of Z n+m−1 in terms of those X n and Y m are

Lie Subalgebras of Schouten Algebra
For the zeroth-order tensors, that is, for the space of smooth functions F (Q), Schouten concomitant reduces to the trivial Poisson bracket of functions on Q. So that F (Q) is a subalgebra of T(Q). For the first order tensors, that is, for the space of smooth vectors X (Q), the concomitant turns out to be the Jacobi-Lie bracket of vector fields. For instance, when X 1 and X 2 are vector fields in the classical sense, then the coefficient function becomes which is the Jacobi-Lie bracket (commutator) of the two vector fields. So, X (Q) is another subalgebra of T(Q). A semidirect product of the subalgebras is another subalgebra of T(Q). In this case, the Schouten concomitant becomes where L X ρ is the directional derivative of the function ρ in the direction of X and [X, Y] is the Jacobi-Lie bracket of vector fields X and Y.
Note also that the complement of the vector space s in TQ is closed under the Schouten concomitant as well. We denote this subalgebra by To see this, herein, we only recall the graded character of the Schouten bracket. The least order tensors in n are of order 2, and the Schouten bracket of two such tensors is a third-order tensor, so it is in n. We record here the decomposition TQ ∼ = s ⊕ n. That is, we can write any generalized tensor field X in the Lie algebra TQ in the form of where the first factor (ρ, X) lives in the subalgebra s whereas the second factor X = ∑ ∞ n=2 X n is in the Lie subalgebra n.

Remark 1.
This splitting can be interpreted in the context of Grad hierarchy of kinetic theory. The s subalgebra represents fluid mechanics (compressible and isentropic), whereas subalgebra n represents the higher moments of the Grad hierarchy. The splitting then means that closed evolution equations can be formulated within s, within n or within the whole TQ. The first case corresponds to fluid mechanics (Euler equations), the second to dynamics of higher moments only (the reducing dynamics approaching fluid mechanics), and the third case is equivalent to solving the whole Vlasov (or Boltzmann) equation, see the work by the authors of [62] for details.

Matched Pair (Bicross Product) Realization of TQ
Let us examine algebraic foundation of the decomposition of TQ ∼ = s ⊕ n presented in the previous paragraph. To this end, we find the possible actions of s and n onto each other by computing the Schouten bracket where X k is a symmetric tensor field of order k greater then or equal to 2. Here, L X X k is the Lie derivative of X k in the direction of X, whereas Therefore, if k is strictly greater than 2, then the order of the tensor [X k , ρ] SC is strictly greater than 1 and it is an element of n. Otherwise, that is for k = 2, [X k , ρ] SC is a first order tensor field (that is a vector field in the classical sense), so that [X 2 , ρ] SC is in s. If the algebraic structure of the decomposition TQ ∼ = s ⊕ n were a direct product, then two components (ρ, X) and X = ∑ ∞ k=2 X k , c.f. Equation (80), would be orthogonal with respect to Schouten bracket. However, we have proved that this is not the case for TQ. If the algebraic structure of the decomposition TQ ∼ = s ⊕ n were a semidirect product, then the Schouten bracket of an element of n and an element of s would lie in one of the subspaces. But we have shown, this is not the case for TQ either. Instead, the bracket [X, (ρ, X)] SC results in some terms lying in s and some other terms lying in n. In light of the work by the authors of [63], this observation manifests that the total space TQ can be realized a matched pair of its subalgebras s and n. A matched pair is a generalization of the semidirect product in the sense that there exists mutual nontrivial actions of s and n on to the each other. If there is only a one-sided action of s on n or vice versa, the matched pair becomes a semidirect product. The mutual actions are computed from the identity [X, (ρ, X)] SC = X (ρ, X) ⊕ X (ρ, X).
Here, the first term X (σ, X) is the left action of X ∈ n on (σ, X) ∈ s, and the second term X (σ, X) is the right action of (σ, X) on X. Accordingly, we compute the mutual actions as follows, : s ⊗ n → n, We summarize this matched pair decomposition of TQ in the following proposition which says that TQ is a matched pair Lie algebra. Proposition 3. The space TQ of symmetric contravariant tensor fields can be written as a matched pair product of its Lie subalgebras s and n exhibited in (77) and (79), that is, where the mutual actions are in Equations (83) and (84).

Remark 2.
Kolmogorov cascade. Simple fluids are fluids with an internal structure that remains unchanged during the time evolution. Experience shows that for such fluids the level of classical fluid mechanics (i.e., the level on which the hydrodynamic fields play the role of state variables) is autonomous. This experimental observation is compatible with the result (Equation (85)). The question now is whether there are autonomous mesoscopic levels with a larger, but finite, number of Grad fields. The result (Equation (85)) indicates that the answer is negative. There are two additional observations supporting the negative answer and thus indirectly also (Equation (85)). First, it is the dissipation added by Boltzmann to the Hamiltonian kinetic equation. The Boltzmann dissipation, which indeed drives solutions to the level of fluid mechanics, is supported by an independent and a very strong physical argument, namely, that the principal culprit of the disorder created in the gas particle trajectories are the binary collisions. There does not seem to be any other physical process in the gas that would drive solutions to a higher order mesoscopic level.
The second is the observation of turbulent flows. When the external force driving the laminar flow increases the macroscopic order of the laminar flow starts to break up. Vortices start to emerge. The vortices can be regarded as an internal structure characterized by higher order Grad moments (i.e., by the fields that have the physical interpretation of higher order velocity correlations). Fluids subjected to a turbulent flow can be thus regarded as complex fluids (i.e., fluids in which the time evolution of the internal structure is coupled to the time evolution of the hydrodynamic fields). Observations of turbulent flows show that the break up continues into smaller and smaller vortices until they completely disappear and become a part of the molecular motion. This observation is known as Kolmogorov cascade. If there was an autonomous mesoscopic level with n higher order Grad moments then, when the vortices would reach the size corresponding to the n-moment, the turbulent flow would become a flow that, in the context of the n -Grad moment level, would appear laminar. We would then expect that the continuation of the break up would appear as an onset of an n-order turbulence, i.e., a turbulence emerging on the n-Grad moment level. In other words, the Kolmogorov cascade would have a more complex dependence on the driving force than the one observed.
Due to the mutual interactions existing in matched pair products, they can be considered as a generalization of the semidirect products. The Lagrangian and Hamiltonian dynamics on these systems are available [64,65], and the discrete matched pair dynamics of Lie groupoids are discussed in the work by the authors of [66].

Lie Group Underlying s
Let us denote the right and left actions of group Di f f (Q) of diffeomorphisms on Q by respectively. Infinitesimal generators of the diffeomorphisms ψ R,L g are vector fields X R,L on Q, respectively. By fixing a point q ∈ Q in (86), we obtain induced mappings ψ R,L q from the group Di f f (Q) to the manifold Q. Actions of Di f f (Q) on space F (Q) of smooth functions are given by means of pull back operation, that is, Infinitesimal generators X R,L F (Q) of the transformations σ R,L g are vector fields on F (Q), that is linear transformations on F (Q). Explicitly, we compute the generators as follows, where Tφ is the tangent mapping of φ. L X R,L denote the Lie derivative operator and in this case, they are directional derivatives of φ in the directions of X R,L (q), respectively. Actions σ R,L of Di f f (Q) on F (Q) in Equation (87) define respective semidirect product structures at (g, φ) are given by two-tuples (X R,L g , ν). Here, X R,L g are material velocity fields satisfying the identities τ Q • X R,L g = ψ R,L g , where τ Q is the tangent bundle projection and ψ R,L g are diffeomorphisms in Equation (86). The second term ν is a function on Q, that is an element of F (Q). Tangent space T (e,0) S R,L at the identity (e, 0) is the product space X (Q) × F (Q). It is the underlying vector space for both of the Lie algebras s R,L induced from the group structures S R,L , respectively. The Lie algebra bracket on s R is the subalgebra structure given in (78) whereas the bracket for s L needs a minus sign.

Tensors to Functions
Let X k be a contravariant tensor field of order k. Due to the canonical inclusion T k Q → T k (T * Q), we may assume X k as a tensor field on the cotangent bundle T * Q. Using the canonical one-form θ = p i dq i , we define a mapping from T k Q to the space F (T * Q) of smooth functions on T * Q by contracting the contravariant tensor X n ∈ T n Q with n-th tensor power θ n T * Q = θ T * Q ⊗ ... ⊗ θ T * Q of the canonical one form θ T * Q , that is, We extend the operation given in Equation (91) to the product space TQ. For X = ⊕ X n ∈ TQ we define a function h X on T * Q as the sum Ref. [67]. This infinite sum may be considered as the Taylor expansion of the function h X in terms of the p−polynomials. A straightforward calculation proves the following lemma. For this result, we refer to the work by the authors of [47].

Lemma 6.
The map X → h X is a Lie algebra anti-homomorphism, that is, where the bracket at the left hand side is the Schouten concomitant of contravariant tensors whereas the bracket at the right hand side is the canonical Poisson bracket of functions on T * Q.

Generalized Complete Cotangent Lift (GCCL)
We further define the following operation, which we call the generalized complete cotangent lift (and denoted as GCCL) from X k , to the space of Hamiltonian vector fields which take a contravariant tensor field X n on Q to Hamiltonian vector field, corresponding to the Hamiltonian function defined by Equation (91), c.f., the work by the authors of [68]. In Darboux' coordinates q i , p i , the Hamiltonian function h X k is a p−polynomial, and GCCL of X k is thus Notice that GCCL of a vector field X 1 = X is exactly the same as the complete cotangent lift of (21). GCCL is indeed a generalization of the complete cotangent lift. More generally, the generalized complete cotangent lift of X = ⊕ X n is defined as Using the identity in (47) and Lemma (6), we arrive at the following equalities, which enable us to state the following proposition.

GCCL on the Subalgebras of T(Q)
We have presented four subalgebras of the the Lie algebra T(Q) equipped with Schouten bracket. They are the space of smooth functions F (Q), the space X(Q) of vector fields, their semidirect product s in (77) , and the complement of the space s denoted by n in (79). On F (Q), the mapping in Equation (91) reduces to the natural inclusion of F (Q) into the space F (T * Q), h ρ = ρ, see also [30]. GCCL then takes the particular form Note that the Jacobi-Lie algebra bracket of vector fields of the form X ρ on T * Q vanishes. In the following section, we show that X ρ generates the gauge invariance of the canonical Hamiltonian structures on T * Q. For X(Q), GCCL reduces to the complete cotangent lift X c * given in (21). Moreover, GCCL X c * is the infinitesimal generator of the right action of the diffeomorphism group Di f f (Q) on T * Q. Image of (ρ, X) in s under the GCCL is which is the sum of the vector fields X c * and X ρ in Equation (97). Notice that GCCL(ρ, X) is a Hamiltonian vector field with the Hamiltonian function, see the work by the authors of [30], This result can be seen for instance in the work by the authors of [48]. Moreover, the Lie algebra identity in Equation (96) gives that where the bracket on the left hand side is minus the Jacobi-Lie bracket of vector fields on T * Q, whereas the bracket on the right hand side is the semidirect product structure on s given in Equation (78).
The space s is a subalgebra of g. The product vector field [(ρ 1 , X 1 ) , (ρ 2 , X 2 )] s is a Hamiltonian vector field on T * Q with Hamiltonian function We can regard Equation (99) as a mapping, and, due to the identity in Equation (101), it is an embedding of the algebra s into F (T * Q). In other words, F s (T * Q) is an isomorphic copy of s in the space F (T * Q).

The Dual of TQ and Kuperschmidt-Manin Bracket
The dual T * Q of TQ is the direct sum ⊕ ∞ n=0 T * n Q of symmetric covariant tensor fields T * n Q of all order [47,49]. In coordinates q i , an element of T * Q is given by where A n ∈ T n Q is a symmetric covariant tensor field of order n. The pairing between T * Q and TQ is given by the infinite sum where d 3 q is a volume form on Q. As the dual of a Lie algebra, T * Q has a Lie-Poisson structure called the Kuperschmidt-Manin bracket. Note that since the result of differentiation is a multivector field, the Schouten concomitant is needed. Here, F and H are functionals on T * Q, and the reflexivity assumption takes the particular form δF/δA ∈ TQ. The bracket inside the integral is the Schouten concomitant, and the pairing inside the integral defined in relation (103) [45].

The Dual of s and Compressible Fluid Bracket
On the dual space s * = Den (Q) × X * (Q) consisting of densities Den (Q) and one-form densities X * (Q), the Lie-Poisson structure in Equation (2) takes the particular form where CF stands for compressible fluids.
Ref. [69][70][71]. H and F are two functionals on s * and reflexivity condition is assumed, that is, δF/δM ∈ X (Q) and δF/δρ ∈ F (Q). To obtain the equations governing the dynamics of isentropic compressible fluid, we choose the Hamiltonian functional which is the total energy of the continuum consisting of a kinetic term and a potential term with internal energy w = w (ρ).

The Dual of g and Momentum-Vlasov Bracket
The Lie-Poisson bracket on the dual space is given by where the bracket [•, •] inside the integral is the Jacobi-Lie bracket and dµ is the symplectic volume.

Lifts of Actions to Cotangent Bundle
The cotangent lifts of the left and right actions of Di f f (Q) on Q are the right and left actions on the cotangent bundle T * Q, respectively. Both actions, φ R,L g , are canonical, which means they respect canonical symplectic two-form Ω T * Q on T * Q [51]. Thus, for all g ∈ G, transformations φ R,L g are elements of group Di f f can (T * Q) of canonical diffeomorphims on T * Q. Infinitesimal generators of the actions φ R,L g are vector fields on T * Q and computed by d dt where X R,L are vector fields generating ψ R,L g . The mappings φ R,L z are obtained by fixing a point z ∈ T * Q and they are from Di f f (Q) to T * Q. X R,L c * are complete cotangent lifts of X R,L , as described in Equation (21).
We define an action t of additive group F (Q) to cotangent bundle T * Q by momentum translations. Explicitly, action of φ ∈ F (Q) to an element z ∈ T * Q over the point q = π Q (z) is In a Darboux' chart, an element z of T * Q is represented by (q, p) and t φ (q, p) = q, p − ∇ q φ . The canonical symplectic structure Ω Q is invariant under the momentum translations, which is the gauge symmetry of Hamiltonian dynamics. In other words, the transformation t φ is canonical, hence t φ is an element of Di f f can (T * Q). Infinitesimal generator X φ (q, p) = −∇ q φ · ∇ p is a Hamiltonian vector on T * Q with the Hamiltonian function φ = φ (q) regarded as an element of F (T * Q). The mapping F (Q) → X ham (T * Q) : ρ → X ρ is a Lie algebra isomorphism into in Equation (97). The following lemma shows that φ R,L g and t ρ commute up to the actions σ R,L g of Di f f (Q) on ρ [34].

Lemma 7.
The composition of the actions φ R,L g in Equation (108) and t φ in Equation (110) on T * Q are intertwining, that is, where σ R,L g are the actions of Di f f (Q) on F (Q) given in Equation (87).
Proof. Let us consider z ∈ T * Q over the point q = π Q (z). Then, we have This lemma enables us to define two possible embeddings of the semidirect product group S into the group of canonical diffeomorphisms G, given by where the actions t and φ are in Equation (110) and Equation (109), respectively. On the Lie algebra level, this turns out be the mapping which is the one in Equation (98). The infinitesimal version of this lemma is where X L (X R,L ) φ are Hamiltonian vector fields for the functions obtained by the Lie derivations of φ in the directions of X R,L . What we derive in Equation (114) is a particular case of the algebra in Equation (78).

The Dual Mapping of GCCL
In the previous subsection we showed that the GCCL, mapping TQ to g, is a Lie algebra homomorphism, see Equation (96). Therefore, its dual mapping Φ : g * → T * Q is a momentum and a Poisson mapping [30]. Taking explicitly Π f = Π i dq i + Π i dp i ∈ g * , the dual operation becomes where θ n−1 T * Q is the (n − 1)-th tensor power of the canonical one form θ T * Q and ϑ is a one-form on T * Q, given explicitly by The definition stems from the duality Left hand side of this equation is the n-th component of the image of Φ while the right hand side can be explicitly calculated from the definition of GCCL.
The image of Π f under the dual mapping Φ gives the moments of the momentum-Vlasov dynamics. The n−th moment of Π f is given by Note that the substitution of the momentum map Π f → f in Equation (52), we have the kinetic moments of the Vlasov dynamics [45,47]. Indeed, the n-th moment reads explicitly where the Ω-divergence is interpreted as the one-particle distribution function f (q, p) = ∂Π j /∂q j − ∂Π j /∂p j . The A n moment is thus the standard n-th moment in kinetic theory (up to some geometrical prefactors dq). In particular, the zero-th moment reads A 0 = dp∂Π i /∂q i , whereas the first moment is A 1 = dp Π i + p i ∂Π j /∂q j dq i . The following proposition summarizes the situation.
holds for all functionals F and H on T * Q, see also the work by the authors of [48].
Proof. To prove this fact, we consider a linear functional F X on T * Q in form Due to linearity, we have that δF X /δA = X. The pull-back Φ * F X of F X to g * via the momentum map Φ in Equation (115) is where Ω 3 Q is the symplectic volume form on the cotangent bundle T * Q. Variation of Φ * F X with respect to its argument Π f is where X h X is the Hamiltonian vector field corresponding to the Hamiltonian function h X in Equation (92). The momentum-Vlasov bracket is on g * , where the bracket inside the integral is minus the Jacobi-Lie bracket of vector fields satisfying Hence, the Poisson map relation in Equation (118) follows the direct substitutions.

M-Vlasov to Fluid Map
We have established that semidirect product s R is a subalgebra of the space TQ of symmetric contravariant tensor fields. It was also shown that the generalized complete cotangent lift in Equation (93) reduces to injective homomorphism s → s ⊂ g in Equation (98) when restricted to the subalgebra s R . The dual Φ : g * → s * : Π f → (ρ, M) of this Lie algebra homomorphism is the first two moments of momentum-Vlasov dynamics given in Equation (115). In Darboux' chart where Π f = Π i dq i + Π i dp i , the momentum mapping Φ is explicitly given by where ρ is a real valued function on Q and M = M i (q) dq i is a differential one-form on Q, which are the zero-th and first moments A 0 and A 1 , respectively. Hence, we arrive the following lemma. We call the operation in Equation (120) as "m-Vlasov to fluid map". The substitution of g * → F (T * Q) : Π f → f (q, p) gives plasma to fluid map in the work by the authors of [72].

Remark 3 (TQ represents conjugate variables).
The function on the manifold Q discussed above can be thought of as the conjugate density in the energetic representation, i.e., δE δρ with ρ being density of matter.
The function on Q can be thought of as chemical potential, usually denoted by µ. Similarly, the function can also stand for the conjugate entropy density, T = δE δs , which is the field of temperature. The vector field above can be thought of as the conjugate momentum density, v = δE δM , which is the velocity field.

Momentum Formulation of Compressible Fluid Flow
In this section, we shall show how some of the physical theories fit the geometrization procedure presented in Section (2). For this end, we start with a generic Hamiltonian vector field, on the canonical symplectic manifold (T * Q, Ω Q ). Then the complete cotangent lift of X h is a vector field on iterated cotangent bundle T * T * Q, which can be computed in the Darboux' coordinates q i , p i ; Π i , Π i as follows, We use that Π f = Ω Q Π f is the image of a one-form Π f under the musical isomorphism induced by the canonical symplectic two-form Ω Q on T * Q. Π f is given locally by Therefore, the action Π f ∂h/∂q i in (122) is simply the action of the vector field Π f on the real valued function ∂h/∂q i . It is interesting to note that X c * h is a Hamiltonian vector field on the symplectic manifold (T * T * Q, Ω T * Q ) with the Hamiltonian function Π, X h , that is, The decomposition of the complete cotangent lift X c * h into the sum its vertical representative VX c * h and its holonomic part HX c * h are computed to be

Momentum-Euler Equations
It was shown in the previous section that the generalized complete cotangent lift determines an embedding s → s ⊂ g, as given in Equation (98). The image (X, φ) is a Hamiltonian vector field on T * Q. The complete cotangent lift of (X, φ) is the Hamiltonian vector field on T * T * Q satisfying the Hamilton's equations in (124) with the Hamiltonian function The vertical representative of the cotangent lift (X, φ) c * is a generalized vector field of order 1 and is given by the following abbreviated formula, where the coefficient functions are locally in the forṁ We call the system of equations given in (128) the momentum-Euler equations. In the density variable these system of equations reduces to by the substitution of the momentum map in Equation (52) into Equation (128). Note that X i can be though of as the i-th component of the fluid velocity, and φ can be thought of as the chemical potential. Equation (129) can be interpreted physically as dynamics of fluctuations around mean velocity field X i and field of chemical potential φ. Geometrization of the right hand side of Equation (127) can also be described as follows. Vertical lift of the one-form Π f is a vector field: where T * π T * Q is the cotangent lift of the projection π T * Q : T * T * Q → T * Q and Ω T * Q is the musical isomorphism induced from the canonical symplectic form Ω T * Q on T * T * Q. Hence, momentum-Euler equations can be written asΠ If we change coordinates to in system (131b), we obtain the equations for compressible fluids in standard formulation. The first of this substitution in Equation (132) is simple relation between velocity and momentum, and the second one is related to Bernoulli's theorem for isentropic fluid flows.
Here, h (ρ) = ρw + w is the enthalpy function and w = w (ρ) is the internal energy of the continuum. Yet another form of Equation (131) iṡ where φ = δE/δρ and X k = δE/δm k are chemical potential and velocity. These are the usual equations for fluid mechanics in absence of entropy (or isentropic), see, e.g., the work by the authors of [13]. Entropy density can be added as follows. In kinetic theory entropy density is defined as where η( f ) is a real function of real variable f , e.g., −k B f (ln(h 3 f ) − 1) for ideal gases and k B and h are the Boltzmann and Planck constants, respectively. Evolution of this field is then given by Equation (129), which is the usual law of entropy conservation. However, to recover the antisymmetric coupling between s and M i , one should add the term −s∂ i E s to the evolution equation for M i . The evolution equations are then completely equivalent to the evolution equations of fluid mechanics coming from the underlying Poisson bracket, e.g., the work by the authors of [13].

The 10-Moment Approximation
In this section, we present a generalization of the momentum-Euler equations to ten kinetic moments (1 density + 3 momentum densities + 6 second moments). The procedure is analogical to construction of the momentum-Euler equations with the only difference that along the φ and X i fields there is a R ij tensor field on the base manifold Q.

Double GCCL of the Second Order Tensor Field
Consider a second order symmetric contravariant tensor field on Q given by representing chemical potential, velocity, and the conjugate pressure tensor. Using the mapping (92), we define the following Hamiltonian function, on T * Q, The GCCL of X is then which is a vector field on T * Q. The subsequent GCCL ofX giveŝ which is a vector field on T * T * Q.

Vertical Representative
The vertical representative V(X) of the second GCCL of X is whereX v =X −X. This vector field has only components in the directions of Π i and Π i , and the components are then interpreted as evolution equations for Π i and Π i , which contain as a special case R ij = 0 the momentum-Euler equations (128).
The distribution function f (q, p) = ∂Π i ∂q i − ∂Π i ∂q i then evolves as This equation can be interpreted physically as kinetic theory of fluctuations around mean fields of velocity, chemical potential and conjugate pressure tensor.

Projection to Moments
Subsequent projection to density, momentum density, and kinetic stress tensor, ρ(q) = dp f (q, p) (144a) M i (q) = dpp i f (q, p) (144b) P ij (q) = dpp i p j f (q, p) (144c) leads to evolution equations for the projected variables (by chain rule) where Q ijk = dpp i p j p k f (q, p) are the higher higher moments. The evolution equation for the stress tensor is thus not in a closed form, which is typical in the Grad hierarchy, e.g., the work by the authors of [62].
Let us now seek an appropriate closure, i.e., specification of the Q ijk terms. Note first that the fields φ, X i and R ij can be interpreted as the corresponding derivatives of energy with respect to the kinetic moments, In other words, φ is chemical potential, X i is velocity, and R ij is the conjugate variable to the second kinetic moments. We now seek the closure so that energy is conserved regardless the choice of energy, which requires equations (145) to possess antisymmetric structure. The coupling between the evolution of ρ and evolution of M i is antisymmetric as can be directly verified by construction of the generating Poisson bracket. To make the coupling between m i and Q ijk antisymmetric as well, we have to set Q ijk = 0, which is the sought closure. Besides automatic energy conservation, Jacobi identity is then satisfied for the evolution equations as shown, for instance, in the work by the authors of [73]. The closure can be thus referred to as a Hamiltonian closure.

Adding Entropy
As in the case of momentum-Euler equations, entropy density can be defined as s(q) = dpη( f ).
It evolves due to Equation (143) as is the first entropic moment (considered also in the work by the authors of [62]. Interestingly, the evolution for the first entropic moment is (using again Equation (143)) where is the tensor of second entropic moments. This way a hierarchy of entropic moments coupled to the hierarchy of kinetic moments can be constructed (similar to kinetic theory of non-ideal gases [74]).
The coupling to the kinetic moments is made antisymmetric by adding complementary terms (in the case of s among the state variables) to the equations for M i and P ij .

Central Kinetic Moments
In the evolution equations for the kinetic moments (145), it is interesting that the density evolves not solely due to the advection by velocity X i . This might appear strange at first sight, but it is actually due to the definition of P ij as the second kinetic moments, not the central second kinetic moments, The transformation of variables from (ρ, M, P) to (ρ, M,P) turns derivatives of energy to The evolution equation for density (145a) then becomes which already has the usual form.

Discussion and Conclusions
An injective Lie algebra homomorphism called generalized complete cotangent lift (GCCL) was defined in Equation (95), mapping the space TQ of generalized symmetric contravariant tensor fields on a manifold Q to the space g = X ham (T * Q) of Hamiltonian vector fields on the cotangent bundle T * Q. It has been shown that kinetic moments in Equation (115) of the momentum-Vlasov equations represent Poisson mappings obtained by the dualization of this homomorphism. The configuration space of compressible isentropic fluids is the semidirect product group S = Di f f (Q) F (Q), and its Lie algebra, s = X (Q) F (Q), is a subalgebra of TQ. Restriction of GCCL on s gives that embedding s → g, whereas the intertwining lemma and Equation (111) establish this embedding on the group level, that is, S → G = Di f f can (T * Q). The dual of the embedding s → g, called m-Vlasov to fluid map in Equation (120), is a Poisson and momentum mapping relating m-Vlasov bracket in Equation (107) and compressible fluid bracket in Equation (105). Generalized complete cotangent lift (X, φ) of a pair (X, φ) ∈ s is a Hamiltonian vector field on T * Q. We have introduced momentum-Euler equations, by taking the vertical representative of the complete lift of (X, ρ). It is shown that after the substitution of "m-Vlasov to fluid map" and the coordinate transformation in Equation (132), momentum-Euler equations reduce to compressible fluid equation in the classical form. Thus, it is achieved to arrive Euler's equation starting from the particle motion and applying pure geometric operations. Note that this geometrization procedure does not need any Hamiltonian functional or Poisson structure. Finally, the reversible evolution equations for ten kinetic moments are found by the same procedure, and a hierarchy of entropic moments coupled to the kinetic moments is identified including entropy. The approach to mechanics presented in this paper can lead to another geometrization of the GENERIC framework [13,17], as it does not rely on Poisson brackets although being purely geometric.