Relations de réciprocité d'Onsager
# Relations de réciprocité d'Onsager

En thermodynamique hors équilibre, les relations de réciprocité d'Onsager ou relations de réciprocité d'Onsager-Casimir caractérisent les coefficients phénoménologiques qui apparaissent dans les relations qui relient les flux de variables extensives caractérisant le système considéré aux affinités thermodynamiques correspondantes.

Elles ont été établies par Lars $Onsager^,$ en 1931 et précisées par Hendrik Casimir en 1945.

Comme d'autres relations de ce type, par exemple le principe de Curie, elles sont l'expression des propriétés d'invariance ou de symétrie de ces systèmes. Dans ce cas précis, il s'agit de l'invariance par renversement du temps.

## Exemple

Les exemples de relations flux-affinités sont très nombreux dans la physique. Ces relations peuvent être directes comme la relation flux de chaleur - gradient de température dans un milieu quelconque (conduction thermique, associée à la loi de Fourier), flux de masse - gradient de concentration des espèces dans un fluide (diffusion, associée à la loi de Fick). Mais il peut aussi y avoir des relations croisées flux de chaleur - gradient de concentration de masse (effet Dufour) et flux de masse - gradient de température (effet Soret).

Le problème est de savoir si tous les coefficients phénoménologiques qui apparaissent dans ces équations sont ou non liés entre eux.

## Potentiels thermodynamiques, affinités et flux

### Dans l'exemple précédent

Le potentiel thermodynamique fondamental est l'énergie interne $u$. Dans un système fluide, $u$ dépend de la densité de matière $r$ et de la densité d'entropie $S$ de la façon suivante :

$\mathrm {d} u=T\,\mathrm {d} S+m\,\mathrm {d} r$

où $T$ est la température et $m$ le potentiel chimique. Nous pouvons écrire :

$\mathrm {d} S={1 \over T}\,\mathrm {d} u-{m \over T}\,\mathrm {d} r$

les quantités extensives $u$ et $r$ sont conservées et leurs flux satisfont aux équations de conservation :

${\begin{array}{rcl}\partial _{t}u+\nabla .\,\mathbf {J} _{u}&=&0\\[0.6em]\partial _{t}r+\nabla .\,\mathbf {J} _{r}&=&0\end{array}}$

où $\partial _{t}$ indique la dérivée partielle par rapport à la variable t et $\nabla \cdot \mathbf {J}$ la divergence du flux de densité $\mathbf {J}$.

Les gradients des variables conjuguées de $u$ et de $r$, qui sont $1/T$ et $-m/T$ sont les affinités thermodynamiques (en anglais thermodynamic forces) et elles provoquent l'apparition de flux des variables extensives correspondantes, mais également des termes croisés dans la relation entre affinités et flux :

${\begin{array}{rcl}J_{u}&=&L_{uu}\nabla \left({1 \over T}\right)+L_{ur}\nabla \left(-{m \over T}\right)\\[0.6em]J_{r}&=&L_{ru}\nabla \left({1 \over T}\right)+L_{rr}\nabla \left(-{m \over T}\right)\end{array}}$

La relation de réciprocité de Onsager établit la relation existante entre les coefficients cinétiques $L_{ur}$ et $L_{ru}$.

### Formulation générale

Par la suite on adopte la convention de sommation d'Einstein.

Soit $E_{i}$ un jeu de variables extensives dont dépend une entropie $S$. Dans l'analyse suivante, ces symboles réfèreront à la densité de ces quantités thermodynamiques. Alors :

    $\mathrm {d} S=I_{i}\,\mathrm {d} E_{i}\,,\quad I_{i}:={\frac {\partial S}{\partial {E_{i}}}}$

On définit ainsi les quantités intensives $I_{i}$ conjuguées des quantités extensives $E_{i}$. Les gradients de ces quantités sont les affinités :

    $\mathbf {F} _{i}=\nabla I_{i}$

Les flux $\mathbf {J} _{i}$ et les variables extensives $E_{i}$ sont liés par les équations de continuité :

    $\partial _{t}E_{i}+\nabla \cdot \mathbf {J} _{i}=0$

L'entropie étant une fonction régulière, on peut écrire pour des petits écarts à l'équilibre thermodynamique :

    $\mathbf {J} _{i}=L_{ij}\,\mathbf {F} _{j}$

La production d'entropie s'écrit :

    $\sigma =L_{jk}\,\mathbf {F} _{j}\cdot \mathbf {F} _{k}\geq 0$

On en déduit que la matrice des $L_{jk}$ est définie positive et que les termes diagonaux $L_{jj}$ sont positifs ou nuls.

## Relations d'Onsager-Casimir

La relation la plus générale est obtenue dans le cadre de la physique statistique hors d'équilibre et s'exprime sous la forme suivante :

