Este trabalho apresenta uma metodologia para a simulação numérica de eventos de capotamento de veículos automotores. Pois apesar de o rolamento veicular não ser o mais comum tipo de acidente, ele é o de maior importância com respeito a lesões e traumas provocados nos ocupantes. Assim, reduzir a incidência de mortes e lesões graves tem aumentado a importância de simulações computacionais e testes físicos. Inicia-se com o desenvolvimento de um modelo computacional de multicorpos para entendimento e avaliação dos parâmetros da suspensão que influenciam no fenômeno de capotamento. Os resultados deste modelo de multicorpos são comparados com testes experimentais. Na segunda etapa são desenvolvidos modelos numéricos computacionais de veículos automotivos para simulação e comparação de diferentes normas e procedimentos que avaliam a resistência da estrutura do teto em suportar as cargas decorrentes do fenômeno de capotamento. Este estudo também apresenta modelos computacionais de elementos finitos para simular eventos de capotamento e predizer possíveis lesões ocasionadas à cabeça, pescoço, tórax e coluna. Modelos numéricos de um veículo de passageiros e de um veículo utilitário esportivo (SUV) são simulados com manequim antropomórfico para representar o condutor. São verificados os riscos de lesões e traumas ao condutor se o mesmo estiver usando ou não cinto de segurança. O risco de lesões ao condutor de um veículo de passageiros é comparado ao risco de um condutor de um SUV. A metodologia computacional desenvolvida mostrou-se bastante eficiente na avaliação da estrutura dos veículos em eventos de capotamento. Resultados computacionais demonstram que algumas normas e procedimentos, que atualmente são amplamente utilizados na indústria automotiva para avaliar a estrutura do teto dos veículos, podem não prevenir a ocorrência de lesões e traumas aos ocupantes. Além disso, os resultados das simulações sugerem que a carga utilizada para avaliar a resistência do teto dos veículos seja aumentada. Demonstrando que métodos numéricos podem ser utilizados na avaliação da estrutura dos veículos e verificação de riscos de lesões ao condutor, reduzindo a quantidade de testes experimentais a serem realizados. / This work presents a methodology for the numerical simulation of vehicle’s rollover events. Even though the rollover is not the most frequent type of accident, it is of the greatest significance with respect to injury and trauma caused to the vehicle occupants. The need to reduce death incidence and serious injuries has increased the importance of computational simulations and prototype testing. Firstly, a multi-body model is developed for understanding and evaluating the parameters of the suspension that most influence in the rollover phenomenon. The results of this numerical model are compared with experimental tests. In the second part, numerical models of automotive vehicles are developed for simulation and comparison of different standards and procedures that evaluate the roof resistance with respect to withstand the rollover loads. This study also presents finite element models to simulate rollover events and to predict possible injuries caused in the head, neck, thorax and cervical spine. Numerical models of passenger's vehicle and a sport utility vehicle (SUV) are simulated including anthropomorphic dummy to represent the driver. The injury risks and traumas are verified to the driver considering belted and unbelted dummies. The driver's risk of injury when driving a passenger car is also compared with the risks when driving a SUV. The computational methodology developed proved to be efficient for the evaluation of the vehicle's roof structure in rollover events. Regarding the computational results, some standards and procedures which are broadly used by the automotive industry to evaluate the vehicle's roof resistance, cannot prevent the occurrence of injuries and traumas to the occupants. Moreover, according to the results of numerical simulations, it is suggested that the load required evaluating the car's roof resistance should be increased. It is shown that numerical methods can be used to assess the structure of the vehicle and verifying the driver's injury risk, reducing the amount of experimental tests to be performed.

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.

