Esta tese apresenta uma nova metodologia para animação de objetos líquidos em imagens. Contrariamente às técnicas existentes, este método é baseado em um modelo físico, o que proporciona efeitos realísticos. A perspectiva da imagem é obtida com a intervenção do usuário, por um esquema simples de calibração da câmera, o qual permite a projeção da camada da imagem a ser animada sobre um plano horizontal no espaço tridimensional. As equações de águas rasas conduzem a simulação e as informações de altura são projetadas de volta ao espaço da imagem utilizando traçado de raios. Além disso, efeitos de refração e iluminação são aplicados durante este estágio, resultando em animações realísticas e convincentes / This work presents a new methodology for animating liquid objects depicted in a still image. In contrast to existing techniques, the proposed method relies on a physical model to accomplish the animation, resulting in realistic effects. Image perspective is handled through a simple user assisted camera calibration scheme which allows one to project the image layers to be animated onto the horizontal plane in the three-dimensional space. Shallow-Water equations drive the simulation and the resulting height field is projected back to the image space via ray-tracing. Refraction and lighting effects are also accomplished during the ray-tracing stage, resulting in realistic and convincing animations

No presente trabalho, foi desenvolvido um modelo para o escoamento em reator de membrana do tipo contator ativo, através da abordagem de dinâmica dos fluidos computacional (CFD), utilizando o código comercial ANSYS FLUENT. O modelo incluiu todo o módulo de membrana constituído por uma membrana tubular e um casco metálico. A reação modelo estudada foi a hidrogenação parcial de 1,5-ciclooctadieno, realizada pelo bombeamento da mistura reacional, dissolvida em n-heptano, através da membrana, a partir das extremidades do tubo. Como catalisador, considerou-se a presença de nanopartículas de Pd impregnadas na membrana. O meio poroso foi aproximado por leito granular representado pela equação de Ergun, tendo como parâmetros a porosidade e o tamanho de grão da membrana de alfa-Al2O3 . O valor para o tamanho de grão foi adotado como equivalente ao diâmetro de partícula determinado com uso do código aberto de estereologia ImageJ, do instituto Nacional de Saúde dos Estados Unidos. O modelo de turbulência adotado foi o RNG k-epsilon. Um estudo de sensibilidade incluiu simulações comparando escoamento desprezando reações como escoamento reativo, variação da velocidade, alteração da saída do fluxo e ativação de modelo de turbulência no meio poroso. Foram realizadas simulações de defeitos estruturais na membrana, correspondendo a regiões de porosidade alterada, com e sem perda de sua uniformidade azimutal. Conclui-se que a presença de defeitos estruturais que afetem a uniformidade azimutal da membrana pode resultar em sensível alteração do escoamento em CMRs. / This study focused on the development of a model for the flow in a reactor membrane of the type active contactor, approached through computational fluid dynamics (CFD), using the commercial code ANSYS FLUENT. The model included the entire membrane module, consisted of a tubular membrane and a metal shell. The model reaction studied was the partial hydrogenation of 1,5-cyclooctadiene initiated by the pumping of the reaction mixture, dissolved in n-heptane, through the membrane, from the ends of the tube. As a catalyst, the study considered the presence of impregnated Pd nanoparticles in the membrane. The porous medium was approximated by a granular bed as represented by the Ergun equation, having as parameters the porosity and the grain size of the alfa-Al2O3 membrane. The value for the grain size was adopted as equivalent to particle diameter determined through the open source stereology software ImageJ, of the National Institute of Health USA. The turbulence model used was the RNG k-epsilon. A sensitivity study included simulations of flow neglecting and including reactions, speed variation, change the flow outlet and activation of turbulence model in the porous media. Simulations of structural defects in the membrane were performed, defining regions of porosity changes with and without loss of azimuthal uniformity. The conclusion was that the presence of structural defects that affect the azimuthal uniformity of the membrane can result in marked alteration of the flow regime in CMRs.

