Solução analítica da equação da energia estacionária e bidimensional para simulação de escoamento plenamente desenvolvido em placa paralela pelo método da giltt

Neste trabalho e apresentado um avanço na tecnica GILTT(Generalized Integral and Laplace Transform Technique) solucionando analiticamente um sistema de EDO's(Equações Diferenciais Ordinarias) de segunda ordem resultante da transformação pela GITT(Generalized Integral Transform Technique). Este tipo de problema usualmente aparece quando esta tecnica é aplicada na solução de problemas bidimensionais estacionários. A principal idéia consiste na redução de ordem do problema transformado em outro sistema de EDO's lineares de primeira ordem e a solução analítica deste problema, pela técnica da transformada de Laplace. Como exemplo de aplicação é resolvida a equação da energia linear bidimensional e estacionária. São apresentadas simulações numéricas e comparações com resultados disponíveis na literatura.

Métodos multigrid paralelos em malhas não estruturadas aplicados à simulação de problemas de dinâmica de fluidos computacional e transferência de calor

Fenômenos naturais, tecnológicos e industriais podem, em geral, ser modelados de modo acurado através de equações diferenciais parciais, definidas sobre domínios contínuos que necessitam ser discretizados para serem resolvidos. Dependendo do esquema de discretização utilizado, pode-se gerar sistemas de equações lineares. Esses sistemas são, de modo geral, esparsos e de grande porte, onde as incógnitas podem ser da ordem de milhares, ou até mesmo de milhões. Levando em consideração essas características, o emprego de métodos iterativos é o mais apropriado para a resolução dos sistemas gerados, devido principalmente a sua potencialidade quanto à otimização de armazenamento e eficiência computacional. Uma forma de incrementar o desempenho dos métodos iterativos é empregar uma técnica multigrid. Multigrid são uma classe de métodos que resolvem eficientemente um grande conjunto de equações algébricas através da aceleração da convergência de métodos iterativos. Considerando que a resolução de sistemas de equações de problemas realísticos pode requerer grande capacidade de processamento e de armazenamento, torna-se imprescindível o uso de ambientes computacionais de alto desempenho. Uma das abordagens encontradas na literatura técnica para a resolução de sistemas de equações em paralelo é aquela que emprega métodos de decomposição de domínio (MDDs). Os MDDs são baseados no particionamento do domínio computacional em subdomínios, de modo que a solução global do problema é obtida pela combinação apropriada das soluções obtidas em cada um dos subdomínios Assim, neste trabalho são disponibilizados diferentes métodos de resolução paralela baseado em decomposição de domínio, utilizando técnicas multigrid para a aceleração da solução de sistemas de equações lineares. Para cada método, são apresentados dois estudos de caso visando a validação das implementações. Os estudos de caso abordados são o problema da difusão de calor e o modelo de hidrodinâmica do modelo UnHIDRA. Os métodos implementados mostraram-se altamente paralelizáveis, apresentando bons ganhos de desempenho. Os métodos multigrid mostraram-se eficiente na aceleração dos métodos iterativos, já que métodos que utilizaram esta técnica apresentaram desempenho superior aos métodos que não utilizaram nenhum método de aceleração.

Proteção passiva de estruturas e sistemas elétricos contra incêndio aplicada a uma unidade de refino de petróleo

Nesta dissertação visa-se estudar e propor alternativas de solução para a proteção de estruturas e sistemas elétricos contra fogo numa unidade de craqueamento catalítico de uma refinaria de petróleo, por meio de proteção passiva. A proteção passiva tem por finalidade garantir a integridade das estruturas sujeitas a incêndio, durante um determinado período de tempo, para possibilitar, no caso da refinaria, a realização de procedimentos de parada da unidade de forma segura e controlar o incêndio a fim de diminuir a possibilidade de propagação do fogo para outras áreas. Com base em técnicas de análise de riscos fez-se a identificação de zonas potencialmente sujeitas a cenários de acidente envolvendo jato de fogo e/ou incêndio em poça. A delimitação das áreas onde haveria necessidade de proteção passiva foi realizada com base em modelos para jatos de fogo e incêndio em poça já estabelecidos na literatura. O dimensionamento da proteção passiva de estruturas e sistemas elétricos com o uso de diversos materiais usados comercialmente para este fim foi estimado com base em equações empíricas desenvolvidas por Jeanes, 1980, Stanzak, 1973 e PABCO, 1984, e, para alguns casos particulares foi feita uma verificação por solução numérica da equação da condução do calor em meio sólido.Assim, foram determinados quais os materiais mais adequados em cada caso de aplicação e qual a espessura em que deve ser aplicado para que a temperatura no elemento estrutural ou no sistema elétrico não atinja a sua determinada temperatura crítica em um período de tempo pré-determinado. Para os casos de elementos estruturais como colunas de sustentação da unidade de seção cilíndrica, o principal material para proteção passiva é a argamassa projetada e para perfil I, é o emprego de placas de gesso. Já para o caso de sistemas elétricos, podem ser utilizadas tanto tintas intumescentes quanto as mantas reforçadas com fibras minerais, esta escolha depende da geometria do sistema em que será empregado. Da comparação entre estes dois métodos pode-se concluir que o dimensionamento da proteção passiva fazendo o uso das correlações empíricas é menos conservativo que para o caso do uso da equação da difusão do calor resolvida por método numérico. Porém, os resultados diferem dentro de um limite considerado aceitável (em torno de 15%) levando-se em consideração os erros embutidos em cada método de cálculo. É importante mencionar que as correlações empíricas são de mais simples aplicação por possuir apenas operações matemáticas básicas. Usando as correlações empíricas para os perfis cilíndricos de aço (diâmetro de 0,1524 m e espessura de parede de 0,0254 m), a espessura de revestimento estimada com o uso das correlações empíricas necessária para garantir que a temperatura na interface entre os dois materiais não atinja 550°C em duas horas seria de 13,5 mm para argamassa projetada, 19,7 mm para vermiculita com silicato de sódio e 34,5 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Cabe ressaltar que com a realização desta dissertação busca-se uma integração entre o mestrado acadêmico e o meio empresarial com o desenvolvimento de trabalhos de natureza acadêmica que tenham aplicação direta na prática. Espera-se assim permitir que Universidade dê retorno à sociedade que a mantém e propiciar que setores da sociedade possam usufruir da capacidade disponível na academia.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Existência de soluções para equações elípticas semilineares com crescimento singular crítico

