Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Desenvolvimento de um conjunto de modelagem de alta resolução para os métodos eletromagnéticos indutivos

A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre.

Modelagem paramétrica linear aplicada à faixa de sinal de televisão digital para cidades densamente arborizadas

Este trabalho apresenta uma modelagem paramétrica (auto-regressiva) linear aplicável a estudos de propagação de televisão digital e telefonia celular para cidades densamente arborizadas. A modelagem proposta apresenta um forte embasamento estatístico e depende apenas de dados provenientes de medição, no caso dados relativos a potência recebida e o valor de PSNR (Peak Signal-to-Noise Ratio). Um algoritmo genético é utilizado no cálculo dos parâmetros de ajuste do modelo a um conjunto de dados. O trabalho foi realizado na faixa de televisão digital e foram analisadas duas variáveis: a potência recebida do sinal e o valor de PSNR. Foram executadas campanhas de medição na cidade de Belém. Nestas medições foram coletados dados de potência e gravados vídeos da programação diária de uma emissora de televisão. Os resultados podem ser aplicados no planejamento de serviços de telecomunicações.

Ensaios sobre a contribuição do ar nos resultados da modelagem do modo transverso-elétrico do método magnetotelúrico

Apresentamos a solução numérica do modelo da falha infinita para o modo TE utilizando o método dos elementos finitos para investigar o efeito do ar na solução TE. A comparação de nossa solução para o campo elétrico normalizado com a solução de Weaver evidenciou uma discrepância esperada devido ao efeito do ar, que foi intencionalmente negligenciado por Weaver, com o intento de facilitar a solução do problema. O problema analítico, levando em conta a presença do ar, é substancialmente mais difícil. Uma solução do problema real, parcialmente analítica e parcialmente numérica, foi apresentada por Sampaio. Devido a uma discrepância inesperada entre a nossa solução por elementos finitos e a solução exata de Sampaio, recalculamos a sua solução. Iterando a parte numérica da solução de Sampaio até ordens mais altas pudemos verificar que esta apresenta comportamento divergente. Para examinar o efeito do ar na resistividade aparente, comparamos três perfis de resistividade aparente obtidos a partir de nossa solução por elementos finitos com os mesmos perfis obtidos pela solução de Weaver. A discrepância encontrada entre estas soluções foi muito pequena, situando-se na faixa de erro instrumental de campo.

Estudo de geração de campo uniforme para modelagem eletromagnética

Para estudar os problemas de prospecção geofísica eletromagnética através de modelagem analógica, as condições em escala natural são representadas, no laboratório, em escala reduzida de acordo com a teoria da similitude. Portanto, para investigar os problemas de técnicas VLF, AFMAG e MT, frequentemente é necessário criar campo uniforme no arranjo experimental. Vários sistemas físicos para geração de campos uniformes são analisados teoricamente nesta tese. Os sistemas estudados aqui são a bobina circular, bobina de Helmholtz, solenóide, um plano de corrente e dois planos paralelos de correntes. As equações analíticas foram obtidas para campo magnético num ponto do espaço e subsequentemente as condições de campo uniforme. Nos casos em que as condições para o campo uniforme não puderam ser obtidas analiticamente, a porcentagem de desvio do campo em relação a um ponto pré-selecionado foi calculada. Contudo, os mapas de campo magnético, assim como o mapa de porcentagem de desvio, estão presentes para todos os sistemas estudados aqui. Também, foram calculados as áreas e os volumes espaciais de vários desvios de porcentagem do campo uniforme. Um estudo comparativo desses sistemas mostra que o solenóide é a maneira mais eficiente para criar um campo uniforme, seguido pelo sistema de bobinas de Helmholtz. Porém, o campo criado em um solenóide está em um espaço fechado onde é difícil colocar modelos e substituí-los para executar experimentos. Portanto, recomenda-se o uso de bobinas de Helmholtz para criar um campo uniforme. Este último sistema fornece campo uniforme com espaço aberto suficiente, o que facilita o experimento.

Estudos dos efeitos do manto de intemperismo sobre anomalias de VLF geradas por corpos inclinados: modelagem numérica através de programa eGs

