A presente discussão e uma extensão do modelo de cinética pontual de nêutrons, onde a reatividade e decomposta em termos de contribuição de curtas e longas escalas de tempo. A primeira representa o controle operacional do reator, enquanto a segunda e devido a alteração da composição química do combustível nuclear, como consequência do burn-up. E um primeiro passo em uma nova direção, uma vez que considera os efeitos dos principais venenos na cinética de nêutrons, ou seja, Xenônio-135 e Sam ario-149. O modelo proposto consiste em um sistema de equações não-lineares acoplado para a densidade de nêutrons, para os precursores de nêutrons atrasados e para as cadeias de decaimento dos venenos produtos de fissão. O sistema de equações e resolvido através de um método de decomposição, que expande os termos não-lineares em uma série infinita, obtendo um sistema recursivo, onde a inicialização da recursão e uma equação linear homogênea e os passos de recursão subsequentes consideram contribuições não-lineares como termo fonte constru dos em passos de recursão anteriores. A construção hierárquica do modelo também e realizada, onde graus espaciais de liberdade são considerados. São apresentados casos de estudos com várias estruturas temporais afim de mostrar a robustez da abordagem atual para este tipo de problema. / The present discussion is an extension to Neutron point kinetics models, where the reactivity is decomposed in a short and a long term contribution. The rst one represents operational reactor control, whereas the second one is due to the change of the chemical composition of the nuclear fuel as a consequence of burn-up. This is a rst step into a new direction where we consider only the e ects of the principal neutron poisons on neutron kinetics, i.e, Xenon-135 and Samarium-149. The proposed model consists in a system of coupled nonlinear equations for the neutron density, the delayed neutron precursors and the neutron poison decay chains. The equation system is solved using a decomposition method, which expands the non-linear terms in an in nit series, obtaining a recursive system, where the recursion initialization is a homogeneous linear equation and the subsequent recursion steps consider the non-linear contributions as source terms constructed from previous recursion steps. A hierarchical construction of the model is also performed, where spatial degrees of freedom are considered. We present case studies with severe time structure in order to show the robustness of the present approach for this kind of problems.

A pesquisa voltada para a retroanálise de taludes de mineração está se tornando cada vez mais importante. As dimensões das atuais cavas de mineração visam sempre obter a melhor relação entre custo, benefício e segurança. Isso exige que por vezes o trabalho ocorra próximo destes limites. Neste trabalho o estudo de retroanálise consiste em um talude misto composto por aluvião e rocha em uma mineração de carvão localizada no município de Arroio dos Ratos, RS. O objetivo foi realizar um estudo de retroanálise onde se utilizou dois métodos de análise método dos Elementos Finitos (MEF) e o método do Equilíbrio Limite (EL). A metodologia adotada foi dividida em três fases: uma etapa de campo, para coleta do material e dados geométricos do talude, uma etapa de laboratório para caracterizar todo o material e a última etapa computacional onde foram gerados os modelos e simulações numéricas através dos softwares comerciais da empresa RocScience, sendo eles o SLIDE e o RS2. Os resultados encontrados mostraram que a porção do talude composta por aluvião apresentava seis camadas intercaladas de areias bem graduadas e argilas de média a alta plasticidade. Em relação as simulações, observou uma diferença entre o MEF e EL devido ao grande número de camadas presentes no talude e a interação entre elas. Foram criados quatro diferentes cenários para observar o comportamento do talude, com presença e ausência de água. O MEF mostrou uma superfície de ruptura próxima a encontrada porem com a presença de água o Fator de Segurança (FS) sofreu uma grande depreciação. O EL considerou uma superfície de ruptura satisfatória, nas simulações com presença de água FS sofreu uma depreciação mínima mesmo com máxima saturação de água no talude. Apesar das diferenças ambos os métodos são confiáveis mas devem ser usados de forma adequado ao tipo de problema existente. Situações mais simples e com taludes mais homogêneos o EL trona-se mais indicado, para problemas de maior complexidades e interações entre camadas o MEF se mostrou satisfatório e com boa resposta. O resultados através do MEF se aproximaram mais da situação real encontrada em campo. / A research focus on retro analysis on mine slopes is becoming increasingly important. The dimensions of current mining caves aim to always obtain the best cost benefit ratio as well as the highest safety level. This demands that the work sometimes is done within these limits. In this paper the study of retro analysis is made up of a mixed slope composed of alluvial and rocks, in a coal mine in Arroio dos Ratos, RS. The goal was to carry out a retro analysis study using two analysis methods: The Finite Element Method (FEM) and Limit Equilibrium analysis (LE). The methodology adopted was split into three phases: a field stage for the sample collection and slopes geometric , a laboratory stage to classify all the material and a final computacional stage were model and numerical simulations were created by way of commercial software, Slide and RS2, by RocScience. The results showed that the portion of the slope composed of alluvium had six layers interspersed with well graded sands and high plasticity clays. The simulations observed a difference between the FEM and LE due to the number of layers present in the slope and the interaction among them. Four different scenarios were created in order to observe the behavior of the slope. The FEM presents surface failure similar to one found in the field. The LE always takes into consideration the overall critical surface which is not observed in the field, it also show high security factor levels above one 1.3 that does not reflect a real situation since the slope is collapsed. In spike of the differences, both methods are reliable but they must be used according to the existing problems. For simpler situations with homogeneous slopes LE is recommended while move complex problems and interactions among lawyers calls for FEM.

