Caso de usuario
El siguiente ejemplo describe los experimentos en túnel de viento llevados a cabo por el Laboratorio de Túnel de Viento Ambiental (Environmental Wind Tunnel Laboratory, EWTL) de la Universidad de Hamburgo [1] como caso de validación en la Parte 9.4 del WTG-Merkblatt M3. Vamos a utilizar los campos de velocidad medidos del modelo de la Ciudad de Michel (Caso BL3-3) para validar las simulaciones numéricas de CFD en estructuras urbanas complejas. El ejemplo puede pertenecer al Grupo 1, según Figura 2.2 en el WTG-Merkblatt-M3, basado en la investigación del valor promedio de la velocidad del viento:
- G2: Valores absolutos con requisitos de precisión media: El área de aplicación puede incluir parámetros o estudios preliminares cuando se planifican investigaciones posteriores con mayor precisión (p. ej., examen en túnel de viento de clase G3).
- R2: Aislado: todas las direcciones de viento relevantes con una resolución direccional suficientemente fina.
- Z2: Valores medios estadísticos y desviaciones estándar: siempre que impliquen procesos de flujo estacionario para los cuales una verificación estadística de las fluctuaciones con un factor de pico sea suficiente.
- S1: Efectos estáticos: Son suficientes para representar el modelo estructural con el detalle mecánico necesario, pero sin propiedades de masa y amortiguamiento.
Descripción
La investigación se centra en un modelo de ciudad idealizado pero geométricamente detallado, colocado en un flujo de capa límite atmosférica. Las mediciones en el túnel de viento se realizaron en la instalación WOTAN, que cuenta con una sección de prueba de 18 m de largo, 4 m de ancho y 2.75-3.25 m de alto. El campo de rugosidad correspondiente se caracterizó por una longitud de rugosidad de z0=1.53 m y un exponente de perfil α=0.27, lo que representa condiciones de terreno "muy rugoso". Se registraron un total de 1,838 puntos de medición para varias configuraciones de techo. Los componentes de velocidad horizontal dependientes del tiempo u y v, incluyendo valores medios, varianzas, correlaciones y espectros, se obtuvieron con un Anemómetro Láser-Doppler (LDA) 2D a 500-600 Hz. Los puntos de medición se distribuyeron en perfiles verticales y horizontales, en calles estrechas y en ubicaciones de repetibilidad definidas. El conjunto de datos de la Ciudad de Michel sirve como caso de validación de referencia (C5) según la Directriz VDI 3783 Parte 9 [2]. Para la validación, además de la tasa de aciertos, se aplican una desviación relativa D=0.25 y una desviación absoluta W=0.08 para tener en cuenta la repetibilidad y la incertidumbre de la medición. Este conjunto de datos ha sido verificado y adoptado por varias instituciones (p. ej., KalWin [3]) para fines de validación de CFD y comparación de modelos.
Requisito de precisión del WTG-Merkblatt M3
El WTG-Merkblatt M3 proporciona dos métodos clave para validar los resultados de la simulación. El método de la Tasa de Aciertos (Hit Rate) evalúa cuántos de los valores simulados Pi coinciden correctamente con los valores de referencia Oi dentro de una tolerancia definida, utilizando un enfoque de clasificación binaria (acierto o fallo). Este enfoque evalúa la fiabilidad de la simulación calculando una tasa de aciertos q, similar a las funciones de confianza utilizadas en la teoría de la fiabilidad. Por el contrario, el método del Error Cuadrático Medio Normalizado (e2) ofrece una evaluación de la precisión más detallada al cuantificar la desviación cuadrática media entre los valores simulados y de referencia, normalizada para tener en cuenta las diferencias de escala. Juntos, estos métodos proporcionan medidas tanto cualitativas como cuantitativas para la validación de la simulación.
Resultados y discusión
La comparación entre los valores de velocidad normalizada (U/Uref) obtenidos de las simulaciones de RWIND y las mediciones experimentales demuestra un nivel moderado de concordancia en todo el conjunto de datos investigado. Se analizaron un total de 43 puntos de validación en los que los valores de desviación oscilan entre aproximadamente el 2% y casi el 50%, lo que indica que, si bien la simulación captura la magnitud y la tendencia generales del campo de velocidad, persisten discrepancias locales significativas en ciertas regiones. El análisis de la tasa de aciertos destaca aún más este comportamiento: solo el 18.60% de los puntos de datos se encuentran dentro de una tolerancia estricta de ±10%, aumentando al 37.21% cuando la tolerancia se relaja a ±20%. El valor del error medio normalizado 𝑒2=0.2498 confirma una desviación global moderada entre los valores predichos y los medidos.
Espacialmente, se observan desviaciones más bajas en regiones probablemente caracterizadas por condiciones de flujo más estables o adherido, mientras que surgen discrepancias mayores donde se esperan gradientes más fuertes o fenómenos de flujo complejos, como zonas de separación o regiones de estela. Estas diferencias pueden atribuirse a limitaciones inherentes al modelado de turbulencia RANS estacionario, las suposiciones de las funciones de pared, los efectos de la resolución de la malla o la sensibilidad a las condiciones de entrada. A pesar de estas limitaciones, la simulación demuestra un rendimiento adecuado para análisis a nivel de ingeniería y predicción de tendencias. Sin embargo, para aplicaciones que requieren una mayor precisión local, puede ser necesario un mayor refinamiento del modelado de turbulencia, la resolución de la malla o la definición de las condiciones de contorno.
Tabla 1: Comparación de la velocidad normalizada (U/Uref) entre RWIND y datos experimentales
| X | Y | Z | U/Uref – RWIND | U/Uref – Experimental | Desviación (%) |
|---|---|---|---|---|---|
| -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 |
Tabla 2: Métricas de validación para la comparación de velocidad normalizada (U/Uref)
| Métrica | Valor |
|---|---|
| Número total de puntos de datos (N) | 43 |
| Tasa de Aciertos q (tolerancia 10%) | 18.60% |
| Tasa de Aciertos q (tolerancia 20%) | 37.21% |
| Forma de desviación del error medio (e²) | 0.2498 |