Este trabalho tem por objetivo simular e analisar uma usina termelétrica a carvão em várias condições de funcionamento. A usina simulada neste trabalho é a AVV 1 localizada em Copenhague, Dinamarca. A AVV 1 é uma usina de geração de potência e aquecimento distrital que pode funcionar em várias condições de carga. A simulação da usina supracitada foi tema de um concurso de simuladores proposto no congresso ECOS 2003 realizado em Copenhague, Dinamarca. Para a realização deste trabalho foi construído um programa na linguagem FORTRAN 90. Cada componente da usina é modelado através de equações de balanço de massa e energia, e o sistema completo tem sua solução obtida pelo método de substituição sucessiva. Para viabilizar essa solução é necessário também implementar uma rotina de cálculo de propriedades do fluido de trabalho. No caso estudado, o fluido de trabalho da usina é a água e a formulação utilizada para o cálculo de suas propriedades nos diversos estados é a IAPWS IF-97. A usina é simulada em dois modos de operação: modo de condensação, onde ocorre apenas geração de eletricidade, e em modo de contrapressão, onde há geração de eletricidade e aquecimento distrital, conforme nomenclatura sugerida pela organização do concurso No modo de condensação, são feitas quatro séries de simulações variando a carga de 100% a 40%. Cada série contém um conjunto de hipóteses quanto à variação das eficiências isoentrópicas e pressões das turbinas em função da vazão mássica. No modo de contrapressão, a usina é simulada funcionando com 100% da carga. O programa desenvolvido calcula as propriedades para qualquer ponto de trabalho ao longo da planta, assim como a eficiência da mesma, a potência gerada, e todas as vazões mássicas pertinentes. Além disso, é feita também uma análise exergética da planta. A simulação demonstrou que a planta possui uma eficiência global de 42,02% com uma geração de 250,2 MW em 100% de carga no modo de condensação. Nessas mesmas condições, do ponto de vista exergético, a eficiência encontrada é de 37,21%. No modo de contrapressão, a usina apresenta uma eficiência exergética de 40,19% com um aproveitamento energético de 90,55%. Por fim, é possível também verificar a comportamento da eficiência da planta e a variação de água de resfriamento do condensador com a carga. Os resultados gerados são próximos àqueles encontrados pelos diversos pesquisadores que abordaram o problema.

Este trabalho apresenta uma sistemática para realizar a otimização numérica de pré-formas e de matrizes em problemas de forjamento axissimétricos e em estado plano de deformações. Para este fim, desenvolveu-se um código computacional composto basicamente de três módulos: módulo de pré-processamento, módulo de análise e módulo de otimização. Cada um destes foi elaborado acrescentando rotinas em programas comerciais ou acadêmicos disponíveis no GMAp e no CEMACOM. Um programa gerenciador foi desenvolvido para controlar os módulos citados no processo de otimização. A abordagem proposta apresenta uma nova função objetivo a minimizar, a qual está baseada em uma operação booleana XOR (exclusive or) sobre os dois polígonos planos que representam a geometria desejada para o componente e a obtida na simulação, respectivamente. Esta abordagem visa eliminar possíveis problemas geométricos associados com as funções objetivo comumente utilizadas em pesquisas correlatas. O trabalho emprega análise de sensibilidade numérica, via método das diferenças finitas. As dificuldades associadas a esta técnica são estudadas e dois pontos são identificados como limitadores da abordagem para problemas de conformação mecânica (grandes deformações elastoplásticas com contato friccional): baixa eficiência e contaminação dos gradientes na presença de remalhamentos. Um novo procedimento de diferenças finitas é desenvolvido, o qual elimina as dificuldades citadas, possibilitando a sua aplicação em problemas quaisquer, com características competitivas com as da abordagem analítica Malhas não estruturadas são tratadas mediante suavizações Laplacianas, mantendo as suas topologias. No caso de otimização de pré-formas, o contorno do componente a otimizar é parametrizado por B-Splines cujos pontos de controle são adotados como variáveis de projeto. Por outro lado, no caso de otimização de matrizes, a parametrização é realizada em termos de segmentos de reta e arcos de circunferências. As variáveis de projeto adotadas são, então, as coordenadas das extremidades das retas, os raios e centros dos arcos, etc. A sistemática é fechada pela aplicação dos algoritmos de programação matemática de Krister Svanberg (Método das Assíntotas Móveis Globalmente Convergente) e de Klaus Schittkowski (Programação Quadrática Sequencial – NLPQLP). Resultados numéricos são apresentados mostrando a evolução das implementações adotadas e o ganho de eficiência obtido.

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

