Neste trabalho foram utilizadas imagens do satélite meteorológico GOES-8, em órbita geoestacionária, tomadas em três bandas do infravermelho, para caracterizar a absorção de radiação nestes comprimentos de onda, em função da variação espacial de parâmetros atmosféricos. Um dos principais objetivos foi o desenvolvimento de métodos capazes de apontar, em regiões préselecionadas por sua altitude, zonas de menor absorção atmosférica no infravermelho, caracterizando-as como sítios potencialmente apropriados para observações astronômicas nestes comprimentos de onda. A área de estudo faz parte dos Andes Peruanos, contendo regiões com mais de 2500 m de altitude. A metodologia do trabalho baseou-se na escolha dos canais 3, 4 e 5 do satélite GOES devido a sua sensibilidade à presença de vapor d'água na atmosfera. Em especial, no canal 3 é possível detectar a emissão de radiação infravermelha pelo vapor d'água presente nos níveis médio e superior da troposfera, enquanto no canal 4 a atmosfera quase não absorve radiação infravermelha, concentrando-se este efeito nas nuvens. Considerou-se que na ausência de nuvens pode-se observar a radiação proveniente do solo permitindo estimar sua temperatura. O canal 5 situa-se na porção do espectro contaminada pelo vapor d'água, isto é, onde o vapor d'água absorve parte da radiação proveniente da baixa troposfera. Como critério de qualidade dos sítios estudados foram escolhidos a altitude, a visibilidade da superfície e a quantidade de vapor d'água presente na atmosfera. Os resultados mostraram que as regiões ao sul do Peru concentram os sítios mais adequados à instalação de um observatório astronômico. Comparações com resultados de outros pesquisadores reforçam nossas conclusões.

Neste trabalho foram utilizadas imagens do satélite meteorológico GOES-8, em órbita geoestacionária, tomadas em três bandas do infravermelho, para caracterizar a absorção de radiação nestes comprimentos de onda, em função da variação espacial de parâmetros atmosféricos. Um dos principais objetivos foi o desenvolvimento de métodos capazes de apontar, em regiões préselecionadas por sua altitude, zonas de menor absorção atmosférica no infravermelho, caracterizando-as como sítios potencialmente apropriados para observações astronômicas nestes comprimentos de onda. A área de estudo faz parte dos Andes Peruanos, contendo regiões com mais de 2500 m de altitude. A metodologia do trabalho baseou-se na escolha dos canais 3, 4 e 5 do satélite GOES devido a sua sensibilidade à presença de vapor d'água na atmosfera. Em especial, no canal 3 é possível detectar a emissão de radiação infravermelha pelo vapor d'água presente nos níveis médio e superior da troposfera, enquanto no canal 4 a atmosfera quase não absorve radiação infravermelha, concentrando-se este efeito nas nuvens. Considerou-se que na ausência de nuvens pode-se observar a radiação proveniente do solo permitindo estimar sua temperatura. O canal 5 situa-se na porção do espectro contaminada pelo vapor d'água, isto é, onde o vapor d'água absorve parte da radiação proveniente da baixa troposfera. Como critério de qualidade dos sítios estudados foram escolhidos a altitude, a visibilidade da superfície e a quantidade de vapor d'água presente na atmosfera. Os resultados mostraram que as regiões ao sul do Peru concentram os sítios mais adequados à instalação de um observatório astronômico. Comparações com resultados de outros pesquisadores reforçam nossas conclusões.

Neste trabalho foram utilizadas imagens do satélite meteorológico GOES-8, em órbita geoestacionária, tomadas em três bandas do infravermelho, para caracterizar a absorção de radiação nestes comprimentos de onda, em função da variação espacial de parâmetros atmosféricos. Um dos principais objetivos foi o desenvolvimento de métodos capazes de apontar, em regiões préselecionadas por sua altitude, zonas de menor absorção atmosférica no infravermelho, caracterizando-as como sítios potencialmente apropriados para observações astronômicas nestes comprimentos de onda. A área de estudo faz parte dos Andes Peruanos, contendo regiões com mais de 2500 m de altitude. A metodologia do trabalho baseou-se na escolha dos canais 3, 4 e 5 do satélite GOES devido a sua sensibilidade à presença de vapor d'água na atmosfera. Em especial, no canal 3 é possível detectar a emissão de radiação infravermelha pelo vapor d'água presente nos níveis médio e superior da troposfera, enquanto no canal 4 a atmosfera quase não absorve radiação infravermelha, concentrando-se este efeito nas nuvens. Considerou-se que na ausência de nuvens pode-se observar a radiação proveniente do solo permitindo estimar sua temperatura. O canal 5 situa-se na porção do espectro contaminada pelo vapor d'água, isto é, onde o vapor d'água absorve parte da radiação proveniente da baixa troposfera. Como critério de qualidade dos sítios estudados foram escolhidos a altitude, a visibilidade da superfície e a quantidade de vapor d'água presente na atmosfera. Os resultados mostraram que as regiões ao sul do Peru concentram os sítios mais adequados à instalação de um observatório astronômico. Comparações com resultados de outros pesquisadores reforçam nossas conclusões.

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.

