O enchimento de reservatórios hidrelétricos provoca impactos ambientais, interferindo nos ecossistemas e nos modos de vida das populações envolvidas, pois durante este processo a biomassa terrestre é decomposta lançando substâncias que se concentram nos volumes dos diferentes compartimentos do reservatório. Nestas condições, algumas regiões do reservatório passam por períodos em que os teores de oxigênio dissolvido e a concentração de matéria orgânica comprometem o equilíbrio da flora e fauna locais, afetando tanto a qualidade das águas como também os ecótonos vizinhos. Neste trabalho, será tratado o problema do transporte das concentrações das espécies constituintes, da temperatura e da turbulência, 2D e 2DH, durante o enchimento de compartimentos de reservatórios / The wadding of hidroelectric reservoirs provokes ambient impacts, intervening with the ecosystems and the ways of life of the involved populations, because during this process the terrestrial biomass is decomposed launching substances that stay concentrated in the volumes of the different compartments of the reservoir. In these conditions, some regions of the reservoir pass for periods where the texts of dissolved oxygen and the concentration of organic substance compromise the balance of the local flora and fauna, affecting as much the quality of waters as the neighboring ecotones. In this work, will be treated the problem of the transport of the concentrations of the constituent species, of the temperature and the turbulence, 2D and 2DH, during the wadding of compartments of reservoirs

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The present work is dedicated to perform the mathematical modeling and DNS and LES simulations of a three-dimensional, temporally evolving incompressible plane wake are performed, seeking to evidence differences in stability, transition and onset of both coherent and small scale structures, when the flow is subjected to random perturbations of different amplitudes. The perturbations are generated using the Random-Flow-Generator (RFG) technique, being imposed to the flow as initial conditions. The Navier-Stokes equations are solved in a prismatic domain, with periodic boundary conditions in all directions, using Fourier pseudospectral method. The invariants of the velocity gradient tensor, Q and R, are analyzed for random perturbations with magnitudes 10−3, 10−4 and 10−5, showing the onset of their characteristic teardrop correlation map. Moreover, maps of the second and third invariants of the rate-of-strain tensor, QS and RS, are shown, in order to evidence the differences in local flow strain and topological characteristics of the dissipation of kinetic energy. Isosurface plots of Q and QW, as well as vorticity contours are shown, allowing visual identification of the coherent structures and confirming patterns predicted by the invariant maps. / O presente trabalho se dedica a modelagem matemática e a simulações numéricas DNS e LES de uma esteira tridimensional, incompressível, em desenvolvimento temporal, buscando evidenciar diferenças na estabilidade, transição e no desenvolvimento de estruturas coerentes e de pequena escala, quando o escoamento é submetido a perturbações randômicas de diferentes amplitudes. As perturbações são geradas utilizando-se a técnica Random Flow Generator (RFG), sendo sobrepostas à condição inicial do escoamento. As equações de Navier-Stokes são resolvidas em um domínio prismático, com condições de contorno periódicas em todas as direções, utilizando-se o método pseudoespectral de Fourier. Os invariantes do tensor gradiente de velocidade, Q e R, são analisados para perturbações de magnitude 10−3, 10−4 and 10−5, mostrando a formação de uma correlação no formato de gota, característica da resolução das equações de Navier-Stokes. Além disso, são apresentados mapas do segundo e terceiro invariante do tensor taxa de deformação, QS e RS, a fim de evidenciar as diferenças locais no escoamento e as características topológicas na taxa de dissipação de energia cinética. Isosuperfícies de Q e QW, bem como contornos de vorticidade são apresentados, possibilitando a identificação visual das estruturas coerentes, e confirmando os padrões de estruturas previstos pelos mapas de invariância. / Mestre em Engenharia Mecânica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this work a single column model is implemented, aiming primarily to simulate the Stable Boundary Layer (SBL). The model is composed of prognostic equations for the potential temperature, specific humidity and the wind components. A prognostic equation for the turbulent kinetic energy (TKE) is also used. The turbulent fluxes are approximated from diffusion coefficients that are directly proportional to TKE and to the turbulent mixing length, with a adjustable proportionality factor. To estimate the turbulent mixing lengths, four different parametrizations from the literature are considered. First, a stable case in the Arctic is simulated, then the model is adjusted to reproduce previously published Large Eddy Simulations (LES) results for the same case. After the calibration, the model is validated through a comparison with observed nocturnal data, from tethered balloon soundings performed at a deforested area in the Amazon region, in the months of July and October 2001 and November 2003. The temporal evolution of potential temperature and specific humidity vertical profiles are simulated and compared to the observations from 13 different nights. The performance of the different formulations is evaluated through a statistic analysis for all nights. Potential temperature profiles were better represented by the model than specific humidity ones. In spite being able to reproduce different stability conditions, the model represented the most stable nights better. Finally, the model was coupled to a surface energy budget scheme and for these simulations the initial conditions from the night of 28 July 2001 were used. The scheme was used to evaluate how the simulated atmospheric variables depend on the wind intensity at the top of the vertical domain for the different formulations and with or without the energy budget scheme being added. / Neste trabalho é implementado um Modelo de Coluna Simples, com o intuito de simular a Camada Limite Estável (CLE). O modelo é composto por equações prognósticas para a temperatura potencial, umidade específica e para as componentes horizontais do vento, além de uma equação prognóstica para a energia cinética turbulenta (ECT). Os fluxos turbulentos são aproximados por coeficientes de difusão que são diretamente proporcionais à ECT e ao comprimento de mistura turbulenta, de forma que o coeficiente de proporcionalidade é um fator ajustável. Para estimar os comprimentos de mistura turbulentos são utilizadas quatro diferentes parametrizações encontradas na literatura. Primeiramente, é simulado um caso estável no Ártico, sendo o modelo ajustado de forma a produzir resultados próximos aos da Simulação dos Grandes Turbilhões (LES), obtidos por outros autores para o mesmo caso. Realizados os ajustes, o modelo é validado através da comparação com dados noturnos, obtidos por sondagens de balão cativo realizadas em um sítio experimental localizado numa área desmatada da região Amazônica, nos meses de julho e outubro de 2001 e de novembro de 2003. São simulados os perfis verticais de temperatura potencial e de umidade específica e comparados com as observações de 13 diferentes noites. A avaliação do desempenho das diferentes formulações utilizadas é feita através de uma análise estatística, utilizando dados de todas as noites. Observou-se, que os perfis verticais de temperatura potencial foram melhor representados do que os de umidade específica. Embora o modelo tenha sido capaz de reproduzir as diferentes condições de estabilidade, as noites mais estáveis foram melhor simuladas. Finalmente, é acrescentado um esquema que resolve o balanço de energia em superfície ao modelo, sendo utilizadas para este fim as condições iniciais da noite de 28 de julho de 2001. Este esquema foi utilizado para avaliar como as variáveis atmosféricas simuladas dependem da intensidade do vento no topo do domínio vertical, para as diferentes formulações analisadas e com ou sem a utilização do balanço de superfície.

