Lattice Boltzmann modeling for conduction heat transfer based on the Qadyan method

Résumé technique : Modélisation de Boltzmann sur réseau pour la conduction thermique basée sur la méthode de Qadyan

Énoncé du problème
L'article traite de la simulation du transfert de chaleur par conduction en utilisant la méthode de Boltzmann sur réseau (LBM). Bien que la LBM soit efficace pour la dynamique des fluides, son application à la conduction thermique se heurte à des limitations concernant l'uniformité de la grille et la stabilité numérique. Les schémas LBM standards, particulièrement ceux basés sur l'approximation de Bhatnagar–Gross–Krook (BGK), sont souvent limités aux grilles uniformes avec des pas spatiaux égaux. De plus, les méthodes de discrétisation spatiale standard, telles que les schémas de différences finies de second ordre, peuvent induire des oscillations numériques dans la solution. L'étude cherche à évaluer des stratégies de discrétisation spatiale alternatives pour améliorer la précision et la stabilité sans introduire de complexité supplémentaire au cadre standard de la LBM ou recourir à des grilles non uniformes.

Méthodologie
La recherche implémente et compare trois approches distinctes pour résoudre l'équation de Boltzmann sur réseau (LBE) pour la conduction thermique :

1. LBM standard : Utilise des modèles de réseau conventionnels (D2Q9, D2Q8 et D2Q5) avec des processus de collision-streaming standards. L'étude note que pour les problèmes purement diffusifs, les modèles avec moins d'ensembles de vitesses (comme D2Q8 et D2Q5) sont suffisants, car la diffusion thermique ne présente pas de préférence directionnelle.
2. Méthode de Boltzmann sur réseau par différences finies (FDLBM) : Cette approche discrétise le terme convectif de la LBE en utilisant des schémas de différences finies. Plus précisément, elle emploie une différence directe de premier ordre pour les dérivées temporelles et un schéma amont (upwind) de premier ordre pour les dérives spatiales afin d'éviter les oscillations, bien que cela limite la précision spatiale au premier ordre près des frontières.
3. La méthode de Qadyan : Une nouvelle approche proposée par Ghadyani et al. qui combine des schémas semi-discrets avec la méthode de quadrature différentielle (DQM). Dans cette méthode, les dérivées spatiales de la LBE sont discrétisées en utilisant des coefficients de pondération DQM. Pour reproduire l'effet amont inhérent aux problèmes d'advection-diffusion, la méthode modifie les coefficients de pondération pour attribuer un poids plus important aux nœuds en amont qu'aux nœuds en aval. Cela permet de résoudre la LBE sous la forme d'un système d'équations différentielles ordinaires sans être lié à un type spécifique de modèle de vitesse discrète.

L'étude valide ces méthodes en utilisant la conduction thermique transitoire unidimensionnelle dans une plaque avec des conditions aux limites de Dirichlet (température constante) et de Neumann (isolée/flux nul). Des solutions analytiques dérivées de l'équation de la chaleur servent de références. Des normes d'erreur ($L_1$, $L_2$ et $L_\infty$) sont calculées pour évaluer la précision, et un problème de diffusion dans une plaque 2D est résolu pour démontrer l'applicabilité de la méthode aux cas multidimensionnels.

Contributions clés et résultats

• Précision vs Stabilité : Les résultats indiquent que si la LBM traditionnelle offre une grande précision, la méthode Qadyan proposée constitue une alternative robuste qui améliore la précision par rapport à l'approche FDLBM, particulièrement pour les conditions de Neumann. Pour les conditions de Dirichlet, la LBM D2Q5 standard a présenté les normes d'erreur les plus faibles, mais la méthode Qadyan a nettement surpassé les formulations FDLBM.
• Indépendance de la grille : La méthode Qadyan réussit à élargir la plage d'exactitude et de stabilité des solutions, augmentant la précision sans introduire de complexité supplémentaire à la méthode de Boltzmann sur réseau standard ou recourir à des grilles non uniformes.
• Coût computationnel : Une conclusion critique est l'équilibre entre précision et efficacité. La méthode Qadyan implique un coût de calcul plus élevé et un temps de simulation plus long par rapport à la LBM standard et à la FDLBM, en raison du nombre accru d'opérations de calcul requises pour les coefficients de pondération DQM.
• Sensibilité aux paramètres : L'étude établit que pour la LBM et la FDLBM standards, les valeurs de référence pour la vitesse de réseau dépendent du type de condition aux limites (température vs flux). En revanche, la méthode Qadyan utilise une valeur de vitesse de réseau constante ($c=0,015$) quel que soit le type de condition aux limites, simplifiant ainsi la sélection des paramètres.
• Validation 2D : Dans une simulation de plaque carrée bidimensionnelle, la méthode Qadyan a produit des contours isothermes et des distributions de température sur la ligne centrale qui concordent étroitement avec les résultats de la méthode des différences finies, démontrant sa capacité pour les problèmes de transfert de chaleur multidimensionnels.

Signification et affirmations
L'article affirme que la méthode de Qadyan est une approche appropriée pour une large gamme de problèmes de transfert de chaleur par conduction. Sa principale importance réside dans sa capacité à accroître la précision numérique sans introduire la complexité de grilles non uniformes ou abandonner le cadre de la LBM. Les auteurs soutiennent que la méthode est robuste et simple à implémenter. Cependant, l'article adopte une position modeste concernant l'efficacité, reconnaissant que la méthode est plus coûteuse en termes de calcul que les approches traditionnelles. L'étude conclut que bien que la méthode de Qadyan nécessite un temps de calcul plus long, son amélioration de la précision et de la stabilité pour des conditions aux limites spécifiques justifie son application, offrant une alternative compétitive aux formulations LBM standards pour les simulations de transfert de chaleur par conduction.