Les principes de l'affinement des mailles dans le CFD

Le processus consiste à augmenter la densité cellulaire dans les régions où les caractéristiques du flux nécessitent une résolution plus élevée pour capturer les gradients physiques avec précision. Les équations Navier-Stokes, qui régissent le mouvement du fluide, sont discrétées sur une grille de calcul composée de milliers à des centaines de millions de cellules. La qualité de la solution dépend fondamentalement de la façon dont cette grille résout les variations spatiales de vitesse, de pression et de température qui conduisent à la montée et à la descente.

Dans l'aérodynamique externe, le débit est rarement uniforme. De grands gradients de pression apparaissent aux bords de tête et de piste des aérofoils, les couches limites se développent le long de surfaces solides et les sillages se forment en aval des corps. Un maillage grossier peut représenter une aile entière avec seulement quelques cellules couvrant son accord, ce qui rend impossible de capturer la courbure simplifiée, la croissance de la couche limite ou le début de la séparation du flux.

Les ingénieurs identifient les zones avec des gradients importants attendus : points de stagnation, couches limites, bouts d'ailes et régions sensibles à l'éviction du vortex. Dans ces endroits, la taille des cellules est réduite jusqu'à ce que la solution devienne insensible à l'amélioration ultérieure — condition vérifiée par une étude d'indépendance du maillage. Sans une telle étude, les résultats du CFD ne peuvent être considérés comme entièrement fiables pour les décisions d'ingénierie.

Types de raffinage de mesh

Les trois familles classiques sont le h-affinement, le p-affinement et le r-affinement, chacune ayant des forces et des faiblesses distinctes qui influencent leur applicabilité à la prédiction de la force aérodynamique.

  • h-refinement: La stratégie la plus courante, le h-refinement ajoute plus de cellules dans les zones critiques tout en maintenant l'ordre polynôme des fonctions de forme fixes. Les outils modernes de maillage appliquent le h-refinement isotrope ou anisotrope selon que les caractéristiques de flux nécessitent une résolution égale dans toutes les directions ou seulement dans des orientations spécifiques, comme la direction mur-normale dans les couches limites.
  • p-refinement: Augmente l'ordre polynôme des fonctions d'interpolation au sein de chaque élément au lieu d'ajouter des cellules. Cela peut fournir une convergence exponentielle pour des problèmes lisses et est souvent utilisé dans les méthodes d'éléments finis à haut ordre. p-refinement est attrayant pour résoudre les caractéristiques d'écoulement à haut ordre sans gonfler la taille du maillage, mais exige des résolveurs robustes et peut lutter contre les discontinuités telles que les ondes de choc.
  • r-affinement: Relocalise les nœuds de mailles existants vers des régions à forte erreur sans changer le nombre total de cellules. Bien que calculablement moins cher que l'ajout de cellules, la qualité des cellules peut se dégrader si les cellules deviennent trop déformées ou tendues.

Dans les applications aérodynamiques, le raffinage h domine, souvent combiné avec des couches de prisme proches de la paroi qui se transforment en un noyau tétraédrique ou hexaédrique. Le mélange de cellules structurées et non structurées permet aux ingénieurs de capturer les couches limites en détail tout en maintenant le nombre global de cellules gérables. Le choix entre le raffinement isotrope et anisotrope dépend de la physique du flux; pour les couches limites attachées, les cellules anisotropes étirées le long de la direction de la paroi sont efficaces, tandis que pour les régions séparées, il peut être nécessaire de cellules isotropes plus nombreuses pour résoudre les structures vortiques tridimensionnelles.

Lier la résolution de mesh à la prévision de la force aérodynamique

Même de petites erreurs dans ces distributions peuvent se combiner en erreurs significatives dans les coefficients de force. Un maillage grossier pourrait endommager le pic d'aspiration sur la surface supérieure d'une aile, réduisant ainsi le levage prévu. Il peut déplacer le point de séparation, ce qui peut faire apparaître le sillage artificiellement large et augmenter la traînée de pression. Pour les flux transoniques, une onde de choc mal résolue peut conduire à un nombre de Mach complètement incorrect de traînée, ce qui peut entraîner une perte de performance de la conception d'un aéronef par une large marge.

