Recentemente, surgiu na literatura uma solução analítica das equações de cinética pontual que considera a reatividade como função do tempo, utilizando o método da decomposição. O presente trabalho dá um passo a frente, considerando as equações de cinética pontual em conjunto com efeitos de realimentação de temperatura. Mas, primeiro, uma breve introdução do modelo de cinética pontual e dos aspectos relevantes desta abordagem são apresentados. O trabalho prossegue acrescentando a realimentação de temperatura como uma perturbação de primeira ordem na equação reatividade, e através de manipulações algébricas, o conjunto de equações de cinética pontual passa a ser expresso como uma equação diferencial não linear de segunda ordem. Esta equação é, então, resolvida pelo método de decomposição, ou seja, expandindo as variáveis dependentes como séries infinitas, construindo-se então um sistema recursivo que permite calcular cada um dos termo destas séries. A não linearidade é tratada utilizando os polinômios Adomian. Os resultados aqui obtidos são comparados com a literatura, apresentando variações percentuais máximas da ordem de [0,1%]. Faz-se uma breve análise da convergência e da estabilidade da solução, usando um método baseado no critério de Lyapunov. / An analytical solution of the point kinetics equations to calculate reactivity as a function of time by the decomposition method has recently appeared in the literature. The present work goes one step forward, by considering the neutron point kinetics equations together with temperature feedback effects. But first, a brief introduction of the point kinetics model and the relevant aspects of this approach are presented. The work proceeds adding the temperature feedback as a first order perturbation in the reactivity equation, and using algebraic manipulation, the set of point kinetics equations becomes expressed as a non linear second order differential equation. This equation is then solved by the decomposition method, that is, expanding the dependent variables as infinite series, building a recursive system that allows the evaluation of each term of these series. The non linearity is treated using the Adomian polynomials. The results obtained are compared with literature, with maximum percentage changes of about [0,1%]. A brief analysis of the convergence and stability of the solution is made, using a method based on the Lyapunov criterion.

