Search

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior The PBL is the region on the troposphere characterized by existence of turbulence (BUSINGER, 1981). Turbulent flow is that when the fluctuations of the variables that describe the flow are amplified in time and space (LANDAU; LIFSHITZ, 1959). The study of the planetary boundary layer (PBL) through the computational fluid dynamics, particularly the large eddie simulation (LES) is an important tool and driven by constant improvement of computing resources. The use of numerical simulations has some advantages over other techniques because it gives a full description of the vertical extent of the PBL. However, it also shows its limitations by simplifications in the equations solved and approximations of numerical methods. The technique employed in the construction of LES models utilizes the volumeavarege to solve the equations of motion for a particle of fluid. With this procedure is expects to solve these equations for large eddies, in which it is believed to be the most turbulent kinetic energy. We solved this scale. Therefore, we solve the scales associated to the large eddies. The rest of the turbulent scales, that is, the smaller scales are approximated by a subfilter or subgrid model. The separation between the scales resolved and subfilter is made by a process of filtering. The width of this filter is directly related to resolution of the grid computing. In this work we establish a comparison between the results obtained from the Moeng s LES model (MOENG, 1984) with experimental data obtained in the PBL. In addition, are accomplished comparisons with results from other models. Futhermore, were generated turbulent statisticals, which describe the flow in approch based on the turbulent statistical theory. To perform this task, were used for initial conditions and forcing the surface to a good weather day at the experimental site of Candiota. The performance of the LES model to reproduce the experimental data was considered satisfactory. However, the model was less effective in simulating the turbulent night. For this period is necessary a higher resolution grid. A CLP é a região da troposfera caracterizada pela existência da turbulência (BUSINGER, 1981). Um escoamento turbulento é aquele no qual as flutuações das variáveis que descrevem o escoamento são amplificadas no tempo e no espaço (LANDAU; LIFSHITZ, 1959). O estudo da camada limite planetária (CLP) através da dinâmica de fluidos computacional, sobretudo a simulação dos grandes turbilhões (LES), é uma ferramenta bastante empregada e impulsionada pelo constante aprimoramento dos recursos computacionais. O uso de simulações numéricas apresenta algumas vantagens em relação a outras técnicas pois fornece uma descrição completa da extensão vertical da CLP. No entanto, também apresenta suas limitações através de simplificações nas equações resolvidas e aproximações dos métodos numéricos. A técnica empregada na construção do modelo LES consiste na utilização de médias de volume para resolver as equações de movimento para uma partícula de fluido. Com isto espera-se resolver estas equações para os grandes turbilhões, nos quais acredita-se estar a maior parte da energia cinética turbulenta. Chamamos esta escala de resolvida. O restante das escalas, ou seja, as menores escalas do escoamento são aproximadas por um modelo de subfiltro ou subgrade. A separação entre as escalas resolvidas e de subfiltro é realizada por um processo de filtragem. A largura deste filtro está diretamente relacionada a resolução da grade computacional. Neste trabalho propõem-se uma comparação entre os resultados obtidos a partir do modelo LES de Moeng (MOENG, 1984) com dados experimentais obtidos na CLP. Além disto, são realizadas comparações com resultados de outros modelos. Adicionalmente, foram geradas estatísticas turbulentas, que descrevem o escoamento em uma aproximação baseada na teoria estatística da turbulência. Para isso foram utilizadas condições iniciais e forçantes de superfície para um dia de bom tempo do sítio experimental de Candiota. O desempenho do modelo LES em reproduzir os dados experimentais foi considerado satisfatório. Entretanto o modelo se mostrou menos eficaz em simular a turbulência noturna. Para tal período se faz necessária uma maior resolução de grade.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior Throughout this paper a spectral analysis is conducted on the superficial limit layer in complex terrains. The data used were collected in the central part of the state of Rio Grande do Sul, in the Valley of Jacuí River. A 15-meter tower, with fast-response sensors and slowresponse ones, collected data at frequencies of 10 Hz and 1 Hz, in July and August 2000, respectively. The time series were numeric analyzed through a piece of program developed in Fortran. The calculated spectra were classified according to stability class, intensity and direction of the wind speed. The spectra of the wind speed vertical component have a well-defined peak for all the analyzed conditions, except for the night series in which the wind direction is transverse to the valley axis. This same spectrum is in accordance with -5/3 Kolmogorov s law, with the beginning of the inertial sub-range in f ≈ 2 for winds that are parallel to the valley axis and any wind intensity. For winds transverse to the valley, the beginning of the inertial sub-range is in f ≈ 3. The frequencies associated with the spectral maxima are inferior to those observed in the Kansas experiment. Being superior in stable conditions when compared to convective ones in parallel conditions. In transverse conditions there is higher frequencies scattering for stable conditions. For parallel winds in stable conditions, the maximum no dimensional spectra are approximately 0.4 independently of z / L , and for convective conditions these maximum vary from 0.4 to 0.6. For transversal mean winds in stable conditions the high frequency vary from 0.4 to 0.5 and for convective conditions they are approximately 0.7. The spectra of the lateral components of velocity higher than 1m/s under stable conditions showed a cut frequency of ≈ 0.06 as initial number of this frequency, for all the cases in the low-frequency region. In the high-frequency region, the aliasing for all the parallel conditions starts at ≈ 5.0 and for transversal cases the starting point was at ≈ 10.0 . The consequences of the mesoescala movements are more important in the nocturnal boundary layer. The average time of ≈ 30 minutes renders contaminations in the computed flux measurements, since it captures mesoescala movements. The lateral spectra for convective conditions show an only spectral peak, highlighting that the importance of the thermal and mechanical effects have the same magnitude. In the region of low frequencies there is a scattering of data that could be explained by the factors associated with topographical influences. In the region of high frequencies the spectrum is in accordance with Kolmogorov s law, indicating, still under not homogeneous conditions, the presence of isotropic eddies. In convective and stable conditions with winds slower than 1m/s, the Kolmogorov s law is not applicable to most of the series; therefore such conditions were not analyzed in this study. For the different classes of z / L it is shown that the reason between vertical and horizontal spectra w u S / S fast increases until it reaches its isotropic number. For the parallel condition, we have a frequency of about f ≅ 2 and for the transverse case, about f ≅ 3. Neste trabalho faz-se uma análise espectral da turbulência na camada limite superficial em terrenos complexos. Os dados usados neste trabalho foram coletados no centro do estado do RS, no Vale do Rio Jacuí. Uma torre de 15 metros com sensores de resposta rápida e resposta lenta coletou dados com freqüências de 10 Hz e 1 Hz, respectivamente em julho e agosto de 2000. As séries temporais foram submetidas a um tratamento numérico a partir de um programa desenvolvido em Fortran. Os espectros calculados foram então classificados, conforme a classe de estabilidade, intensidade e direção da velocidade do vento. Os espectros da componente vertical da velocidade do vento, possuem um pico bem definido para todas as condições analisadas, excetuando-se as séries noturnas nas quais a direção do vento é transversal ao eixo do vale. Este mesmo espectro obedece a lei de -5/3 de Kolmogorov, com o início do sub intervalo inercial em f ≈ 2 para ventos paralelos ao eixo do vale e qualquer intensidade do vento. Para ventos transversais ao vale o início do subintervalo inercial ocorre em f ≈ 3. As freqüências associadas com os máximos espectrais são inferiores aquelas observadas no experimento e Kansas. E são maiores em condições estáveis do que em condições convectivas para condições paralelas. Para condições transversais existe um maior espalhamento nas freqüências para as condições estáveis. Para ventos paralelos em condições estáveis, os máximos espectrais adimensionais, são aproximadamente iguais a 0,4 independente de z / L , e para condições convectivas, estes máximos variam de 0,4 a 0,6. Para ventos médios transversais em condições estáveis, os máximos variam entre 0,4 e 0,5 e para condições convectivas, são de aproximadamente 0,7. Os espectros das componentes laterais da velocidade com ventos maiores que 1 m/s sobre condições estáveis mostraram uma freqüência de corte de ≈ 0,06 como valor de início desta freqüência para todos os casos na região de baixas freqüências. Na região de altas freqüências, o aliasing , para o caso de condições paralelas tem início em ≈ 5,0 e para os casos transversais o início foi em ≈ 10,0 . A conseqüência dos movimentos de mesoescala são mais importantes na camada limite noturna. O tempo de média de ≈ 30 minutos leva à contaminações do fluxo computado por capturar movimentos de mesoescala. Os espectros laterais para condições convectivas mostraram um único pico espectral, evidenciando-se que a importância dos efeitos térmicos e mecânicos são de mesma magnitude. Na região de baixas freqüências existe espalhamento de pontos que são explicáveis por fatores associados às influências topográficas. Na região de altas freqüências o espectro segue a lei de Kolmogorov indicando, ainda sobre condições de inomogeneidade, a presença de vórtices isotrópicos. Para condições convectivas e estáveis com ventos médios menores que 1 m/s, a lei de Kolmogorov não é verificada para a maioria das séries, e por esta razão não foi analisada neste estudo. Para as diferentes classes de z / L mostra-se que a razão entre os espectros vertical e horizontal w u S / S cresce rapidamente até atingir seu valor isotrópico. Para a condição paralela, temos uma freqüência em torno de f ≅ 2 e para o caso transversal, em torno de f ≅ 3.

