Équation de Boltzmann

L' équation de Boltzmann ou équation de transport de Boltzmann décrit le comportement statistique d'un système thermodynamique hors état d' équilibre, conçue par Ludwig Boltzmann en 1872. L'exemple classique d'un tel système est un fluide avec des gradients de température dans l'espace provoquant un flux de chaleur des régions les plus chaudes vers les plus froides, par le transport aléatoire mais orienté des particules composant ce fluide. Dans la littérature moderne, le terme équation de Boltzmann est souvent utilisé dans un sens plus général, se référant à toute équation cinétique décrivant le changement d'une quantité macroscopique dans un système thermodynamique, comme l'énergie, la charge ou le nombre de particules, que ce soit dans un cadre classique, semi-classique ou relativiste.

Simulation microscopique d'un système de N (ici petit) particules dans un container, pour lequel l'équation de Boltzmann s'applique. Sans forces externes, les collisions binairs font évoluer le système. Certaines particules sont colorées en rouge pour un meilleur suivi.

L'équation de Boltzmann est une des équations les plus importantes de la physique par l'étendue de ses applications (par exemple en mécanique des fluides) et par ses implications épistémologiques, en donnant un caractère universel au concept d'entropie introduit par Clausius en thermodynamique. Elle jette un pont entre la physique microscopique (comportement individuel des composants d'un système) et la physique macroscopique (propriétés collectives comme l'énergie thermique, la viscosité, la conductivité électrique, etc...)

D'un point de vue mathématique, c'est une équation intégro-différentielle dont l'étude des propriétés et la recherche de solutions analytiques a fortement mobilisé les mathématiciens (Seiji Ukai dans les années 1970 ; ou plus récemment les médaillés Fields Pierre-Louis Lions et Cédric Villani). Si l'existence de solutions est prouvée, même loin de l'équilibre, leur régularité et unicité reste un problème ouvert très important. Mais dans la pratique, l'équation est résolue numériquement par différentes méthodes (lagrangienne, eulérienne, particule-in-cell, etc.).

L'espace de phase et la fonction de densité

L'équation de Boltzmann détermine l'évolution d'un système de N particules (N très grand) dans l'espace des phases, espace qui contient l'ensemble des valeurs possibles des positions r et impulsions p des particules. Chaque particule est représentée par 6 nombres $(\mathbf {r} \;,\mathbf {p} )=(x,y,z\;,p_{x},p_{y},p_{z})$. L'espace des phases complet est donc a priori un espace à 6 N dimensions, l'évolution au cours du temps du système étant déterminée par une trajectoire dans l'espace des phases $F(\mathbf {r} _{1},\mathbf {p} _{1}\;;\mathbf {r} _{2},\mathbf {p} _{2}\;;...;\mathbf {r} _{N},\mathbf {p} _{N}\;;t)$, fonction (dite à N-corps) qui contient une quantité phénoménale d'information ; non calculable, elle n'a pas d'intérêt si on s'intéresse aux propriétés globales du système.

L'équation de Boltzmann est une équation à un corps qui considère $f(\mathbf {r} \,,\mathbf {p} \;;t)$ la distribution de probabilité pour la position et l'impulsion d'une particule quelconque du système ; c'est-à-dire la probabilité qu'une particule occupe à un instant donné, une très petite région donnée de l'espace des phases, de volume $\mathrm {d} ^{3}{\bf {r}}$ autour de la position ${\bf {r}}$, et de volume $\mathrm {d} ^{3}{\bf {p}}$,autour de l'impulsion ${\bf {p}}$.

${\text{d}}^{3}\mathbf {r} \,{\text{d}}^{3}\mathbf {p} ={\text{d}}x\,{\text{d}}y\,{\text{d}}z\,{\text{d}}p_{x}\,{\text{d}}p_{y}\,{\text{d}}p_{z}.$

Le cœur de l'équation de Boltzmann est donc l'évolution au cours du temps de la fonction de densité de probabilité $f(\mathbf {r} \,,\mathbf {p} \;;t)$ définie :

${\text{d}}N=f(\mathbf {r} \,,\mathbf {p} ;t)\,{\text{d}}^{3}\mathbf {r} \,{\text{d}}^{3}\mathbf {p}$

de sorte que $dN$ est le nombre de toutes les molécules dont les positions et les impulsions sont situées dans l'élément de volume centré autour de $(\mathbf {r} ,\mathbf {p} )$. Cette fonction ne dépend plus de $\mathbf {r}$ dans un système homogène dans l'espace ; dans un système homogène dans l'espace et isotrope en impulsion (la probabilité des vitesses est la même dans toutes les les directions), le seul degré de liberté est l'énergie cinétique des particules $\epsilon ={\frac {\|\mathbf {p} \|^{2}}{2m}}$, on obtient

$${\displaystyle {\text{d}}N=g(\epsilon ;t){\text{d}}\mathbf {r} \,{\text{d}}\epsilon \;{\text{avec}}\;g(\epsilon )=4\pi {\sqrt {(2m^{3}\epsilon )}}}\;f(\epsilon )\;;$$

forme que l'on retrouve pour la distribution de probabilité des systèmes en équilibre thermique (voir infra).

L'intégration sur une région de l'espace des positions et de l'espace des impulsions donne le nombre total de particules qui ont des positions et des impulsions dans cette région (intégrale sextuple):

${\begin{aligned}N&=\int \limits _{\mathrm {momenta} }{\text{d}}^{3}\mathbf {p} \int \limits _{\mathrm {positions} }{\text{d}}^{3}\mathbf {r} \,f(\mathbf {r} ,\mathbf {p} ,t)\\[5pt]&=\iiint \limits _{\mathrm {momenta} }\quad \iiint \limits _{\mathrm {positions} }f(x,y,z,p_{x},p_{y},p_{z},t)\,{\text{d}}x\,{\text{d}}y\,{\text{d}}z\,{\text{d}}p_{x}\,{\text{d}}p_{y}\,{\text{d}}p_{z}.\end{aligned}}$

La fonction à un corps f correspond au niveau le plus bas dans la hiérarchie BBGKY des équations du problème à N corps ; ce qui implique que l'on ne prend pas en compte les corrélations entre paires de particules, et donc :

$f(\mathbf {r} _{1},\mathbf {p} _{1};\mathbf {r} _{2},\mathbf {p} _{2})\approx f(\mathbf {r} _{1},\mathbf {p} _{1})f(\mathbf {r} _{2},\mathbf {p} _{2}).$

On suppose que les particules du système sont identiques (donc chacune a une masse m identique). Pour un mélange de plusieurs espèces chimiques, une distribution est nécessaire pour chacune d'entre elles, voir ci-dessous.

Principe de base

La fonction de distribution évolue au cours du temps en fonction de 3 facteurs : les forces exercées sur les particules par une influence externe (le modèle peut s'étendre aux forces s'exerçant entre les particules elles-mêmes dans l'hypothèse du champ moyen), la diffusion ("diff") des particules liée à leur mouvement dans l'espace, et leur éventuelles collisions directes ("coll"). L'équation générale peut alors s'écrire

${\frac {df}{dt}}=\left({\frac {\partial f}{\partial t}}\right)_{\text{force}}+\left({\frac {\partial f}{\partial t}}\right)_{\text{diff}}+\left({\frac {\partial f}{\partial t}}\right)_{\text{coll}}.$

Les expressions détaillées de chacune de ces contributions sont fournies ci-dessous.

Note : certains auteurs utilisent la vitesse des particules v au lieu de la quantité de mouvement p ; grandeurs liées par la définition de la quantité de mouvement par $\mathbf {p} =m\mathbf {v}$.

Termes de force, de diffusion et de collision

Dans un système de particules décrit par la distribution $f$, soumis à une force externe F, à un instant donné t, considérons le sous-ensemble de celles situées à la position r dans l'élément $d^{3}{\bf {r}}$ avec la quantité de mouvement p dans $d^{3}{\bf {p}}$. Au temps $t+\Delta t$, leur position sera $\mathbf {r} +\Delta \mathbf {r} =\mathbf {r} +{\frac {\mathbf {p} }{m}}\,\Delta t$ et leur impulsion $\mathbf {p} +\Delta \mathbf {p} =\mathbf {p} +\mathbf {F} \Delta t$. En l'absence de collisions, la conservation du nombre de particules implique que la densité dans l'espace des phases est conservée, soit :

$f\left(\mathbf {r} +{\frac {\mathbf {p} }{m}}\,\Delta t,\mathbf {p} +\mathbf {F} \,\Delta t,t+\Delta t\right)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} =f(\mathbf {r} ,\mathbf {p} ,t)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} .$