No presente trabalho, é apresentado um modelo numérico para análise tridimensional de túneis escavados em maciços rochosos reforçados por tirantes aderidos ao longo de todo o seu comprimento. O método numérico utilizado é o método dos elementos finitos. A abordagem usada na modelagem da interação maciço/tirante é o método de homogeneização de meios periódicos. Essa abordagem substitui o meio compósito constituído de rocha e tirantes por um meio homogeneizado equivalente, dispensando a tarefa de discretizar cada tirante separadamente. O comportamento do maciço é considerado elastoplástico perfeito, caracterizado pelo critério de Drucker-Prager, com regra de fluxo associada. O revestimento é considerado elástico, enquanto os tirantes são tratados como inclusões lineares em tração/compressão com lei de comportamento elastoplástica. Com esse modelo foi implementado um código, tridimensional, que simula o processo de escavação pelo método da ativação/desativação de elementos. Esse código é capaz de levar em conta diferentes situações de profundidade, geometria e estado inicial de tensões. Com a finalidade de avaliar o desempenho da abordagem de homogeneização, foi implementado um segundo código, que modela a interação maciço/tirante por meio do método incorporado. Esse método, muito usado em análises de concreto armado, tem em comum com o método de homogeneização a capacidade de lidar com diferentes configurações de reforços utilizando a mesma malha base. Ao final, os dois códigos são utilizados em modelagens de um caso real, o túnel experimental Kielder (Inglaterra), para o qual se dispõe de dados medidos in situ. A comparação destes dados com os resultados das duas modelagens mostra uma aproximação muito boa. / A numerical model for three-dimensional analysis of tunnels excavated in reinforced rock mass by fully grouted bolts is presented in this work. The numeric method utilized is the finite element method and the approach used in the modeling of the interaction mass/bolt is the homogenization of periodic media method. This approach replaces the composite medium formed by rock and bolts by a homogenized equivalent medium eliminating the task of discretizing each bolt separately. The rockmass behavior is considered elastoplastic perfect caracterized by the Drucker-Prager criterion with associated flow rule. The lining is elastic, while the bolts are treated as linear inclusions in tension / compression with elastoplastic behavior law. A three-dimensional code was implemented with this model that simulates the excavation process by means of the activation-deactivation method. This code is able to consider different situations of depth, geometry and initial state of stress. In order to evaluate the performance of the homogenization approach, a second code was implemented to model the interaction rockmass/bolts by means of the embedded method. This method, widely used in analysis of reinforced concrete, has in common with the approach of the homogenization the capacity to deal with different configurations of reinforcements with the same base mesh. Finally, the two codes are used to model a real case, the Kielder experimental tunnel (England), for which in situ measured data are available. The comparison of these data with the results of the two modeling shows a very good approximation.

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Elétrica, 2012. / Submitted by Alaíde Gonçalves dos Santos (alaide@unb.br) on 2013-01-25T12:20:01Z No. of bitstreams: 1 2012_LeandroCarisioFernandes.pdf: 11244387 bytes, checksum: 56968dbdb42928aca7dfd8945faf01cd (MD5) / Approved for entry into archive by Guimaraes Jacqueline(jacqueline.guimaraes@bce.unb.br) on 2013-02-01T12:59:34Z (GMT) No. of bitstreams: 1 2012_LeandroCarisioFernandes.pdf: 11244387 bytes, checksum: 56968dbdb42928aca7dfd8945faf01cd (MD5) / Made available in DSpace on 2013-02-01T12:59:34Z (GMT). No. of bitstreams: 1 2012_LeandroCarisioFernandes.pdf: 11244387 bytes, checksum: 56968dbdb42928aca7dfd8945faf01cd (MD5) / A demanda por serviços de comunicações móveis cresce ano após ano. A implantação de sistema que ofereça esse tipo de serviço requer mecanismos de predição de cobertura de sinal, a fim de otimizar parâmetros usados na fase de planejamento. E importante considerar que, em um enlace de comunicação via rádio, a propagação do sinal não ocorre em situações ideais, na maioria das vezes. A quantidade de fenômenos que influenciam a propagação do sinal é proporcional a complexidade do ambiente. Em regra, em áreas urbanas existe um número maior de obstáculos do que em localidades rurais, o que torna mais complicado estimar a perda de propagação em cidades e metrópoles. Melhorar a estimativa da perda de propagação pode assegurar mais qualidade ao serviço oferecido por operadoras de telecomunicações. Este trabalho mostra como o ambiente urbano pode ser considerado no cálculo da perda de propagação em duas faixas de freqüência. Para a faixa de 900 MHz, são propostos métodos que analisam situações em que a altura da antena transmissora é menor do que a altura média das construções. Os resultados demonstram uma melhora signicativa quando comparados com aqueles obtidos por meio do modelo COST-Walfisch-Ikegami. Em ondas médias, na faixa de 1MHz, propõe-se uma correção da Recomendação ITU-R P.368 de forma a se incluir informações sobre a densidade de área ocupada por construções na localização do receptor. Os valores estimados do campo elétrico estão mais próximos dos dados de medidas do que os obtidos usando-se apenas a recomendação do ITU. O estudo da propagação em ambientes indoor também e um problema de grande interesse prático. Nesse caso, devido a complexidade do ambiente, o uso de equações analíticas torna-se impraticável. Uma forma de contornar essa restrição e usando métodos numéricos. Nesse sentido, será mostrado também como problemas de propagação em ambientes indoor podem ser analisados via métodos numéricos. _______________________________________________________________________________________ ABSTRACT / The demand for mobile communications services grows year after year. Path loss prediction models are necessary to deploy any system that o ers this service. In a wireless communication link, the propagation of the signal occurs mostly over a non line of sight path. The more complex the environment, the greater the number of phenomena that a ect the propagation of the signal. In urban areas, the number of obstacles is greater than in rural areas, and therefore it is more di cult to estimate the path loss in those sites. However, improve estimates of path loss in urban environment can also improve the quality of service o ered by telecommunications operators. Thus, this thesis shows how the urban environment can be considered to calculate the path loss in two frequency bands. Two path loss models are proposed at 900 MHz, when transmitting antenna is lower than the average height of the buildings. The results show a signi cant improvement when compared with COST-Wal sch-Ikegami model. For medium waves, at 1 MHz, a correction of Recommendation ITU-R P.368 is proposed to include information of area occupied by buildings in the receiver's location. The estimated values of the electric fi eld are better than those obtained using the Recommendation ITU-R P.368. Propagation in indoor environments is also a problem with practical interest. Due to the complexity of the environment, the use of analytical equation is impractical. Numerical methods solve this limitation. Thus, it will be shown how to analyse indoor environments using numerical methods.

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.

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.

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.