La résistance visqueuse à la traînée est particulièrement sensible à la résolution des mailles. La traînée de frottement de la peau provient du gradient de vitesse à la paroi, et la précision dépend de l'existence d'un maillage suffisamment fin dans la direction de la paroi normale. La distance sans dimension y+ est la mesure standard : pour la simulation de grande taille par construction (LES) ou la simulation numérique directe (DNS), une y+ d'environ 1 est nécessaire, tandis qu'une approche de navigation-réveil (RANS) avec un modèle de turbulence à faible nombre de réynolds vise souvent y+dessous 5. Si le premier centroïde cellulaire est trop éloigné du mur, le résolveur compte sur des fonctions de paroi qui ne représentent pas le véritable flux si la couche limite n'est pas entièrement attachée.

La vitesse de la traînée de pression est tout aussi exigeante. La bonne capture de la séparation de débit autour des corps de bluff, du train d'atterrissage ou de la partie arrière d'un wagon dépend de la capacité de la grille à résoudre le gradient de pression défavorable et l'éviction subséquente du vortex. Une maille trop grossière dans le sillage entraînera la déstabilisation numérique de la bulle de séparation ou sera complètement supprimée, modifiant radicalement la pression de base et le coefficient de traînée de pression.

Comment différents éléments de la force réagissent à l'affinement

Le levage est largement régi par des forces de pression et tend à être moins sensible à la résolution du treillis que la traînée, à condition que la circulation générale autour du corps soit correctement capturée. Cela conduit parfois les praticiens à conclure prématurément qu'un treillis est adéquat parce que le soulèvement prévu correspond aux données du soufflet, alors que la traînée reste mal résolue. La traînée induite, liée au vortex à bout d'aile, nécessite une résolution de travée suffisante et un réveil correctement résolu pour calculer correctement le lavage du fond. Si la grille dans le sillage est trop grossière, la structure vortique diffuse artificiellement et la traînée induite peut être sous-estimée de plusieurs pour cent. La traînée parasitaire, qui comprend la friction cutanée et la traînée de pression due à la séparation, est très dépendante du mesh et nécessite généralement la meilleure résolution.

Un choc doit être capté à travers le moins de cellules possible pour maintenir la netteté, mais ces cellules doivent être suffisamment petites pour éviter les frottis qui élargiraient artificiellement le choc et modifieraient la distribution de la pression. Le raffinement adaptatif du maillage (RMA) excelle souvent ici, agglomérant automatiquement des cellules près du choc une fois qu'une solution préliminaire est obtenue. Sans ces capacités d'adaptation, les ingénieurs doivent placer manuellement des zones de raffinement basées sur l'expérience, qui peut être longue et sujette à erreur.

Études de cas : Amélioration des mailles en action

L'un des plus cités est le modèle de recherche commun (CRM) de la NASA, conçu pour valider les prédictions aérodynamiques fondées sur le RANS. Dans les ateliers de prévision de la dérive de l'AIAA, les participants ont présenté des résultats CFD pour une configuration transonique du corps des ailes à l'aide de diverses mailles. La dispersion dans la traînée totale prévue entre les groupes — même ceux utilisant le même modèle de solveur et de turbulence — était souvent des centaines de comptages de traînées, particulièrement près du nombre de Mach de la dérive. L'analyse a révélé que des parties substantielles de cette dispersion pouvaient être directement tracées par des différences de densité de mailles dans la région de choc et l'espacement près du mur.

Dans le sport automobile, le corps ahmède, forme automobile simplifiée à surface arrière inclinée, est un repère classique pour l'aérodynamique du corps bluff. Le débit sur le dos incliné est dominé par la séparation et les tourbillons longitudinaux forts. Lorsque le maillage est trop grossier dans la région inclinée, la paire de vortex contre-rotatif est souvent complètement oubliée ou fortement affaiblie, ce qui entraîne une surestimation de la pression de base et une sous-prédiction de la traînée de 20 %. Le affinage du maillage sur le maillage et dans le sillage proche révèle la structure du vortex complet et apporte le coefficient de traînée en ligne avec les expériences.

