Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Nocturnal turbulent flux data are analysed, initially for two experimental sites in Brazil, and later for a large variety of biomes. In such case, a wide dataset, originated from more than a hundred flux stations scattered over Europe and the Americas is used. Emphasis is given to important aspects associated to the occurrence of intermittent turbulence and its implications for CO2 nocturnal flux estimates. Intermittency factors (IF) for CO2 and sensible heat fluxes are determined for the set of stations. In most cases, there was agreement between the results for CO2 and sensible heat, specially for those stations that were less intermittent, on average. A larger number of stations showed low IF for CO2 than for sensible heat fluxes. CO2 flux intermittency was significantly larger over higher canopies than over savannas or crops, while the same difference is smaller for the sensible heat flux. An investigation on the CO2 flux dependence on the temporal variability of the turbulence intensity is carried on . The main purpose of such an anlysis is identifying whether ther is CO2 storage during low turbulence periods and if such accumulated CO2 affects the fluxes in subsequent periods. In the vast majority of the stations, the classical pattern of NEE decreasing as u* tends to zero is observed. The hypothesis of the present study is that such NEE decrease under calm situations does not necessarily imply on flux lost by the eddy covariance technique, but that it may have not been properly captured by the CO2 storage term, being later measured in subsequent turbulent situations. Therefore, the present study has as the main goal identifying a potential problem in the correcting procedure commonly applied for low turbulent periods. Specifically, it is important to analyse whether the periods whose data are being replaced are not suceeded by others for which there is flux excess, characterizing, therefore, an excessive correction. Results indicate that there is a large quantity os stations for which the fluxes are independent on turbulence intensity variability, with evidences that the observed reduction during calm periods is compensated when it gets turbulent. In such cases, it is possible that corrections commonly used for calm conditions are, indeed, excessive. . / Dados de fluxos turbulentos noturnos são analisados, primeiramente, para dois sítios experimentais no Brasil e uma posterior expansão dessas análises é realizada para uma grande variedade de biomas, utilizando para tanto, um conjunto amplo de dados originados de mais de uma centena de estações de fluxos, espalhadas nas Américas e na Europa. A ênfase é dada aos aspectos importantes associados à ocorrência de turbulência intermitente e suas implicações para as estimativas de fluxos turbulentos noturnos de CO2. Fatores de intermitência (FI) dos fluxos de CO2 e de calor sensível são definidos e quantificados para o conjunto de estações. Na maioria dos casos, houve concordância entre os comportamentos obtidos para o fluxo de CO2 e de calor sensível, especialmente nas estações que na média apresentaram condições pouco intermitentes. Um número maior de estações apresentaram valores reduzidos de FI para fluxo de CO2 do que para o de calor sensível. A intermitência dos fluxos de CO2 se mostrou significativamente maior sobre estações de maior altura do dossel (florestas) do que sobre estações de savana ou agrícolas, enquanto que a mesma diferença é muito mais reduzida para os fluxos de calor sensível. Uma investigação da dependência dos fluxos de CO2 em relação à variabilidade temporal de turbulência ao longo da noite foi realizada. O propósito principal da análise é identificar se há acúmulo de CO2 em períodos pouco turbulentos e se este CO2 acumulado afeta os fluxos em períodos subsequentes. Na imensa maioria das estações, se observou um padrão clássico de NEE diminuindo conforme u* tende a zero. É hipótese do presente trabalho que o decréscimo de NEE em situações calmas não necessariamente indica fluxo não contabilizado pelo método da covariância dos vórtices, mas que este pode não ter sido propriamente capturado pelo termo de armazenamento, sendo posteriormente capturado pelo sensor em situações turbulentas subsequentes. Desta forma, o presente estudo tem como precípuo intento, contribuir no sentido de identificar um potencial problema na utilização das técnicas de correção de dados utilizadas para os períodos com pequena intensidade turbulenta. Especificamente, é importante analisar se os períodos que estão tendo seus dados substituídos não são sucedidos por outros em que houve excesso de fluxo, caracterizando, desta forma, uma correção excessiva. Os resultados indicam que há uma grande quantidade de estações para as quais os fluxos são independentes da variabilidade da intensidade turbulenta, havendo indicações que a redução observada em períodos calmos seja compensada em intervalos turbulentos. Nesses casos, é possível que correções comumente usadas para a correção de fluxos em condições calmas sejam excessivas.

