O presente trabalho apresenta uma comparação das eficiências energética e exergética de três alternativas diferentes de montagem de sistemas de co-geração de energia que empregam o gás natural como combustível. O caso estudado é uma situação real de um curtume que utiliza lenha como combustível para geração de sua energia térmica e compra energia elétrica da distribuidora. A primeira alternativa é composta por um motor de combustão interna a ciclo Otto, a segunda emprega uma turbina a gás a ciclo Brayton e a terceira usa um gerador de vapor superaquecido juntamente com uma turbina a vapor em contrapressão. Para simulação das plantas é utilizado um software comercial que foi alterado, ficando capacitado para calcular exergias físicas e químicas de plantas de utilidades, cálculos exergoeconômicos, tais como taxas de custo das correntes e eixos, e também de determinar para cada equipamento os valores das variáveis usadas para a avaliação e o melhoramento do desempenho de plantas. Após as modificações o programa foi testado resolvendo um problema conhecido e cuja solução encontrada foi comparada com o da bibliografia. Foram obtidos resultados muito próximos ao padrão usado para comparação, os desvios encontrados são desprezíveis para os objetivos deste trabalho. A planta com melhor eficiência energética e exergética é aquela com motor ciclo Otto, seguida pela montagem que emprega a turbina a gás e finalmente, o conjunto gerador de vapor e turbina a vapor. São calculados os custos de insumos: investimentos necessários para construção das plantas, os custos de operação e manutenção e os custos com combustível. A partir destes valores são resolvidos os sistemas de equações para determinar as taxas de custos das correntes e eixos. Para as três plantas as variáveis de decisão são modificadas buscando aumentar a eficiência exergética e diminuir os custos dos produtos de co-geração. A metodologia para avaliação e melhoramento do desempenho das plantas é empregada e observa-se que é na terceira configuração é onde as modificações tiveram maior impacto, seguida pela primeira planta e por último a segunda alternativa.

Tubos de aço inoxidável duplex (AID) são utilizados na indústria petroquímica para transporte de óleo e gás em ambientes agressivos e em condições de corrosão sob tensão. O processo de soldagem destes tubos impõem grandes dificuldades, como longos períodos de soldagem e formação de intermetálicos devido a baixas taxas de resfriamento. Visando investigar alternativas para este processo, o Laboratório de Metalurgia Física (LAMEF) desenvolveu um equipamento (MASF 1500) capaz de realizar o processo de Soldagem de Tubos por Atrito com Anel Rotativo (STAAR). Este processo é mais rápido e não prejudica as propriedades mecânicas como a soldagem convencional, porém apresenta o viés da formação de rebarbas no tubo que precisam ser removidas para transporte de fluidos. Para investigar a formação de rebarba, os campos de deformação e a distribuição térmica na peça, um modelo utilizando o método dos elementos finitos (MEF) foi desenvolvido a partir de resultados experimentais. O modelo termomecânico desenvolvido com o software ABAQUS se baseia apenas nos parâmetros de entrada do equipamento, e se mostrou capaz de reproduzir os resultados de soldagem de força e deslocamento com erros abaixo 15% e 4% respectivamente. A formação de rebarba interna dos tubos foi calculada com erros abaixo de 5% e as distribuições de temperatura e deformação plástica foram analisadas. Considerando que a formação de fases intermetálicas é um dos principais problemas na soldagem convencional dos AID, as temperaturas calculadas com o modelo foram utilizadas para estimar a formação de fase sigma no tubo, tendo valores calculados abaixo de 1%. / Duplex stainless steel pipes are used in oil and gas transport while in aggressive environment and stress corrosion cracking conditions. The common welding process of these pipes is difficult due to long welding periods and intermetallic phase formation by slow cooling rate. As an alternative for this process, the Physical Metallurgy Laboratory (LAMEF) at UFRGS University developed a machine (MASF 1500) capable of performing tube friction welding with a rotary ring process (STAAR – Soldagem de Tubos por Atrito com Anel Rotativo). This process is faster and does not diminish original mechanical properties of the component, however the process presents flash formation which must be removed for fluid transportation. To investigate flash formation, strain fields and thermal fields in the welded samples, a FEM model was developed based on experimental results. The thermo-mechanical model developed in ABAQUS environment is solely based in process input parameters, showing results of force and displacement with errors smaller than 15% and 4% respectively. The internal flash formation calculation showed less than 5% error and the thermal distribution and plastic deformation were analyzed. Considering that intermetallic phases are one of the biggest concerns on duplex stainless steel arc welds, sigma phase formation was calculated from model’s thermal output, showing less than 1% sigma formation.