No desenvolvimento de projetos de iluminação de ambientes, um dos objetivos que se destaca é a busca pelo conforto visual, que emprega metodologias de resolução conhecidas, como o Método dos Lumens e o Método Ponto a Ponto. A luz visível está contida no espectro da radiação térmica e, portanto, o fluxo luminoso pode ser relacionado ao fluxo de radiação térmica. Determinar as posições e as competências das fontes de luz necessárias na superfície de projeto ganha importância quando o comportamento, em termos de uniformidade ou de fluxo radiante, é especificado. O presente trabalho visa a estabelecer diferentes valores de fluxo em duas regiões distintas da superfície do projeto. Por meio do posicionamento das fontes de luz, é estabelecido um fluxo maior na região denominada principal e um fluxo menor na região denominada secundária. A modelagem matemática da radiação térmica (Método das Radiosidades) é aplicada ao projeto de iluminação, considerando-se as características da visão humana e o comportamento das fontes de luz. Na modelagem, é considerada uma cavidade retangular tridimensional com superfícies cinza e com condição de parede fria, na qual o poder emissivo das paredes é nulo. As fontes de luz são representadas por unidades de malha no teto. A relação de equações é resolvida por metodologia inversa, usando o algoritmo de Otimização Extrema Generalizada (GEO). Este algoritmo é classificado como um método de otimização estocástica de busca global para a resolução de sistemas considerados inicialmente mal condicionados. A posição e a potência das fontes luminosas são determinadas pela resolução do sistema de equações, de forma a proporcionar um fluxo de radiação duas vezes maior na região principal em relação à região secundária. A função objetivo do processo consiste em minimizar a diferença entre o fluxo desejado e os valores de fluxo de radiação incidente nas duas regiões da superfície de projeto. Em virtude das características de simetria do problema, a relação é estabelecida para apenas um quarto da cavidade. Assim, por exemplo, aplicar a metodologia com 9 fontes de luz a um quarto da região resulta em 36 fontes de luz em toda a cavidade. Os resultados mostram que é possível encontrar um arranjo de fontes de luz preestabelecendo-se duas condições de potência. / In the development of environmental illumination projects, one of the main goals to be achieved is the visual comfort, which is usually done by known methodologies, like the Lumens Method and the Point by Point Method. Since the visible light is contained in the spectrum of thermal radiation, the luminous flux can be related to the thermal radiation flux. The determination of the position and power of the light sources required by the design surface gains an higher importance whenever a behavior is specified, should it be in terms of uniformity or in therms of radiant flux. In this work, we describe a method that allows the establishment of different flux values in two distinct regions of the design surface, which are referred by the names main region and secondary region. Through the spatial arrangement of the light sources, the method sets a more intense flux in the main region and a less intense one in the secondary region. The mathematical model of thermal radiation, known as Radiosity Method, is applied to the illumination design, along with the characteristics of the human vision and the behavior of light sources. In this model, a rectangular three-dimensional cavity is considered. It has gray surfaces and exhibits the conditions of a cold wall, in which the emissivity power of the walls is null. The light sources are represented by a mesh unit in the ceiling. The system of equations is solved by inversemethodology, using the Generalized Extremal Optimization (GEO) algorithm. This algoritm is classified as being a stochastic optimization method of global search to solve systems that are initially considered ill-conditioned. By solving this system, the position and power of light sources can be determined, and this is done in such a way that the flux radiation in the main region is twice more intense then the one in the secondary region. The target function of the whole process is to minimize the difference between the desired flux and the incident flux radiation values for each one of the two design surface regions. We further explore the problem symmetry, solving the equation system for only a quarter of the cavity. This way, if the methodology is applied with nine light sources into a quarter of the region, the entire cavity will behave as if it has 36 light sources. Our results show that, given two prescribed conditions of power, it is possible to find an arrangement of light sources.