A estimativa acurada do fluxo de calor radiativo na região próxima à chama do tipo jato não pré-misturada é necessária para garantir a segurança de pessoas e equipamentos em caso de vazamentos ou processo de descarte na indústria de petróleo e gás. A simulação computacional dos fenômenos físicos envolvidos na transferência de calor e combustão do processo tem alto custo computacional. No presente trabalho é apresentado um estudo para o emprego do modelo de múltiplas fontes ponderadas para a estimativa do fluxo de calor radiativo no campo próximo à chama. O modelo matemático simplificado tem baixo custo computacional e consiste em representar a transferência de calor radiativa por fontes pontuais distribuídas no eixo central da chama. Cada fonte tem um peso proporcional à contribuição de cada região discretizada da chama na transferência de calor por radiação. Para determinar o peso de cada fonte foi utilizada a análise inversa pelo método da Otimização Extrema Generalizada, no qual o fluxo de calor é dado de entrada enquanto o peso de cada fonte é dado de saída. Como dado de entrada foi utilizado o fluxo de calor radiativo medido experimentalmente de um conjunto de 12 chamas, com potência entre 0,139 e 0,554 kW. A análise inversa foi utilizada para recuperar os pesos, e a fração radiante, que geram o fluxo de calor radiativo com maior compatibilidade com os dados experimentais em três abordagens. A primeira abordagem consiste em aplicar a análise inversa em cada chama separadamente, e depois de obtidos os pesos correlaciona-los com a potência da chama. Na aplicação do método em cada chama individualmente, o desvio máximo do resultado do modelo com os dados experimentais é de 5%. Em uma segunda abordagem, foi utilizada a análise inversa para obter diretamente os coeficientes de uma função entre os pesos do modelo e a potência da chama e seu comprimento estimado, o desvio máximo encontrado é de 18,6%. Na terceira abordagem, a análise inversa foi utilizada novamente para encontrar os coeficientes da função que correlaciona a potência da chama com parâmetros do modelo, e neste caso foi utilizado o comprimento experimental da chama para posicionar as fontes pontuais. Para o caso em que sete fontes foram posicionadas a 2,25 vezes o comprimento da chama medido experimentalmente, o desvio máximo observado foi de 8,6% e o desvio médio de 2,9%. / The accurate approximation of radiative heat flux on non-premixed flame in the region adjacent to the flame is required to guarantee the staff and machinery safety on oil and gas industry in cases of fuel leaking or disposal process. The physical phenomena involved on heat transfer and combustion has high computational cost in numerical simulation, in accidental leaking cases it is prohibitive. In the present study it is presented an application of Weighted Multi-Point Sources model to predict the radiative heat flux in the region adjacent to the flame. The simplified mathematical model has low computational cost and simulates the radiative heat transfer by punctual sources placed at the flame axis. Each source has the weight proportional to heat transfer contribution by each flame portion. To set each source's weight, it was used the inverse analysis by the Generalized Extremal Optimization. In inverse analysis the radiative heat fluxes are the input data while the weight of each source is the sought variable. As input data are used the radiative heat fluxes measured from a set of 12 flames, with power between 0.139 and 0.554 kW. The inverse analysis was used to recover the weights, and the fraction of heat radiated, which generate radiative heat fluxes with greater compatibility with the experimental data on three approaches. The first approach is to apply the inverse analysis in each flame separately, obtain the weights and then correlates them with the flame power. In method application for each flame, the maximum deviation between model outcome and experimental data was less than 5%. In a second approach, it was used the inverse analysis to directly obtain the coefficients of a function between the model's parameters and the flame’s power. When used just the flame power as model input parameter to distributes the sources and predict the heat flux, the maximum deviation is 18.6%. In the third approach, the inverse analysis was again used to find the coefficients of the function which correlates the flame power with model parameters, in this case was used the experimental flame length for positioning the point sources. For the case in which seven sources are positioned at 2.25 times the experimental length of the flame, the maximum deviation observed was 8.6% and the mean deviation of 2 9%.