Os processos de conversão de energia tendem a considerar cada vez mais restrições econômicas e ambientais, tornando-se necessário o entendimento da interação entre combustão e turbulência. Este trabalho tem como objetivo o desenvolvimento de soluções analíticas para a fração de mistura de uma chama difusiva, sob forma de um jato turbulento axissimétrico. Foi desenvolvida, também, uma metodologia analíticonumérica para a determinação das frações mássicas dos componentes, considerando uma reação de combustão de dois passos. Os resultados foram comparados com dados experimentais encontrados na literatura para uma chama de hidrogênio H2=N2 (50/50% em volume). De modo geral, os resultados obtidos foram satisfatórios frente aos dados experimentais, sendo a principal limitação o fato das expressões analíticas obtidas não serem capazes de representar o jato próximo à saída do bocal, sendo válidas a partir de, aproximadamente, x=d > 10, onde x é a coordenada ao longo do comprimento do jato. A principal vantagem do método empregado neste trabalho é a diminuição da complexidade do sistema de equações a ser resolvido numericamente. / The energy conversion processes tend to consider even more economical and environmental constraints, making it necessary to understand the interaction between combustion and turbulence. This study aims at the development of analytical solutions for the mixture fraction of a diffusive flame in the form of an axisymmetric turbulent jet. It was also considered an analytical-numerical approach for the determination of the mass fractions of the compounds, for a two-step reaction. The results were compared with data found in literature for a hydrogen flame H2=N2 (50/50 % by volume). Overall, the results were satisfactory when compared with the experimental data, however the principal limitation was the fact that the analytical expressions were not able to represent the jet near the nozzle exit, being the solution valid from, approximately, x=d > 10, where x is the coordinate along the jet length. The main advantage of the method employed in this work is the decrease in the complexity of the equations system to be solved numerically.

O estudo de parâmetros cinéticos para a queima de carvão se dá através de utilização de equipamentos laboratoriais como o forno de queda livre (DTF). O presente trabalho trata do desenvolvimento e análise de um modelo numérico unidimensional de um DTF. A modelagem é realizada a partir de dados da literatura, e as equações discretizadas que caracterizam os fenômenos do escoamento, da transferência de calor e da combustão do carvão foram programadas na linguagem FORTRAN 90. Os resultados são comparados com um modelo de referência e com dados experimentais, e apresentaram boa concordância. Um estudo de sensibilidade é realizado para entender o comportamento da queima de carvão frente a mudanças de alguns parâmetros de operação do DTF. A análise de sensibilidade mostrou a versatilidade do modelo e, consequentemente, do experimento. A partir da variação da concentração de O2, da temperatura de operação e das vazões de entrada pode-se obter uma gama de resultados. O uso do modelo aqui apresentado concomitante com resultados experimentais resulta em um estudo aprofundado dos parâmetros cinéticos de combustão de carvão. / The study of kinetic parameters for coal combustion occurs through use of laboratory equipment such as drop tube furnace (DTF). The present work shows the development and analysis of a one-dimensional numerical model of a DTF. The modeling is developed from literature works, and the discrete equations that characterize the flow, heat transfer and coal combustion are programmed in FORTRAN 90 language. The results are compared with a reference model and experimental data, and showed good agreement. A sensitivity study is performed to understand the behavior of coal combustion due to changes of some working parameters of the DTF. The sensitivity analysis showed the versatility of the model and thereby the experiment. From the variation of the O2 concentration, operating temperature and input flow rates a range of results can be obtained. Using the model presented here together with experimental results leads to a detailed study of the kinetic parameters of coal combustion.