Um processo competitivo no mercado atual exige a união de qualidade e produtividade. Além disto, o aumento da eficiência dos processos vem ao encontro da atual preocupação com os impactos ambientais das atividades produtivas, principalmente no que concerne ao uso de combustíveis fósseis. O objetivo principal deste estudo foi verificar a possibilidade de aumento da produtividade de um forno de reaquecimento, passando de 42 para 45 toneladas por hora de tarugos de aço com qualidade similar à AISI/SAE 52100. Para esta finalidade, foram efetuadas simulações numéricas dos ciclos de reaquecimento utilizando um aplicativo computacional. Inicialmente foi efetuada uma revisão da bibliografia referente ao assunto, seguindo-se de um trabalho prático que constou de um procedimento de instrumentação em um tarugo com qualidade similar à de um aço AISI/SAE 52100 para obtenção de dados reais do processo de reaquecimento e sua utilização para verificação do comportamento do software utilizado e, por fim, simulações numéricas com dois diferentes tempos de reaquecimento, relativos a uma produção de 42 e 45 toneladas por hora. Considerando-se o perfil térmico do tarugo de aço, obtido a partir das simulações, verificou-se que o equipamento em questão encontra-se próximo aos limites de produtividade. Diminuições no tempo de reaquecimento, ou seja, aumentos na taxa de produção, provavelmente comprometeriam a qualidade do processo e do produto. / A productive process to be competitive requires the union of quality and productivity. Furthermore, the increased efficiency of processes encounters the current concern with the environmental impacts of production activities, particularly with regard to fossil fuels consumption and emissions. The main objective of this study was to investigate the possibility of increasing the productivity of a reheating furnace from 42 to 45 tons per hour of AISI/SAE 52100 steel billets. For this purpose, numerical simulations of reheating cycles were performed using a computer application. Initially a literature review on the subject was performed, followed by a practical work that consisted of a instrumentation procedure in a AISI / SAE 52100 billet to obtain real data from reheating process. This data was used to validate the numerical simulation and, finally, two different reheating times relating to a production of 42 and 45 tons per hour were simulated. Considering the thermal profile obtained from the simulations, it is likely that the furnace is within its productivity limits. Reductions in the reheat time, thus increasing the production rate, are likely to compromise the quality of the final product.

O estudo do escoamento ao redor de corpos submersos tem sido foco atual de muitas pesquisas. Em especial, escoamentos ao redor de obstáculos rombudos constituem sistemas de grande interesse prático em vários campos da engenharia e outras ciências como, por exemplo, na engenharia automobilística, na concepção de edificações submetidas à ação do vento, no estudo de correntes fluviais e marítimas sobre pilares de pontes, em risers, cabos de ancoragem e plataformas de captação de petróleo. O presente trabalho teve como objetivo principal o estudo por Simulação Numérica Direta do escoamento ao redor de um cilindro com dois graus de liberdade, submetido a Vibração Induzida por Vórtices (VIV), buscando o esclarecimento a cerca dos fenômenos tridimensionais de formação e evolução do escoamento em função da interação fluido-cilindro. As equações de Navier-Stokes e da continuidade são resolvidas numericamente usando esquemas de diferenças finitas compacto de sexta ordem na discretização espacial, em associação com um método de representação das fronteiras do cilindro, denominado Método de Fronteiras Imersas. Inicialmente foram realizadas simulações para os casos de cilindro fixo e em movimento forçado, onde os resultados mostraram que o código computacional consegue modelar satisfatoriamente os fenômenos físicos de acordo com os resultados encontrados na literatura. Para as simulações com cilindro fixo, em movimento forçado e em VIV, nota-se que as estruturas do escoamento tridimensional apresentam modos com comprimentos de onda alterados em função do movimento imposto. Os resultados das simulações em VIV com dois graus de liberdade, para Re = 1250 mostraram que as estruturas tridimensionais do escoamento sofrem significativas alterações espaciais e temporais em função dos parâmetros do modelo dinâmico. Foram identificados diferentes fenômenos de formação e emparelhamento dos vórtices que são modificados em função das condições impostas ao modelo massa-mola amortecimento. / The study of flow around submerged bodies has been the focus of many researches. In special, the flow around bluff bodies are systems of great practical interest in many fields of engineering and other sciences, e.g., the flow around vehicles, buildings submitted to the action of the wind, river and sea currents over pillars of the bridges, risers and cables in the oil platforms. The present work aims to study through Direct Numerical Simulation (DNS) the flow dynamics around a circular cylinder under Vortex-Induced Vibration with two degrees of freedom, searching for identify the changes in the three-dimensional structures of the street wake imposed by the interactions between the flow and the immersed cylinder. The simulations were performed using a sixth-order compact finite differences scheme associated to a Immersed Boundary Method (IBM) to represent the cylinder surface. Firstly, simulations of the flow around fixed and forced movement were accomplished. The computational code Incompact3d presented results in agreement with the literature. The validation of the code shown that is possible to obtain satisfactory results relate to the street wake behind cylinder and the force coefficients are in agreement with others experimental and numerical works. In case of free oscillating cylinder, several changes in the three-dimensional vortex shedding modes are observed as the parameters of the oscillating dynamic model was changed. It was observed that different three-dimensional structures are formed in the wake depending on the initial conditions of the mass-stiffness damping model. The results of the simulations at Reynolds numbers equal to 1250 presenting some hydrodynamic phenomena as paring and turbulent structures in the flow that are modifying during the time and in the computational domain.

