Estudo da camada limite atmosférica em regiões metropolitanas costeiras com simulações de brisa marítima

O principal objetivo deste trabalho foi identificar e caracterizar a evolução diária da Camada Limite Atmosférica (CLA) na Região da Grande Vitória (RGV), Estado do Espírito Santo, Brasil e na Região de Dunkerque (RD), Departamento Nord Pas-de-Calais, França, avaliando a acurácia de parametrizações usadas no modelo meteorológico Weather Research and Forecasting (WRF) em detectar a formação e atributos da Camada Limite Interna (CLI) que é formada pelas brisas marítimas. A RGV tem relevo complexo, em uma região costeira de topografia acidentada e uma cadeia de montanhas paralela à costa. A RD tem relevo simples, em uma região costeira com pequenas ondulações que não chegam a ultrapassar 150 metros, ao longo do domínio de estudos. Para avaliar os resultados dos prognósticos feitos pelo modelo, foram utilizados os resultados de duas campanhas: uma realizada na cidade de Dunkerque, no norte da França, em Julho de 2009, utilizando um sistema light detection and ranging (LIDAR), um sonic detection and ranging (SODAR) e dados de uma estação meteorológica de superfície (EMS); outra realizada na cidade de Vitória – Espírito Santo, no mês de julho de 2012, também usando um LIDAR, um SODAR e dados de uma EMS. Foram realizadas simulações usando três esquemas de parametrizações para a CLA, dois de fechamento não local, Yonsei University (YSU) e Asymmetric Convective Model 2 (ACM2) e um de fechamento local, Mellor Yamada Janjic (MYJ) e dois esquemas de camada superficial do solo (CLS), Rapid Update Cycle (RUC) e Noah. Tanto para a RGV quanto para a RD, foram feitas simulações com as seis possíveis combinações das três parametrizações de CLA e as duas de CLS, para os períodos em que foram feitas as campanhas, usando quatro domínios aninhados, sendo os três maiores quadrados com dimensões laterais de 1863 km, 891 km e 297 km, grades de 27 km, 9 km e 3 km, respectivamente, e o domínio de estudo, com dimensões de 81 km na direção Norte-Sul e 63 km na Leste-Oeste, grade de 1 km, com 55 níveis verticais, até um máximo de, aproximadamente, 13.400 m, mais concentrados próximos ao solo. Os resultados deste trabalho mostraram que: a) dependendo da configuração adotada, o esforço computacional pode aumentar demasiadamente, sem que ocorra um grande aumento na acurácia dos resultados; b) para a RD, a simulação usando o conjunto de parametrizações MYJ para a CLA com a parametrização Noah produziu a melhor estimativa captando os fenômenos da CLI. As simulações usando as parametrizações ACM2 e YSU inferiram a entrada da brisa com atraso de até três horas; c) para a RGV, a simulação que usou as parametrizações YSU para a CLA em conjunto com a parametrização Noah para CLS foi a que conseguiu fazer melhores inferências sobre a CLI. Esses resultados sugerem a necessidade de avaliações prévias do esforço computacional necessário para determinadas configurações, e sobre a acurácia de conjuntos de parametrizações específicos para cada região pesquisada. As diferenças estão associadas com a capacidade das diferentes parametrizações em captar as informações superficiais provenientes das informações globais, essenciais para determinar a intensidade de mistura turbulenta vertical e temperatura superficial do solo, sugerindo que uma melhor representação do uso de solo é fundamental para melhorar as estimativas sobre a CLI e demais parâmetros usados por modelos de dispersão de poluentes atmosféricos.

Interacção de um corpo cúbico assente numa superfície muito rugosa com uma camada limite turbulenta

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia Mecânica na especialidade de Aerodinâmica pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Solução da equação de difusão unidimensional transiente para o estudo da dispersão de poluentes na camada limite planetária

Neste trabalho apresenta-se uma solução analítica para a dispersão vertical turbulenta em uma Camada Limite Convectiva e em uma Camada Limite Estável. A equação analisada considera a difusão com velocidades finitas, o que representa o transporte turbulento fisicamente correto. Considerando o caráter não-local, adicionam-se na equação que representa uma fonte área instantânea, termos como: o tempo de relaxação, a assimetria, a escala de tempo Lagrangeana e a velocidade turbulenta vertical. A solução é obtida utilizando-se a técnica da Transformada de Laplace. Os parâmetros que encerram a turbulência são derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade. Foram utilizados coeficientes de difusão especáficos para cada uma das camadas. A transformada inversa é obtida através do esquema numérico de quadratura Gaussiana. São apresentadas várias simulações para diferentes alturas de fonte área e obtém-se o valor da concentração para alturas próximas ao solo e próximas ao topo da Camada Limite Planetária. A inserção do termo de contra-gradiente na equação resultou em uma pequena influência na concentração de poluentes, observada de forma mais expressiva na Camada Limite Convectiva.

