Neste trabalho apresenta-se uma análise numérica tridimensional dos campos de temperatura e velocidade em reservatórios térmicos cilíndricos. Estuda-se o fenômeno da convecção natural laminar e a estratificação vertical da temperatura, que tem grande importância no desempenho de diversos sistemas térmicos, tais como os sistemas de aquecimento de água por energia solar. O código computacional escrito em linguagem FORTRAN resolve as equações da quantidade de movimento e da energia nas três dimensões pelo Método dos Volumes Finitos em malha estruturada. As simulações tridimensionais transientes do processo de resfriamento mostram a formação e evolução do perfil estratificado de temperatura ao longo do tempo. Observa-se a presença de recirculações na região superior do reservatório e a formação das camadas limites térmica e hidrodinâmica próximo das paredes. Através de estudo experimental realizado em paralelo, determinou-se os perfis de temperatura para diversos instantes do processo de resfriamento de um reservatório térmico cilíndrico. Obteve-se também perfis experimentais de temperatura para o caso de consumo de água quente e ingresso de água fria. Em todos os ensaios, além dos perfis de temperatura, mediu-se também as variações da temperatura ambiente. Utilizou-se os dados experimentais colhidos como condição inicial e de contorno para as simulações, possibilitando assim a confrontação numérico versus experimental para a validação da solução numérica. A partir das diversas simulações realizadas, obteve-se correlações para o coeficiente de transferência de calor e para o grau de estratificação térmica no interior do reservatório. Realizou-se ainda a simulação de casos considerando o reservatório térmico interligado a coletores solares e investigou-se tanto a influência da posição do jato de entrada quanto o uso de um anteparo defronte a este jato na preservação da estratificação térmica. Os resultados mostram que, para certas condições do jato de entrada, tanto a escolha adequada da posição do jato de entrada quanto o uso de um anteparo defronte a este jato permitem construir perfis de temperatura com melhor estratificação térmica. / This work presents a three-dimensional numerical analysis of the temperature and velocity fields in cylindrical thermal storage tanks. It is studied the phenomenon of the laminar natural convection and the vertical stratification of the temperature, that has great importance in the performance of diverse thermal systems, such as solar energy water heating systems. The numerical solution is obtained through the Finite Volume Method in structured mesh. The transient three-dimensional simulations of the cooling process show the formation and evolution of the profile of temperature along the time. It is observed the presence of recirculations in the top of the tank and the formation of walls thermal and hydrodynamic boundary layers. Through an experimental study, it was determined the profiles of temperature for several instants of time of the cooling process of a cylindrical thermal storage tank. It was also obtained experimental profiles of temperature for the case of outlet hot water and inlet cold water, as well as measured the ambient air temperature. In some cases the numerical simulation applied experimental data as initial and boundary condition, in order to compare numerical and experimental results, for the validation of the developed computational code. From the simulations, it was obtained correlations for the heat transfer coefficient and degree of thermal stratification in the interior of the reservoir. Others simulations were made with the thermal thank connected to solar collectors, investigating the influence of the position of the inlet water jet, and the use of an inner flat plate in front of this jet, in the preservation of the thermal stratification. The results show that, for certain conditions of the jet, the adequate choice of its position, as well as the use of a flat plate, allows to obtain more stratified temperature profiles.

