Multiscale thermodynamics of charged mixtures

TL;DR

This paper develops a comprehensive multiscale thermodynamic framework for charged fluid mixtures, integrating continuum mechanics, electrodynamics, and thermodynamics using geometric and Hamiltonian methods, with applications to engineering models.

Contribution

It introduces a novel multiscale theory combining geometric mechanics and thermodynamics for charged mixtures, unifying reversible and irreversible processes within the GENERIC framework.

Findings

Derived a new Poisson bracket for the Hamiltonian system.

Implemented gradient dynamics for irreversible thermodynamics.

Recovered standard engineering models as low-level limits.

Abstract

A multiscale theory of interacting continuum mechanics and thermodynamics of mixtures of fluids, electrodynamics, polarization and magnetization is proposed. The mechanical (reversible) part of the theory is constructed in a purely geometric way by means of semidirect products. This leads to a complex Hamiltonian system with a new Poisson bracket, which can be used in principle with any energy functional. The thermodynamic (irreversible) part is added as gradient dynamics, generated by derivatives of a dissipation potential, which makes the theory part of the GENERIC framework. Subsequently, Dynamic MaxEnt reductions are carried out, which lead to reduced GENERIC models for smaller sets of state variables. Eventually, standard engineering models are recovered as the low-level limits of the detailed theory. The theory is then compared to recent literature.

Tables

Table 1. Table 1: List of the basic quantities, symbols and SI units.

Quantity	Symbol	SI units
magnetic field	$𝐁$	$\mathrm{kg}/s^2/A$
magnetic field intensity	${𝐁}^{†},𝐇$	$A{m}^{-1}$
displacement field	$𝐃$	$C/m^2$
electric induction field	$𝓓$	$C/m^2$
electric field	$𝐄,{𝐃}^{†},{𝓓}^{†}$	$V{m}^{-1}$
electric potential	$\phi$	$V$
polarization density	$𝐏$	$C/m^2$
polarization momentum	$𝝅$	$Vs{m}^{-1}$
magnetization density	$𝐌$	$A{m}^{-1}$
partial mass density	${\rho }_{\alpha }$	$\mathrm{kg}/m^3$
chemical potential	${\rho }_{\alpha }^{†}$	$J{\mathrm{kg}}^{-1}$
mass momentum density	$𝐮$	$\mathrm{kg}/m^2/s$
total momentum density	$𝐦,\widehat{𝐦}$	$\mathrm{kg}/m^2/s$
entropy density	$s$	$J/K/m^3$

Table 2. Table 2: List of momenta densities and the corresponding levels of description.

momentum density	level of description
$𝐮$	$({\rho }_{\alpha },s,𝐮)$
$𝐦$	$({\rho }_{\alpha },s,𝐦,𝐃,𝐁)$
$\widehat{𝐦}$	$({\rho }_{\alpha },s,\widehat{𝐦},𝐏,𝝅,𝓓,𝐁,𝐌)$

Equations

$$\displaystyle\dot{F}(\mathbf{x})=\Big{\{}F,E\Big{\}}+\left\langle\frac{\delta F}{\delta\mathbf{x}},\frac{\delta\Xi}{\delta\mathbf{x}^{*}}\middle|_{\mathbf{x}^{*}=\frac{\delta S}{\delta\mathbf{x}}}\right\rangle,$$

$$\left(\frac{\delta F}{\delta \mathbf{x}}\right)\equiv F_{\mathbf{x}}$$

$$\left\{F,E\right\}={\left(\dot{F}\right)}_{\text{rev}}={\left({\partial }_{t}F\right)}_{\text{rev}}.$$

$$⟨F_{\mathbf{x}},{\Xi }_{{\mathbf{x}}^{\ast }}{|}_{{\mathbf{x}}^{\ast }=S_{\mathbf{x}}}⟩={\left(\dot{F}\right)}_{\text{irr}}={\left({\partial }_{t}F\right)}_{\text{irr}}.$$

$${\partial }_{t}F={\left({\partial }_{t}F\right)}_{\text{rev}}+{\left({\partial }_{t}F\right)}_{\text{irr}}.$$

$$(\dot{\mathbf{x}}{)}_{\text{rev}}=\{\mathbf{x},E\},$$

$$\{\bullet ,\bullet \}:\mathcal{F}\times \mathcal{F}⟶\mathcal{F}$$

$$\{FG,H\}=\{F,H\}G+F\{G,H\},$$

$$\{\{F,G\},H\}+\{\{G,H\},F\}+\{\{H,F\},G\}=0,$$

$$\{F,H{\}}^{T^{\ast }V}(a,a^{\ast })=⟨F_a,H_{a^{\ast }}⟩-⟨H_a,F_{a^{\ast }}⟩,$$

$$[\bullet ,\bullet ]:\mathfrak{g}\times \mathfrak{g}⟶\mathfrak{g}$$

$$\{F,H{\}}^{{\mathfrak{g}}^{\ast }}(\mu )=\pm ⟨\mu ,[F_{\mu },H_{\mu }]⟩$$

$$\{F,H{\}}^{P_1\times P_2}=\{F,H{\}}^{P_1}+\{F,H{\}}^{P_2}$$

$$\{F,H{\}}^{{\mathfrak{g}}_1^{\ast }\times {\mathfrak{g}}_2^{\ast }}=\{F,H{\}}^{{\mathfrak{g}}_1^{\ast }}+\{F,H{\}}^{{\mathfrak{g}}_2^{\ast }},$$

$$\rho :\mathfrak{g}\times V\mapsto V.$$

$$\left[\left({\xi }_1,v_1\right),\left({\xi }_2,v_2\right)\right]=\left(\left[{\xi }_1,{\xi }_2\right],{\rho }_{{\xi }_1}(v_2)-{\rho }_{{\xi }_2}(v_1)\right).$$

$$\{F,H{\}}^{{\mathfrak{g}}^{\ast }\times V^{\ast }}(\mu ,a^{\ast })=\pm ⟨\mu ,\left[F_{\mu },H_{\mu }\right]⟩\pm ⟨a^{\ast },{\rho }_{F_{\mu }}(H_{a^{\ast }})-{\rho }_{H_{\mu }}(F_{a^{\ast }})⟩,$$

$$\{F,H{\}}_{(DP)}^{T^{\ast }V\times {\mathfrak{g}}^{\ast }}=\{F,H{\}}^{T^{\ast }V}\pm \{F,H{\}}^{{\mathfrak{g}}^{\ast }}.$$

$$⟨a^{\ast },{\rho }_{\xi }(a)⟩=⟨{\rho }_{\xi }^{\ast }(a^{\ast }),a⟩$$

$$\{F,H{\}}_{(SDP)}^{T^{\ast }V\times {\mathfrak{g}}^{\ast }}=$$

$$\pm ⟨a^{\ast },{\rho }_{F_{\mu }}(H_{a^{\ast }})-{\rho }_{H_{\mu }}(F_{a^{\ast }})⟩\pm ⟨a,{\rho }_{F_{\mu }}^{\ast }(H_a)-{\rho }_{H_{\mu }}^{\ast }(F_a)⟩.$$

$$\{F,G{\}}^{(\text{FM})}(\rho ,\mathbf{u},s)$$

$$+\int \mathrm{d}\mathbf{r}{u}_i\left({\partial }_j{F}_{u_i}{G}_{u_j}-{\partial }_j{G}_{u_i}{F}_{u_j}\right)$$

$$+\int \mathrm{d}\mathbf{r}s\left({\partial }_i{F}_s{G}_{u_i}-{\partial }_i{G}_s{F}_{u_i}\right),$$

$$(\dot{F}{)}_{\text{rev}}=$$

$$=$$

$$+\int \mathrm{d}\mathbf{r}{F}_{u_i}\left(-\rho {\partial }_i{\rho }^{†}-u_j{\partial }_i{u}^{†j}-s{\partial }_i{s}^{†}-{\partial }_j(u_i{u}^{†j})\right)$$

$$+\int \mathrm{d}\mathbf{r}{F}_s\left(-{\partial }_i(s{u}^{†i})\right),$$

$$\dot{F}=\int \mathrm{d}\mathbf{r}\left(F_{\rho }{\partial }_t\rho +F_{u_i}{\partial }_t{u}_i+F_s{\partial }_{t}s\right)$$

$$({\partial }_t\rho {)}_{\text{rev}}=$$

$$({\partial }_t{u}_i{)}_{\text{rev}}=$$

$$({\partial }_{t}s{)}_{\text{rev}}=$$

$$E^{\text{Euler}}\left(\rho ,\mathbf{u},s\right)=\int \mathrm{d}\mathbf{r}\left(\frac{{\mathbf{u}}^2}{2\rho }+\epsilon (\rho ,s)\right),$$

$$\{F,G\}({\rho }_{\beta },{\mathbf{u}}^{\beta },s_{\beta })=\sum _{\alpha =0}^n\{F,G{\}}^{{\text{(FM)}}_{\alpha }}({\rho }_{\alpha },{\mathbf{u}}^{\alpha },s_{\alpha }).$$

$$\mathbf{u}=\sum _{\alpha =0}^n{\mathbf{u}}^{\alpha }\text{mbox}ands=\sum _{\alpha =0}^n{s}_{\alpha }.$$

Full text

1 Introduction

Theoretical electrochemistry aims to describe and predict behavior of chemically reacting systems of charged substances. The modeling methods vary according to the characteristic times, lengths and details of the observed electrochemical systems. This paper aims to develop a hierarchy of continuum models on different levels of description using the framework of the General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) Grmela1997 ; ottinger1997dynamics ; ottinger2005minimal ; pavelka2018multiscale .

Let us first briefly recall GENERIC. Consider an isolated system described by state variables $\mathbf{x}$. The state variables can be for instance position and momentum of a particle, field of probability density on phase space, fields of density and momentum density, electromagnetic fields, etc. Evolution of functionals $F(\mathbf{x})$ of the state variables is then expressed as

[TABLE]

where the former term on the right hand side stands for a Poisson bracket of the functional $F$ and the energy $E$ while the latter for scalar product of gradient of $F$ and gradient of the dissipation potential $\Xi$. Conjugate variables (derivatives of the entropy $S$ in the entropic representation) are denoted by $\mathbf{x}^{*}$. The Poisson bracket is antisymmetric, which leads to automatic energy conservation, and satisfies Jacobi identity, which expresses intrinsic compatibility of the reversible evolution. The irreversible term yields a generalized gradient flow driven by gradient of entropy and ensures the second law of thermodynamics. Many successful models in non-equilibrium thermodynamics have been formulated in the GENERIC structure (1), and many new thermodynamically consistent models have been obtained by seeking that structure, see e.g. ottinger2005minimal ; pavelka2018multiscale .

We shall start the present work by recalling the semi-direct product coupling of reversible fluid mechanics and electrodynamics in vacuum marsden1982hamiltonian . In other words, we let the electrodynamics be advected by fluid mechanics as in pavelka2018multiscale . To go beyond, we shall add also the field of polarization density and a canonically coupled momentum of polarization. It is important to express the behavior of dipole moment of molecules in interaction with electromagnetic field and the overall motion. To this picture, we shall add magnetization (the famous Landau & Lifshitz model) advected by fluid mechanics. This way we shall build a hierarchy of levels of description with appropriate Poisson brackets expressing kinematics on the levels, see Figure 1.