Orientador : Joaquim Severino Paiva Netto / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia de Campinas / Made available in DSpace on 2018-07-15T08:43:43Z (GMT). No. of bitstreams: 1 Park_KilJin_D.pdf: 2331857 bytes, checksum: c1cd37e1ec0d4ca6c47d47aec01099f8 (MD5) Previous issue date: 1979 / Resumo: Neste trabalho foram estudadas as correlações de transferências de calor e de quantidade de movimento em um trocador, munido de promotores estacionários de turbulência......Observação: O resumo, na íntegra, poderá ser visualizado no texto completo da tese digital / Abstract: Heat and momentum transfer correlations were studied in exchanger fitted with mationless terbulence promoters...Note: The complete abstract is available with the full electronic digital thesis or dissertations / Doutorado / Doutor em Engenharia Mecânica

Para vegetação esparsa e de porte alto a determinação dos parâmetros de rugosidade é comprometida pela dificuldade em se observar condições que satisfaçam a lei logarítmica da velocidade do vento. Estimou-se o comprimento de rugosidade (z0) e o deslocamento do plano zero (d) por alguns métodos com medidas micrometeorológicas e da estrutura física de arbustos esparsos em região semi-árida, durante o experimento HAPEX-Sahel. A velocidade do vento foi medida em quatro alturas acima da superfície (3,0; 4,1; 5,3 e 8,5 m), e os fluxos determinados por correlações dos turbilhões a 9m de altura. Métodos baseados no perfil logarítmico foram aplicados em condições de atmosfera neutra. A altura média da vegetação era h = 2,06 ± 0,47 m. O método convencional (ajuste estatístico) resultou em estimativas satisfatórias de d e z0 em condições nas quais a validade do perfil logarítmico foi satisfeita. Com uma única altura de medida localizada acima da subcamada inercial as estimativas resultaram em valores ou fisicamente inconsistentes ou que não caracterizam a rugosidade da superfície. Quando se utilizou a velocidade de fricção dada pela correlação dos turbilhões na solução do perfil logarítmico, as estimativas melhoraram. A combinação do perfil logarítmico com a relação z0 = λ (h - d) proporcionou estimativas satisfatórias para os valores de λ = 0,188 e 0,190 determinados em função da estrutura física da vegetação, o que não foi observado para o valor médio da literatura (0,166). Relações entre a estrutura física da vegetação e o transporte de momentum estimaram apropriadamente d e z0. A rugosidade da área foi melhor descrita por d = 0,95 m = 0,46 h e z0 = 0,204 m = 0,1 h, sendo λ = 0,185. As velocidades horizontal do vento e de fricção foram mais sensíveis a variações em z0 do que em d. / For sparse and tall vegetation the estimate of roughness parameters is compromised by the difficulty in observing conditions that satisfy the windspeed logarithmic law. The roughness length (z0) and the zero-plane displacement (d) were estimated by some methods with micrometeorological measurements and the physical structure of sparse shrubs in semi-arid region, during the HAPEX-Sahel experiment. The wind speed was measured at four heights above of surface (3.0, 4.1, 5.3 and 8.5 m), and the turbulent flows determined by eddy correlations at the height of 9m. Methods based on the logarithmic profile have been applied in neutral atmosphere conditions. The average height of the vegetation was h = 2.06 ± 0.47 m. The conventional method (statistical fit) resulted in good estimates of d and z0 only under conditions of validity of the logarithmic law. Only one height of measurement located above of the inertial sublayer is enough to result in physically inconsistent values. When the friction velocity, given by eddy correlation, was used in the logarithmic law, the estimates improved. The combination of the logarithmic law with z0 = λ (h - d) provided satisfactory estimates of the surface roughness for λ = 0.188 and 0.190 determined in function of the physical structure of the vegetation; but for λ = 0.166, the average value of literature, the estimates where not good. Relationships between the physical structure of the vegetation and the momentum transfer estimated appropriately d and z0. The area roughness was better described by d = 0.95 m = 0.46 h and z0 = 0.204 m = 0.1 h, being λ = 0.185. Wind speed and friction velocity were more sensible to variations in z0 than in d.