A estimativa acurada do fluxo de calor radiativo na região próxima à chama do tipo jato não pré-misturada é necessária para garantir a segurança de pessoas e equipamentos em caso de vazamentos ou processo de descarte na indústria de petróleo e gás. A simulação computacional dos fenômenos físicos envolvidos na transferência de calor e combustão do processo tem alto custo computacional. No presente trabalho é apresentado um estudo para o emprego do modelo de múltiplas fontes ponderadas para a estimativa do fluxo de calor radiativo no campo próximo à chama. O modelo matemático simplificado tem baixo custo computacional e consiste em representar a transferência de calor radiativa por fontes pontuais distribuídas no eixo central da chama. Cada fonte tem um peso proporcional à contribuição de cada região discretizada da chama na transferência de calor por radiação. Para determinar o peso de cada fonte foi utilizada a análise inversa pelo método da Otimização Extrema Generalizada, no qual o fluxo de calor é dado de entrada enquanto o peso de cada fonte é dado de saída. Como dado de entrada foi utilizado o fluxo de calor radiativo medido experimentalmente de um conjunto de 12 chamas, com potência entre 0,139 e 0,554 kW. A análise inversa foi utilizada para recuperar os pesos, e a fração radiante, que geram o fluxo de calor radiativo com maior compatibilidade com os dados experimentais em três abordagens. A primeira abordagem consiste em aplicar a análise inversa em cada chama separadamente, e depois de obtidos os pesos correlaciona-los com a potência da chama. Na aplicação do método em cada chama individualmente, o desvio máximo do resultado do modelo com os dados experimentais é de 5%. Em uma segunda abordagem, foi utilizada a análise inversa para obter diretamente os coeficientes de uma função entre os pesos do modelo e a potência da chama e seu comprimento estimado, o desvio máximo encontrado é de 18,6%. Na terceira abordagem, a análise inversa foi utilizada novamente para encontrar os coeficientes da função que correlaciona a potência da chama com parâmetros do modelo, e neste caso foi utilizado o comprimento experimental da chama para posicionar as fontes pontuais. Para o caso em que sete fontes foram posicionadas a 2,25 vezes o comprimento da chama medido experimentalmente, o desvio máximo observado foi de 8,6% e o desvio médio de 2,9%. / The accurate approximation of radiative heat flux on non-premixed flame in the region adjacent to the flame is required to guarantee the staff and machinery safety on oil and gas industry in cases of fuel leaking or disposal process. The physical phenomena involved on heat transfer and combustion has high computational cost in numerical simulation, in accidental leaking cases it is prohibitive. In the present study it is presented an application of Weighted Multi-Point Sources model to predict the radiative heat flux in the region adjacent to the flame. The simplified mathematical model has low computational cost and simulates the radiative heat transfer by punctual sources placed at the flame axis. Each source has the weight proportional to heat transfer contribution by each flame portion. To set each source's weight, it was used the inverse analysis by the Generalized Extremal Optimization. In inverse analysis the radiative heat fluxes are the input data while the weight of each source is the sought variable. As input data are used the radiative heat fluxes measured from a set of 12 flames, with power between 0.139 and 0.554 kW. The inverse analysis was used to recover the weights, and the fraction of heat radiated, which generate radiative heat fluxes with greater compatibility with the experimental data on three approaches. The first approach is to apply the inverse analysis in each flame separately, obtain the weights and then correlates them with the flame power. In method application for each flame, the maximum deviation between model outcome and experimental data was less than 5%. In a second approach, it was used the inverse analysis to directly obtain the coefficients of a function between the model's parameters and the flame’s power. When used just the flame power as model input parameter to distributes the sources and predict the heat flux, the maximum deviation is 18.6%. In the third approach, the inverse analysis was again used to find the coefficients of the function which correlates the flame power with model parameters, in this case was used the experimental flame length for positioning the point sources. For the case in which seven sources are positioned at 2.25 times the experimental length of the flame, the maximum deviation observed was 8.6% and the mean deviation of 2 9%.