Subsequently, we shall introduce dissipation on the most detailed levels of description, we shall present reductions to less detailed (lower) levels down to the level of mechanical equilibrium. On this level the evolution is governed by the generalized Poisson-Nernst-Planck equations. We believe that such a complete and geometric picture of continuum thermodynamics of matter coupled with electrodynamics (including polarization and magnetization) was missing in the literature.

1.1 Mathematical aspects, terminology and further remarks

A description in terms of fields describing physical quantities, e.g. mass density, magnetic field, is developed throughout the manuscript. A mathematically rigorous formulation of the Poisson brackets and the dissipation potentials for the fields would require deployment of an advanced functional analysis, see e.g. marsden2013introduction . We opt-out from discussing the setting of the function spaces, convergence of the integrals and other mathematical aspects that are not pertinent to the discussed physics at this stage of the development.

Only the bulk equations are investigated in this manuscript, and all boundary terms are thus neglected in the calculations below. This means that the considered systems are assumed to be either boundary-less, isolated or (infinitely) large with a reasonable decay of the fields at infinity. This is a standard assumption in geometric field theories, see e.g. ottinger2005minimal ; pavelka2018multiscale ; marsden2013introduction .

A level of description is given by a set of variables. If two sets of variables, belonging to the respective levels of description, can be ordered by inclusion, then we say that the level of description with smaller set of variables is less detailed level of description and vice versa.

This brings us to an another interpretation of the multiscale description not necessarily connected to different space-time scales. For instance in the Grad hierarchy grad ; miroslav-grad , where kinetic moments are considered struch , the scale (or rather level of description) is given by the number of moments considered. The advantage of the latter meaning is that it is independent of the particular non-objective and observer-dependent spatio-temporal scale.

The term reversible is used through out the article in connection with the structure of the evolution equations in the sense of the time-reversal transformation pavelka2014time . The Hamiltonian systems may often be cast in to the first-order, symmetric hyperbolic partial differential equation form peshkov2018continuum . Assuming smooth initial data, solutions of the hyperbolic systems may be constructed for short times using the Cauchy-Kowalevskaya theorem evansPDE . Such solutions are then short-time reversible in the usual sense of partial differential equations.

1.2 Notation $(\partial_{t}\mathbf{x})_{\text{rev/irr}}$ and $\mathbf{x}^{\dagger}$, $\mathbf{x}^{*}$

A problem under consideration usually admits multiple equivalent choices of state variables. Let us denote the family of the admissible, equivalent, levels of description as $\mathbb{X}$. For a given level of description, characterized by state variables $\mathbf{x}\in\mathbb{X}$, the presented theory aims to describe evolution of functionals depending on $\mathbf{x}$, i.e. we are interested in $\dot{F}(\mathbf{x})$, see (1). Therefore, functional derivatives

[TABLE]

are often used in this manuscript.

The part of the evolution generated by the Poisson bracket $\{\cdot,\cdot\}$ and the energy $E(\mathbf{x})$, see (1), will be called Hamiltonian or reversible and denoted as

[TABLE]

The functional derivative of the energy w.r.t. variable $\mathbf{x}$, also called the energy-conjugate to the variable $\mathbf{x}$, will be briefly denoted by $\dagger$ in the superscript, i.e. $E_{\mathbf{x}}=\mathbf{x}^{\dagger}$. Poisson brackets are usually conveniently expressed in the energetic representation, i.e. when entropy density $s$ is amongst the state variables, $(s,\boldsymbol{\xi})\in\mathbb{X}$, where $\boldsymbol{\xi}$ denotes the state variables other than the entropy density.

The part of the evolution generated by the dissipation potential $\Xi$ and the entropy $S(\mathbf{x})$, see (1), will be called irreversible and denoted as

[TABLE]

The symbol $\mathbf{x}^{*}$, called entropy-conjugate, denotes the variables on which the dissipation potential $\Xi(\mathbf{x}^{*})$ depends. The notation $\mathbf{x}^{*}$ is going to be overloaded in Section 3 because it will be also used to denote the value of the entropy derivative, i.e. $\mathbf{x}^{*}=S_{\mathbf{x}}$. Since the irreversible evolution is expressed in terms of the derivatives of the dissipation potential w.r.t. $\mathbf{x}^{*}$, presence of energy density $e$ amongst the variables makes it simple to control the energy conservation, i.e., $(e,\boldsymbol{\xi})\in\mathbb{X}$. Such set of variables is called entropic representation callen .

The relations between the energy-conjugates in the energetic representation $(s^{\dagger},\boldsymbol{\xi}^{\dagger})$ and the entropy-conjugate in the entropic representation $(e^{*},\boldsymbol{\xi}^{*})$ are explained in pavelka2018multiscale ; CMAT2018 .

The overall evolution of a functional $F$ is then composed from the reversible and irreversible parts,

[TABLE]

Symbols $(\partial_{t}F)_{\text{rev/irr}}$ be interpreted as the reversible and irreversible (w.r.t TRT) parts of the right hand side of the evolution equations, or as the reversible and irreversible vector fields, see e.g. pavelka2014time .

2 Hierarchy of Poisson brackets

2.1 Hamiltonian dynamics and semidirect product theory

A physical system described by state variables $\mathbf{x}$ is said to be Hamiltonian if the differential equations governing its dynamics can be written in the form of

[TABLE]

where dot stands for (typically partial) time derivative, and $E$ is a Hamiltonian (energy) function marsden2013introduction ; de2011methods ; Marsden1983 . Here, the (Poisson) bracket

[TABLE]

is a (skew-symmetric and bilinear) mapping on the space of smooth functions depending on the state variables $\mathbf{x}$. It is required that a Poisson bracket must satisfy both the Leibniz and the Jacobi identities weinstein1998poisson that is

[TABLE]

for any real valued functions $F$, $G$ and $H$. A space endowed with such a bracket is called a Poisson space abraham1978foundations ; Marsden-Ratiu-book ; libermann2012symplectic . Skew-symmetry of the bracket manifests conservation of the Hamiltonian function all along the motion. This structure implies in particular to conservation of energy. This observation is fundamental for reversible character of the Hamiltonian systems. More on time-reversibility and Onsager-Casimir reciprocal relations can be found in pavelka2018multiscale .

In order to write a system of equations in the Hamiltonian form (6), two tasks must simultaneously be accomplished after deciding the state variables abraham1978foundations . These are (1) to choose a Hamiltonian function and (2) to construct a Poisson bracket. Even though there is no algorithmic way to accomplish these tasks, there are some techniques to perform Hamiltonian analysis of a given system. Let us depict here one of them that we use in the present work. Start with the most basic form of the system (by removing all possible field extensions and symmetries). After the Hamiltonian analysis of the basic model is established, field extensions can be added by applying pure geometric and algebraic techniques. For example, consider a classical hydrodynamical system, start initially with analysis of incompressible and isentropic fluid flow. After this is done, let the mass density and the entropy vary in time. Additions of these scalar fields can be done in purely abstract framework by using the semi-direct product theory marsden-ratiu-weinstein ; esen2017hamiltonian . Accordingly, in this work, we shall start with Hamiltonian realizations of two basic models, namely the one-component compressible fluid mechanics. Then we extend these models with the additions of some field theories such as electromagnetic fields, polarization and magnetization.

For the basic models (without field extensions and symmetries), in the literature, there exist some popular Poisson structures to begin with. Let us comment on these Poisson structures one by one in order to fix the notation and to make the present work more complete.

Canonical Poisson bracket.

Consider a vector space $V$ and its linear algebraic dual $V^{*}$. In finite dimensions, the duality (pairing) $\langle\bullet,\bullet\rangle$ between $V$ and $V^{*}$ can be considered as the Euclidean scalar (dot) product whereas, in infinite dimensions, (that is if $V$ is a function space or the space of some fields) the duality can be an $L^{2}$-pairing (that is simply multiply-and-integrate form). The cotangent bundle $T^{\ast}V$ is isomorphic to the product space $V\times V^{\ast}$ and carries a canonical Poisson structure. Let us exhibit this geometry. The elements of the cotangent bundle $T^{*}V$ are pairs $(a,a^{*})$ in $T^{\ast}V$. Consider two function(al)s $F$ and $H$ on $T^{*}V$, that depend on $(a,a^{*})$. By referring to the duality between $V$ and $V^{*}$, the canonical Poisson bracket is defined to be

[TABLE]

where, for example, $F_{a}$ denotes the partial (for infinite cases, Fréchet) derivative of $F$ with respect to $a$. We assume that the vector space is reflexive that is $V^{**}=V$. This enables us consider $F_{a}$ as an element of $V^{*}$ for $a$ in $V$, and $F_{a^{*}}$ as an element of $V$ for $a^{*}$ in $V^{*}$.

Lie-Poisson bracket.

Let us start once more with a vector space $\mathfrak{g}$. In this case, we are not coupling $\mathfrak{g}$ with its dual as in the previous case, but we consider $\mathfrak{g}$ as a Lie algebra. A vector space $\mathfrak{g}$ is called a Lie algebra if it admits a bilinear skew-symmetric mapping (Lie bracket)

[TABLE]

satisfying the Jacobi identity. In this case, linear algebraic dual $\mathfrak{g}^{*}$ of $\mathfrak{g}$ carries a Poisson bracket called as the Lie-Poisson bracket. More concretely, if $F$ and $H$ are two function(al)s on $\mathfrak{g}^{*}$ then, for $\mu$ in $\mathfrak{g}^{*}$,

[TABLE]

where $F_{\mu}$ is assumed to be an element of $\mathfrak{g}$ whereas $\langle\bullet,\bullet\rangle$ is the pairing between $\mathfrak{g}^{*}$ and $\mathfrak{g}$. Notice that there exist two Lie-Poisson brackets according to the sign of the bracket (10). The choice of the sign is not arbitrary for a chosen physical system. It is minus for the rigid body dynamics, since it admits a symmetry due to the left action that is, matrix multiplication of the rotation group $SO(3)$ on the configuration space $\mathbb{R}^{3}$, see Fecko ; Marsden-Ratiu-book . In this present work, we employ the plus Lie-Poisson bracket since, in the continuum theories there exists a right symmetry, physically called as the particle relabeling symmetry MaWe81 ; Marsden1983 . We will continue to keep plus minus notation in this introductory subsection. After having presented the canonical and Lie-Poisson structures, let us now briefly explain (both direct and semi-direct) couplings and extensions of these geometries.

Direct product of two Poisson structures.

The most direct way to couple two Poisson structures is to add the associated Poisson brackets directly without any cross terms and interactions. Consider for example two Poisson spaces $P_{1}$ and $P_{2}$ equipped with Poisson brackets $\{\bullet,\bullet\}^{P_{1}}$ and $\{\bullet,\bullet\}^{P_{2}}$, respectively. Then the direct product Poisson structure is

[TABLE]