O incêndio florestal é uma complexa combinação da energia liberada na forma de calor devido à combustão dos produtos oriundos da pirólise da vegetação e o transporte dessa energia para o ar e para a vegetação à sua volta. O primeiro é o domínio da química e ocorre na escala de moléculas e o segundo é o domínio da física e ocorre em escalas de até quilômetros. É a interação desses processos sobre uma ampla gama de escalas temporais e espaciais envolvidas no incêndio florestal que faz a modelagem do seu comportamento uma tarefa tão difícil. A propagação do incêndio através de vegetação rasteira e folhas mortas foi simulada numericamente usando a formulação física do WFDS. A abordagem utilizada foi tridimensional e transiente, e baseada em uma descrição dos fenômenos físicos que contribuem para a propagação de um incêndio de superfície através de uma camada de combustível. Neste cenário de incêndio, existem duas regiões: vegetação e ar, cada uma com suas propriedades físicas e químicas e, embora elas precisem ser integradas no mecanismo de solução, há diferentes fenômenos que ocorrem em cada uma. Na região de vegetação, a abordagem é representá-la como partículas submalha cercadas de ar. O caráter heterogêneo da vegetação, como sua natureza, folhagens, pequenos galhos, etc. foi levado em conta usando propriedades físicas médias características da floresta amazônica. Os fenômenos na região de vegetação são a evaporação da sua umidade, a pirólise e a transferência de calor por radiação e por convecção. Na região do ar, a combustão com chama ocorre em um ambiente turbulento, onde as transferências de calor por radiação e por convecção desempenham um papel significativo. Para incorporar a radiação dos gases de combustão, o modelo físico emprega o método de volumes finitos, que resolve a equação de transferência de calor por radiação como uma equação de transporte para um número finito de discretos ângulos sólidos, e que pode ser usado em uma ampla faixa de espessuras óticas e meios participantes. A combustão turbulenta para a fase gasosa é modelada com base no modelo Eddy Dissipation Concept (EDC). O modelo de combustão turbulenta adota a hipótese de reação química infinitamente rápida entre o combustível e o ar e é controlado apenas pela velocidade de mistura desses reagentes. Esse modelo representa bem a física de incêndios em ambientes ventilados, como é o caso dos incêndios florestais. Para incluir os efeitos do transporte turbulento é utilizado o método Large Eddy Simulation (LES), que calcula explicitamente as grandes estruturas turbulentas, mas trata a dissipação e a cascata inercial em escalas menores usando aproximações na escala submalha. As regiões de vegetação e ar trocam massa e energia. O comportamento da mistura gasosa resultante da degradação térmica da vegetação e das reações de combustão é regido pelas equações de Navier-Stokes. As equações que regem os modelos físicos são formuladas como equações diferenciais parciais que são resolvidas por métodos numéricos. O método utilizado para discretização das equações é o método de diferenças finitas em malha deslocada. O modelo numérico utilizado resolve as equações de Navier-Stokes para fluidos compressíveis usando o filtro de Favre. A dissipação de energia cinética é obtida através de um fechamento simples para a tensão turbulenta: o modelo de coeficiente constante de Deardorff. O transporte turbulento de energia e massa é contabilizado pelo uso, respectivamente, de números de Prandtl e de Schmidt turbulentos constantes. Os resultados das simulações do modelo físico descrito foram comparados aos dados experimentais obtidos em campo para a propagação do incêndio na floresta amazônica. Apesar da idealização das condições de combustível, vento e as incertezas dos dados experimentais, as previsões do modelo estão na mesma ordem de grandeza dos experimentos. As taxas de propagação do incêndio experimentais variam de 0,12 +/-0,06 a 0,35+/-0,07 m/min. Mesmo considerando-se o desvio padrão da taxa de propagação do incêndio experimental, os valores das taxas simuladas ficaram dentro do erro experimental somente em dois de sete casos. As simulações mostraram que os parâmetros importantes para o modelo são a área superficial por volume da vegetação, sua massa específica aparente e sua umidade. Como o coeficiente de absorção por radiação é função direta da massa específica aparente e da área superficial por volume da vegetação, esses parâmetros afetam o comportamento numérico do incêndio de superfície. De acordo com os resultados das simulações numéricas, a umidade da vegetação também tem importância no incêndio de superfície. A temperatura inicial da vegetação e a umidade do ar na faixa de variação analisada não influenciam a taxa de propagação do incêndio. As simulações também mostraram que o processo de radiação é muito importante, e afeta diretamente todos os demais processos e a taxa de propagação do incêndio. A convecção tem importância muito menor que a radiação na condição de ausência de vento externo. A coerência das taxas de propagação do incêndio experimental e numérica em função da massa específica aparente de material combustível e da umidade da vegetação foi investigada. O modelo numérico é coerente em todas as nove combinações de casos. Já o experimento é coerente em quatro combinações. Com base nas comparações entre cada dois casos experimentais e as respectivas simulações numéricas, nota-se que as taxas de propagação a partir das simulações numéricas foram mais coerentes que as experimentais. / Forest fire is a complex combination of energy released as heat due to the combustion of the products from the vegetation pyrolysis and the transport of this energy to the surrounding air and vegetation. The first is the domain of chemistry and occurs on the molecular scale, and the second is the domain of physics and occurs at scales up to kilometers. It is the interaction of these processes on a wide range of temporal and spatial scales involved in forest fires that makes modeling its behavior such a challenging task. The spread of fire through small plants and dead leaves was simulated numerically using WFDS physical formulation. The approach used was three-dimensional and transient, based on a description of the physical phenomena that contribute to the spread of a surface fire through a layer of fuel. In this fire scenario, there are two regions: vegetation and air, each one with its physical and chemical properties and, although they need to be integrated into the solution mechanism, there are different phenomena that occur in each one. In the vegetation region, the approach is to represent it as subgrid particles surrounded by air. The heterogeneity of the vegetation, such as its nature, leaves, twigs, etc. was taken into account by using average physical properties that are representative of the Amazon forest. The phenomena in the vegetation region are the evaporation of its moisture, pyrolysis, heat transfer by radiation and convection. In the air region, the flaming combustion occurs in a turbulent environment, and heat transfer by radiation and convection play a significant role. To incorporate the radiation from the combustion gases, the physical model employs the finite volumes method, solving the radiation transfer equation as a transport equation for a finite number of discrete solid angles, which can be used in a wide range of optical thicknesses and participating media. Turbulent combustion for the gaseous phase is modeled using the Eddy Dissipation Concept (EDC) model. The mixing controlled turbulent combustion model adopts the assumption of infinitely fast chemical reaction between the fuel and air. This model represents well the fire physics in ventilated areas, as is the case of forest fires. To include the turbulent flow effects, it is used the Large Eddy Simulation (LES) method, which explicitly calculates the large turbulent structures, but models the dissipation and inertial cascade using approximations in the sub-grid scale. The vegetation and air regions exchange mass and energy. The behavior of the gas mixture resulting from the vegetation thermal degradation and combustion reactions is governed by the Navier-Stokes equations. The equations governing the physical model are formulated as partial differential equations, which are solved by numerical methods. The method used for discretization of the equations is the finite difference method on a staggered grid. The numerical model solves the Navier-Stokes equations for compressible fluids using the Favre filter. Dissipation of kinetic energy is achieved through a simple closure for the turbulent stress: the constant coefficient Deardorff model. The turbulent transport of heat and mass is accounted for by use of constant turbulent Prandtl and Schmidt numbers, respectively. The physical model simulation results were compared to experimental data obtained in the field for the spread of fire in the Amazon forest. Despite of the idealized conditions of fuel, wind and the uncertainties of the experimental data, the model predictions and the experiments are in the same order of magnitude. Experimental rate of spread range from 0.12 +/- 0.06 to 0.35 +/- 0.07 m/min. Even considering rate of spread experimental standard deviation, simulated rate values were within experimental error only in two of seven cases. The simulations showed that the important parameters for the model are the vegetation surface area to volume ratio, its bulk density and moisture. As the radiation absorption coefficient is a direct function of vegetation bulk density and surface area to volume ratio, these parameters affect the numeric behavior of the surface fire. According to the numerical simulations results, vegetation moisture is also important in the surface fire scenario. Vegetation initial temperature and air humidity in the range analyzed does not influence the rate of spread. The simulations also showed that the radiation process is very important and directly affects all other processes and rate of spread. Convection heat transfer has much less significance than radiation heat transfer in the absence of external wind. The consistency of the experimental and numerical rate of spread, as a function of combustible material bulk density and vegetation moisture was investigated. The numerical model is consistent in all nine case combinations. The experiment is consistent in four cases. Based on comparisons between each two experiments and their numerical simulations, it is noted that the rate of spread variation from the numerical simulation is more consistent than the experimental one.