A presente dissertação contempla o desenvolvimento e validação de um aplicativo para efetuar simulações de hadronterapia. Este aplicativo foi empregado no estudo das reações nucleares de fragmentação que ocorrem para as energias e partículas envolvidas no tratamento. Estas reações representam a chave para o planejamento, tratamento, e pós-análise de precisão. Inicialmente é realizada uma revisão a respeito do tratamento de tumores profundos utilizando prótons e íons. Os aspectos físicos, biológicos e técnicos da utilização de feixes energéticos destas partículas são discutidos. A situação atual do tratamento e perspectivas são mencionadas. Faz-se também referências às di culdades do planejamento devido a ausência de ferramentas precisas para uma estimativa el da deposição de energia. A sequência do texto contém uma breve revisão dos códigos de Monte Carlo que podem ser adotados para a simulação do tratamento. Uma ferramenta que empregue o método de Monte Carlo é uma viável solução para levar em conta todos os efeitos físicos provocados pelas interações das partículas do feixe com os tecidos biológicos. Neste mesmo contexto, são discutidos os processos físicos relevantes que precisam ser considerados nas simulações de hadronterapia. O código de Monte Carlo Geant4 é, então, apresentado. São descritos os recursos que são utilizados para as simulações do tratamento. Este toolkit foi empregado para o desenvolvimento de um aplicativo, sendo este o primeiro item de originalidade do presente trabalho. O aplicativo desenvolvido é discutido e suas capacidades e funcionalidades são relatadas. Posteriormente, são apresentadas simulações utilizando o aplicativo supracitado. Os resultados são confrontados com dados experimentais e as limitações dos modelos físicos disponíveis são analisadas. O aplicativo demonstra, pela avaliação dos resultados, que: consiste em uma possível ferramenta para as estimativas precisas necessárias para a hadronterapia; é versátil, podendo ser empregado para um amplo intervalo de energias e para diferentes partículas do feixe; as funcionalidades podem ser facilmente estendidas para se moldar à proposta da aplicação devido ao paradigma da orientação a objetos; o problema de tempo computacional para o emprego do método de Monte Carlo no planejamento de tratamento pôde ser contornado com a implementação de processamento paralelo. Além disso, as análises demonstram que os modelos físicos necessitam modi cações na faixa de energia da hadronterapia para realizar uma melhor descrição do processo de fragmentação nuclear. Os resultados originais da dissertação abrem caminho para tornar as simulações em hadronterapia um instrumento de precisão no combate contra tumores situados em localizações que não permitem procedimentos invasivos. / The present master thesis comprises the development and validation of an application to perform a simulation of hadrontherapy. This application was employed in the study of nuclear fragmentation reactions that occur in the energy range and with particles involved in cancer treatment. These reactions are the key to acquire high precision in the planning, treatment and post-analysis of the treatment. A review is initially performed with regard to the treatment of deep seated tumors with beams of protons and ions. The physical, biological and technical aspects of the utilization of energetic beams using these particles are discussed. The current situation of the treatment and perspectives are indicated. In addition, reference is made to the di culties in planning due to the lack of precise tools for an accurate estimative of the energy deposition. The sequence of the text contains a brief review of the Monte Carlo codes that may be used for simulation of the treatment. A tool which would employ the Monte Carlo method shall be a reliable way to determine all the physical e ects caused by the interaction of the particles from the beam with the biological tissues. In this context, relevant physical processes which may be considered in the hadrontherapy simulation are discussed. The Monte Carlo code Geant4 is presented. The resources used for simulation of treatment are described. The toolkit is applied for the development of an application and that is novelty of the present work. The developed application is described and its capabilities and functionalities are listed. At a later moment simulations are presented by utilizing the application. The results are compared to experimental data and the limitations of the available physics models are analyzed. Evaluation of the results indicates that the application represents a useful tool for precise estimations that are needed for the hadrontherapy; further it is versatile and it can be used for a wide range of energies and di erent particles from the beam; moreover the functionalities can be easily extended in order to mold the terms of the application due to the object-oriented paradigm; last the issue regarding computation time for the use of the Monte Carlo method at treatment planning is solved with the implementation of parallel processing. Besides, the analyses indicate that the physics models need to be modi ed in the energy range suited for hadrontherapy, in order to achieve a better description of the process of nuclear fragmentation. These novel results shown in the present thesis open a pathway to make simulation in hadrontherapy a precise instrument to ght tumors located in deep seated areas where an invasive procedure is not be possible.

