Historia de usuario
El siguiente ejemplo describe experimentos en túnel de viento realizados por el Laboratorio de Túnel de Viento Ambiental (EWTL) de la Universidad de Hamburgo [1] como un caso de validación en la Parte 9.4 del WTG-Merkblatt M3. Vamos a utilizar los campos de velocidad medidos y los datos de rugosidad del modelo de la Ciudad de Michel (Caso BL3-3) para validar simulaciones numéricas de CFD en estructuras urbanas complejas. El ejemplo puede pertenecer al Grupo 2, de acuerdo con Figura 2.2 en WTG-Merkblatt-M3, basándose 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 planean investigaciones posteriores con mayor precisión (por ejemplo, examen en túnel de viento de clase G3).
- R2: Solitario: todas las direcciones de viento relevantes con 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 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 ensayo 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 que representa condiciones de terreno "muy rugoso". Se registraron un total de 1,838 puntos de medición para varias configuraciones de techos. Se obtuvieron los componentes de velocidad horizontal dependientes del tiempo u y v, incluyendo valores promedio, varianzas, correlaciones y espectros, con un Anemómetro Láser-Doppler 2D (LDA) a 500-600 Hz. Los puntos de medición se distribuyeron en perfiles verticales y horizontales, en cañones de calles y en ubicaciones de repetibilidad definidas. El conjunto de datos de la Ciudad de Michel 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 tasa de aciertos, 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 de medición. Este conjunto de datos ha sido verificado y adoptado por varias instituciones (por ejemplo, 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 simulación. El método de Tasa de Aciertos 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 error). 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 fiabilidad. En contraste, el método de Error Cuadrático Medio Normalizado (e2) ofrece una evaluación de precisión más detallada al cuantificar la desviación cuadrática media entre los valores simulados y los 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 simulaciones.
Resultados y discusión
La comparación entre los valores de velocidad normalizados (U/Uref) obtenidos de simulaciones RWIND y medidas 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 cuales los valores de desviación varían aproximadamente del 2% a casi el 50%, indicando que aunque la simulación captura la magnitud y la tendencia general del campo de velocidad, las discrepancias locales siguen siendo 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 caen dentro de una estricta tolerancia de ±10%, aumentando al 37.21% cuando se relaja la tolerancia a ±20%. El valor del error cuadrático medio normalizado 𝑒2=0.2498 confirma una desviación general moderada entre los valores predichos y los medidos.
Espacialmente, se observan desviaciones menores en regiones probablemente caracterizadas por condiciones de flujo más estables o adjuntas, mientras que ocurren discrepancias mayores donde se esperan gradientes más fuertes o fenómenos de flujo complejos, como zonas de separación o áreas de estela. Estas diferencias pueden atribuirse a limitaciones inherentes al 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 el análisis a nivel de ingeniería y la predicción de tendencias. Sin embargo, para aplicaciones que requieren una mayor precisión local, puede ser necesario un refinamiento adicional del modelado de turbulencia, la resolución de malla o la definición de condiciones de frontera.
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 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 |