Caso de Utilizador
O exemplo seguinte descreve os ensaios em túnel de vento realizados pelo Laboratório de Túnel de Vento Ambiental (EWTL) da Universidade de Hamburgo [1] como um caso de validação na Parte 9.4 do Caderno WTG M3. Vamos utilizar os campos de velocidade medidos do modelo Michel City (Caso BL3-3) para validar simulações numéricas de CFD em estruturas urbanas complexas. O exemplo pode pertencer ao Grupo 1, de acordo com Figura 2.2 no Caderno WTG M3, com base na investigação do valor médio da velocidade do vento:
- G2: Valores absolutos com requisitos de precisão média: A área de aplicação pode incluir estudos paramétricos ou preliminares quando estão planeadas investigações posteriores com maior precisão (por exemplo, exame em túnel de vento da classe G3).
- R2: Isolado: todas as direções de vento relevantes com resolução direcional suficientemente fina.
- Z2: Valores médios estatísticos e desvios-padrão: desde que envolvam processos de escoamento estacionário para os quais uma verificação estatística das flutuações com um fator de pico seja suficiente.
- S1: Efeitos estáticos: São suficientes para representar o modelo estrutural com o detalhe mecânico necessário, mas sem propriedades de massa e amortecimento.
Descrição
A investigação foca-se num modelo de cidade idealizado, mas geometricamente detalhado, colocado num escoamento de camada limite atmosférica. As medições em túnel de vento foram realizadas na instalação WOTAN, que possui uma secção de teste com 18 m de comprimento, 4 m de largura e 2,75-3,25 m de altura. O campo de rugosidade correspondente foi caracterizado por um comprimento de rugosidade de z0=1,53 m e um expoente de perfil α=0,27 representando condições de terreno "muito rugoso". Foram registados um total de 1.838 pontos de medição para várias configurações de cobertura. As componentes horizontais da velocidade dependentes do tempo, u e v, incluindo valores médios, variâncias, correlações e espectros, foram obtidas com um Anemómetro Laser-Doppler (LDA) 2D a 500-600 Hz. Os pontos de medição foram distribuídos em perfis verticais e horizontais, em desfiladeiros urbanos e em locais definidos de repetibilidade. O conjunto de dados Michel City serve como um caso de validação de referência (C5) de acordo com a Diretriz VDI 3783 Parte 9 [2]. Para validação, além da taxa de acerto, são aplicados um desvio relativo D=0,25 e um desvio absoluto W=0,08 para ter em conta a repetibilidade e a incerteza de medição. Este conjunto de dados foi verificado e adotado por várias instituições (por exemplo, KalWin [3]) para fins de validação de CFD e comparação de modelos.
Requisito de Precisão do Caderno WTG M3
O Caderno WTG M3 fornece dois métodos principais para validar resultados de simulação. O método da Taxa de Acerto (Hit Rate) avalia quantos dos valores simulados Pi correspondem corretamente aos valores de referência Oi dentro de uma tolerância definida, usando uma abordagem de classificação binária (acerto ou erro). Esta abordagem avalia a fiabilidade da simulação calculando uma taxa de acerto q, semelhante às funções de confiança utilizadas na teoria da fiabilidade. Em contraste, o método do Erro Quadrático Médio Normalizado (e2) oferece uma avaliação de precisão mais detalhada ao quantificar o desvio quadrático médio entre os valores simulados e de referência, normalizado para ter em conta as diferenças de escala. Juntos, estes métodos fornecem medidas qualitativas e quantitativas para a validação da simulação.
Resultados e Discussão
A comparação entre os valores de velocidade normalizada (U/Uref) obtidos das simulações RWIND e as medições experimentais demonstra um nível moderado de concordância em todo o conjunto de dados investigado. Foram analisados um total de 43 pontos de validação nos quais os valores de desvio variam de aproximadamente 2% a quase 50%, indicando que, embora a simulação capte a magnitude geral e a tendência do campo de velocidade, permanecem discrepâncias locais significativas em certas regiões. A análise da taxa de acerto destaca ainda mais este comportamento: apenas 18,60% dos pontos de dados se enquadram dentro de uma tolerância estrita de ±10%, aumentando para 37,21% quando a tolerância é relaxada para ±20%. O valor do erro médio normalizado 𝑒2=0,2498 confirma um desvio global moderado entre os valores previstos e medidos.
Espacialmente, desvios mais baixos são observados em regiões provavelmente caracterizadas por condições de escoamento mais estáveis ou colado, enquanto discrepâncias maiores ocorrem onde são esperados gradientes mais fortes ou fenómenos de escoamento complexos, como zonas de separação ou regiões de esteira. Estas diferenças podem ser atribuídas a limitações inerentes à modelação de turbulência RANS estacionária, pressupostos de funções de parede, efeitos de resolução da malha ou sensibilidade às condições de entrada. Apesar destas limitações, a simulação demonstra um desempenho adequado para análises de nível de engenharia e previsão de tendências. No entanto, para aplicações que exigem maior precisão local, pode ser necessário um maior refinamento da modelação de turbulência, da resolução da malha ou da definição das condições de fronteira.
Tabela 1: Comparação da Velocidade Normalizada (U/Uref) Entre Dados Experimentais e RWIND
| X | Y | Z | U/Uref – RWIND | U/Uref – Experimental | Desvio (%) |
|---|---|---|---|---|---|
| -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 |
Tabela 2: Métricas de Validação para a Comparação da Velocidade Normalizada (U/Uref)
| Métrica | Valor |
|---|---|
| Número total de pontos de dados (N) | 43 |
| Taxa de Acerto q (tolerância de 10%) | 18,60% |
| Taxa de Acerto q (tolerância de 20%) | 37,21% |
| Forma de Desvio do Erro Médio (e²) | 0,2498 |