Conselho Nacional de Desenvolvimento Científico e Tecnológico Serial mesoscale convective system (MCS) events with lifetimes over 18 h and up to nearly 70 h are routinely observed over southeastern South America from infrared satellite imagery during the spring and summer. These events begin over the southern La Plata River basin, with individual convective systems generally moving eastward with the cloud-layer mean wind. However, an important and common subset of these serial MCS events shows individual MCSs moving to the east or southeast, yet the region of convective development as a whole shifts upstream to the north or northwest. Analyses of the composite mean environments from 10 of these upstream-propagating serial MCS events using NCEP/NCAR reanalysis data events indicates that the synoptic conditions resemble those found in mesoscale convective complex environments over the United States. The serial MCS events form within an environment of strong low-level warm advection and strong moisture advection between the surface and 700 hPa from the Amazonas region southward. One feature that appears to particularly influence the lowlevel flow pattern at early times is a strong surface anticyclone located just off the coast of Brazil. At upper-levels, the MCSs develop on the anticyclonic side of the entrance region to an upperlevel jet. Mean soundings show that the atmosphere is moist from the surface to near 500 hPa, with values of convective available potential energy above 1200 J kg-1 at the time of system initiation. System dissipation and continued upstream propagation to the north and northwest occurs in tandem with a surface high pressure system that crosses the Andes Mountains from the west. The 72 h simulated idealized numerical experiment reproduces the main conditions synoptic patterns in a upstream propagating mesoscale convective system, and generates a unique 18 h Upstream propagating MCS, which produces a significant precipitation event with more than 80mm/18h in some regions, it represents about 50% of the total monthly summer precipitation.A transient anticyclone drives the large scale low level circulaton at the early times, it produces moisture an war advection from Amazonia southward. A low thermal pressure system close to 22oS is observed and its interaction between with the south Atlantic anticyclone circulaton produces a moisture corridor transport promoted into a northerly-enhanced Low Level Jet. The numerical experiment evidences the importance of local storms to the development and triggering in the mesoscale convective environments, since it forms the first convective cluster what becomes organized and grow in size. After the MCS becomes mesoscale organized, its outflows dominates the convective environment interacting with the gravity waves pre-existing region, but the gravity waves propagation upstream still exerting a important hole since it generate new convective areas ahead of the outflow boundary. During the MCS life two main characteristics can be observed: Gravity waves propagate out ahead of the line and may initiate new convective development upstream of the existing convection. Once this new convection develops downdrafts, the density-current mechanism again dominates its evolution. Sistemas convectivos de mesoescala (SCM) com tempos de vida superiores a 18 h e alguns com ate 70 h são rotineiramente observados sobre a região sul da América do Sul através de imagens de satélite principalmente durante o período de primavera e verão. Estes eventos iniciam sobre o sul da Bacia do Rio da Prata, com sistemas convectivos geralmente movendo-se para leste ou sudoeste, mas a região convectiva principal propaga-se contra o escoamento no sentido norte ou noroeste. A análise de compostos médios de dados de reanálise para 10 destes SCM com propagação serial contra o escoamento indica que as condições sinóticas assemelham-se àquelas encontradas nos ambientes de desenvolvimento de complexos convectivos e mesoescala sobre os Estados Unidos.SCMs com propagação serial ocorrem em ambientes de forte advecção quente em baixos níveis e forte advecção úmida na camada entre a superfície e 700 hPa. Esta advecção transporta calor e umidade da região Amazônica em direção ao sul da América do Sul. Um padrão que particularmente influencia o escoamento em baixos níveis é um forte anticiclone localizado na costa do Brasil que atua principalmente nos momentos que antecedem o SCM. Em níveis superiores, o SCM se desenvolve no lado anticiclônico da região de entrado do jato de altos níveis. Sondagens médias mostram que a atmosfera é umedecida desde a superfície até aproximadamente 500 hPa, com valores de CAPE > 1200 J kg-1 no momento de iniciação do sistema. A dissipação do sistema e a continuada propagação para norte e noroeste ocorre em fase com um sistema de alta pressão que cruza os Andes vindo de oeste. O experimento de simulação numérica de 72h reproduz os principais padrões sinóticos observados nos dados de reanálise e gera um SCM com propagação contra o escoamento de 18 h de vida, o qual produz um significante volume de precipitação com que chegam a mais de 80mm/18h em algumas regiões, isto representa aproximadamente 50% do total de um mês de verão. Um anticiclone transiente governa o escoamento de larga escala nos momentos iniciais da simulação, produzindo advecção quente e úmida do Equador em direção ao Pólo Sul. Uma baixa térmica próxima de 22oS interage com a Alta do Atlântico Sul e produz um corredor de umidade associado à um jato de baixos níveis. O experimento numérico evidencia a importância de tempestades locais no desenvolvimento e iniciação de sistemas convectivos de mesoescala, uma vez que elas tornam-se organizadas em aglomerados convectivos de mesoescala que dão origem aos SCM. A atividade de ondas de gravidade nos momentos inicial do sistema dá suporte para a iniciação convectiva, após o SCM tornar-se organizado e produzir os primeiros outflows , correntes de densidade dominam o ambiente convectivo interagindo com as ondas de gravidade pré-existentes e com o jato em baixos níveis. Durante a vida do SCM duas características principais de sua evolução podem ser observadas: Num primeiro momento, ondas de gravidade propagam-se em bandas gerando e organizando a convecção. Uma vez que a convecção torna-se organizada e produz os primeiros downdrafts , correntes de densidade tornam-se o principal mecanismo no desenvolvimento convectivo.