Dois dos principais objetivos no estudo da aerodinâmica de veículos comerciais são a redução no consumo de combustível e o aumento na eficiência da refrigeração do motor. Esses objetivos podem ser alcançados através do desenvolvimento de dispositivos que modifiquem o escoamento do ar ao redor do veículo e também através da alteração da forma das superfícies externas. A inclinação das superfícies da parte traseira de um ônibus, por exemplo, tem grande influência sobre a esteira turbulenta que se forma atrás do veículo. O uso de ferramentas computacionais permite uma redução de custo e maior flexibilidade na análise aerodinâmica de autoveículos. Ainda é preciso, no entanto, que o resultado dessas ferramentas computacionais seja verificado com o maior número possível de casos para que se possa escolher e ajustar o modelo matemático de forma adequada. O objetivo do presente trabalho é a verificação dos resultados computacionais e experimentais no desenvolvimento de metodologias que visem à redução no consumo de combustível e aumento na eficiência da refrigeração do motor. Foram comparados resultados experimentais e computacionais do escoamento sobre um modelo de um ônibus comercial em escala 1:17,5. Para a realização do experimento foi utilizado um túnel de vento de seção aberta, onde foram analisadas as distribuições de pressão nas superfícies da carroceria e o arrasto aerodinâmico. Para o teste computacional, foi utilizado um software de dinâmica dos fluidos computacional em que as equações de Navier-Stokes com média de Reynolds são resolvidas pelo método dos volumes finitos usando um modelo de turbulência RNG 'capa' - 'épsilon' / Two main objectives in the study of commercial vehicle aerodynamics are the reduction in fuel consumption and the improvement in engine refrigeration efficiency. These objectives can be achieved through development of devices which vary the flow characteristics around the vehicle and also through modification of the shape of external surfaces. The slope of rear surfaces, for instance, has large influence over the wake turbulence which forms behind the vehicle. The use of computational tools yields cost reduction and greater flexibility in automotive aerodynamic analysis. There is still a need, however, for verification of results, generated by these computational tools, with the largest possible number of test cases so that the mathematical model is adequately chosen and adjusted. The objective of the present work is the verification of experimental and computational results in the development of methodologies aiming at reduction of fuel consumption and improvement in engine refrigeration efficiency. Experimental and computational aerodynamic results were compared for a commercial bus model with a 1:17.5 scale. The experiments were conducted in an open section wind tunnel where pressure distribution and aerodynamic drag were analyzed. The numerical analysis was conducted using computational fluid dynamics software which solves the Reynolds Averaged Navier-Stokes equations using the finite volume method with a RNG 'capa' - 'épsilon' turbulence model

