No presente trabalho estudam-se diferentes técnicas explícitas, em diferenças finitas, no emprego de algoritmos do tipo velocidade-pressão para simulação de escoamentos incompressíveis. O método de resolução da pressão de maneira direta e explícita, introduzido por Bravo e Claeyssen [BRA 97a], é analisado. Faz-se uma aproximação para o erro causado por esta técnica, e verifica-se como isto afeta a equação da continuidade. As simulações são realizadas na cavidade quadrada, comparando-se os diferentes métodos e validando as aproximações realizadas no estudo do método de resolução da pressão. Além disso, simula-se o escoamento em cavidades profundas e rasas, observandose a formação de vórtices e distribuição de energia cinética. Simulações do escoamento na cavidade cúbica também são apresentadas. / In this work different explicit technics in finite defferences in the application of velocity-pressure algorithm to simulate incompressible flows have been studied. The direct and explicit method of pressure resolution, introduced by Bravo anel Claeyssen [BRA 97a] is analyzed. An approximation to the error caused by this method is made, anel how this affects the continuity equation is verified. The simulations are maele in a square cavity, comparing the differents methoels anel valielating the approximations maele in the study of the pressure resolution method. Besieles this, flow in eleep and shallow cavities is simulateel, observing the formation of vortices and kinetic energy distribution. Simulations of the flow in the cubical cavity are also considered.

Submitted by Liliane Ferreira (ljuvencia30@gmail.com) on 2018-09-13T14:03:50Z No. of bitstreams: 2 Dissertação - Jaqueline Rodrigues Ferreira - 2018.pdf: 3324369 bytes, checksum: 7bb3737ade608223e8320020b0b39180 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2018-09-14T10:23:57Z (GMT) No. of bitstreams: 2 Dissertação - Jaqueline Rodrigues Ferreira - 2018.pdf: 3324369 bytes, checksum: 7bb3737ade608223e8320020b0b39180 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2018-09-14T10:23:57Z (GMT). No. of bitstreams: 2 Dissertação - Jaqueline Rodrigues Ferreira - 2018.pdf: 3324369 bytes, checksum: 7bb3737ade608223e8320020b0b39180 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2018-06-29 / The use of steel pipe piles as a foundation element brings a new focus to the contributions of soil plug to the behavior of open-ended pipe piles. Therefore, considering the necessity to understand the plug formation, its influence on pipe pile load capacity and the limited of research developed in real scale models, this work aims to numerically simulate the process jacking open ended steel pipe piles in sands, seeking to investigate the soil properties and piles geometry interference on the jacking process. For the development of these simulations an explicit finite difference program, FLAC (Fast Lagrangian Analysis of Continua), was used. The methodology developed started with the finite difference mesh, jacking velocity and interface parameters (normal and shear modulus) calibration. This process was validated with calibration chamber tests results developed in other researches. Henceforth, parametric analysis involving the pipe geometry – internal diameter, wall thickness and tip angle – and the soil mechanical parameters – friction angle, density and elastic modulus – were conducted. As result, among all the factors analyzed, soil density and pile tip angle do not interfere in soil plug formation. Furthermore, the program selected for the simulations was able to model plugging, exhibiting similar behavior to calibration chamber tests results. / A crescente utilização de estacas metálicas tubulares de ponta aberta como elemento de fundação traz novo enfoque às contribuições do tampão de solo ao comportamento destes elementos de fundação. Assim, considerando a necessidade da compreensão do processo de formação deste tampão, sua influência na capacidade de carga das estacas e a pequena quantidade de pesquisas desenvolvidas com modelos em escala real, o objetivo deste trabalho é simular numericamente o processo de cravação estática de estacas metálicas tubulares de ponta aberta em areias, buscando estudar a interferência das propriedades do solo e da geometria das estacas na dinâmica do processo executivo. Para o desenvolvimento destas simulações utilizou-se um programa de diferenças finitas explicito, FLAC (Fast Lagrangian Analysis of Continua). A metodologia desenvolvida iniciou-se com a calibração da malha de diferenças finitas, da velocidade de cravação e dos parâmetros das interfaces – módulos normal e cisalhante. Sendo esta calibração validada com a utilização dos resultados de ensaios em câmara de calibração realizados em outras pesquisas. Posteriormente, desenvolveu-se uma análise paramétrica onde foram variáveis a geometria do problema – diâmetro interno, espessura da parede e inclinação da ponta das estacas, e os parâmetros mecânicos do solo – ângulo de atrito, peso específico e módulo de elasticidade. Como resultado, têm-se que, o peso específico do solo e a inclinação da ponta das estacas não interferem na formação do tampão de solo, enquanto que o diâmetro interno e a espessura da parede das estacas apresentam grande influência nos resultados. Além disso, o programa selecionado para as simulações foi capaz de modelar o fenômeno do tamponamento, exibindo comportamento semelhante ao encontrado nos ensaios realizados em câmara de calibração.

