A caracterização aerodinâmica de foguetes de sondagem é requisito fundamental para o projeto e operação destes veículos. O conhecimento das forças atuantes determina o carregamento a que sua estrutura estará submetida e assim orientará o seu projeto estrutural. Da mesma forma, para se determinar a sua trajetória, faz-se necessário o conhecimento das características aerodinâmicas em todas as fases do vôo. Uma das maneiras de se obter as características aerodinâmicas de um foguete é por meio de simulações numéricas do escoamento sobre o foguete. Para isso, a região em torno do foguete é discretizada sob a forma de uma malha computacional que será então empregada para a solução numérica das equações de governo do escoamento. Os objetivos da presente dissertação consistem em gerar malhas computacionais sobre configurações de foguetes, e avaliar a sua adequação para a simulação do escoamento nos diferentes regimes de velocidade. Estas simulações serão feitas empregando-se o programa de cálculo NS3D da subdivisão de Aerodinâmica da ASE, que resolve numericamente as equações de Navier-Stokes 3D de escoamentos turbulentos. Estas simulações serão validadas pela comparação dos resultados dos cálculos com dados experimentais.

A study of the flowfield around axisymmetric launch vehicles in different flight conditions and configurations is performed in this work. Particularly, the VLS second stage flight configuration is analyzed considering the case with and without a propulsive jet in the vehicle base. Among the different cases considered for this configuration there are inviscid as well as viscous turbulent flows. The study is performed using a finite volume cell centered formulation on unstructured grids. Different spatial discretization schemes are compared, including a centered and an upwind scheme. The upwind scheme is a second-order version of the Liou flux vector splitting scheme and a simple scalar advection test problem is used to assess the influence of different reconstruction and limiting methods in the second-order extension of the scheme. Turbulence effects are accounted for using two one-equation turbulence closure models, namely the Baldwin and Barth and the Spalart and Allmaras models. An agglomeration multigrid algorithm is used to accelerate the convergence to steady state of the numerical solutions. Mesh refinement procedures as well as hybrid and adaptive meshes are discussed.

O difusor de seção anular é frequentemente empregado em turbinas a gás de motores aeronáuticos, dada sua característica de alcançar elevada pressão em espaço físico limitado. Guias e normas que auxiliem no projeto otimizado desse tipo de difusor não contemplam as diversas condições de contorno e regimes de escoamento nos quais o dispositivo pode operar. Para tanto, técnicas de modelagem numérica têm sido aplicadas para predizer o comportamento do escoamento nesse dispositivo e o seu desempenho. Este trabalho propõe a aplicação de ferramenta CFD (Dinâmica dos Fluidos Computacional) para análise dos efeitos no desempenho do difusor tipo anular de núcleo paralelo devido às variações do seu ângulo de divergência e dos números de Reynolds e de Mach do escoamento. O programa FLUENT 13.0 é utilizado para execução das simulações de escoamentos em regimes incompressível e compressível, empregando-se o modelo de turbulência - realizável e tratamento de parede melhorado. As equações de conservação de massa, de quantidade de movimento e de energia são resolvidas pelo emprego do Método dos Volumes Finitos (MVF). Adotou-se a formulação bidimensional axissimétrica para o domínio computacional, assim como as hipóteses de escoamento de ar turbulento, subsônico, em regime permanente, ausência de swirl e ação gravitacional desprezível. Nas simulações de escoamento incompressível, as propriedades do ar são constantes e condição de contorno de velocidade é atribuída na entrada do domínio. Em regime compressível, o fluido é tratado como gás ideal e condição de vazão mássica e temperatura total são impostas na entrada. Na saída da geometria é especificada condição da atmosfera padrão. As paredes são tratadas como adiabáticas, não-deslizantes e não-rugosas. Estudo de refinamento de malha é realizado para garantir a independência da malha nos resultados. O método numérico é validado pela comparação dos resultados obtidos para o coeficiente de recuperação de pressão estática ( ) com dados de estudos experimentais disponíveis na literatura. Os resultados são apresentados para os difusores com ângulos de divergência de , e . Dos resultados, percebe-se que as recirculações se acentuam com o aumento do ângulo e que o fenômeno de separação de escoamento pode ocorrer de forma cada vez mais precoce, seguido de perda acentuada de pressão estática. Além disso, variações nas propriedades do fluido próximo à parede do difusor, quando sujeito a escoamento compressível, contribuem com o melhor desempenho do mesmo se comparado aos resultados do difusor sob regime incompressível.