A estimativa acurada do fluxo de calor radiativo na região próxima à chama do tipo jato não pré-misturada é necessária para garantir a segurança de pessoas e equipamentos em caso de vazamentos ou processo de descarte na indústria de petróleo e gás. A simulação computacional dos fenômenos físicos envolvidos na transferência de calor e combustão do processo tem alto custo computacional. No presente trabalho é apresentado um estudo para o emprego do modelo de múltiplas fontes ponderadas para a estimativa do fluxo de calor radiativo no campo próximo à chama. O modelo matemático simplificado tem baixo custo computacional e consiste em representar a transferência de calor radiativa por fontes pontuais distribuídas no eixo central da chama. Cada fonte tem um peso proporcional à contribuição de cada região discretizada da chama na transferência de calor por radiação. Para determinar o peso de cada fonte foi utilizada a análise inversa pelo método da Otimização Extrema Generalizada, no qual o fluxo de calor é dado de entrada enquanto o peso de cada fonte é dado de saída. Como dado de entrada foi utilizado o fluxo de calor radiativo medido experimentalmente de um conjunto de 12 chamas, com potência entre 0,139 e 0,554 kW. A análise inversa foi utilizada para recuperar os pesos, e a fração radiante, que geram o fluxo de calor radiativo com maior compatibilidade com os dados experimentais em três abordagens. A primeira abordagem consiste em aplicar a análise inversa em cada chama separadamente, e depois de obtidos os pesos correlaciona-los com a potência da chama. Na aplicação do método em cada chama individualmente, o desvio máximo do resultado do modelo com os dados experimentais é de 5%. Em uma segunda abordagem, foi utilizada a análise inversa para obter diretamente os coeficientes de uma função entre os pesos do modelo e a potência da chama e seu comprimento estimado, o desvio máximo encontrado é de 18,6%. Na terceira abordagem, a análise inversa foi utilizada novamente para encontrar os coeficientes da função que correlaciona a potência da chama com parâmetros do modelo, e neste caso foi utilizado o comprimento experimental da chama para posicionar as fontes pontuais. Para o caso em que sete fontes foram posicionadas a 2,25 vezes o comprimento da chama medido experimentalmente, o desvio máximo observado foi de 8,6% e o desvio médio de 2,9%. / The accurate approximation of radiative heat flux on non-premixed flame in the region adjacent to the flame is required to guarantee the staff and machinery safety on oil and gas industry in cases of fuel leaking or disposal process. The physical phenomena involved on heat transfer and combustion has high computational cost in numerical simulation, in accidental leaking cases it is prohibitive. In the present study it is presented an application of Weighted Multi-Point Sources model to predict the radiative heat flux in the region adjacent to the flame. The simplified mathematical model has low computational cost and simulates the radiative heat transfer by punctual sources placed at the flame axis. Each source has the weight proportional to heat transfer contribution by each flame portion. To set each source's weight, it was used the inverse analysis by the Generalized Extremal Optimization. In inverse analysis the radiative heat fluxes are the input data while the weight of each source is the sought variable. As input data are used the radiative heat fluxes measured from a set of 12 flames, with power between 0.139 and 0.554 kW. The inverse analysis was used to recover the weights, and the fraction of heat radiated, which generate radiative heat fluxes with greater compatibility with the experimental data on three approaches. The first approach is to apply the inverse analysis in each flame separately, obtain the weights and then correlates them with the flame power. In method application for each flame, the maximum deviation between model outcome and experimental data was less than 5%. In a second approach, it was used the inverse analysis to directly obtain the coefficients of a function between the model's parameters and the flame’s power. When used just the flame power as model input parameter to distributes the sources and predict the heat flux, the maximum deviation is 18.6%. In the third approach, the inverse analysis was again used to find the coefficients of the function which correlates the flame power with model parameters, in this case was used the experimental flame length for positioning the point sources. For the case in which seven sources are positioned at 2.25 times the experimental length of the flame, the maximum deviation observed was 8.6% and the mean deviation of 2 9%.