Para o projeto de estruturas com perfis de aço formados a frio, é fundamental a compreensão dos fenômenos da instabilidade local e global, uma vez que estes apresentam alta esbeltez e baixa rigidez à torção. A determinação do carregamento crítico e a identificação do modo de instabilidade contribuem para o entendimento do comportamento dessas estruturas. Este trabalho avalia três metodologias para a análise linear de estabilidade de perfis de aço formados a frio isolados, com o objetivo de determinar os carregamentos críticos elásticos de bifurcação e os modos de instabilidade associados. Estritamente, analisa-se perfis de seção U enrijecido e Z enrijecido isolados, de diversos comprimentos e diferentes condições de vinculação e carregamento. Determinam-se os carregamentos críticos elásticos de bifurcação e os modos de instabilidade globais e locais por meio de: (i) análise com o Método das Faixas Finitas (MFF), através do uso do programa computacional CUFSM; (ii) análise com elementos finitos de barra baseados na Teoria Generalizada de Vigas (MEF-GBT), via uso do programa GBTUL; e (iii) análise com elementos finitos de casca (MEF-cascas) por meio do uso do programa ABAQUS. Algumas restrições e ressalvas com relação ao uso do MFF são apresentadas, assim como limitações da Teoria Generalizada de Viga e precauções a serem tomadas nos modelos de cascas. Analisa-se também a influência do grau de discretização da seção transversal. No entanto, não é feita avaliação em relação aos procedimentos normativos e tampouco análises não lineares, considerando as imperfeições geométricas iniciais, tensões residuais e o comportamento elastoplástico do material. / For the design of cold formed steel members, it is essential to understand the effects of local and global instability, since these members typically have a high slenderness and low torsion stiffness. The determination of critical loads and the associated buckling modes contribute to understand the behavior of these members. This work performs a evaluation of three methods for linear stability analysis of isolated cold-formed steel members in order to determine the elastic critical loads and the corresponding buckling modes. Specifically, Ue and Ze shape members were studied with various length, different boundary conditions and loads. The elastic critical loads and buckling modes are determined by means of: (i) analysis with the Finite Strip Method (FSM), by the computer program CUFSM, (ii) beam finite element analysis based on the Generalized Beam Theory (FEM-GBT), by GBTUL program, and (iii) Finite Element Method with shell analysis using ABAQUS program. Some restrictions and warnings regarding the use of the FSM are presented, as well as limitations of the Generalized Beam Theory and precautions to be taken in the shell models. It is also analyzed the influence of the degree of discretization of the cross section. In the present study, no evaluation was made with respect to normative procedures neither nonlinear analyses considering the initial geometric imperfections, residual stresses and elastoplastic behavior of the material.