Estudo micrometeorológico da camada limite planetária estável na região de Porto Alegre

A análise micrometeorológica da camada limite planetária implica na obtenção de parâmetros de escala. Estas escalas típicas dão a noção da estrutura espacial e temporal da Camada Planetária, que estão diretamente ligada aos forçantes na superfície. Para tanto, necessita-se dos dados de perfis que caracterizam estes fluxos, os dados utilizados foram obtidos em 1995 e 1996. Estas radiossondagens foram realizadas no Aeroporto Salgado Filho na cidade de Porto Alegre, RS. Os perfis da temperatura e do cisalhamento do vento foram obtidos com uma resolução na vertical em torno de 50 metros. Com estes dados estimou-se a altura da camada de inversão através do perfil de temperatura potencial. A camada de mistura turbulenta estimou-se pelo cálculo do perfil do Número de Richardson. Como também, analisa-se o comportamento do jato de nível baixo, sua altura e direção. Procura-se reconhecer o tipo de resfriamento que é significante, radiativo ou turbulento. Com estes parâmetros típicos procura-se caracterizar a dinâmica da dispersão de poluentes dentro da camada.

Evolução de camada limite planetária para dispersão de poluentes pelo método da GILTT

O objetivo deste trabalho é obter os parâmetros turbulentos para o crescimento da camada limite planetária (CLP), durante a realizaçãoo do experimento Olad (Overland along wind dispersion experiment), conduzido na transição da noite para o dia. Nesta hora a CLP exibe uma altura, geralmente, pequena, disponibilizando pouco volume para a dispersão dos poluentes. Assim, concentrações superficiais elevadas podem ocorrer, atacando materiais, plantas e a saúde da população. Logo, conhecer os parâmetros do crescimneto é de fundamental importância para o correto modelamento da dispersão atmosférica ao amanhecer. A validação dos parâmetros é realizada a partir da solução da equação da difusão-advecção bidimensional, pelo método da GILTT (Generalized Integral Laplace Transform Technique). São empregados coeficientes de difusão turbulenta (problema de fechamento) dependentes da estabilidade atmosférica. As concentrações superficiais tridimensionais são obtidas através do espalhamento lateral da pluma com distribuição gaussiana. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com os dados experimentais. O modelo proposto mostrou-se aceitável em relação aos dados do experimento.

Análise quali-quantitativa do deslocamento da camada-limite em carrocerias veiculares

Dois dos principais objetivos no estudo da aerodinâmica de veículos é a segurança e o desempenho. Esses objetivos podem ser alcançados através do desenvolvimento de dispositivos que modifiquem o escoamento de ar ao redor do veículo e também através de alterações da forma das superfícies externas. O perfil mais baixo da frente do veículo, por exemplo, tem grande influência sobre o coeficiente de penetração aerodinâmica (Cx) e posteriormente sobre grande parte da distribuição de pressão ao longo da superfície do veículo. O objetivo deste trabalho foi de analisar o comportamento aerodinâmico que incide sobre 2 (dois) tipos de carrocerias de veículos, baseado em experimentação em túnel aerodinâmico, visando determinar a eficiência aerodinâmica de cada uma delas. Foram comparados dois modelos de veículos, sendo um deles com configuração familiar, muito comum nas grandes cidades e outro modelo com configuração esportiva mais voltado para o desempenho, ambos com escala 1:18. Para a realização do experimento foi construído um túnel aerodinâmico do tipo soprador com circuito aberto, onde foram analisadas as distribuições de pressão nas superfícies inferiores e superiores da carroceria dos modelos e o deslocamento da camada-limite. Tomando-se por base toda a teórica deste trabalho, com relação à aerodinâmica de carroceria, considerou-se nesta dissertação que essa filosofia de pesquisa é capaz de gerar reflexos positivos, que podem tornar os veículos cada vez mais eficientes e influenciar positivamente nos níveis de segurança e desempenho

Simulação numérica da camada limite planetária na região de Iperó, SP-Brasil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo comparativo da estrutura da turbulência atmosférica na camada limite superficial acima da floresta primária e região desmatada da Amazônia: análise de ondeletas aplicadas à detecção de estruturas correntes