[ES] La situación energética y medioambiental actual requiere de un cambio de paradigma energético, sustituyendo las fuentes de energía primarias basadas en combustibles fósiles (carbón, petróleo y gas natural) por otras más limpias y renovables. Dentro de este cambio, la fotocatálisis, simulando a la fotosíntesis natural, surge como una alternativa para la producción de combustibles a partir de la molécula de CO2. Los hidróxidos dobles laminares (LDH) son una familia de materiales con propiedades prometedoras para su aplicación en el campo de la fotocatálisis, dada la posibilidad de modular su composición química y morfología dependiendo del método de síntesis. En concreto, los hidróxidos dobles laminares que contienen Ti en su estructura (un elemento bien conocido por el uso del TiO2 en fotocatálisis) pueden tener aplicación en este campo, aunque debido a la dificultad de incorporación de un metal tetravalente en la estructura LDH, no han sido muy estudiados. En este contexto, la presente Tesis Doctoral ha probado la aplicabilidad de estos LDH conteniendo Ti en su composición en procesos de transformación fotocatalítica de la molécula de CO2 para dar lugar a combustibles solares y otras moléculas de interés para la industria. Además, se han estudiado los beneficios de la incorporación de un número elevado de elementos metálicos en un mismo material LDH (hasta siete en el caso de la presente Tesis) en la actividad fotocatalítica de dicho material. / [CA] La situació energètica i mediambiental actual requereix d'un canvi de paradigma energètic, substituint les fonts d'energia primàries basades en combustibles fòssils (carbó, petroli i gas natural) per altres més netes i renovables. Dins d'este canvi, la fotocatàlisis, simulant a la fotosíntesi natural, sorgeix com una alternativa per a la producció de combustibles a partir de la molècula de CO2. Els hidròxids dobles laminars (LDH) són una família de materials amb propietats prometedores per a la seua aplicació en el camp de la fotocatàlisis, donada la possibilitat de modular la seua composició química i morfologia depenent del mètode de síntesi. En concret, els hidròxids dobles laminars que contenen Ti en la seua estructura (un element ben conegut per l'ús del TiO2 en fotocatàlisis) poden tindre aplicació en este camp, encara que a causa de la dificultat d'incorporació d'un metall tetravalent en l'estructura LDH, no han sigut molt estudiats. En este context, la present Tesi Doctoral ha provat l'aplicabilitat d'estos LDH, contenint Ti en la seua composició, en processos de transformació fotocatalítica de la molècula de CO2 per a donar lloc a combustibles solars i altres molècules d'interès per a la indústria. A més, s'han estudiat els beneficis de la incorporació d'un nombre elevat d'elements metàl·lics en un mateix material LDH (fins a set en el cas de la present Tesi) en l'activitat fotocatalítica d'este material. / [EN] The current energy and environmental problematics call for a shift in the energy paradigm, replacing primary energy sources based on fossil fuels (coal, oil, and natural gas) with cleaner and renewable alternatives. Within this transition, photocatalysis, mimicking natural photosynthesis, emerges as an alternative for fuel production from CO2 molecules. Layered double hydroxides (LDHs) constitute a family of materials with promising properties for application in the field of photocatalysis, given the ability to modulate their chemical composition and morphology depending on the synthesis method. Specifically, LDHs containing Ti in their structure (a well-known element for the TiO2 usage in photocatalysis) may find application in this field, although the incorporation of a tetravalent metal into the LDH structure has been less explored due to its difficulty. In this context, the present Doctoral Thesis has demonstrated the applicability of these Ti-containing LDHs in photocatalytic transformation processes of CO2 molecules to produce solar fuels and other molecules of interest for the industry. Additionally, the Thesis has explored the benefits of incorporating a high number of metallic elements into a single LDH material (up to seven in this case) on the photocatalytic activity of the resulting material. / Molina Muriel, M. (2024). Preparación de materiales basados en hidróxidos dobles laminares como fotocatalizadores [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/203889

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.