As redes ópticas de comunicação têm possibilitado, cada vez mais, o oferecimento de serviços do tipo faixa larga. A rede de comunicação totalmente óptica está se tornando a meta tecnológica mais ambiciosa. Grandes esforços têm sido concentrados no desenvolvimento e aperfeiçoamento de materiais e dispositivos não-recíprocos, como isoladores e circuladores constituem uma importante classe de dispositivos ópticos. Os isoladores são utilizados em sistemas ópticos para prevenir a reflexão de luz em lasers e amplificadores. Os circuladores são empregados em esquemas de derivação de sinais que utilizam multiplexação em comprimento de onda (WDM). O funcionamento destes dispositivos é baseado nas propriedades de materiais magnetoópticos. As metas desta dissertação são apresentar as principais características dos materiais magnetoópticos, explorando as características não-recíprocas dos modos TM. Guias planares e tridimensionais são analisados neste trabalho. Para tanto, são obtidas as expressões das componentes dos campos e as equações características dos modos de interesse em estruturas planares por meio da técnica de matriz de transferência (TMT). A análise de propagação de onda em guias planares com materiais magnetoópticos é feita por meio de simulação numérica empregando o método de propagação de feixe (BPM) baseado em diferenças finitas (FD) e o esquema de Crank-Nicholson (CN) na discretização da solução da equação de onda. A condição de fronteira transparente (TBC) é incorporada ao algoritmo FD-BPM com a finalidade de se evitar reflexões de ondas eletromagnéticas para dentro da janela computacional. O método do índice efetivo é empregado na análise de guias de onda tridimensionais do tipo rib. / Optical communication networks have allowed a continuous increase of broadband services offer. The all-optical communication networks are becoming the most ambitious technological goal. Great efforts have been concentrated on the materiaIs and devices development and improvement to make it possible. Nonreciprocal devices, such as isolators and circulators constitute an important class of optical devices. Isolators are used in optical systems to avoid reflection of light in lasers and amplifiers. Circulators are used in signal derivation schemes that use wavelength division multiplexing (WDM). The operation of these devices is based on the properties of magnetooptic materiaIs. The purposes of this dissertation are to present the main features of the magnetooptic materiaIs as well as to analyze the eletromagnetic wave propagation in magnetooptic waveguides, exploring nonreciprocal features of TM modes. Planar and three-dimensional waveguides are analysed in this present study. Therefore expressions of electromagnetic field components and characteristic equations of the modes of interest in planar structures are obtained by using transfer matrix technique (TMT). The wave propagation analysis in planar magnetooptic waveguides is realized by using the finite-difference beam propagation method (FD-BPM) and Crank-Nicholson scheme (CN) applied to wave equation solution discretization. In order to avoid electromagnetic wave reflection into computational window, the transparent boundary condition (TBC) is incorporated to the FD-BPM formalism. The effective index method (EIM) is used in the analysis of three-dimensional rib magnetooptic waveguides.