Made available in DSpace on 2016-06-02T19:55:37Z (GMT). No. of bitstreams: 1 5587.pdf: 3066852 bytes, checksum: 17d88d3ab913d61bb315842b57ec5a74 (MD5) Previous issue date: 2013-09-30 / Universidade Federal de Sao Carlos / Cyclones are used to separate suspended particles from gas streams. The analysis of flow inside the cyclone is complex due to the large number of parameters and operating variables that influence the dynamics of the system. Due to its versatility of use and its robustness, results in computational fluid dynamic are an important alternative often used to study the dynamics of flow in the cyclone. In the literature, several works make use of computational fluid dynamics in order to study the pressure drop of cyclones. However, the works about the efficiency of collection are fewer. For the cyclone, computational fluid dynamics allows an accurately simulation of the collection efficiency of particles with diameter greater than 5 µm; but for particles with smaller diameters the simulated results diverge from experimental values. In this work, Stairmand s cyclone was numerically simulated with injected particles with diameters from 1 to 5 µm. The results were verified using available experimental data reported in the literature. The experimental data taken from literature belong to Zhao and Shen Kang (2004) and Zhao (2005). The simulation mesh was analyzed before and after the simulations. The turbulence models used in the simulation of flow in the cyclone were: Reynolds Stress Model (RSM) and Large Eddy Simulation (LES). The dispersed phase was simulated considering one and two way couplings. The equation of motion of the particles was integrated using the integration schemes: implicit, analytic, trapezoidal, and Runge-Kutta. The results showed that the methodology was adequate to reproduce the behavior of the flow in the cyclone. The error percentages obtained in pressure drop were under 5% and the average collection efficiency was reproduced with good accuracy for diameters of 3, 4 and 5 µm. / O ciclone é um equipamento utilizado para separar partículas suspensas em correntes gasosas. A análise do escoamento dentro do ciclone é complexa devido ao grande número de parâmetros e variáveis operacionais que influenciam na dinâmica do sistema. Devido a sua versatilidade de uso e a sua robustez nos resultados a fluidodinâmica computacional é uma alternativa importante utilizada com frequência no estudo da dinâmica do escoamento no ciclone. Na literatura existem vários trabalhos que fazem uso da fluidodinâmica computacional no estudo da queda de pressão. Entretanto, no estudo da eficiência de coleta, os artigos são em menor quantidade. No ciclone, a fluidodinâmica computacional permite simular a eficiência de coleta de partículas com diâmetro maiores que 5 µm com bastante precisão; a modelagem da eficiência de coleta para partículas com diâmetros menores ainda resultam em valores diferentes dos valores experimentais. Neste trabalho foram realizadas as simulações numéricas no ciclone tipo Stairmand; foram injetadas partículas de 1 a 5 µm. Para a verificação foram tomados dados experimentais disponíveis na literatura. Os dados experimentais foram obtidos por Zhao, Shen e Kang (2004) e Zhao (2005). A malha utilizada nas simulações, foi analisada antes e depois das simulações. Os modelos de turbulência utilizados para simular o escoamento no ciclone foram: Reynolds Stress Model (RSM) e Large Eddy Simulation (LES). A fase dispersa foi simulada considerando acoplamento de uma via e acoplamento de duas vias. Para integrar as equações de movimento das partículas, foram utilizados os esquemas de integração: implícito, analítico, trapezoidal e Runge-Kutta. Os resultados mostraram que a metodologia empregada mostrou-se adequada para reproduzir o comportamento do escoamento no ciclone. As porcentagens de erro obtidos na queda de pressão foram em média menores que 5% e na eficiência de coleta foram reproduzidas com boa precisão para os diâmetros de 3, 4 e 5 µm.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / The Amazon river, located in northernBrazil, discharges between 80,000 and 250,000 m3s-1 of water onto the adjacent shelf, creating a plume of brackish water that extends hundreds of kilometers away from the river mouth. This river also carries a large amount of fine sediments to the ocean where fluid mud has been found in the topset and upper foreset layers of the subaqueous delta formed on the mid-shelf. One of the main goals of this dissertation is to describe how turbulence and suspended sediment concentration vary along the Northern Channel of the Amazon river. Water column measurements were carried out in October 2008 at six anchor stations (P1, P3, P5, P6, P8 e P9) located seaward of the river mouth; P1 and P9 were 125 km apart. Each station was occupied during 13 hours during which current speed and direction were continuously sampled with a 600 kHz Teledyne-RDI ADCP; hourly profiles of temperature, salinity, turbidity and depth were also obtained. Water samples were collected for determination of Suspended Particulate Matter (SPM) concentration and calibration of the turbidity sensor. Current speed reached values above 1.5 m s1 in the along-channel direction (NE-SW); a remarkable ebb-flood asymmetry was observed and flows were strongly ebb-dominated. Throughout the water column, SPM concentration at stations P1 and P3 varied between 100 and 300 mg L1 in association with the presence of freshwater. In contrast, a strong salinity gradient was observed between stations P6 and P9, coinciding with the occurrence of concentrations of SPM above 10 g L-1 (fluid mud). At stations P3, P5 and P6, interface between freshwater from the Amazon river and salt water from the continental shelf, shear stresses wereestimated through four diferents methods: Reynolds, Turbulent Kinetic Energy (TKE), modified TKE and Quadratic Law; in the nearbed region (3 mab) the computed values varied between 0 and 3 Pa. At the three stations (P3, P5 and P6) the lowest and the highest shear stress values were obtained through, respectively, the Reynolds and the TKE methods. Over the whole water column turbulence intensity was estimated through the standard deviation of the turbulent component of the along-channel current velocity (root-mean square of u); from these values, it was estimated the turbulent dissipation of energy (G), whose values at 3 mab varied between zero and 20 s1.