for two function(al)s on $P_{1}\times P_{2}$. In particular, consider two Lie algebras $\mathfrak{g}_{1}$ and $\mathfrak{g}_{2}$, and the dual spaces $\mathfrak{g}^{*}_{1}$ and $\mathfrak{g}^{*}_{2}$, respectively. Then by employing the isomorphism $(\mathfrak{g}_{1}\times\mathfrak{g}_{2})^{*}=\mathfrak{g}_{1}^{*}\times\mathfrak{g}_{2}^{*}$, we construct the Lie-Poisson bracket on the (direct product) dual space by

[TABLE]

where the Lie-Poisson brackets on the right hand side are the ones on $\mathfrak{g}_{1}^{*}$ and $\mathfrak{g}_{2}^{*}$, respectively.

Semi-direct product extension.

We now extend the Lie-Poisson bracket. Consider a Lie algebra $\mathfrak{g}$ and a vector space $V$. Assume that $\mathfrak{g}$ acts on $V$, that is there exist bilinear mapping

[TABLE]

By fixing an element $\xi$ in $\mathfrak{g}$, we define a linear mapping $\rho_{\xi}$ on $V$ Hall-LA-book . This enables us to define a Lie algebra bracket on the product space $\mathfrak{g}\times V$ given by

[TABLE]

By referring to the Lie-Poisson structure presented in the previous paragraph, the dual space $\mathfrak{g}^{*}\times V^{*}$ turns out to be a Poisson space equipped the Lie-Poisson bracket. For two function(al)s $F$ and $H$ depending on $(\mu,a^{*})$ in $\mathfrak{g}^{*}\times V^{*}$, the Lie-Poisson bracket is computed to be

[TABLE]

where the first pairing on the right hand side is the one between $\mathfrak{g}^{*}$ and $\mathfrak{g}$ whereas the second pairing on the right hand side is the one between $V^{*}$ and $V$.

Direct product coupling of the canonical Poisson bracket and the Lie-Poisson bracket.

Consider two Poisson structures, the canonical Poisson bracket (8) on $T^{*}V$ and the Lie-Poisson bracket (10) on $\mathfrak{g}^{*}$. If a physical system has a phase space as the product of these two spaces $\mathfrak{g}^{*}\times T^{*}V$, then one way to define a Poisson bracket on the product space by simply adding two brackets (8) and (10). This is, for two function(al)s $F$ and $H$ depending on $(a,a^{*},\mu)$, given by

[TABLE]

Semi-direct product coupling of the canonical Poisson bracket and the Lie Poisson bracket.

If the quantities in $V$ are advected by the quantities in $\mathfrak{g}$, that is if there is an action of $\mathfrak{g}$ as given in (13), then the Poisson bracket in (16) fails to be the correct Poisson structure for the system. In this case, the cross terms (see, for example, the second term in (15)) must be added as in the semi-direct product theory exhibited in the previous paragraph. In this case, one is equipped, in addition to the action (13), with a (dual) action of $\rho^{*}_{\xi}$ on the dual space $V^{*}$. This is due to the preservation of the duality between $V$ and $V^{*}$ under the action,

[TABLE]

for all $a$ in $V$ and $a^{*}$ in $V^{*}$. Note that the dual action $\rho^{*}$ is computed as linear algebraic dual of the linear mapping $\rho_{\xi}$. A direct computation shows that for two function(al)s $F$ and $H$, under the existence of the actions, the Poisson bracket is

[TABLE]

Notice that the first line is the same with the direct product case (non interaction) in (16), and the second line is a manifestations of the advection. This is a particular case of the matched pair of Hamiltonian dynamics presented in ogul2016matched2 .

2.1.1 Matched pairs

Only one sided couplings and semidirect products are considered in this work. Only recently the two-sided actions have been developed in classical physics, called Matched pairs ogul2016matched2 . In the Galilean setting the one-sided actions are usually enough, but it is indeed possible that matched pairs would turn necessary for the fully relativistic and quantum (see works by Majid Majid-book ) treatment. So far the only physically important matched pair in classical physics has been identified in Grad hierarchy, where coupling of fluid mechanics and higher kinetic moments plays an important role momentum-Euler .

2.2 Hamiltonian fluid mechanics of mixtures

The Poisson bracket expressing kinematics of fluid mechanics has been long known arnold ; dv ; GrPD ; MG84 ; Morrison-brackets . The Poisson bracket can be easily extended to mixtures with multiple densities, momenta and entropies (i.e. temperatures), see e.g. hierarchy .

2.2.1 Classical fluid mechanics

The state variables of classical fluid mechanics are fields of density, momentum density and entropy density, $\mathbf{x}=(\rho,\mathbf{u},s)$. The Poisson bracket generating one-component compressible fluid mechanics (hydrodynamic Poisson bracket) is

[TABLE]

where $\rho$, $\mathbf{u}$ and $s$ are mass density, mass momentum density and volumetric entropy density, respectively. The Poisson bracket in (2.2.1) is also an example of a semidirect product bracket obeying the form exhibited in (15). This well-know analysis can be found, for example, in arnold ; marsden-ratiu-weinstein ; be .

For an arbitrary energy $E(\rho,\mathbf{u},s)$, the reversible evolution of a functional $F$ of the state variables reads

[TABLE]

where integration by parts was used several times111The energy-conjugate to variable $\mathbf{x}$, i.e. functional derivative of the energy E w.r.t. $\mathbf{x}$, is going to be denoted by $\mathbf{x}^{\dagger}$, i.e. $E_{\mathbf{x}}=\mathbf{x}^{\dagger}$. The notation is summarized in section 1.2.. Boundary terms disappear as we assume isolated (e.g. periodic) system.

By comparing with the chain rule

[TABLE]

we can read the reversible evolution equations for fluid mechanics,

[TABLE]

For the usual choice of energy,

[TABLE]

the compressible non-isothermal Euler equations for ideal fluids are obtained from the system (22). The energy conjugates $\rho^{\dagger}$, $\mathbf{u}^{\dagger}$ and $s^{\dagger}$ represent the chemical potential, barycentric velocity and temperature, respectively.

2.2.2 Hamiltonian fluid mechanics of mixtures

Consider now a mixture of $n+1$ species, each of which is described by its own density, momentum density and entropy density. The Poisson bracket expressing kinematics of state variables $\mathbf{x}=(\rho_{\alpha},\mathbf{u}^{\alpha},s_{\alpha})$, $\alpha\in\{0,1,\cdots,n\}$, is

[TABLE]

This Poisson bracket can be derived for instance by projection from the Liouville equation hierarchy . It consists of the sum of $n+1$ Poisson brackets (2.2.1), each expressed in terms of variables of mixture component $\alpha$. The Poisson bracket in (24) is a direct product Poisson bracket fitting the abstract framework in (11).

Poisson bracket (24) depends on $n+1$ momenta and $n+1$ entropies, each for one component of the mixture, which is a rather detailed description allowing for independent motion of the constituents and for different temperatures of them (as in cold plasma, where electrons have different temperature than ions). We are, however, often interested in less detailed description, keeping only densities of the species, the total mass momentum and the total entropy,

[TABLE]

By letting the arbitrary functional depend only on state variables $\mathbf{x}=(\rho_{\alpha},\mathbf{u},s)$, bracket (24) becomes

[TABLE]

which is referred to as the classical mixture hydrodynamic bracket. It generates the reversible part of Classical Irreversible Thermodynamics (CIT) degroot1984nonequilibrium .

The descriptions of the fluid and fluid mixture dynamic considered in the remainder of the paper will be based on the brackets (2.2.1) and (26), respectively. Hence, the description considering the distinguished momenta of species (24) is further avoided.

2.3 Electrodynamics in vacuum

The reversible evolution of electromagnetic fields is generated by the canonical Poisson bracket, see marsden1982hamiltonian ; holm1986hamiltonian ; pavelka2018multiscale ,

[TABLE]

where A stands for the vector potential and $\mathbf{Y}=-\mathbf{D}$ denotes negative of the electric displacement field (either in variables $(\mathbf{A},\mathbf{Y})$ or $(\mathbf{A},\mathbf{D})$). Poisson bracket in (27) is the canonical Poisson bracket as in (8). Here, the vector space $V$ is consisting of vector potentials whereas the dual space $V^{*}$ is the space of electric displacement fields.

Let us define the magnetic field $\mathbf{B}$ as

[TABLE]

In order to express the bracket (27) in terms of magnetic field, we assume that the functionals depend only on the curl of $\mathbf{A}$. Bracket (27), transformed in terms of $(\mathbf{D},\mathbf{B})$, see esen2017hamiltonian , becomes

[TABLE]

This is the Poisson bracket expressing kinematics of electromagnetic fields $\mathbf{D}$ and $\mathbf{B}$.

For an arbitrary energy $E$ the evolution equations of the electromagnetic field given by (29) are

[TABLE]

The conjugates are interpreted as electric field and magnetic intensity in the laboratory frame,

[TABLE]

A concrete energy will be specified later after matter is added to the system so that the theory becomes Galilean invariant.

Applying divergence to (30) gives the following evolution equations:

[TABLE]

The first equality holds due to (28), the second equality demonstrates the absence of charges in a vacuum. Hence, the usual constraints–Gauß’s laws  jackson1998classical –hold true if satisfied by the initial condition, see marsden1982hamiltonian . Equations (32) are thus a consequence of evolution equations (30).

2.4 Electromagnetic field advected by charged fluids

The purpose of this section is to formulate coupled kinematics of fluids and electromagnetic fields. We employ the theory of semidirect product to find such coupling, and then we perform a transformation unveiling the usual form of the Lorentz force.

2.4.1 Semidirect product extension

Let us recall the Lie-Poisson structure in (26) expressing kinematics of Classical Irreversible Thermodynamics, i.e. for state variables $(\rho,\mathbf{u},s)$, and the canonical Poisson bracket in (29) for the electromagnetic field theory in terms of the fields $(\mathbf{D},\mathbf{B})$. In this Subsection, we couple these two in the light of the abstract framework in (18) assuming that the action is minus the Lie derivative $-\mathcal{L}$ with respect to the velocity field $\mathbf{u}^{\dagger}$. The physical meaning is that a quantity is Lie-dragged or advected by velocity of the continuum (conjugate to the total momentum), see e.g. arnoldbook ; LagEul . We obtain Poisson bracket

[TABLE]

where $\mathbf{m}$ denotes total momentum density (of matter and electromagnetic field),

[TABLE]

This standard relation, see e.g. Landau2 , can be either verified by a direct calculation pavelka2018multiscale or inferred using physical arguments, see Appendix D. Note that the notation introduced in (33) allows to write briefly

[TABLE]

Bracket (33) expresses kinematics of a CIT mixture and electromagnetic field with state variables $\mathbf{x}=(\rho_{\alpha},\mathbf{m},s,\mathbf{D},\mathbf{B})$ and was found (for the single species case) in ADER ; CMAT2018 ; holm1986hamiltonian .

2.4.2 Transformation to mass momentum

Let us suppose that each mixture component carries charge $e_{0}\frac{z_{\alpha}}{m_{\alpha}}$ proportional to the free charge density, which is defined as