Este trabalho tem por objetivo analisar a aplicação do Método Reaction Diffusion Manifold (REDIM) na simulação de chamas axisimétricas, difusivas e laminares. Modelos detalhados de cinética química são usualmente computacionalmente proibitivos para a simulação de fluxos reativos complexos e, portanto, modelos reduzidos são necessários. Ferramentas automáticas de redução de modelos normalmente exploram a estrutura multi-escalar natural dos sistemas de combustão. A nova técnica REDIM aplica o conceito de variedades invariantes para tratar, também, a influência dos processos de transporte no modelo reduzido. Assim, ela supera um problema fundamental na redução de modelos em negligenciar o acoplamento do transporte molecular com o processo termodinâmico. Para o estudo, considera-se como caso teste uma chama de metano/ar (65%CH4 + 35%N2), difusiva, laminar, axisimétrica, a pressão atmosférica para validar a metodologia aqui apresentada. Primeiramente, REDIM uni e bidimensionais são computadas e tabuladas em tabelas lookup. Então, o sistema completo de equações governantes é projetado na REDIM e implementado no código Gascoigne através de uma nova biblioteca addon para tratar as tabelas REDIM. O sistema projetado de equações governantes é discretizado pelo Método de Elementos Finitos (MEF) e resolvido por uma iteração GMRES precondicionada por um método multigrid geométrico. Refinamento de malha local, malha adaptativa e paralelismo foram aplicados para garantir eficiência e precisão. Os resultados numéricos obtidos usando a técnica REDIM mostram boa concordância com os dados de simulação numérica detalhada e dados experimentais reportados na literatura. / The goal of this work is to analyze the use of automatically reduced chemistry by the Reaction-Diffusion-Manifold (REDIM) method in simulating axisymmetric laminar coflow diffusion flames. Detailed chemical kinetic models are usually computationally prohibitive for simulating complex reacting flows and, therefore, reduced models are required. Automatic reduction models approaches usually exploit the natural multiscale structure of combustion systems. The novel REDIM approach applies the concept of invariant manifolds to treat also the influence of the transport processes on the reduced model, which overcomes a fundamental problem of model reduction in neglecting the coupling of molecular transport with thermochemical processes. We have considered a previously well studied atmospheric pressure nitrogen-diluted methane/air flame as a test case to validate the methodology presented here. First, one-dimensional and two-dimensional REDIMs were computed and tabulated in lookup tables. Then, the full set of governing equations are projected on the REDIM and implemented in the object-oriented C++ Gascoigne code with a new add-on library to deal with the REDIM-tables. The projected set of governing equations has been discretized by the Finite Element Method (FEM) and solved by a GMRES iteration preconditioned by a geometric multigrid method. Local grid refinement, adaptive mesh and parallelization are applied to ensure efficiency and precision. The numerical results obtained using the REDIM approach have shown very good agreement with detailed numerical simulations and experimental data.