Desde a exploração até o planejamento de lavra, metodologias para quantificação do risco associado ao empreendimento mineiro podem melhorar drasticamente o processo de tomada de decisão. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo. Dentro deste contexto, na mineração, a incerteza geológica é vista como um dos principais fatores que contribuem para o insucesso do projeto. A medição e gerenciamento do risco geológico na avaliação do projeto e tomada de decisão pode ser realizada em varias etapas do empreendimento mineiro, desde a quantificação dos recursos, sequenciamento de lavra até a exaustão das reservas. Além disso, a incerteza geológica pode ser traduzida no risco financeiro associado, proporcionando uma análise sobre a viabilidade econômica do empreendimento. A necessidade da quantificação da incerteza geológica já era reconhecida pela indústria de mineração desde a década de 1970, onde estimativas globais e locais eram consideradas insuficientes para a otimização do planejamento de produção, sequenciamento de lavra ou estratégias de homogeneização de pilhas. Desta forma, a simulação estocástica vem desempenhando um papel crucial na construção de modelos de incerteza de teores em depósitos minerais, proporcionando uma ferramenta para obter uma análise do risco associado. Este trabalho tem como objetivo o desenvolvimento de uma metodologia que contribua para a construção de cenários de planejamento de lavra e avaliação do impacto provocado pelas incertezas geológicas fornecidas a partir de modelos simulados. As respostas buscadas, usando a abordagem sugerida, compreendem a definição de uma metodologia para controle e planejamento de lavra incorporando a incerteza local associada ao minério, avaliando a sensibilidade do retorno financeiro e o impacto no planejamento de lavra, considerando variações e incertezas relacionadas ao modelo geológico e estimativa de teores. / From exploration to mine planning, methodologies for quantifying the risk associated with mining enterprise can dramatically improve decision making process. Small variations in the boundary conditions of a project may have significant impact on the final return of the same. Within this context, in mining, geological uncertainty is seen as a major factor contributing to the failure of the project. The need to quantify and manage the geological risk for project evaluation and decision making can be translated into the need to assess the geological risk of any parameter at all stages, from the resources quantification, mining sequencing to reserves exhaustion. In addition, the geological uncertainty can be translated into financial risk, providing an analysis of the economic feasibility of the project. However, the need for quantification of geological uncertainty was already recognized by the mining industry since the 1970, where local and global estimates were considered insufficient for the optimization of production planning, sequencing mining, strategies of homogenization. Thus, the stochastic simulation has played a crucial role in the construction of grade uncertainty models in mineral deposits, providing a tool to perform risk analysis. This study aims to develop a methodology able to contribute in building scenarios of mine planning and evaluation of the impact caused by geological uncertainties supplied from simulated models. Answers searched using the suggested approach include understand the definition of a methodology to control and mine planning incorporating the local uncertainty associated with the ore, assessing the sensitivity of financial return and mine planning considering variantions and uncertainty related to the geological model and grade estimation.