Na área de simulação de processos, existe uma visível tendência da migração das ferramentas seqüenciais modulares, que hoje são as mais amplamente utilizadas, para as baseadas em equações. Uma das principais vantagens do paradigma baseado em equações ou simultâneo é que este se mostra eficiente na solução de problemas de simulação, otimização, estimação de parâmetros e reconciliação de dados, todos baseados em um mesmo conjunto de modelos, evitando retrabalho de modelagem. Porém, a tecnologia simultânea também apresenta algumas deficiências, onde destacam-se os problemas de robustez tanto na modelagem quanto na obtenção de resultados numéricos. Este trabalho tem como objetivo reunir e desenvolver técnicas que permitam reduzir estas deficiências. Para tanto, as técnicas conhecidas para depuração de sistemas de equações que representam problemas estacionários e dinâmicos foram estudadas em detalhe. Pôde-se observar que para o caso estático os métodos disponíveis para depuração de modelos, já se apresentam em um nível bem desenvolvido. Já para o caso dinâmico, onde há uma maior complexidade, as técnicas conhecidas encontram-se em um nível de desenvolvimento muito menor. Neste ponto encontram-se as principais contribuições deste trabalho. / In the field of process simulation the movement from the sequential modular tools, which are currently the most widely used, to the equation based approach is clear. One of the key advantages of the equation based or simulatneous approach is that using a single model one can solve simulation, optimization, parameter estimation, and optimization problems. This fact avoids modeling rework for each application. However, the simultaneous technology has problems regarding modeling and solving robustness. This work aims to group and develop methods capable of minimize these deficiencies. In order to achieve this goal, available debugging approaches for both steady-state and dynamic system of equations were studied in detail. For the steady-state case well stablished debugging techniques are known. For dynamic models, where the complexity is higher, the analysis and debugging methods are much less mature. This was the source for the major contributions of this work.

O desenvolvimento de técnicas que deem sobrevida a navios com mais de 20 anos de navegação interior é crucialmente importante para um país com a demanda de mobilidade como a que tem o Brasil. Nesse sentido, uma técnica que pode ser citada é a determinação da vida útil potencial (VUP), que propõe uma forma de avaliar o aumento da vida útil de embarcações de navegação interior. Na busca por índices globais que identifiquem as condições do navio, considera-se o cálculo do momento fletor último da seção mestra do navio. O presente trabalho apresenta, de forma detalhada, a aplicação dos principais métodos de cálculo do momento fletor último em três casos de estudo: um perfil caixão, o perfil reduzido de embarcação de fundo duplo (estudada por Nishihara, 1983) e a seção mestra de uma embarcação típica utilizada para navegação na Bacia do Sudeste (NT Guarita). Os métodos empregados para o cálculo do momento fletor último são: métodos direto e iterativo propostos pela Associação Internacional de Sociedades Classificadoras (IACS 2008), Método da Resistência Direta, Método da Largura Efetiva e análises não lineares através do Método dos Elementos Finitos. Com o intuito de obter uma resposta mais coerente, é proposta uma modificação no Método Iterativo [IACS, 2008]. O método modificado é então utilizado no estudo de sensibilidade do momento fletor último frente à redução de espessura das chapas de convés (considerando a seção mestra da embarcação NT Guarita). Os resultados obtidos, bem como as conclusões a respeito, são apresentados ao final do trabalho. Adicionalmente, é proposto um índice para uso na avaliação do nível de não linearidade da resposta dos elementos locais, quando submetidos a reduções de espessura. / The increasing of ships life extension can only be provided by new methodologies development. A new method to determine the potential ship shelf life is being proposed, being part of the new method, in this document the ultimate bending momentum of the master ship section is study. The study presents the application of the most popular methods to calculate the ultimate bending for three different cases: box section, double hull ship studied by Nishihara, 1983 and master section of a typical boat used in Bacia do Sudeste, NT Guarita. The methods applied to calculate the ultimate bending momentum are: direct and iterative methods proposed by the International Association of Classification Societies, IACS 2008, Direct Strength Method, Method of Effective Width and also non-linear analysis using the Finite Element Method. In order to obtain a more coherent response, a change in the Iterative Method [IACS, 2008] has been proposed. The new method is used to study the sensitivity of the ultimate momentum to the reduction of deck thickness applied the NT Guarita. The results and conclusions will be presented in this study. Additionally, will be proposed an index to evaluate the level of non-linearity of the local elements response, when exposed to thickness reduction.