Orientadores: Jose Manoel Balthazar, Paulo R. G. Kurka / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-09T22:56:56Z (GMT). No. of bitstreams: 1 Pereira_DaniloCarlos_D.pdf: 3711639 bytes, checksum: 2413d33f619760c04be8cc5320c0a84b (MD5) Previous issue date: 2007 / Resumo: Neste trabalho analisou-se a dinâmica não linear de uma aeronave em vôo longitudinal. Efetuou-se a análise do comportamento bifurcacional da aeronave F-8 ¿Cruzader¿. Na análise bifurcacional foi estudado o comportamento topológico desta aeronave tomando-se dois parâmetros de controle: a deflexão do profundor e a alteração da massa da referida aeronave. Ante a pesquisa desenvolvida, foi proposto um projeto de controle linear ótimo com o objetivo de estabilizar as oscilações do ângulo de ataque, considerando-se regiões criticas do comportamento não linear da aeronave. Adicionalmente, incluiu-se no modelo matemático a variação da velocidade longitudinal da aeronave, visto tratar-se de simulações numéricas em um túnel de vento virtual / Abstract: In this work it was analyzed the non linear dynamic of an aircraft taken onto longitudinal flight. It was done analysis of the bifurcacional behavior of the aircraft F-8 ¿Cruzader¿. In the bifurcational analysis was studied the topological behavior of this aircraft taken into account two parameters of control: the deflection of the elevator and the alteration of the mass of the related aircraft. In the face of the developed research, an optimum linear control project was proposed with the objective of stabilizing the oscillations of the angle-of-attack. Additionally, the variation of the longitudinal speed of the aircraft was included in the mathematic model in order to simulate the oscillatory movement of the aircraft considered, in a tunnel of virtual wind / Doutorado / Materiais e Processos de Fabricação / Mestre em Engenharia Mecânica