Diversas aplicações industriais utilizam escoamentos de fluidos viscoelásticos, e em muitos casos é necessário saber se os escoamentos propagam-se no estado laminar ou no turbulento. Embora a hidrodinâmica de fluidos viscoelásticos seja fortemente afetada pelo balanço entre forças inerciais e elásticas no escoamento, o efeito da elasticidade sobre a estabilidade de escoamentos inerciais não foi completamente estabelecido. Neste trabalho, estuda-se o que ocorre durante a transição laminar-turbulenta, investigando a convecção de ondas de Tollmien-Schlichting para o escoamento incompressível, para um fluido viscoelástico, entre placas paralelas, utilizando a equação constitutiva Giesekus. Para isto, adotou-se a simulação numérica direta para verificar a estabilidade dos escoamentos à perturbações não estacionárias deste fluido. Experimentos computacionais para verificação do código foram realizados. Com os resultados numéricos obtidos, foi possível verificar e analizar a estabilidade de escoamentos utilizando-se o modelo não newtoniano Giesekus. / Several industrial applications use viscoelastic fluid flows, and it is necessary to know if the flows propagate in the laminar or turbulent state. Although the hydrodynamics of viscoelastic fluids is strongly affected by the balance between inertial and elastic forces in the flow, the effect of elasticity on the stability of inertial flows has not been completely established. In this work we study what happens during the laminar-turbulent transition, investigating the convection of Tollmien-Schlichting waves for the incompressible flow, for a viscoelastic fluid, between parallel plates, using the constitutive equation Giesekus. For this, the direct numerical simulation was used to verify the stability of the flows to the non-stationary perturbations of this fluid. Computational experiments to verify the code were performed. With the numerical results obtained, it was possible to verify and analyze the stability of flows modelled by Giesekus non-newtonian model.

Dois dos principais objetivos no estudo da aerodinâmica de veículos comerciais são a redução no consumo de combustível e o aumento na eficiência da refrigeração do motor. Esses objetivos podem ser alcançados através do desenvolvimento de dispositivos que modifiquem o escoamento do ar ao redor do veículo e também através da alteração da forma das superfícies externas. A inclinação das superfícies da parte traseira de um ônibus, por exemplo, tem grande influência sobre a esteira turbulenta que se forma atrás do veículo. O uso de ferramentas computacionais permite uma redução de custo e maior flexibilidade na análise aerodinâmica de autoveículos. Ainda é preciso, no entanto, que o resultado dessas ferramentas computacionais seja verificado com o maior número possível de casos para que se possa escolher e ajustar o modelo matemático de forma adequada. O objetivo do presente trabalho é a verificação dos resultados computacionais e experimentais no desenvolvimento de metodologias que visem à redução no consumo de combustível e aumento na eficiência da refrigeração do motor. Foram comparados resultados experimentais e computacionais do escoamento sobre um modelo de um ônibus comercial em escala 1:17,5. Para a realização do experimento foi utilizado um túnel de vento de seção aberta, onde foram analisadas as distribuições de pressão nas superfícies da carroceria e o arrasto aerodinâmico. Para o teste computacional, foi utilizado um software de dinâmica dos fluidos computacional em que as equações de Navier-Stokes com média de Reynolds são resolvidas pelo método dos volumes finitos usando um modelo de turbulência RNG 'capa' - 'épsilon' / Two main objectives in the study of commercial vehicle aerodynamics are the reduction in fuel consumption and the improvement in engine refrigeration efficiency. These objectives can be achieved through development of devices which vary the flow characteristics around the vehicle and also through modification of the shape of external surfaces. The slope of rear surfaces, for instance, has large influence over the wake turbulence which forms behind the vehicle. The use of computational tools yields cost reduction and greater flexibility in automotive aerodynamic analysis. There is still a need, however, for verification of results, generated by these computational tools, with the largest possible number of test cases so that the mathematical model is adequately chosen and adjusted. The objective of the present work is the verification of experimental and computational results in the development of methodologies aiming at reduction of fuel consumption and improvement in engine refrigeration efficiency. Experimental and computational aerodynamic results were compared for a commercial bus model with a 1:17.5 scale. The experiments were conducted in an open section wind tunnel where pressure distribution and aerodynamic drag were analyzed. The numerical analysis was conducted using computational fluid dynamics software which solves the Reynolds Averaged Navier-Stokes equations using the finite volume method with a RNG 'capa' - 'épsilon' turbulence model