Calamidades naturais são cada dia mais frequentes e a maior parte delas está relacionada aos danos ocasionados por tempestades severas. Desta forma, são necessários estudos que auxiliam tanto na compreensão desses fenômenos quanto na previsão dos mesmos para que sejam desenvolvidas estratégias mitigatórias e edificações cada vez mais preparadas para as adversidades meteorológicas. Intensas rajadas de vento provocam um dos principais danos provenientes de tempestades severas e podem ser originadas a partir de fenômenos tais como downbursts ou tornados. Downbursts são fenômenos meteorológicos que apresentam distinto perfil de velocidade de vento e turbulência que causam diversos acidentes no mundo e, por meio deste trabalho, verificou-se que são a causa de diversos acidentes também no Brasil, em especial na região da Bacia Amazônica, na Região Sul e Sudeste brasileira e na região próxima ao litoral do nordeste brasileiro. A previsibilidade de downbursts também é abordada nessa pesquisa; buscou-se verificar a relação existente entre altos valores do índice CAPE e baixos valores de CINE na ocorrência de ventos provenientes de downbursts, tendo sido constatado que esses índices não são eficientes se analisados de forma isolada na previsão de downbursts. Por fim, verificou-se que estudos mais aprofundados a respeito das ações de downbursts sobre edificações fazem-se necessários para o Brasil. Por fim desenvolveu-se uma revisão de literatura a respeito da simulação física e numérica do fenômeno, buscando gerar um embasamento teórico para futuros trabalhos nesta área. / Natural disasters are becoming more frequent every day and most of them are related to damages caused by severe storms, thus studies that help both in understanding these phenomena and in predicting are necessary to the development of mitigation strategies and to the development of building increasingly prepared to the weather adversity. Intense wind gusts are responsible for the major damage from severe thunderstorms and can arise from phenomena such as downbursts or tornadoes. Downbursts are meteorological phenomena that exhibit distinctive profiles of wind speed and turbulence causing many accidents in the world and through this work are found to be the cause of many accidents in Brazil too, especially in the Amazon Basin region in the South and Southeast Brazilian’s Regions and in the region next to the coast of Brazilian’s Northeastern. Predictability of downbursts is also addressed in this study, it was valued the relationship between high values of CAPE index and low values of CINE index for the occurrence of winds from downbursts, it was found that these index are not efficient when analyzed in isolation in predicting downbursts. Finally, it was concluded that more intensive studies on the effects of downbursts on structures are up needed in Brazil. Finally it was developed an investigation on literature about the physical and numerical simulation of this phenomena, seeking to generate a theoretical basis for future work in this area.

As estruturas sanduíche são formadas por duas faces externas (geralmente feitas com fibra sintética e resina) e um núcleo no meio. Estas estruturas possuem elevada rigidez em flexão, resistência ao cisalhamento e à compressão, e um peso muito baixo. São muitas vezes produzidos por moldagem líquida, como a infusão. Neste trabalho, foi utilizada resina éster‑vinílica, manta e tecido de fibra de vidro e núcleos de PET e PVC, para obter estruturas sanduíche pelo processo de RTM Light. O principal objetivo desta tese é validar a utilização desse processo para a fabricação dessas estruturas utilizando combinações dos núcleos e da manta ou tecido nas faces. Buscou‑se também estudar numericamente, utilizando o software OpenFOAM, o processo de RTM Light para estruturas sanduíches, sobre o qual há pouco estudo na literatura. A resina éster‑vinílica apresentou propriedades adequadas para utilização nestes compósitos, como viscosidade intermediária e tempo de cura apropriado para o processo de RTM Light, sendo que sua pós‑cura teve um efeito importante nas propriedades finais dos compósitos. Os dois núcleos demonstraram propriedades similares, como resistência à compressão transversal e longitudinal, mas o PVC apresentou uma resistência ao cisalhamento superior ao PET. Quanto ao tipo de reforço, também foram obtidos resultados similares, mas o tecido apresentou uma maior resistência à flexão que a manta. A simulação numérica do processo possibilitou prever o comportamento da resina dentro do molde de tal maneira a identificar padrões de preenchimento, bem como regiões mais propícias à formação de vazios. Sendo, portanto, uma ferramenta importante na construção de moldes, definição dos parâmetros de processo e na escolha das matérias‑primas. / The sandwich structures are comprised of two external faces (usually made of synthetic fiber and resin) and a core in between. These structures have high stiffness in bending, shear and compression strength, with a very low weight. They are often produced by liquid molding, such as infusion. In this work, vynil‑ester resin, glass mats and cloths and PET and PVC cores were used to obtain sandwich structures by Light RTM process. The main objective of this thesis is to validate the use of this process for the manufacture of these structures produced with a combination of cores and mats/clothes on the faces. The study was also based on a numerical study of the RTM Light process using the OpenFOAM software, for sandwiches structures, about which very little is available in the literature. The vinyl ester resin showed suitable properties for use in these composites, such as intermediate viscosity and adequate curing time for the Light RTM process, and its post‑cure had an important effect on the final properties of the composites. Both cores showed similar properties, such as resistance the transverse and longitudinal compression, but PVC showed greater shear strength than PET. The two types of reinforcement also produced similar results, but the use of cloths increased flexural strength. Numerical simulation of the process enabled prediction of the behavior of the resin inside the mold in a way to identify the filling pattern, as well as regions more prone to void formation. Thus, being an important tool for mold construction, definition of process parameters and the selection of raw materials.