Dans cette expression on utilise le fait que l'élément de volume de l'espace des phases $d^{3}{\bf {r}}$ $d^{3}{\bf {p}}$ est constant au cours du temps (conséquence du théorème de Liouville). La seule déviation de cette identité est causée par les collisions :

${\displaystyle {\begin{aligned}dN_{\mathrm {coll} }&=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }\,\Delta t\;d^{3}\mathbf {r} \,d^{3}\mathbf {p} \\[5pt]&=f\left(\mathbf {r} +{\frac {\mathbf {p} }{m}}\Delta t,\mathbf {p} +\mathbf {F} \Delta t,t+\Delta t\right)d^{3}\mathbf {r} d^{3}\mathbf {p} -f(\mathbf {r} ,\mathbf {p} ,t)\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} \\[5pt]&=\Delta f\,d^{3}\mathbf {r} \,d^{3}\mathbf {p} \;;\end{aligned}}}$

où $\Delta f$ est la variation totale de$f$. En passant aux limites $\Delta t\rightarrow 0$ et $\Delta f\rightarrow 0$, on obtient :

${\frac {df}{dt}}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }$

(2)

Par ailleurs, la conservation de la densité dans l'espace implique que, sans collisions, la différentielle totale suivante soit nulle :

$$df={\frac {\partial f}{\partial t}}dt+\nabla _{r}f.d\mathbf {r} +\nabla _{p}.d\mathbf {p} =0$$