[TABLE]

Suppose functionals $F(\rho_{\alpha},\mathbf{m},s,\mathbf{D},\mathbf{B})=\widetilde{F}(\rho_{\alpha},\mathbf{u},s,\mathbf{D},\mathbf{B})$, the relation of the total and mass momenta (34) implies the following transformation rules between their derivatives, it reads

[TABLE]

see Appendix D for further details. Poisson bracket (33) can be now transformed by means of relations (37) to the mass momentum $\mathbf{u}$ instead of the total momentum $\mathbf{m}$. The calculation was carried out in esen2017hamiltonian and pavelka2018multiscale and leads to Poisson bracket

[TABLE]

which is the Poisson bracket expressing evolution of a CIT mixture coupled with electromagnetic field (using the mass momentum $\mathbf{u}$).222the transformation crucially depends on $\operatorname{div}\mathbf{B}=0$

Energy $E(\rho_{\alpha},\mathbf{m},s,\mathbf{D},\mathbf{B},\mathbf{E},\mathbf{H})$ when considering $\mathbf{m}$ as a state variable is equal to energy $\widetilde{E}(\rho_{\alpha},\mathbf{u},s,\mathbf{D},\mathbf{B})$ when considering $\mathbf{u}$ as a state variable. However, to distinguish between the conjugates, we introduce the notation $\widetilde{E}_{\mathbf{x}}=\mathbf{x}^{\widetilde{\dagger}}$. It follows from the formulae (37) that the chemical potentials, temperature and velocity are unaffected by transformation (34), i.e.

[TABLE]

respectively. Note that the electromagnetic fields on the left hand side are seen from the fluid frame while those on the right hand side from the inertial lab frame, see C for more details.

Energy can be chosen as

[TABLE]

c.f. van-magnetic . The first part represents kinetic energy of the center of mass. The second part is the electromagnetic energy, see e.g. feynman2 , where all the fields are measured in the laboratory frame. Finally, the last part is the internal energy. In order to express the energy in terms of the state variables $(\rho_{\alpha},\mathbf{v},s,\mathbf{D},\mathbf{B})$, we have to use material relations

[TABLE]

which are valid in the frame co-moving with the center of mass of the matter, i.e. with velocity $\mathbf{v}$ with respect to the laboratory frame. Symmetric positive definite tensors $\boldsymbol{\epsilon}$ and $\boldsymbol{\mu}$ represent the electric permittivity and magnetic permeability. The co-moving fields $\mathbf{D}^{\prime}$, $\mathbf{B}^{\prime}$, $\mathbf{E}^{\prime}$ and $\mathbf{H}^{\prime}$ are then transformed to the laboratory fields by transformation (130). This leads to energy

[TABLE]

This energy has the right behavior with respect to Galilean transformations, see e.g. leblond1973galilean .

The evolution equations implied by bracket (38) and energy (42) are

[TABLE]

The terms on the last line represent the Lorentz force acting upon the free charge density. Equations (43) are invariant with respect to Galilean transformations (130). Therefore, they represent a Galilean invariant form of electrodynamics coupled with matter. Note that Galilean invariance would be impossible without the explicit presence of evolution equations for matter.

The relation between formulation (38) and (33) can also be understood in terms of surface balances of the electric induction flux and magnetic induction flux for static and moving surfaces as it is shown in (marsik1998biotermodynamika, , Eqns. 2.95, 2.96)

Note that the bracket formulated with the total momentum (33) satisfies the Jacobi identity unconditionally. This is not the case for the bracket with the mass momentum (38), where the validity of the Jacobi identity cannot be established without $\operatorname{div}\mathbf{B}=0$.

2.4.3 Gauß’s law for electric charge and non-existence of magnetic monopoles

Equations (32) represent the Gauß’s laws. Let us now consider the dynamics of coupled matter and electromagnetic field generated by brackets (38) and (33). After applying divergence on Eq. (43d) we obtain that

[TABLE]

where $\rho_{f}$ stands for the free charge. The equality holds if it holds at some initial time (which can be for instance vacuum, so it indeed holds). Applying divergence on Eq. (43e), we obtain again that $\operatorname{div}\mathbf{B}=0$. The Gauß’s law is implied by the evolution equations even in the case of electrodynamics coupled with mixtures of fluids.

2.5 Polarization

The reversible evolution of a charged mixture in electromagnetic field is described by one of the Poisson brackets in the previous section and a choice of energy. But such description does not, in general, capture the intrinsic dipole moments of the molecules, i.e. polarization. Indeed, a fluid of dipoles can be charge neutral while still electromagnetically interacting. An additional bound charge is present due to internal dipole density of the matter on top of the modeled free charge.

Description of the bound charge depends profoundly on the chosen variables, the time/space scales and the internal structure of the assumed matter. The classical treatment on the macroscopic level, see e.g. purcell1985electricity , resorts to the definition of polarization vector $\mathbf{P}$. The divergence part of $\mathbf{P}$ is set equal to the density of bound charge. The time derivative of $\mathbf{P}$ represents current, therefore it is added to the left-hand side of Ampere’s law (43d).

Russakoff in russakoff1970derivation acquired the polarization as a consequence of averaging of microscopic Maxwell’s equations with point charges and subsequent expansion with respect to spatially correlated charges. This approach leads, compared to the Purcell’s, to a definition of the polarization related with magnetization proportional to the averaged relative velocities of the correlated charges.

2.5.1 Intrinsic dynamics of polarization

Density of polarization $\mathbf{P}$ represents a vector field just as the displacement field $\mathbf{D}$. Advection of the vector field by the fluid mechanics is then expressed by a semidirect product as in Sec. 2.4.1. Similarly as in that Section, we can add (besides $\mathbf{P}$) the conjugate momentum variable (to be denoted by $\boldsymbol{\pi}$). The Poisson bracket expressing advection of the pair $(\boldsymbol{\pi},-\mathbf{P})$ by fluid mechanics is shown in the following section.

The canonical Poisson bracket (8) of state variables $\boldsymbol{\pi}$ and $\mathbf{P}$ reads

[TABLE]

Bracket (45) can represent a continuum of elementary dipoles with fixed centers of mass, but changing length and orientation. Indeed, covector field $\boldsymbol{\pi}$ can be interpreted as proportional to the relative momentum of particles forming the dipole, c.f. variable t in (127) and (128). The divergence part of $\mathbf{P}$ represents the bound charge density,

[TABLE]

2.5.2 Advected cotangent bundle $(\boldsymbol{\pi},-\mathbf{P})$

When the dipoles are not fixed in space, but advected by a fluid, the interaction is captured by coupling of bracket (45) and fluid dynamics bracket (26) by semidirect product,

[TABLE]

see bracket (18) later comments for more details. This Poisson bracket expresses kinematics of state variables $(\rho_{\alpha},{\,{}^{\boldsymbol{\pi}}{}}\mathbf{u},s,\mathbf{P},\boldsymbol{\pi})$. Note that the total momentum of the coupled system is denoted as ${\,{}^{\boldsymbol{\pi}}{}}\mathbf{u}$. This Poisson bracket expresses kinematics of fluid mechanics advecting the polarization density with its conjugate momentum (relative momentum of the intrinsic dipole charges).

Note that bracket (2.5.2) may be briefly expressed as

[TABLE]

using the definition of $\{F,G\}^{\textrm{(SP)}}$ from (33) and $\{F,G\}^{\textrm{(SP)}_{\mathbf{A}}}$ from (2.5.2).

The evolution equations given by (2.5.2) are as follows,

[TABLE]

Note that the semi-direct product theory does not fully determine the relation of the total momentum ${\,{}^{\boldsymbol{\pi}}{}}\mathbf{u}$ and the mass momentum. We only know that the velocities coincide, which means that

[TABLE]

This relation can be ascertained by further assumptions on the physical properties of the interaction.

Let us assume ${\,{}^{\boldsymbol{\pi}}{}}\mathbf{u}=\mathbf{u}$, i.e., there appears no extra momentum on top of the mass momentum. The assumption corresponds with the nature of the transformation (129), where motion of an elementary dipole is described by two de-coupled canonical brackets, one for the center of mass and one for the relative quantities, and the center-of-mass momentum of the dipole affects neither $\mathbf{P}$ nor $\boldsymbol{\pi}$ directly.

2.5.3 Coupling to the electromagnetic field – total charge density

Having coupled fluid mechanics with polarization density and its conjugate momentum, let us finally make the coupling to the electromagnetic field, i.e. having state variables $(\rho_{\alpha},{\,{}^{\boldsymbol{\pi}}{}}\mathbf{u},s,\mathbf{D},\mathbf{B},\mathbf{P},\boldsymbol{\pi})$. Both pairs $(\boldsymbol{\pi},\mathbf{P})$ and $(\mathbf{D},\mathbf{B})$ were coupled to fluid mechanics by semidirect product. Advection of both pairs by fluid mechanics can be thus expressed (using  (29), (33), (2.5.2)) by Poisson bracket

[TABLE]

which, however, does not contain any direct coupling between $\mathbf{P}$ and $\mathbf{D}$. The coupling appears after the transformation to the field of electric induction $\boldsymbol{\mathcal{D}}$,

[TABLE]

for which it holds that

[TABLE]

where the second equality is due to (46) and (44). The divergence of $\boldsymbol{\mathcal{D}}$ represents the total charge density and is given as a sum of the bound charge density and the free charge density.

For a functional $\widehat{F}(\mathbf{D},\mathbf{P},\dots)=F(\boldsymbol{\mathcal{D}},\mathbf{P},\dots)$ then holds that

[TABLE]

Poisson bracket (51) can be reformulated using the electric induction (52) and the transformation (54) as

[TABLE]

Surprisingly, just the evolution equations of $\boldsymbol{\mathcal{D}}$ and of polarization momentum $\boldsymbol{\pi}$ are affected by the transformation. The transformed evolution equations of the electric induction and polarization momentum $\boldsymbol{\pi}$ then become

[TABLE]

Note that the Gauß’s law (53) is again compatible with the evolution of the $\operatorname{div}\boldsymbol{\mathcal{D}}$, given by (56a).

The evolution of the bound charge density, given by bracket (2.5.3) and formally identical to (49d), reads

[TABLE]

The divergent part of (57) evolves as a convected density which also be changed by $\boldsymbol{\pi}^{\dagger}$. The term $\boldsymbol{\pi}^{\dagger}$ can be interpreted as a reversible polarization current. This reversible current appears in the evolution of $\boldsymbol{\mathcal{D}}$, see Eq.  (56a).

The material relation for electric intensity now changes to

[TABLE]

which is the usual relation in the frame co-moving with matter. Transformation (130) brings these relation to the laboratory frame.

Total momentum.

In order to obtain the Lorentz force, the total momentum ${\,{}^{\boldsymbol{\pi}}{}}\mathbf{m}$ has to be related to the mass momentum $\mathbf{u}$. We choose

[TABLE]

