Aller au contenu principal
fr

Gradients et limiteurs sur un maillage médiocre

Dans un solveur centré sur les cellules, chaque flux repose sur une estimation de gradient. C'est là que la qualité du maillage entre en premier dans la solution.

Dans la méthode des volumes finis, les inconnues sont stockées aux centres des cellules, mais les flux sont nécessaires aux faces. Le pont entre les deux est l'estimation du gradient, et c'est là que la qualité du maillage entre en premier dans la solution — sur des cellules obliques, une erreur de gradient devient directement une erreur de flux.

Deux méthodes courantes

Moindres carrés. On ajuste un plan sur les centres des cellules voisines. Plus robuste sur maillages irréguliers ; elle faiblit lorsque les voisins sont regroupés dans une seule direction.

Green-Gauss. On intègre sur les valeurs aux faces. Économique et exacte sur maillage régulier ; sur des cellules obliques, l'interpolation de la valeur de face se dégrade et l'erreur grandit.

Le choix dépend de la qualité du maillage et fait partie des points à vérifier sur une solution qui ne converge pas. Voir la note sur les réglages du solveur.

Ce que protègent les limiteurs, et ce qu'ils coûtent

La reconstruction du second ordre produit des dépassements non physiques aux gradients raides : énergie cinétique turbulente négative, fractions volumiques au-delà de 100 %, températures sous le zéro absolu. Un limiteur les écrête.

Le coût est le suivant : partout où le limiteur intervient, la solution retombe localement au premier ordre et de la diffusion numérique est générée. Plus le limiteur est agressif, plus la solution est stable et floue.

Le point pratique

Regardez où le limiteur intervient. S'il intervient à une transition réellement raide — une onde de choc ou une surface libre — c'est normal. S'il intervient au milieu d'un écoulement régulier, c'est un symptôme de qualité de maillage : il y a dans cette zone des cellules obliques ou des sauts de volume brusques.

La plupart des cas où la qualité du maillage paraît correcte en moyenne alors que la solution se dégrade viennent exactement de là — ce qu'il faut examiner, ce sont les pires cellules. Voir la note sur le type de cellule.

Références

  1. Ferziger, J. H., Perić, M. & Street, R. L. (2020). Computational Methods for Fluid Dynamics, 4th ed. Springer. DOI: 10.1007/978-3-319-99693-6
  2. Patankar, S. V. (1980). Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing. ISBN 978-0-89116-522-4.

Discutons de votre projet

Définissons ensemble le système à analyser et la question à laquelle répondre.

Discutons de votre projet