Nesse trabalho temos como objetivo simular diversos estágios da ameba Dictyostelium discoideum durante o seu ciclo de vida social. Tal organismo apresenta um ciclo de vida típico no qual se comporta como uma ameba Eucarionte, alimentando-se de bactérias e dividindo-se por mitose. É quando as fontes de nutrientes no ambiente se tornam escarças que tem início o ciclo de vida social no qual as amebas começam a emitir um sinal químico que guiará as mesmas com objetivo de formar agregados com aproximadamente 50 mil amebas. Dentro de tal agregado as amebas se diferenciam e formam um ser pluricelular na forma de uma lesma que tem como objetivo transportar as amebas para uma região com maior probabilidade de sobrevivência. Para simular tal organismo iremos mesclar dois modelos diferentes para reproduzir a sinalização química entre as amebas e o seu movimento. O modelo utilizado para descrever a sinalização química do composto Adenosina Monofosfato cíclico 3'5' foi desenvolvido de acordo com a análise dos processos bioquímicos envolvidos e os parâmetros foram obtidos experimentalmente. A movimentação das amebas é implementada através de um modelo de animóides com velocidade fixa onde o ângulo de movimento varia de acordo com regras simples que já provaram ser eficientes na simulação desde grupos de animais a processos de segregação celular. Nessa tese mostraremos como o método utilizado para resolver o modelo de sinalização química não é somente mais eficiente numericamente, como também é capaz de apresentar diversas das propriedades da ameba vistas nos experimentos até o momento da agregação, entre elas: sincronização; adaptação; formação de espirais; padrões de reação difusão. Ao implementarmos o movimento das amebas verificamos outras propriedades como: caminhos preferenciais de deslocamento até o centro do agregado (streaming); formação de mound e movimento coletivo circular. Nesse trabalho além de conseguirmos unir diversas propriedades em um único modelo, também conseguimos demonstrar que a origem do movimento helicoidal coletivo das amebas dentro da lesma pode ser proveniente de simples propriedades físicas ao contrário de algumas vertentes que afirmam ser proveniente da quimiotaxia em um sistema com ondas de sinal químico helicoidal. / In this work we aim to simulate different stages of the amoeba Dictyostelium discoideum through its social cycle. This organism has a typical life cycle in which it behaves like an Eukaryotic amoeba, feeding on bacteria and reproducing through mitosis. It is when the sources of nutrients in the environment become scarce that the social cycle begins. In this cycle the amoebae start to emit a chemical signal that will guide them to form aggregates with approximately 50 thousand amoebae. Within such aggregate the amoebae differentiate themselves to create a multicellular being in the form of a slug, which aims to transport the amoebae to a region with a higher probability of survival. In order to simulate such organism we have merged two different models to reproduce the chemical signaling among amoebae and their movement. The model chosen to describe the chemical signaling of the compound cyclic Adenosine Monophosphate 3'5' has ben developed according to the biochemical processes and the parameters obtained experimentally. The amoebae movement is implemented through a boids model with fixed speed where the angle of movement varies according to simple rules that have proven effective in other simulations of animal groups and cellular segregation. In this thesis we show how the method used to solve the model chemical signaling is not only numerically more efficient, but also capable of reproducing several properties of the amoeba seen in experiments until the aggregation stage, including: synchroniiation, adaptation, spiral formation and reaction-diffusion patterns. Upon implementation of amoeba movement we verified other properties such as streaming, mound formation and collective rotational movement. In this thesis we in addition to uniting several properties in a single model, were also able to demonstrate that the origin of the helicoidal movement within the slug can come from simple physical properties, in contrast to some studies which claim the origin to be from chemotaxis in a system with scroll waves of chemical signaling.

