O objetivo desta tese é introduzir uma nova metodologia para a cinética bidimensional multi- grupo de difusão de nêutrons em reatores nucleares. A metodologia apresentada usa uma aproximação polinomial em um domínio homogêneo retangular com condições de contornos não homogêneas. Como ela consiste em uma série de Taylor truncada, sua estimativa de erro varia de acordo com o tamanho do retângulo. Os coeficientes são obtidos principalmente pelas suas relações com o termo independente, que _e determinado pela equação diferencial. Estas relações são obtidas apenas pelas condições de contorno, e é demonstrado serem linearmente independentes. Um esquema numérico é feito para assegurar uma rápida convergência. Estes procedimentos feitos para um retângulo homogêneo são feitos para construir soluções para problemas de autovalor e dependentes do tempo de geometria ortogonal global com parâmetros seccionalmente constantes pelo método iterativo SOR. O autovalor dominante e sua autofunção são obtidos pelo método da potência no problema de autovalor. A solução para casos dependentes do tempo usam o método de Euler modificado na variável tempo. Quatro casos-teste clássicos são considerados para ilustração. / The objective of this thesis is to introduce a new methodology for two{dimensional multi{ group neutron diffusion kinetics in a reactor core. The presented methodology uses a polyno- mial approximation in a rectangular homogeneous domain with non{homogeneous boundary conditions. As it consists on a truncated Taylor series, its error estimates varies with the size of the rectangle. The coefficients are obtained mainly by their relations with the independent term, which is determined by the differential equation. These relations are obtained by the boundary conditions only, and these relations are proven linear independent. A numerical scheme is made to assure faster convergence. The procedures done for one homogeneous rectangle are used to construct the solution of global orthogonal geometry with step{wise constant parameters steady state and time dependent problems by the iterative SOR algo- rithm. The dominant eigenvalue and its eigenfunction are obtained by the power method in the eigenvalue problem. The solution for the time dependent cases uses the modi ed Euler method in the time variable. Four classic test cases are considered for illustration.

Taludes escavados em filitos muitas vezes estão sujeitos à instabilidade, principalmente quando há outras famílias de descontinuidades presentes na estrutura. Este problema está presente na mina de Candiota, onde foi identificado, em um dos taludes, potencialidade para ruptura por tombamento. Este trabalho aborda uma discussão sobre métodos analíticos e numéricos aplicado a tombamento primário, com objetivo de avaliar o fator de segurança do talude em estudo por meio de método numérico por elementos finitos, utilizando o Phase2 da Rocscience, e métodos analíticos. Também são apresentadas doze simulações para avaliar se o movimento de tombamento em modelagens numéricas obedece a um efeito de escala. Com uso de métodos numéricos dois modelos foram construídos: um contínuo equivalente e outro contínuo (mas com as descontinuidades inseridas no modelo). No modelo com as descontinuidades a rocha foi tratada como elástica por meio do critério de Hoek-Brown e as descontinuidades foram tratadas como elasto-plásticas por meio do critério de Mohr-Coulomb. Os resultados das modelagens da mina de Candiota foram semelhantes e mostraram fatores de segurança em níveis de estabilidade, sendo que no modelo contínuo (mas com descontinuidades inseridas) se observou de forma mais clara as tensões cisalhantes induzidas nas descontinuidades quando ocorre o movimento de tombamento. O uso dos métodos analíticos mostrou-se ineficientes para o problema apresentado. Os resultados das doze modelagens sugerem que o aumento do espaçamento das descontinuidades influencia nos fatores de resistência dos taludes e na forma do movimento de tombamento. / Slopes in phyllites are often subject to instability, especially when there are other families of discontinuities in the structure. This problem is present in Candiota mine, which was identified in one of the slopes potential to toppling. This work show a discussion of analytical and numerical methods applied to toppling, to evaluate the slope safety factor being studied by numerical method finite element using the Phase2 of Rocscience, and analytical methods. Also shown are twelve simulations to evaluate if the toppling in numerical modeling follows a scaling effect. With use of numerical methods two models were built: an equivalent continuous and another solid (but with discontinuities inserted in the model). In the model with the discontinuities the rock was treated as elastic by the Hoek-Brown criterion and discontinuities were treated as elastic-plastic by the Mohr-Coulomb criterion. The results of Candiota mine were similar and showed safety factors in stability levels, and in the continuous model (with discontinuities) was observed more clearly the shear stress induced in discontinuities when the movement of overturning occurs. The use of analytical methods proved inefficient for the problem presented. The results of the twelve modeling suggests increasing the spacing of the discontinuities of slope influences the resistance factor and the form of the toppling movement.

Nesta tese desenvolvemos um método numérico para diagonalizar o Hamiltoniano de Heisenberg iterativamente. O método consiste basicamente em diagonalizar cadeias de spins, cada vez maiores, em que cada passo da diagonalização corresponde à adição de um novo spin à cadeia. A base de vetores para calcular o Hamiltoniano de uma cadeia de N spins, HN, é construída por meio do produto direto dos autovetores do Hamiltoniano Hn-1 da rede diagonalizada no passo anterior, pelos autoestados correspondentes ao N-ésimo spin adicionado. Além de usar a comutação do Hamiltoniano com a componente azimutal do spin total, Sz, prática comum em outros métodos, usufruímos da conservação com o quadrado do spin total, S2. Para uma classe específica de redes também implementamos a simetria de reflexão. Obtemos o espectro completo de energia de cadeias de spins 1/2 com até 20 sítios, para as quais mostramos resultados da dependência com a temperatura da susceptibilidade magnética e do calor específico, para redes com impurezas tipo spin substitucionais, com defeitos nas ligações ou com efeitos de bordas, isto é, para sistemas sem invariância translacional. Usualmente essa restrição impõe enormes dificuldades em métodos tradicionais. Para diagonalizar cadeias com um número maior de sítios, implementamos um procedimento que seleciona os estados de mais baixa energia para serem usados na base de vetores do passo seguinte. Com esse tipo de truncamento de estados, fomos capazes de obter o estado fundamental e alguns estados de baixa energia de cadeias com mais de uma centena de sítios, com precisão de até cinco algarismos significativos. Nossos resultados reproduzem os da literatura para os casos conhecidos, em geral sistemas homogêneos. As aproximações desenvolvidas recentemente no contexto da Teoria do Funcional da Densidade, aplicada ao modelo de Heisenberg, e que também se aplicam a sistemas inomogêneos, estão em conformidade com nossos resultados numericamente exatos. Generalizamos o método para diagonalizar escadas de spins 1/2. Calculamos o estado fundamental e o gap de energia desse sistema, onde variamos a razão entre os acoplamentos ao longo das pernas da escada e ao longo dos degraus da mesma; nossos resultados são comparados com os da literatura. Apresentamos também a implementação do método iterativo no modelo de Hubbard, que descreve um sistema de spins itinerantes. Sabe-se que no regime de alta repulsão Coulombiana entre os spins e densidade um (número de spins igual ao número de sítios da cadeia), esse modelo é mapeado no modelo de Heisenberg, resultado que é verificado numericamente em nosso procedimento por meio do cálculo de energias de ambos os modelos em um regime paramétrico apropriado. / In this Thesis we develop a numerical method to diagonalize the Heisenberg model iteratively. In essence, we diagonalize spin chains in steps, each one corresponding to an addition of a spin to a smaller chain. The basis vectors to calculate the Hamiltonian of a N-spin chain, HN, is built by means of the direct product of the eigenvectors of the (N-1)-spin Hamiltonian, diagonalized on the previous step, by the eigenstates of the N-th added spin. Besides the common use of the conservation of the z-component of the total spin, Sz, we also exploit the conservation of the squared total spin, S2. For a specific class of spin systems we also implemented the reflection symmetry. We obtain the entire energy spectrum of spin-1/2 chains up to 20 sites, for which we show the temperature dependence of the magnetic susceptibility and specific heat, for systems with substitutional impurity spins, bond defects, border effects, i.e., for systems without translational invariance. This normally imposes enormous restrictions in many traditional methods. In order to diagonalize chains with a larger number of sites we implemented a procedure that selects lower energy states to be used in the basis vector on the next step. Using this truncation scheme, we are able to obtain low-lying energy states for chains with more than a hundred sites, up to five significant figures of accuracy. Our results reproduce those of the literature for the known cases, in general homogeneous systems. The approaches recently developed in the context of Density Functional Theory to the Heisenberg model, which also apply to inhomogeneous systems, are consistent with our numerical results. We generalize the method to diagonalize spin-1/2 ladders. We calculate the ground-state and the energy gap of this system, for arbitrary ratio of the couplings along the lags or over the rungs of the ladder. We also present the implementation of our iterative method to the Hubbard model, which describes a system of itinerant spins. It is known that in the regime of high Coulomb repulsion between the spins and unitary density (number of spins equal to the number of sites in the chain), this model is mapped onto Heisenberg one, a result which is verified numerically in our procedure by calculating the energy spectrum of both models in na appropriated parametric regime.

nÃo hà / Para determinar o potencial eÃlico de uma regiÃo à de fundamental importÃncia que seja realizado um estudo sobre as caracterÃsticas intrÃnsecas do vento do local. Este trabalho aborda alguns mÃtodos numÃricos a serem empregados no cÃlculo dos parÃmetros da distribuiÃÃo de Weibull que auxilia no estudo da velocidade do vento, para que haja uma correta definiÃÃo das condiÃÃes naturais existentes. O prÃvio conhecimento destas informaÃÃes coopera no processo de tomada de decisÃo sobre a viabilidade tÃcnica na instalaÃÃo de novos parques eÃlicos industriais. à realizada uma anÃlise estatÃstica entre sete sistemas matemÃticos conhecidos da literatura para estimar os parÃmetros ( k ) de forma e ( c ) de escala da curva de distribuiÃÃo de frequÃncias de Weibull. SÃo utilizados dados de velocidade do vento e de potÃncia eÃlica de duas cidades litorÃneas do Estado do Cearà pertencentes à regiÃo Nordeste do Brasil, Icapuà e Camocim. Os mÃtodos apurados no desenvolvimento desta pesquisa sÃo: MÃtodo GrÃfico, MÃtodo da MÃxima VerossimilhanÃa, MÃtodo da MÃxima VerossimilhanÃa Modificado, MÃtodo EmpÃrico, MÃtodo do Momento, MÃtodo da Energia PadrÃo e MÃtodo da Energia Equivalente. A realizaÃÃo da anÃlise comparativa de eficiÃncia e exatidÃo entre estes, compreende a aplicaÃÃo dos seguintes testes estatÃsticos: AnÃlise de VariÃncia (R2 ) , Raiz Quadrada dos Erros QuadrÃticos MÃdios (RMSE) e Teste do Qui-quadrado (X 2 ) . / To determine the wind potential of a region is of paramount importance that a study be conducted on the intrinsic characteristics of the wind site. This paper address some numerical methods to be used in calculating the parameters of the Weibull distribution the aids in the study of wind speed, so there is a correct definition of natural conditions existing. The previous knowledge of this information assist in the decision-making process on the technical feasibility of installing industrial wind farms. It perform a statistical analysis of seven mathematical systems known in the literature for estimating the parameters ( k ) form and ( c ) scale of the frequency distribution curve of Weibull. Data are used for wind speed and wind power from two coastal cities of Cearà State belonging to Northeast Brazil, Icapuà and Camocim. The methods to be employed in the development of this research are: Graphical Method, Maximum Likelihood Method, Maximum Likelihood Modified Method, Empirical Method, Moment Method, Energy Pattern Factor Method and the Equivalent Energy Method. The realization of the comparative analysis of efficiency and accuracy among these include the application of the following statistical tests: ANOVA(R2 ) , Square Root of Average Quadratic Errors (RMSE) and Chi-square (X 2 ).

Este trabalho estuda a técnica de reconstrução de imagens conhecido como tomografia por impedância elétrica em um domínio bidimensional. Esta técnica consiste na alocação de eletrodos na fronteira do volume e uma fonte injeta padrões de corrente através dos eletrodos e medem-se as voltagens resultantes na fronteira. Com estes dados estima-se a condutividade (ou resistividade) do interior do domínio criando-se uma imagem do mesmo. A tomografia por impedância elétrica é um problema inverso e mal posto no sentido de Hadamard. Estudam-se diversos métodos de solução para resolver o problema direto usando métodos numéricos como diferenças finitas e volumes finitos. Proporemos os métodos numéricos a serem aplicados na solução do problema direto os quais serão testados com problemas onde a solução analítica é conhecida. Posteriormente aplicaremos os métodos propostos ao problema especifico. Uma questão importante na reconstrução de imagens é propor a maneira como aproximar o Jacobiano (ou matriz de sensibilidade) do problema, assim desenvolvemos uma técnica para a aproximação do mesmo usando os dados fornecidos pelo problema direto. / In this work is studied the technique of reconstruction of images known as electrical impedance tomography for a two-dimensional domain. This technique consists in the allocation of electrodes on the border of the volume and a source injects patterns of current through the electrodes and then measuring voltages through the other electrodes. With these data it is estimated the conductivity (or resistivity) on the interior of the domain and an image is create of it. The electrical impedance tomography is an inverse and ill conditioned problem in the Hadamard sense. In this work, is studying some numerical methods to solve the direct problem and are applied numerical methods such as the finite difference method and the finite volume method. It is proposed some numerical methods to solve the direct problem which will be tested with analytical problems where the solution is known. Later, apply the methods proposed to the specific issue. An important issue in the reconstruction problems is about the Jacobian (or sensitivity matrix) aproximation, thus proposing a technique for the calculation of even using the data provided by the direct problem. Keywords:

Neste trabalho um método de partículas para simular a dinâmica de sólidos elásticos e interação fluido estrutura e implementado. O método e baseado no Moving Particle Semi-implicit (MPS), originalmente desenvolvido para escoamentos incompressíveis com superfície livre. A estratégia principal do MPS e substituir os operadores diferenciais das equações governantes por operadores diferenciais discretos em uma distribuição de nos irregulares, derivados de um modelo de interação entre partículas. Inicialmente são apresentados os detalhes da formulação do método e modelos constitutivos utilizados. Uma condição simplificada de fragmentação e proposta, assim como um algoritmo de detecção de contato, permitindo a fragmentação entre vários sólidos. No caso da interação fluido-estrutura, as partículas de superfície do solido são tratadas como partículas de fluido e as pressões destas partículas são calculadas pela resolução da equação de Poisson para a pressão, tal como as partículas de fluido. Desta forma, o acoplamento entre solido e fluido e realizado utilizando o deslocamento e velocidade do solido elástico, como condições de contorno do fluido, e a pressão na interface, obtida pela resolução do movimento do fluido, e aplicada ao movimento do solido elástico. São apresentados e detalhados os algoritmos de solido elástico, fragmentação, colisão e acoplamento fluido-estrutura. Validações qualitativas e quantitativas do método são realizadas para casos estáticos e dinâmicos sujeitos a diferentes condições de contorno, comparando os resultados numéricos obtidos pelo MPS, outros métodos numéricos, soluções analíticas e medições experimentais presentes na literatura. / In this work a particle method to simulate the dynamics of elastic solids and fluid-structure interaction is implemented. It is based on the Moving Particle Semi-implicit Method (MPS), which was originally developed for incompressible flows with free surface. The main strategy of the MPS is to replace the differential operators of the governing equations by discrete differential operators on irregular nodes, which are derived from a model of interaction between particles. Initially details of the method and constitutive equations are shown. A simplified condition of fragmentation and collision between solids are proposed to allow the investigation of fragmentation amount multiple solids. In case of fluid-structure interaction, the solid\'s surface particles are treated as a fluid particle and the pressures of the surface particles are computed by solving Poisson equation for the pressure, just as the fluid particles. Therefore, the coupling between solid and fluid is done by using the displacement and velocity of elastic solid as the boundary conditions of the fluid, and the pressure at the interface, which is obtained when solving the fluid motion, is used to calculate the motion of the elastic solid. The algorithms for elastic solid, fragmentation, collision and fluid-structure interaction are presented and detailed. The qualitative and quantitative validations of the method are carried out herein considering static and dynamic cases subjected to deferent boundary conditions by comparing the numerical results from MPS with other numerical, analytical and experimental results available in the literature.

O trabalho desenvolvido tem como objetivo aplicar metodologias de otimização de suspensão para veículos de passageiro e comerciais leves, baseando-se em parâmetros de dirigibilidade e conforto veicular, através do uso de simulação numérica computacional. São apresentadas métricas objetivas utilizadas para a avaliação de um veículo com relação ao seu desempenho em termos de dirigibilidade e conforto e é proposta uma nova métrica global conjunta. São desenvolvidos no trabalho modelos analíticos para quantificar as métricas de dirigibilidade. Os resultados obtidos foram comparados com medições experimentais e resultados de modelos multicorpos mais complexos, atingindo o nível de correlação necessário para os propósitos deste trabalho. Alguns dos modelos analíticos desenvolvidos são contribuições inovadoras deste trabalho, não tendo correspondente anterior na literatura. É apresentada a modelagem dinâmica vertical que possibilita avaliar as métricas de conforto, cuja aplicação foi feita em conjunto com os modelos analíticos de dirigibilidade desenvolvidos, obtendo-se uma avaliação global conjunta. Três ferramentas de otimização são apresentadas e avaliadas na aplicação ao problema de otimização global. A análise comparativa dos resultados dos métodos de otimização permite identificar qual método mais adequado com relação ao desempenho computacional, praticidade de uso e disponibilidade de informação. Finalmente, demonstra-se que a aplicação da otimização numérica proporciona resultados efetivos para melhoria do produto, trazendo ganhos de tempo e custo no desenvolvimento de um novo projeto. / This work is intended to apply suspension set-up optimization methodologies to passenger vehicles and pick-up trucks based on their ride and handling parameters through the use of numerical computational simulation. Metrics used to evaluate a vehicle regarding its performance in terms of ride and handling are shown and a new global single metric for ride and handling is proposed. Analytical models are developed to quantify the vehicle handling metrics. The results obtained were compared with experimental measurements and the results of more complex multibody models, achieving the correlation level required for the purposes of this work. Some of the analytical models herein developed are new contributions from this work and were not previously available in the literature. The vertical dynamic model that allows the computation of the ride comfort metrics is shown, and its application simultaneously with the analytical handling models developed allows the calculation of the proposed global single ride and handling metric. Three different optimization techniques are presented and studied in order to compare their performance for the proposed problem. The comparative analysis among the optimization results allows determining where each method is more adequate with respect to computational performance, usage friendliness and information availability. Finally, it is shown that the application of numerical optimization is effective to improve the product performance, with gains in terms of development time and cost for a new project.

Esta dissertação apresenta uma revisão histórica dos trabalhos acerca de superfícies minimais, ressaltando a pertinência da analogia entre a busca de superfícies de mínima área e a busca de formas de membranas estruturais sujeitas a um estado de tensões superficiais, homogêneo e isótropo. São colocados alguns conceitos geométricos das superfícies parametrizáveis, com base na geometria diferencial, a fim de realizar o equilíbrio diferencial de membranas e determinar as suas equações de equilíbrio. Além disso, é apresentada uma metodologia puramente geométrica para a determinação de superfícies minimais, baseada na minimização do funcional da área, dado pela soma das áreas das facetas triangulares nas quais a superfície é discretizada. O trabalho discute a formulação matemática do problema e apresenta resultados obtidos tanto por meio das rotinas implementadas no software MATLAB quanto por meio daquelas da biblioteca de otimização deste mesmo software. Finalmente, são realizados alguns exemplos e um teste de convergência, comparando as superfícies resultantes dos métodos numéricos com suas respectivas respostas analíticas. A geometria final de um dos exemplos é verificada por meio da analogia dos filmes de sabão, realizando-se uma análise não-linear de equilíbrio através do software Ansys. As soluções foram bastante satisfatórias, resultando em formas muito próximas das analíticas e com pequenos erros relativos das áreas. O teste de convergência também comprovou que o refinamento da discretização leva a uma solução mais próxima da desejada. Portanto, os procedimentos apresentados podem ser empregados no processo de busca da forma de membranas estruturais. / This dissertation presents a historical review on the theoretical developments on minimal surfaces, highlighting the important analogy between the problems of finding minimal area surfaces and finding membrane surfaces with homogeneous and isotropic stress fields. Some geometric concepts of the parametric surfaces are placed, on the basis of differential geometry, in order to do the differential equilibrium of membranes and to achieve its equilibrium equations. Moreover, a purely geometric methodology for the determination of minimal surfaces is presented, based on the minimization of the area functional, which is computed by the simple addition of a finite number of triangular facet areas in which the surface is divided. It discusses the mathematical formulation of the problem as well as some results obtained with the algorithms implemented in MATLAB and others obtained with the aid of MATLAB optimization routines. Finally, some examples and a convergence test are produced, comparing their analytical and numerical results. The final geometry of one of examples is verified by means of the soap film analogy, with a nonlinear equilibrium analysis through Ansys. The solutions have been sufficiently satisfactory, resulting forms very close to the analytical ones and with small areas relative errors. Convergence test also confirm that the method lead to numerical solutions as close to the analytical one as required, as long as the triangular facets mesh is refined. Therefore, the presented procedures can be used in structural membranes form finding.

A utilização de perfis formados a frio na construção civil tem sido motivada pela elevada eficiência estrutural, expressa pela relação entre capacidade resistente e peso, e pela facilidade de fabricação, caracterizada pela possibilidade de produção de elementos com diferentes seções transversais. Devido à alta esbeltez dos elementos que os constituem, o projeto desses perfis, à temperatura ambiente ou em situação de incêndio, é governado pelos fenômenos de instabilidade local, distorcional e global. O objetivo dessa Tese é o desenvolvimento de ferramentas computacionais que possibilitem a avaliação do comportamento estrutural de pilares de aço formados a frio em situação de incêndio. Para esse fim são desenvolvidos dois programas de computador. O primeiro programa, denominado de ATERM, tem o objetivo de determinar o campo de temperaturas em regime transiente, de estruturas bidimensionais formadas por qualquer material, submetidas a qualquer curva temperatura-tempo de incêndio, com base no método dos elementos finitos. A esse programa acopla-se o programa ATERM-DIM para o dimensionamento em regime plástico de vigas de aço continuamente travadas, em situação de incêndio. Resultados obtidos do ATERM são comparados aos obtidos do programa sueco de computador Super Tempcalc. O segundo programa, denominado de INSTAB, realiza análises lineares e não lineares de estabilidade de perfis formados a frio, considerando as instabilidades local, global e distorcional, empregando o método das faixas finitas splines, para material elastofrágil, considerando a redução das propriedades mecânicas devido à temperatura. Valores de esforços resistentes obtidos pelo programa INSTAB são comparados aos determinados por meio do programa comercial de elementos finitos ANSYS, o qual considera o comportamento elastoplástico do material, e a resultados obtidos por meio de um método simplificado de dimensionamento de perfis formados a frio em situação de incêndio, proposto pelo autor para fins de normatização brasileira. Essa comparação permite avaliar o efeito da elastoplasticidade nos perfis axialmente comprimidos. Pretende-se, com base neste estudo, fornecer subsídios para o desenvolvimento de procedimentos para a verificação estrutural de perfis formados a frio em situação de incêndio. / The use of cold-formed steel profile in construction has increased because of its high efficiency, expressed as the ratio between load capacity and weight, and ease of manufacturing, characterized by the possibility of production of elements with different cross sections. Due the high slenderness ratio of the sections elements, the design of these profiles, either at room temperature or in case of fire, is determined by local, distortional and global buckling phenomena. The aim of this Thesis is to develop computational tools that allow the assessment of the structural behavior of cold-formed steel columns in fire. For this purpose, two softwares are developed. The first one, called ATERM, allows determining the temperature field under transient analysis of two-dimensional structures formed by any material, subjected to any time-temperature fire curve, and is based on the finite element method. This software interacts with another program, ATERM-DIM, used for plastic design of lateral restrained beams steel in fire. Results of ATERM are compared to those obtained from the Swedish software Super Tempcalc. The second program, named INSTAB, can perform linear and nonlinear stability studies of cold formed profiles, taking into account the local, distortional and global buckling, by means of the splines finite strip method for elastofragile material, considering the reduction of mechanical properties caused by the increase in temperature. Resistant values obtained by the INSTAB software are compared with results from the commercial finite element program ANSYS, which considers the plastic behavior of the material and also with results obtained by means of a simplified method for design of cold-formed profiles in fire, proposed by the author for the Brazilian fire standard. This comparison allows analyzing the effect of elastoplasticity in columns of cold-formed steel. Another objective of this study is to provide background for the development of procedures for the structural analysis of cold formed profiles in fire.

O Tanque de Provas Numérico (TPN) é um laboratório pioneiro em hidrodinâmica aplicada e fruto de uma colaboração entre a indústria brasileira de petróleo (PETROBRAS S.A.) e as principais instituições de pesquisa do país. Seu principal objetivo é atuar como parceiro da indústria offshore e de petróleo, colaborando para a obtenção da autossuficiência da produção nacional de petróleo como uma poderosa ferramenta para projeto e análise de sistemas flutuantes de produção de óleo e gás. O coração do TPN é um cluster de computadores SMP que é hoje um dos maiores agrupamentos do Brasil para fins de pesquisa. Um dos focos de atenção do TPN é a aplicação do Moving Particles Semi-implicit Method (MPS) na exploração de soluções para muitos problemas de Engenharia. Por trabalhar sem a necessidade do uso de malhas (método tradicional Euleriano), o método tem diversas aplicações na simulação de corpos flutuantes e na hidrodinâmica aplicada, sendo atualmente utilizado para realização de estudos sobre a influência do movimento de ondas em navios; simulações de fenômenos que envolvem fragmentações; superfícies livres; grandes deformações; dinâmica de fluidos em condições extremas, como é o caso em processos prospecção do petróleo onde muitas vezes é difícil e economicamente inviável fazer ensaios físicos. Devido ao grande número de partículas utilizadas na simulação de sistemas complexos pelo método MPS, é necessário aproveitar de forma eficiente os recursos computacionais disponíveis para a análise de modelos com o refinamento adequado às aplicações práticas. Com tera-FLOPS disponíveis na rede cluster do TPN para modelagem computacional, há uma grande necessidade de uma solução computacional paralela altamente escalonável que, além disto, seja fácil de manutenção e extensibilidade. Dentro desta linha de pesquisa, foi desenvolvida uma solução com essas características através do emprego de modernas técnicas de engenharia de software. / The Numerical Offshore Tank (TPN) is a pioneer laboratory in applied hydrodynamics and result of collaboration between the Brazilian oil (Petrobras SA) and the major research institutions in the country. Its main purpose is to act as a partner of industry and offshore oil, contributing to the achievement of self-sufficiency of domestic oil production as a powerful tool for design and analysis of floating production systems for oil and gas. The heart of TPN is a cluster of SMP computers that is now one of the largest groupings of Brazil for research purposes. One focus of attention of TPN is the application of Moving Particles Semi-implicit Method (MPS) in exploring solutions to many engineering problems. By working without the use of mesh (Eulerian traditional method), the method has several applications in the simulation of floating bodies and applied hydrodynamics, currently being used for studies on the influence of the movement of ships in waves; simulations of phenomena involving fragmentation; free surfaces, large deformations; fluid dynamics in extreme conditions, as is the case in processes where petroleum exploration is often difficult and uneconomical to do physical tests. Due to the high number of particles used in the simulation of complex systems by the MPS method, it is necessary to efficiently take advantage of the computational resources available for the analysis of models with the refinement suitable for practical applications. With tera-FLOPS available in the TPN network cluster for computational modeling, there is a great need for a parallel highly scalable solution which, moreover, must be easy maintenance and extensibility. Within this line of research, we developed a solution with these characteristics through the use of modern software engineering techniques.