Un solveur CFD ne calcule pas de forces. Il calcule des champs de pression et de vitesse ; les forces en sont déduites, par intégration sur les surfaces. Comprendre cette étape change la façon de lire un résultat de traînée portance CFD — et permet d'attraper les erreurs les plus courantes.

• Prérequis : un calcul terminé, sur une aile ou une hélice ;
• Durée : une demi-heure de lecture ;
• Ce que vous en tirez : savoir ce que recouvre exactement chaque nombre du tableau de résultats.

• Déroulé du didacticiel :

1. De la pression à la force
2. Pression et frottement, deux parts à distinguer
3. Cas aile : portance, traînée, Cx et Cz
4. Cas hélice : poussée et couple
5. Lire le tableau des forces
6. Les erreurs d'interprétation classiques

1 : De la pression à la force

Sur chaque petite face du maillage qui touche l'objet, le fluide exerce deux actions : une pression, perpendiculaire à la face, et une contrainte de cisaillement, tangente à la face. Multipliées par l'aire de la face, elles donnent une force élémentaire.

La force totale est la somme de ces contributions sur toutes les faces d'une frontière donnée. Le moment se calcule de la même façon, en multipliant chaque force élémentaire par son bras de levier par rapport à l'axe choisi.

Cette mécanique a une conséquence pratique : la précision des forces dépend directement de la qualité du maillage de surface. Une pale décrite par des faces trop grandes intègre une pression moyennée, et une couche limite mal maillée fausse la contrainte de cisaillement — donc toute la traînée de frottement.

Sur chaque face : une pression perpendiculaire et un cisaillement tangent. Les forces sont la somme de tout cela

2 : Pression et frottement, deux parts à distinguer

Heliciel PRO CFD/BEM sépare, pour chaque frontière, la part due à la pression et la part visqueuse. Cette séparation n'est pas cosmétique : elle est le premier outil de diagnostic d'un résultat.

• Sur une pale, la poussée doit venir massivement de la pression. C'est la différence entre intrados et extrados qui travaille ;
• La part visqueuse représente le frottement. Elle est petite sur une pale bien profilée, et elle est toujours résistante ;
• Sur un moyeu ou une paroi de conduit, une part visqueuse importante signale un frottement parasite. Si elle domine le moment axial, la puissance qu'elle consomme n'est pas de la propulsion, c'est de la perte.

Le frottement du moyeu et celui de la paroi de conduit ne se règlent pas depuis l'interface : Heliciel décide selon le type de cas — moyeu glissant en conduit fermé, moyeu frottant sur un rotor libre — et inscrit son choix dans le heliciel_benchmark.txt du cas, où il reste lisible. Un cas rouvert conserve la valeur qu'il avait. C'est un exemple typique de ce que les préréglages prennent en charge : une décision dont une erreur fausserait le couple sans jamais faire échouer le calcul.

3 : Cas aile — portance, traînée, Cx et Cz

Pour un cas aile, l'écoulement va vers +X, et la décomposition est directe : la traînée est la composante de la force selon X, la portance est la composante selon Z. Les coefficients s'en déduisent en divisant par la pression dynamique et la surface de référence.

La traînée elle-même se décompose en deux termes que la CFD sépare naturellement :

• la traînée de pression, ou traînée de forme, liée à la distribution de pression autour du profil et à d'éventuels décollements ;
• la traînée de frottement, intégrale du cisaillement pariétal.

Attention à la convention : cette décomposition en traînée et portance n'a de sens que pour les cas aile, seuls consommateurs de ces grandeurs. Pour une hélice, l'écoulement va vers −X et la lecture est différente.

Les coefficients d'une aile, avec la part de pression et la part de frottement séparées

4 : Cas hélice — poussée et couple

Pour une hélice, deux grandeurs comptent : la force axiale selon l'axe de rotation, qui est la poussée ou la traînée du rotor selon qu'il propulse ou capte ; et le moment autour de ce même axe, qui est le couple.