Foram analisados dados turbulentos de resposta rápida, coletados em torres micrometeorológicas de 54 metros de altura em Caxiuanã-Pa, de 60 metros de altura na Rebio-Jarú-Ro e em torre de 4,04 metros de altura na Fazenda Nossa Senhora-Ro. Foi desenvolvido estudo comparativo de detecção de estruturas coerentes em dados obtidos acima de florestas e pastagem, portanto sob diferentes condições de rugosidade superficial. Aplicou-se a Transformada em Ondeletas de Morlet para identificar a escala de ocorrência das estruturas coerentes nos sinais de temperatura. Comparou-se a duração média e o comprimento médio das estruturas coerentes para duas.superfícies com rugosidades superficiais diferentes: Caxiuanã e Fazenda Nossa Senhora. Parâmetros de estabilidade atmosférica foram calculados para propiciar um estudo da variabilidade das estruturas coerentes em função das condições de estabilidade atmosférica. Foi calculado sobre a Rebio-Jarú a altura do ponto de inflexão do perfil vertical da velocidade do vento médio e a partir do valor do mesmo observou-se uma nítida correlação deste com a escala de ocorrência das estruturas coerentes. Finalmente, foram calculados os coeficientes de correlação para os sinais brutos e na escala de ocorrência das estruturas coerentes com a finalidade de melhor compreender os resultados de um ponto de vista da ação das estruturas coerentes na mistura turbulenta. Os resultados mostram um comportamento marcadamente diferente, no que diz respeito à floresta e à pastagem. Observou-se que em superfícies rugosas, tais com florestas, em condições próximas da neutralidade, predominam padrões de oscilações associados as escalas temporais da ordem de 80 s, as quais provavelmente, são manifestações de estruturas do tipo "rolos". No entanto, à medida que a altura do ponto de inflexão se aproximou do topo do dossel ocorreu uma diminuição na duração das estruturas coerentes. Além disso, foi observado que o coeficiente de correlação calculado na escala de ocorrência das estruturas coerentes foi consideravelmente maior que o mesmo calculado para o sinal bruto. Tais resultados sugerem mudanças significativas na estrutura da turbulência em regiões desmatadas comparativamente àquelas em que foram mantidas as florestas primárias da Amazônia.

Análise de um cilindro submetido à excitação acústica por camada limite turbulenta

A launching rocket and its payload are submitted to strong acoustic loads in some moments of its flight like lift off, during transonic flight and in the instant of maximum dynamic pressure. These loads could affect the payload and other rocket internal equipment. It must be taken into account that in the higher stages, where delicate control equipment and the payload are placed, the acoustic excitations are dominant. The knowledge of the acoustic loads is of great importance in order to provide estimated loads for the design phase, define intensity level to specify qualification and acceptation acoustic tests to which subsystems and elements of the launcher and its payload shall be submitted and to design an acoustic … (Complete abstract click electronic access below)

Influência de efeitos não-locais na dispersão de poluentes na camada linite planetária

A equação de difusão-advecção é muito utilizada no campo de estudos da poluição atmosférica na determinação da concentração de poluentes. Uma maneira de solucionar o problema de fechamento desta equação está baseada na hipótese de transporte por gradiente que, em analogia com a difusão molecular, assume que o fluxo turbulento de concentração é proporcional à magnitude do gradiente de concentração média. Neste trabalho, diferentemente do modo tradicional, utiliza-se uma equação genérica para a difusão turbulenta considerando-se que o fluxo mais a sua derivada são proporcionais ao gradiente médio. Desta forma, obtém-se uma equação que leva em conta a assimetria no processo de dispersão de poluentes atmosféricos. Portanto, a proposta do presente trabalho é a obtenção da solução analítica desta nova equação utilizando-se a técnica da Transformada de Laplace, considerando-se a Camada Limite Planetária (CLP) como um sistema multicamadas. Os parâmetros que encerram a turbulência sâo derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade convectiva válidos para grandes tempos de difusão. Finalmente, na avaliação da performance deste modelo que considera a assimetria no processo de dispersão de poluentes atmosféricos, utilizam-se os dados experimentais de Copenhagen e Prairie Grass.

Transição à turbulência na camada de mistura estavelmente estratificada utilizando simulação numérica direta e simulação de grandes escalas