since it is compatible with the choices made for the polarization-fluid coupling (50) and the EMHD coupling (34). The bracket (2.5.3) may be transformed using (59) and transformation rules analogous (39b) so that the description with the mass momentum $\mathbf{u}$ is obtained. The calculation is analogous to the transformation in the case of EMHD brackets (33) to (38). Apart from the presence of $(\mathbf{P},\boldsymbol{\pi})$ and its coupling to the fluid mechanics, the only actual difference is the $\{\cdot,\cdot\}^{\mathbf{P}\boldsymbol{\pi}}(\boldsymbol{\mathcal{D}},\boldsymbol{\pi})$ appearing in (2.5.3). Relation (59) implies the following form of the Lorenz force density (right hand side of the evolution equation for momentum $\mathbf{u}$)

[TABLE]

where $\mathbf{x}^{\widetilde{\dagger}}$ denotes derivative of energy dependent on $\mathbf{u}$ with respect to $\mathbf{x}$, analogously as in Section 2.4.2 transformation (39b).

Let us analyze (60). The first term contains the electric field $\mathbf{E}=\boldsymbol{\mathcal{D}}^{\widetilde{\dagger}}$ and $\mathbf{v}\times\mathbf{B}$ acting upon the total charge density $\operatorname{div}\boldsymbol{\mathcal{D}}$. The main difference to the Lorentz force appearing in (43), is the total charge density appearing in front of the bracket. The second term in (60) is a product of the magnetic field $\mathbf{B}$ and the reversible rate of change of polarization $\boldsymbol{\pi}^{\widetilde{\dagger}}$. Hence, it can be interpreted as a force exerted by the magnetic field due to the current of the polarized charge. Due to the form of (59) it holds $\boldsymbol{\pi}^{{\dagger}}=\boldsymbol{\pi}^{\widetilde{\dagger}}$ and $\mathbf{P}^{{\dagger}}=\mathbf{P}^{\widetilde{\dagger}}$.

Energy (42) turns to

[TABLE]

Note that for pure substances the tensor $\boldsymbol{\epsilon}$ in this formula for energy becomes the vacuum permittivity (multiplied by the inverse metric) because then the usual formula for dielectrics is recovered, see Sec. 3.3.2, after the polarization has relaxed to equilibrium. However, for a mixture one should consider a polarization for each constituent, and some of the polarizations can be relaxed while some not. We keep the tensor $\boldsymbol{\epsilon}$ in the formula for energy in order to cover also the mixtures.

In summary, reversible evolution equations for a mixture coupled with electromagnetic fields, polarization and its conjugate momentum were constructed in a geometric way (semidirect products).

2.5.4 Polarization waves and stress

Let us demonstrate the dynamics of the $\mathbf{P}$, $\boldsymbol{\pi}$ couple on the canonical polarization bracket (45), bringing new equations for propagation of polarization. For the sake of simplicity, we assume state variables $(\mathbf{P},\boldsymbol{\pi})$ in this Section. Let us assume that the energy contains a contribution corresponding to a dipole-dipole interaction that is proportional to the square of the polarization gradient $\nabla\mathbf{P}$. Hence, we have

[TABLE]

where $\alpha^{ij}$ and $\beta_{ij}$ are assumed to be symmetric positive definite constant tensors. The Hamiltonian evolution of $\mathbf{P}$ and $\boldsymbol{\pi}$ due to bracket (45) with energy (62) reads

[TABLE]

Assuming isotropy, $\alpha^{ij}=\alpha g^{ij}$ and $\beta_{ij}=\beta g_{ij}$, where $\alpha$ and $\beta$ are positive numbers and $g_{ij}$ the Eucledian metric, Eqs. (63) become

[TABLE]

which is a polarization wave equation. Speed of the polarization sound is $(\alpha\beta)^{-1/2}$.

Let us now consider the level of description with state variables $(\rho_{\alpha},s,{\,{}^{\boldsymbol{\pi}}{}}\mathbf{m},\mathcal{D},\mathbf{B},\mathbf{P},\boldsymbol{\pi})$, which is equipped with bracket (51) and where (62) contributes to the energy. The evolution equation for total momentum ${\,{}^{\boldsymbol{\pi}}{}}\mathbf{m}$ is given by the bracket (2.5.3) also determines the evolution of the mass momentum using transformation (59). Thus we have

[TABLE]

without specifying the energy dependence of the other variables. Note that we no longer pay attention the upper and lower indexes in this section for simplicity. The last line is the Lorentz force acting upon the total charge, c.f. (60). The terms with $\mathbf{P}$ and $\boldsymbol{\pi}$ on the second line of the right hand side of (65) can be written in a divergence form as

[TABLE]

which is the stress contribution due to the dipole-dipole interaction. Momentum is conserved.

The third term in (66) is not symmetric w.r.t. change of $i$ and $j$, but conservation of the angular momentum is not violated. Indeed, the cross product of the discussed term and the position vector $\mathbf{r}$ can be written in a divergence form. We have, omitting $\beta$,

[TABLE]

where the last term is simultaneously symmetric and antisymmetric in $i$ and $j$, hence it is zero. Similar observation was made in extra-mass-flux . The angular momentum is conserved for the proposed energy of the dipole-dipole interaction (62).

2.6 Magnetization

The magnetization of matter, see e.g. (purcell1985electricity, , Sec. 11), is due to the orientation of spins. Perhaps due to the resemblance with dynamics of rigid body rotation, the pioneering model of magnetization by Landau & Lifshitz LL1935 was based on that dynamics. In the following text we first recall the Hamiltonian formulation of rigid body dynamics and then we let the rigid body dynamics be advected by the fluid (using again the semidirect product theory).

2.6.1 Intrinsic dynamics of magnetization

The configuration manifold of rigid body rotations is the Lie group $SO(3)$. The standard machinery of differential geometry, e.g. arnold ; Fecko or  (pavelka2018multiscale, , Eq. 3.69), concludes that the Lie algebra dual, where angular momentum $\mathbf{M}$ seen from the body reference frame plays the role of state variable, is equipped with a Poisson bracket. Local version of the bracket (see also dv ) is

[TABLE]

where $\gamma$ is the gyromagnetic ratio and $\mathbf{M}(\mathbf{r})$ is the local angular momentum covector field. Bracket (68) is also called the spin bracket, see e.g. lakshmanan2011fascinating . The evolution equation implied by this bracket is (by the same procedure as in Sec. 2.2.1)

[TABLE]

$E$ being energy of the rotation. Derivative of energy with respect to $\mathbf{M}$ is the angular velocity $\boldsymbol{\omega}$.

2.6.2 Advection by fluid mechanics

Similarly as in the case of electromagnetic field or polarization, we will now construct the Poisson bracket expressing advection of magnetization and its dynamics by fluid mechanics (using again semidirect product). Now, however, the advected structure is not a cotangent bundle, but a Lie algebra dual, i.e. having state variables $(\rho,{\,{}^{\mathbf{M}}}\mathbf{m},s,\mathbf{M})$. The general formula for semidirect product of two Lie algebra duals, see  (esen2017hamiltonian, , Eq. 40), then gives Poisson bracket

[TABLE]

where ${\,{}^{\mathbf{M}}}\mathbf{m}$ is the new total momentum of the coupled system and $\mathbf{M}$ denotes the magnetization. Note that the advected electrodynamics are kept in the bracket for completeness. The right action of velocity field $F_{\mathbf{m}}$ on the vector field $G_{\mathbf{m}}$ is defined as negative of Lie derivative as usually,

[TABLE]

Using (71), bracket  (70) becomes

[TABLE]

which is the explicit form of Poisson bracket expressing kinematics of magnetization advected by fluid mechanics.

Let us choose

[TABLE]

which, as in the polarization case, means that there is no excess momentum and ${\,{}^{\mathbf{M}}}\mathbf{m}=\mathbf{m}$.

The evolution equation equations implied by (LABEL:LL03) are

[TABLE]

These evolution equations show how magnetization is advected by fluid mechanics, and how such advection affects the fluid motion itself. Moreover, magnetization keeps its intrinsic rigid-body-like dynamics. Finally, note that there is no explicit coupling to the electromagnetic field just as in the original LL1935 paper. The coupling is achieved implicitly later by letting energy depend on both $\mathbf{M}$ and $\mathbf{B}$.

2.6.3 Spin-spin interaction stress

Let us now suggest novel dynamics of magnetization advected by fluid mechanics. Let now $(\rho_{\alpha},\mathbf{m},s,\mathbf{D},\mathbf{B},\mathbf{M})$ be the considered level of description governed by (70). The spin-spin interaction contribution to the energy can, on the continuum level, be modeled by the gradient of magnetization $\nabla\mathbf{M}$ as in LL1935 ,

[TABLE]

where $\alpha^{ik}$ is a constant positive definite tensor. The mass momentum $\mathbf{u}$, given by the evolution of the total momentum (74a), relation (73) and relations analogous to (37), evolves as

[TABLE]

where only the derivatives of energy w.r.t. $\mathbf{M}$ are expressed explicitly. Note that we have put $\alpha^{ij}=\alpha g^{ij}$ and that we stopped distinguishing between upper and lower indexes for simplicity. The terms on the last line represent the Lorentz force acting on the free charge. The terms on the second line in (76) are the stress due to the spin-spin interaction. The terms can be, similarly to the dipole-dipole interaction (66), written in the divergence form as

[TABLE]

The non-symmetric term in (77) is equivalent to the one already discussed in the polarization case (67), therefore, it also does not violate the conservation of the angular momentum.

2.7 General level $(\rho_{\alpha},\widehat{\mathbf{m}},s,\mathcal{D},\mathbf{B},\mathbf{P},\boldsymbol{\pi},\mathbf{M})$

The hierarchy of the brackets built in the preceding paragraphs will be completed on the level of description containing all the discussed variables. The General bracket reads

[TABLE]

The first line is due to the dynamics of the CIT-mixture. The second line of (2.7) accounts for the electromagnetism and its coupling to continuum. The third line contains the polarization bracket and its coupling to continuum. The fourth line of (2.7) is due to the magnetization dynamics and its coupling to continuum. Brackets $\{\cdot,\cdot\}^{\textrm{(SP)}}$, $\{\cdot,\cdot\}^{\textrm{(SP)}_{\mathbf{A}}}$ and $\{\cdot,\cdot\}^{\textrm{(SP)}_{\mathbf{M}}}$ were found due to the semidirect product theory.

The evolution equations implied by the General bracket are

[TABLE]

When energy is among the state variables (thus transforming first to $\mathbf{u}=\rho\mathbf{v}$ to $\mathbf{v}$), the evolution equations read

[TABLE]

The bracket (2.7) can be projected to the levels of description upon which it was built. One can simply evaluate the bracket (2.7) on a set of functionals independent of a certain variables, see hierarchy . In the rest of this paper we enrich the reversible equations by irreversible terms in order to reduce this rather detailed description to the common continuum models coupling matter and electromagnetic field.

Consistently with the previously devised relations of the total momenta and mass momentum $\mathbf{u}$, we choose

[TABLE]

The Lorentz force is equal to the one given by (60).

3 Continuum thermodynamics and reductions