A interação turbulência-radiação (TRI, do inglês Turbulence-Radiation Interaction) resulta do acoplamento altamente não linear entre flutuações da intensidade de radiação e flutuações da temperatura e da composição química do meio, e tem-se demonstrado experimentalmente, teoricamente e numericamente que este é um fenômeno relevante em diversas aplicações envolvendo altas temperaturas, especialmente em problemas reativos. Neste trabalho, o TRI é analisado em um escoamento interno não reativo de um gás participante que se desenvolve em um duto de seção transversal quadrada, para diferentes intensidades de turbulência do escoamento e considerando duas espécies distintas para a composição do fluido de trabalho (dióxido de carbono e vapor de água). O objetivo central é avaliar como a inclusão ou não da variação espectral das propriedades radiativas do meio no cálculo influencia a magnitude do TRI. Isso é feito através de simulações numéricas no código de dinâmica dos fluidos computacional Fire Dynamics Simulator (FDS), que resolve, através do método dos volumes finitos, as equações fundamentais que regem o problema – isto é, os balanços de massa, de quantidade de movimento e de energia e a equação de estado – em uma formulação adequada para baixos números de Mach, utilizando um algoritmo de solução explícito e de segunda ordem no tempo e no espaço. A turbulência é modelada através da simulação de grandes escalas (LES, do inglês Large Eddy Simulation), empregando-se o modelo de Smagorinsky dinâmico para o fechamento dos termos submalha; para a radiação térmica, o método dos volumes finitos é utilizado na discretização da equação da transferência radiativa e os modelos do gás cinza e da soma-ponderada-de-gases-cinza (WSGG, do inglês Weighted-Sum-of-Gray-Gases) são implementados como forma de desconsiderar e de incluir a dependência espectral das propriedades radiativas, respectivamente. A magnitude do TRI sobre o problema é avaliada através de diferenças entre as médias temporais dos fluxos de calor superficiais e do termo fonte radiativo obtidas em cálculos que consideram os efeitos do fenômeno e cálculos que os negligenciam. Em geral, a interação turbulência-radiação mostrou ser pouco importante em todos os casos considerados, o que concorda com resultados de outros estudos sobre o tema em escoamento não reativos. Com o modelo WSGG, as contribuições do fenômeno foram maiores do que com a hipótese do gás cinza, evidenciando que a inclusão da variação espectral na solução do problema radiativo tem um impacto sobre a magnitude dos efeitos do TRI. Além disso, é feita uma discussão, em parte inédita no contexto do TRI, sobre diferentes metodologias para a análise do fenômeno. Finalmente, é proposto um fator de correção para o termo fonte radiativo médio no modelo WSGG, que é validado através de sua implementação nos casos simulados. Em estudos futuros, uma análise de sensibilidade sobre os termos constituintes desse fator de correção pode levar a um melhor entendimento de como as flutuações de temperatura se correlacionam com o fenômeno da interação turbulência-radiação. / Turbulence-radiation interaction (TRI) results from the highly non-linear coupling between fluctuations of radiation intensity and fluctuations of temperature and chemical composition of the medium, and its relevance in a number of high-temperature problems, especially when chemical reactions are included, has been demonstrated experimentally, theoretically, and numerically. In the present study, the TRI is analyzed in a channel flow of a non-reactive participating gas for different turbulence intensities of the flow at the inlet and considering two distinct species for the medium composition (carbon dioxide and water vapor). The central objective is to evaluate how the inclusion or not of the spectral variation of the radiative properties of a participating gas in the radiative transfer calculations affects the turbulence-radiation interaction. With this purpose, numerical simulations are performed using the computational fluid dynamics Fortranbased code Fire Dynamics Simulator, that employs the finite volume method to solve a form of the fundamental equations – i.e., the mass, momentum and energy balances and the state equation – appropriate for low Mach number flows, through an explicit second-order (both in time and in space) core algorithm. Turbulence is modeled by the large eddy simulation approach (LES), using the dynamic Smagorinsky model to close the subgrid-scale terms; for the thermal radiation part of the problem, the finite volume method is used for the discretization of the radiative transfer equation and the gray gas and weighted-sum-of-gray-gases (WSGG) models are implemented as a way to omit and consider the spectral dependence of the radiative properties, respectively. The TRI magnitude in the problem is evaluated by differences between values for the time-averaged heat fluxes at the wall (convective and radiative) and for the time-averaged radiative heat source calculated accounting for and neglecting the turbulence-radiation interaction effects. In general, TRI had little importance over all the considered cases, a conclusion that agrees with results of previous studies. When using the WSGG model, the contributions of the phenomenon were greater that with the gray gas hypothesis, demonstrating that the inclusion of the spectral variance in the solution of the radiative problem has an impact in the TRI effects. Furthermore, this paper presents a discussion, partly unprecedented in the context of the turbulence-radiation interaction, about the different methodologies that can be used for the TRI analysis. Finally, a correction factor is proposed for the time-averaged radiative heat source in the WSGG model, which is then validated by its implementation in the simulated cases. In future studies, a sensibility analysis on the terms that compose this factor can lead to a better understanding of how fluctuations of temperature correlate with the turbulence-radiation interaction phenomenon.