A transição à turbulência em uma camada de mistura estavelmente estratificada é de grande interesse para uma variedade de problemas geofísicos e de engenharia. Esta transição é controlada pela competição entre o cisalhamento do escoamento de base e as forças de empuxo, devido à estratificação em densidade do ambiente. Os efeitos do empuxo atuam no escoamento reduzindo a taxa de crescimento das perturbações e retardando a transição à turbulência, enquanto o cisalhamento fornece energia cinética ao escoamento. O presente trabalho investiga a natureza da transição à turbulência em uma camada de mistura temporal estavelmente estratificada através de Simulação Numérica Direta (DNS) e Simulação de Grandes Escalas (LES). O propósito da investigação é analisar o efeito da estratificação estável no desenvolvimento da instabilidade de Kelvin-Helmholtz (K-H) e na formação dos vórtices longitudinais, que se formam após a saturação dos turbilhões primários de K-H. Além deste propósito, é examinado, utilizando de DNS e LES, o desenvolvimento das instabilidades secundárias de K-H na camada baroclínica. Os testes numéricos tridimensionais são realizados com diferentes tipos de condições iniciais para a flutuação de velocidade transversal, enquanto uma condição forçada é usada para as outras duas componentes de flutuação de velocidade. Em particular, o efeito do comprimento transversal do domínio de cálculo é testado empregando diferentes comprimentos, enquanto são usadas as mesmas dimensões para a direção longitudinal e vertical. As simulações bidimensionais mostram que o aumento da estratificação inibe o processo de emparelhamento, reduz a troca de energia entre os turbilhões de K-H e o escoamento, atenua a instabilidade de K-H e diminui o fluxo vertical de massa. A instabilidade secundária do tipo K-H é identtificada na camada baroclínica para Re ¸ 500 quando há o processo de emparelhamento dos vórtices simulados. Na simulação a Re = 500, a instabilidade secundária de K-H aparece tanto para Ri = 0.07 (fraca estratificação) como para Ri = 0.167 (forte estratificação). Os resultados tridimensionais demonstram que os vórtices longitudinais são claramente formados na camada a Ri = 0. Por outro lado, nos casos estratificados os vórtices são enfraquecidos, devido ao gradiente longitudinal de densidade, que diminui a vorticidade nos turbilhões de K-H enquanto aumenta na região entre eles.

Uma solução da equação difusão-advecção com o termo contragradiente

Neste trabalho, apresenta­se uma solução para equação de difusão­advecção considerando o  termo  contragradiente  que  é um termo  adicional. Esse termo  adicional  contém informações  sobre a assimetria, escala de tempo Lagrangeana e velocidade turbulenta vertical. A solução  da equação foi obtida pela utilização da técnica de Transformada de Laplace, considerando a  Camada  Limite  Planetária  (CLP)  como  um  sistema  de  multicamadas.  Os  parâmetros  turbulentos foram derivados da teoria de difusão estatística de Taylor, combinada com a teoria  da similaridade. Assim, são apresentadas simulações para diferentes valores de assimetria, o  que  propiciou  a  obtenção  de  uma  concentração  de  contaminantes  em  diferentes  alturas,  em  uma  camada  limite  convectiva.  A  avaliação  do  desempenho  do  modelo,  considerando  a  assimetria  no  processo  de  dispersão  de  poluentes  atmosféricos, foi realizada  através  de  um  experimento  de  tanque  convectivo  tradicional.  Nesse  experimento,  o  termo  contragradiente  influenciou a concentração de poluentes para uma camada limite convectiva. Entretanto, com  as  parametrizações  utilizadas,  o  modelo  não  conseguiu  captar  de  forma  eficiente  o  comportamento da concentração em pontos mais distantes da fonte.

Interacção térmica de edíficios altos com o vento atmosférico

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do Grau de Mestre em Engenharia Mecânica.

Estudo laboratorial de escoamentos em leitos com vegetação emersa e rígida

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do grau de Mestre em Engenharia Civil,Perfil de Construção

Efeito do Vento como Ação Dinâmica em Edifícios Altos

Este trabalho aborda uma série de conceitos base no que concerne à ação do vento sobre edifícios altos, começando por ser estabelecidas algumas considerações fundamentais acerca da circulação do vento na camada limite atmosférica bem como acerca da sua interação com as estruturas. É feita uma análise da metodologia proposta pelo Eurocódigo 1 para quantificação de tal ação sobre os edifícios, bem como é elaborada uma comparação da metodologia proposta por este com a metodologia ainda vigente na regulamentação portuguesa. Foram modelados computacionalmente, com recurso a um programa de cálculo estrutural automático, três edifícios altos com diferente secção geométrica em planta que servirão de caso de estudo. Para estes mesmos edifícios são aplicados os dois regulamentos considerados com vista à determinação de esforços e deslocamentos. Sendo os edifícios altos um género de estruturas capazes de ser excitadas dinamicamente perante a ação do vento, adota-se uma metodologia para quantificação desta ação de forma dinâmica na direção do escoamento. Assim, é obtida a resposta dinâmica ao longo do tempo em termos de deslocamentos e acelerações para o caso de estudo considerado e é feita uma comparação da resposta do edifício quadrangular sob a ação dinâmica do vento com a resposta estática regulamentar.