Há diversas metodologias empregadas para a determinação das cargas de rajadas usadas em dimensionamentos estruturais da aeronave, e que incluem métodos determinísticos e métodos estocásticos. As cargas de resposta à turbulência continua são, hoje, obtidas pelo método estocástico (Random Process Theory - RPT) em conformidade com o estabelecido nos requisitos aeronáuticos em vigor (FAR). Este trabalho visa implementar uma alternativa ao RPT na determinação das cargas de turbulência contínua, baseada na teoria dos filtros correlacionados (Matched Filter Theory - MFT). Historicamente, o MFT foi originalmente utilizado na obtenção mde sinais maximizados de radares. Pototzky (POTOTZKY, 1997) demostrou que o MFT é aplicável também para a solução de problemas dinâmicos e aeroelásticos genéricos, especificamente para cálculos de cargas correlacionadas de rajadas. Será demonstrado aqui que as cargas correlacionadas geradas pelo MFT são fortemente similares às geradas pelo RPT. É apresentado neste documento uma descrição detalhada da metodologia de cálculo de cargas baseada em RPT e em MFT, e os resultados numéricos obtidos.

The influence of physical properties on heat transfer between solid and fluid phases is investigated for laminar and turbulent flows in a channel filled with a moving porous material. Concurrent, counterflow and crossflow configurations are analyzed. To simulate flow and heat transfer between phases, a two-energy equation model using a thermal non-equilibrium condition is applied. Transport equations are discretized using the control volume method and the system of algebraic equations is relaxed via the SIMPLE algorithm. Validations are made for laminar model under concurrent and counterflow configurations. Effects of thermal and hydrodynamic properties on heat transfer for several conditions are analyzed and compared with analytical results in the literature. For concurrent laminar flow, simulations indicate that, when the speed of the solid approaches that of the fluid, the strong axial convection of the solid phase, as well as the reduction of the relative velocity, cause an increase in the axial length needed for thermal equilibrium between phases to occur. Longer thermal developing lengths are also found for higher permeability and porosity. Results for a counterflow moving bed indicate that motion of the solid material, contrary to the direction of the fluid, enhances heat transfer between phases. The same effect is observed for smaller Darcy number and porosity, as well as for higher solid-to-fluid thermal capacity and thermal conductivity ratios. In the case of crossflow, where there are two fluid inlets, more energy is convected into the system in both longitudinal and transversal directions .The fluid temperature reaches the highest values in the symmetry region of the channel. This occurs mainly for high velocity, high thermal capacity and high thermal conductivity ratios between fluid and solid phases. These behaviors were observed for laminar and turbulent flows, in both fully filled and half filled channels. The studies presented here might have applications to problems involving engineering equipment in which a moving porous bed is identified.