Análises ecotoxicológicas vêm sendo amplamente empregadas no monitoramento de amostras ambientais, efluentes industriais e substâncias complexas. Com o objetivo de avaliar a toxicidade de fluidos de usinagem usados na indústria de peças metálicas, foram feitos ensaios de toxicidade aguda com espécies de três diferentes níveis tróficos: Vibrio fischeri, Daphnia similis, Daphnia laevis e Danio rerio. As amostras dos fluidos foram também analisadas quanto aos valores de DQO, fenol, pH, cor, densidade e surfactantes. Os parâmetros físicos e químicos encontram-se dentro dos limites da legislação CONAMA 357 (D.O.U., 2005). Os resultados dos bioensaios revelaram que os fluidos de usinagem apresentaram toxicidade aos organismos empregados neste estudo, e que o tratamento das amostras com radiação gama não foi eficiente na diminuição da toxicidade dessas matrizes. Foi observado que a radiação gama potencializou o efeito tóxico dos fluidos estudados. A biodegradação em solo demonstrou ser efetiva para os fluidos de usinagem e bactérias endêmicas foram identificadas e isoladas para possível utilização no tratamento de áreas contaminadas com fluidos de usinagem. O monitoramento e gestão dos resíduos de fluido de usinagem são necessários para a preservação dos organismos aquáticos devido a sua toxicidade elevada. / Ecotoxicological analyses have being used to monitor environmental samples, industrial effluents and complex substances. With the objective to analyze the toxicity of cutting fluids used in the machinery industry, acute toxicity test with species of three different throphic levels: Vibrio fischeri, Daphnia similis, Daphnia laevis e Danio rerio, were performing. The samples of fluids were analyzed by COD, phenol, pH, color, density and surfactants. The physical and chemical parameters are the according with the brazilian law, CONAMA 357 (D.O.U. 2005}. The results of the toxicity tests showed that the cutting fluids have high toxicity to the organisms used in this study and the gamma radiation treatment was not efficient to decrease the matrix. The biodegradation in soil demonstrated be effective to the cutting fluids and the indigenous bacteria were identified and isolated to possible treatment of soils contaminated with these kinds of substances. The monitoring and management of residues of cutting fluids are necessary to preservation of aquatic live, in consequence of their high toxicity.

Neste trabalho é proposta uma metodologia numérica como ferramenta para avaliação de ambientes térmicos com manequins. A simulação de CFD do ambiente térmico em simulador comercial é integrada à simulação do sistema térmico do corpo humano realizada em código acadêmico. As soluções dos fluxos de calor e temperaturas nas peles são retroalimentadas e a transferência de informações é realizada via arquivo. A geometria do ambiente térmico é simplificada para minimizar os efeitos de problemas com a malha computacional na simulação de CFD, permitindo uma melhor análise do método interativo proposto. O manequim é separado em 15 segmentos cilíndricos representando cabeça, pescoço, tronco, braços, antebraços, mãos, coxas, pernas e pés. Cada segmento é subdividido em quatro quadrantes, totalizando 60 zonas, para capturar assimetrias térmicas e aerodinâmicas. Foi conduzido estudo de validação geométrica de manequim virtual pela comparação de resultados dos coeficientes térmicos obtidos com CFD e resultados obtidos de ensaios experimentais da literatura. A qualidade da malha e o tratamento de parede são discutidos. Os resultados tornam evidente que uma geometria simplificada do manequim é suficiente para estudos e avaliações de ambiente térmico e de conforto térmico quando se utiliza técnicas numéricas de CFD. Os resultados a partir da integração dos simuladores mostram que o método numérico pode ser instável nos segmentos com baixo metabolismo e baixa vazão de sangue, como nos pés e mãos. Pretende-se introduzir na metodologia proposta algum mecanismo que identifique automaticamente este fenômeno, para evitar a divergência do método e tornar a ferramenta mais robusta. / In this work is proposed a new numerical methodology as a tool for thermal comfort evaluation. This method promotes the interaction of the thermal environment simulation and the thermal system of the human body simulation. The commercial CFD simulator FLUENT R and an academician code for human body simulation are used. The solutions are fed back and the transfer is made by file. The geometry of the room is simplified to minimize the effects of problems with the computational mesh in the CFD simulation, allowing a better analysis of the proposed interactive method. The dummy is separated into 15 cylindrical segments representing head, neck, trunk, arms, forearms, hands, thighs, legs and feet. Each segment is subdivided into four quadrants, totaling 60 zones, to capture asymmetries in the heat flux field and temperature field. Was conducted a geometric validation of virtual dummy by comparing results of heat transfer coefficients from literature and CFD simulation. The mesh quality and near wall treatment are discussed. The results show that a simplified geometry of the dummy is sufficient for thermal environment studies and evaluations in CFD simulations. The results from coupled simulations show that the numerical method can be unstable in the segments with low metabolism and low blood flow, as the feet and hands. So, its intended to introduce mechanisms in the methodology to automatically identify this phenomenon and to avoid the divergence of the method to make more robust this methodology.