Neste trabalho apresenta-se uma combinação de formulações numéricas para a análise da interação estaca-solo com ou sem blocos de capeamento rígido, sujeita à carga horizontal e vertical. Nestas formulações as estacas são representadas pelo método das diferenças finitas (MDF) ou pelo método dos elementos finitos (MEF) e o solo é representado pelo método dos elementos de contorno (MEC). Na utilização do MEF, para a análise das estacas, os deslocamentos e as forças de interação foram representados por várias funções polinomiais chegando-se a um elemento finito final considerado eficiente e constituído por quatro pontos nodais, 14 parâmetros nodais, sendo quatro para deslocamentos lineares em cada uma das direções (X1, X2 e X3) e mais dois parâmetros referentes as rotações do topo da estaca em torno dos eixos X1 e X2. Os deslocamentos transversais ao longo da estaca foram representados por uma função polinomial do 4º grau e os deslocamentos axiais foram representados por uma função cúbica. Para as forças da interface nas direções X1 e X2 são utilizados funções polinomiais cúbicas. As forças de superfície cisalhantes que ocorrem ao longo do fuste da estaca são representadas por um polinômio quadrático e a tensão normal à seção da extremidade inferior da estaca é suposta constante. O maciço de solos é modelado pelo MEC como um meio contínuo, elástico-linear, semi-infinito, isótropo e homogêneo. Combinando-se estes métodos de análise, obtém-se um sistema de equações lineares representando o problema de interação estaca-solo. Após a resolução deste sistema, são obtidos os deslocamentos e rotações nos nós do elemento e as tensões de contato estaca-solo. Vários exemplos envolvendo as formulações propostas são analisados e os resultados obtidos são concordantes com os de outros autores. / This work presents a mixed numerical formulation of analysis of the pile-soil interaction, with or without rigid caps, under horizontal and vertical loads. In these formulations the piles are modeled by finite difference method (FDM) or finite element method (FEM) and the soil is represented by the boundary element method (BEM). In the finite element method for the pile analysis the displacements and interactions forces with the soil were represented by several polynomials functions leading to an efficient element with four nodal points, 14 nodal parameters, where four parameters are the linear displacements for each directions (X1, X2 e X3) and two parameters for the rotations at the top of the pile. The transversal displacements along the pile were represented by a fourth degree polynomial and the axial displacements were represented by cubic functions. The interface forces in the X1 and X2 directions are represented by cubic functions. The shear contact forces along of the pile surface are approximated by a second degree polynomial and the normal tractions in the pile tip is assumed to be constant over its crosssection. The soil is modeled by BEM as an isotropic, homogeneous, semi-infinite and linear-elastic continuum. Combining these analysis methods a complete system representing a pile-soil interaction problem can be obtained and from it the displacements and the pile-soil contact tractions can be achieved. Various examples are presented and the results closely agree with others authors.

Este trabalho tem como objetivo desenvolver um método numérico para simular escoamentos incompressíveis, isotérmicos, confinados ou com superfícies livres, de fuidos viscoelásticos governados pelos modelos integrais de Maxwell e K-BKZ (Kaye-Bernstein, Kearsley e Zapas). A técnica numérica apresentada é uma extensão do método GENSMAC (Tomé McKee - J. Comp. Phys., (110), pp 171--186, 1994 ) para a solução das equações de conservação, juntamente com as equações constitutivas integrais de Maxwell e K-BKZ. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada. O tensor de Finger, B_t\'(t) é calculado com base nas idéias do método de campos de deformação (Peters et al. - J. Non-Newtonian Fluid Mech. (89), de maneira que não há a necessidade de seguir a trajetória da partícula de fuido para descrever a história de deformação da partícula. Uma abordagem diferente para a discretização do instante passado é utilizada e o tensor de Finger e o tensor das tensões são calculados utilizando um método de segunda ordem. A validação do método numérico descrito nesse trabalho foi feita utilizando o escoamento em um canal bidimensional e a solução numérica obtida para a velocidade e para as componentes de tensão com o modelo de Maxwell foram comparadas com as respectivas soluções analíticas no estado estacionário, mostrando excelente concordância. Os resultados numéricos para a simulação do escoamento em uma contração planar 4 : 1 mostraram bons resultados, tanto qualitativos quanto quantitativos, quando comparados com os resultados experimentais de Quinzani et al. ( J. Non-Newtonian Fluid Mech. (52), pp 1?36, 1994 ). Além disso, utilizando os modelos Maxwel e K-BKZ, o escoamento em uma contração planar 4 : 1 foi simulado para vários números de Weissenberg e os resultados obtidos estão de acordo com os encontrados na literatura. Resultados numéricos de escoamentos com superfícies livres modelados pelas equações integrais de Maxwell e K-BKZ são apresentados. Em particular, a simulação numérica do jato oscilante para diferentes números de Weissenberg e diferentes números de Reynolds é apresentada. / The aim of this work is to develop a numerical technique for simulating incompressible, isothermal, free surface (also con¯ned) viscoelastic flows of fuids governed by the integral models of Maxwell and K-BKZ (Kaye-Bernstein, Kearsley and Zapas). The numerical technique described herein is an extension of the GENSMAC method (Tome and McKee, J. Comput. Phys., 110, pp. 171-186, 1994) to the solution of the momentuum and mass conservation equations together with the integral constitutive Maxwell and K-BKZ equations. The governing equations are solved by the finite difference method on a staggered grid using a Marker-and-Cell approach. The fluid is represented by marker particles on the fluid surface only. This provides the visualization and location of the fluid free surface so that the free surface stress conditions can be applied. The Finger tensor Bt0(t) is computed using the ideias of the deformation fields method (Peters et al. J. Non-Newtonian Fluid Mech., 89, pp. 209-228, 2001) so that it is not necessary to track a fluid particle in order to calculate its deformation history. However, in this work modifcations to the deformation fields method are introduced: the past time is discretized using a different formula, the Finger tensor Bt0(x; t) is obtained by a second order method and the stress tensor ? (x; t) is computed by a second order quadrature formula. The numerical method presented in this work is validated by simulating the flow of a Maxwell fluid in a two-dimensional channel and the numerical solutions of the velocity and the stress components are compared with the respective analytic solutions providing a good agreement. Further, the flow through a 4:1 planar contraction of a specific fuid studied experimentally by Quinzani et al. (J. Non-Newtonian Fluid Mech., 52, pp. 1-36, 1994) was simulated and the numerical results were compared qualitatively and quantitatively with the experimental results and very good agreement was obtained. The Maxwell and the K-BKZ models were applied to simulate the 4:1 planar contraction problem using various Weissenberg numbers and the numerical results were in agreement with those published in the literature. Finally, numerical results of free surface flows using the Maxwell and K-BKZ integral constitutive equations are presented. In particular, the numerical simulation of jet buckling using several Weissenberg numbers and various Reynolds numbers are presented