Neste trabalho estudamos uma equação diferencial parcial elíptica semilinear contendo uma singularidade e um termo de crescimento crítico. A existência de soluções depende da dimensão do espaço e do coeficiente da singularidade. Através da caracterização variacional e com o uso de seqüências de Palais-Smale provamos que o problema possui soluções não triviais.

Tratamento com calor no controle de Botryosphaeria dothidea (Moug.) Ces. & De Not.) em maçãs cv. Fuji

A necessidade de se prolongar o período pós-colheita fez com que a utilização de tratamentos em pós-colheita com o uso de defensivos agrícolas aumentasse nas últimas décadas. No entanto, a preocupação de consumidores quanto a resíduos as deficiências de controle observadas incentivaram a procura por tratamentos alternativos, entre os quais os tratamentos com calor vem sendo testados em frutos de várias espécies. No presente trabalho estudou-se o efeito do tratamento com calor no crescimento do micélio e germinação dos conídios de Botryosphaeria dothidea e, avaliou-se a viabilidade da termoterapia com calor, isolada ou em combinação com carbonato de sódio no controle da podridão branca e na manutenção da qualidade das maçãs ‘Fuji’ em diferentes estádios de maturação, mantidas ou não sob armazenamento refrigerado. Foi objetivo também verificar a expressão de genes que são regulados durante a interação maçã-Botryosphaeria dothidea, após tratamento com calor. Os experimentos foram conduzidos testando-se combinações de diferentes temperaturas e períodos de exposição ao calor sobre o patógeno B. dothidea e sobre frutos inoculados e não inoculados artificialmente com o patógeno. Os resultados obtidos sobre o patógeno mostraram que, independente da concentração de caldo batata-dextrose (BD) utilizada na suspensão, 45% dos conídios de B. dothidea germinaram nos frutos após 24h e 91% após 48h de incubação a 26°C. Constatou-se um aumento significativo da germinação dos conídios quando acrescentado 2,4% de caldo BD a suspensão. Sobre o efeito da temperatura no patógeno, as temperaturas inferiores a 62°C, em qualquer período de exposição não foram capazes de inibir o micélio do patógeno. Os conídios submetidos à temperatura de 58°C somente foram inativados quando o tratamento foi estabelecido por 60s. A termoterapia com calor associada ou não ao carbonato de sódio reduziu a incidência e severidade da doença podridão branca. A severidade da doença foi influenciada pelo estádio de maturação dos frutos verificando pelo aumento no número de lesões/fruto quanto mais avançada a maturação das maçãs. Os resultados obtidos no RT-PCR evidenciaram que os protocolos utilizados não produziram RNA de alta qualidade e desta forma inviabilizaram a análise da expressão diferencial dos genes na interação maçã-B.dothidea.

Estudo ergonômico sobre o calor

O calor, no ambiente onde o ser humano vive ou trabalha, atribui-se a: -fatores físicos: relação entre temperatura, radiação, umidade e movimento do ar; fatores humanos: o ser humano atua como fonte de energia através de seu metabolismo e atividade física. o calor caracterizado por um determinado meio ambiente é o resultado da atuação de diferentes variáveis, tais como: sistema de construção, situação geográfica do ambiente físico, climatização artificial, etc.; idade, sexo, capacidade física, estado de aclimatação, vestuário, tipo, carga e regime de trabalho, etc. A partir do momento em que o indivíduo for introduzido num determinado meio ambiente térmico, todos estes fatores vão influenciar a transmissão de calor entre ele e o ambiente. Na pretensão de haver equilíbrio térmico no meio ambiente quente, constata-se a necessidade de providenciar medidas de proteção a nível do sujeito e do ambiente, para que prevaleçam situações ambientais "confortáveis", ou pelo menos, "tolerantes". Desse modo, através da: definição das condições térmicas tolerantes e de conforto, parte-se para projetar meios ambientes de trabalho, que tornem praticáveis um isolamento térmico do calor exterior, assim como a perda de calor de dentro para fora. Atuando-se sobre variáveis individuais e ambientais estaremos incidindo diretamente sobre os meios de transmissão, procurando-se diminuir a quantidade de calor que o organismo produz e/ou recebe e aumentar a possibilidade de dissipá-lo.