A elaboração de um projeto de obra de edificação é fruto da interação entre profissionais de diversas disciplinas desenvolvendo, simultaneamente, suas opções e decisões com relação ao mesmo. É necessário considerar a qualidade do processo de elaboração do projeto dependente de uma efetiva comunicação entre os membros da equipe. A indústria da construção civil está sendo beneficiada pelos recursos oferecidos pela tecnologia da informação como, por exemplo, os sistemas colaborativos (extranets de projeto). Estes sistemas corroboram na integração e comunicação entre os membros de um projeto. O objetivo da presente pesquisa foi a apresentação do desenvolvimento de modelo validado de ferramenta computacional projetada para complementar um sistema colaborativo pré-existente. Esta ferramenta foi denominada Sistema de Programação e Controle do Processo de Projeto (SIPROCON/PP). A finalidade do SIPROCON/PP é auxiliar no processo de tomada de decisão dos coordenadores através do monitoramento do cronograma e do fluxo do intercâmbio de informações das atividades geradoras de informações neste processo. O processo de desenvolvimento da ferramenta foi delineado em harmonia com a organização da pesquisa. Neste sentido, é assinalado o uso da prototipação e da simulação como estratégias de pesquisa, utilizados respectivamente na construção e validação do modelo da ferramenta. Ante a necessidade do experimento, foi criada uma atividade lúdica através da qual foi reproduzido o intercâmbio de informações de um processo de projeto denominada Sistemática de Simulação do Processo de Projeto (SS/PP) para validação do modelo da ferramenta. Logo, por meio do uso das simulações foram obtidos os dados correspondentes a validação da modelagem e aplicação do SIPROCON/PP e, ao mesmo tempo, necessários à conclusão da pesquisa. Diante disso, a partir da análise dos dados coletados é apontada a melhoria na qualidade do processo de elaboração do projeto considerando a formação de uma equipe equilibrada de projeto. Isso, fruto do conhecimento mais acurado do coordenador sobre o desempenho de cada projetista.

Este trabalho apresenta um modelamento computacional dos diversos componentes de um ciclo de Rankine regenerativo em operação transiente, utilizando as leis fundamentais de balanço de massa e energia, além de algumas equações constitutivas. O cálculo das propriedades termodinâmicas foi realizado através de rotinas desenvolvidas no Grupo de Estudos Térmicos e Energéticos (GESTE) da UFRGS. Os componentes com comportamento inercial são o condensador e o tanque de água de alimentação. Os demais respondem instantaneamente. O método de solução das equações algébricas é seqüencial e o avanço no tempo é totalmente implícito. Primeiramente o sistema é colocado em regime permanente, a partir do qual é introduzida uma mudança de carga na turbina em um determinado intervalo de tempo. Algumas variáveis prescritas foram obtidas a partir de dados reais de uma usina termelétrica. A partir desse modelamento foi possível investigar o comportamento dinâmico do condensador e do tanque de alimentação e de que maneira estes efeitos eram percebidos pelos demais componentes. Os resultados obtidos demonstram que ambos possuem um comportamento dinâmico, sendo que o condensador apresenta uma constante de tempo menor que o tanque de alimentação. O condensador, por rejeitar o calor necessário, acaba por perceber rapidamente as oscilações introduzidas no ciclo Além disto é apresentado um método para determinação de coeficientes globais de troca térmica em condensadores e regeneradores. Os resultados obtidos neste trabalho foram utilizados na modelagem desses componentes durante a construção de um simulador em regime permanente da UTE Presidente Médici - Fase B. Também foi realizada uma pesquisa sobre os diversos programas simuladores de plantas de potência, assim como uma discussão das vantagens introduzidas por simuladores dinâmicos e em regime permanente nas decisões de engenharia.