Neste trabalho estudamos as propriedades estruturais e eletrônicas do CO2 supercrítico, iniciando com a avaliação de campos de força balizados por aplicações anteriores de simulação quântica do tipo Dinâmica Molecular de Born-Oppenheimer (BOMD). A aplicação principal é a descrição do espectro de absorção da paranitroanilina (pNA) em CO2 supercrítico. O CO2 supercrítico pode ser considerado como uma ``alternativa verde para os solventes orgânicos convencionais e a busca por solventes mais seguros, juntamente com a crescente consciência sobre a questão ambiental, tem levado a uma ``química verde com o intuito de se buscar soluções sustentáveis. A princípio estudamos três campos de força tradicionais para o CO2, aplicados na região supercrítica. Estes campos de força podem ser validados por meio de simulação de primeiros principios. Iniciamos considerando a condição supercrítica para o CO2 como T = 315 K, = 0.81 g/cm³ e o campo de força clássico de Zhang e Duan. Depois fizemos uma análise consistindo de uma alteração de cargas e também da geometria do CO2, que seria um caso não linear no qual foi considerado um ângulo (O-C-O) = 176° . O estudo do solvatocromismo da pNA em CO2 supercrítico foi feito considerando todas estas situações descritas para o campo de força, avaliando os resultados experimentais e teóricos já existentes. A simulação gera estruturas usando Monte Carlo e são usadas em cálculos de Mecânica Quântica do tipo DFT (CAM-B3LYP). Por fim, para verificar a importância da geometria do sistema, ou seja, a propriedade estrutural, consideramos uma outra geometria para a pNA, diferente da geometria que utilizamos a princípio nas simulações com o CO2 supercrítico. Essa ``geometria modificada\" da pNA foi obtida de uma simulação existente de Born-Oppenheimer e a utilizamos numa simulação Monte Carlo com o caso não linear para o CO2 supercrítico. Os resultados de todas essas simulações nos indicaram que a alteração das cargas e por consequência a alteração da polarização do solvente, não possui muita importância na mudança do espectro de absorção da pNA. Ao se considerar o CO2 não linear, obtivemos resultados um pouco melhor, mas não muito, comparados com a previsão teórica. Mas os resultados mais significativos são os obtidos para a situação em que utilizamos a geometria modificada da pNA. Uma parte do deslocamento do máximo da banda de absorção no espectro da pNA vem com a contribuição eletrostática da interação soluto-solvente e a outra parte vem da mudança estrutural. / In this work we study the structural and electronic properties of CO2 supercritical starting with the evaluation of force fields based on previous ab initio Born-Oppenheimer molecular dynamics (BOMD). The main application is the description of the absorption spectrum of paranitroanilina (pNA) in supercritical CO2. The supercritical CO2 is considered a ``green alternative\" to conventional organic solvents and the search for safer solvents, along with the increasing awareness of environmental issues has led to the interest in ``green chemistry\", seeking sustainable solutions. At first we studied three traditional force fields for CO2, applied in the supercritical region. These force fields can be validated by first principles simulation. We considered the supercritical condition for CO2 as T=315K, =0.81g/cm³ and the classical force field of Zhang and Duan. We also did an analysis consisting of a change of the atomic point charges and the geometry of CO2, including a non-linear case in which an angle (O-C-O)=176° was considered. The study of the solvatochromism of pNA in supercritical CO2 was made considering all these situations, evaluating the theoretical outcome and the experimental results. The simulation generates structures using Monte Carlo and are used in quantum mechanics calculations of DFT (CAM-B3LYP). To verify the importance of geometry in the system, that is, the structural property, we considered another geometry for the pNA geometry different from that we used initially in the simulations with supercritical CO2. This ``modified geometry\" of pNA was obtained from a previous Born-Oppenheimer simulation and was used in a Monte Carlo simulation with the non-linear case for supercritical CO2. The results of all these simulations indicated that the alterations of charge and thus the change in the polarization of the solvent, has no great importance in the change of the absorption spectrum of the pNA. When considering the nonlinear CO2, we obtained slightly better results. But the most significant results are obtained for the situation in which we use the modified geometry of pNA. Part of the shift in the absorption spectrum of the pNA comes with the electrostatic contribution of solute-solvent interaction and the other part comes from the structural change.