Neste trabalho desenvolvemos um PSE (Ambiente para Soluca-o de Problemas) para gera,c-ao de malhas O programa chamado MAKEGRID fornece uma interface amig avel para a definica-o da regi-ao f 1sica a ser discretizada Esta regi-ao pode ser simplesmente ou multiplamente coneXa V arios geradores de malha foram implementados entre m etodos alg ebricos e el 1pticos in- cluindo o gerador TTM n-ao-homogAeneo para atrac-ao da malha para um ponto ou linha coordenada espec 1fica O programa tamb emoferece ferramentas graficas para a an alise qualitativa da malha gerada Utilizando as id eias do m etodo MAC introduzimos o uso de uma malha de estados para a identifica,cao das c elulas em uma malha generalizada Assim cada ponto da malha possui um valor associado informando se este ponto pertence a um corte uma c elula fict 1cia ou a condi,c-ao de contorno que deve ser aplicada neste ponto Pode-se especificar uma condi,c-ao de contorno diferente para cada ponto da malha para cada vari avel primitiva utilizada O programa retorna arquivos contendo a malha de estados e a discretiza,c-ao das condic es de contorno requeridas Os arranjos co-localizado e diferenciado s-ao apresentados como poss 1veis escolhas para a discretizaca-o dos operadores gradiente divergente e Laplaciano Estes s-ao implementados para malhas cartesianas e generalizadas / A PSE (Problem-Solving Environment) on Mesh Generation was developed The program entitled MAKEGRID provides a friendly interface for the definition of physical region boundaries which can define simply or multiply-connected domains Several options for mesh generators among algebraic and elliptical methods were implemented including the inhomogeneous TTM method for attracting a mesh to a grid point or coordinate line Besides the programoffers graphical tools for a qualitative analysis of the generated grid Based on the ideas of the MAC method we introduce the use of a state grid for identification of cells in a general mesh Each grid point receives a associated value which indicates if this point belongs to a cut branch to a fictitious cells or what boundary condition is required for this point The used methodology permits the definition of a different boundary condition for each grid cell The user can define the boundary cells and the conditions to be applied on these cells in a practical way receiving a state grid and a file containing the discrete formof the specified boundary conditions as result The co-localized and staggered schemes are presented as possible choices for the gradient divergent and Laplacian operators discrete form .

Este trabalho descreve uma metodologia para a geração de programas CNC para otimizar um sistema produtivo. Quanto às informações geométricas, o sistema caracteriza-se por utilizar os ciclos de usinagem que o comando da máquina disponibiliza, por possibilitar a aplicação de ciclos desenvolvidos com recursos de programação parametrizada e por gerar novos ciclos, através de uma linguagem genérica, proposta neste trabalho, em um editor de programação parametrizada. Quanto aos dados tecnológicos, ou seja, as condições de usinagem definidas no código de programação, o sistema permite a otimização dos parâmetros de corte através de um banco de dados com informações realimentadas com resultados obtidos da produção. A análise da viabilidade da metodologia proposta é realizada através da implementação de um protótipo de um sistema dedicado à programação de tornos CNC. O protótipo é composto de aplicativos computacionais desenvolvidos e orientados à otimização do sistema produtivo. A metodologia para geração de funções CAD/CAM foi implementada e testada mostrando resultados eficientes. O protótipo implementado como um conjunto de ferramentas aplicativas integradas é de fácil operação e utilização unindo todos os módulos projetados. / This study describes a methodology for the generation of CNC programs to optimize a productive system. Regarding geometrical information, the system is characterized by the use of machining cycles that the machine tool control provides because a allows the application of cycles developed with parameterized programming and by generation of new cycles, through a generic language, proposed in this study on a parameterized programming editor. As for technological data, i.e. the machining conditions set out in the programming code, the system allows the optimization of cutting parameters through a data base with information fed back with results of production. The feasibility analysis of the proposed methodology is accomplished through the implementation of a prototype of a system dedicated to programming CNC lathes. The prototype consists of computer applications developed and oriented to the optimization of the production system. The methodology for the generation of CAD/CAM functions has been implemented and tested showing effective results. The prototype implemented as an integrated set of applicative tools is of easy operation and use connecting all designed modules.