As referências apresentadas e utilizadas através de métodos indutivos de propagação eletromagnética E.M., são utilizadas na prospecção de minerais, entre outros, haja vista que em regiões tropicais existem camadas superficiais condutivas (manto intemperismo) que levam a influenciar, evidentemente, na interpretação do resultado real. Em nossa região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo. O objetivo deste trabalho é estudar os efeitos dos mantos sobre anomalias VLF, devidos a corpos tabulares inclinados em contato ou não com o manto (manto galvânico e manto indutivo). Deste modo, foram realizadas umas séries de experimentos em modelagens numéricas através da eGs, considerando-se diferentes parâmetros de resposta para o manto e para o corpo. Os trabalhos de modelagem numérica, bem como seus resultados foram apresentados por diversas situações, como: variação da profundidade do topo do condutor, variação do mergulho do condutor, variação da condutância do manto, entre outros. Para interpretação dos resultados, foi gerado um conjunto de curvas considerando-se valores pico-a-pico das anomalias de Tilt Angle e de Elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. A medida em que se aumenta a profundidade, o pico-a-pico e a forma do pico da anomalia do Tilt Angle e da Elipsidade diminuem na presença do manto, ou na ausência do mesmo, fazendo com que este se afaste do ponto de “Cross-Over”. Para um condutor da mesma espessura, o aumento da condutância causa um aumento nas anomalias de Tilt Angle e da Elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de Tilt Angle e Elipsidade. O aumento da condutância do manto de intemperismo causa uma diminuição nas anomalias de Tilt Angle e de Elipsidade. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de condutância do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e, tendendo, a ser menos condutivo.

O uso de Thin-Plate Splines na transformação de coordenadas com modelagem de distorções entre realizações de referenciais geodésicos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação computacional da adsorção dos poluentes benzeno, tolueno e p-xileno sobre carvão ativado

Os maiores problemas de contaminação de aquíferos e solos são atribuídos aos hidrocarbonetos monoaromáticos, que são os constituintes mais solúveis e mais móveis da fração de algumas substâncias, como por exemplo, da gasolina. Para a remoção destes contaminantes, a adsorção por carvão ativado é o método mais utilizado, pois o carvão apresenta uma habilidade significativa para adsorver componentes orgânicos de baixo peso molecular, como o benzeno, o tolueno e o p-xileno. Neste trabalho, verificou-se a adsorção dos mesmos sobre carvão ativado via simulação computacional. Como base, utilizou-se o modelo postulado de carvão preparado por Bourke et al. (2007). Várias etapas foram concluídas desde os modelos das estruturas do carvão e dos poluentes até as simulações de dinâmica molecular. Para a análise conformacional da estrutura do carvão, foi utilizado o método semi-empírico PM3 e para o processo de dinâmica, o campo de força AMBER FF99SB. A estrutura passou por um aquecimento, à pressão constante, até alcançar uma temperatura final de 298K (25ºC), sendo suas informações coletadas a cada 50ps. Posteriormente, a estrutura foi submetida a equilíbrio de sistema, à temperatura constante de 298K (25ºC), por 500ps para então suas informações serem analisadas. Por fim, o sistema foi então submetido à dinâmica molecular durante 15 ns. Após análise dos resultados, constatou-se que os grupos éter, lactona e carbonila (cetona) presentes na estrutura de carvão ativado conferem caráter ácido à mesma e devido a isto e à sua consequente carga superficial negativa, a adsorção tornou-se viável uma vez que os poluentes apresentavam carga superficial positiva, o que corrobora o entendimento que já se tem a respeito desse tipo de fenômeno.

Equilíbrio líquido-líquido dos sistemas biodiesel de castanha do Brasil + metanol + (glicerina ou água): determinação experimental e modelagem termodinâmica