* Le coefficient cinétique $L_{jk}$ est tel que, en présence d’un champ magnétique B et d’une vitesse angulaire ω, le coefficient cinétique traduisant l’influence de l'affinité $F_j$ sur le flux $J_k$ est relié à celui qui traduit l’influence de l'affinité $F_k$ sur le flux $J_j$ par :
* $L_{kj}(\mathbf {B} ,{\boldsymbol {\omega }})=\epsilon _{j}\epsilon _{k}L_{jk}(-\mathbf {B} ,-{\boldsymbol {\omega }})$
* où $\epsilon =\pm 1$ est la signature de la variable correspondante dans un retournement du temps : ε=1 si la variable est inchangée par l'opération $t\mapsto -t$, ε=-1 sinon.

La démonstration des relations d'Onsager repose sur l'analyse des fluctuations cinétiques du milieu. Elle sera faite dans l'hypothèse où il n'y a ni rotation, ni champ magnétique.

### Les coefficients cinétiques

Soit $x_i$ l'écart à l'équilibre du aux fluctuations statistiques des N variables extensives $E_i$ et $dS(x_i)$ l'écart à l'équilibre correspondant $S(E_i)$ pour l'entropie. Cet écart est supposé suffisamment faible pour que l'on puisse linéariser les équations. La probabilité d'un ensemble d'écarts ${x_i}$ est proportionnelle au nombre de micro-états, donné par la loi de Boltzmann :

$p(x_{i})=Ce^{\frac {S}{k}}$

S est une fonction positive, concave, continûment différentiable et l'équilibre est défini par $\left.{\frac {\partial S}{\partial E_{i}}}\right|_{S=S_{0}}=0\,,\;\;\;\forall i$

On peut donc écrire un développement en série de Taylor conduisant à une distribution gaussienne :

$p\simeq Ce^{\frac {S_{0}}{k}}e^{-{\frac {1}{2}}\beta _{ik}x_{i}x_{k}}\,,\;\;\;\beta _{ik}=-{\frac {1}{k}}{\frac {\partial ^{2}S}{\partial x_{i}\partial x_{k}}}$

la matrice $β=(β_{ik})$ est symétrique définie positive et la normalisation de p permet d'obtenir la constante :

$\int _{-\infty }^{\infty }...\int _{-\infty }^{\infty }p\,\mathrm {d} x_{1}...\mathrm {d} x_{N}=1\;\;\;\Rightarrow \;\;\;Ce^{\frac {S_{0}}{k}}={\frac {\sqrt {\det {\boldsymbol {\beta }}}}{(2\pi )^{\frac {N}{2}}}}$

Les quantités thermodynamiques conjuguées sont :

$Y_{i}=-{\frac {\partial S}{\partial x_{i}}}=\beta _{ik}x_{k}$

On peut exprimer sous forme de développement les dérivées temporelles et faire apparaître les coefficients cinétiques (β est inversible):

${\dot {x}}_{i}=\lambda _{ik}x_{i}x_{k}=-L_{ik}Y_{k}\,,\;\;\;L_{ik}=\lambda _{ij}({\boldsymbol {\beta }}^{-1})_{jk}$

### Corrélations temporelles

Les relations ci-dessus permettent de calculer par intégration les corrélations temporelles entre $x_i$ et $Y_k$ :

${\overline {x_{i}Y_{k}}}=\delta _{ik}$

Les $x_i$ sont des quantités microscopiques invariantes par renversement du temps (microréversibilité). On peut donc écrire leurs moyennes temporelles de la façon suivante :

${\overline {x_{i}(t)x_{k}(t+\tau )}}={\overline {x_{i}(t+\tau )x_{k}(t)}}$

En dérivant par rapport à τ et en faisant τ=0 dans l'expression obtenue il vient :

${\overline {x_{i}{\dot {x}}_{k}}}={\overline {{\dot {x}}_{i}x_{k}}}$

En remplaçant les dérivées par leurs valeurs ci-dessus on obtient :

${\overline {x_{i}L_{kj}Y_{j}}}={\overline {L_{ij}Y_{j}{x}_{k}}}$

or :

${\overline {x_{i}L_{kj}Y_{j}}}=L_{kj}{\overline {Y_{j}{x}_{j}}}=L_{kj}\delta _{ij}=L_{ki}$

En faisant de même sur le second membre on obtient la relation de réciprocité d'Onsager :

$L_{ki}=L_{ik}$

Dans le cas où l'une des fluctuations $x_i$ change de signe par retournement du temps (microréversibilité), par exemple parce qu'elle est proportionnelle à une vitesse, la relation de corrélation devient :

${\overline {x_{i}{\dot {x}}_{k}}}=-{\overline {{\dot {x}}_{i}x_{k}}}$

D'où :

$L_{ki}=-L_{ik}$

On peut généraliser simplement cela à un certain nombre de variables en introduisant la notion de signature par rapport à l'action de retournement du temps.

Lorsqu'un champ magnétique où une rotation est présente l'invariance n'est réalisée qu'à condition d'inverser le signe des vecteurs correspondants.