O presente trabalho abordou a simulação numérica de problemas relacionados a um tópico de grande importância na área médica: o comportamento de tecidos ósseos. Foram utilizadas ferramentas computacionais, fundamentos de Mecânica dos Sólidos e teorias baseadas em estudos e experimentos, para simular o comportamento desses tecidos vivos. Aqui se desenvolveu uma metodologia para criar modelos computacionais de ossos fazendo uso de leituras tomográficas, designando as propriedades em função do valor das densidades. A metodologia foi aplicada juntamente com algoritmos genéticos e o método dos elementos finitos para obter a melhor posição de implantes ortodônticos, visando uma homogeneização da tensão no osso na região de colocação do implante. Também foi possível fazer uma comparação entre modelos utilizados na literatura e o desenvolvido aqui. Em particular, foi demonstrado que o uso de materiais isotrópicos homogêneos e materiais com propriedades variáveis de acordo com a densidade levam a resultados diferentes para a otimização da orientação do implante. Foi criado um modelo computacional de remodelamento ósseo com isotropia transversal, com base no modelo de crescimento ósseo isotrópico de Stanford e utilizando as relações constitutivas dos tecidos ósseos. Foi realizada uma comparação entre os modelos isotrópico e isotrópico transversal na simulação do processo de remodelamento ósseo, sendo possível destacar as principais diferenças entre eles. A seguir foi realizada a aplicação do remodelamento ósseo para a otimização de próteses femorais. Foram construídos modelos computacionais integrando o método dos elementos finitos, um modelo matemático do comportamento adaptativo do material ósseo e o método dos algoritmos genéticos, visando a minimização da perda de massa óssea após a inserção de próteses em um osso hospedeiro. Foram realizadas três otimizações para combinações diferentes de parâmetros de projeto de uma prótese femoral. Na primeira, foi realizada a otimização da geometria e do módulo de Young da prótese. Na segunda, foi otimizado apenas o módulo de Young da prótese. Na terceira, foi proposto um modelo de prótese bimaterial, no qual os módulos de Young e a posição da divisão entre os materiais são os parâmetros de projeto escolhidos. A partir dos resultados obtidos, foi possível fazer comentários e sugestões sobre o material utilizado para construção de próteses femorais, bem como sobre a aplicabilidade da metodologia empregada na simulação de problemas de biomecânica. / The present work is addressed to the numerical simulation of problems related to an important topic in the medical area, namely, bone tissues. Computational methods, solid mechanics and theoretical/experimental studies were used to simulate the behavior of live bone tissues. A methodology was developed to create computational models using tomography data, assigning the mechanical properties from the corresponding densities map. The methodology was applied with genetic algorithms and the finite element method to predict the optimal position of orthodontic implants, aiming the homogenization stress in the bone. It was possible to compare some models found in literature and the one developed here. In particular, it was showed that homogeneous isotropic materials and materials with variable density result in different predictions in the optimization of the implant orientation. A computational scheme for bone remodeling was developed using three material behavior laws: transversally isotropic, isotropic Stanford model and classical constitutive relations for bone tissues. A comparison between the isotropic model and the transversally isotropic in the simulation of the bone remodeling was carried out, and the main differences between them were highlighted. Then the proposed bone remodeling scheme was applied to optimize femoral prosthesis. The computational models integrated the finite element method, a mathematical model to describe the behavior of adaptive bone material, and the genetic algorithm. The optimization aimed reduction of bone absorption after the insertion of a prosthesis in a bone host. Three cases of cemented femoral prosthesis were optimized. In the first case, the Young Modulus and geometry of the prosthesis were selected as design variable. In the second case, only the Young Modulus was used as design variable. In the third case, a bi-material prosthesis was proposed was proposed and optimized. Conclusions about the applicability of the proposed formulation for the simulation of biomechanics problems are drawn.

Durante a realização de projetos de sistemas de refrigeração assistidos por sistemas de armazenagem de energia, na forma latente, é necessário que se saiba de antemão as taxas de transferência de energia dos mesmos. Na literatura soluções para este problema são escassos e não abordam todas as geometrias utilizadas para o armazenamento de energia latente, inexistindo correlaçõesmatemáticas que permitam aos projetistas prever o comportamento da fusão e da solidificação do material de mudança de fase (MMF) e as taxas de transferência de energia térmica. Na solução destes problemas são usualmente adotadas hipóteses simplificativas para amenizar o grau de complexidade da descrição matemática do problema Uma destas hipóteses constitui-se em desconsiderar o processo advectivo na transferência de calor no MMF, em sua fase líquida. Outra hipótese adotada, quando o MMF é a água pura, consiste em utilizar um comportamento linear para a massa específica ou um comportamento médio, consistindo em um afastamento da realidade. Este trabalho descreve esforços no sentido de aproximar estas soluções cada vez mais da realidade fisica do problema. São apresentadas simulações para a mudança de fase da água pura levando em conta a inversão da massa específica, em geometria polar, na forma transiente, utilizando o método de Volumes Finitos para solução das equações governantes. As soluções apresentadas, otimizadas quanto às malhas espacial e temporal, envolvem condições de contorno de primeira e terceira espécies Como resultado das soluções são apresentados o Nusselt local, ao longo da parede do cilindro além do Nusselt médio. Correlações matemáticas para o número de Nusselt médio durante a fusão e para o volume de MMF fundido, são propostos e apresentados. A evolução geométrica da fronteira de mudança de fase também é avaliada para as duas condições de contorno e diversos números de Rayleigh. Comparações entre o volume fundido considerando apenas o processo difusivo e o volume fundido utilizando a convecção são realizadas. Para a verificação da independ~encia de malha aplica-se o índice de convergência de malha (GCI) . Algumas comparações qualitativas com soluções experimentais disponíveis na literatura são apresentadas.