In this work, a DNS code was developed to investigate problems on transition in subsonic compressible boundary layer on a flat plate. Code validation tests were performed for linear and nonlinear stages of transition, on incompressible and compressible regimes. The focus of the present work is to investigate natural transition in subsonic boundary layers modeled by wave packets; and perform a preliminary study of transition induced by white noise. Three main problems were considered, namely, a DNS simulation and analysis of the ex- periment (MEDEIROS; GASTER, 1999b) of wave packet evolution on incompressible boundary layer, the influence of compressibility on wave packet evolution at subsonic Mach numbers and finally, a preliminary study of the evolution of a white noise perturbation in the boundary layer at Mach 0.2 and Mach 0.9. Comparisons between numerical and experimental results show remarkably good agreement in the linear and nonlinear stages, in both, spatial and Fourier spaces. A numerical simulation of this experiment and the analysis carried out is not available in the literature for wave packets in the incompressible boundary layer. The nonlinear modal analysis performed established the existence of tuned fundamental and subharmonic resonance of H-type and K-type in the packet. Influence of compressibility in the wave packet evolution was here investigated in boundary layers at Mach 0.7 and Mach 0.9. There are no works reported in the literature on wave packets in compressible subsonic boundary layer. In the linear regime, the oblique modes were the most unstable for Mach > 0.7, as expected by the results of the literature. In the nonlinear regime, strong streaks were observed, associated with low frequency modes that eventually decay downstream. An isolated wave packet at Mach 0.9 showed nonlinear amplification only in the subharmonic band, which may be associated to H-type or detuned resonance. However this packet has a relatively stable character. On the other hand, at Mach 0.9 spanwise interaction of wave packet pairs were more unstable than the isolated case, because stable modes for the isolated packet evolution becomes unstable in the wave packet interaction. This scenario evidenced the presence of oblique transition. Finally, the nonlinear evolution of the same white noise disturbance at Mach 0.2 and Mach 0.9 were observed to be completely different. In the incompressible boundary layer localized lambda vortex structures were observed, that could be associated to the local presence of H-type and/or K-type resonance. In the compressible regime, longitudinal vortex structures distributed across the entire domain seemed to be linked to oblique transition. In the white noise evolution, compressibility seems to have a stronger effect than in the wave packet evolution. In the conditions considered, the wave packet interaction appear to be a better representation of white noise compressible transition scenario. / No presente trabalho, um código DNS (Direct Numerical Simulation) foi desenvolvido para abordar problemas de transição para turbulência em camada limite subsônica compressível em uma placa plana. Foram realizados testes de validação de código , nos regimes linear e não linear do processo de transição, nos regimes incompressível e compressível. O foco do presente trabalho é estudar transição natural modelada por meio de pacotes de onda em camada limite compressível subsônica, e realizar uma análise preliminar da transição induzida por ruído branco. Três assuntos principais foram considerados: uma simulação DNS e uma análise comparativa com o experimento (MEDEIROS; GASTER, 1999b) sobre a evolução de um pacote de ondas em camada limite incompressível, a influência da compressibilidade na evolução de pacotes de ondas no regime subsônico, e por último, um estudo preliminar da transição induzida por ruído branco em Mach 0.2 e Mach 0.9. As comparações realizadas entre a solução numérica e os dados experimentais mostram uma boa concordância, nos regimes linear e não linear, tanto no espaço físico quanto no espaço de Fourier. A simulação numérica deste experimento e a análise realizada neste trabalho, não são encontradas na literatura para o regime incompressível. A análise modal não linear aplicada aos resultados, permitiu identificar a presença das ressonâncias tipo H e tipo K no pacote de ondas. A influência da compressibilidade na evolução dos pacotes de onda foi estudada em Mach 0.7 e Mach 0.9. Na literatura não há trabalhos sobre pacotes de ondas no regime sub- sônico. No regime linear da transição, os modos oblíquos resultam ser os mais instáveis para Mach > 0.7, como era de esperar, de acordo com os resultados da literatura. No regime não linear, foram observadas estrias de moderada amplitude, associadas com modos de baixa frequência que acabam decaindo. O pacote de ondas em Mach 0.9 apresentou amplificação não linear somente na banda subharmônica, que pode ser associada com transição tipo H ou ressonância dessintonizada. No entanto, o comportamento geral neste regime é estabilizante. Por sua vez, a interação entre pacotes de ondas em Mach 0.9 mostrou um comportamento desestabilizante, pois a interação acaba gerando amplificação não linear em modos que decaem no pacote isolado. Os modos amplificados sugerem a presença do mecanismo de transição oblíqua. Finalmente, a evolução da mesma perturbação constituída por ruído branco em Mach 0.2 e Mach 0.9, resultaram ser completamente diferentes. Na camada limite incompressível foram observados vórtices tipo lambda, que poderiam ser gerados pela presença localizada das ressonâncias tipo H e/ou tipo K. No regime compressível foram observados vórtices distribuidos em todo o domínio, o que sugere a presença da transição oblíqua. Na transição gerada por ruído branco a compressibilidade teve uma influência maior que no pacote de ondas. Nas condições estudadas, a interação entre pacotes de ondas parece ser uma melhor representação do ruído branco no regime compressível.