Simulações acopladas do escoamento reativo e dos processos de transferência de calor para o estudo de chamas são problemas dispendiosos computacionalmente. A transferência de calor por radiação em processos de combustão, devido às elevadas temperaturas, é o processo de troca energética dominante. Ainda, o comportamento altamente irregular do coeficiente de absorção com o comprimento de onda se constitui em uma dificuldade adicional na modelagem da transferência radiativa em meios participantes. Para contornar essa dificuldade modelos espectrais foram desenvolvidos com o objetivo de simular o comportamento de um gás real. Dentre esses modelos destacam-se o gás cinza (GG: Gray Gas), o mais simples, que negligencia o comportamento espectral do coeficiente de absorção, e o modelo da soma-ponderada-de-gases-cinza (WSGG: Weighted-Sum-of-Gray-Gases) onde a integração sobre todo o espectro é substituída por um número finito de gases cinza. Com o avanço de ferramentas computacionais, principalmente códigos CFD (Computational Fluid Dynamics), abordagens computacionais se tornaram atrativas frente ou em complemento às abordagens experimentais. Desse modo, o presente trabalho tem por objetivo a aplicação dos modelos WSGG e GG com novas correlações na simulação detalhada de chamas difusivas laminares de metano diluído com dióxido de carbono e nitrogênio com o código CFD comercial ANSYS/Fluent. Foram desenvolvidas rotinas de usuário (UDF: User-Defined Functions) para o acoplamento dos modelos espectrais ao código CFD. A verificação das rotinas de usuário foi realizada comparando os resultados obtidos via simulação Fluent com dados obtidos pelo modelo WSGG com um código FORTRAN próprio desenvolvido pelo grupo de pesquisa do Laboratório de Radiação Térmica (LRT/UFRGS) para o problema unidimensional de superfícies negras e infinitas preenchidas por um meio não-isotérmico e não-homogêneo. Os erros encontrados para o fluxo de calor radiativo nas superfícies e para o termo fonte radiativo ao longo do meio foram da ordem de 1% indicando o funcionamento correto das rotinas UDF acopladas ao Fluent. Por fim, as rotinas foram aplicadas na simulação numérica para chamas de potência constante com diluição dos reagentes e os dados obtidos com a solução numérica foram comparados com dados experimentais para a fração radiante e fluxo de calor radiativo. Os desvios médios encontrados para o fluxo de calor radiativo ficaram em torno de 10% para todas as chamas, excetuando as chamas com diluição de CO2 de 30%, 40% e 50%, em volume, para as quais os desvios médios ficaram em torno de 15%. O termo fonte para as chamas apontou para a predominância da emissão do meio em relação à absorção. Todas as chamas estudadas se encontram no regime opticamente fino (optically thin) para o qual, segundo apontam estudos da literatura, a escolha do modelo espectral possui impacto pequeno em resultados globais da chama como a temperatura e a concentração das espécies na mistura. Nesse aspecto os resultados encontrados concordaram com a previsão da literatura, entretanto para a transferência radiativa, o modelo GG se mostrou sensivelmente menos preciso em comparação ao modelo WSGG, principalmente para a fração radiante e para o fluxo radiativo na região da pluma aquecida, indicando a dependência do modelo espectral adotado. / Coupled simulations of the reactive flow with the heat transfer processes for flame studying are computationally demanding problems. The radiative transfer in combustion processes is the main heat transfer mechanism due to the high temperatures involved. However, the highly irregular behavior o f the absorption coefficient with the wavenumber composes in an additional difficulty on modeling the radiative transfer in participating media. In order to overcome this issue, spectral models were developed with the objective of simulate the behavior of real gases. Some of the most known models are the gray gas (GG) for which the spectral behavior of the radiative properties of the medium is neglected and the weighted-sum-of-gray-gases (WSGG) for which the integration over the entire spectrum is replaced by a summation over a finite number of gray gases with constant absorption coefficients. With the development of computational tools, mainly Computational Fluid Dynamics (CFD) codes, numerical approaches became attractive instead or in complement of experimental set ups. In this way, the present work aims to couple the WSGG and the GG models with new correlations in a detailed simulation of diffusive laminar flames of methane diluted with carbon dioxide and nitrogen with the commercial CFD code ANSYS/Fluent. User-defined functions (UDF) were developed to the coupling of the spectral models. The verification was carried out through the WSGG model by comparing the Fluent solution with a solution obtained with a FORTRAN code developed by the Thermal Radiation Laboratory (LRT/UFRGS) research group for the one-dimensional system of black surfaces filled with a non-homogeneous and non-isothermal medium. The deviations for the radiative heat flux for the walls and the radiative heat source along the domain were of 1% or less, indicating the correct coupling between the UDF routines and the CFD code. Finally, the UDF were applied in the solution of constant power flames with fuel diluted with carbon dioxide and nitrogen. The obtained data was then compared with experimental measurements for the radiant fraction and the radiative heat flux along the flame axis. The average deviations found were in order of 10% for all flames, except for the flames with 30%, 40% and 50% of CO2 dilution, in volume, for which the deviatioms found were in order of 15%. The radiative heat source was plotted and indicated for the medium emission predominance in comparison with the medium absorption. All flames studied were optically thin flames for which, studies pointed, the spectral model have minor impact over global results as flame temperature and mixture concentration. For this aspect the results found showed agreement with the literature studies predictions, however the GG model showed itself less accurate in comparison with the WSGG model for the radiant fraction and the radiative heat flux computations. So the spectral models have influence on the radiative transfer even if its effect on flame structure can be negligible.