No presente trabalho estudam-se diferentes técnicas explícitas, em diferenças finitas, no emprego de algoritmos do tipo velocidade-pressão para simulação de escoamentos incompressíveis. O método de resolução da pressão de maneira direta e explícita, introduzido por Bravo e Claeyssen [BRA 97a], é analisado. Faz-se uma aproximação para o erro causado por esta técnica, e verifica-se como isto afeta a equação da continuidade. As simulações são realizadas na cavidade quadrada, comparando-se os diferentes métodos e validando as aproximações realizadas no estudo do método de resolução da pressão. Além disso, simula-se o escoamento em cavidades profundas e rasas, observandose a formação de vórtices e distribuição de energia cinética. Simulações do escoamento na cavidade cúbica também são apresentadas. / In this work different explicit technics in finite defferences in the application of velocity-pressure algorithm to simulate incompressible flows have been studied. The direct and explicit method of pressure resolution, introduced by Bravo anel Claeyssen [BRA 97a] is analyzed. An approximation to the error caused by this method is made, anel how this affects the continuity equation is verified. The simulations are maele in a square cavity, comparing the differents methoels anel valielating the approximations maele in the study of the pressure resolution method. Besieles this, flow in eleep and shallow cavities is simulateel, observing the formation of vortices and kinetic energy distribution. Simulations of the flow in the cubical cavity are also considered.

