Histoire d'utilisateur
L'exemple suivant décrit les expériences en soufflerie menées par le laboratoire 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 de la ville de Michel (Cas BL3-3) pour valider les simulations numériques CFD dans 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 enquêtes ultérieures avec une plus grande précision sont prévues (par exemple, examen 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 qu'ils impliquent des processus d'écoulement stationnaires pour lesquels une vérification statistique des fluctuations avec un facteur de pointe est suffisante.
- S1 : Effets statiques : Ils sont 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'étude se concentre sur un modèle idéalé mais géométriquement détaillé d'une ville placé dans un écoulement de couche limite atmosphérique. Les mesures en soufflerie ont été effectuées dans l'installation WOTAN, avec une section de test 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 toit. 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 [2]. Pour la validation, en plus du taux de réussite, une déviation relative D=0,25 et une déviation absolue W=0,08 sont appliquées pour prendre en compte la répétabilité 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
Le 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 à l'intérieur d'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 la déviation quadratique moyenne 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 d'accord 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 de déviation varient d'environ 2% à presque 50%, indiquant que bien que la simulation capte 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 une déviation globale modérée entre les valeurs prédites et mesurées.
Spatialement, des déviations plus faibles sont observées 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 stationnaire, aux hypothèses de fonction de paroi, 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 tendance. Cependant, pour des applications nécessitant une précision locale plus élevée, un affinement plus poussé de la modélisation de la turbulence, de la résolution du maillage ou de la définition des conditions aux limites peut être nécessaire.
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 de déviation de l'erreur moyenne (e²) | 0.2498 |