Multispecies Bhatnagar-Gross-Krook models and the Onsager reciprocal relations

L'article démontre que la plupart des modèles existants de type Bhatnagar-Gross-Krook indépendants de la vitesse pour les systèmes multispecies ne satisfont pas les relations de réciprocité d'Onsager, ce qui complique l'étalonnage de leurs propriétés de transport afin de les faire correspondre à des fluides spécifiques.

L'essentiel

Imaginez que vous essayez de prédire comment deux types différents de billes (disons rouges et bleues) se mélangent lorsque vous secouez une boîte. Les scientifiques disposent d'un « code de règles » décrivant le comportement de ces billes, appelé l'équation de Boltzmann. Elle est incroyablement précise, mais si complexe que la résoudre revient à compter chaque grain de sable d'une plage pendant qu'une ouragan souffle.

Pour simplifier les choses, les scientifiques ont créé une version simplifiée de ce code de règles, appelée le modèle BGK (du nom de Bhatnagar, Gross et Krook). Considérez le modèle BGK comme une « fiche d'aide » ou un « raccourci » qui approxime le comportement complexe des billes sans effectuer tous les calculs lourds. Il est utilisé depuis des décennies pour simuler tout, depuis l'écoulement de l'air sur une aile jusqu'au plasma d'une étoile.

Le Problème : Une Boussole Cassée

Cet article, écrit par E. S. Benilov, met en évidence une faille majeure dans les versions les plus courantes de cette « fiche d'aide » lorsqu'il s'agit de mélanges de gaz différents (comme l'oxygène et l'azote, ou la vapeur d'eau et l'air).

L'auteur découvre que ces modèles BGK populaires violent une loi fondamentale de la physique connue sous le nom de relations de réciprocité d'Onsager.

Voici une analogie simple pour les relations d'Onsager :
Imaginez une rue à double sens où la circulation s'écoule entre deux villes.

• Règle A : Si vous construisez une colline dans la Ville A, cela affecte la vitesse des voitures dans la Ville B.
• Règle B : Si vous construisez une colline dans la Ville B, cela affecte la vitesse des voitures dans la Ville A.

Les relations d'Onsager stipulent que ces deux effets doivent être parfaitement « accordés » l'un à l'autre. Si la colline dans la Ville A ralentit le trafic dans la Ville B de 10 %, alors la colline dans la Ville B doit ralentir le trafic dans la Ville A d'une quantité mathématiquement liée. C'est une règle de symétrie ; l'univers exige que ces interactions s'équilibrent.

Ce Que L'Article A Découvert

Benilov a testé la « fiche d'aide » BGK standard contre cette règle. Il a constaté que :

1. Le Modèle est Asymétrique : Dans le modèle BGK, la « colline » dans la Ville A (un changement de température) a aucun effet sur le trafic dans la Ville B (écoulement de masse). Cependant, la « colline » dans la Ville B (un changement de densité) affecte le trafic dans la Ville A (flux de chaleur).
2. Le Désaccord : Parce qu'un côté de l'équation est nul et l'autre non, la symétrie est brisée. Le modèle est comme une balance qui est définitivement penchée d'un côté.
3. La Conséquence : Parce que le modèle enfreint cette règle fondamentale, il est impossible de le « calibrer ». Le calibrage consiste à accorder une radio pour obtenir un signal clair. Si vous essayez d'ajuster les boutons (paramètres) du modèle BGK pour qu'il corresponde aux données réelles d'un fluide spécifique, vous ne pouvez pas. Le modèle est fondamentalement défectueux d'une manière qui l'empêche d'être jamais parfaitement précis, peu importe comment vous le modifiez.

L'Exception « Vapeur d'Eau » (et pourquoi elle ne sauve pas la mise)

Vous pourriez penser : « Eh bien, peut-être que cela ne compte que pour des gaz étranges. Qu'en est-il des mélanges courants comme la vapeur d'eau et l'air ? »