Existe uma variedade de modelos de dispersão de poluentes, em geral, os modelos gaussianos são usados em todo o mundo por agências ambientais com intuito de regulação. O modelo CALPUFF é um deles. Neste estudo, _e avaliada a influência de escalas de tempo de descorrelação no sistema de modelagem CALPUFF, sob condições atmosféricas neutras. Para fazer isso uma nova parametrização das escalas de tempo descorrelação é proposta. Uma distribuição espectral de um perfil de velocidade euleriana e uma for- mutação da evolução das escalas temporais de descorrelação Lagrangiana são utilizadas como os mecanismos forçantes na dispersão turbulenta, numa camada limite dominada pelo cisalhamento do vento. O desempenho do modelo foi estabelecido com a comparação das concentrações superficiais do experimento Over-Land Alongwind Dispersion. Emissões de fontes em linha foram avaliadas com o modelo CALPUFF com distintas formas de inicialização. Um segundo modelo também foi testado, normalmente utilizado para estudar e prever o impacto ambiental e validar parametrizações turbulentas. É o modelo estocástico de partículas Lagrangiano LAMBDA. Também neste estudo, os resultados do modelo LAMBDA e CALPUFF enfatizam a capacidade da nova derivação de escalas de tempo em representar o comportamento estocástico desconhecido do fenômeno da dispersão de poluentes. / There exists a variety of pollution of dispersion models and in general, Gaussian models are used worldwide by environmental agencies in regulatory applications. The CALPUFF model is one of them. In this study, the in uence of decorrelation time scales in the CALPUFF modeling system under neutral conditions is evaluated. To do this a new parameterization of decorrelation time scales is proposed. A spectral distribution of an Eulerian velocity pro le and a formulation of the evolution of the Lagrangian decorrela- tion timescales are used as the forcing mechanisms (shear-dominated boundary layer) for the turbulent dispersion. The model performance was established by comparing ground- level concentrations with Over-Land Alongwind Dispersion experimental results. Line source emissions was evaluated using the CALPUFF model with different forms of the initialization. A second model was also tested, normally used to study and predict the environmental impact and validate turbulent parameterizations. Is the stochastic La- grangian dispersion model LAMBDA (Ferrero and Anfossi, 1998). Also in this study the model LAMBDA and CALPUFF results emphasized the ability of the new derivation of decorrelation time scales to represent the unknown stochastic behavior.

Nos escoamentos em turbomáquinas temos como principais características a tridimensionalidade, possível ocorrência de separação da camada limite, relaminarização, transição laminar-turbulenta, dentre outros efeitos físicos. De acordo com alguns estudos experimentais em turbinas observouse que a transição laminar-turbulenta pode se estender por até 60% da corda de uma pá de turbina. Uma boa estimativa para se prever corretamente o local da transição é indispensável para que seja obtida uma melhoria na eficiência das turbinas. Escoamentos sobre superfícies côncavas estão sujeitos à instabilidade centrífuga, podendo dar origem a vórtices longitudinais, conhecidos como vórtices de Görtler. Esses vórtices são responsáveis por gerar distorções fortes nos perfis de velocidade e consequentemente nos perfis de temperatura. O presente estudo tem por objetivo estudar a influência da variação da curvatura de uma superfície côncava, e os efeitos do comprimento de onda transversal no processo de transição, e sua influência nas taxas de transferência de calor. Para tal, um código de simulação numérica paralelizado, com alta ordem de precisão, foi utilizado para resolver numericamente as equações de Navier-Stokes. Este código é validado através de comparações entre resultados obtidos com uso da teoria de estabilidade linear, e com resultados de simulações numéricas não lineares. Resultados obtidos evidenciam a influência da variação da curvatura, e os efeitos causados pelo comprimento de onda transversal nas instabilidades de Görtler, e secundária. Tais evidências comprovam que a variação da curvatura pode ser útil no controle do processo de transição laminar-turbulenta, e que as taxas de transferência de calor de um escoamento de Görtler desenvolvido em superfícies de curvatura variável podem ser intensificadas, atingindo valores superiores aos obtidos em escoamentos turbulentos. / Some characteristics of flows in turbomachinery are the three-dimensionality, possible occurrence of separation of the boundary layer, relaminarization, laminar-turbulent transition, among other physical effects. According to some experimental observations in turbines, it has been observed that the laminar-turbulent transition can extend over 60% chord of a turbine blade. A good estimate to correctly predict the location of the transition is essential for an improvement in the efficiency of turbines. Flow over concave surfaces is subjected to centrifugal instability, which may lead to formation of longitudinal vortices, known as the Görtler vortices. These vortices are responsible for generating strong distortions in the velocity profiles and hence the temperature profiles. The current goal aims to study the influence of the curvature variation of a concave surface and the effects of spanwise wavelength on the transition process and its influence on the heat transfer rates. For this, a parallel numerical simulation code, with a high order of precision, was used to numerically solve the Navier-Stokes equations. This code is validated through comparisons between results obtained using linear stability theory, and nonlinear numerical simulations results. Results obtained show the influence of the curvature variation, and the effects caused by the spanwise wavelength on the Görtler and secondary instabilities. This evidence proves that the curvature variation can be useful in the control of the laminar-turbulent transition process, and that heat transfer rates of a Görtler flow developed on variable curvature surfaces can be intensified, and reach values higher than these achieved in turbulent flows.