A avaliação da solicitação produzida por explosões, assim como da resposta de estruturas, são temas de muito interesse na engenharia atualmente, tanto pela quantidade crescente de acidentes relacionados com explosões quanto pelas ações terroristas muitas vezes vinculadas a estes tipos de ações. Neste contexto, o presente trabalho tem por objetivo explorar técnicas de análise tanto na modelagem da excitação quanto na resposta de estruturas consideradas como alvos. Para isto, são utilizadas metodologias de diferentes tipos: implementações baseadas em sistema de elementos finitos comerciais como Ansys [2000] e LS-Dyna [2003] e técnicas simplificativas que permitem realizar uma avaliação preliminar. As aplicações consideradas são indicadas a seguir: • Análise da Resposta de Estruturas Laminares Submetidas à Ação de Cargas Explosivas: determina-se a pressão produzida por explosivos sólidos a certa distância do epicentro, através de métodos simplificados, determinando a resposta esperada em placas retangulares; • Efeito da Pressão Interna em Vasos de Pressão (Extintores de Incêndio): comparando resultados numéricos e experimentais verifica-se a influência da pressão interna nas propriedades dinâmicas do sistema; • Estudo de Um Vaso Esférico de GLP Sob Ação de Uma Carga Explosiva: aplica-se a ação de uma onda explosiva produzida por um gás inflamável pesado sobre uma estrutura de vaso de pressão esférico com fluido e gás em seu interior, determinando sua resposta, avaliando também a influência de diferentes quantidades de líquido e pressão interna na resposta da estrutura; • Modelamento de uma Camada de Solo / Propagação das Ondas: verifica-se o comportamento da propagação de ondas em um meio elástico, comparando valores encontrados em testes experimentais. Analisa-se numericamente o efeito da inserção de uma valeta na atenuação de tais ondas; • Simulação Numérica Completa de uma Explosão: modela-se um semi-espaço submetido à ação de um explosivo sólido sobre sua superfície, avaliando os campos de pressão gerados. Ao final de cada aplicação são apresentadas conclusões parciais obtidas e as possibilidades de trabalhos futuros vislumbrados. Finalmente, conclui-se que as técnicas empregadas para as simulações são extremamente eficientes, considerando conhecidos todos os parâmetros envolvidos em cada modelo. Neste ponto é fundamental o trabalho do engenheiro, utilizando-se de seus conhecimentos técnicos para transformação do evento real em um modelo numérico, considerando e selecionando as simplificações necessárias.

Este trabalho tem por objetivo estudar a transferência de calor em tubos capilares cilíndricos utilizados na técnica de separação de moléculas denominada Eletroforese Capilar. Esta técnica é usada, por exemplo, na análise de biomoléculas e no sequenciamento de DNA, onde o controle da temperatura está diretamente ligado ao desempenho destes métodos e à qualidade dos resultados. Para empregar esta técnica, tensões elétricas da ordem de 20 kV são aplicadas entre as extremidades dos tubos capilares, que possuem normalmente 50 cm de comprimento, 350 µm de diâmetro externo e 50 µm de diâmetro interno, preenchidos por uma solução aquosa. Tais tensões geram uma corrente elétrica na solução, provocando aquecimento distribuído por Efeito Joule. Os tubos capilares são construídos em quartzo amorfo e protegidos por uma camada de material polimérico (poliimida). Para implementar o controle da temperatura, os tubos capilares são colocados em contato com um fluido de resfriamento. Num primeiro momento, os estudos são realizados por simulação numérica, empregando o Método dos Volumes Finitos em rotinas escritas em FORTRAN. São simulados casos onde os tubos são recobertos por camadas cilíndricas de materiais com uma condutividade térmica relativamente boa, com o objetivo de aumentar a superfície de troca de calor com o fluido de resfriamento. Como resultado, obtêm-se curvas da temperatura no centro dos tubos capilares em função do coeficiente de transferência de calor por convecção. Um caso de interesse é quando os tubos capilares são posicionados excentricamente ao recobrimento cilíndrico Num segundo momento, é utilizado o software de simulação numérica ANSYS CFX®, onde é simulado o resfriamento dos mesmos tubos capilares expostos a um escoamento transversal de ar a 15°C. Neste caso, também são aplicados os recobrimentos cilíndricos e, além disso, opta-se por simular o resfriamento de um arranjo de vários tubos (sistema multicapilar) dispostos entre placas de vidro, no formato de um sanduíche. Como resultados mais importantes salientam-se: a) o aumento do raio do recobrimento resulta no aumento da transferência de calor, fazendo com que a temperatura no núcleo do capilar fique estacionada em valores baixos que não comprometem as separações/análises; b) chegou-se a um valor de raio crítico da ordem de 10 mm para a condição de operação mais típicas na área da Eletroforese Capilar; c) as montagens com o tubo capilar concêntrico e excêntrico ao recobrimento não apresentam diferenças significativas no perfil de temperatura da solução tampão; e finalmente d) observa-se que o uso de duas placas de material dielétrico com os capilares posicionados em forma de sanduíche entre elas permite uma eficiente dissipação do calor gerado na solução tampão.