Um algoritmo para a simulação numérica de escoamentos laminares incompressíveis tridimensionais de fluidos viscosos é proposto neste trabalho. A equação de conservação da energia, assim como as forças de flutuação devido à ação da gravidade, calculadas através da aproximação de Bousinessq, são implementadas para analisar problemas não-isotérmicos. Emprega-se o esquema de Taylor-Galerkin de um passo para a dicretização no tempo e no espaço, respectivamente. O elemento finito empregado é o hexaédrico isoparamétrico de oito nós. Adota-se um enfoque de pseudo-compressibilidade o que permite escrever a equação da continuidade em termos da pressão, a qual é resolvida implicitamente através do Método dos Gradientes Conjugados com Precondicionamento, enquanto as equações de conservação de quantidade de movimento e de energia são resolvidas explicitamente. A integração das matrizes de elemento é efetuada analiticamente considerando o elemento não distorcido e usando um único ponto de integração para avaliar a matriz Jacobiana. O programa é codificado na linguagem FORTRAN 90 e vetorizado para aproveitar as características dos processadores vetoriais existentes nos atuais supercomputadores. / The formulation and application of an algorithm to simulate numerically 3-D incompressible laminar flows is presented in this work. The energy equation and buoyancy forces, which are calculated by Bousinessq approximation, are also implemented to analyse non-isothermic problems. The one-step Taylor-Galerkin scheme is used in order to obtain a useful numerical model. The eight node isoparametric hexahedron finite element is employed. A pseudocompressibility approach is adopted, allowing to write the continuity equation in terms of the pressure. The continuity equation is solved implicitly by the Conjugate Gradient Method with a preconditioning technique, whereas the Equation of Momentum and Energy Conservation are solved explicitly. Element matrices are calculated analytically considering undistorted elements and using one integration point to evaluate the Jacobian Matrix. The code is written in FORTRAN 90 language, and is vectorized in order to take advantage of fast vectorial processors existing in modern supercomputers.

Neste trabalho apresenta-se uma análise numérica tridimensional dos campos de temperatura e velocidade em reservatórios térmicos cilíndricos. Estuda-se o fenômeno da convecção natural laminar e a estratificação vertical da temperatura, que tem grande importância no desempenho de diversos sistemas térmicos, tais como os sistemas de aquecimento de água por energia solar. O código computacional escrito em linguagem FORTRAN resolve as equações da quantidade de movimento e da energia nas três dimensões pelo Método dos Volumes Finitos em malha estruturada. As simulações tridimensionais transientes do processo de resfriamento mostram a formação e evolução do perfil estratificado de temperatura ao longo do tempo. Observa-se a presença de recirculações na região superior do reservatório e a formação das camadas limites térmica e hidrodinâmica próximo das paredes. Através de estudo experimental realizado em paralelo, determinou-se os perfis de temperatura para diversos instantes do processo de resfriamento de um reservatório térmico cilíndrico. Obteve-se também perfis experimentais de temperatura para o caso de consumo de água quente e ingresso de água fria. Em todos os ensaios, além dos perfis de temperatura, mediu-se também as variações da temperatura ambiente. Utilizou-se os dados experimentais colhidos como condição inicial e de contorno para as simulações, possibilitando assim a confrontação numérico versus experimental para a validação da solução numérica. A partir das diversas simulações realizadas, obteve-se correlações para o coeficiente de transferência de calor e para o grau de estratificação térmica no interior do reservatório. Realizou-se ainda a simulação de casos considerando o reservatório térmico interligado a coletores solares e investigou-se tanto a influência da posição do jato de entrada quanto o uso de um anteparo defronte a este jato na preservação da estratificação térmica. Os resultados mostram que, para certas condições do jato de entrada, tanto a escolha adequada da posição do jato de entrada quanto o uso de um anteparo defronte a este jato permitem construir perfis de temperatura com melhor estratificação térmica. / This work presents a three-dimensional numerical analysis of the temperature and velocity fields in cylindrical thermal storage tanks. It is studied the phenomenon of the laminar natural convection and the vertical stratification of the temperature, that has great importance in the performance of diverse thermal systems, such as solar energy water heating systems. The numerical solution is obtained through the Finite Volume Method in structured mesh. The transient three-dimensional simulations of the cooling process show the formation and evolution of the profile of temperature along the time. It is observed the presence of recirculations in the top of the tank and the formation of walls thermal and hydrodynamic boundary layers. Through an experimental study, it was determined the profiles of temperature for several instants of time of the cooling process of a cylindrical thermal storage tank. It was also obtained experimental profiles of temperature for the case of outlet hot water and inlet cold water, as well as measured the ambient air temperature. In some cases the numerical simulation applied experimental data as initial and boundary condition, in order to compare numerical and experimental results, for the validation of the developed computational code. From the simulations, it was obtained correlations for the heat transfer coefficient and degree of thermal stratification in the interior of the reservoir. Others simulations were made with the thermal thank connected to solar collectors, investigating the influence of the position of the inlet water jet, and the use of an inner flat plate in front of this jet, in the preservation of the thermal stratification. The results show that, for certain conditions of the jet, the adequate choice of its position, as well as the use of a flat plate, allows to obtain more stratified temperature profiles.

