Histoire d’utilisateur
L’exemple suivant décrit les essais en soufflerie menés par le Environmental Wind Tunnel Laboratory (EWTL) de l’Université de Hambourg [1] en tant que cas de validation dans la Partie 9.4 de la WTG-Merkblatt M3. Nous allons utiliser les champs de vitesse mesurés du modèle nommé « Michel City » (Cas BL3-3) pour valider les simulations numériques CFD pour des structures urbaines complexes. L’exemple peut appartenir au Groupe 1, selon Figure 2.2 dans WTG-Merkblatt-M3, basé sur l’étude de la valeur moyenne de la vitesse du vent :
- G2 : Valeurs absolues avec exigences de précision moyennes : le domaine d’application peut inclure des paramètres ou des études préliminaires lorsque des investigations ultérieures avec une plus grande précision sont prévues (par exemple, essai en soufflerie de classe G3).
- R2 : Solitaire : toutes les directions du vent pertinentes avec une résolution directionnelle suffisamment fine.
- Z2 : Valeurs moyennes statistiques et écarts-types : à condition d’impliquer des processus de flux stationnaires pour lesquels une vérification statistique des fluctuations avec un facteur de pointe est suffisante.
- S1 : Effets statiques : suffisants pour représenter le modèle structurel avec le détail mécanique nécessaire, mais sans propriétés de masse et d’amortissement.
Description
L’investigation porte sur un modèle idéalisé mais géométriquement détaillé d’une ville placé dans un écoulement de couche limite atmosphérique. Les essais en soufflerie ont été effectuées dans l’installation WOTAN, avec une section d’essai de 18 m de long, 4 m de large et 2,75-3,25 m de haut. Le champ de rugosité correspondant était caractérisé par une longueur de rugosité de z0=1,53 m et un exposant de profil α=0,27 représentant des conditions de terrain « très rugueuses ». Un total de 1 838 points de mesure ont été enregistrés pour plusieurs configurations de toiture. Les composantes de vitesse horizontale dépendantes du temps u et v, y compris les valeurs moyennes, les variances, les corrélations et les spectres, ont été obtenues avec un anémomètre Laser-Doppler 2D (LDA) à 500-600 Hz. Les points de mesure étaient répartis dans des profils verticaux et horizontaux, dans des canyons de rue, et à des emplacements de répétabilité définis. Le jeu de données de la ville de Michel sert de cas de validation de référence (C5) selon la directive VDI 3783 Partie 9 [1]. Pour la validation, en plus du taux de réussite, un écart relatif D=0,25 et un écart absolu W=0,08 sont appliqués pour prendre en compte la reproductibilité et l’incertitude de mesure. Ce jeu de données a été vérifié et adopté par plusieurs institutions (par exemple, KalWin [3]) pour la validation et la comparaison de modèles CFD.
Exigence de précision WTG-Merkblatt M3
WTG-Merkblatt M3 fournit deux méthodes principales pour valider les résultats de simulation. La méthode du taux de réussite évalue combien des 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 (e2) 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ée pour tenir compte des différences d’échelle. Ensemble, ces méthodes fournissent des mesures qualitatives et quantitatives pour la validation de simulation.
Résultats et discussion
La comparaison entre les valeurs de vitesse normalisées (U/Uref) obtenues à partir des simulations RWIND et des mesures expérimentales démontre un niveau de concordance modéré sur l’ensemble du jeu de données étudié. Un total de 43 points de validation ont été analysés, où les valeurs d’écart varient d’environ 2 % à presque 50 %, indiquant que bien que la simulation capture l’ordre de grandeur et la tendance générale du champ de vitesse, des écarts locaux restent significatifs dans certaines régions. L’analyse du taux de réussite souligne davantage ce comportement : seuls 18,60 % des points de données se trouvent dans une tolérance stricte de ±10 %, augmentant à 37,21 % lorsque la tolérance est relâchée à ±20 %. La valeur de l’erreur moyenne normalisée 𝑒2=0,2498 confirme un écart global modéré entre les valeurs prévues et mesurées.
Spatialement, des écarts plus faibles sont observés dans des régions probablement caractérisées par des conditions d’écoulement plus stables ou attachées, tandis que des écarts plus importants apparaissent là où des gradients plus forts ou des phénomènes d’écoulement complexes sont attendus, tels que des zones de séparation ou des régions de sillage. Ces différences peuvent être attribuées à des limitations inhérentes à la modélisation de la turbulence RANS des flux stationnaires, aux hypothèses de fonction de voile, aux effets de résolution de maillage ou à la sensibilité aux conditions d'entrée. Malgré ces limitations, la simulation démontre une performance adéquate pour une analyse au niveau de l’ingénierie et la prédiction de tendances. Cependant, pour des applications nécessitant une précision locale plus élevée, un raffinement plus important de la modélisation de la turbulence, de la résolution du maillage ou de la définition des conditions aux limites peut être requis.
Tableau 1 : Comparaison des vitesses normalisées (U/Uref) entre RWIND et les données expérimentales
| X | Y | Z | U/Uref – RWIND | U/Uref – Expérimental | Déviation (%) |
|---|---|---|---|---|---|
| -67.25 | -41.76 | 2 | 0.113 | 0.0942 | 16.64 |
| -70.41 | -27.84 | 2 | 0.120 | 0.0787 | 34.42 |
| -73.58 | -13.92 | 2 | 0.096 | 0.0693 | 27.81 |
| -76.75 | 0 | 2 | 0.098 | 0.0959 | 2.14 |
| -63.23 | 0 | 2 | 0.144 | 0.1305 | 9.38 |
| -61.71 | -12 | 2 | 0.134 | 0.1285 | 4.10 |
| -59.44 | -21.99 | 2 | 0.113 | 0.1399 | 23.81 |
| -57.17 | -31.98 | 2 | 0.124 | 0.1385 | 11.69 |
| -56.89 | -44.45 | 2 | 0.158 | 0.1186 | 24.94 |
| -46.54 | -47.15 | 2 | 0.149 | 0.1393 | 6.51 |
| -48.19 | -23.16 | 2 | 0.117 | 0.156 | 33.33 |
| -48.19 | -12 | 2 | 0.105 | 0.154 | 46.67 |
| -49.71 | 0 | 2 | 0.177 | 0.1347 | 23.90 |
| -36.19 | 0 | 2 | 0.117 | 0.1547 | 32.22 |
| -36.19 | -16.61 | 2 | 0.210 | 0.1633 | 22.24 |
| 31.81 | 127.75 | 2 | 0.091 | 0.086 | 5.49 |
| 58.81 | 119.75 | 2 | 0.096 | 0.093 | 3.13 |
| 72.31 | 127.75 | 2 | 0.065 | 0.0558 | 14.15 |
| 72.31 | 115.54 | 2 | 0.082 | 0.0423 | 48.41 |
| 85.81 | 118.91 | 2 | 0.149 | 0.0748 | 49.80 |
| 85.81 | 127.75 | 2 | 0.192 | 0.0978 | 49.06 |
| -149.89 | -124.58 | 2 | 0.103 | 0.1115 | 8.25 |
| -156.52 | -106.63 | 2 | 0.632 | 0.4036 | 36.15 |
| -142.43 | -99.38 | 2 | 0.658 | 0.334 | 49.24 |
| -141.58 | -112.44 | 2 | 0.476 | 0.4192 | 11.93 |
| -130.02 | -143.76 | 2 | 0.117 | 0.1723 | 47.26 |
| -119.53 | -120.2 | 2 | 0.313 | 0.384 | 22.68 |
| -127.31 | -105.09 | 2 | 0.634 | 0.3833 | 39.54 |
| -114.27 | -84.87 | 2 | 0.119 | 0.0894 | 24.87 |
| -105.26 | -112.85 | 2 | 0.466 | 0.3084 | 33.82 |
| -100.19 | -77.62 | 2 | 0.174 | 0.1187 | 31.78 |
| -36.27 | -94.25 | 2 | 0.483 | 0.2455 | 49.17 |
| -35.24 | -109.91 | 2 | 0.249 | 0.2782 | 11.73 |
| -48.97 | -103.35 | 2 | 0.366 | 0.2337 | 36.15 |
| -55.52 | -120.02 | 2 | 0.220 | 0.284 | 29.09 |
| -69.23 | -113.45 | 2 | 0.265 | 0.2553 | 3.66 |
| -103.49 | -58.25 | 2 | 0.270 | 0.1459 | 45.96 |
| -118.54 | 21.4 | 2 | 0.024 | 0.0224 | 6.67 |
| -121.46 | 20.73 | 2 | 0.062 | 0.0361 | 41.77 |
| -120.54 | 30.17 | 2 | 0.117 | 0.0952 | 18.63 |
| -117.61 | 30.84 | 2 | 0.139 | 0.0896 | 35.53 |
| -122.53 | 38.95 | 2 | 0.170 | 0.0907 | 46.65 |
| -125.45 | 38.28 | 2 | 0.158 | 0.1128 | 28.61 |
Tableau 2 : Métriques de validation pour la comparaison des vitesses normalisées (U/Uref)
| Métrique | Valeur |
|---|---|
| Nombre total de points de données (N) | 43 |
| Taux de réussite q (tolérance 10 %) | 18,60 % |
| Taux de réussite q (tolérance 20 %) | 37,21 % |
| Forme d’écart de l’erreur moyenne (e²) | 0,2498 |