Este trabalho objetivou a ampliação dos conhecimentos sobre a dinâmica de chamas de difusão turbulenta. A abordagem escolhida foi o método numérico da simulação de grandes escalas turbulentas (LES) acoplado ao modelo de combustão do tipo folha de chama ou \"Flamesheet\". A modelagem dos efeitos de dissipação das pequenas escalas turbulentas, foi realizada pelo modelo de Smagorinsky (1963). Foi feita a verificação e validação do código implementado pelo método das soluções manufaturadas, descrito por Roache (2002), e pelas soluções de engenharia do tipo tampa móvel em cavidade laminar, escoamento de Poiseuille e jato laminar. Neste trabalho foi estudada a chama de difusão turbulenta do tipo D, padronizada pelo \"SANDIA National Laboratories\". Foi verificado que, quando estudado o jato isotérmico, este obedece ao critério de jato auto-sustentável, definido por Hussein et al. (1994). Já os resultados obtidos para o jato reativo, quando comparados aos dados experimentais obtidos por Masri et al. (1996) e Barlow e Frank (1998), apresentaram boa concordância para os campos de velocidades e razoável correspondência para os campos de fração de mistura e temperatura, indicando que maiores estudos são necessários. / The mean goal of this work was to better understand the dynamics of turbulent diffusion flames. The turbulent flow was resolved appling the methodology of Large Eddy Simulation, coupled with the Flamesheet model for reactive systems. The effects of the turbulent subgrid dissipation were accounted with the Smagorinsky (1963) subgrid model. The numerical code verification and validation were accomplished with the manufactured solutions method, described by Roache (2002), and with the method of engineering solutions of the type of moveable wall in laminar cavity flow, Poiseuille flow and laminar jet flow. In this work was studied the turbulent flame of the type D flow, standardized by the SANDIA National Laboratories. It was verified that, when the isotermal jet flow was studied, the criterium of a self-sustained jet, defined by Hussein et al. (1994), was achieved. The results for the reactive flow, when compared against experimental data obtained by Masri et al. (1996) e Barlow e Frank (1998), presented resonable agreement for the velocity field and satisfactory correspondence for mixture fraction and temperature fields, indicating that further studies are necessary.