Esta tese apresenta um estudo do comportamento térmico de um coletor solar acumulador e desenvolve uma metodologia para medir a sua eficiência diária. O coletor solar acumulador está instalado na face norte do prédio de Energia Solar da UFRGS e possui cerca de 26 m2. É constituído de uma massa espessa de concreto com uma superfície absorvente feita de tijolos, possuindo uma cobertura dupla de vidros colocada de modo a deixar um espaço para a circulação de ar. Os raios solares atravessam a cobertura de vidro e aquecem a massa absorvente de tijolo, a qual aquece o ar que é introduzido no interior da construção por efeito de termossifão. Uma das principais características do coletor solar acumulador consiste no fato de que a resposta do coletor é defasada no tempo. Este fenômeno permite que o coletor entregue calor ao ambiente mesmo após o término da radiação solar. Essa defasagem dos picos de energia térmica ocorre devido ao baixo valor da difusividade térmica do concreto. Este trabalho foi dividido em duas etapas. A primeira etapa consistiu na montagem de um calorímetro para controle e monitoração das variáveis envolvidas. No interior do calorímetro foram instaladas 36 garrafas com água. As temperaturas dos conteúdos das garrafas, do coletor solar e as radiações envolvidas foram monitoradas através de 26 sensores de temperatura de CI, 8 sensores resistivos PT100 e dois sensores de radiação fotovoltaicos. Para obter as medidas dos sensores instalados foi feita a montagem de um sistema de aquisição de dados interfaceado a um microcomputador A segunda etapa consistiu na produção de um programa computacional, escrito em linguagem Fortran 90, para simular o comportamento térmico dos diversos elementos constituintes do coletor, determinar a potência térmica do coletor solar e sua eficiência diária. Para a simulação numérica do coletor solar acumulador, adotou-se um modelo simplificado bidimensional do mesmo. Foi integrada, através do Método dos Volumes Finitos, a equação de difusão de calor transiente em 2 dimensões. Na formulação das equações lineares optou-se pelo emprego das diferenças centrais no espaço e formulação explícita no tempo. Ao todo foram produzidas 4 malhas computacionais, com distintos refinamentos e foi realizado o estudo da estabilidade numérica das diversas malhas. Através da montagem experimental obtiveram-se várias características térmicas do comportamento do sistema, entre as quais, a transmitância da cobertura, curvas de temperatura do ar fornecido ao calorímetro e curva da eficiência diária do coletor solar . Através da simulação numérica foi possível determinar a potência térmica que o coletor entrega para o laboratório, a eficiência do coletor, os campos de temperatura e a vazão mássica nos diversos canais interiores do coletor solar.