L'article vérifie cela aussi. Même si l'effet de la température sur l'écoulement de masse est minuscule (ce qui est le cas pour la vapeur d'eau et l'air), le modèle échoue toujours. Pour faire fonctionner le modèle dans ce cas spécifique, vous devriez tourner un cadran à l'infini, ce qui brise le modèle entièrement en arrêtant tout mouvement. Ainsi, le modèle échoue aussi bien pour les mélanges complexes que pour les simples.

Existe-t-il de Bons Modèles ?

L'article note qu'il existe quelques autres modèles BGK, plus complexes, qui respectent les règles, mais ils ont leurs propres problèmes (comme enfreindre d'autres lois de la physique, telles que le « théorème H », qui garantit que l'entropie augmente toujours).

L'auteur conclut que, actuellement, aucun modèle BGK existant n'est parfait. Un modèle parfait devrait :

• Conserver la masse, la quantité de mouvement et l'énergie.
• Respecter les lois de la thermodynamique (entropie).
• Traiter les particules identiques équitablement.
• Maintenir les températures et les concentrations positives.
• Permettre aux scientifiques de l'ajuster pour correspondre à n'importe quel fluide réel.
• Et obéir aux relations de réciprocité d'Onsager (la règle de symétrie).

Pour l'instant, chaque modèle que nous possédons échoue à au moins l'un de ces tests.

La Conclusion

L'article est un avertissement aux scientifiques qui utilisent ces modèles. Si vous utilisez le modèle BGK standard pour simuler des mélanges de gaz, vous utilisez un outil fondamentalement « désaccordé » par rapport aux lois de la physique. Il pourrait vous donner une idée approximative de ce qui se passe, mais vous ne pouvez pas lui faire confiance pour fournir des résultats précis et calibrés pour des fluides réels, car il viole une règle de symétrie fondamentale de la nature. L'auteur espère qu'à l'avenir, quelqu'un construira un modèle « parfait » qui résoudra tous ces problèmes.

Résumé technique

Résumé technique : Modèles multispecies de Bhatnagar–Gross–Krook et relations de réciprocité d'Onsager

Énoncé du problème
Le modèle de Bhatnagar–Gross–Krook (BGK) est une approximation phénoménologique largement utilisée de l'équation cinétique de Boltzmann, en particulier pour les mélanges gazeux. Bien que le modèle original mono-espèce et ses extensions aux mélanges (par exemple, Gross et Krook [2]) soient appréciés pour leur justesse qualitative et leur adaptabilité, ils n'ont pas été rigoureusement testés quant à leur conformité avec les relations de réciprocité d'Onsager. Ces relations imposent des contraintes thermodynamiques fondamentales sur les coefficients de transport, reliant spécifiquement le coefficient du gradient de température dans le flux de masse (thermodiffusion ou effet Soret) au coefficient du gradient de densité dans le flux de chaleur (effet Dufour). Les modèles dérivés de premiers principes satisfont automatiquement ces relations, tandis que les modèles phénoménologiques peuvent ne pas le faire. Le non-respect de ces relations par un modèle jette un doute sur sa pertinence physique et, de manière cruciale, empêche l'étalonnage des paramètres du modèle pour correspondre aux propriétés de transport de fluides réels spécifiques.

Méthodologie
L'article formule le modèle BGK multispecies le plus général pour un mélange binaire de gaz monoatomiques. Cette formulation inclut des fréquences de collision ($\nu_{ij}$) et des paramètres ($\alpha_i, \beta_i, \gamma_i$) qui définissent la relaxation vers des distributions de Maxwell mixtes ($M_{ij}$) avec des vitesses et des températures dépendantes de l'espèce. L'auteur s'assure que le modèle conserve la masse, la quantité de mouvement et l'énergie en imposant des contraintes spécifiques sur ces paramètres.

Pour analyser les propriétés de transport, l'article utilise une approximation de diffusion plutôt que l'approximation hydrodynamique standard. Cela implique :

