No presente trabalho propõe-se empregar algoritmos genéticos em associação com o design construtal para a otimização de geometrias em problemas de transferência de calor. O objetivo principal de todos os estudos deste trabalho é minimizar a máxima temperatura que ocorre no domínio computacional. Investigou-se, inicialmente, uma cavidade isotérmica em forma de Y inserida em um sólido retangular com geração de calor uniforme a uma taxa volumétrica constante, onde foi feita uma comparação e validação do algoritmo genético frente à busca exaustiva para poucos graus de liberdade. Após, foi feita uma otimização usando somente algoritmos genéticos considerando todos os quatro graus de liberdade do problema e diferentes valores para suas restrições geométricas. O estudo seguinte foi feito considerando a mesma geometria anteriormente discutida, porém considerou-se as paredes da cavidade Y com uma condição de contorno convectiva. Da mesma forma anterior, foi feita uma validação do algoritmo genético frente à busca exaustiva e na sequência uma otimização de todos os quatro graus de liberdade e diferentes valores do parâmetro convectivo a, empregando somente algoritmos genéticos. No terceiro caso, estudou-se um caminho assimétrico em forma de V de um material de alta condutividade. A geometria tem sua base recebendo um fluxo de calor constante e o remove através das extremidades de dois braços ligados a um sumidouro de calor. Otimizou-se a forma pelo método exaustivo considerando quatro graus de liberdade e uma restrição constante . Após, usou-se algoritmos genéticos para otimizar a geometria considerando os mesmos graus de liberdade e diferentes valores para a restrição de ocupação do material condutivo. Similarmente ao caso da cavidade convectiva em forma de Y, por fim, estudou-se a otimização geométrica de um corpo cilíndrico onde cavidades convectivas retangulares com dois pares de braços são inseridas. Realizaram-se otimizações de até sete graus de liberdade e também se estudou a influência de um parâmetro convectivo e das frações de ocupação das áreas do corpo e braços da cavidade. Deste estudo, concluiu-se que quanto maior o número de cavidades, menores são as máximas temperaturas que ocorrem no domínio. Destaca-se, também, a dependência do parâmetro convectivo, que influenciou na forma da melhor geometria encontrada. Para todos os estudos feitos, os resultados mostraram que a busca por meio de algoritmos genéticos levou a uma redução significativa do número de simulações necessárias para obter a geometria ótima com resultados concordantes aos obtidos com busca exaustiva. Além disso, foi possível estender o estudo para problemas com mais graus de liberdade, restrições e propriedades térmicas. Conclui-se que o melhor design é altamente dependente dos graus de liberdade e restrições, este sendo alcançado de acordo com o princípio construtal da ótima distribuição das imperfeições. / In this work, we propose employing genetic algorithms in association with constructal design for geometry optimization in heat transfer problems. The main objective of all studies is to minimize the maximum temperature that occurs in the computational domain. It was investigated initially an isothermal Y-shaped cavity intruded into a rectangular solid conducting wall with heat generation uniformly at a volumetric rate, where a comparison and validation of genetic algorithm against exhaustive search for few degrees of freedom was made. Then, an optimization is performed by means of genetic algorithms considering all four degrees of freedom of the problem and different values for geometric constraints. The following study has been done considering the same geometry as previously discussed, but it is considered the walls of the Y-cavity with a convective boundary condition. Thus, a dimensionless heat transfer parameter to study (a) was added. Similarly, foregoing study, a genetic algorithm validation was performed comparing to the exhaustive search. After, all four degrees of freedom and different values of a parameter only using genetic algorithms were optimized. In the next investigation, an asymmetric V-shaped pathway of high conductivity material was studied. This geometry receives a constant heat transfer rate in its base and removes it by the end of the two branches that are in touch with the heat sink. The shape was optimized by exhaustive approach considering four degree of freedom and a constraint. After, we used genetic algorithms to optimize the geometry considering the same degrees of freedom and different values for the restriction. Finally, similar to the case of the Y-shaped convective cavity, rectangular convective cavities with two pairs of arms inserted into a cylindrical solid body were optimized. Optimizations of up to seven degrees of freedom were performed and the influence of the convective parameter and of the area fractions of the body and arms of the cavity, were also investigated. From this study, it was concluded that the higher the number of cavities, the lower the maximum temperatures occurring in the domain. Also, the dependence of the convective parameter, influenced in the form of the best geometry, is highlighted. For all studies carried out, the results showed that the search using genetic algorithms led to a significant reduction of the number of simulations required to obtain the optimal geometry. Moreover, it was possible to extend the study where it was considered other degrees of freedom, constraints and thermal properties. We concluded that the best design is highly dependent of degrees of freedom and constraints, and this has been achieved according to the constructal principle of optimal distribution of imperfections.