Diversas aplicações industriais utilizam escoamentos de fluidos viscoelásticos, e em muitos casos é necessário saber se os escoamentos propagam-se no estado laminar ou no turbulento. Embora a hidrodinâmica de fluidos viscoelásticos seja fortemente afetada pelo balanço entre forças inerciais e elásticas no escoamento, o efeito da elasticidade sobre a estabilidade de escoamentos inerciais não foi completamente estabelecido. Neste trabalho, estuda-se o que ocorre durante a transição laminar-turbulenta, investigando a convecção de ondas de Tollmien-Schlichting para o escoamento incompressível, para um fluido viscoelástico, entre placas paralelas, utilizando a equação constitutiva Giesekus. Para isto, adotou-se a simulação numérica direta para verificar a estabilidade dos escoamentos à perturbações não estacionárias deste fluido. Experimentos computacionais para verificação do código foram realizados. Com os resultados numéricos obtidos, foi possível verificar e analizar a estabilidade de escoamentos utilizando-se o modelo não newtoniano Giesekus. / Several industrial applications use viscoelastic fluid flows, and it is necessary to know if the flows propagate in the laminar or turbulent state. Although the hydrodynamics of viscoelastic fluids is strongly affected by the balance between inertial and elastic forces in the flow, the effect of elasticity on the stability of inertial flows has not been completely established. In this work we study what happens during the laminar-turbulent transition, investigating the convection of Tollmien-Schlichting waves for the incompressible flow, for a viscoelastic fluid, between parallel plates, using the constitutive equation Giesekus. For this, the direct numerical simulation was used to verify the stability of the flows to the non-stationary perturbations of this fluid. Computational experiments to verify the code were performed. With the numerical results obtained, it was possible to verify and analyze the stability of flows modelled by Giesekus non-newtonian model.

