Données de Mesure pour Charges de Vent sur Structures Cylindriques

User Story

Cet exemple présente des mesures expérimentales des forces aérodynamiques et de la distribution de pression sur un cylindre circulaire, qui sont largement utilisées comme données de référence pour valider les simulations CFD en ingénierie du vent. L'écoulement autour d'un cylindre circulaire représente un problème aérodynamique classique où la séparation de l'écoulement, la formation du sillage et les effets du nombre de Reynolds influencent fortement les forces aérodynamiques. En raison de ces phénomènes, les surfaces courbes telles que les cylindres sont particulièrement difficiles à appréhender pour les simulations numériques.

Les études expérimentales montrent que les coefficients aérodynamiques d'un cylindre varient considérablement avec le nombre de Reynolds et la rugosité de la surface. Aux nombres de Reynolds élevés typiques des écoulements de vent atmosphérique, les mesures présentent souvent une dispersion considérable, indiquant que les résultats dépendent non seulement du nombre de Reynolds mais aussi des caractéristiques de la surface et des conditions de turbulence. L'exemple peut appartenir au **Groupe 1**, selon Figure 2.2 du WTG-Merkblatt-M3, sur la base de l'investigation de la valeur moyenne de la vitesse du vent :

• G1 : Valeurs qualitatives avec des exigences de faible précision pour une utilisation dans l'investigation de base ou la conception préliminaire. L'effort et les exigences relatives au niveau de détail sont réduits, car souvent toutes les conditions aux limites ne sont pas complètement clarifiées.
• R1 : Solitaire (sans bâtiments environnants), analyse des directions de vent importantes individuelles.
• Z1 : Valeurs moyennes statistiques, à condition qu'elles concernent des processus d'écoulement stationnaires où les fluctuations (par exemple, dues à la turbulence de l'écoulement approchant) peuvent être suffisamment capturées par d'autres mesures.
• S1 : Effets statiques. Il suffit de représenter le modèle structural avec le détail mécanique nécessaire, mais sans les propriétés de masse et d'amortissement.

Description

Cet exemple de validation considère l'écoulement autour d'un cylindre circulaire fini monté verticalement sur une plaque plane sous une vitesse d'entrée uniforme, fournissant un banc d'essai simplifié et bien contrôlé pour les simulations CFD. Le cylindre a un diamètre de D = 30 mm et une hauteur de L = 180 mm, et il est positionné dans un domaine de type soufflerie où l'écoulement entrant est aligné avec la direction de l'écoulement (X), tandis que Y et Z désignent respectivement les directions latérale et verticale. Contrairement aux simulations de couche limite atmosphérique (CLA), la condition aux limites d'entrée est définie comme un profil de vitesse constant, U(z)=U_o, avec une vitesse d'écoulement libre de 𝑈_o=10, ce qui signifie qu'aucun cisaillement vertical ou gradient de vitesse n'est présent. Cette hypothèse supprime la complexité associée au développement de la couche limite et permet une investigation ciblée de la physique fondamentale de l'écoulement autour du cylindre.

Dans les conditions d'écoulement données, le nombre de Reynolds basé sur le diamètre du cylindre est d'environ 𝑅𝑒≈20 000, plaçant l'écoulement dans un régime subcritique où la séparation et le détachement tourbillonnaire sont attendus. L'écoulement entrant uniforme conduit à une région de stagnation symétrique à l'avant du cylindre, suivie d'une séparation de l'écoulement le long des côtés et de la formation d'un sillage caractérisé par une recirculation et un détachement tourbillonnaire périodique (dans les simulations transitoires). L'interaction entre le cylindre et la plaque plane introduit une complexité supplémentaire par le biais des effets de proche paroi et du développement d'une couche limite au sol, en fonction de la résolution numérique et de l'approche de traitement de paroi.

1 Montage expérimentale et système de coordonnées (mm)

Tableau 1 : Données d'entrée du cylindre circulaire

Étude de maillage computationnel

La figure 2 présente une analyse de sensibilité de maillage d'un modèle cylindrique dans RWIND. Le coefficient de force calculé (C_f) diminue légèrement de 0,76 à une densité de maillage de 15 % à 0,71 à 25 %, et ensuite à 0,70 à 35 %. Cette réduction graduelle indique que la solution se stabilise à mesure que le maillage est raffiné. La faible variation de C_f à des densités de maillage plus élevées démontre une convergence globale, suggérant qu'un raffinement supplémentaire n'a qu'un impact mineur sur les résultats.

2 Analyse de sensibilité du maillage du modèle cylindrique montrant la convergence du coefficient de force (Cf) avec l'augmentation de la densité du maillage

De plus, l'étude de maillage computationnel doit être effectuée conformément au lien suivant :

• Étude de maillage numérique

Exigence de précision du WTG-Merkblatt M3

Le WTG-Merkblatt M3 fournit deux méthodes clés pour valider les résultats de simulation. La méthode du **Taux de réussite (Hit Rate)** évalue combien de valeurs simulées Pi correspondent correctement aux valeurs de référence Oi dans une tolérance définie, en utilisant une approche de classification binaire (succès ou échec). Cette approche évalue la fiabilité de la simulation en calculant un taux de réussite q, similaire aux fonctions de confiance utilisées en théorie de la fiabilité. En revanche, la méthode de l'**Erreur quadratique moyenne normalisée** **(e²)** offre une évaluation de précision plus détaillée en quantifiant l'écart quadratique moyen entre les valeurs simulées et de référence, normalisé pour tenir compte des différences d'échelle. Ensemble, ces méthodes fournissent des mesures à la fois qualitatives et quantitatives pour la validation de la simulation.

• Validation des résultats de simulation à l’aide de la méthode du taux de réussite définie dans la WTG-Merkblatt M3

• Évaluation de la précision des résultats à l'aide de l'erreur quadratique moyenne normalisée définie dans la WTG-Merkblatt M3

Résultats et discussion

La figure 3 compare le coefficient de pression moyen (Cp) de RWIND et les données expérimentales pour différentes directions du vent. Les deux montrent une forte concordance dans les tendances globales, y compris la stagnation à 0°, la chute de pression due à la séparation de l'écoulement et la récupération à des angles plus élevés. Le Cp minimum se produit autour de 60°–70° et est bien prédit par RWIND. Cependant, RWIND sous-estime légèrement la dépression dans la région du sillage à des directions de vent plus élevées. Dans l'ensemble, les résultats indiquent une bonne précision globale avec des écarts mineurs.

3 Coefficient de pression moyen (Cp) vs direction du vent

La comparaison entre RWIND et les coefficients de pression expérimentaux montre un très bon accord pour les directions de vent entre 0° et 100°, où les écarts restent inférieurs à 10 % et tous les points de données satisfont aux critères de taux de réussite de 10 % et 20 %. Les écarts les plus importants se produisent dans la région du sillage (110°–180°), où RWIND sous-estime systématiquement la dépression, conduisant à des écarts plus élevés allant jusqu'à environ 43 %.

Dans l'ensemble, 58 % des résultats se situent dans une plage d'écart de 10 % (et 20 %). Les métriques d'erreur globales (e² = 0,05, MAE = 0,12, RMSE = 0,16) indiquent une précision prédictive acceptable, l'écart étant principalement concentré dans la région d'écoulement séparé.

Tableau 2 : Comparaison du coefficient de pression (Cp) entre RWIND et les données expérimentales

0x

0x

Exemple WTG 9.5 – Structure cylindrique