Neste trabalho, é apresentada a simulação de escoamentos compressíveis turbulentos no entorno de corpos móveis rígidos ou deformáveis empregando-se técnicas adaptativas. As simulações numéricas são conduzidas utilizando-se o método dos elementos finitos. A discretização espaço-temporal é desenvolvida através do método das linhas ou direções características (Characteristic-Based Split - CBS) e a modelagem da turbulência é feita através de um modelo de Simulação de Grandes Escalas (SGE, ou na terminologia em inglês, Large Eddy Simulation – LES) com o coeficiente de Smagorinsky variável no tempo e espaço (SGE ou LES dinâmico). A análise estrutural de corpos deformáveis imersos no fluido é realizada através de um modelo de elementos finitos triangulares para análise de placas e cascas com não linearidade geométrica, usando materiais elásticos com comportamento linear. Conjuntamente, um método de adaptação anisotrópica transiente de malhas é empregado para obter resultados com boa resolução a baixos custos computacionais. A consideração do movimento relativo de corpos imersos no escoamento é feita através de um método híbrido de movimento da malha que emprega interpolação com funções de base radial. Exemplos bidimensionais e tridimensionais são apresentados de forma a validar cada uma das metodologias desenvolvidas. Por fim, exemplos de simulações complexas são investigados, comparando-se os resultados obtidos com resultados experimentais e numéricos presentes na literatura. / In this work, the simulation of compressible turbulent flows around rigid and flexible moving bodies is presented using adaptative techniques. The numerical simulations are solved employing the finite element method. The space-time discretization is performed using the Characteristic-Based Split scheme (CBS) and turbulence is modelled with Large Eddy Simulation (LES) and a dynamic Smagorinsky sub-grid model. The structural analysis of deformable bodies immersed on the flow is performed using a triangular finite element model for the analysis of geometrically non-linear elastic plates and shells. An anisotropic mesh adaptation algorithm for transient simulations is coupled with the solver to achieve results with good resolution and low computational costs. The consideration of the relative movement of immersed bodies on the flow is performed employing an hybrid method of mesh movement based on radial basis function interpolation. Twodimensional and three-dimensional examples are presented in order to validate the proposed methodologies. Finally, complex simulations are investigated, where results are compared with experimental and numerical data available in the literature.