O objetivo deste trabalho consiste no desenvolvimento e aplicação de um software de simulação para o processo de lingotamento contínuo que seja capaz de simular com precisão e confiabilidade a solidificação de tarugos de aço, possibilitando a otimização do processo de fabricação dos produtos siderúrgicos. Para a validação do software, utilizou-se de resultados obtidos experimentalmente de medições de temperaturas em planta e de amostras de aços lingotados. As medidas de temperatura foram realizadas com a inserção de termopares ao longo da parede de moldes e com uso de pirômetro óptico posicionado em diferentes regiões da máquina de LC. O modelo desenvolvido mostrou-se bastante coerente e preciso quando seus resultados foram comparados com os dados experimentais do processo. Foram analisados lingotes de seção quadrada de 240 mm, 180 mm e 150 mm. Com a simulação de uma seqüência de corridas monitoradas foi possível obter resultados para ajuste de parâmetros térmicos como os coeficientes de transferência de calor nas interfaces metal/molde e metal/ambiente. Posteriormente, aferições foram realizadas com dados reais de processo e com outros programas numéricos da literatura aplicados ao processo de LC. Com o programa de simulação devidamente aferido, correlações numéricas e experimentais foram realizadas para verificar a influência de parâmetros térmicos (ρ, c, k, TL, TS), parâmetros de processo (TV, VL, CQ) e parâmetros estruturais (λ1, λ2, TCE) na evolução da casca solidificada e perfil térmico dos tarugos. Dessa forma, equações empíricas que correlacionam esses parâmetros foram analisadas e posteriormente adotadas no programa numérico. / The development and application of simulation software for the continuous casting process is proposed in this work. The heat transfer and solidification models were validated with experimental results of temperature measurements in an industrial plant as well as with steel billets samples. Thermocouples embedded in the mold walls were used to measure temperatures along the mold length and non-contact pyrometers positioned at different locations along the machine were used to monitored the strand surface temperatures. The model seemed to predict quite accurately the continuous casting simulation when calculated results were compared with experimental measurements. In this study, steel billet with square section of 240 mm, 180 mm and 150 mm were analyzed. Some of the monitored heats were used to adjust thermal parameters as the heat transfer coefficients. After, the model results were compared with results from a heat transfer model of literature. With the simulation software properly calibrated, numerical and experimental correlations were performed to analyze the influence of thermal parameters (ρ, c, k, TL, TS), process parameters (TV, VL, CQ) and structural parameters (λ1, λ2, TCE) in the evolution of solidified shell and thermal profile of the billets. Therefore, empirical equations that correlate these parameters were analyzed and subsequently adopted in the numerical model.

Ce travail porte sur l’étude du transfert de chaleur entre un milieu poreux et un milieu libre en utilisant une approche multi-échelle. Dans un premier temps, nous dérivons les conditions limites à imposer à une interface libre-poreux dans le cas des transferts de chaleur à l’équilibre thermique local puis dans le cas du déséquilibre thermique local. Pour les transferts de chaleur turbulent, une simulation numérique directe est réalisée afin d’apporter une meilleur compréhension de la physique existant à l’interface libre poreux. Puis, nous déterminons un modèle turbulent avec des conditions de saut. Ces études répondent à des questions fondamentales liées à la compréhension physique des conditions de saut, des valeurs des paramètres des sauts et de la position de l’interface dans le cadre des transferts de chaleur. / This work deals with the study of heat transfer between a porous medium and a free medium, using multi scale approaches. First, we derive the boundary conditions that must be applied at a free-porous interface for laminar heat transfer at local thermal equilibrium and, then, at local thermal non-equilibrium. For turbulent heat transfer, a direct numerical simulation is performed supplying a better understanding of the physic at the free-porous interface. Then, we determine a turbulent model with associated jump conditions. These studies answer fundamental questions regarding the physical meaning of the jump conditions, the values of the jump parameters and the location of the interface for heat transfer.

