Historia del Usuario
El siguiente ejemplo describe experimentos de túnel de viento realizados por el Laboratorio de Túnel de Viento Ambiental (EWTL) en la Universidad de Hamburgo [1] como un caso de validación en la Parte 9.4 del WTG-Merkblatt M3. Vamos a usar los campos de velocidad medidos del modelo de Michel City (Caso BL3-3) para validar simulaciones numéricas CFD en estructuras urbanas complejas. El ejemplo puede pertenecer al Grupo 1, según Figura 2.2 en WTG-Merkblatt-M3, basado en la investigación del valor promedio de 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 (por ejemplo, examen de túnel de viento de la clase G3).
- R2: Solitaria: todas las direcciones del viento relevantes con una resolución direccional suficientemente fina.
- Z2: Valores medios estadísticos y desviaciones estándar: siempre que involucren procesos de flujo estacionario para los cuales una verificación estadística de 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 amortiguación.
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 del 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 del perfil de α=0.27, representando 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, incluidos valores medios, varianzas, correlaciones y espectros, se obtuvieron con un Anemómetro Laser-Doppler 2D (LDA) a 500-600 Hz. Los puntos de medición se distribuyeron en perfiles verticales y horizontales, en cañones urbanos y en ubicaciones de repetibilidad definidas. El conjunto de datos de Michel City sirve como un caso de validación de referencia (C5) según la Directriz VDI 3783 Parte 9 [2]. Para la validación, además de la razón de impacto, se aplica una desviación relativa D=0.25 y una desviación absoluta W=0.08 para tener en cuenta la repetibilidad y la incertidumbre en la medición. Este conjunto de datos ha sido verificado y adoptado por varias instituciones (por ejemplo, KalWin [3]) para la validación CFD y propósitos de 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 simulación. El Método de la Tasa de Acierto 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 confiabilidad de la simulación calculando una tasa de acierto q, similar a las funciones de confianza utilizadas en la teoría de la confiabilidad. En contraste, el método del Error Cuadrático Medio Normalizado (e2) ofrece una evaluación de precisión más detallada al cuantificar la desviación cuadrática promedio 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 simulación.
Resultados y Discusión
La comparación entre los valores de velocidad normalizados (U/Uref) obtenidos de simulaciones RWIND y mediciones experimentales demuestra un nivel moderado de acuerdo en el conjunto de datos investigado. Se analizaron un total de 43 puntos de validación donde los valores de desviación varían aproximadamente del 2% a casi el 50%, indicando que, aunque la simulación captura la magnitud y tendencia general del campo de velocidad, las discrepancias locales siguen siendo significativas en ciertas regiones. El análisis de la tasa de acierto destaca aún más este comportamiento: solo el 18.60% de los puntos de datos caen dentro de una estricta tolerancia de ±10%, aumentando a 37.21% cuando la tolerancia se relaja a ±20%. El valor del error medio normalizado 𝑒2=0.2498 confirma una desviación general moderada entre los valores predichos y medidos.
Espacialmente, se observan desviaciones menores en regiones probablemente caracterizadas por condiciones de flujo más estables o adheridas, mientras que ocurren mayores discrepancias 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 en el modelado de turbulencia RANS estacionario, suposiciones de funciones de pared, efectos de resolución de malla o 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 malla o la definición de condiciones de contorno.
Tabla 1: Comparación de 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 Acierto q (tolerancia del 10%) | 18.60% |
| Tasa de Acierto q (tolerancia del 20%) | 37.21% |
| Forma de Desviación del Error Medio (e²) | 0.2498 |