${\displaystyle {\begin{aligned}df&={\frac {\partial f}{\partial t}}\,dt+\left({\frac {\partial f}{\partial x}}\,dx+{\frac {\partial f}{\partial y}}\,dy+{\frac {\partial f}{\partial z}}\,dz\right)+\left({\frac {\partial f}{\partial p_{x}}}\,dp_{x}+{\frac {\partial f}{\partial p_{y}}}\,dp_{y}+{\frac {\partial f}{\partial p_{z}}}\,dp_{z}\right)\\[5pt]&={\frac {\partial f}{\partial t}}dt+\nabla _{\mathbf {r} }f\cdot d\mathbf {r} +{\frac {\partial f}{\partial \mathbf {p} }}\cdot d\mathbf {p} \\[5pt]&={\frac {\partial f}{\partial t}}dt+\nabla _{\mathbf {r} }f\cdot {\frac {\mathbf {p} }{m}}dt+{\frac {\partial f}{\partial \mathbf {p} }}\cdot \mathbf {F} \,dt\end{aligned}}}$

avec le produit scalaire ( . ), et l'opérateur gradient $\nabla _{r}$ usuel pour les positions, qui se formule de manière analogue pour les moments :

$\nabla _{\mathbf {p} }f=\mathbf {\hat {e}} _{x}{\frac {\partial f}{\partial p_{x}}}+\mathbf {\hat {e}} _{y}{\frac {\partial f}{\partial p_{y}}}+\mathbf {\hat {e}} _{z}{\frac {\partial f}{\partial p_{z}}},$

avec ${\widehat {e}}_{x},{\widehat {e}}_{y},{\widehat {e}}_{z}$ les vecteurs unitaires en coordonnées cartésiennes.

Forme finale

Regroupant tous les termes, on obtient l'équation de Boltzmann :

${\frac {\partial f}{\partial t}}+{\frac {\mathbf {p} }{m}}\cdot \nabla f+\mathbf {F} \cdot {\frac {\partial f}{\partial \mathbf {p} }}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }.$

Dans ce contexte, $\mathbf {F} (\mathbf {r} ,t)$ est le champ de force agissant sur les particules du système, et m est la masse des particules. Lorsque les forces dérivent d'un potentiel $V({\mathbf {r} })$, on peut exprimer l'équation de Boltzmann en fonction de l'hamiltonien $H(\mathbf {r} ,\mathbf {p} )={\frac {\mathbf {p} ^{2}}{2m}}+V(\mathbf {r} )$ avec le crochet de Poisson :

${\frac {\partial f}{\partial t}}+\{f,H\}=\left({\frac {\partial f}{\partial t}}\right)_{\mathrm {coll} }$

Lorsque les forces entre particules sont répulsives à longue portée, comme l'interaction coulombienne, les collisions directes sont impossibles ; on obtient l'équation de Boltzmann sans collision, généralement appelée équation de Vlassov. L'équation de Boltzmann quantique est obtenue lorsque la densité dans l'espace des phases et les collisions entre particules tiennent compte des effets du principe d'exclusion de Pauli pour les fermions.