Cette thèse aborde la question de la modélisation des transferts thermiques dans les milieux poreux soumis à de hautes températures. Un modèle macroscopique hors équilibre thermique local entre phases est obtenu par changement d’échelle. Cette procédure tient compte à l’échelle locale du couplage entre rayonnement et autres modes de transfert. Le modèle de rayonnement repose sur l’équation de transfert radiatif généralisée (GRTE) et, à la limite courante d’un milieu macroscopiquement optiquement épais, sur la loi de Fourier radiative. L’originalité de cette procédure réside dans l’application de la méthode de prise de moyenne volumique (VAM) aux équations de bilan local, dans lesquelles les transferts radiatifs sont inclus. Cette homogénéisation couplée soulève trois difficultés : - les différents transferts sont de natures différentes. Le système matériel (site de transferts conductifs et convectifs) coexiste avec le champ des photons qui est homogénéisé par une méthode statistique reposant sur la caractérisation des propriétés radiatives au moyen de fonctions de distribution continûment définies sur le volume du milieu poreux. - les échelles de longueur mises en jeu dans la procédure de changement d’échelle doivent être compatibles entre elles. On établit que la séparation des échelles, requise par la prise de moyenne volumique, est compatible avec l’échelle de longueur caractéristique de l’homogénéisation statistique radiative, seulement limitée par la résolution d’une tomographie du milieu. - le phénomène d’émission dépend de la température de la matière. Cette température spécifique au calcul radiatif est obtenue en appliquant un opérateur de prise de moyenne à la température de la matière sur une échelle locale représentative. En pratique, c’est la résolution de cette prise de moyenne qui définit l’échelle des couplages locaux avec la méthode VAM. Le modèle macroscopique résultant est appliqué à la résolution d’un problème unidimensionnel et stationnaire. Dans ce cas simple, le rôle du couplage avec le rayonnement à l’échelle locale est mis en évidence. / This thesis deals with the modeling of heat transfer in porous media subjected to high temperature. An upscaling procedure yields a macroscopic model based on local thermal non-equilibrium. This procedure accounts for local scale coupling effects between radiation and other transfer modes. Radiation modeling uses the generalized radiation transfer equation (GRTE) and, at the commonly encountered limit of a macroscopically optically thick medium, the radiative Fourier law. An original feature of this procedure is the application of the volume averaging method (VAM) to local conservation equations in which radiation transfer is included. This raises three major challenges: - the physical natures of the various transfer modes involved are different. The material system (in which conduction and convection occur) coexists alongside with the photon field, which is homogenized using a statistical method based on the characterization of the radiative properties through statistical distribution functions, continuously defined over the whole volume of the porous medium. - the length scales involved in the upscaling procedure must be compatible with each other. The compatibility of the scale separation constraint (required by the VAM) with the length scale of the radiative homogenization technique (which is limited by the resolution of a tomography of the medium) is established. - the emission phenomenon depends on the temperature of matter. This temperature, specific to the radiation calculation, is obtained by applying a dedicated averaging operator. This operator is associated with an averaging volume whose length scale has to be representative at the local scale. In practice, the resolution of this averaging procedure defines the scale of the coupling between the VAM and the radiation model. The resulting macroscopic model is applied to a one-dimensional, steady case. The solving of this simple case shows the influence of local coupling effects.