Quantificação de incerteza e risco tem aplicação direta no projeto de limites de cava e na programação de produção. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo na ordem de milhões de dólares. Preço de mercado do bem mineral, o custo de capital, taxas de atratividade, etc, são fatores usualmente testados para analisar a viabilidade de um empreendimento, porém, raramente é considerada a incerteza relacionada ao atributo geológico em questão. O propósito de um planejamento de lavra tem sido prover subsidio para o engenheiro de minas decidir sobre a capacidade de minerar determinadas unidades de lavra do depósito a partir de teores estimados. Salienta-se porém que existe, a partir dos dados amostrais, incertezas a respeito do modelo geológico que devem ser consideradas ao se projetar uma cava e desenvolver a lavra. A simulação geoestatistica tem a capacidade de produzir múltiplos modelos equiprováveis, os quais podem ser avaliados independentementecomo possíveis cenários do depósito mineral. Simulação condicional, ao contrário de técnicas de interpolação e estimativa, provê respostas para questões associadas a risco devido à variações nos teores do modelo geológico. Ao gerar múltiplos cenários tem-se acesso à probabilidade e conseqüentemente às variações de retorno financeiro e rotas de extração de minério de um projeto. o presente trabalho tem como objetivo a investigação de novas metodologias que contribuam para a construção de cenários de planejamento de lavra e avaliação do impacto provocado nestes pela incerteza fornecida a partir de modelos simulados. As respostas buscadas dentro da abordagem sugerida por esse trabalho, compreendem a definição de uma metodologia para controle e planejamento de lavra a médio e longo prazo por incorporação da flutuabilidade local associada ao minério, avaliando a sensibilidade do retorno financeiro e o traçado de lavra em relação ao método de geração do modelo geológico de teores. Para solucionar o problema em questão, sugere-se a geração de modelos estocásticos e a alimentação de múltiplos modelos selecionados por critérios considerados relevantes para a delimitação do espaço de incerteza. A aplicação de funções de transferência tais como mudança de suporte e seqüenciamento de produção sobre os modelos simulados é a chave para a obtenção de respostas a respeito do impacto sobre a lucratividade de um projeto. Ao alimentar à essas funções modelos equiprováveis que preservam as características estatísticas e a conectividade espacial dos dados amostrais mas que invariavelmente possuem diferentes distribuições de teores tem-se a dimensão em termos econômicos do risco associado à incerteza geológica. Foi confrontado o retorno financeiro produzido por modelos simulados e verificou-se para o caso especifico que os métodos de simulação geoestatistica empregados não produziram diferenças significativas quando comparados casos selecionados segundo os mesmos critérios. Porém, entre cenários extremos gerados pelo mesmo algoritmo de simulação foram verificadas desigualdades relevantes. Verificou-se também a transferência de minério entre diferentes classes de qualidade ao proceder-se com variações nas dimensões dos blocos de lavra.