Made available in DSpace on 2016-12-23T14:36:45Z (GMT). No. of bitstreams: 1 dissertacao.pdf: 1896132 bytes, checksum: dc313d94261c073031be0aad2e3bffbf (MD5) Previous issue date: 2005-06-30 / Discute a implementação do algoritmo numérico para simulação de escoamento de fluidos incompressíveis, baseado no método de diferenças finitas, projetado para plataformas de processamento paralelo com memória distribuída, particularmente para clusters de estações de trabalho. O algoritmo de solução para as equações de Navier-Stokes utiliza um esquema explicito para pressão e um esquema implícito para as velocidades. A implementação paralela é baseada na decomposição do domínio, onde o domínio computacional do problema é decomposto em vários blocos, sendo um ou mais destinados a nós de processamento distintos. Todos os nós então processam em paralelo as tarefas de computação sobre os blocos a eles designados. O processamento paralelo inclui inicialização, cálculo de coeficientes, solução linear nos subdomínios, e comunicação entre os nós. A troca de informação entre os processos referentes a cada subdomínio é realizada utilizando a biblioteca message passing interface (MPI), o que assegura portabilidade entre diferentes plataformas computacionais, abrangendo desde máquinas maciçamente paralelas (MPP) até clusters de estações de trabalho. Para melhorar os níveis de desempenho obtidos pelo algoritmo, foram investigadas técnicas para a redução do volume de comunicação entre processadores e utilização mais eficiente da memória cache dos microprocessadores. Para avaliar o desempenho do algoritmo desenvolvido e analisar as diferentes estratégias de paralelização foram executadas simulações com cluster de 2 a 56 processadores, nas quais foram avaliados o tempo de execução, speedup e eficiência paralela. Os resultados experimentais mostram que as otimizações relacionadas aos fatores de comunicação melhoram o speedup em até 165%, e a técnica de utilização mais eficiente da memória cache pode melhorar o speedup em mais 40% acima da otimização da comunicação. / This work discusses the implementation of a numerical algorithm for simulating incompressible fluid flows, based on the finite difference method, and designed for parallel computing platforms with distributed-memory, particularly for clusters of workstations. The solution algorithm for the Navier-Stokes equations utilizes an explicit scheme for pressure and an implicit scheme for velocities. The parallel implementation is based on domain decomposition, where the original calculation domain is decomposed into several blocks, each of which given to a separate processing node. All nodes then execute computations in parallel, each node on its associated sub-domain. The parallel computations include initialization, coefficient generation, linear solution on the sub-domain, and inter-node communication. The exchange of information across the sub-domains, or processors, is achieved using the message passing interface standard, MPI. The use of MPI ensures portability across different computing platforms ranging from massively parallel machines to clusters of workstations. Three different optimization strategies were evaluated in order to improve the computational performance of the algorithm, which include techniques exploring a reduction in the communication volume between processors and a more efficient utilization of the microprocessor s cache memory. In order to evaluate the performance levels obtained, and to analyze the effectiveness of the optimization strategies adopted, simulations using a 64 nodes cluster were executed. The simulations were performed using 2 to 56 processors, where execution time and speed-up were measured. The results indicate that the optimizations related to communication factors can improve the speed-up obtained up to 165%, while the cache memory optimization technique used can improve the speed-up obtained in further 40%.

Orientador: Hugo Enrique Hernandez Figueroa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-05T08:10:41Z (GMT). No. of bitstreams: 1 Santos_CarlosHenriquedaSilva_M.pdf: 1752834 bytes, checksum: 8ed5b0902bb130762ff802db03187fbb (MD5) Previous issue date: 2005 / Resumo: Esse trabalho tem por objetivo desenvolver soluções computacionais de alto desempenho a um baixo custo, seguindo as propostas incentivadoras do Governo Federal para adoção de software livre. Essas soluções possibilitam simular, de maneira eficiente, os domínios computacionais de médio e grande porte utilizados no eletromagnetismo computacional. Os bons resultados obtidos nesse trabalho mostram a importância e eficiência da computação massivamente paralela utilizando cluster Beowulf para o processamento do método FDTD aplicado em estruturas complexas, porém a um baixo custo financeiro. O desempenho desse sistema ficou comprovado na realização de experimentos para analisar a SAR na cabeça humana e estudar os efeitos de estruturas metamateriais / Abstract: This work has as objective to develop high performance computational solutions to a low cost, following the stimulated proposals of the Federal Government for adoption of free software. They make possible to simulate, in efficient way, the computational domains of middle and high size useful on the computational electromagnetism. The good results gotten in these work showed the importance and efficiency of the massive parallel computation using the Beowulf cluster for the process the FDTD method applied on complex structures, however to a low financial cost. The performance of this system was proved in the realization of experiments to analyze the SAR on the human head and to study the effects of metamarial structures / Mestrado / Telecomunicações e Telemática / Mestre em Engenharia Elétrica