Simulações acopladas do escoamento reativo e dos processos de transferência de calor para o estudo de chamas são problemas dispendiosos computacionalmente. A transferência de calor por radiação em processos de combustão, devido às elevadas temperaturas, é o processo de troca energética dominante. Ainda, o comportamento altamente irregular do coeficiente de absorção com o comprimento de onda se constitui em uma dificuldade adicional na modelagem da transferência radiativa em meios participantes. Para contornar essa dificuldade modelos espectrais foram desenvolvidos com o objetivo de simular o comportamento de um gás real. Dentre esses modelos destacam-se o gás cinza (GG: Gray Gas), o mais simples, que negligencia o comportamento espectral do coeficiente de absorção, e o modelo da soma-ponderada-de-gases-cinza (WSGG: Weighted-Sum-of-Gray-Gases) onde a integração sobre todo o espectro é substituída por um número finito de gases cinza. Com o avanço de ferramentas computacionais, principalmente códigos CFD (Computational Fluid Dynamics), abordagens computacionais se tornaram atrativas frente ou em complemento às abordagens experimentais. Desse modo, o presente trabalho tem por objetivo a aplicação dos modelos WSGG e GG com novas correlações na simulação detalhada de chamas difusivas laminares de metano diluído com dióxido de carbono e nitrogênio com o código CFD comercial ANSYS/Fluent. Foram desenvolvidas rotinas de usuário (UDF: User-Defined Functions) para o acoplamento dos modelos espectrais ao código CFD. A verificação das rotinas de usuário foi realizada comparando os resultados obtidos via simulação Fluent com dados obtidos pelo modelo WSGG com um código FORTRAN próprio desenvolvido pelo grupo de pesquisa do Laboratório de Radiação Térmica (LRT/UFRGS) para o problema unidimensional de superfícies negras e infinitas preenchidas por um meio não-isotérmico e não-homogêneo. Os erros encontrados para o fluxo de calor radiativo nas superfícies e para o termo fonte radiativo ao longo do meio foram da ordem de 1% indicando o funcionamento correto das rotinas UDF acopladas ao Fluent. Por fim, as rotinas foram aplicadas na simulação numérica para chamas de potência constante com diluição dos reagentes e os dados obtidos com a solução numérica foram comparados com dados experimentais para a fração radiante e fluxo de calor radiativo. Os desvios médios encontrados para o fluxo de calor radiativo ficaram em torno de 10% para todas as chamas, excetuando as chamas com diluição de CO2 de 30%, 40% e 50%, em volume, para as quais os desvios médios ficaram em torno de 15%. O termo fonte para as chamas apontou para a predominância da emissão do meio em relação à absorção. Todas as chamas estudadas se encontram no regime opticamente fino (optically thin) para o qual, segundo apontam estudos da literatura, a escolha do modelo espectral possui impacto pequeno em resultados globais da chama como a temperatura e a concentração das espécies na mistura. Nesse aspecto os resultados encontrados concordaram com a previsão da literatura, entretanto para a transferência radiativa, o modelo GG se mostrou sensivelmente menos preciso em comparação ao modelo WSGG, principalmente para a fração radiante e para o fluxo radiativo na região da pluma aquecida, indicando a dependência do modelo espectral adotado. / Coupled simulations of the reactive flow with the heat transfer processes for flame studying are computationally demanding problems. The radiative transfer in combustion processes is the main heat transfer mechanism due to the high temperatures involved. However, the highly irregular behavior o f the absorption coefficient with the wavenumber composes in an additional difficulty on modeling the radiative transfer in participating media. In order to overcome this issue, spectral models were developed with the objective of simulate the behavior of real gases. Some of the most known models are the gray gas (GG) for which the spectral behavior of the radiative properties of the medium is neglected and the weighted-sum-of-gray-gases (WSGG) for which the integration over the entire spectrum is replaced by a summation over a finite number of gray gases with constant absorption coefficients. With the development of computational tools, mainly Computational Fluid Dynamics (CFD) codes, numerical approaches became attractive instead or in complement of experimental set ups. In this way, the present work aims to couple the WSGG and the GG models with new correlations in a detailed simulation of diffusive laminar flames of methane diluted with carbon dioxide and nitrogen with the commercial CFD code ANSYS/Fluent. User-defined functions (UDF) were developed to the coupling of the spectral models. The verification was carried out through the WSGG model by comparing the Fluent solution with a solution obtained with a FORTRAN code developed by the Thermal Radiation Laboratory (LRT/UFRGS) research group for the one-dimensional system of black surfaces filled with a non-homogeneous and non-isothermal medium. The deviations for the radiative heat flux for the walls and the radiative heat source along the domain were of 1% or less, indicating the correct coupling between the UDF routines and the CFD code. Finally, the UDF were applied in the solution of constant power flames with fuel diluted with carbon dioxide and nitrogen. The obtained data was then compared with experimental measurements for the radiant fraction and the radiative heat flux along the flame axis. The average deviations found were in order of 10% for all flames, except for the flames with 30%, 40% and 50% of CO2 dilution, in volume, for which the deviatioms found were in order of 15%. The radiative heat source was plotted and indicated for the medium emission predominance in comparison with the medium absorption. All flames studied were optically thin flames for which, studies pointed, the spectral model have minor impact over global results as flame temperature and mixture concentration. For this aspect the results found showed agreement with the literature studies predictions, however the GG model showed itself less accurate in comparison with the WSGG model for the radiant fraction and the radiative heat flux computations. So the spectral models have influence on the radiative transfer even if its effect on flame structure can be negligible.

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.