The fow around bluff bodies is an essential topic in fuid dynamics. This fow is characterized by large vortical fow regions separating from the surface of the bluff body, and they cause oscillating drag and lift forces on the structure. The fow around an infinite long cylinder is a well-known case being studied in the literature. However, a cylinder with low aspect-ratio piercing the free surface has not been studied much although such an arrangement can be found in many floating offshore structures. In this thesis the results of CFD calculations are presented for a fixed, free surface piercing cylinder with an aspect-ratio L/D equals to 2. The Reynolds number was equal to 4:3 x 104 indicating that the flow is in the sub-critical turbulent regime. An extensive methodology of verication and validation is followed to investigate the reliability of the results. To investigate the effect of the free surface on the calculated hydrodynamic loads, two approaches have been considered: a double-body symmetry condition and an interface capturing Volume-of-Fluid (VoF) method. Additionally, two turbulence models are investigated: a two-equation turbulence model; a non-linear Explicit Algebraic Reynolds Stress Model (EARSM); and the Improved Delayed Detached Eddy Simulation (IDDES) turbulence model. The results are presented in terms of integral results (drag and lift coefficients) and flow visualizations. Based on the results of the cases in which the free surface was modeled as a double body symmetry boundary condition, it is concluded that the model is not suitable for this type of flow as the model damps out the flow dynamics due to over-production of eddy-viscosity. Hence, the characteristic oscillating lift forces are not captured using this turbulence model. However, this turbulence model showed good agreements regarding the flow fields in comparison with experimental PIV measurements. Results of the case modeled with EARSM turbulence model shows better agreement with the experimental results compared with the turbulence model. In the cases where the free-surface is considered, results with the EARSM turbulence model show similar results for the drag forces whereas the lift uctuations were one order of magnitude smaller, compared with the double body case. Lastly, the results using the IDDES turbulence model and free-surface VoF modeling are shown to produce the best comparison with the experimental results, regarding both, integral values and flow field results. / O escoamento ao redor de corpos rombudos é um tópico essencial na dinâmica de fluidos. O escoamento é caracterizado por regiões com grande vorticidade que se separam do corpo e causam oscilações das forças de arrasto e sustentação sobre a estrutura. O escoamento ao redor de cilindros longos é um tema que tem sido extensivamente estudado com muitos trabalhos encontrados na literatura. Entretanto, o cilindro com baixa razão de aspecto perfurante à superfície livre é um caso pouco estudado, apesar desta estrutura ser encontrada em várias estruturas oceânicas flutuantes. Esta dissertação apresenta cálculos numéricos para o escoamento ao redor de um cilindro fixo, que trespassa a superfície livre com razão de aspecto L/D igual a 2. O problema é estudado em um regime subcrítico de turbulência, com número de Reynolds igual a 4:3 x 104. Uma vasta metodologia de verificação e validação foi seguida para avaliar a confiabilidade dos resultados obtidos numericamente. Para investigar os efeitos da superfície livre nas cargas hidrodinâmicas, duas abordagens s~ao consideradas: condição de simetria de duplo corpo e um método de captura de interface Volume of Fluid. Além disso, dois modelos de turbulência foram investigados: o modelo não linear Explicit Algebraic Reynolds Stress Model (EARSM), e o modelo de turbulência Improved Delayed Detached Eddy Simulation (IDDES). Os resultados relacionados aos coeficientes de arrasto e sustentação são apresentados a partir de análise estatística, complementados através de ilustrações que permitem visualizar os campos de escoamento e pressão. Com base nos resultados de casos em que a superfície livre é modelada com uma condição de contorno de simetria, conclui-se que o modelo de turbulência não é adequado para este tipo de escoamento, pois o modelo amortece a dinâmica do escoamento devido à superprodução de viscosidade turbulenta. Consequentemente, as oscilações na força de sustentação não são capturadas usando este modelo. Entretanto, resultados dos campos médios do escoamento mostram-se concordantes com imagens experimentais obtidas com técnicas de PIV - Particle Image Velocimetry. Resultados do caso modelado com o modelo de turbulência EARSM mostram melhores concordâncias na comparação dos parâmetros estatísticos com experimentos do que o modelado com o modelo EARSM. Nos casos em que a superfície livre é modelada com o método VoF, o modelo de turbulência EARSM mostra resultados semelhantes para o arrasto, enquanto as flutuações da sustentação apresentam-se uma ordem de grandeza menores, quando comparadas ao caso de duplo corpo. Resultados usando o modelo de turbulência IDDES e VoF apresentam melhores comparações aos resultados experimentais, tanto para os parâmetros estatísticos quando para as visualizações do escoamento.