A metodologia LES (Large Eddy Simulation) é uma alternativa viável para a solução numérica de escoamentos de interesse prático em virtude da limitação computacional imposta pela resolução direta de todas as escalas presentes em escoamentos turbulentos. Entretanto, a compreensão detalhada do fenômeno da turbulência é ainda uma tarefa desafiadora em consequência do seu comportamento não linear e alta sensibilidade às condições iniciais e de contorno. Dessa forma, o sucesso de simulações LES está associado à utilização de um código computacional eficiente, com modelagem submalha que represente corretamente a dinâmica do escoamento, juntamente com a especificação de condições iniciais turbulentas fisicamente consistentes. Nesse contexto, o presente trabalho tem como objetivo o desenvolvimento de um código LES de alta ordem aliado a um método de geração de perturbações para o estudo de escoamentos turbulentos em camada limite sobre superfície plana. Foi adotada a formulação vorticidadevelocidade. A metodologia numérica baseia-se no método de diferenças finitas em malhas colocalizadas, onde as derivadas nas direções longitudinal e normal ao escoamento são aproximadas usando diferenças compactas de alta ordem. Esse estudo assume periodicidade na direção transversal do escoamento e então um método espectral é adotado nessa direção. A integração temporal é feita através do método Runge-Kutta de 4a ordem e a solução da equação de Poisson se dá por meio de um método multigrid. Para a modelagem submalha é adotado o modelo WALE (Wall-Adapting Local Eddy-viscosity). O método RFG (Random Flow Generation) foi responsável pela geração das flutuações de velocidade. Os resultados obtidos mostraram-se em boa concordância com os dados DNS (Direct Numerical Simulation) e LES presentes na literatura. / LES methodology is a viable alternative for the numerical solution of practical interest flows due to the computational limitations imposed by the direct resolution of all scales presented in turbulent flow. However, the detailed understanding of the turbulence phenomenon is still a challenging task as a result of its non-linear behavior and high sensitivity to initial and boundary conditions. Thus, the success of LES simulations is associated with the use of an efficient computational code, wherein the subgrid scale modeling accurately represents the flow dynamics, together with the specification of realistic inicial boundary conditions. In this context, this study aims to develop a high-order LES code combined with a method for generating velocity fluctuations to compute turbulent boundary layer flows over a flat plate. The vorticity-velocity formulation was adopted. The numerical scheme is based on the finite difference method in collocated grid, where the derivatives in the streamwise and wall-normal are approximated using high order compact finite difference schemes. We also assume periodicity in spanwise direction therefore it is adopted a spectral method in this direction. The method chosen for the temporal evolution is the 4th order Runge-Kutta method and the solution of Poisson equation solution is accessed via a multigrid algorithm. For subgrid modelling it is adopted the Wall-Adapting Local Eddy-viscosity (WALE) model. The RFG (Random Flow Generation) method was responsible for the generation of unsteady turbulent velocity signal. The results obtained were in good agreement with DNS (Direct Numerical Simulation) and LES from the literature.

[pt] Entender o processo de produção de fuligem é crucial para o projeto de novos queimadores, como os de fornos industriais. Estes queimadores, que utilizam processos de combustão turbulenta, dependem de transferência de calor via radiação das partículas de fuligem para as paredes do forno para seu bom funcionamento. A fuligem formada na região de radiação deve ser oxidada para evitar problemas de saúde e meio ambiente. Mesmo tendo havido significativo progresso no decorrer das duas últimas décadas em relação às chamas laminares, a interação entre a turbulência e a produção de fuligem ainda é um problema em aberto. Este trabalho apresenta resultados experimentais recentes da distribuição instantânea e média da distribuição de fuligem em chamas turbulentas de etileno/ar não prémisturadas estabilizadas em um queimador do tipo bluff-body. A intensidade de turbulência na região de esteira deste queimador é muito alta, levando a uma presença de fuligem intermitente e a estruturas de fuligem altamente distorcidas. A distribuição de fuligem é medida usando incandescência induzida por laser (LII), com uma excitação em 266 nm a 10 Hz e fluência de 0,8 J/cm2 e medição em 400 nm por uma câmera intensificada. Os resultados da técnica LII são comparados à técnica clássica de extinção da luz. Resultados da distribuição de hidrocarbonetos aromáticos policíclicos também são apresentados. Os resultados obtidos permitem caracterizar a distribuição da função de densidade de probabilidade de fuligem. Nas situações de escoamento onde a turbulência da esteira é controlada pelo escoamento de ar, demonstra-se que a PDF da fração volumétrica de fuligem corresponde a uma distribuição lognormal. / [en] Understanding the soot production process is crucial to the design of new burners, such as those in industrial furnaces. Indeed, these burners, which use turbulent combustion processes, rely on radiative heat transfer from the soot particles to the furnace walls to operate properly. The soot formed within the radiation region must the be oxidized in order to avoid health and environment issues. Although there has been significant progress over the past two decades in relation to laminar flames, the interaction between turbulence and soot production is still an open problem. This works presents recent experimental results of the instantaneous and mean soot distribution on non-premixed turbulent ethylene/air flames stabilized at a bluff-body burner. The turbulence intensity in the wake region of this burner is very high, leading to a soot intermittent presence and to highly distorted soot structures. The soot distribution is measured using laserinduced incandescence (LII), with 266 nm excitation at 10 Hz, 0.8 J/cm2 fluence and collected at 400 nm by an intensified camera. The results of the LII technique are compared to those of a classical of light extinction technique. Polyciclic aromatic hydrocarbon distribution results are also presented. The results obtained allow to characterize the soot probability density function distribution. In flow situations where the wake turbulence is controlled by the air flow, the soot volume fraction PDF is shown to correspond to a lognormal distribution.