Made available in DSpace on 2014-12-17T14:57:46Z (GMT). No. of bitstreams: 1 MirtaniaAL.pdf: 1813042 bytes, checksum: 3d26a66f383220285205f479f11b4656 (MD5) Previous issue date: 2009-09-27 / This research is based, at first, on the seeking of alternatives naturals reinforced in place of polymeric composites, also named reinforced plastics. Therein, this work starts with a whole licuri fiber micro structural characterization, as alternative proposal to polymeric composites. Licuri fiber is abundant on the Bahia state flora, native from a palm tree called Syagrus Coronata (Martius) Beccari. After, it was done only licuri fiber laminar composite developing studies, in order to know its behavior when impregnated with thermofix resin. The composite was developed in laminar structure shape (plate with a single layer of reinforcement) and produced industrially. The layer of reinforcement is a fabric-fiber unidirectional of licuri up in a manual loom. Their structure was made of polyester resin ortofit?lica (unsaturated) only reinforced with licuri fibers. Fiber characterization studies were based on physical chemistry properties and their constitution. It was made by tension, scanning electron microscopy (SEM), x-ray diffraction (RDX) and thermal analyses (TG and DTA) tests, besides fiber chemistry analyses. Relating their mechanical properties of strength and hardness testing, they were determined through unit axial tension test and flexion in three points. A study in order to know fiber/matrix interface effects, in the final composites results, was required. To better understand the mechanical behavior of the composite, macroscopic and microscopic optical analysis of the fracture was performed / Esta investiga??o baseia-se, em principio, na busca de refor?os naturais alternativos aos comp?sitos polim?ricos, tamb?m denominados de pl?sticos refor?ados. Neste sentido, o presente trabalho de pesquisa se inicia com um completo estudo de caracteriza??o micro estrutural da fibra de licuri, essa como proposta alternativa de refor?o aos comp?sitos polim?ricos. A fibra de licuri ? abundante na flora do estado da Bahia e origin?ria da palmeira denominada Syagrus Coronata (Martius) Beccari. Em seguida, foi realizado um estudo do desenvolvimento de um comp?sito laminar somente a base de fibras de licuri, de forma a se conhecer o comportamento da fibra quando impregnada com resina termofixa. O comp?sito foi desenvolvido na forma de estrutura laminar (placa com uma ?nica l?mina de refor?o) e fabricado industrialmente. A l?mina de refor?o ? de um tecido unidirecional de fibras de licuri confeccionado em um tear manual. O comp?sito tem em sua estrutura resina poli?ster ortofit?lica (insaturada) refor?ada somente com fibras de licuri. O estudo da caracteriza??o da fibra teve como base a determina??o das propriedades f?sicas, qu?micas e de constitui??o. As mesmas foram realizadas atrav?s dos ensaios de: tra??o uniaxial, an?lise por microscopia eletr?nica de varredura (MEV), ensaio de difra??o de raios X (RDX), an?lises t?rmicas (TG e DTA), al?m da an?lise qu?mica das fibras. Com rela??o ?s propriedades mec?nicas de resist?ncia e rigidez do comp?sito, estas foram determinadas atrav?s dos ensaios de tra??o uniaxial e flex?o em tr?s pontos. Um estudo com o objetivo de conhecer a influ?ncia da ader?ncia da interface fibra/matriz, na resposta final do comp?sito, se fez necess?rio. Para melhor entendimento do comportamento mec?nico do comp?sito, an?lises macrosc?picas e microsc?picas ?pticas da fratura foram realizadas

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.