Made available in DSpace on 2015-04-14T13:59:11Z (GMT). No. of bitstreams: 1 462338.pdf: 6932658 bytes, checksum: 002a8671d2441b25bba6d1226ae0b910 (MD5) Previous issue date: 2014-08-29 / Gravity currents begin when two fluid bodies, with different densities, collides one against other, starting a relative movement in both, being these processes responsible for the transport of sediment from the continental shelf into the deep ocean and the marine construction relief. When the density of sediments transported by rivers exceeds the density of the ambient fluid, it can plunge and form a submerged current, called hyperpycnal plume. The purpose of this research is to investigate, using Direct Numerical Simulation (DNS), the plunge point dynamics of these currents, being employed the Immersed Boundary Method (IBM) to simulate the slope of the continental shelf. Are investigated the influence of the fall velocity of the particle, the initial concentration of suspended sediment, the declivity of the slope in plunge point formation and consequently in the dynamics of particle sedimentation. These results were compared qualitatively with experimental case found in the literature. Quantitative discrepancies were found arising from the use of two-dimensional configuration, and also uncertainty as to the experiment data. The two-dimensional cases showed that the initial concentration of particles and the slope of the platform directly affects the position of plunge point and so little the sediment deposits profile, as the fall velocity shows opposite behavior. The 2D configuration was inefficient for low Ri numbers, showing the necessity of three-dimensional simulations. / Correntes de gravidade tem in?cio quando corpos fluidos de diferentes densidade se encontram, dando in?cio a um movimento relativo entre eles, sendo estes processos respons?veis pelo transporte de sedimento das plataformas continentais para o oceano profundo e constru??o do relevo marinho. Quando a densidade dos sedimentos trasportados pelos rios excede a densidade do fluido ambiente, a corrente gerada pode afundar, gerando uma corrente submersa, denominada pluma hiperpicnal. A inten??o desta pesquisa ? investigar, atrav?s da Simula??o Num?rica Direta (DNS), a din?mica do mergulho destas correntes, sendo empregado o m?todo das fronteiras imersas (IBM), para simular a inclina??o da plataforma continental. S?o investigadas a influ?ncia da velocidade de queda da part?cula, da concentra??o inicial de sedimentos em suspens?o, da declividade da encosta na forma??o do ponto de mergulho e consequentemente na din?mica de sedimenta??o das part?culas. Estes resultados foram comparados qualitativamente com um caso experimental encontrado na literatura. Discrep?ncias quantitativas foram encontradas decorrentes da limita??o bidimensional empregada aos casos, e tamb?m das incertezas quanto aos dados do experimento. Os casos bidimensionais mostraram que a concentra??o inicial de part?culas e a inclina??o da plataforma afeta diretamente a posi??o do ponto de mergulho e pouco o perfil de dep?sito de sedimentos, j? a velocidade de queda mostra comportamento contr?rio. A configura??o 2D mostrou-se ineficiente para baixos n?meros de Ri, mostrando a necessidade da realiza??o de simula??es tridimensionais.