Les transferts radiatifs jouent un rôle important dans les chambres de combustion des installations industrielles. En effet, il existe un couplage fort entre la combustion turbulente et le rayonnement. Dans le but d’étudier ce couplage, le code Rainier est développé pour les calculs de pertes par rayonnement dans un écoulement réactif dans des géométries complexes. Ce code repose sur des simulations aux grandes échelles (LES) de la combustion turbulente. Il est basé sur les maillages tétraédriques non structurés. Le modèle de rayonnement appliqué à la modélisation des propriétés radiatives des gaz est le modèle CK (Correlated-k). La méthode statistique de Monte-Carlo (ERM) est utilisée pour résoudre l’équation de Transfert du Rayonnement (ETR). Le code de rayonnement est parallélisé et il montre une réponse linéaire en fonction du nombre de processeurs très proche de la réponse idéale. Une méthode de couplage de code de rayonnement avec le code de combustion LES est développée. Chacun des codes a sa propre logique d’architecture et de développement. En conséquence, le couplage entre les deux domaines d’étude est réalisé de telle façon que les échanges des données et les synchronisations entre eux soient assurés. Les résultats obtenus à partir du couplage des sur une chambre de combustion d’hélicoptère sont présentés. Nous avons montré que le rayonnement modifie les champs instantanés de température et d’espèces à l’intérieur de la chambre de combustion. / Radiation plays an important role in industrial combustion chambers. In fact, there is a strong coupling between combustion and heat transfer in these turbulence chambers. The Rainier code is developed for the calculation of radiation in a reactive flow within a complex geometry. This code is dedicated to Large Eddy Simulation (LES) of turbulent reactive flows. The code is based on unstructured tetrahedral mesh. The correlated k-distribution method (CK) is applied to the modelling of radiative properties of gases. A statistical method of Monte Carlo (ERM) is used to resolve the Radiation Transfer Equation (RTE). The code is parallelized and it shows a linear response to the number of processors, which is very close to the ideal response. A coupling method of the Rainier code and the turbulent combustion code LES is implemented. Each code has its own logic architecture and development. For this, the coupling between these two different fields of study is achieved in such a way that the data exchange and synchronization between them is assured. The results obtained by applying the coupled codes to the combustion chamber of helicopter are presented. We have shown that the radiative transfer modifies the temperature and species fields inside the combustion chamber.

O trabalho consiste em uma análise experimental comparativa do processo de combustão e de transferência de calor ao longo de uma fornalha flamotubular. Os combustíveis utilizados são o óleo diesel e o biodiesel proveniente de óleo de fritura. Uma outra parcela deste trabalho é investigar os níveis de emissões gasosas no processo de combustão do biodiesel em fornalha flamotubular. Não seria possível realizar esta pesquisa sem uma fornalha flamotubular dotada de câmaras calorimétricas, janelas de observação e extração de gás, e totalmente instrumentada. Assim, foi efetuado em paralelo o projeto e a construção de uma fornalha calorimétrica flamotubular. Os parâmetros utilizados para se analisar os combustíveis na fornalha são: rendimento térmico, taxa de fluxo mássico, temperatura de chama e perda de carga. Observou-se que o óleo diesel apresentou maior taxa de transferência de calor para a maioria das seções expostas à chama. Na região onde o corpo da chama não está presente, a taxa de transferência de calor do biodiesel torna-se ligeiramente superior. Em relação às emissões gasosas observaram-se baixos níveis de poluentes / This work consists in a comparative experimental evaluation of the process of combustion and heat transfer along a flame tube furnace. The fuels used are diesel oil and ethyl ester from waste vegetable oil (biodiesel). Another part of this work investigates the percentage of pollutants and emissions of gases in the flame tube furnace. It would not be possible to carry out this research without a flame tube furnace endowed with calorimetric chambers, windows of observation and extraction of gases, and totally instrumented. Thus, the project and the construction of a flame tube calorimetric furnace were done parallel. The parameters used to analyze the fuels in the furnace are: thermal efficiency, mass flow rate, temperature of flame and draught loss. It was observed that diesel oil presented higher heat transfer rate in most parts of the furnace that was near the flame. In the region where the body of flame is not present, the heat transfer rate of biodiesel becomes a little higher. In relation to emission of gases, it was observed low indexes of pollutants