Submitted by Edisangela Bastos (edisangela@ufpa.br) on 2013-08-22T12:06:12Z No. of bitstreams: 2 Dissertacao_EstudoInfluenciaGranulometria.pdf: 4489674 bytes, checksum: dee436227a77f4053e8ba002d10daad7 (MD5) license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) / Approved for entry into archive by Ana Rosa Silva(arosa@ufpa.br) on 2013-08-23T14:32:24Z (GMT) No. of bitstreams: 2 Dissertacao_EstudoInfluenciaGranulometria.pdf: 4489674 bytes, checksum: dee436227a77f4053e8ba002d10daad7 (MD5) license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) / Made available in DSpace on 2013-08-23T14:32:24Z (GMT). No. of bitstreams: 2 Dissertacao_EstudoInfluenciaGranulometria.pdf: 4489674 bytes, checksum: dee436227a77f4053e8ba002d10daad7 (MD5) license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Previous issue date: 2011 / O presente trabalho apresenta estudos de caracterização reológica e determinações de fatores de atrito em três polpas de bauxita, originárias da Mina de Miltônia, Paragominas – Pará, resultantes de condições operacionais com graus diferentes de moagem. Deste modo, a diferenciação básica entre as mesmas reside na distribuição granulométrica. O objetivo é fazer uma proposta de revisão da especificação granulométrica do produto (polpa de bauxita para o mineroduto), permitindo uma comparação com as polpas do projeto (planta piloto) e da operação atual (usina de beneficiamento) quanto a resposta a alterações, principalmente, da viscosidade e fatores de atrito. A polpa proposta incorpora um percentual maior de finos (< 10 microns) em relação as demais polpas. Os benefícios possibilitarão ganhos de recuperação com a incorporação desta parcela de finos, a qual atualmente é descartada para a bacia de rejeitos. O material foi caracterizado por análises granulométricas da série Tyler, físicas e físico-químicas, microscopia eletrônica de varredura (MEV) e espectrometria de raios-x dispersiva em energia (EDS). Utilizando-se de um viscosímetro de cilindros coaxiais, tipo Searle, abordo use os aspectos reológicos destas polpas, em função da sua distribuição granulométrica e da variação da concentração de sólidos. As propriedades reológicas avaliadas foram a viscosidade, a tensão de escoamento e também os comportamentos que mostram a dependência da viscosidade com o tempo de aplicação de uma taxa de cisalhamento constante, verificando-se por meio de curva de histerese, a existência ou não da tixotropia e da reopexia. Os resultados experimentais mostraram que para as polpas ensaiadas, a granulometria e a concentração de sólidos são parâmetros relevantes na determinação da viscosidade. Os testes com estas três polpas de bauxita, de granulometrias diferentes, foram realizados a concentrações pré-determinadas, de modo a determinar o modelo reológico cujos parâmetros apresentaram os melhores coeficientes de correlação (R2), sendo que os melhores ajustes encontrados foram segundo o modelo de Herschel-Bulkley. Os fatores de atrito foram obtidos por formulações que envolvem uma relação logarítmica, utilizando-se da Lei da Parede como método de obtenção das expressões aplicáveis a estes fluidos não- Newtonianos viscoplásticos, sendo comparados com os calculados pelas correlações de Dodge e Metzner e verificou-se que de uma maneira geral há a concordância entre o modelo estudado nas três polpas. A análise dos fatores de atrito na concentração de 50% de sólidos, concentração de bombeio por mineroduto, mostra que o comportamento das polpas estudadas é semelhante, com os melhores ajustes para a faixa de 10.000 a 100.000 Reynolds (faixa média de trabalho para o escoamento turbulento de polpas de minério, homogêneas e heterogêneas). / This paper presents studies of the rheological characterization and determination of friction factors in three bauxite slurries originating in Miltonia Mine, Paragominas - Para, resulting from operating conditions with different degrees of grinding. Thus, the basic difference between them lies in the size distribution. The goal is to make a proposal to revise the specification of the product particle size (for the bauxite slurry pipeline), allowing a comparison with the slurries of the project (pilot plant) and the current operation (processing plant) and the response to changes, especially, viscosity and friction factors. The slurry proposal incorporates a higher percentage of fines (<10 microns) compared with other slurries. The benefits can earn recovery with the incorporation of this portion of fines, which currently is discarded to the tailings pond. The material was characterized by size analysis of the series Tyler, physical and physical-chemical, electron microscopy (SEM) and X-ray spectrometry energy dispersive (EDS). Testes were conducted in coaxial cylinder viscometer, Searle type, in order to determine the rheological behavior of these slurries, depending on the size distribution and varying the concentration of solids. The rheological properties were measured viscosity, yield stress and also behaviors that show the dependence of viscosity with the time of application of a constant shear rate, verifying through the hysteresis curve, the presence or absence of thixotropy and reopexia. The experimental results showed that for the slurries tested, the particle size and solids concentration are important parameters in determining the viscosity. Tests with these three bauxite slurries of different sizes, were performed at predetermined concentrations, to determine the model whose rheological parameters showed the best correlation coefficients (R2), and the best fits were found using the Herschel-Bulkley model. The friction factors were obtained for formulations involving a logarithmic relationship, using the Law of the Wall as a method of obtaining these expressions apply to non-Newtonian viscoplastic fluid, and compared with those calculated by Dodge and Metzner correlation and found that in general there is agreement between the model studied in the three slurries. The factor analysis of friction in the concentration of 50% solids, which is the concentration of pumping by pipeline, shows that the behavior of the slurries studied is similar, with the best settings for the range from 10,000 to 100,000 Reynolds (average work range for the turbulent flow of the ore pulp, homogeneous and heterogeneous).

