Dynamique des fluidesProcédés et eau

Mécanique des fluides numérique

Résoudre une équation aux dérivées partielles sur une grille — un champ de vitesse dans l’espace et le temps. Un solveur Lattice-Boltzmann autonome reproduit le résultat le plus emblématique de la dynamique des fluides : l’allée de tourbillons de von Kármán émise par un écoulement autour d’un cylindre, validée sur un écoulement en canal, le nombre de Strouhal, la traînée et la portance, et une cavité entraînée par un couvercle.

DjiniousLabModélisation et simulation

Un champ de vorticité montrant des tourbillons rouges et bleus alternés se détachant d’un cylindre en une allée de von Kármán

erreur de l’écoulement en canal contre l’analytique: 0.27%erreur de l’écoulement en canal contre l’analytique
nombre de Strouhal (fréquence de détachement): 0.20nombre de Strouhal (fréquence de détachement)
coefficient de traînée du cylindre: 0.94coefficient de traînée du cylindre
la traînée oscille au double de la fréquence de portance: 2×la traînée oscille au double de la fréquence de portance
notebooks de conception: 10notebooks de conception
exigences RÉUSSIES: 5 / 5exigences RÉUSSIES

Un champ de vitesse entier, pas une poignée d’états.

Chaque autre modèle de ce portefeuille suit quelques variables dans le temps. L’écoulement fluide est différent : c’est un champ — vitesse et pression en tout point de l’espace, couplés par une équation aux dérivées partielles, évoluant ensemble. Le résoudre sur une grille, c’est la mécanique des fluides numérique, le domaine de COMSOL et d’ANSYS Fluent. Ce programme construit un solveur autonome à partir de zéro et l’utilise pour reproduire le résultat le plus reconnaissable du domaine — l’allée de tourbillons qui se détache d’un cylindre dans un écoulement — puis le valide face à des références de manuel.

Bien comprendre la physique sur un problème que vous pouvez vérifier.

Le solveur est une méthode de Lattice-Boltzmann — une voie cinétique, à distribution de particules, vers les équations de Navier-Stokes incompressibles, qui est explicite, locale et robuste là où une résolution directe de la pression est délicate. Avant de lui faire confiance sur quoi que ce soit de complexe, le programme le valide sur un écoulement en canal sous gradient de pression, où la réponse exacte est une parabole : le profil simulé correspond à l’analytique à 0,27 % près. Un solveur qui réussit le cas que vous pouvez vérifier sur papier est un solveur auquel vous pouvez croire sur le cas que vous ne pouvez pas vérifier.

Un profil de vitesse d’écoulement en canal simulé superposé à la parabole analytique — L’ancrage de justesse : l’écoulement en canal sous gradient de pression (Poiseuille) a un profil de vitesse parabolique exact, et le solveur Lattice-Boltzmann le reproduit avec une erreur L2 de 0,27 %. Validez d’abord face à la forme close, puis lâchez le solveur sur les écoulements qui n’en ont pas.

Le champ de vorticité d’une allée de tourbillons de von Kármán — des tourbillons rouges et bleus alternés se détachant d’un cylindre — Le héros : un écoulement autour d’un cylindre au-delà d’un nombre de Reynolds critique détache des tourbillons alternés — l’allée de tourbillons de von Kármán. Le champ de vorticité montre le motif caractéristique de tourbillons contrarotatifs (rouge et bleu) se détachant du cylindre et convectés en aval. Le résultat le plus reconnaissable de la dynamique des fluides, résolu sur une grille.

Le sillage ne se détache pas toujours — jusqu’à ce qu’il le fasse.

À faible nombre de Reynolds, l’écoulement derrière le cylindre est une paire stable et symétrique de bulles de recirculation. Augmentez le nombre de Reynolds au-delà d’une valeur critique et cette symétrie se brise : le sillage devient instable et commence à détacher des tourbillons à une fréquence périodique nette — une bifurcation de Hopf dans l’écoulement lui-même. Le programme balaye le nombre de Reynolds et capture la transition, l’amplitude d’oscillation du sillage passant d’essentiellement nulle à un détachement pleinement développé.