Este trabalho tem como propósito, investigar a validade da lei de Hess-Murray, através da experimentação numérica, aplicando o método do Design Construtal associado ao método de otimização da busca exaustiva, no escoamento interno em estruturas em forma de T com seção circular. Variação do número de Reynolds, escoamento de fluidos newtonianos e não newtonianos, estrutura em forma de T com paredes impermeáveis e permeáveis, foram as principais características avaliadas para confrontar a lei de Hess-Murray. O estudo proposto assume escoamento tridimensional, laminar, incompressível, regime permanente e propriedades fluidodinâmicas constantes, sendo o regime de escoamento governado pelo número de Reynolds (Re). O objetivo principal consiste em determinar as configurações ótimas que facilitem o escoamento de fluido ou minimizem as resistências ao escoamento, quando a área ocupada pelos dutos (A) e o volume ocupado pelos dutos (V) são mantidos constantes, variando-se as razões de diâmetros (aD) e comprimentos (aL). As equações de conservação de massa e quantidade de movimento, foram resolvidas através do método de volumes finitos. A geometria foi discretizada através de uma malha tridimensional composta por aproximadamente 1.950.000 elementos. Como resultados, obteve-se as geometrias ótimas que apresentaram resistências ao escoamento até 30 vezes menor do que as demais configurações. Além disso, foi possível verificar que a lei de Hess-Murray nem sempre é válida, visto que o sistema adapta sua geometria ótima para cada condição de escoamento, a fim de proporcionar a melhor arquitetura de escoamento para atender ao objetivo de minimizar as resistências ao escoamento em acordo com a Lei Construtal. Esta dissertação avançou no presente estado da arte, pois desenvolveu um modelo tridimensional sem simplificações, aplicado ao sistema de escoamento de fluidos em estrutura em forma de T utilizando o método do Design Construtal, validando os resultados analíticos apresentados na bibliografia e apresentando novas referências que permitem ampliar a complexidade dos sistemas de escoamento bem como a implementação de métodos de otimização mais avançados. / This work investigates, through the numerical experimentation together with the Construtal Design method, the Hess-Murray Law in the internal flow in T-shaped structures with a circular section for the laminar flow of Newtonian and non-Newtonian fluids with impermeable and permeable walls, determining the optimal configurations that facilitate fluid flow or minimize flow resistance. The geometric global constants, the volume occupied by the ducts (V) and the area occupied by the ducts (A), delimit the space occupied by the T-shaped structure and the degrees of freedom, the ratio between the diameter of the parent duct and daughter (aD) and the ratio between parent duct length and daughter (aL), are the main geometric parameters to be evaluated. The proposed study is assumed three-dimensional, laminar, incompressible, permanent and constant fluidodynamic properties being the flow regime governed by Reynolds number (Re). Construtal Design method, associated with the exhaustive search, was used to determine the global geometric constants, degrees of freedom and objective function in the geometric evaluation of the system. The numerical solution of the mass conservation and momentum equations is solved based on the finite volume method. The geometries and mesh of the computational domain was discretized through a three-dimensional composed of approximately 1,950,000 elements. The results show that the optimal geometries that presented resistance to the flow up to 30 times smaller than the other configurations. In addition, it was possible to verify that the Hess-Murray Law is not always valid, since the system adapts its optimal geometry to each flow condition, in order to provide a better flow architecture to meet the objective of minimizing resistance to flow in agreement with the Constructal Law. This work advanced in the present state of the art, since it developed a three-dimensional model without simplifications, applied to the fluid flow system in T-shaped structure using the Construtal Design method, validating the analytical results presented in the bibliography and presenting new references that allow increase the complexity of flow systems as well as the implementation of more advanced optimization methods.