In this thesis, I discuss the results of radio polarization observations of selected Compact Steep Spectrum sources made with the Very Long Baseline Array using selected frequencies in the 18 cm, 21 cm, 13 cm and 6 cm observing bands. An analysis of the polarized emission for 3C216, C147, 3C119 and 3C43 allowed us to conclude that they show a hook-shaped structure suggesting that it may be colliding with a dense environment since both jet brightness and fractional polarization increases at the hook. The average rotation measures in the hook are about 30% higher than the values seen in superluminal knots nearest to the core component with rotation measures exceeding 1500 ± 200 rad m2. The polarized flux suggests that although the average Faraday depth < neBd > increases towards the core within the first 80 pc of the jet, it becomes higher as the jet bends. The source components which are nearest to the cores are generally unpolarized even at such short wavelength, suggesting that the spread in the λ2 polarization spectra is much narrower than in extended sources and that the wavelengths at which the fractional polarization falls by a half of the value at zero depolarization may be in the submillimeter wavelengts. The properties of the Faraday depolarization and rotatation seen in the selected objects sources are briefly discussed, as well as hypothesis to interpret them. Nesta tese, discuto resultados de observações em polarização de radiofontes compactas selecionadas do grupo CSS, feitas com o Very Long Baseline Array em freqüências selecionadas nas bandas de observação em 18 cm, 21 cm, 13 cm e 6 cm. Uma análise da emissão polarizada para os quasares 3C216, 3C147, 3C119 e 3C43 permite concluir que todos mostram um perfil em forma de gancho, o que sugere que podem estar colidindo com um meio ambiente denso, o que explicaria a distribuição de brilho do jato e o aumento da polarização fracional na direção da curva do gancho. As medidas de rotação nas curvas são aproximadamente 30% mais altas que os valores vistos nas regiões dos nós superluminais mais próximas às componentes nucleres, onde a medida de rotação excede 1500 ± 200 rad m2. Uma análise da despolarização < neBd > que aumenta na direção do núcleo dentro dos primeiros 80 pc dos jatos torna-se mais alta à medida que o jatos tornan-se curvos. As componentes das radiofontes que são mais distantes aos núcleos são geralmente despolarizadas, em relação às componentes mais próximas, mesmo em curtos comprimentos de onda, sugerindo que a dispersão no espectro de polarização com λ2 é bem mais estreita que em radiofontes estendidas e que os comprimentos de onda nos quais a polarização fracional cai pela metade do valor esperado em despolarização nula podem estar localizados nas bandas submilimétricas. Essas componentes, portanto, sofrem despolarização com a distância. As propriedades da despolarização e da rotação de Faraday, vistas nessa seleção de radiofontes, são discutidas brevemente, assim como as hipótese para interpretá-las.