Made available in DSpace on 2016-12-12T20:25:12Z (GMT). No. of bitstreams: 1 Diego da Silva.pdf: 2316377 bytes, checksum: 8cdc8aef19b579ce777e88c236323f74 (MD5) Previous issue date: 2015-12-10 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In practical applications related to the drying process is absolutely essential to understand the basic echanism of mass transfer of chemical species by convection. Once understood this physical phenomenon, it is possible to propose techniques to increase the convective mass transfer coefficient of the chemical species in order to enhance or control the drying. Within this context, the main objective of this study was to develop, validate and apply a numerical simulation tool for the study of the turbulent flow of a binary mixture (air-water vapor) inside the channels. The numerical method was developed using finite difference equations being discretized by central differences in co-located mesh. Artificial dissipation terms have been introduced to control the roblem of odd-even decoupling and minimize numerical instabilities. Because the flow is turbulent, the turbulence model type high Reynolds_ &#603;, since the industrial applicability of the model was used. The physical study was to evaluate drying important parameters, namely: (i) determining the rates of mass transfer of chemical species on the surface of the channels (solid / gas interface) in problems of forced convection; (ii) analyzing the effects of parameters / important characteristics of the flow of the mixture over a drying process. The results established that the mass transfer of the chemical species to the flow channel wall is directly affected by parameters considered, particularly in the entrance region and the channel region resulting in recirculation. / Em aplicações práticas relacionadas ao processo de secagem é absolutamente essencial compreender o mecanismo básico da transferência de massa da espécie química por convecção. Uma vez compreendido este fenômeno físico, é possível propor técnicas para aumentar o coeficiente convectivo de transferência de massa de forma a intensificar ou controlar a secagem. Dentro deste contexto, o principal objetivo deste trabalho foi desenvolver, validar e aplicar uma ferramenta de simulação numérica para o estudo do escoamento turbulento de uma mistura binária (ar-vapor d água) no interior de canais. O método numérico foi desenvolvido em diferenças finitas, sendo as equações de governo discretizadas em diferenças centrais em malha co-localizada. Foram introduzidos termos de dissipação artificial para controlar o problema do desacoplamento par-ímpar e minimizar instabilidades numéricas. Devido o escoamento ser turbulento, foi utilizado o modelo de turbulência do tipo k &#8722; &#603; alto Reynolds, visto a aplicabilidade industrial do modelo. O estudo físico realizado consistiu em avaliar parâmetros importantes da secagem, a saber: (i) determinação das taxas de transferência de massa da espécie química na superfície dos canais (interface sólido/gás) em problemas de convecção forçada; (ii) análise dos efeitos de parâmetros/características importantes do escoamento da mistura sobre um processo de secagem. Os resultados obtidos estabeleceram que a transferência de massa da espécie química das paredes do canal é diretamente afetada pelos parâmetros considerados, principalmente na região de entrada do canal e em regiões que resultam em recirculação.