Le couple multiplié par la vitesse angulaire donne la puissance à l'arbre. C'est le point de jonction avec le dimensionnement mécanique : moteur, réducteur, arbre, structure de pale.

Le signe distingue les deux familles de machines. Une hélice propulsive consomme un couple ; un rotor capteur — éolienne, hydrolienne, turbine — en produit un. Une inversion de signe inattendue s'explique presque toujours par la convention d'axes, jamais par la physique : rappelez-vous que le flux va vers −X pour une hélice tandis que l'axe de rotation reste +X.

Un dernier point sur ce qui est intégré : la CFD somme sur toutes les frontières solides déclarées — pales, moyeu, et selon le cas carénage ou objet inséré. Le BEM d'Heliciel, lui, rend les forces des pales. Ne comparez pas un total à une partie.

Les forces d'un rotor, frontière par frontière : pales, moyeu, carénage. Le total n'est pas la pale

5 : Lire le tableau des forces

Le bouton Matrix Me! de l'onglet Résultats ouvre le tableau numérique du cas : les grandeurs aérodynamiques, frontière par frontière, avec leur décomposition.

En Phase 2, une précision de méthode s'impose. Les forces instationnaires oscillent à la fréquence de passage des pales ; la valeur qui fait résultat est la moyenne sur le dernier tour, pas la valeur du dernier pas de temps. Citer un instant, c'est citer le sommet ou le creux d'une oscillation.

Cette même convergence des forces sert de critère d'arrêt au calcul : le moniteur compare la moyenne sur deux fenêtres glissantes et arrête proprement quand la variation relative passe sous son seuil. Forces et convergence sont donc le même sujet, vu de deux côtés — le didacticiel n°7 traite l'autre côté.

Les forces CFD en regard des forces BEM : ce que chaque total recouvre doit être su avant de comparer

6 : Les erreurs d'interprétation classiques

Quatre pièges, tous rencontrés en pratique :

• Comparer un total CFD à une pale BEM. Le moyeu traîne et frotte : son inclusion suffit à créer un écart qu'on attribuera à tort à la méthode ;
• Oublier la surface de référence. Un coefficient sans sa surface de référence n'est pas un coefficient, c'est un nombre ;
• Citer une force instantanée de Phase 2 en la présentant comme un résultat stationnaire ;
• Ignorer la diffusion numérique. Un schéma décentré du premier ordre — celui du préréglage Stable, et celui de la Phase 2 — adoucit les gradients et tend à surestimer la traînée. C'est le prix de la stabilité, il faut le savoir quand on compare deux calculs conduits avec des schémas différents.

La suite naturelle : le didacticiel n°12, qui applique tout ceci à un profil isolé et le confronte aux polaires de la base d'Heliciel.

La collection des vingt didacticiels CFD

« Précédent : Lire une carte de pression  |  Suivant : Vérifier une polaire de profil »

1. Première simulation CFD d'hélice
2. Lire le dossier de cas OpenFOAM
3. Choisir son préréglage de maillage
4. La zone MRF et la Phase 1
5. La Phase 2 en maillage glissant
6. Confronter le BEM et la CFD
7. Lire la convergence d'un calcul
8. Les modèles de turbulence
9. Couche limite et y+
10. Lire une carte de pression
11. D'où viennent les forces (vous êtes ici)
12. Vérifier une polaire de profil
13. Hélice marine et cavitation
14. Éolienne et hydrolienne
15. Ventilateur en conduit fermé
16. Hélice d'avion, croisière et point fixe
17. La boucle d'optimisation
18. Se caler sur une référence
19. De la carène à l'hélice
20. Le dossier de calcul livrable

Cette série accompagne les didacticiels de conception Heliciel, qui traitent la partie BEM : conception de la pale, choix du régime, courbes de performance. La CFD vient après eux, pour vérifier et pour voir.