Transferência de calor e controle de temperatura em tubos capilares utilizados em eletroforese capilar

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.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Otimização da produção de frango de corte em condições de estresse por calor

Foram realizados dois experimentos (EXP) com a finalidade de apresentar alternativas para reduzir os efeitos do estresse por calor (EPC), aumentando a produtividade em épocas quentes. O primeiro EXP teve como objetivos verificar o efeito de dietas com mais gordura (2,4 vs 4,0%) e menos proteína bruta (19,5 vs 18,5%) na metabolizabilidade, desempenho, rendimento de carcaça, parâmetros morfológicos (baço, bursa, coração, intestino e fígado), bioquímicos (proteínas totais, glicose, fructosamina, albumina e globulinas) e hematológicos (heterófilos. linfócitos, eosinófilos, basófilos, monócitos e relação heterófilo:linfócito) das aves aos 42 dias submetidas a EPC cíclico (25-32C). O efeito direto do EPC no desempenho e na metabolizabilidade da dieta, na situação de consumo pareado (comparado ao ambiente termoneutro) Também foi estudada. O segundo EXP teve como objetivos verificar o efeito de dietas suplementadas com vitaminas C e E (100UI de vitamina E/kg de ração e 300 ppm de vitamina C/kg de ração) e minerais orgânicos Zn e Se (40 ppm Zn e 0,3 ppm Se/kg de ração) no desempenho de frangos de corte submetidos a EPC e nos parâmetros morfológicos, bioquímicos e hematológicos das aves aos 35 dias. Foi observado que a queda de desempenho verificada no calor está relacionada principalmente com a diminuição no consumo e à redução na metabolizabilidade da matéria seca. A dieta com 1,6% a mais de gordura e 1% a menos proteína proporcionou às aves em EPC melhora na conversão alimentar (CA), mas não interferiu no rendimento de carcaça e cortes. Esta dieta também proporcionou um melhor peso relativo de bursa, a diminuição no número de linfócitos, heterófilos, relação H/L e monócitos das aves em EPC, indicando que as alterações metabólicas ocasionadas pelo estresse foram atenuadas pela dieta. No segundo EXP, a suplementação vitamínica e/ou mineral melhorou o desempenho das aves em função de um menor consumo que resultou em melhor CA, independentemente do ambiente. O EPC reduziu o peso absoluto e relativo dos órgãos linfóides, os valores de hemoglobina, aumentou o número de heterófilos e a relação H/L . As suplementações vitamínico e/ou mineral não influenciaram os parâmetros bioquímicos séricos e hematológicos dos frangos em EPC. A relação H/L foi um bom indicador de estresse, aumentando sempre que as aves estiveram estressadas, independentemente do tipo de estresse. O peso dos órgãos linfóides também pode ser considerado um bom indicador de EPC desde que não estejam sofrendo outro tipo de estresse.

Análise experimental e numérica do comportamento térmico de um coletor solar acumulador

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.

Utilização do método de Monte Carlo na resolução de problemas de transferência de calor por radiação em cavidades que contém meio participante

Um estudo foi realizado com a finalidade principal de investigar de que maneira a teoria de transferência de calor por radiação em meios participantes pode ser melhor aplicada na resolução de problemas de troca de energia no interior de cavidades que se assemelhem a fornalhas industriais. Verificou-se que as superfícies internas da maioria das fornalhas podem ser consideradas cinzentas e difusoras, o que facilita a aplicação de modelos de soma ponderada de gases cinzas. Sendo assim, as maiores dificuldades aparecem quando a geometria da cavidade é complexa. Com o método de Monte Carlo, é possível resolver problemas com características geométricas bastante complexas, bem como considerar efeitos espectrais e direcionais, quando necessário. Este método foi, então, aplicado, em conjunto com o método da zona, na elaboração de três códigos computacionais, que são capazes de determinar áreas de troca totais em cavidades onde as paredes e o gás são cinzas. Um deles pode ser aplicado a cavidades paralelepipédicas; outro, a cilíndricas. Tanto no primeiro, quanto no segundo, as dimensões podem assumir qualquer valor. Com o terceiro código computacional, é possível obter-se áreas de troca totais entre zonas de cavidades que têm geometria genérica formada a partir de combinações de cubos. Resultados obtidos através dos três códigos computacionais foram comparados com outros disponíveis na literatura.