Un troisième cas instructif provient de simulations de rotors d'hélicoptère. Le vortex à pointe de pale, qui interagit avec la lame suivante causant une charge impulsive et du bruit, nécessite des mailles exceptionnellement fines dans le sillage pour éviter une diffusion prématurée. Les ingénieurs utilisent souvent des mailles surdimensionnées ou un raffinement adaptatif qui suit le vortex à pointe, ajoutant des cellules seulement lorsque nécessaire.

Meilleures pratiques pour la génération et le raffinage des meshs

Compte tenu de l'effet profond de la qualité des mailles sur les prédictions de force, un workflow discipliné est essentiel. Le point de départ est une étude d'indépendance des mailles, où le même cas est exécuté sur des grilles systématiquement raffinées. La méthodologie de l'indice de convergence des grilles (ICG), proposée par Roache, fournit une mesure quantitative de la distance entre une solution discrète et la valeur asymptotique.

Pour réaliser une étude d'indépendance de maillage appropriée, il faut au moins trois mailles : grossières, moyennes et fines. Chaque maille doit être affinée par un facteur constant dans chaque dimension, idéalement autour de √2 (1.41) pour un raffinement isotrope, ou par une réduction équivalente du volume cellulaire. Le résolveur doit être exécuté avec des réglages identiques sur tous les maillages, et les coefficients de force résultants sont tracés sur une mesure de l'espacement des mailles (comme la racine cube du nombre de cellules). L'extrapolation à la taille zéro des cellules donne la valeur asymptotique autour de laquelle les solutions doivent se regrouper si le maillage se situe dans la plage asymptotique.

L'étude d'indépendance du maillage doit examiner non seulement les coefficients de force globaux, mais aussi les distributions locales de la pression et de la friction cutanée. L'élévation et la traînée peuvent apparaître presque constantes alors que la distribution de la pression à une station spandwise spécifique continue de se déplacer avec raffinement. Ces variations locales peuvent signaler une séparation non résolue ou un mouvement de choc qui pourrait devenir critique dans des conditions hors-conception.

Mesures de qualité à surveiller