Amplitude d’oscillation du sillage en fonction du nombre de Reynolds — proche de zéro à faible Re, importante au-dessus de la valeur critique — La transition de détachement : en dessous d’un nombre de Reynolds critique (~47), le sillage est stable et l’amplitude d’oscillation est essentiellement nulle ; au-dessus, le sillage se détache et l’amplitude bondit. Une bifurcation dans l’écoulement, capturée en balayant un seul nombre sans dimension.

Le détachement laisse une empreinte sur les forces.

L’allée de tourbillons n’est pas seulement esthétique — elle est mesurable. Sondez le sillage et sa vitesse oscille à une fréquence qui, rendue sans dimension, est le nombre de Strouhal ≈ 0,2, la valeur que cite chaque manuel de mécanique des fluides. Et le détachement secoue le cylindre : la force de portance oscille à la fréquence de détachement tandis que la traînée ondule à exactement le double — la signature de von Kármán, et la raison pour laquelle cheminées, câbles et tabliers de pont doivent être conçus contre les vibrations induites par les tourbillons.

Traînée et portance sur le cylindre au fil du temps — la portance oscillant à la fréquence de détachement, la traînée au double — La signature de force : les tourbillons alternés poussent le cylindre d’un côté puis de l’autre, si bien que la portance oscille à la fréquence de détachement et que la traînée ondule à exactement le double (coefficient de traînée ≈ 0,94). C’est la vibration induite par les tourbillons — la raison pour laquelle les structures élancées dans un écoulement sont conçues précisément contre cela.

L’autre test canonique, également réussi.

La cavité entraînée par un couvercle — un fluide dans une boîte entraîné par un couvercle supérieur mobile — est l’autre cas de validation CFD standard, avec des décennies de solutions de référence publiées. Le solveur produit le tourbillon de recirculation principal attendu, son centre se plaçant exactement là où la littérature (Ghia et al.) le situe. Deux références canoniques réussies, à partir d’un seul solveur autonome écrit dans le notebook.

Un écoulement en cavité entraînée par un couvercle — un tourbillon de recirculation principal avec des flèches de vitesse circulant dans la boîte — La cavité entraînée par un couvercle : un couvercle supérieur mobile entraîne le fluide en un seul tourbillon de recirculation principal, son centre se plaçant là où les solutions de référence le situent. L’autre référence classique du domaine CFD — même solveur, une boîte fermée au lieu d’un sillage ouvert.

Chaque chiffre redérivé à la validation finale.

Le notebook de V&V reconstruit le solveur depuis zéro et redérive chaque exigence face à la théorie et aux références, en imprimant un tableau RÉUSSI/ÉCHEC.

Résultat

Profil de Poiseuille contre analytique : 0,27 %
Amplitude de détachement tourbillonnaire / U₀ : 0,87
Nombre de Strouhal : 0,20
Coefficient de traînée / signature : 0,94 ; traînée = 2× portance
Centre du tourbillon principal de cavité : (0,61, 0,72)

Exigence

Profil de Poiseuille contre analytique : R-01 < 2 %
Amplitude de détachement tourbillonnaire / U₀ : R-02
Nombre de Strouhal : R-03 ≈ 0,15–0,2
Coefficient de traînée / signature : R-04
Centre du tourbillon principal de cavité : R-05

Un véritable solveur de champ, à fidélité pédagogique.

Le livrable est dix notebooks et le dossier — un champ EDP discrétisé, pas un réseau acausal localisé, donc pas de bloc personnalisé ni de canevas. Le solveur est un Lattice-Boltzmann 2D incompressible, à faible nombre de Mach, avec une collision à temps de relaxation unique ; la résolution de grille est modeste, le cylindre utilise une frontière de rebond en escalier, les coefficients de traînée sont indicatifs plutôt que validés pour publication, et il n’y a pas de modèle de turbulence. Ce que le programme démontre est la capacité qui manquait le plus au portefeuille — résoudre une équation aux dérivées partielles comme un champ spatial sur une grille — avec un solveur autonome (sans dépendance CFD externe) qui reproduit l’allée de tourbillons, le nombre de Strouhal, la signature de force, et deux références canoniques, chaque chiffre réexécutable.