After having constructed a hierarchy of Poisson brackets for fluid mechanics of mixtures advecting electrodynamics, polarization and its conjugate momentum and magnetization, let us now enrich that detailed reversible dynamics by dissipative irreversible terms. This allows to see relaxation of fast mesoscopic variables and the effects on dynamics of less detailed variables. For instance we let the conjugate polarization momentum $\boldsymbol{\pi}$ relax to recover the standard Single Relaxation Time (SRT) model, which is widely used for comparison with experiments. We also let the magnetization $\mathbf{M}$ dissipate to recover the full Landau & Lifshitz model not only evolving in the laboratory frame, but being advected by the fluid. Finally, we approach the level of mechanical equilibrium, where neither momentum nor velocity is among the state variables. In particular, a generalized Nernst-Planck-Poisson equations are formulated on that level of description. In summary, a comprehensive multiscale thermodynamic construction of fluid mixtures equipped with electrodynamics, polarization and magnetization is provided.

3.1 Gradient dynamics

Before adding dissipative terms to the actual evolution equations, let us recall the general framework of gradient dynamics, where irreversible evolution is generated by derivatives of a dissipation potential Landau-Ginzburg ; gyarmati . Sound statistical arguments for gradient dynamics based on the large deviations principle was found in mielke2014 ; mielke-peletier ; Alberto-FDT . The following paragraph closely follows (pavelka2018multiscale, , Sec. 4.5, 4.6).

3.1.1 Dissipation potential

Consider a set of state variables $\mathbf{x}$ and let energy, entropy and mass of the system be denoted by $E(\mathbf{x})$, $S(\mathbf{x})$ and $M(\mathbf{x})$, respectively. A dissipation potential $\widetilde{\Xi}(\mathbf{x},\mathbf{x}^{*}):\mathbf{x}^{*}\rightarrow\mathbb{R}$ is a family of functionals of conjugate variables $\mathbf{x}^{*}$ parametrized by $\mathbf{x}$. We require every parametrization $\Xi[\mathbf{x}^{*}]=\widetilde{\Xi}(\mathbf{x})[\mathbf{x}^{*}]$ to satisfy:

Positiveness: $\Xi[\mathbf{x}^{*}]\geq 0$ and $\Xi[\mathbf{x}^{*}=0]=0$. 2. 2.

Monotonicity of derivative: $\left\langle\Xi_{\mathbf{x}^{*}},\mathbf{x}^{*}\right\rangle\geq 0\quad\forall\mathbf{x}^{*}$. 3. 3.

Convexity near $\mathbf{x}^{*}=0$. 4. 4.

Degeneracy with respect to mass $\left\langle\Xi_{\mathbf{x}^{*}},M_{\mathbf{x}}\right\rangle=0$. 5. 5.

Degeneracy with respect to energy $\left\langle\Xi_{\mathbf{x}^{*}},E_{\mathbf{x}}\right\rangle=0$. 6. 6.

$\Xi[\mathbf{x}^{*}]$ be even with respect to time-reversal transformation, see pavelka2014time .

The irreversible part of the evolution of a functional $F(\mathbf{x})$ is then given as

[TABLE]

Gradient dynamics automatically satisfies the second law of thermodynamics (growth of entropy in isolated systems). This is guaranteed for instance for convex dissipation potentials, but also non-convexity far from the origin (equilibrium) be taken into account nonconvex . Moreover, it is in close relation to the method of entropy production maximization Adam-EPM . Gradient dynamics plays a key role when formulating dissipation in the GENERIC framework.

Entropy-conjugate variables $\mathbf{x}^{*}$ are most suitable for formulation of the gradient dynamics. However, it is often more straightforward to use the energy-conjugate variables $\mathbf{x}^{\dagger}$, especially when combining the irreversible evolution with the Hamiltonian evolution (see e.g. be ; ADER ; ader-em ; godr ). Details on the transformation between the conjugate variables can be found in pavelka2018multiscale ; CMAT2018 .

A simple notorious example of the dissipation potential expressed in the energy-conjugate variables is

[TABLE]

This is a prototype of dissipation potential, since any general dissipation potential can be approximated by a quadratic one due to the convexity near equilibrium and flatness at equilibrium. For further details see (pavelka2018multiscale, , Sect. 4.5).

3.2 Dynamic maximum entropy principle

The principle of maximum entropy (MaxEnt), where unknown value of a variable is determined by finding the maximum value of entropy subject to constraints given by declared knowledge, has been successfully applied in many fields (information theory, thermodynamics, etc.) jaynes . However, in non-equilibrium thermodynamics the problem is not only to find value of a fast variable that has relaxed, but also to find the vector field along which the fast variable evolves, its evolution equation, when only less detailed variables are among the state variables (observables). To this end we recall the method of Dynamic MaxEnt (DynMaxEnt) grmela2013-cmat ; pavelka2018multiscale . Extensive treatment of the DynMaxEnt principle in the context of the continuum thermodynamics can be found in klika2019dynmaxent .

3.3 Relaxation of conjugate polarization momentum $\boldsymbol{\pi}$

In subsection 2.5 polarization was equipped with its conjugate momentum $\boldsymbol{\pi}$. Inspired by the relaxation of the damped particle in subsection 3.2, we shall let the conjugate momentum relax to recover dissipative evolution of polarization. Let us choose state variables $(\rho_{\alpha},\mathbf{m},s,\boldsymbol{\mathcal{D}},\mathbf{B},\mathbf{P},\boldsymbol{\pi})$, represented in (2.5.3). The energy depends on all the variables unless assumed otherwise.

3.3.1 Polarization relaxation via $\boldsymbol{\pi}$

Let us assume a dissipation potential quadratic in $\boldsymbol{\pi}^{*}$,

[TABLE]

where $\Lambda_{(\boldsymbol{\pi})}$ is a positive definite symmetric constant tensor. Assuming also energy quadratic in $\boldsymbol{\pi}$, the MaxEnt value of $\boldsymbol{\pi}$ is zero. Using (56b), the DynMaxEnt relaxation of $\boldsymbol{\pi}$ can be formulated as

[TABLE]

which is the constitutive relation to be plugged into the remaining evolution equations. In particular, evolution equation for polarization (79f) becomes

[TABLE]

which is an equation for polarization equipped with both reversible and irreversible terms.

Note that the irreversible terms on the right hand side of equation (86) are generated by dissipation potential $\Xi^{(\boldsymbol{\pi})}$ evaluated at constitutive relation (85),

[TABLE]

Similar ideas were presented in entrpo . This is the dissipation potential generating irreversible evolution of $\mathbf{P}$ and $\boldsymbol{\mathcal{D}}$. We may also pursue further relaxation of the polarization. If DynMaxEnt value of $\mathbf{P}$ is zero, e.g. the total energy is quadratic in polarization, all the momentum coupling terms from (86) vanish and we obtain

[TABLE]

as the ultimate value of $\mathbf{P}^{\dagger}$.

3.3.2 Linear isotropic dielectric

Let us now demonstrate the results of the previous paragraph on a particular choice of energy (61) restricted to state variables $(\rho_{\alpha},\mathbf{v},s,\boldsymbol{\mathcal{D}},\mathbf{B},\mathbf{P})$. Assuming that the internal energy be independent of $\mathbf{P}$, the MaxEnt value of $\mathbf{P}$ is zero so that (88) holds. We are going to determine the energy on the level where $\mathbf{P}$ is relaxed. Since the following considerations are simpler with velocity $\mathbf{v}$ instead of momentum $\widehat{\mathbf{m}}$, transformation (39b) is carried out. Recalling (88), we have

[TABLE]

Let us first observe that the relation (89)right can be also reformulated, using (37b), as

[TABLE]

Energy (61) then becomes

[TABLE]

The relations (52) and (90) can be combined to obtain,

[TABLE]

Finally, for either negligible magnetic field $\mathbf{B}=0$ or negligible velocity $\mathbf{v}=0$ we obtain that

[TABLE]

which leads to the usual form of Coulomb’s law for linear isotropic dielectrics purcell1985electricity ,

[TABLE]

since then $\boldsymbol{\mathcal{D}}^{\dagger}=\boldsymbol{\mathcal{D}}^{\widetilde{\dagger}}=\boldsymbol{\epsilon}_{0}\mathbf{E}$.

3.3.3 Single Relaxation Time model

Consider again equation (86) and assume further that $\mathbf{v}=0$ and $\mathbf{B}^{\dagger}$=0 and that the internal energy is negligible. Then, for energy

[TABLE]

the evolution of polarization (86) becomes

[TABLE]

which represents a dissipative evolution of polarization subject to electric field.

Let us first analyze the evolution equation by applying harmonic electric field $\mathbf{E}={\mathbf{E}}_{0}\exp(i\omega t)$. Equation (96) then gives the Single Relaxation Time (SRT) model of polarization, see e.g. boettcher1979dielectric2 ,

[TABLE]

provided that $\mathbf{P}=\mathbf{P}_{0}\exp(i\omega t)$.

In summary, by letting the conjugate polarization momentum $\boldsymbol{\pi}$ relax, a dissipative evolution of polarization is obtained (86). If the mechanical equilibrium is further assumed, this dissipative evolution is compatible with the SRT model widely used for comparison to experiments. Finally, the equilibrium of the dissipative evolution is the linear relation between polarization and electric intensity known from electrostatics.

3.4 Relaxation of magnetization

Let us now discuss relaxation of magnetization $\mathbf{M}$ inspired by the Landau & Lifshitz model LL1935 . We consider dissipation potential

[TABLE]

where $\Lambda_{(\mathbf{M})}$ is a positive definite symmetric tensor field. Let us, for simplicity, assume that $\Lambda_{(\mathbf{M})}$ is a constant multiple of the constant Eucledian metric and let us not distinguish between the upper and lower indexes in this section. The derivative of (98) w.r.t $\mathbf{M}^{\dagger}$ is

[TABLE]

Hence, the irreversible evolution of $\mathbf{M}$ due to (98) is

[TABLE]

which is compatible with the Landau & Lifshitz model of magnetization once suitable energy is provided,  (brown1963micromagnetics, , Sec. 3.7).

Having recovered the Landau & Lifshitz model, let us also formulate its generalized version advected by fluid mechanics and interacting with magnetic field. Let us assume the following energy

[TABLE]

The energy-conjugate to the magnetization $\mathbf{M}$ then reads

[TABLE]

The field of magnetic intensity, which is conjugate to $\mathbf{B}$, satisfies

[TABLE]

which is the usual relation between $\mathbf{H}$, $\mathbf{B}$ and $\mathbf{M}$ in the frame co-moving with matter. It holds, moreover, that

[TABLE]

Combining (80h) and (100), we obtain the evolution of magnetization, using also (102),

[TABLE]

Equation (105), supplied with the rest of Eqs. (79f), is the generalized Landau-Lifshitz magnetization relaxation model333To fully recover the model, one should also add a weakly non-local term $\frac{1}{2}\alpha(\nabla\mathbf{M})^{2}$ to the energy density, see Sec. 2.6., where magnetization relaxes, interact with electromagnetic field and where it is advected by the fluid.

3.4.1 Linear isotropic magnetizable medium