Este trabalho apresenta a simulação numérica do escoamento turbulento em torno de uma turbina eólica de eixo vertical de pequeno porte, Savonius tipo helicoidal com torção de 180° nas pás. Com o intuito de avaliar a metodologia computacional empregada os resultados numéricos obtidos são comparados com os resultados experimental e numérico contidos no estado da arte. Também, compara-se o coeficiente de toque da turbina Savonius helicoidal com a turbina Savonius convencional. As simulações numéricas são baseadas no Método de Volumes Finitos, e para tal emprega-se o programa Fluent /Ansys versão 13.0 que resolve as equações da continuidade e as equações de Navier-Stokes com médias de Reynolds, juntamente com o modelo de turbulência . As simulações são desenvolvidas empregando diferentes malhas computacionais em estudos transientes, tridimensionais, com a turbina estacionária. A avaliação da qualidade da malha é realizada através do método de Índice de Convergência de Malha (GCI) o qual analisa o quão longe os resultados estão da solução assintótica para a malha utilizada. Após a análise da qualidade de malha, realizam-se simulações com a turbina em rotação as quais fazem uso da malha contendo uma região móvel possibilitando a imposição de uma velocidade angular ao rotor. O coeficiente de torque é obtido nas simulações e a partir dele calcula-se o coeficiente de potência. Além da análise do desempenho do rotor realiza-se uma análise qualitativa das características do escoamento sobre a turbina. A turbina Savonius helicoidal apresenta um valor de coeficiente de potência de 0,175 para a razão de velocidade de ponta de 0,58 considerando correção do efeito de bloqueio. Os resultados obtidos apresentam boa concordância com os resultados publicados por outros autores. / This dissertation presents the numerical simulation of the turbulent flow around of a small sized vertical axis wind turbine, consisting in a helical Savonius type with a 180° degree of blade twist. In order to evaluate the used methodology the obtained results are compared with the state of the art numerical and experimental data. It will be also presented the comparison between the torque coefficient of the conventional Savonius turbine and the helical Savonius turbine. The numerical simulations are based on the Finite Volume Method (FVM), using the commercial code Fluent/ANSYS version 13.0, which solves the continuity and Navier-Stokes through the Reynolds time-averaged methodology, including the turbulence model. The simulations are developed using different computational meshes for transient and three-dimensional studies with the stationary turbine. The evaluating the quality of the mesh is performed by of Grid Convergence Index (GCI) method which analyzes how far the results are the asymptotic solution to the mesh used. After the evaluation of the mesh quality, it was simulated a case considering the rotor motion using the moving mesh configuration, allowing the imposition of an angular velocity to the turbine. In the post-processing stage, it is possible to obtain the torque coefficient on the rotor shaft, allowing the calculation of the power coefficient for the turbine. In addition to the performance analysis, it is also made a qualitative analysis of the flow characteristics over the turbine rotor and in both cases presenting a good correspondence with the results in the literature. The helical Savonius turbine presents a value of power coefficient of 0.175 to a tip speed ratio of 0.58 whereas blocking effect correction.