A simulação paralela de eventos é uma área da computação que congrega grande volume de pesquisas, pela importância em facilitar o estudo de novas soluções nas mais diferentes áreas da ciência e tecnologia, sem a necessidade da construção de onerosos protótipos. Diversos protocolos de simulação paralela podem ser encontrados, divididos em dois grandes grupos de acordo com o algoritmo empregado para a execução em ordem dos eventos: os conservadores e os otimistas; contudo, ambos os grupos utilizam trocas de mensagens para a sincronização e comunicação. Neste trabalho, foi desenvolvido um novo protocolo de simulação paralela, fazendo uso de memória compartilhada, o qual foi implementado e testado sobre um ambiente de estações de trabalho, realizando, assim, simulação paralela com uso de memória compartilhada distribuída. O protocolo foi desenvolvido tendo como base de funcionamento os protocolos conservadores; utilizou diversas características dos mesmos, mas introduziu várias mudanças em seu funcionamento. Sua execução assemelha-se às dos protocolos de execução síncrona, utilizando conceitos como o lookahead e janelas de tempo para execução de eventos. A principal mudança que o novo protocolo sofreu foi proporcionada pelo acesso remoto à memória de um LP por outro, produzindo diversas outras nas funções relativas à sincronização dos processos, como o avanço local da simulação e o agendamento de novos eventos oriundos de outro LP. Um ganho adicional obtido foi a fácil resolução do deadlock, um dos grandes problemas dos protocolos conservadores de simulação paralela. A construção de uma interface de comunicação eficiente com uso de memória compartilhada é o principal enfoque do protocolo, sendo, ao final da execução de uma simulação, disponibilizado o tempo de simulação e o tempo de processamento ocioso (quantia utilizada em comunicação e sincronização). Além de uma implementação facilitada, propiciada pelo uso de memória compartilhada ao invés de trocas de mensagens, o protocolo oferece a possibilidade de melhor ocupar o tempo ocioso dos processadores, originado por esperas cada vez que um LP chega a uma barreira de sincronização. Em nenhum momento as modificações efetuadas infringiram o princípio operacional dos protocolos conservadores, que é não possibilitar a ocorrência de erros de causalidade local. O novo protocolo de simulação foi implementado e testado sobre um ambiente multicomputador de memória distribuída, e seus resultados foram comparados com dois outros simuladores, os quais adotaram as mesmas estratégias, com idênticas ferramentas e testados em um mesmo ambiente de execução. Um simulador implementado não utilizou paralelismo, tendo seus resultados sido utilizados como base para medir o speedup e a eficiência do novo protocolo. O outro simulador implementado utilizou um protocolo conservador tradicional, descrito na literatura, realizando as funções de comunicação e sincronização através de trocas de mensagens; serviu para uma comparação direta do desempenho do novo protocolo proposto, cujos resultados foram comparados e analisados.

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

O presente trabalho apresenta uma comparação entre resultados computacionais, obtidos através do "software" EnergyPlus e experimentais do comportamento térmico de um ambiente condicionado e não condicionado. Para tanto, monitoraram-se os dados climáticos de radiação, velocidade do vento e temperatura, no período de 11 a 20 de janeiro de 2002 e produziu-se um arquivo climático. Simultaneamente, fez-se a aquisição das temperaturas de uma sala-teste, localizada no terceiro pavimento de um prédio na cidade de Porto Alegre, bem como das salas adjacentes. As temperaturas do ar de insuflamento e de retorno dos condicionadores de ar, localizados na sala-teste, foram medidas durante o dia, em seis dias do período de monitoramento. Mediu-se também a velocidade do ar de retorno e determinou-se a potência sensível de refrigeração. Os ganhos de calor interno da sala também foram medidos e utilizaramse as formulações apresentadas pela ASHRAE, 2001, para determiná-los em relação às pessoas e à infiltração. Tais dados foram declarados ao programa como variáveis de entrada. As simulações do comportamento térmico da sala-teste foram implementadas informando-se ao programa a temperatura das salas ou os coeficientes de uma equação representativa das mesmas. Por considerar que a primeira representava melhor as condições térmicas das salas contíguas, utilizou-se esta modalidade para análise As simulações foram conduzidas, alterando-se opções fornecidas pelo programa: modelo de céu isotrópico e anisotrópico, coeficiente de convecção simples e detalhado. Os resultados da carga térmica e temperatura ambiente da sala-teste obtidos nas simulações foram comparados com os dados experimentais do período de monitoramento. A melhor concordância foi obtida para o modelo anisotrópico, coeficiente de convecção detalhado. Constatou-se uma grande discrepância entre as cargas térmicas, para o modelo de convecção simples. Assim, conclui-se que o "software" EnergyPlus representa bem o comportamento térmico de uma edificação "termicamente pesada" para coeficiente de convecção detalhado, necessitando pesquisa para as demais edificações.