Let us demonstrate some properties of energy $E^{(\mathbf{M})}$ in the context of MaxEnt. Introducing the MaxEnt value of the magnetization,

[TABLE]

into the energy $E^{(\mathbf{M})}$ gives

[TABLE]

for which it holds that

[TABLE]

In particular, $\mathbf{M}^{\dagger}=0$ we obtain $\mathbf{M}=\chi_{m}\mathbf{H}$. The MaxEnt value of the magnetization is used to recover the usual relation between the magnetic field and its intensity for linear isotropic materials.

3.5 Electro-diffusion – dissipation of $\boldsymbol{\mathcal{D}}$ and $\rho_{\alpha}$

Let us now consider only state variables $(\rho_{\alpha},\boldsymbol{\mathcal{D}})$, which are essential for electrodiffusion. Poisson bracket (2.7) for mixtures may be endowed with a weakly non-local electro-diffusion dissipation potential expressing friction between components of the mixture and the zero-th species (e.g. solvent),

[TABLE]

where $M_{\alpha\beta}$ is a symmetric, positive definite matrix of binary diffusion coefficients. Note that $\nabla\rho^{*}_{\alpha}-e_{0}\frac{z_{\alpha}}{m_{\alpha}}\mathcal{D}_{i}^{*}$ is proportional to the gradient of electrochemical potential of species $\alpha$. This dissipation potential conserves both mass and energy.

The irreversible part of the evolution equations can be expressed in the energy-conjugate variables as

[TABLE]

Note that the divergence of (110c) is equal to the sum of (110b) and (110a) weighted by the charge per mass of the respective species. Therefore, the irreversible evolution given by (109) is compatible with Gauß’s law given by (44). In other words

[TABLE]

Equations (110) express irreversible evolution of densities of the constituents and electric displacement under mutual interaction. The right hand sides of the equations for densities can be regarded as gradients of respective electrochemical potentials. Evolution for the $\boldsymbol{\mathcal{D}}$ field leads to relaxation of the field consistent with the Gauß law and with the Poisson equation.

Note that a dissipation potential introducing dissipation of the partial mass densities (i.e. containing $\rho_{\alpha}^{*}$) is required to contain corresponding terms with $\boldsymbol{\mathcal{D}}^{*}$ otherwise the irreversible evolution would not be compatible with the Gauß law, c.f. (111).Therefore, if the validity of Gauß’s law for the free charge is required, then the form of dissipation potential involving $\boldsymbol{\mathcal{D}}^{\dagger}$ is constrained.

On short enough distances magnetic field effects are usually negligible compared with effects of the electric field, in contrast to long distances, where electric field usually does not play any relevant role due to screening. Let us analyze the former case, paving the way towards electrochemical problems. The electro-diffusion due to (109) introduces dissipative fluxes identical to those of the generalized Planck-Nernst-Poisson systems (gPNP) presented in dreyer2018bulksurface . Assuming that no magnetic field is present and that its evolution equation (80e) is satisfied, it follows that

[TABLE]

see also (92). The considered level of description consists of $(\rho_{\alpha},\mathbf{u},s,\varphi)$. Let us consider a further relaxation of the mass momentum $\mathbf{u}$. We can proceed in two ways. Either we suppose for that energy does not depend on the mass momentum and the barycentric velocity vanishes, or we can think of a dissipation of the mass momentum, e.g. viscosity, and use the DynMaxEnt principle in order to determine the value of the barycentric velocity $\mathbf{u}^{\dagger}$. Let us now follow the latter route. The energy is quadratic in $\mathbf{u}$, see (92), hence, its MaxEnt value is zero, leaving the remainder of the momentum equation as an constitutive equation for the velocity field. For the dissipation potential generating the irreversible part of the Navier-Stokes equations, see e.g. (pavelka2018multiscale, , Eqns. 4.74, 4.76). The viscous dissipative terms are then added to the reversible balance of the mass momentum. The whole gPNP-Stokes systems eventually reads

[TABLE]

Note that the third equation is not violated if the $curl$-operator is applied on it.

4 Discussion

W. Dreyer, C. Guhlke and R. Müller in dreyer2018bulksurface published a comprehensive analysis of fluid mixtures coupled with electromagnetic fields, including polarization and magnetization which will be further referred as the DGM approach. Their treatment of surfaces, c.f. ondra-surface , as independent thermodynamic systems interacting with the bulk, being beyond the scope of the presented work, have elucidated many electrochemical problems using non-equilibrium thermodynamics, for instance a unified theory of the Helmholtz and Stern layers, a derivation of Butler-Volmer equations, or useful asymptotic techniques, see fuhrmann2016numerical ; dreyer2016new ; dreyer2018bulksurface ; dreyer2014mixture ; dreyer2014mixture ; guhlke2015theorie . Since our goal is in close relation to that works, let us compare the two approaches in detail.

An important conceptual difference between the two approaches lies in the treatment of the state variables, i.e., the levels of the description. In our approach, the level of description is always defined first by declaring the set of state variables. Energy and entropy of the system can then, in principle, depend on all the state variables. Each state variable has its evolution given by the GENERIC equation (1), where the reversible and the irreversible parts of the evolution are separated. In contrast, the DGM approach develops a system of general balance equations which outlines the relations between the physical quantities and their fluxes and productions. The dependence of the fluxes and production upon the variables is determined using the entropy principle.

The DGM approach is formulated in an inertial frame of reference for which holds that (i) mass center not subjected to external forces moves with a constant velocity and (ii) the Lorentz-Maxwell-aether (or material) relations 41 are valid. This is similar, although not equivalent, to our construction. The choice of the energy in DGM,

[TABLE]

is analogous to the energy 42. Moreover, the formulation of the Maxwell equation in DGM is, in vacuum, formally equivalent to our approach.

In the non-vacuum case, the total charge is for both formulations given by $\operatorname{div}\boldsymbol{\mathcal{D}}$. The $\operatorname{div}\mathbf{P}$ also bears the same meaning, i.e. the density of the bounded charge. The DGM evolution of the polarization $\mathbf{P}$ is given as formal solution to the bound charge balance. It reads

[TABLE]

where $\mathbf{M}^{\text{D}}$ is the Lorentz magnetization and $J^{\textrm{P}}$ denotes the non-convective part of the polarization current, cf. equation (79f). The DGM formulation contains no balance equation concerning $\mathbf{M}^{\text{D}}$, an evolution equation is later on found for $\widehat{\mathbf{M}}^{\text{D}}=\mathbf{M}^{\text{D}}+\mathbf{v}\times\mathbf{P}$ as a consequence of the closure. In contrast, the polarization bracket (45) contains the polarization momentum $\boldsymbol{\pi}$ as a variable and therefore implies its reversible evolution. The polarization momentum $\boldsymbol{\pi}$ can be projected to a reduced variable $\boldsymbol{\mathcal{M}}=\operatorname{curl}\boldsymbol{\pi}$, would lead to the appearance of $\operatorname{curl}\boldsymbol{\mathcal{M}}^{\dagger}$ in the equation for polarization. This is perhaps the closest the two approaches can get in this respect.

In DGM the balance equations for the densities of the partial masses $\rho_{\alpha}$, mass momentum $\mathbf{u}$, total momentum $\mathbf{u}+\boldsymbol{\mathcal{D}}\times\mathbf{B}$, total energy $e^{\text{D}}$, internal energy $\epsilon^{\text{D}}$, total charge $\rho_{\text{f}}-\operatorname{div}\mathbf{P}$, bound charge $-\operatorname{div}\mathbf{P}$, flux of the magnetic field, and the entropy $s^{\text{D}}$ are formulated, see dreyer2018bulksurface , which is translated in our approach to the respective evolution equations. In particular, the total energy and the total momentum $\mathbf{u}+\boldsymbol{\mathcal{D}}\times\mathbf{B}$ are, in the absence of external forces, conserved quantities in DGM as in our approach. Coupling between the charged fluid and electromagnetic field is given by the choice of Lorentz’s force (dreyer2018bulksurface, , Eqn. 36a):

[TABLE]

as a source term in the mass momentum balance and by the choice of Joule heating

[TABLE]

as a source term in the internal energy balance, cf. Appendix D. Symbols $J^{\mathrm{F}}$ denotes the non-convective part of the free charge current. This is similar to our formula (60), at least in the absence of irreversible terms.

As it is shown in Appendix D, the presented coupling between the motion of mass and the electromagnetic field can be formulated using a similar choice of Joule’s heating, see (138), and restricts the form of the total momentum to (142). The dissipation potentials for the irreversible, conductive, currents are formulated in terms of total momentum $\widehat{\mathbf{m}}$, so that $\boldsymbol{\mathcal{D}}^{\dagger}$ is the co-moving electric field. Hence, the Lorentz force acting upon the conductive currents can be found when the dissipation potentials are transformed into the variables with the mass momentum $\mathbf{u}$.

The choice of Lorentz’s force (119) tells that the electric field acts upon the total charge and the magnetic field upon convective current density, $\mathbf{J}^{\mathrm{F}}$ and $\mathbf{J}^{\text{P}}$ within the DGM theory. This is equivalent for the total charge and the conductive free charge $\mathbf{J}^{\mathrm{F}}$. The difference appears for $\mathbf{J}^{\text{P}}$, since it is not purely irreversible and contains reversible terms in the sense of the time-reversal transform.

Our evolution equation for polarization consists of reversible and irreversible parts. The former is given by Eq. (80f), and the terms of the right hand side containing velocity be rewritten as the Lie derivative of $\mathbf{P}$, $\mathcal{L}_{\mathbf{v}}\mathbf{P}$. Polarization is thus simply advected by the fluid (apart from interacting with the electromagnetic fields and apart from relaxation processes). On the other hand, polarization is not simply advected in DGM as the equation for $\mathbf{P}$ contains other velocity-dependent terms therein. This is a difference between our approach and DGM.

Authors of DGM assumed that the entropy density $s^{\text{D}}$ depends only on a specific subset of the variables

[TABLE]

thus obtaining a specific definition of the temperature $T^{\text{D}}$, see (dreyer2018bulksurface, , Eqn. 48a, 49a)444the dependence on the bulk deformation gradient is omitted. This choice allowed them to find a reasonably simple closure of the equations using the entropy principle, c.f. jou2010extended . On the other hand, temperature is in our approach defined as derivative of the internal energy density with respect to entropy as is usual degroot1984nonequilibrium .

Let us now focus on some features of the closure, especially, the non-convective flux of the bound charge that. They derived the following evolution of the polarization $\mathbf{P}$ and magnetization $\widehat{\mathbf{M}}^{\text{D}}$,

[TABLE]

The relaxation of the polarization, i.e. the last bracket in (122a), is equivalent to the dissipation derived using the DynMaxEnt principle given by potential $\Xi^{\mathbf{P}}$, see (87). The conductive flux of free charge $\mathbf{J}^{\textrm{F}}$ is equivalent to the one generated by the dissipation potential $\Xi^{\textrm{D}}$ in (110c), which also holds for the mass diffusion fluxes.