Sans collisions, l'équation se réduit au formalisme classique de la mécanique hamiltonienne ; le calcul du terme de collisions est donc essentiel et repose sur le type de statistique qui s'applique aux particules du système (Maxwell-Boltzmann, Fermi-Dirac ou Bose-Einstein).

Solutions à l'équilibre

Lorsque le système est homogène et n'est pas soumis à des forces externes, Boltzmann a démontré que la distribution f tend asymptotiquement vers une limite, qui est la distribution de probabilité de Maxwell-Boltzmann. Comme vu plus haut, le seul degré de liberté des particules étant leur énergie, la loi donne le nombre de particules dans un petit domaine d'énergie $dn_{E}$

$$\mathrm {d} n_{E}={\frac {N}{\int g(\varepsilon )\mathrm {e} ^{-\varepsilon /k_{B}T}\mathrm {d} \varepsilon }}~g(E)\mathrm {e} ^{\frac {-E}{k_{B}T}}\,\mathrm {d} E\;;$$

où N est le nombre totale de particule, T est la température d'équilibre, $k_{B}$ la constante de Boltzmann et $g(E)$ la densité de niveau d'énergie. Cette distribution correspond au maximum de l'entropie du système (voir Théorème H).

Le terme de collision (Stosszahlansatz) et le chaos moléculaire

Terme de collision à deux corps

Une idée clé de Boltzmann a été de calculer le terme de collision en le limitant aux collisions à deux corps entre des particules, supposées ponctuelles et non corrélées avant la collision. Cette hypothèse désignée comme le " Stosszahlansatz " et est également connue sous le nom de " hypothèse du chaos moléculaire ". Sous cette hypothèse, et parce que l'interaction est locale dans l'espace des positions, le terme de collision peut être écrit comme une intégrale dans l'espace des impulsions sur le produit des fonctions de distribution à une particule ;

$\left({\frac {\partial f}{\partial t}}\right)_{\text{coll}}=\iint g\,I(g,\Omega )[f(\mathbf {r} ,\mathbf {p} '_{A},t)f(\mathbf {r} ,\mathbf {p} '_{B},t)-f(\mathbf {r} ,\mathbf {p} _{A},t)f(\mathbf {r} ,\mathbf {p} _{B},t)]\,d\Omega \,d^{3}\mathbf {p} _{B},$

où $\mathbf {p} _{A}$et $\mathbf {p} _{B}$ sont les impulsions de deux particules quelconques avant la collision (étiquetées A et B pour plus de commodité) ; $\mathbf {p} '_{A}$et $\mathbf {p} '_{B}$ sont les impulsions après la collision ; la conservation des impulsions (chocs élastiques) implique que $g$, la valeur de la quantité de mouvement relative, soit conservée :

$g=|\mathbf {p} _{B}-\mathbf {p} _{A}|=|\mathbf {p} '_{B}-\mathbf {p} '_{A}|.$

$I(g,\Omega )$ est la section efficace différentielle de la collision ; la collision se traduit donc par une rotation d'un angle θ dans l'angle solide $d\Omega$ de la quantité de mouvement relative des particules.

Simplifications du terme de collision

Le principal défi dans la résolution de l'équation de Boltzmann provient du terme de collision. Désormais, la puissance des ordinateurs permet de faire les calculs explicites pour des milliers de particules (voir beaucoup plus, voir simulation d'un système à N-corps) en interaction. Historiquement, des tentatives ont été faites pour «modéliser» et simplifier le terme de collision. L'équation modèle la plus connue est l'équation BGK due à Bhatnagar, Gross et Krook  ; elle est obtenue en faisant l'hypothèse que l'effet des collisions moléculaires est de forcer une fonction de distribution hors d'équilibre à revenir à la fonction d'équilibre (fonction maxwellienne dans le cas classique) en tous points de l'espace des positions, et que la vitesse à laquelle cela se produit est proportionnelle à la fréquence des collisions . L'équation de Boltzmann est donc modifiée sous la forme BGK :

${\frac {\partial f}{\partial t}}+{\frac {\mathbf {p} }{m}}\cdot \nabla f+\mathbf {F} \cdot {\frac {\partial f}{\partial \mathbf {p} }}=\nu (f_{0}-f),$

où $\nu$ est la fréquence de collision, et $f_{0}$ est la fonction de distribution locale (maxwellienne dans le cas classique, voir ci-dessus) compte tenu de la température du gaz en ce point de l'espace.

Équation générale (pour un mélange)