Conselho Nacional de Desenvolvimento Científico e Tecnológico According to the Kolmogorov's fully developed turbulence theory established in 1941, there is an equilibrium range in which phenomenological quantities associated to the turbulence flow are described by laws of universal character. The main objection to this theory was given by Landau, who critized the universality of these laws arguing that an universal statistical equilibrium state would not be reached since there is a dependence on energy transference between the eddies belonging to the inertial range with the largest ones, which are external to this interval, thus non-universal. From the discussion between Kolmogorov and Landau, emerged in 1962 a theory proposed to be a correction to the 1941's theory. Following the new ideas, an intermittent energy cascade should exist in the inertial range thus, as a consequence, corrections to the scaling laws prediced by 1941's theory should to be taken in account. Since that, intermittency has became the central problem in fully developed turbulence. Another question which still remains from the discussion between Kolmogorov and Landau is that concerning to the universality of the second order structure constant, known as the Kolmogorov constant. Some theoreticians claim that this constant should be Reynolds' number dependent, however experimental evidences contradict themselves. 14 In this work, based on structure function of several orders of velocity, temperature and a passive scalar, anisotropy is related to a strange behavior of velocity odd-order structure functions. The Kolmogorov constant is determined in diferent flows and appear to be universal with the Reynolds' number. In addiction, the scatter and the dependence with the Reynolds' number found in the literature are explained as a consequence of anisotropy in the intertial subrange and imprecise determination of the dissipation rate. On other hand, experimental verification of intemittency in inertial subrange through the correlation of the rate dissipation R²,²(r) is explained to be an undesired efect of anisotropy on the equations whose derivation is based on the isotropy assumption. The existence of intermittency correction coeficent μ of the second order structure funtion GR2 is attributed to be a consequence of an arbitrary choice of the inertial subrange, nite Reynolds' number efect and another singularities. These conclusions lead to the idea that the classical theory of 1941 is correct in its original form and the discrepancies observed are due the nonexistence of a real fully developed turbulence state. A teoria de turbulência completamenta desenvolvida de Kolmogorov, elaborada em 1941, prevê a existência de um intervalo de equilíbrio, onde leis de caráter universal descrevem as quantidades fenomenológicas associadas ao escoamento turbulento. Entretanto Landau criticou a universalidade destas leis alegando que um equilíbrio estatístico universal não seria alcançado, desde que há uma dependência na transferência de energia entre os turbilhões que pertencem ao intervalo de equilíbrio e os turbilhões externos a este intervalo. Desta discussão entre Kolmogorov e Landau, surgiu em 1962 uma teoria de correção à de 1941. Esta teoria prevê a existência de uma cascata de energia intermitente entre os turbilhões pertencentes ao intervalo inercial de equilíbrio. Desde então intermitência tornou-se o centro das pesquisas em turbulência bem desenvolvida. Outra questão que ainda é motivo de discussão é a universalidade da constante da função estrutura de segunda ordem prevista por Kolmogorov. Alguns teóricos argumentam que esta constante seria dependente do número de Reynolds do escoamento, entretanto as evidências experimentais são contraditórias. Na presente análise, baseado em funções estrutura de diversas ordens de velocidade, temperatura e um escalar passivo, a anisotropia é relacionada à um estranho comportamento das funções estrutura de velocidade de ordem ímpar. Além do 12 mais, o estudo aponta que a constante de Kolmogorov, determinada em diferentes escoamentos, apresenta um comportamento universal com o número de Reynolds. Adicionalmente, o espalhamento desta constante em torno das médias e a dependência com o número de Reynolds obtidos nos trabalhos encontrados na literatura são explicados como um efeito devido à anisotropia no sub-intervalo inercial e a determina ção imprecisa da taxa de dissipação. Por outro lado, a verificação experimental da intermitência no sub-intervalo inercial através da correlação da taxa de dissipação R²,²(r) é explicada ser um efeito indesejado de anisotropia nas equações baseadas na hipótese de isotropia. A existência de um valor significativo para o coeficiente de correção de intermitência μ da função estrutura de segunda ordem GR2 é atribuído à uma escolha arbitrária do sub-intervalo inercial e a efeitos de número de Reynolds bonito e outras singularidades. Desta forma conclui-se que a teoria clássica de 1941 está correta na sua forma original e as discrepâncias observadas originalmente são devido a não existência de um estado de turbulência completamente desenvolvida.