A interação turbulência-radiação (TRI, do inglês Turbulence-Radiation Interaction) resulta do acoplamento altamente não linear entre flutuações da intensidade de radiação e flutuações da temperatura e da composição química do meio, e tem-se demonstrado experimentalmente, teoricamente e numericamente que este é um fenômeno relevante em diversas aplicações envolvendo altas temperaturas, especialmente em problemas reativos. Neste trabalho, o TRI é analisado em um escoamento interno não reativo de um gás participante que se desenvolve em um duto de seção transversal quadrada, para diferentes intensidades de turbulência do escoamento e considerando duas espécies distintas para a composição do fluido de trabalho (dióxido de carbono e vapor de água). O objetivo central é avaliar como a inclusão ou não da variação espectral das propriedades radiativas do meio no cálculo influencia a magnitude do TRI. Isso é feito através de simulações numéricas no código de dinâmica dos fluidos computacional Fire Dynamics Simulator (FDS), que resolve, através do método dos volumes finitos, as equações fundamentais que regem o problema – isto é, os balanços de massa, de quantidade de movimento e de energia e a equação de estado – em uma formulação adequada para baixos números de Mach, utilizando um algoritmo de solução explícito e de segunda ordem no tempo e no espaço. A turbulência é modelada através da simulação de grandes escalas (LES, do inglês Large Eddy Simulation), empregando-se o modelo de Smagorinsky dinâmico para o fechamento dos termos submalha; para a radiação térmica, o método dos volumes finitos é utilizado na discretização da equação da transferência radiativa e os modelos do gás cinza e da soma-ponderada-de-gases-cinza (WSGG, do inglês Weighted-Sum-of-Gray-Gases) são implementados como forma de desconsiderar e de incluir a dependência espectral das propriedades radiativas, respectivamente. A magnitude do TRI sobre o problema é avaliada através de diferenças entre as médias temporais dos fluxos de calor superficiais e do termo fonte radiativo obtidas em cálculos que consideram os efeitos do fenômeno e cálculos que os negligenciam. Em geral, a interação turbulência-radiação mostrou ser pouco importante em todos os casos considerados, o que concorda com resultados de outros estudos sobre o tema em escoamento não reativos. Com o modelo WSGG, as contribuições do fenômeno foram maiores do que com a hipótese do gás cinza, evidenciando que a inclusão da variação espectral na solução do problema radiativo tem um impacto sobre a magnitude dos efeitos do TRI. Além disso, é feita uma discussão, em parte inédita no contexto do TRI, sobre diferentes metodologias para a análise do fenômeno. Finalmente, é proposto um fator de correção para o termo fonte radiativo médio no modelo WSGG, que é validado através de sua implementação nos casos simulados. Em estudos futuros, uma análise de sensibilidade sobre os termos constituintes desse fator de correção pode levar a um melhor entendimento de como as flutuações de temperatura se correlacionam com o fenômeno da interação turbulência-radiação. / Turbulence-radiation interaction (TRI) results from the highly non-linear coupling between fluctuations of radiation intensity and fluctuations of temperature and chemical composition of the medium, and its relevance in a number of high-temperature problems, especially when chemical reactions are included, has been demonstrated experimentally, theoretically, and numerically. In the present study, the TRI is analyzed in a channel flow of a non-reactive participating gas for different turbulence intensities of the flow at the inlet and considering two distinct species for the medium composition (carbon dioxide and water vapor). The central objective is to evaluate how the inclusion or not of the spectral variation of the radiative properties of a participating gas in the radiative transfer calculations affects the turbulence-radiation interaction. With this purpose, numerical simulations are performed using the computational fluid dynamics Fortranbased code Fire Dynamics Simulator, that employs the finite volume method to solve a form of the fundamental equations – i.e., the mass, momentum and energy balances and the state equation – appropriate for low Mach number flows, through an explicit second-order (both in time and in space) core algorithm. Turbulence is modeled by the large eddy simulation approach (LES), using the dynamic Smagorinsky model to close the subgrid-scale terms; for the thermal radiation part of the problem, the finite volume method is used for the discretization of the radiative transfer equation and the gray gas and weighted-sum-of-gray-gases (WSGG) models are implemented as a way to omit and consider the spectral dependence of the radiative properties, respectively. The TRI magnitude in the problem is evaluated by differences between values for the time-averaged heat fluxes at the wall (convective and radiative) and for the time-averaged radiative heat source calculated accounting for and neglecting the turbulence-radiation interaction effects. In general, TRI had little importance over all the considered cases, a conclusion that agrees with results of previous studies. When using the WSGG model, the contributions of the phenomenon were greater that with the gray gas hypothesis, demonstrating that the inclusion of the spectral variance in the solution of the radiative problem has an impact in the TRI effects. Furthermore, this paper presents a discussion, partly unprecedented in the context of the turbulence-radiation interaction, about the different methodologies that can be used for the TRI analysis. Finally, a correction factor is proposed for the time-averaged radiative heat source in the WSGG model, which is then validated by its implementation in the simulated cases. In future studies, a sensibility analysis on the terms that compose this factor can lead to a better understanding of how fluctuations of temperature correlate with the turbulence-radiation interaction phenomenon.