O processamento térmico contínuo usando trocadores de calor é uma forma muito comum de industrialização de alimentos líquidos. Atrelada à inativação de micro-organismos e enzimas pela alta temperatura, está a perda de qualidade do produto; portanto, o estudo e avaliação deste tipo de processo é fundamental para garantir a inocuidade e a melhora da qualidade do produto. Neste trabalho foi desenvolvida a modelagem matemática de um processo térmico contínuo em trocadores bitubulares para determinação da distribuição da temperatura média do produto e da letalidade sobre um atributo de segurança ou qualidade. No modelo, foram considerados os trocadores de aquecimento e de resfriamento, bem como o tubo de retenção. Admitiu-se regime permanente, escoamento pistonado e leva-se em conta a troca de energia com o ambiente. Para a aplicação do modelo, foi adotado um sistema em escala de laboratório com vazão de processamento entre 10 e 50 L/h, trocadores com quatro módulos de troca térmica e área total de 0,13 m² e sistema de indicação e aquisição de dados de temperatura. Como fluidos de trabalho foram usados: mistura de glicerina/água 80 % (Newtoniano) e solução 1 % de carboxi-metil-celulose (CMC, não-Newtoniano). Fluidos e utilidade foram água quente pressurizada e água gelada. Os parâmetros de troca térmica do modelo foram ajustados com sucesso através de ensaios experimentais. Foram ajustados: 1) coeficientes de convecção natural do ar sobre as seções de aquecimento e de resfriamento; 2) parâmetros da correlação de Nusselt x Reynolds para o coeficiente de convecção dos fluidos de trabalho; 3) coeficiente global de troca térmica com o ambiente no tubo de retenção. A validação do modelo ajustado foi realizada através da comparação das distribuições de temperatura experimental com a predição do modelo. A letalidade, levando em conta o tempo espacial no percurso do produto foi calculada e avaliada, indicando importante contribuição da seção de aquecimento e forte influência da elevação de temperatura na entrada do tubo de retenção para compensar as perdas para o ambiente. O modelo proposto mostra-se útil para a avaliação de processos térmicos em sistemas tubulares. / The continuous thermal processing using heat exchangers is a very usual form of industrialization of liquid foods. Linked to the inactivation of microorganisms and enzymes by high temperature is the loss of product quality; thus, the study and evaluation of this type of process is critical to ensure the safety and improve product quality. In this work, it was developed a mathematical model of a continuous process with tubular heat exchangers to determine the distribution of the average temperature of the product and the lethality considering a safety or quality attribute. In the model, are considered the exchangers for heating and cooling and the holding tube. It was assumed steady state plug-flow and it was taken into account the energy exchange with the ambient air. For the application of the model, it was adopted a laboratory system with processing flow rate between 10 and 50 L/h, exchangers with four heat transfer modules and total area of 0.13 m² and a temperature indication and acquisition system. The working fluids used were: a mixture of glycerol/water 80 % (Newtonian) and a 1 % solution of carboxymethylcellulose (CMC, non-Newtonian). Utility fluids were pressurized hot water and cold water. The heat transfer parameters of the model were adjusted successfully using the experimental data. Were adjusted: 1) the coefficients of natural convection of the air over the heating and cooling sections; 2) the parameters of the Nusselt x Reynolds correlation for the convective coefficient of the working fluid; 3) the overall heat exchange coefficient between the product and the ambient air in the holding tube. The validation of the fitted model was performed by comparing the experimental temperature distributions with the prediction from the model. Lethality, taking into account the space-time in the path of the product was calculated and evaluated, indicating a significant contribution of the heating section and a strong influence on the temperature rise in the entrance of the holding tube to compensate for losses to the surroundings. The model proves useful for the evaluation of thermal processes in tubular systems.