Au-delà du comptage cellulaire, plusieurs mesures de qualité influencent directement la précision des prédictions de force.

  • Orthogonalité et inclinaison: Les cellules près de la paroi devraient être aussi proches que possible de rectangulaires; une grande inclinaison introduit la diffusion numérique et dégrade la prédiction de la contrainte de cisaillement. Un angle d'orthogonalité inférieur à 10 degrés est souhaitable, en particulier dans la couche limite.
  • Rapport d'aspect: Dans la couche limite, les cellules doivent avoir des rapports d'aspect très élevés (souvent supérieurs à 10 000) pour résoudre le gradient mur-normal sans exploser le nombre de cellules dans le sens du cours d'eau. La stabilité du solveur pour les cellules à haut rapport d'aspect doit être vérifiée, car certains schémas deviennent rigides ou nécessitent des nombres de Courant inférieurs.
  • Ratio de transition de taille de cellule:[ Un saut soudain dans la taille des cellules peut créer des réflexions fallacieuses ou des erreurs d'interpolation. Les cellules adjacentes ne devraient pas différer en volume de plus de 20 % pour maintenir des gradients lisses.
  • Changement de volume :[ Les cellules volumétriques négatives ou les éléments fortement déformés causent une défaillance du résolveur. Un contrôle approfondi avant l'exécution est obligatoire, en utilisant des diagnostics intégrés dans le logiciel de maillage.
  • La variation du rapport d'aspect:[ Même si les cellules individuelles ont des rapports d'aspect acceptables, une variation rapide d'une cellule à l'autre peut introduire des erreurs numériques. La transition doit être progressive, en particulier au-delà du bord de la couche limite.

Ces paramètres devraient être vérifiés pour le maillage final et, idéalement, surveillés pendant le processus de raffinement. De nombreux générateurs commerciaux de maillage fournissent des rapports de qualité qui indiquent les régions problématiques pour la correction manuelle.

Erreur quantifiante : l'indice de convergence de grille en pratique

L'indice de convergence des grilles (ICG) est une méthode largement utilisée pour rapporter les études de convergence des grilles. Il a été développé par Patrick Roache et est recommandé par l'American Society of Mechanical Engineers (ASME) et l'AIAA pour la vérification CFD. L'ICG fournit une façon uniforme de comparer les erreurs de discrétisation de différentes études. Il est calculé à partir des résultats de trois mailles à l'aide de la formule:

GCIfin = Fs

où ε est le changement relatif entre les solutions fines et moyennes, r est le rapport de raffinement, p est l'ordre de précision observé, et Fs est un facteur de sécurité (généralement 1,25 pour les études à trois mailles). Un ICG inférieur à 2–5% indique que le maillage est suffisamment raffiné pour la quantité d'intérêt. Toutefois, l'ordre de précision observé s'écarte souvent de l'ordre théorique en raison de la qualité du maillage, de la non-linéarité de la solution ou de l'affinement insuffisant du maillage.

Dans la pratique, les ingénieurs devraient indiquer les valeurs de l'ICG pour les coefficients de levage et de traînée dans toute publication ou révision de conception. Cette transparence permet à d'autres d'évaluer la crédibilité des résultats.

Coût de calcul vs Fidelité: Équilibrer

Chaque augmentation de la densité de mailles est coûteuse. Dans une simulation RANS, le nombre de cellules influence directement la taille du système linéaire et la mémoire requise. Alors que les résolveurs modernes s'échellent bien sur les grappes HPC, le doublement du nombre de cellules peut augmenter le temps de calcul de plus d'un facteur de deux en raison de la nature non linéaire du processus de solution itérative. Pour les simulations dépendantes du temps comme la simulation Eddy Détachée (DES) ou la simulation Grand Eddy (LES), le coût est encore plus élevé car l'étape de temps doit être réduite proportionnellement à la plus petite taille de cellules pour maintenir un nombre de Courant-Friedrichs-Lewy (CFL) sous l'unité.

La tâche de l'ingénieur consiste à trouver le maillage le plus grossier qui donne encore suffisamment de précision pour la décision à venir. Ce seuil est spécifique à l'application. Dans la conception conceptuelle, où des centaines de configurations doivent être évaluées, un maillage qui produit des coefficients de traînée dans un rayon de 5% de la valeur attendue peut être acceptable.

Les processus de travail en scripts permettent d'appliquer la même topologie de maille à différentes géométries, assurant la cohérence et réduisant le temps sur la réparation manuelle de surface. Combinés à une estimation d'erreur basée sur l'interface, ces outils indiquent où le perfectionnement ultérieur permettra d'améliorer les prévisions de force, guidant l'ingénieur vers un maillage optimal avec un minimum d'essai et d'erreur. Dans bien des cas, un maillage de 20 millions de cellules qui est affiné de façon optimale peut fournir une meilleure précision qu'un maillage uniforme de 50 millions de cellules à une fraction du coût.

Techniques avancées de raffinage du mesh adaptatif

Cependant, le raffinement statique du maillage (AMR) – où le résolveur modifie automatiquement la grille en fonction d'un indicateur d'erreur dérivé de la solution – gagne rapidement en traction tant dans le milieu universitaire que dans l'industrie. Le RMA peut être basé sur des caractéristiques, des cellules de raffinage où l'amplitude du gradient de vitesse, de la vorticité ou de la pression dépasse un seuil. Ou bien, il peut être basé sur l'adjoint, ciblant directement l'erreur dans une sortie fonctionnelle telle que le levage ou le traîné. Le RMA basé sur l'adjoint est particulièrement puissant parce qu'il se concentre sur les régions qui influencent la quantité d'intérêt, évitant les cellules gaspillées dans des zones qui n'affectent pas la force finale. Cette approche peut réduire le nombre de cellules d'un ordre de grandeur tout en maintenant la même précision dans la prédiction du traînée.

Dans la pratique, l'AMR a été utilisé avec succès pour résoudre les jets supersoniques, les interactions entre les couches de choc et les couches de vortex instables avec des réductions spectaculaires du nombre total de cellules par rapport aux mailles uniformément fines. Pour une aile transonique, un parcours AMR à base d'un joint peut automatiquement détecter et affiner le pied de choc, la bulle de séparation et le vortex à bout d'aile, tout en laissant le champ éloigné grossier. Le résultat est une prédiction de traînée qui rivalise avec une grille plus fine uniforme mais à un dixième du coût de calcul.

Le rôle de la modélisation de la turbulence dans la sensibilité au mesh

Dans les simulations RANS, le modèle de turbulence fournit une fermeture pour les contraintes de Reynolds, et la précision dépend de la capacité du modèle à représenter correctement la production, la dissipation et le transport de turbulences. Le maillage compte parce que de nombreux modèles, en particulier ceux basés sur la formulation k‐ш (comme SST k‐ш), sont extrêmement sensibles à la résolution près du mur. Si le premier centre cellulaire se trouve dans la couche tampon (y+ entre 5 et 30) au lieu de la sous-couche visqueuse (y+ en dessous de 5), le modèle passe à une approche de fonction murale qui peut ne pas être appropriée pour les flux avec des gradients de pression ou une séparation forte.

Les méthodes hybrides RANS-LES, telles que la simulation de l'Eddy à distance (DDES) ou la simulation de l'Eddy à stress (SBES), imposent des exigences encore plus grandes sur le maillage. Ces méthodes reposent sur la grille pour déclencher la transition du mode RANS au mode ERP dans des régions séparées. Si le maillage n'est pas suffisamment fin dans ces régions, le modèle reste en mode RANS et les grandes structures turbulentes ne sont pas résolues, niant l'avantage de l'approche hybride.

Dans tous les cas, un maillage optimisé pour un modèle de turbulence peut ne pas convenir à un autre. Lors de la comparaison des résultats de différents modèles — ou lorsque l'analyse comparative des données expérimentales — la sensibilité au maillage devrait être revue. Une bonne pratique consiste à consulter les meilleures pratiques communautaires, telles que la NASA Langley Turbulence Modeling Resource[, qui fournit des conseils sur les exigences en maillage pour divers modèles.

Conclusion et perspectives d'avenir

L'efficacité du raffinement calculal du maillage dans la simulation des forces de levage et de traînée n'est pas simplement « plus de cellules égale plus de précision ». Il s'agit d'un jeu nuancé entre la résolution, l'erreur numérique, la modélisation de turbulences et le budget de calcul. Lorsqu'il est appliqué avec une méthodologie claire – en commençant par une étude d'indépendance du maillage utilisant GCI, en respectant les meilleures pratiques pour la discrétisation des couches limites, la surveillance des mesures de qualité et l'exploitation des techniques d'adaptation lorsque c'est possible – le raffinement fournit constamment des prédictions de plus haute fidélité qui permettent aux ingénieurs de repousser les limites de la performance aérodynamique.

Plusieurs tendances promettent de réduire le fardeau de la génération et du raffinement des mailles. Les indicateurs d'erreur basés sur l'apprentissage automatique peuvent un jour prédire la distribution optimale des cellules avec une entrée minimale de l'utilisateur, en tirant des leçons des simulations précédentes pour guider le raffinement. Les résolveurs accélérés GPU permettront d'exécuter des mailles de centaines de millions de cellules sur des postes de travail de classe bureau, réduisant ainsi le temps nécessaire aux études de haute fidélité.

En attendant, les fondamentaux restent inchangés. Un maillage bien construit est le fondement sur lequel se fondent des prédictions aérodynamiques fiables. Négligence de cette fondation — en coupant des coins sur la résolution près du mur, en ignorant l'indépendance du maillage, ou en sur-relyant sur le maillage automatisé sans contrôle de qualité — conduira à des prédictions de force qui semblent raisonnables mais sont physiquement erronées. L'investissement dans le raffinement soigneux maillage est remboursé à plusieurs reprises sous la forme d'avions qui répondent aux garanties de performance, de voitures de course qui restent collées à la piste, et de turbines qui fonctionnent efficacement pendant des décennies.