Neste trabalho, é apresentada a simulação de escoamentos compressíveis turbulentos no entorno de corpos móveis rígidos ou deformáveis empregando-se técnicas adaptativas. As simulações numéricas são conduzidas utilizando-se o método dos elementos finitos. A discretização espaço-temporal é desenvolvida através do método das linhas ou direções características (Characteristic-Based Split - CBS) e a modelagem da turbulência é feita através de um modelo de Simulação de Grandes Escalas (SGE, ou na terminologia em inglês, Large Eddy Simulation – LES) com o coeficiente de Smagorinsky variável no tempo e espaço (SGE ou LES dinâmico). A análise estrutural de corpos deformáveis imersos no fluido é realizada através de um modelo de elementos finitos triangulares para análise de placas e cascas com não linearidade geométrica, usando materiais elásticos com comportamento linear. Conjuntamente, um método de adaptação anisotrópica transiente de malhas é empregado para obter resultados com boa resolução a baixos custos computacionais. A consideração do movimento relativo de corpos imersos no escoamento é feita através de um método híbrido de movimento da malha que emprega interpolação com funções de base radial. Exemplos bidimensionais e tridimensionais são apresentados de forma a validar cada uma das metodologias desenvolvidas. Por fim, exemplos de simulações complexas são investigados, comparando-se os resultados obtidos com resultados experimentais e numéricos presentes na literatura. / In this work, the simulation of compressible turbulent flows around rigid and flexible moving bodies is presented using adaptative techniques. The numerical simulations are solved employing the finite element method. The space-time discretization is performed using the Characteristic-Based Split scheme (CBS) and turbulence is modelled with Large Eddy Simulation (LES) and a dynamic Smagorinsky sub-grid model. The structural analysis of deformable bodies immersed on the flow is performed using a triangular finite element model for the analysis of geometrically non-linear elastic plates and shells. An anisotropic mesh adaptation algorithm for transient simulations is coupled with the solver to achieve results with good resolution and low computational costs. The consideration of the relative movement of immersed bodies on the flow is performed employing an hybrid method of mesh movement based on radial basis function interpolation. Twodimensional and three-dimensional examples are presented in order to validate the proposed methodologies. Finally, complex simulations are investigated, where results are compared with experimental and numerical data available in the literature.