Orientador: Jose Roberto Nunhez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-11T08:15:27Z (GMT). No. of bitstreams: 1 Lopes_GabrielaCantarelli_M.pdf: 1710566 bytes, checksum: efab0c6308734536b15c2ef91d4d8337 (MD5) Previous issue date: 2008 / Resumo: Processos que envolvem escoamentos de fluidos no interior de dutos estão presentes em muitas aplicações industriais. Na indústria petroquímica um desses processos que vem se tornando cada vez mais importante é o craqueamento catalítico de frações pesadas do petróleo, já que as converte em frações leves e mais nobres. Por causa da crescente demanda mundial por gasolina e GLP e à sua alta rentabilidade para uma refinaria de petróleo, unidades de craqueamento catalítico em leito fluidizado (FCC) estão presentes em todo o mundo. Testes experimentais costumam ser usados no estudo dos fenômenos envolvidos nesse processo. Porém esse tipo de análise possui alto custo, que pode ser reduzido com o uso de simulações computacionais em seu estudo inicial. Assim, o objetivo deste trabalho foi desenvolver um modelo tridimensional e elíptico, em linguagem Fortran, capaz de fornecer dados para a análise preliminar de escoamentos no interior de reatores de FCC. Na modelagem desses problemas são usadas equações diferenciais parciais, e essas não possuem solução analítica conhecida, sendo necessário o emprego de métodos numéricos para esse fim. Neste trabalho foi usado o Método dos Volumes Finitos, que tem a função de substituir as equações diferenciais parciais por equações algébricas aplicadas a pequenos volumes de controle finitos pertencentes ao domínio. Uma das maiores dificuldades encontradas no tratamento numérico de escoamentos incompressíveis é a determinação de um campo de pressão que satisfaça a Equação da Continuidade. Esse problema foi resolvido fazendo-se uso da abordagem acoplp.da de solução. Para análise do modelo foram obtidos perfis numéricos de velocidade e pressão para fluidos escoando em regimes laminar e turbulento, que foram validados usando-se os dados obtidos com a solução analítica das equações, por correlações (semi-) empíricas ou por dados experimentais, conforme cada um dos casos. Notou-se que o modelo representa muito bem casos laminares, e gue 11.oS casos turbulentos foi necessário um maior refino da malha próximo a parede do tubo. Também foram feitas simulações para que se pudessem observar as características tridimensionais, elípticas e transientes da modelagem do escoamento. De maneira geral o modelo se mostrou bastante rápido, convergindo em poucas iterações. Palavras-Chave: Fluidodinâmica computacional; dutos cilíndricos; modelo tridimensional e elíptico; método dos Volumes Finitos; solução acoplada; turbulência / Abstract: Processes involving fluid flow in tubes are present in many industrial applications. In petrochemical industry one of these processes that are becoming more and more important is the fluid catalytic cracking of heavy petroleum fractions. This fact is due to the process capacity to convert heavy fractions in light and valuable ones. Because of the increasingly worldwide demand for gasoline and LPG and its high yield for a petroleum refinery, fluid catalytic cracking (FCC) units are present in the whole world. Experimental tests are used in the study of the phenomena involved in this processo However this kind of analysis has high cost which can be reduced by using computational simulations in its initial study. Thus, the aim of this work was the development of a three-dimensional and elliptical mo dei in Fortran language in order to provide data for fluid flow preliminary analysis in FCC reactors. Partial differential equations were used in the modeling of these problems. These equations do not have known analytical solution, being necessary therefore the use of numerical methods. In this work the Finite Volume Method were applied with this purpose. This method has as a role to substitute the partial differential equations of the mo dei for algebric equations applied to small finite control volumes of the domain. One of the biggest difficulties found in the numerical treatment of incompressible fluid flows is the determination of apressure field that satisfies the Continuity Equation. This problem was solved using the coupled solution approach. For model analysis, numerical velocity and pressure proJ:iles for laminar and turbulent flows were obtained, that had been validated using the data obtained through the analytical solution of the equations, by empirical correlations or by experimental data, according to each one of the cases. The model represented well laminar cases, and in the turbulent ones the mesh had to be more refined near the tube wall. Other simulations were performed, in aQalyzing the three-dimensional, elliptical and transient model characteristics. In general, the mo dei was very fast, converging in a few interations. Keywords: Computational fluid dyn~mic; cylindrical ducts; three-dimensional and elliptical model; Finite Volume Method; coupled solution; turbulence / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Neste trabalho foi realizado um estudo utilizando os modelos de turbulência Simulação das Grandes Escalas, Large Eddy Simulation (LES), e Simulação dos Vórtices Desprendidos, Detached Eddy Simulation (DES), para simular o escoamento compressível interno em um duto contendo válvulas controladoras na saída dos gases de combustão de um reator de craqueamento catalítico fluido, com o objetivo de comparar o desempenho numérico e computacional de ambas as técnicas. Para isso foi utilizado um programa comercial de dinâmica dos fluidos computacional, Computational Fluid Dynamics (CFD), que possui em seu código os dois modelos de turbulência. / In the present work a study was conducted using Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) turbulence models in order to simulate the internal compressible flow in a duct containing the flue gas discharge control valves of a fluid catalytic cracking reactor so as to compare the numerical and computational behavior of both techniques. A commercial Computational Fluid Dynamics (CFD) software, which includes these turbulence models in its code, was used.