Pour un mélange d'espèces chimiques étiquetées par les indices $i=1,2,3,...,n$, l'équation pour une espèce i est

${\frac {\partial f_{i}}{\partial t}}+{\frac {\mathbf {p} _{i}}{m_{i}}}\cdot \nabla f_{i}+\mathbf {F} \cdot {\frac {\partial f_{i}}{\partial \mathbf {p} _{i}}}=\left({\frac {\partial f_{i}}{\partial t}}\right)_{\text{coll}},$

où le terme de collision s'écrit :

$\left({\frac {\partial f_{i}}{\partial t}}\right)_{\mathrm {coll} }=\sum _{j=1}^{n}\iint g_{ij}I_{ij}(g_{ij},\Omega )[f'_{i}f'_{j}-f_{i}f_{j}]\,d\Omega \,d^{3}\mathbf {p'} ,$

avec les notations simplifiées $f_{i}=f(\mathbf {r} _{i},\mathbf {p} _{i},t)$ et $f'_{i}=f(\mathbf {r} '_{i},\mathbf {p} '_{i},t)$ ; l'amplitude des moments relatifs est

$g_{ij}=|\mathbf {p} _{i}-\mathbf {p} _{j}|=|\mathbf {p'} _{i}-\mathbf {p'} _{j}|,$

et $I_{ij}$ est la section efficace différentielle entre les particules des espèces i et j. La somme des intégrales décrit l'entrée et la sortie des particules d'espèces i dans ou hors de l'élément de l'espace des phases, sous l'effet des collisions avec les particules de toutes les espèces.

Applications et extensions

Équations de conservation

L'équation de Boltzmann peut être utilisée pour dériver les lois de conservation de la dynamique des fluides pour la masse, la charge, la quantité de mouvement et l'énergie. Pour un fluide constitué d'un seul type de particules, la densité numérique $n$ (nombre de particules par unité de volume d'espace) est donnée par

$n=\int f\,d^{3}\mathbf {p} .$

La valeur moyenne locale de toute fonction $A(\mathbf {r} ,\mathbf {p} )$ est :

$\langle A\rangle ={\frac {1}{n}}\int Af\,d^{3}p.$$\langle A\rangle ={\frac {1}{n}}\int Af\,d^{3}\mathbf {p} .$

Étant donné que les équations de conservation impliquent des tenseurs, la convention de sommation d'Einstein sera utilisée : les indices répétés dans un produit indiquent une sommation sur ces indices. Ainsi $\mathbf {x} \mapsto x_{i}$ et $\mathbf {p} \mapsto p_{i}=mv_{i}$, où $v_{i}$ est le vecteur vitesse des particules. Considérons une quantité $A(p_{i})$, une fonction de la quantité de mouvement $p_{i}$ seulement, qui est conservée dans une collision.

Supposons également que la force $F_{i}$ est une fonction de la position uniquement, et que la distribution $f$ est nulle pour $p_{i}\to \pm \infty$ . En multipliant l'équation de Boltzmann par A et en intégrant sur la quantité de mouvement, on obtient quatre termes qui, en utilisant l'intégration par parties, peuvent être exprimés comme :

$\int A{\frac {\partial f}{\partial t}}\,d^{3}p={\frac {\partial }{\partial t}}(n\langle A\rangle ),$

$\int {\frac {p_{j}A}{m}}{\frac {\partial f}{\partial x_{j}}}\,d^{3}p={\frac {1}{m}}{\frac {\partial }{\partial x_{j}}}(n\langle Ap_{j}\rangle ),$

$\int AF_{j}{\frac {\partial f}{\partial p_{j}}}\,d^{3}p=-nF_{j}\left\langle {\frac {\partial A}{\partial p_{j}}}\right\rangle ,$

$\int A\left({\frac {\partial f}{\partial t}}\right)_{\text{coll}}\,d^{3}p=0,$