1. Mise à l'échelle : Introduction d'un petit paramètre $\epsilon$ pour étirer les coordonnées spatiales ($\nabla \to \epsilon \nabla$) et le temps ($\partial_t \to \epsilon^2 \partial_t$), approprié pour des écoulements de diffusion lents et faibles.
2. Développement : Développement des fonctions de distribution, des vitesses macroscopiques, des températures et des densités en puissances de $\epsilon$.
3. Dérivation : Substitution de ces développements dans les équations BGK redimensionnées pour résoudre la correction du premier ordre de la fonction de distribution ($f^{(1)}_i$).
4. Calcul des flux : Intégration de $f^{(1)}_i$ pour dériver des expressions explicites des flux de masse ($J_i$) et du flux de chaleur ($Q$) en fonction des gradients de densité et de température.
5. Comparaison : Comparaison des flux BGK dérivés avec les expressions hydrodynamiques standard (théorie d'Enskog–Chapman) qui satisfont les relations d'Onsager.

Contributions et résultats clés
Le résultat central de l'article est la démonstration que la version la plus courante du modèle BGK multispecies — où les fréquences de collision et les coefficients du modèle sont indépendants de la vitesse moléculaire — ne satisfait pas les relations de réciprocité d'Onsager.

Plus précisément, l'analyse révèle une incohérence fondamentale :

• Selon les relations d'Onsager, le coefficient reliant le gradient de température au flux de masse (thermodiffusivité, $B$) doit être non nul si le coefficient reliant le gradient de densité au flux de chaleur (effet Dufour, $C$) est non nul, et ils sont liés par $C = pB$.
• La dérivation du modèle BGK donne un flux de masse où le terme thermodiffusif est strictement nul ($B=0$) sauf si des limites spécifiques et non physiques sont prises.
• Cependant, le flux de chaleur résultant contient un terme proportionnel au gradient de densité (l'effet Dufour) qui est non nul et proportionnel à la différence de masse des espèces.
• Par conséquent, le modèle BGK prédit un effet Dufour non nul parallèlement à un effet Soret nul, violant la relation d'Onsager $C = pB$.

L'article montre en outre que cette divergence ne peut être résolue par l'étalonnage des paramètres du modèle (tel que le paramètre $\beta$ gouvernant la relaxation de la vitesse relative). Aucune valeur de $\beta$ n'existe qui satisfasse simultanément les équations de transport standard et les relations d'Onsager pour des mélanges généraux. Même dans le cas spécifique de fluides à thermodiffusivité négligeable (où $B \approx 0$), le modèle BGK échoue à reproduire le comportement correct du flux de chaleur sans forcer l'ensemble du flux diffusif à s'annuler.

Signification et affirmations
L'article affirme que ce non-respect pose un problème significatif pour l'application pratique des modèles BGK multispecies. Parce que les coefficients de transport dérivés du modèle ne respectent pas la réciprocité thermodynamique nécessaire, il est impossible d'« étalonner » le modèle — c'est-à-dire de choisir des valeurs de paramètres qui fassent correspondre les propriétés de transport du modèle à celles d'un fluide réel spécifique (tel que la vapeur d'eau et l'air).

L'auteur note que, bien que le modèle échoue à satisfaire exactement les relations d'Onsager, il les satisfait de manière asymptotique dans deux limites spécifiques :

1. Lorsque le rapport de masses $m_1/m_2 \to 0$ (pertinent pour les plasmas ionisés où la masse des électrons est négligeable par rapport à celle des ions).
2. Lorsque le rapport de masses $m_1/m_2 \to 1$ (où le premier terme de l'expression du flux de chaleur s'annule).

L'article conclut en énumérant les exigences d'un modèle cinétique multispecies « parfait », qui incluent les lois de conservation, le théorème H, le principe d'indifférenciation, la positivité de la température/concentration, la capacité à représenter des nombres de Prandtl/Schmidt arbitraires, et la conformité avec les relations d'Onsager. L'auteur déclare qu'à ce jour, aucun modèle de type BGK existant n'a été démontré comme satisfaisant simultanément toutes ces exigences. Les résultats suggèrent que, bien que les modèles BGK actuels soient des approximations utiles, leur incapacité à satisfaire les relations d'Onsager limite leur fidélité physique et leur utilité pour un étalonnage précis dans les applications générales de la dynamique des fluides.