Este trabalho é dedicado à análise e implementação de esquemas "upwind" de alta ordem modernos e o modelo de turbulência k-epsilon padrão no Freeflow-2D; um ambiente integrado para simulação numérica em diferenças finitas de problemas de escoamentos incompressíveis com superfícies livres. O propósito do estudo é a simulação de escoamentos de fluidos newtonianos incompressíveis, bidimensionais, confinados e/ou com superfícies livres e a altos valores do número de Reynolds. O desempenho do código Freeflow-2D atual é avaliada na simulação do escoamento numa expansão brusca e de um jato livre incidindo perpendicularmente sobre uma superfície rígida impermeável. O código é então aplicado na simulação de um jato planar turbulento em uma porção de fluido com superfície livre e estacionário. Os resultados numéricos obtidos são comparados com dados experimentais, soluções analíticas e soluções numéricas de outros trabalhos. / This work is devoted to the analysis and implementation of modern high-order upwind schemes and the standard k-epsilon turbulence model into the Freeflow-2D; a finite difference integrated environment for the numerical simulation of incompressible free surface flow problems. The purpose of this study is the two-dimensional simulation of high-Reynolds incompressible newtonian confined and/or free surface flows. The performance of the current Freeflow-2D code is assessed by applying it to the simulation of flow over a backward facing step and of an impinging free jet onto an impermeable rigid surface. The code is then applied to a turbulent planar jet into a pool. The numerical results are compared with experimental data, analytical solution, and numerical simulations of other works.