We report on three campaigns to measure the structure of the nocturnal boundary layer in cleared and forested areas near Santarém, PA, in 2001 and 2003. The aim was to describe the temporal evolution of wind and thermodynamic profiles in the stable surface layer and to use the accumulation of scalars, such as CO2 in the layer as a way to infer surface fluxes. This was deemed necessary because the frequent very calm conditions limited the validity of eddy covariance flux CO2 measurements at these sites at night. In 2001, tethered balloon observations were made at a pasture site. In 2003, this study was repeated, with a special balloon sonde deployed to measure CO2, temperature, humidity, and pressure. During the convective period just after sunrise, both the eddy covariance and accumulation flux estimation methods are valid and were compared. At night the key is to determine the thickness of the atmospheric layer to which nocturnal respiratory CO2 accumulates. In the forest at night the CO2 accumulation layer was basically constant and had a height of 32 m. The mean CO2 night flux (Fc) was estimated as 0.09(±0.02) mgCO2/m2s. The average CO2 accumulation rate was constant into de forest, decreased linearly with height form canopy height to 100-140 m high. These results are in agreement with the average profiles measured from forest flux towers. Over plowed field we found a larger CO2 accumulation layer at 60-80 m, reaching more then 100 m on some nights. The mean CO2 night flux was estimated Fc=0.20(±0.2) mgCO2/m2s. We found two shape of CO2 accumulation rate near the ground. Two nights the accumulation rate was basically constant with height, other four nights it decreased linearly with height, from the surface to 30-50 m high. The level what the accumulation rate is zero vary from 80 m to 135 m high. Três campanhas de medidas da estrutura da camada limite noturna em área limpa e de floresta próximo de Santarém, PA, realizadas em 2001 e 2003 são apresentadas neste trabalho. O objetivo principal foi descrever a evolução temporal dos perfis de vento e termodinâmicos na camada limite estável para usar o método de acúmulo de escalares, como o CO2, nesta camada, como meio de estimar os fluxos superficiais. Isso se torna necessário porque freqüentemente condições de ventos calmos sobre estes sítios ocorrem durante o período noturno limitando a confiança das estimativas do fluxo de CO2 realizadas pelo método de covariância de vórtices turbulentos. Em 2001, usando balão cativo as observações foram realizadas no sítio de pastagem. Em 2003, estas observações foram repetidas com uma sonda especialmente desenvolvida para medir a concentração de CO2, temperatura, umidade e pressão do ar. Durante os primeiros estágios da convecção, logo após o nascer do sol, ambos métodos de covariância de vórtices e acumulação são validos e foram comparados. Durante a noite o principal é determinar a profundidade da camada atmosférica dentro da qual o CO2 acumula. Na floresta, à noite a camada de acumulação de CO2 basicamente foi constante e a altura foi de 32 m. O fluxo médio de CO2 à noite (Fc) foi estimado igual a 0.09(±0.03) mgCO2/m2s. Na média a taxa de acumulação de CO2 se manteve constante dentro da floresta, decrescendo até aproximadamente zero até níveis que variaram desde 100 m até 140 m. Estes resultados estão de acordo com as observações das torres de fluxo na região de estudo dos dados de longo termo do perfil do gás. Sobre o sítio com o solo arado a camada de acumulação de CO2 à noite em média varia entre 60-80 m, atingindo mais de 100 m em algumas noites. O fluxo noturno médio de CO2 foi estimado em Fc=0.20(±0.2) mgCO2/m2s. A variação vertical da taxa de acumulação apresentou dois padrões distintos próximo a superfície. Duas noites a variação vertical da taxa de acumulação foi aproximadamente constante, em outras quatro noites decresceu rapidamente com a altura, desde a superfície até 30-50 metros de altura. O nível a partir do qual a taxa de acumulação é nula variou desde 80 m a 135 m.