Este trabalho tem por objetivo analisar a aplicação do Método Reaction Diffusion Manifold (REDIM) na simulação de chamas axisimétricas, difusivas e laminares. Modelos detalhados de cinética química são usualmente computacionalmente proibitivos para a simulação de fluxos reativos complexos e, portanto, modelos reduzidos são necessários. Ferramentas automáticas de redução de modelos normalmente exploram a estrutura multi-escalar natural dos sistemas de combustão. A nova técnica REDIM aplica o conceito de variedades invariantes para tratar, também, a influência dos processos de transporte no modelo reduzido. Assim, ela supera um problema fundamental na redução de modelos em negligenciar o acoplamento do transporte molecular com o processo termodinâmico. Para o estudo, considera-se como caso teste uma chama de metano/ar (65%CH4 + 35%N2), difusiva, laminar, axisimétrica, a pressão atmosférica para validar a metodologia aqui apresentada. Primeiramente, REDIM uni e bidimensionais são computadas e tabuladas em tabelas lookup. Então, o sistema completo de equações governantes é projetado na REDIM e implementado no código Gascoigne através de uma nova biblioteca addon para tratar as tabelas REDIM. O sistema projetado de equações governantes é discretizado pelo Método de Elementos Finitos (MEF) e resolvido por uma iteração GMRES precondicionada por um método multigrid geométrico. Refinamento de malha local, malha adaptativa e paralelismo foram aplicados para garantir eficiência e precisão. Os resultados numéricos obtidos usando a técnica REDIM mostram boa concordância com os dados de simulação numérica detalhada e dados experimentais reportados na literatura. / The goal of this work is to analyze the use of automatically reduced chemistry by the Reaction-Diffusion-Manifold (REDIM) method in simulating axisymmetric laminar coflow diffusion flames. Detailed chemical kinetic models are usually computationally prohibitive for simulating complex reacting flows and, therefore, reduced models are required. Automatic reduction models approaches usually exploit the natural multiscale structure of combustion systems. The novel REDIM approach applies the concept of invariant manifolds to treat also the influence of the transport processes on the reduced model, which overcomes a fundamental problem of model reduction in neglecting the coupling of molecular transport with thermochemical processes. We have considered a previously well studied atmospheric pressure nitrogen-diluted methane/air flame as a test case to validate the methodology presented here. First, one-dimensional and two-dimensional REDIMs were computed and tabulated in lookup tables. Then, the full set of governing equations are projected on the REDIM and implemented in the object-oriented C++ Gascoigne code with a new add-on library to deal with the REDIM-tables. The projected set of governing equations has been discretized by the Finite Element Method (FEM) and solved by a GMRES iteration preconditioned by a geometric multigrid method. Local grid refinement, adaptive mesh and parallelization are applied to ensure efficiency and precision. The numerical results obtained using the REDIM approach have shown very good agreement with detailed numerical simulations and experimental data.