A demanda por produtos industrializados que apresentem máxima preservação de suas características naturais têm crescido e feito as indústrias alimentícias re-analisarem seus processos para atingirem essa necessidade do mercado. A abordagem convencional conservadora, utilizada para o dimensionamento do processamento térmico de alimentos, pode levar ao sobre-processamento especialmente no regime laminar, devido às significativas distribuições de temperatura e tempos de residência existentes. Um modelo matemático, composto por equações diferenciais de massa e energia foi elaborado, considerando o processamento térmico de um fluido não-newtoniano, sob regime laminar, escoando em um trocador de calor bitubular. No modelo, se levou em conta as dispersões efetivas de massa e energia associadas com o escoamento laminar não ideal, as trocas de calor com o ambiente, a letalidade que ocorre no aquecimento e resfriamento e o perfil de velocidade. O modelo foi testado por meio de simulações do estudo de caso do processamento térmico de suco de graviola (fluido pseudoplástico) considerando a destruição de bolores e leveduras. Objetivou-se nas simulações avaliar o efeito de distintas considerações do modelo nas variáveis consideradas. Os resultados indicaram que as etapas de aquecimento e resfriamento contribuíram de forma significativa na letalidade do processo, assim como as considerações quanto às dispersões de massa e energia. Como exemplo, as letalidades, considerando a abordagem convencional (tubo de retenção isotérmico com velocidade máxima) e o modelo completo (com todas as considerações) apresentaram valores de 1,46 e 5,74, respectivamente. A flexibilidade do modelo elaborado, assim como os tempos computacionais pequenos necessários para obter os resultados são as principais vantagens do uso do mesmo. Acredita-se que o modelo elaborado pode contribuir de forma importante para o correto dimensionamento e avaliação de processos térmicos em indústrias de alimentos, permitindo que a demanda dos consumidores seja atendida. / The demand for industrialized food with maximum retention of sensorial and nutritional attributes has grown and made the food industries rethink operational conditions to meet this market expectation. The classic conservative approach, used for the design of thermal food processing can lead to over-processing specially in laminar regime, due to the existing significant temperature and residence time distributions. A mathematical model, comprising differential equations for mass and heat transfer was elaborated, considering the thermal processing of a non-Newtonian liquid, under laminar flow in a double-pipe heat exchanger. In the model, it was taken into account the effective mass and energy dispersions associated with the non-ideal laminar flow, the heat exchange with the ambient, the contribution from heating and cooling sections in the lethality and the velocity profile. The model was tested through simulations of a study case of soursop juice processing (pseudoplastic fluid) regarding the destruction of yeast and molds. The objective of the simulations was to evaluate the effect of distinct model assumptions on the variables. The results indicated that the heating and cooling sections and the assumptions regarding the effective mass and energy dispersions had an important contribution to the processing lethality. As an example, the lethality, regarding the conventional approach (isothermal holding tube at the maximum velocity) and the complete model (with all the assumptions) were 1.46 and 5.74, respectively. The model flexibility and the small computational time needed for the results to the obtained are the main advantages of its use. It is expected that the developed model can be an important contribution to the correct design and evaluation of thermal processing in food industries, allowing the consumer demands to the reached.

O objetivo deste trabalho consiste no desenvolvimento e aplicação de um software de simulação para o processo de lingotamento contínuo que seja capaz de simular com precisão e confiabilidade a solidificação de tarugos de aço, possibilitando a otimização do processo de fabricação dos produtos siderúrgicos. Para a validação do software, utilizou-se de resultados obtidos experimentalmente de medições de temperaturas em planta e de amostras de aços lingotados. As medidas de temperatura foram realizadas com a inserção de termopares ao longo da parede de moldes e com uso de pirômetro óptico posicionado em diferentes regiões da máquina de LC. O modelo desenvolvido mostrou-se bastante coerente e preciso quando seus resultados foram comparados com os dados experimentais do processo. Foram analisados lingotes de seção quadrada de 240 mm, 180 mm e 150 mm. Com a simulação de uma seqüência de corridas monitoradas foi possível obter resultados para ajuste de parâmetros térmicos como os coeficientes de transferência de calor nas interfaces metal/molde e metal/ambiente. Posteriormente, aferições foram realizadas com dados reais de processo e com outros programas numéricos da literatura aplicados ao processo de LC. Com o programa de simulação devidamente aferido, correlações numéricas e experimentais foram realizadas para verificar a influência de parâmetros térmicos (ρ, c, k, TL, TS), parâmetros de processo (TV, VL, CQ) e parâmetros estruturais (λ1, λ2, TCE) na evolução da casca solidificada e perfil térmico dos tarugos. Dessa forma, equações empíricas que correlacionam esses parâmetros foram analisadas e posteriormente adotadas no programa numérico. / The development and application of simulation software for the continuous casting process is proposed in this work. The heat transfer and solidification models were validated with experimental results of temperature measurements in an industrial plant as well as with steel billets samples. Thermocouples embedded in the mold walls were used to measure temperatures along the mold length and non-contact pyrometers positioned at different locations along the machine were used to monitored the strand surface temperatures. The model seemed to predict quite accurately the continuous casting simulation when calculated results were compared with experimental measurements. In this study, steel billet with square section of 240 mm, 180 mm and 150 mm were analyzed. Some of the monitored heats were used to adjust thermal parameters as the heat transfer coefficients. After, the model results were compared with results from a heat transfer model of literature. With the simulation software properly calibrated, numerical and experimental correlations were performed to analyze the influence of thermal parameters (ρ, c, k, TL, TS), process parameters (TV, VL, CQ) and structural parameters (λ1, λ2, TCE) in the evolution of solidified shell and thermal profile of the billets. Therefore, empirical equations that correlate these parameters were analyzed and subsequently adopted in the numerical model.