In the study of the dispersion of contaminants in the atmosphere by computational simulations, the Eulerian approach, based in the diffusionadvection equation, is largely employed. In these simulations, concentration values are calculated at each point of a fixed grid in space and time and for this to be done it is necessary to determine an appropriated eddy diffusivity to the planetary boundary layer under consideration. The aim of this work is to describe the decay of turbulence in the residual layer and to derive the corresponding eddy diffusivity for the case of a homogeneous (in the horizontal plane) and isotropic turbulent flow neglecting mechanical and thermal production effects. The decay of turbulence is described as a process in which eddies of a certain size loss their energy by the effect of an eddy viscosity associated with the smallest eddies through a cascade of interactions of local character in the sense of Heisenberg s model. The eddy diffusivities are calculated by considering the dispersive effects of the energy-containing eddies on the basis of Taylor s statistical diffusion model. Considering the hypothesis that the law of decay of the energy spectral density function associated with the ith axis of the reference system is similar to the law of decay of the tridimensional energy spectral density function, the corresponding variance of the turbulent part of the velocity and the corresponding eddy diffusivity are calculated. These quantities, in particular for the vertical axis, are compared with results of large eddy simulations (LES) showing very good agreement. Finally, incorporating the calculated eddy diffusivity in a Eulerian model, the effect of turbulence in the vertical dispersion of an passive scalar contaminant in the residual layer is described for the special case of a instantaneous bidimensional emission source (area source) at a certain height over land and zero mean wind velocity. No estudo da difusão de um contaminante na atmosfera por meio de simulações computacionais é bastante utilizada a perspectiva euleriana, baseada na equação da difusão-advecção. Nessas simulações, os valores de concentração são calculados em cada ponto de um retículo fixo no espaço e no tempo e para que isso possa ser feito, é necessário determinar a difusividade turbulenta (coeficiente de difusão turbulento) para o tipo de camada limite planetária em questão. O objetivo desse trabalho é descrever o decaimento da turbulência na camada residual e derivar a correspondente difusividade turbulenta no caso de um escoamento turbulento homogêneo (no plano horizontal) e isotrópico, desprezando os efeitos mecânicos de cisalhamento e térmicos das correntes de convecção. O decaimento da turbulência é descrito como um processo em que os turbilhões de certo tamanho perdem energia por efeito de uma viscosidade cinemática turbulenta associada aos turbilhões menores através de interações em cascata de caráter local nos moldes do modelo de Heisenberg. A difusividade turbulenta é calculada pelos efeitos dispersivos dos turbilhões energéticos com base no modelo estatístico de Taylor. Com a hipótese de que o decaimento da função densidade espectral de energia unidimensional associada ao i-ésimo eixo do sistema de referência é similar à lei de decaimento da função densidade espectral tridimensional, calcula-se a correspondente variância da componente turbulenta da velocidade e a correspondente difusividade turbulenta e essas grandezas associadas ao eixo vertical, em particular, são comparadas com resultados obtidos por simulação numérica de grandes turbilhões (LES) mostrando uma concordância muito boa. Finalmente, incorporando a difusividade turbulenta calculada num modelo euleriano, descreve-se o efeito da turbulência na dispersão vertical de um contaminante escalar passivo na camada residual no caso especial de uma fonte bidimensional (fonte área) com emissão instantânea situada a uma certa altura acima do solo e na ausência de vento.