O biodiesel é definido como sendo uma mistura de monoésteres de ácidos graxos derivados de gorduras animal ou óleos vegetais, obtido por meio do processo de transesterificação com alcoóis de cadeia curta. Durante a sua produção é utilizada uma quantidade em excesso de álcool e um catalisador para favorecer a reação de formação dos ésteres; desta forma a corrente de saída do reator contém o catalisador, o álcool que não reagiu e os produtos da transesterificação, biodiesel e glicerina, formando um sistema bifásico; dependendo do grau de solubilidade desses compostos, pode haver quantidades de biodiesel na fase rica em glicerina e quantidades de glicerina na fase rica em biodiesel. Durante o processo de purificação do biodiesel é necessário executar uma lavagem com água, para promover a retirada do catalisador e impurezas do produto de interesse. Devido à quantidade de compostos envolvidos na produção e purificação é essencial conhecer os dados de equilíbrio líquido-líquido para poder predizer as proporções em que os compostos coexistem e, posteriormente, proceder com a purificação sob condições adequadas para obtenção do biodiesel com maior rendimento e auxiliar no projeto do reator e sistemas de separação. Os objetivos desse trabalho foram produzir o biodiesel a partir de uma planta oleaginosa nativa da região, a castanha do Brasil (Bertholletia excelsa H. B. K.) e, determinar dados de equilíbrio líquido-líquido (ELL) para os sistemas ternários: biodiesel de castanha do Brasil + metanol + glicerina e biodiesel de castanha do Brasil + metanol + água nas temperaturas de 30°C e 50°C. O biodiesel produzido foi inicialmente caracterizado segundo as normas da Agência Nacional de Petróleo (ANP) utilizando métodos físico-químicos. Os dados obtidos foram posteriormente correlacionados com os modelos termodinâmicos NRTL e UNIQUAC para o calculo do coeficiente de atividade de cada componente na fase liquida, com estimativa de novos parâmetros de interação energética. Os resultados obtidos com a modelagem foram satisfatórios, e foi observado que o modelo NRTL representou melhor os dados experimentais.

Modelagem das exigências de lisina para aves de postura em crescimento

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Avaliação através do desempenho térmico de edificação verticalizada em Belém-PA

A presente dissertação de mestrado estuda o desempenho térmico de materiais de fachada em uma edificação verticalizada na cidade de Belém – PA, Brasil. Este trabalho tem por objetivo contribuir com estudo do desempenho térmico e eficiência de materiais utilizados na envoltória das edificações residenciais verticalizadas, entre os elementos opacos, foram selecionados dois tipos de materiais, bastante utilizados atualmente na construção civil, a pintura e revestimento cerâmico. Também estuda-se a influência da altura da edificação sobre seu desempenho térmico e eficiência energética. O estudo foi desenvolvido com o uso de modelagem em código computacional no programa EnergyPlus que permite avaliar o desempenho térmico da edificação, a influência da altura na eficiência energética da mesma. Compara-se o desempenho térmico dos dois materiais selecionados, verificou-se que o prédio com revestimento cerâmico apresentou a temperatura interna média anual de 0,42 º C acima da temperatura referente ao prédio pintado, portanto com menor desempenho térmico. A influência da altura do pavimento não produziu uma variação significativa na temperatura interna, pois com seu acréscimo, ocorreu pequena diminuição da temperatura interna, com gradiente térmico da ordem de -0,07ºC /m em decorrência da velocidade do vento que é maior quando mais alto é o pavimento. Considera-se também para este estudo a ventilação natural como estratégia passiva para a edificação em estudo, que segundo a análise da Carta Bioclimática de Givoni contribuiria positivamente para a condição de temperatura interna e umidade do ar, e consequentemente do usuário.

Estratégia de identificação paramétrica aplicada à modelagem fenomenológica de um sistema do tipo correia transportadora industrial para fins de detecção de faltas

Os sistemas do tipo correia transportadora, são sistemas essenciais para grandes empresas, embora este equipamento apresente um elevado grau de criticidade, de modo que uma parada não planejada do mesmo, pode gerar perdas imensas ou até mesmo a parada de todo o processo produtivo. Tendo em vista a criticidade do equipamento, torna-se necessário realizar a monitoração adequada do mesmo e detectar com maior antecedência possível a ocorrência de alguma falta ocorrida no sistema. Objetivando reduzir as paradas não planejadas, investiga-se nesta dissertação a modelagem de um sistema do tipo correia transportadora com a finalidade de uso no monitoramento e diagnóstico de faltas neste tipo de sistema. Primeiramente é discutido um modelo fenomenológico do processo, o qual é baseado na aplicação das leis da mecânica e considerando-se os diversos tipos de força de oposição ao movimento de uma correia transportadora. Os principais parâmetros do transportador foram estimados utilizando-se técnicas de identificação baseadas em Mínimos Quadrados Não Recursivos. Em seguida, foi desenvolvido e implementado um algoritmo de detecção de faltas utilizando a teoria de análise intervalar, de modo que possibilite detectar condições inadequadas de funcionamento. Com o intuito de avaliar o desempenho do algoritmo proposto foi projetado e construído um protótipo que emula as condições operacionais típicas de um sistema real do tipo correia transportadora. Os resultados obtidos experimentalmente, confirmam o bom desempenho da metodologia proposta.