Este trabalho tem como propósito investigar como a variação geométrica de determinados parâmetros envolvidos na construção de uma geometria bifurcada de seção circular, em forma de Y, afeta a resistência ao escoamento, tanto de fluidos newtonianos como não newtonianos. As geometrias estudadas foram construídas utilizando-se o princípio do Design Construtal. Os parâmetros variados foram a relação entre os comprimentos dos dutos pais e filhos, a relação entre os diâmetros dos mesmos dutos, e o ângulo central da estrutura em forma de Y. Para as relações geométricas lineares foram utilizados os valores de 0,5; 0,6; 0,7; 0,8; 0,9 e 1, enquanto para os ângulos, foram utilizados os valores de 155°, 135°, 115°, 95°, 75°, 45°, 25° e 10°. Os fluidos utilizados foram do tipo newtoniano e não newtoniano, dentre estes últimos, foram estudados fluidos pseudoplásticos e dilatantes. O trabalho foi realizado através de simulações numéricas, implementadas com a utilização do software comercial Ansys Fluent, o qual resolve as equações governantes através do método dos volumes finitos. As malhas utilizadas foram do tipo poliédrica. Os resultados indicam que há uma diferença em relação ao que se espera da literatura para as relações entre os diâmetros e os comprimentos. A Lei Hess-Murray indica que estas relações ótimas seriam de 2-1/3 para as relações entre os diâmetros e comprimentos. No presente trabalho, foram determinadas relações entre os diâmetros próximas de 0,6, e entre os comprimentos, iguais a 1. Os ângulos ótimos ficaram localizados no intervalo entre 100° e 135°. / This work aims to investigate how the geometric variation of certain parameters involved in the construction of a bifurcated Y-shaped circular cross-section geometry affects the flow resistance of both Newtonian and non-Newtonian fluids. The geometries studied were constructed using the Constructal Design principle. The parameters were the relationship between the lengths of the daughter and parent ducts, the relationship between the diameters of the same ducts, and the central angle of the Y-shaped structure. For the linear geometric relations, values of 0.5; 0.6; 0.7; 0.8; 0.9 and 1 where used, for the angles, the values of 155 °, 135 °, 115 °, 95°, 75 °, 45 °, 25 ° and 10 ° were used. The fluids used were of the Newtonian and non-Newtonian type, among the latter, pseudo plastic and dilatant fluids were studied. The work was carried out through numerical simulations, implemented with the commercial software Ansys Fluent, which solves the governing equations through the finite volume method. The meshes used were of the polyhedral type. The results indicate that there is a difference in relation to what is expected from the literature for the relationships between diameters and lengths. The Hess-Murray Law indicates that these optimal relations would be 2-1/3 for the relationships between diameters and lengths. In the present work, relationships between the diameters close to 0,6 were found and s equal to 1 between the lengths. The optimum angles were located in the range between 100 ° and 135 °.

A simulação numérica de escoamentos reativos, como a combustão, tem um caráter altamente não-linear devido a presença de diversas reações químicas que acontecem entre as espécies que descrevem o processo de oxidação do combustível. Além disso, tais processos ocorrem a nível molecular, tornando o sistema de equações governantes rígido, o que implica na necessidade de esquemas numéricos de alta ordem bem como malhas finas e passo de tempo pequeno, aumentando consideravelmente o custo computacional. Neste sentido, o uso de mecanismos de oxidação detalhados na simulação numérica é proibitivo, e técnicas de redução química são necessárias de modo a desenvolver modelos reduzidos com menos variáveis e rigidez moderado, mantendo a precisão e abrangência do modelo detalhado. O objetivo do presente trabalho é obter uma comparação dos resultados obtidos para duas técnicas de redução química diferentes, Directed Relation Graph - DRG, baseada no desenvolvimento de mecanismos esqueletos, e a Reaction Diffusion Manifolds - REDIM, baseada na separação das escalas de tempo. Como validação dos modelos propostos, simulações numéricas 1D de chamas pré-misturadas e não pré-misturadas, bem como de reatores homogêneos, são desenvolvidas. Além disso, uma estratégia que une as duas técnicas de redução é apresentada, com o objetivo de ser aplicada em mecanismos cinéticos grandes. / Numerical simulation of reactive flows, such as combustion, has a highly non-linear character due to the presence of several chemical reactions that occur among the chemical species that describe the process of fuel’s oxidation. Besides, such processes occur at a molecular level, making the system of governing equations stiff, which implies in the need of high order numerical schemes as well as fine meshes and small time step, enhancing considerably the computational cost. In this sense, the use of detailed oxidation mechanisms in the numerical simulation is prohibitive, and chemical reduction techniques are needed in order to develop reduced models with less variables and moderate stiffness, while maintaining the accuracy and comprehensiveness of the detailed model. The objective of the present works if to obtain a comparison between two chemical reduction techniques, the Directed Relation Graph - DRG, based on the skeletal mechanisms generation, and the Reaction Diffusion Manifolds - REDIM, based on the separation of time scales. As validation of the proposed models, one-dimensional numerical simulations of premixed and non-premixed flames, as well as homogeneous reactors, are carry out. Besides, a coupled methodology between DRG and REDIM is presented, that will provide a useful tool for simulation of fuels with very large detailed kinetic mechanisms.