O presente trabalho apresenta um estudo numérico sobre o desempenho de trocadores de calor com microcanais. Primeiramente, para a verificação da ferramenta numérica, uma comparação com um estudo experimental conhecido é realizado. Em seguida, aperfeiçoamentos disponíveis na literatura são testados e comparados. As simulações numéricas são realizadas com o software Fluent utilizando o Método dos Volumes Finitos (MVF). O resultado numérico da resistência térmica do dissipador do estudo experimental em que a primeira parte desse trabalho foi baseada ficou em 0,097 °C/W (com uma temperatura máxima de 373 K), enquanto o resultado experimental ficou em 0,090 °C/W, o que representa uma diferença de 7,2 %. São realizados aperfeiçoamentos integrando o trocador de calor a um sistema de refrigeração. Isso reduz a temperatura de entrada e, consequentemente, a temperatura máxima encontrada no sistema (350 K para uma mesma potência dissipada) não havendo mudança na resistência térmica, porém trazendo um ganho em confiabilidade. São apresentadas mudanças na geometria do dissipador diminuindo a área de entrada, com uma sensível piora nos resultados da resistência térmica, fazendo com que a resistência aumentasse de 0,097 °C/W para o canal original para 0,272 °C/W para o canal com menor altura. Apresenta-se um trocador de calor cujos microcanais em suas dimensões longitudinais possuem ondulações, porém os efeitos esperados de melhor mistura e maior área de contato para o fluxo de calor não foram observados para as condições de contorno do experimento original; a resistência térmica calculada foi de 0,102 °C/W. A última análise apresenta uma idealização onde o fluxo de calor é distribuído uniformemente no volume de silício, encontrando-se uma resistência térmica de 0,084 °C/W. / This paper presents a numerical study on the performance of heatsinks with microchannels. For the validation of the numerical tool, a comparison with a known experimental study is carried out, then improvements available in the literature on the subject are tested and compared. The numerical simulations are performed with the Fluent software using the Finite Volumes Method (MVF). The numerical result of the thermal resistance of the heatsink from the experimental study which the first part of this paper is based was 0.097 °C/W, whereas the experimental result was 0.090 °C/W (with a maximum temperature of 373 K), representing a difference of 7.2%. An improvement is performed by integrating the heatsink to a HVAC (Heating, ventilation and air conditioning) system, which lowers the inlet temperature and consequently lowers the highest temperature found in the system (350 K for the same dissipated power) without changing the thermal resistance, but providing gains in reliability. Changes in the heatsink geometry are presented by decreasing the inlet area, with a significant worsening in the thermal resistance results, causing the resistance to vary from 0.097 °C/W for the original channel to 0.272 °C/W for the channel that has the lowest height. A heatsink whose microchannels have waves in their longitudinal dimensions is shown, but the expected effects of better mixing and greater contact area for the heat flux were not observed for the boundary conditions of the original experiment, the calculated thermal resistance was of 0.102 °C/W. The last analysis presents an idealization where the heat flux is evenly distributed in the silicon volume, with a thermal resistance of 0.084 °C/W.

White dwarf stars are the most common final stage of stellar evolution, corresponding to 99% of all stars in the Galaxy. White dwarf models can be used to obtain the age of stellar populations, to build an initial to final mass relation to understand the connection between the properties of white dwarfs and their progenitors, determine the upper mass limit that separates white dwarfs progenitors from Type II supernovae, enhance the comprehension of the physical properties of high density matter and derive ages and masses for observed white dwarfs from the cooling tracks. The literature is populated with low mass and intermediate mass white dwarf models, however the massive white dwarfs are often forgotten and the evolutionary sequences are incomplete. In this dissertation, we compute full evolutionary sequences for massive white dwarfs, exploring the evolution of hydrogen-rich and hydrogen-deficient white dwarfs stars with masses between 1.012 and 1.307 M , and initial metallicity of Z = 0.02. These sequences are the result of main sequence stars with masses between 8.8 and 11.8 M . The simulations were performed with the Modules for Experiments in Stellar Astrophysics - MESA code, starting at the zero-age main sequence, through thermally pulsing and mass-loss phases, ending as the white dwarfs at the cooling sequence. Our simulations are full evolutionary, in which we consider the entire evolutionary history of the progenitors. We present reliable nuclear chemical profiles for the whole mass range considered, covering the different expected central compositions, i.e. C/O, O/Ne and Ne/O/Mg, and their dependence with the stellar mass. In addition, we present detailed chemical profiles of hybrid C/O-O/Ne core white dwarfs, found in the mass range between 1.024 and 1.150 M . We present the initial-to-final mass relation, mass-radius relation, and cooling times with improved crystallization limits, considering the effects of atmosphere and core composition.