When the relaxation time of polarization $\tau^{\textrm{P}}$ can be considered large, the dissipative part in (122a) can be neglected. Eventually, the reversible part of the bound charge evolution in a volume $V$ for the two theories read

[TABLE]

The net bound charge in a fixed volume $V$ can thus be reversibly changed by convection or by change of $\boldsymbol{\pi}^{\dagger}$. In contrast to DGM, the bound charge is advected like a scalar field in our approach, which corresponds with that polarization is Lie-dragged in our approach.

Also the evolution of magnetization in DGM, Eq. (122b), is different from our equation (80h). In our case the evolution of magnetization caused by fluid motion is again just the Lie drag. Apart from that, there is a contribution from inertia of the magnetization itself, as in the Landau-Lifshitz model, which brings interaction with the magnetic field when energy depends on both magnetization and magnetic field.

Let us consider the volume and shear viscosity to be vanishing, the stress tensor presented in (dreyer2018bulksurface, , Eqn 62) reads555omitting also the dependence on the volume-preserving deformation gradient

[TABLE]

here $\widehat{\mathbf{E}}=\mathbf{E}+\mathbf{v}\times\mathbf{B}$. The dependence of the stress tensor $\sigma$ on the material properties linked to the polarization $\mathbf{P}$ and magnetization $\widehat{\mathbf{M}}^{\text{D}}$ is provided by the isotropic part $p\delta_{ij}$ the rest of the tensor is linear in $\mathbf{P}$ and $\widehat{\mathbf{M}}^{\text{D}}$, independent of the choice of the energy. An energy weakly non-local in polarization or magnetization, see e.g. (62) and (75), leads in the here presented treatment to a structurally different stress since a non-symmetric components can appear in the evolution equation for the mass momentum $\mathbf{u}$, see (66) or (77). The stress in DGM is different from our formulas for stress, mainly in the off-diagonal part.

The Onsager-Casimir reciprocal relations Onsager1931 ; casimir1945onsager ; degroot1984nonequilibrium are one of the corner-stones of non-equilibrium thermodynamics. Roughly speaking, they say that state variables with the same parity are coupled through an operator symmetric with respect to the simultaneous transposition and time-reversal while variables with opposite parities are coupled through an antisymmetric operator. They are automatically satisfied within the GENERIC framework in a generalized sense (beyond near equilibrium) ottinger2005minimal ; pavelka2018multiscale ; pavelka2014time . In dreyer2018bulksurface they seem to be satisfied as well, but parities of the state variables are determined from the power of seconds in the units of the variables instead of the time-reversal transformation, which is not invariant with respect to changing the physical units. Although the results seem to be all right, one should be careful in principle, perhaps using a precise definition of affine and vector spaces Matolcsi ; tamas-kinematics or the definition based on projections from more detailed levels pavelka2018multiscale .

5 Summary and Conclusion

In the first section a hierarchy of Poisson brackets describing the reversible dynamics of a charged, polarized and magnetized continua coupled with electromagnetic field has been developed by means of differential geometry, see Figure 1. The semidirect product of the fluid mechanics Lie Algebra dual $(\rho,\mathbf{u},s)$ and the electromagnetism cotangent bundle $(\mathbf{A},-\mathbf{D})$ results in the reversible electro-magneto-hydrodynamics already known marsden1982hamiltonian ; holm1986hamiltonian ; esen2017hamiltonian ; pavelka2018multiscale . Newly, cotangent bundle $(\boldsymbol{\pi},-\mathbf{P})$ describing dynamics of bound charge is also coupled to $(\rho,\mathbf{u},s,\mathbf{D},\mathbf{B})$ using the same technique. Finally, the classical spin dynamics represented by the local Lie algebra dual of $\text{SO}(3)$ is coupled to the Lie algebra dual of fluid mechanics, which can be seen as advection of the Landau-Lifshitz model by fluid mechanics. The theory is the Galilean invariant.

The second section is dedicated to introduction of irreversible dynamics using dissipation potentials and subsequent reductions of the before-built levels of description to the less detailed levels. The reduction is carried out by the Dynamic Maximization of Entropy (DynMaxEnt) technique, which forms passage from finer to rougher levels of description in a geometric way. In particular, the polarization momentum $\boldsymbol{\pi}$ is relaxed giving rise to dissipation of the polarization field itself. Further exploitation of the induced dissipation of $\mathbf{P}$ leads to the standard formulas for linear dielectrics and to the Single Relaxation Time model. The dissipation potential for the magnetization $\mathbf{M}$ is found so that the Landau-Lifshitz model of spin relaxation is restored. Finally, the electro-diffusion dissipation potential is introduced, leading to a generalized Poisson-Nernst-Planck-Stokes model.

In this rather complex dynamics, we observe phenomena like advection of polarization and magnetization by the fluid, polarization and magnetization waves and the respective inertial effects. We leave concrete applications of these phenomena for future researchers.

6 Acknowledgment

Authors thank to professor František Maršík for his relentless encouragement to explore the relation of electromagnetism and non-equilibrium mixture theories on a fundamental level. Authors thank also to Clemens Guhlke, Rüdiger Müller and Wolfgang Dreyer for their utter patience and openness during the elucidating discussions on the topic. We are also grateful to Péter Ván and Tamás Fülöp for consulting response to a reviewer regarding electrodynamics and space-time. PV and MP are indebted to Miroslav Grmela for teaching us to geometric thermodynamics.

This work was supported by the German Research Foundation, DFG project no. FU 31614-1. This work was supported by the Czech Science Foundation, project no. 17-15498Y. This work was supported by Charles University Research program No. UNCE/SCI/023.

Appendix A Notation

Appendix B Elementary dipole

Two classical charged mass points are described by their positions, $\mathbf{r}^{1}$ and $\mathbf{r}^{2}$, momenta $\mathbf{p}^{1}$ and $\mathbf{p}^{2}$, carrying charge $q_{1}$ and $q_{2}$, respectively. The dynamics of the particles is governed by the canonical Poisson bracket:

[TABLE]

Assume that $q=q_{1}=-q_{2}$ and consider the following transformation of the variables:

[TABLE]

Bracket (126) then transforms into:

[TABLE]

Appendix C Galilean transformations

Maxwell equations (30) are not in the usual form, e.g. from Landau2 . Usually the equations are derived in the special relativistic setting by first forming the electromagnetic four-cotensor from fields $\mathbf{E}$ and $\mathbf{B}$, and the resulting evolution equations are for $\mathbf{E}$ and $\mathbf{B}$ as well.

The Maxwell equations are naturally Lorentz invariant, as follows from the special relativistic treatment. However, it is often the case that Galilean invariance is sufficient for describing the observed phenomena. Moreover, Galilean transformations can be seen as low-velocity and low-strength limits of Lorentz transformations. Since Maxwell equations adequately describe usual phenomena around us, what is their behavior with respect to the Galilean transformations?

When having the pair $(\mathbf{E},\mathbf{B})$ at hand, as in the usual case, one has to be restricted either to the electric limit (magnetic effects negligible) or magnetic limit (electric effects negligible), as discussed in leblond1973galilean . However, when we have all the fields $(\mathbf{E},\mathbf{B},\mathbf{D},\mathbf{H})$, which is the case of Eqs. (80), the equations are invariant with respect to Galilean transformations

[TABLE]

The reasoning can be found for instance in Matolcsi or leblond1973galilean , where four-vectors and four-covectors are split into the time-like and space-like parts, and the Galilean-invariant Maxwell equations then result from such splitting of the special relativistic four-dimensional structure in the low velocity limit, c.f. Matolcsi-Van .

Maxwell equations in vacuum (30) are not Galilean invariant. However, equations (43), where Maxwell equations are complemented by motion of matter, are Galilean invariant, as noted in leblond1973galilean , albeit without the evolution equations for matter. In leblond1973galilean they found material relations (41) not being Galilean invariant, since they had to be formulated in a special reference frame (e.g. aether). However, we have a natural reference frame at hand (the frame co-moving with the center of mass) where the material relations are valid. This makes the theory fully Galilean invariant.

Appendix D Relation of the mass momentum and the total momentum

Let us define the mass momentum as

[TABLE]

for some given energy $E(\rho,\mathbf{m},s,\boldsymbol{\mathcal{D}},\mathbf{B},\mathbf{P},\boldsymbol{\mu},\mathbf{M})$. This relation thus defines a function $\widetilde{\mathbf{u}}$ such that

[TABLE]

Let us assume that there exists a different representation of energy $\widetilde{E}(\rho,\mathbf{u},s,\boldsymbol{\mathcal{D}},\mathbf{B},\mathbf{P},\boldsymbol{\mu},\mathbf{M})$. However, the energy of the system is invariant w.r.t. the choice of representation. Therefore,

[TABLE]

We require the velocity field to be also invariant of the chosen representation,

[TABLE]

As a consequence of the chain rule and formula (131), we obtain more specific form of $\widetilde{\mathbf{u}}$, that is,

[TABLE]

the extra momentum $\mathbf{f}$ is also assumed to be independent of mass and entropy density and being caused only due to the electromagnetic field.

The form of (135) suggest that the extra momentum should be the momentum of the electromagnetic field. Since the theory of the Maxwell equations is built upon the assumption of the field superposition, we assume $\mathbf{f}$ to be a bilinear function of $\boldsymbol{\mathcal{D}}$ and $\mathbf{B}$, hence we have,

[TABLE]

Using the chain rule, formulas (133) and (136), we find the relations between the derivatives of the energy in the two representations,

[TABLE]

The kinetic energy of the system may be defined using either the velocity or the mass momentum. In general, kinetic energy is not a conserved quantity and may be changed during evolution of the system. Experimental evidence suggests that the change of the kinetic energy of a charged matter, due to the presence of the electromagnetic field, is linearly proportional to the charge density and the electric field, i.e., $\partial_{i}\mathcal{D}^{i}\mathcal{D}^{\dagger}_{j}m^{\dagger j}$. We require this power density to be independent of the chosen representation, i.e., $\partial_{i}\mathcal{D}^{i}\mathcal{D}^{\dagger}_{j}m^{\dagger j}=\partial_{i}\mathcal{D}^{i}\mathcal{D}^{\widetilde{\dagger}}_{j}v^{j}$. Using (137) we have

[TABLE]

To this end, the extra term on the right-hand side must be zero, i.e., it must vanish for arbitrary magnetic and velocity fields. For example, it must hold

[TABLE]

Apparently, the tensor $a_{ijk}$ is antisymmetric and isotropic. Up to a sign, there exists only one such tensor in three dimensions,

[TABLE]

The sign choice is due to the experimentally observed difference in bending of rays of negatively and positively charged particles LL1935 . Let us the calculate time change of the kinetic energy,

[TABLE]

The terms on lines (141a), (141b), (141c) are the changes of the kinetic energy due to the fluid mechanics, dipole mechanics and spin mechanics, respectively.

The terms on line (141d) are the change of the kinetic energy due to the power of the electromagnetic field. In total, the change of kinetic energy is independent of the chosen energy representation, see (137) with the choice of (140).

In summary, the relation between the mass momentum $\mathbf{u}$ and the total momentum $\mathbf{m}$

[TABLE]

is recovered.