A metodologia LES (Large Eddy Simulation) é uma alternativa viável para a solução numérica de escoamentos de interesse prático em virtude da limitação computacional imposta pela resolução direta de todas as escalas presentes em escoamentos turbulentos. Entretanto, a compreensão detalhada do fenômeno da turbulência é ainda uma tarefa desafiadora em consequência do seu comportamento não linear e alta sensibilidade às condições iniciais e de contorno. Dessa forma, o sucesso de simulações LES está associado à utilização de um código computacional eficiente, com modelagem submalha que represente corretamente a dinâmica do escoamento, juntamente com a especificação de condições iniciais turbulentas fisicamente consistentes. Nesse contexto, o presente trabalho tem como objetivo o desenvolvimento de um código LES de alta ordem aliado a um método de geração de perturbações para o estudo de escoamentos turbulentos em camada limite sobre superfície plana. Foi adotada a formulação vorticidadevelocidade. A metodologia numérica baseia-se no método de diferenças finitas em malhas colocalizadas, onde as derivadas nas direções longitudinal e normal ao escoamento são aproximadas usando diferenças compactas de alta ordem. Esse estudo assume periodicidade na direção transversal do escoamento e então um método espectral é adotado nessa direção. A integração temporal é feita através do método Runge-Kutta de 4a ordem e a solução da equação de Poisson se dá por meio de um método multigrid. Para a modelagem submalha é adotado o modelo WALE (Wall-Adapting Local Eddy-viscosity). O método RFG (Random Flow Generation) foi responsável pela geração das flutuações de velocidade. Os resultados obtidos mostraram-se em boa concordância com os dados DNS (Direct Numerical Simulation) e LES presentes na literatura. / LES methodology is a viable alternative for the numerical solution of practical interest flows due to the computational limitations imposed by the direct resolution of all scales presented in turbulent flow. However, the detailed understanding of the turbulence phenomenon is still a challenging task as a result of its non-linear behavior and high sensitivity to initial and boundary conditions. Thus, the success of LES simulations is associated with the use of an efficient computational code, wherein the subgrid scale modeling accurately represents the flow dynamics, together with the specification of realistic inicial boundary conditions. In this context, this study aims to develop a high-order LES code combined with a method for generating velocity fluctuations to compute turbulent boundary layer flows over a flat plate. The vorticity-velocity formulation was adopted. The numerical scheme is based on the finite difference method in collocated grid, where the derivatives in the streamwise and wall-normal are approximated using high order compact finite difference schemes. We also assume periodicity in spanwise direction therefore it is adopted a spectral method in this direction. The method chosen for the temporal evolution is the 4th order Runge-Kutta method and the solution of Poisson equation solution is accessed via a multigrid algorithm. For subgrid modelling it is adopted the Wall-Adapting Local Eddy-viscosity (WALE) model. The RFG (Random Flow Generation) method was responsible for the generation of unsteady turbulent velocity signal. The results obtained were in good agreement with DNS (Direct Numerical Simulation) and LES from the literature.