où le dernier terme est nul, puisque A est conservé dans une collision élastique. Les valeurs de A correspondent aux moments de vitesse $v_{i}$ (et l'impulsion $p_{i}$, car ils sont linéairement dépendants).

Moment zéro

Soit $A=m(v_{i})^{0}=m$, la masse de la particule, l'équation de Boltzmann intégrée devient l'équation de conservation de la masse  :

${\frac {\partial }{\partial t}}\rho +{\frac {\partial }{\partial x_{j}}}(\rho V_{j})=0,$

où $\rho =mn$ est la masse volumique, et $V_{i}=\langle v_{i}\rangle$ est la vitesse moyenne du fluide.

Premier moment

Soit $A=m(v_{i})^{1}=p_{i}$, la quantité de mouvement des particules, l'équation de Boltzmann intégrée devient l'équation de conservation de la quantité de mouvement  :

${\frac {\partial }{\partial t}}(\rho V_{i})+{\frac {\partial }{\partial x_{j}}}(\rho V_{i}V_{j}+P_{ij})-nF_{i}=0,$

où $P_{ij}=\rho \langle (w_{i}-V_{i})(w_{j}-V_{j})\rangle$ est le tenseur de pression (soit le tenseur des contraintes visqueuses plus la pression hydrostatique , voir l'article Équations_de_Navier-Stokes pour leurs définitions).

Deuxième moment

Posant $A={\frac {m(v_{i})^{2}}{2}}={\frac {p_{i}p_{i}}{2m}}$, l'énergie cinétique des particules, l'équation de Boltzmann intégrée devient l'équation de conservation de l'énergie  :

${\frac {\partial }{\partial t}}(u+{\tfrac {1}{2}}\rho V_{i}V_{i})+{\frac {\partial }{\partial x_{j}}}(uV_{j}+{\tfrac {1}{2}}\rho V_{i}V_{i}V_{j}+J_{qj}+P_{ij}V_{i})-nF_{i}V_{i}=0,$

où $u={\tfrac {1}{2}}\rho \langle (w_{i}-V_{i})(w_{i}-V_{i})\rangle$ est la densité d'énergie thermique cinétique, et $J_{qi}={\tfrac {1}{2}}\rho \langle (w_{i}-V_{i})(w_{k}-V_{k})(w_{k}-V_{k})\rangle$ est le vecteur de flux de chaleur.

Mécanique hamiltonienne

En mécanique hamiltonienne, l'équation de Boltzmann s'écrit souvent plus généralement comme

${\hat {\mathbf {L} }}[f]\ =\ \mathbf {C} [f]$

où $\mathbf {C}$ est l'opérateur de collision, opérateur intégral non linéaire.

Théorème de Lanford

Ce théorème étudie la validité de l'équation de Boltzmann lorsque les particules ne sont pas considérées comme ponctuelles.

Limite de Boltzmann-Grad

La limite dite de Boltzmann-Grad consiste à prendre la limite conjointe :

• d'un nombre de particules $N\ \to +\ \infty$ ;
• d'un rayon $r\ \to \ 0$ ;

en maintenant le produit $N\ r^{2}={\text{cte}}$. En particulier, le volume exclu tend vers zéro dans cette limite : $N\ r^{3}\ \to \ 0$.

Théorème de Lanford (1973)

Lanford a démontré rigoureusement qu'un gaz de sphères dures dilué dans $\mathbb {R} ^{3}$ obéit à l'équation de Boltzmann dans la limite de Boltzmann-Grad, au moins pour un temps très court, égal seulement à un cinquième du temps de parcours moyen d'une particule.

En dépit de cette restriction sur la durée, ce théorème mathématique rigoureux est très important conceptuellement, puisque l'équation de Boltzmann entraîne le théorème H, à propos duquel Boltzmann fut accusé de pratiquer des « mathématiques douteuses ». Il n'en demeure pas moins qu'il reste à démontrer que ce résultat reste vrai pour des temps macroscopiques, ainsi que lorsque les atomes sont confinés dans une boîte.

Théorie quantique et violation de la conservation du nombre de particules

Il est possible d'écrire des équations de Boltzmann quantiques (ou plus exactement semi-classiques) qui prennent en compte les contraintes typiquement quantiques sur la densité dans l'espace des phases. L'équation la plus connue est obtenue en appliquant l'approximation du champ moyen pour les interactions moyenne entre particules et un terme de collisions modifié pour interdire les collisions conduisant à des états finaux interdits par le principe de Pauli (le terme de collision de Uehling-Uhlenbeck)

De plus, dans aux hautes énergies, dans les systèmes quantiques relativistes le nombre de particules n'est pas conservé dans les collisions. Cela a plusieurs applications en cosmologie physique, y compris la formation des éléments légers dans la nucléosynthèse du Big Bang, la production de matière noire et la baryogénèse. Il n'est pas établi a priori que l'état d'un système quantique puisse être caractérisé par une densité d'espace de phase classique f. Cependant, pour une large classe d'applications, il existe une généralisation bien définie de f qui est la solution d'une équation de Boltzmann efficace qui peut être dérivée des premiers principes de la théorie quantique des champs.

Relativité générale et astronomie

L'équation de Boltzmann est utile en dynamique galactique. Une galaxie, sous certaines hypothèses, peut être approchée comme un fluide continu ; sa distribution de masse est alors représentée par f ; dans les galaxies, les collisions physiques entre les étoiles sont très rares, et l'effet des collisions gravitationnelles peut être négligé pendant des temps bien plus longs que l' âge de l'univers. L'opérateur relativiste ${\hat {L}}$ s'écrit

${\hat {\mathbf {L} }}_{\mathrm {GR} }=p^{\alpha }{\frac {\partial }{\partial x^{\alpha }}}-\Gamma ^{\alpha }{}_{\beta \gamma }p^{\beta }p^{\gamma }{\frac {\partial }{\partial p^{\alpha }}},$

où $\Gamma _{\;\beta \gamma }^{\alpha }$ est le symbole de Christoffel de seconde espèce (cela suppose qu'il n'y a pas de forces externes, de sorte que les particules se déplacent le long des géodésiques en l'absence de collisions), avec la subtilité importante que la densité est une fonction dans l'espace des phases en mode mixte contravariant-covariant $(x^{i},p_{i})$ par opposition à l'espace des phases entièrement contravariant $(x^{i},p^{i})$ .

En cosmologie physique, l'approche entièrement covariante a été utilisée pour étudier le rayonnement de fond diffus cosmologique.Plus généralement, l'étude des processus dans l' univers primitif tente souvent de prendre en compte les effets de la mécanique quantique et de la relativité générale.Dans le milieu très dense formé par le plasma primordial après le Big Bang, des particules se créent et s'annihilent en permanence. Dans un tel environnement, la cohérence quantique et l'extension spatiale de la fonction d' onde peuvent affecter la dynamique, ce qui rend douteux que la distribution d'espace de phase classique f qui apparaît dans l'équation de Boltzmann soit appropriée pour décrire le système. Dans de nombreux cas, il est cependant possible de dériver une équation de Boltzmann efficace pour une fonction de distribution généralisée à partir des premiers principes de la théorie quantique des champs. Cela inclut la formation des éléments légers dans la nucléosynthèse primordiale, la production de matière noire et la baryogénèse.

Résoudre l'équation

Des solutions exactes aux équations de Boltzmann se sont avérées exister dans certains cas ; cette approche analytique donne un aperçu, mais n'est généralement pas utilisable dans des problèmes pratiques.

Au lieu de cela, les méthodes numériques (y compris les éléments finis et les méthodes de Boltzmann sur réseau) sont généralement utilisées pour trouver des solutions approximatives aux différentes formes de l'équation de Boltzmann. Les exemples d'applications vont de l' aérodynamique hypersonique dans les flux de gaz raréfiés , aux flux de plasma. Une application de l'équation de Boltzmann en électrodynamique est le calcul de la conductivité électrique - le résultat est au premier ordre identique au résultat semi-classique.

Lorsque le système est proche de l'équilibre local, la solution de l'équation de Boltzmann peut être représentée par un développement asymptotique en puissances du nombre de Knudsen (le développement de Chapman-Enskog ). Les deux premiers termes de ce développement donnent les équations d'Euler et les équations de Navier-Stokes. Les termes supérieurs ont des singularités. Le problème du développement mathématique des processus limites, qui conduisent de la vision atomistique (représentée par l'équation de Boltzmann) aux lois du mouvement des continus, est une partie importante du sixième problème de Hilbert.

Limitations et autres utilisations de l'équation de Boltzmann

L'équation de Boltzmann n'est valable que sous plusieurs hypothèses. Par exemple, les particules sont supposées être ponctuelles, c'est-à-dire sans avoir une taille finie. Il existe une généralisation de l'équation de Boltzmann qui s'appelle l'équation de Enskog (voir plus haut).

Les particules sont supposées n'avoir aucun de degré de liberté interne, aussi les collisions sont considérées comme élastiques. S'il existe un degré de liberté interne, l'équation de Boltzmann doit être généralisée pour prendre en compte des collisions inélastiques.

De nombreux fluides réels comme les liquides ou les gaz denses ont en plus des caractéristiques mentionnées ci-dessus des formes de collisions plus complexes, il y aura non seulement des collisions binaires, mais aussi ternaires et d'ordre supérieur. Celles-ci doivent être dérivées en utilisant la hiérarchie BBGKY.

Des équations de type Boltzmann sont également utilisées pour le mouvement des cellules. Puisque les cellules sont des particules composites qui portent un degré de liberté interne, les équations de Boltzmann généralisées correspondantes doivent avoir des intégrales de collision inélastiques. De telles équations peuvent décrire les invasions de cellules cancéreuses dans les tissus, la morphogenèse et les effets liés à la chimiotaxie .

Voir également

 

Théorème H

Théorie cinétique des gaz

Physique statistique hors d'équilibre Théorème H Hiérarchie BBGKY Équation de Vlassov Équation de Poisson-Boltzmann

Équations de Navier-Stokes

Références

Bibliographie

• Anouk Barberousse, La mécanique statistique - De Clausius à Gibbs, Collection Histoire des Sciences, Belin (2002) (ISBN 2-7011-3073-5). Cette collection originale propose une histoire du développement de la théorie cinétique des gaz basée sur des extraits des grands textes fondateurs (traduits en français) mis en perspective contemporaine par une historienne des sciences (CNRS). Accessible dès le niveau premier cycle universitaire.
• Ludwig Boltzmann, Leçons sur la théorie des gaz, Gauthier-Villars (1902-1905). Réédition Jacques Gabay (1987) (ISBN 2-87647-004-7).
• Carlo Cercignani, Ludwig Boltzmann - The man who Trusted Atoms, Oxford University Press (1998) 330 p. (ISBN 0-19-850154-4). Biographie scientifique du professeur Boltzmann, qui a porté la théorie cinétique des gaz à son acmée. Par un professeur de Physique mathématique de l'université de Milan (Italie), spécialiste de l'équation de Boltzmann. Niveau plutôt second cycle universitaire.
• Carlo Cercignani, The Boltzmann Equation & its Applications, Series: Applied Mathematical Sciences, 67, Springer-Verlag (1987) (ISBN 0-387-96637-4).
• Carlo Cercignani, Reinhard Illner et Mario Pulvirenti, The Mathematical Theory of Dilute Gases, Series: Applied Mathematical Sciences, 106, Springer-Verlag (1994) (ISBN 0-387-94294-7).
• Carlo Cercignani et Gilberto M. Kremer, The Relativistic Boltzmann Equation: Theory and Applications, Series: Progress in Mathematical Physics, 22, Birkhäuser (2002) (ISBN 3-7643-6693-1).
• Lawrence Sklar, Physics and Chance - Philosophical Issues in the Foundations of Statistical Mechanics, Cambridge University Press (1995) (ISBN 0-521-55881-6).

• Stewart Harris, An introduction to the theory of the Boltzmann equation, Dover Books, 1971, 221 p. (ISBN 978-0-486-43831-3, lire en ligne). Introduction très peu coûteuse au cadre moderne (à partir d'une déduction formelle de Liouville et de la hiérarchie Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) dans laquelle se place l'équation de Boltzmann). La plupart des manuels de mécanique statistique comme Huang traitent encore le sujet en utilisant les arguments originaux de Boltzmann. Pour dériver l'équation, ces livres utilisent une explication heuristique qui ne fait pas ressortir la plage de validité et les hypothèses caractéristiques qui distinguent les équations de Boltzmann des autres équations de transport comme les équations de Fokker-Planck ou de Landau.

• Arkeryd, « On the Boltzmann equation part I: Existence », Arch. Rational Mech. Anal., vol. 45, n^o 1,‎ 1972, p. 1–16 (DOI 10.1007/BF00253392, Bibcode 1972ArRMA..45....1A, S2CID 117877311)
• Arkeryd, « On the Boltzmann equation part II: The full initial value problem », Arch. Rational Mech. Anal., vol. 45, n^o 1,‎ 1972, p. 17–34 (DOI 10.1007/BF00253393, Bibcode 1972ArRMA..45...17A, S2CID 119481100)
• DiPerna et Lions, « On the Cauchy problem for Boltzmann equations: global existence and weak stability », Ann. of Math., 2^e série, vol. 130, n^o 2,‎ 1989, p. 321–366 (DOI 10.2307/1971423, JSTOR 1971423)

Liens externes

• L'équation de transport de Boltzmann par Franz Vesely
• « Équation de Boltzmann : un pas en avant ! », La Méthode scientifique, France Culture, 5 mai 2022.
• Comportements gazeux de Boltzmann résolus