Em processos de combustão, uma determinação precisa dos parâmetros envolvendo transferência de calor influencia diretamente os demais fenômenos envolvidos. Dentre os mecanismos de transferência de calor presentes na combustão a radiação térmica é predominante, mas sua correta determinação impõe uma elevada complexidade, principalmente quando se trata da solução de meios participantes. O cálculo envolve propriedades de absorção que variam com a temperatura e o comprimento de onda, sendo então necessária a utilização de modelos espectrais para obter bons resultados com um baixo tempo computacional. Para o cálculo da transferência radiante, existem diversos modelos espectrais, desde modelos de simples implementação, como, por exemplo, o GG (gás cinza) e o WSGG (soma-ponderada-dos-gases-cinza), até modelos com um grau elevado de detalhamento, como o SLW (soma-ponderada-dos-gases-cinza baseado em linhas espectrais) e o CW (número de onda cumulativa). Como os modelos com maior grau de detalhamento são de complexa implementação, alguns autores preferem empregar modelos simplistas, como o GG (gás cinza), apenas por questões de conveniência, mesmo em detrimento da qualidade dos resultados. Uma forma de executar o cálculo da radiação térmica sem simplificações é levar em conta as absorções em cada comprimento de onda, sendo esses cálculos denominados integração linha-por-linha (LBL), por executar o cálculo da transferência radiante em cada linha de absorção, o que gera resultados benchmark, podendo ser utilizados para avaliar os diversos modelos existentes. Este trabalho tem por objetivo verificar e sintetizar a aplicação dos modelos espectrais, em configurações envolvendo concentração e temperatura não uniformes, onde são realizados cálculos em um meio contendo CO2, H2O e fuligem. São avaliados os modelos GG, WSGG, SLW e CW. Dentre os modelos avaliados, o que apresenta os melhores resultados para as condições apresentadas é o modelo WSGG. De forma a aprimorar o modelo WSGG, uma nova implementação para a solução de misturas é apresentada, a qual apresenta correlações para o H2O e para o CO2 geradas individualmente, possibilitando misturas com qualquer razão de concentração, mostrando que o modelo apresenta bons resultados em diversas situações e é uma boa opção para a solução de problemas de combustão. / In combustion processes a good determination of the heat transfer parameters are of great importance because of its direct influence in the computation of the chemical reactions rate in the process and, consequently, in the formation of the combustion products. Among the processes of heat transfer in combustion, thermal radiation is predominant, and their determination can be a very complex task, especially with participating medium. The analysis involves absorption properties that vary with the temperature and wavelength, and therefore it is necessary to use spectral models to ensure good results with low computational time. There are several spectral models developed along the years, since the simplistic models such as the GG (gray gas) and WSGG (weighted-sum-of-gray-gases), to more advanced methods such as the SLW (spectral line weighted-sum-of-gray-gases) and CW (cumulative wavenumber). Due advanced models are in general a hard task to implement, the option is to use simplified models, for example the GG, even working with considerably errors. In order to quantify these solutions, for temperature and concentration conditions of the absorbing species, it is necessary to implement the radiation heat transfer taking into account the absorption at each wavelength through line-by-line (LBL) integration, being this solution the exact one, or, the benchmark solution, which it is used to evaluate the spectral models. In this study, the LBL integration is carried out to evaluate some of the existing models in a non-isothermal and inhomogeneous medium containing CO2, H2O and soot. The work involves the GG, WSGG, SLW and CW spectral models. For the presented cases, the best results occur with WSGG model. In order to improve the WSGG model a new implementation for the mixture solution is presented, which solves the correlations for H2O and CO2 generated individually, enabling mixtures containing any concentration ratio, showing the good agreement of the spectral model at any condition, being the WSGG a good option to solve combustion problems.