Submitted by Erika Demachki (erikademachki@gmail.com) on 2016-10-19T16:18:12Z No. of bitstreams: 2 Dissertação - Felipe Resende de Carvalho Sousa - 2016.pdf: 6672203 bytes, checksum: b0d857b45451fe225dfdbd424c35be00 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Erika Demachki (erikademachki@gmail.com) on 2016-10-19T16:18:40Z (GMT) No. of bitstreams: 2 Dissertação - Felipe Resende de Carvalho Sousa - 2016.pdf: 6672203 bytes, checksum: b0d857b45451fe225dfdbd424c35be00 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2016-10-19T16:18:40Z (GMT). No. of bitstreams: 2 Dissertação - Felipe Resende de Carvalho Sousa - 2016.pdf: 6672203 bytes, checksum: b0d857b45451fe225dfdbd424c35be00 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2016-06-10 / Prediction ofmoisture levels inside power transformers shown itselfas an important technique for maintenance engineering, because in addition to being associated with the useful life of equipment, it also provides arguments for decision making about intervention in them. Despite concerns about water in transformers existed for over half a century, the monitoring of this variable is still done in a simply way, given little information to weave a diagnosis, which can lead to superficial interpretations of equipment, and therefore, the misdiagnosis. Given such a scenario, this work aims to obtain an equation that models classical curves relative to moisture in transformers, as well as the development of a mathematical model that simulates water migration inside these equipments. For this purpose, computational techniques and numerich methods are used. Thus, the main contribution of this work is to optimize the above diagnostic process. During the work a literature review of the subject is presented, with further explanations of the mathematical tools used. To prove the effectiveness of the project presented, case studies with data from energized power transformers are analysed. Some values found in other works are also uses to validate the work. Ending with some conclusions, this paper proves that the use of the approach presented here is effective, and its application is fundamental in the area of maintenance engineering. / A predição dos níveis de umidade no interior de transformadores de potência se mostra como uma técnica importante para a Engenharia de Manutenção, pois, além de estar associada à vida útil destes equipamentos, fornece argumentos para a tomada de decisão quanto a intervenções nos mesmos. Apesar da preocupação com a presença de água em transformadores existir há mais de meio século, o acompanhamento de tal variável ainda é realizado de maneira superficial, considerando-se poucas informações para se apresentar um diagnóstico eficiente, o que pode levar a interpretações equivocadas sobre o equipamento e, consequentemente, a falsos diagnósticos. Tendo em vista tal cenário, este trabalho possui como objetivo a obtenção de um equacionamento que modele curvas clássicas relativas à umidade em transformadores, bem como o desenvolvimento de um modelo matemático que simule o comportamento da migração de água no interior destes equipamentos. Para isso, técnicas computacionais e métodos numé- ricos são utilizados. Desta forma, a principal contribuição do trabalho é otimizar o processo de diagnóstico supracitado. Durante o trabalho, é realizada uma análise bibliográfica do tema, com posteriores explicações sobre as ferramentas matemáticas utilizadas. Para comprovar a eficácia do projeto apresentado, são realizados estudos de caso com dados obtidos de transformadores energizados. Alguns valores encontrados em outros trabalhos também são utilizados para validar o trabalho. Finalizando com conclusões sobre o exposto, comprova-se que a utilização da abordagem aqui apresentada é eficaz e de fundamental aplicação na área de Engenharia de Manutenção.