Parâmetros bioecológicos de Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae) e modelagem da dinâmica espaço-temporal da sua interação com Spartocera dentiventris (Berg) (Hemiptera: Coreidae) através da simulação de múltiplos agentes

Estudos referentes à dinâmica espaço-temporal da interação hospedeiro-parasitóide são de fundamental importância para o desenvolvimento de técnicas de manejo de insetos em agroecossistemas. Uma recente abordagem têm sido o estudo de populações através da simulação computacional dos indivíduos que a compõem. Este método tem possibilitado a incorporação de importantes características como individualidade e estrutura espacial nos modelos teóricos. Neste sentido, o presente estudo visou a elaboração de um modelo baseado no comportamento individual para simular a dinâmica espaço-temporal da interação Spartocera dentiventris (Berg) (Hemiptera: Coreidae), inseto associado à cultura do fumo, e de seu inimigo natural, o parasitóide de ovos Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae). Para a elaboração do modelo, estudos sobre a bioecologia de G. gallardoi foram desenvolvidos, de modo a investigar (i) o efeito da temperatura no seu desenvolvimento e viabilidade, (ii) a mortalidade de imaturos em campo, (iii) os parâmetros reprodutivos e a longevidade dos adultos, (iv) os padrões de dispersão em campo e (v) a resposta funcional e a interferência mútua dos parasitóides em diferentes densidades de ovos de S. dentiventris. No intuito de avaliar a capacidade de predição do modelo elaborado, a dinâmica espaço-temporal de ovos sadios e parasitados de S. dentiventris foi acompanhada em um cultivo de fumo, estabelecido em área experimental do Departamento de Fitossanidade da Faculdade de Agronomia da UFRGS, Porto Alegre (30o01’S e 51o13’O), RS, Brasil A viabilidade do desenvolvimento dos parasitóides na faixa de 20 a 30oC não diferiu significativamente, alcançando 98,8%. O tempo de desenvolvimento ovo-adulto de machos e fêmeas foi inversamente proporcional ao aumento da temperatura. Os valores estimados para o limite térmico inferior de desenvolvimento e para a constante térmica foram 15,5oC e 185,19GD para machos e 15,6oC e 192,31GD para fêmeas. Diversos fatores afetam o sucesso de imaturos de G. gallardoi em campo, sendo que o malogro e a predação por sugadores são os principais responsáveis pela emergência de apenas 37,87% dos adultos. Foi observado um período médio de oviposição de 10,1 ± 1,74 dias, com o pico de oviposição no segundo dia, sendo depositados ao longo dos mesmos uma média de 67,5 ± 11,29 ovos por fêmea. Fêmeas de G. gallardoi foram significativamente mais longevas que os machos, vivendo, respectivamente, 13,7 ± 1,94 e 10,6 ± 1,78 dias. A razão sexual total observada foi de 0,79. A dispersão de G. gallardoi em cultivo de fumo não foi sensivelmente influenciada pelo vento, sendo estimada uma capacidade de locomoção diária média para fêmeas de no mínimo 7,62 m. O padrão de parasitismo de G. gallardoi ajustou-se perfeitamente à resposta funcional do tipo II (pseudo-r2= 0,99), sendo obtidos os valores de 0,0557 e 0,9989 h para os componentes a’ e Tm , respectivamente O aumento da densidade de parasitóides acarretou, de maneira geral, uma diminuição no número de ovos parasitados por parasitóide. Foi obtido para o índice de interferência mútua “m” o valor 0,626. A dinâmica temporal de S. dentiventris – G. gallardoi simulada ajustou-se muito bem aos dados obtidos em campo (r2= 0,82 para ovos sadios e r2= 0,72 para ovos parasitados). Da mesma forma, os arranjos espaciais dos ovos e da taxa de parasitismo observados em campo foram satisfatoriamente reproduzidos através das simulações computacionais. Os resultados indicam que G. gallardoi apresenta um bom potencial para controlar populações de S. dentiventris, seja pela sua viabilidade de desenvolvimento em uma larga faixa de temperatura ou pela sua capacidade reprodutiva e padrões de dispersão similares àqueles observados para outras espécies já utilizadas no controle biológico de insetos. Por outro lado, a sensibilidade às baixas temperaturas, a alta mortalidade de imaturos observada em campo e uma resposta funcional do tipo II podem ser fatores que se contrapõe ao sucesso do parasitóide no controle de populações de S. dentiventris. Considerando a independência entre os dados observados e os dados obtidos através de simulações, o modelo elaborado obteve pleno sucesso em simular a dinâmica da interação S. dentiventris – G. gallardoi a partir do comportamento individual dos diferentes agentes Duas importantes características associadas à estabilidade de sistemas hospedeiro-parasitóide foram observadas no estudo de campo e reproduzidas no modelo – um padrão de parasitismo inversamente dependente da densidade dos hospedeiros e uma resposta agregativa dos parasitóides em algumas regiões do sistema. Desta forma, futuros estudos podem avaliar, a partir de um enfoque baseado no indivíduo, quais mecanismos são responsáveis por tais comportamentos e qual é o seu efeito na dinâmica do sistema. Conclui-se, através deste trabalho, que o estudo da dinâmica de populações através de uma abordagem baseada no indivíduo, com base em modelos computacionais, pode ser uma alternativa viável para superar as limitações e a complexidade inerente ao uso de modelos populacionais analíticos. Desta forma, o maior realismo implícito nestes modelos os torna ferramentas muito úteis para diversas áreas da ecologia aplicada, particularmente para a entomologia agrícola.

Simulação e controle de um sistema de suspensão simplificado

As aplicações da mecânica vibratória vêm crescendo significativamente na análise de sistemas de suspensões e estruturas de veículos, dentre outras. Desta forma, o presente trabalho desenvolve técnicas para a simulação e o controle de uma suspensão de automóvel utilizando modelos dinâmicos com um, dois e três graus de liberdade. Na obtenção das equações do movimento para o sistema massa-mola-amortecedor, o modelo matemático utilizado tem como base a equação de Lagrange e a segunda lei de Newton, com condições iniciais apropriadas. A solução numérica destas equações é obtida através do método de Runge-Kutta de 4ª ordem, utilizando o software MATLAB. Para controlar as vibrações do sistema utilizou-se três métodos diferentes de controle: clássico, LQR e alocação de pólos. O sistema assim obtido satisfaz as condições de estabilidade e de desempenho e é factível para aplicações práticas, pois os resultados obtidos comparam adequadamente com dados analíticos, numéricos ou experimentais encontrados na literatura, indicando que técnicas de controle como o clássico podem ser simples e eficientes.

Critérios para locação amostral baseados em simulação estocástica

Em planejamento de lavra, freqüentemente as campanhas de sondagem são complementadas por falta de segurança nas estimativas. Na maioria dos casos, essas sondagens adicionais são locadas sem o uso de qualquer procedimento matemático que indique os locais mais carentes de informação. Essa dissertação propõe uma metodologia fundamentada em simulação estocástica para tratar esse problema. Com a utilização desse método, é possível que seja considerada a variabilidade entre as informações, além da configuração espacial apresentada. Essa metodologia permite gerar várias estimativas equiprováveis do atributo em questão, reproduzindo o histograma e o variograma das informações originais. A variabilidade apresentada pelos conjuntos de n valores equiprováveis, em cada bloco discretizado da jazida, é medida por meio do cálculo de índices de incerteza. Esses índices são calculados após o processo de simulação estocástica por meio da compilação das n estimativas realizadas para cada bloco. Os índices de incerteza considerados nessa dissertação são: (i) variância condicional; (ii) coeficiente de variação condicional e (iii) intervalo entre quartis. A partir da classificação desses índices de incerteza, são identificados os setores que apresentam as maiores variabilidades em suas estimativas. Com essa classificação é possível a confecção de mapas de incerteza que auxiliam na visualização dos setores/blocos mais variáveis da jazida Nesses setores, são adicionadas informações oriundas da escolha aleatória de um dos n valores equiprováveis disponíveis para cada bloco simulado. Por meio dos índices de incerteza e de acordo com a metodologia, são propostas informações adicionais na área considerada, que contribuem na diminuição das incertezas globais e locais. Essa metodologia usa técnicas extremamente robustas, com embasamento teórico muito consistente, atacando diretamente os setores/blocos que apresentam as maiores flutuabilidades nas estimativas. Um estudo de caso em uma jazida de carvão verifica a aplicabilidade do método. Esse estudo, considera a variabilidade nas estimativas de espessura de uma camada carvão.

Utilização do programa energyplus para simulação do conforto térmico em edificações populares em quatro capitais brasileiras

O principal objetivo deste trabalho consiste em utilizar o programa EnergyPlus para simular e analisar o nível de conforto térmico dos moradores de uma edificação térrea simples. Para tanto, escolheu-se e definiu-se uma casa popular padrão que será utilizada para as cidades de Belém, PA, Brasília, DF, São Paulo, SP e Recife, PE. Em relação a esta residência e aos seus moradores foram determinadas diversas características, entre elas, rotinas de ocupação, nível de atividades, fator de vestimenta e número de moradores. Ganhos internos (lâmpadas e equipamentos elétricos) e renovação do ar (infiltração e ventilação) foram, também, considerados neste trabalho. O conforto térmico dos ocupantes da habitação foi simulado e analisado para uma semana representante da época do ano com a temperatura mais elevada, e outra semana que representa o período de frio ou chuvas, ou seja, época do ano com temperaturas mais baixas. Para esta etapa necessitou-se a utilização de arquivos climáticos das quatro capitais brasileiras citadas acima. Os níveis de conforto térmico foram determinados e estudados conforme a metodologia de Fanger, que utiliza os seguintes parâmetros: temperatura, umidade relativa e velocidade do ar, temperatura média radiante do ambiente, vestimenta e metabolismo dos indivíduos. Este estudo foi realizado para um ambiente em evolução livre e analisou-se somente a sala e os quartos da residência. Também não foram considerados neste trabalho os equipamentos elétricos com o intuito de melhorar o conforto térmico como, por exemplo, ventiladores A ferramenta de simulação computacional utilizada para o cálculo do conforto térmico neste trabalho foi o programa EnergyPlus. Para a simulação do conforto necessitou-se calcular a infiltração e a ventilação na residência, e com este intuito utilizou-se uma sub-rotina deste programa baseada no modelo de fluxo de ar multizonal COMIS. Os resultados obtidos para o conforto térmico, através dos valores de PMV, nos dias analisados das quatro cidades em que foi realizado o estudo mostraram-se bastante coerentes com as condições externas e internas da residência, indicando a boa capacidade do programa EnergyPlus nestes casos simulados. Calculado e analisado o conforto térmico das pessoas, melhorias nas estruturas da edificação foram sugeridas a fim da obtenção de condições otimizadas para os ocupantes desta habitação.

Metabolismo da homocisteína e defeitos do tubo neural : um estudo bioquímico e molecular no sul do Brasil

Os defeitos de fechamento de tubo neural constituem uma das malformações mais freqüentes na espécie humana, apresentando alta morbi-mortalidade. Sua etiologia é considerada multifatorial, estando envolvidos fatores genéticos e ambientais. Estes fatores estão relacionados principalmente com o metabolismo da homocisteína. Realizamos um estudo de caso-controle com o objetivo de estudar os fatores bioquímicos e genéticos relacionados ao DTN na nossa população. Em pares de afetados com DTN e suas mães e pares de pacientes normais e suas mães foram avaliados dosagem de folato, vitamina B12, homocisteína e polimorfismos da enzima metileno tetraidrofolato redutase (MTHFR), C677T e A1298C. A dosagem de folato nos casos foi 11,37 ng/mL(±6,72) e nos controles 5,64 ng/mL(±4,16) (p<0,001). O folato sérico das mães foi 7,27 ng/mL (±4,48) e 3,90 ng/mL (±1,77) nas mães controles (p<0,001). A média de dosagem de vitamina B12 foi de 641,88 pg/mL ((±262,21) nos casos e 743,27 pg/mL (±433,52) nos controles (p= 0,205). A média de dosagem de vitamina B12 nas mães dos casos foi 354,75 pg/mL (±142,06) e 465,25 pg/mL (±194,91) nas mães controles (p=0,004). O nível de homocisteína plasmático médio foi 6,89 μmol/L(±4,48) para os casos e 5,41 μmol/L (±2,55) para os controles (p=0,099). Nas mães dos casos a dosagem média de homocisteína foi 7,23 μmol/L (±2,64) e 7,00 μmol/L (±2,24) nas mães controles (p=0,666). Não houve diferença entre a freqüência dos genótipos C677T e A1298C da MTHFR nos casos e controles e suas mães. Para o polimorfismo C677T as freqüências dos alelo C e T foram respectivamente 0,6585 e 0,3414 nos pacientes com DTN; 0,6590 e 0,3410 nos controles; 0,6460 e 0,3540 nas mães dos casos e 0,6136 e 0,3860 nas mães controles. Para o polimorfismo A1298C as freqüências dos alelos A e C foram respectivamente 0,7436 e 0,2564 nos pacientes com DTN; 0,7610 e 0,2390 nos controles; 0,8055 e 0,1945 nas mães dos casos e 0,8065 e 0,1935 nas mães controles. Identificamos que indivíduos homozigotos 677TT apresentam um maior nível de homocisteína e este é inversamente relacionado com os níveis de vitamina B12. Estes achados sugerem que uma alteração metabólica relacionada ao metabolismo da homocisteína e principalmente devido à diminuição da vitamina B12 seja um fator de risco para DTN na nossa população.

Desenvolvimento de um programa de simulação computacional de sistemas de aquecimento solar para água

Esta Tese apresenta uma análise do comportamento térmico de um sistema de aquecimento solar operando por termossifão. Neste tipo de sistema o fluido no coletor solar é circulado por convecção natural, que acontece devido à diferença de massa específica da água ao longo circuito. Nestes sistemas a vazão mássica varia ao longo do dia e do ano, dependendo, dentre outros fatores, da irradiância solar absorvida, do perfil de temperaturas da água no sistema, da geometria, do volume e do perfil de demanda de água quente. Para uma avaliação detalhada do comportamento térmico de aquecedores solares operando por termossifão foram realizados ensaios experimentais e cálculos teóricos. Os resultados dos experimentos concordaram com aqueles apresentados na literatura e sua análise fundamentou o desenvolvimento do aplicativo TermoSim, um programa de simulação computacional do comportamento térmico de sistemas de aquecimento de água com energia solar. O tratamento matemático adotado no TermoSim compreende a modelagem dos coletores solares de acordo com a teoria de Hottel-Bliss-Whillier. O reservatório térmico é modelado com estratificação térmica, convecção e condução entre as camadas. A vazão mássica é obtida a partir do balanço da quantidade de movimento no circuito. Os modelos matemáticos empregados na construção do aplicativo TermoSim foram validados através do confronto dos resultados simulados com medidas experimentais. Foi demonstrado que a utilização destes modelos é adequada e permite reproduzir com precisão o comportamento térmico dos coletores solares e do reservatório térmico. Além do programa TermoSim, foi também desenvolvido o programa TermoDim, que é uma ferramenta para o dimensionamento de sistemas de aquecimento solar, que requer apenas o conhecimento dos parâmetros geométricos do sistema, dados meteorológicos em média mensal e informação a respeito do volume de demanda. O TermoDim é apropriado para estimar o desempenho de aquecedores solares operando por termossifão com tanques verticais e horizontais. O método de dimensionamento do TermoDim é baseado na correlação para a eficiência média mensal obtida neste trabalho a partir de um grande número de simulações.

Modelagem e simulação de um reator trickle-bed para hidrogenação de propeno

Neste trabalho foi modelado um reator catalítico heterogêneo, pressurizado, adiabático e provido de reciclo. O leito deste reator está dividido em duas partes, sendo que a primeira opera como um reator trickle-bed, enquanto a segunda opera como um reator de leito fixo, onde os regentes estão em fase líquida. Um Sistema de Hidrogenação de Propeno em operação no Pólo Petroquímico do Sul foi utilizado para fornecer os dados necessários à validação do modelo proposto. Neste sistema, a reação de hidrogenação de propeno, em fase líquida, ocorre sobre um catalisador níquel/óxido de níquel. Uma equação, que representa a taxa da reação no intervalo de temperatura de interesse para o sistema, foi desenvolvida utilizando-se dados experimentais obtidos em um reator batelada. O comportamento termodinâmico foi representado através de dados de equilíbrio do sistema binário hidrogênio-propano. Considerou-se que a fase líquida escoa em plug-flow e a fase gasosa, devido à existência de um selo líquido no interior do reator, foi considerada estagnada.

Defeitos responsáveis pela isolação de GaAs irradiado com prótons

A Espectroscopia de Transientes de Níveis Profundos (DLTS – Deep Level Transient Spectroscopy) foi, detalhadamente, descrita e analisada. O processo de isolação por implantação em GaAs foi estudado. Sua dependência com a sub-rede, do As ou do Ga, em que o dopante é ativado foi investigada para material tipo-p. Semelhantes doses de implantação de prótons foram necessárias para se tornar semi-isolantes camadas de GaAs dopadas com C ou com Mg possuindo a mesma concentração de pico de lacunas livres. A estabilidade térmica da isolação nestas amostras foi medida. Diferenças no comportamento de recozimento destas apontaram a formação, provavelmente durante a referida etapa térmica, de uma estrutura diferente de defeitos em cada caso. Medidas de DLTS foram realizadas em amostras de GaAs tipo-n e tipo-p implantadas com prótons de 600 keV. A estrutura de picos observada apresentou, além de boa parte dos defeitos introduzidos para o caso de irradiação com elétrons, defeitos mais complexos. Um novo nível, com energia superior em ~0,64 eV ao valor correspondente ao topo da banda de valência, foi identificado nos espectros medidos em material tipo-p. A variação da concentração dos centros de captura introduzidos com diferentes etapas de recozimento foi estudada e comparada com o comportamento previamente observado para a resistência de folha em camadas de GaAs implantadas com prótons. Simulações foram feitas, indicando que a interpretação adotada anteriormente, associando o processo de isolação diretamente à formação de defeitos relacionados a anti-sítios, pode não estar completa.

Simulação da irrigação por inundação e da drenagem nos solos de várzea do Rio Grande do Sul

Os solos de várzea são responsáveis por uma parcela significativa da economia do Rio Grande do Sul. As principais atividades de exploração desses solos são o cultivo do arroz irrigado por inundação e a pecuária. Tendo em vista o método de irrigação adotado, o cultivo do arroz é acompanhado de grandes consumos de água e energia. As más condições de drenagem, associadas à uma série de outros fatores, praticamente limitam a atividade agrícola à essa cultura, e fazem com que a pecuária seja quase uma imposição. Tendo esse quadro como referência foram definidos os objetivos desse trabalho, através dos quais buscaram-se algumas respostas que possam auxiliar na definição de um modelo de exploração mais adequado ao grande potencial que esses solos apresentam. Inicialmente o modelo agroidrológico SWAP foi adaptado às particularidades da irrigação por inundação e dos sistemas de drenagem que têm seus drenos implantados com subsoladores. Criou-se então uma nova versão, denominada SWAP_INU. Na seqüência o modelo foi validado a partir da comparação de seus resultados com observações feitas em áreas da Estação Experimental do Arroz do IRGA. Nessas áreas os perfis de solo têm como característica a presença de horizonte B superficial e impermeável. Uma série de simulações permitiram que fossem confrontadas as lâminas de irrigação e profundidades do lençol freático registradas nos experimentos e calculadas pelo modelo. No caso específico da irrigação, as lâminas observadas e calculadas também foram comparadas com as recomendadas pelo IRGA, quando pode-se constatar que o limite inferior do intervalo recomendado é compatível com as observações e cálculos feitos. Além disso, verificou-se que nesse tipo de perfil as características do horizonte A têm pouca importância sobre os resultados finais. Em outra série de simulações procurou-se relacionar a drenagem com a possibilidade de mecanização. Foram criados cenários com diferentes condições de drenagem superficial e subsuperficial, além da presença de drenos livres. Os resultados mostraram uma dependência muito grande entre os sistemas convencionais, e um desempenho bastante satisfatório da drenagem não convencional. Por fim, os resultados de um modelo simplificado, IRRIGA, foram confrontados com os do modelo SWAP_INU. A proximidade entre os valores calculados indicou um ótimo desempenho do modelo simplificado na estimativa da necessidade total de água para a irrigação por inundação.

Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.

Caracterização metalúrgica do aço SAE 1141 microligado ao nióbio por simulação física

Este trabalho analisa, através do emprego de simulador físico, os efeitos da adição do nióbio no aço SAE 1141 e seu comportamento como refinador de grão e endurecedor por precipitação, nas altas temperaturas de forjamento e tempos de reaquecimento aplicados industrialmente. Analisa também qual a influência das deformações e velocidades de resfriamento sobre o tamanho de grão austenítico, microestrutura, dureza, cinética da precipitação e propriedades finais. Na definição dos parâmetros de simulação foram observadas rotinas de produção de peças forjadas comercialmente, a literatura técnica, e possíveis alterações que otimizassem as propriedades do aço em estudo. Para tanto foram definidas as temperaturas de forjamento 1100, 1200 e 1250°C, a taxa de deformação de 1/s, as deformações de 30, 50 e 70%, e as velocidades de resfriamento de 36, 60, 100 e 150°C/min. Os resultados indicaram que a velocidade de resfriamento tem papel preponderante na microestrutura final, que variou de um agregado de ferrita e perlita para bainita (às vezes com alguma martensita) a medida em que a velocidade de resfriamento aumenta. A simulação térmica indicou que tamanho de grão austenítico aumenta a medida em que se eleva a temperatura de reaquecimento e que existe apenas um pequeno efeito ancorador de grão do nióbio nesse quesito, quando não estiver presente alguma deformação Foi observado, através da simulação física que, com a presença de deformação, principalmente para as temperaturas menores de reaquecimento, a precipitação induzida por deformação parece ancorar o crescimento dos grãos recristalizados, resultando num tamanho de grão austenítico menor. Esse trabalho constitui parte integrante do projeto de desenvolvimento da tecnologia dos aços microligados, num esforço conjunto do Laboratório de Metalurgia Física – LAMEF, através do Grupo de Desenvolvimento de Aços Microligados, da Universidade Federal do Rio Grande do Sul, em parceria com a Aços Finos Piratini e DANA – Albarus, com o apoio da CAPES.

Simulação computacional por dinâmica molecular de propriedades de equilíbrio e espectroscópicas de misturas líquidas CS2/C6H6

Propriedades de equilíbrio e espectroscópicas são determinadas para a mistura líquida CS2/C6H6 na temperatura de 298K e nas frações molares de benzeno iguais a 0,25, 0,50 e 0,75, usando simulação computacional por Dinâmica Molecular. As interações intermoleculares são descritas em relaçãao a todos os átomos por um potencial Lennard-Jones (12/6), com as interações eletrostáticas representadas por quadrupolos pontuais, localizados nos centros de massa das mol´eculas. Um segundo potencial também é considerado, onde a distribuição de cargas na molécula é usada para descrever as interações eletrostáticas. Dados termodinâmicos e propriedades estruturais, descritas através de funções de distribuição radial e de correlação angular, são calculados. Os resultados obtidos com os dois potenciais mostram boa concordância com os dados experimentais. A análise das frações molares locais e das energias internas indica que a mistura líquida CS2/C6H6 apresenta comportamento aproximadamente ideal. As funções de distribuição radial e de correlação angular evidenciam uma configuração preferencialmente paralela entre as moléculas de CS2 e C6H6 a distâncias menores em relação a distância do máximo da função de distribução , correspondente a primeira camada de coordenação . Nesta região, a configuração perpendicularé fracamente favorecida em relação as demais orientações. Funções de correlação temporal da polarizabilidade coletiva e da polarizabilidade molecular, correspondentes aos espalhamentos Rayleigh e Raman, respectivamente, são calculadas a partir do modelo de interação dipolo induzido por dipolo de 1a ordem. As funções de correlação são divididas nas contribuições orientacional, induzida por colisão e cruzada, e nas contribuições dos componentes. Propriedades de equilíbrio como anisotropia efetiva e intensidades integradas são determinadas para as misturas e líquidos puros. Os resultados indicam que a correlação orientacional entre as moléculas de CS2 é a maior responsável pelos espectros Rayleigh e espectros Raman deste componente. A maior participação do benzeno ocorre no espalhamento Raman com uma maior correlação na contribuição induzida por colisão. A separação na escala temporal entre a dinâmica reorientacional e induzida por colisão determinada na simulação é menor em relação ao experimento. Tempos de correlação da contribuição induzida por colisão são maiores nas misturas em relação aos líquidos puros no espalhamento Rayleigh, conforme o experimento, e estão relacionados a contribuições significativas das funções de correlação entre os componentes CS2 e benzeno a tempo longo. Funções de correlação temporal do momento dipolar coletivo são determinadas a partir do modelo dipolo induzido por quadrupolo e analisadas em termos de contribuições isotrópicas, anisotrópicas e cruzadas, como também de contribuições dos componentes. A relação entre estes termos é concordante com os parâmetros moleculares de anisotropia e quadrupolo dos componentes. Os resultados da simulação indicam absorbância de excesso para as misturas em uma ordem semelhante a do experimento. O tempo de correlação maior na mistura de fração molar de benzeno igual a 0,25, também encontrado no experimento, está relacionado a elevada contribuição a tempo longo da função de correlação entre dipolos induzidos em moléculas de CS2 por moléculas de C6H6. Os espectros calculados na simulação e os espectros experimentais do espalhamento de luz despolarizado Rayleigh apresentam boa concordância, validando os modelos de potenciais de interação e de indução usados. No entanto, diferenças são encontradas no infravermelho longínquo. Modi cações no modelo simplificado de indução são propostas.

Geração de modelos de co-simulação distribuída para a arquitetura DCB

O aumento na complexidade dos sistemas embarcados, compostos por partes de hardware e software, aliado às pressões do mercado que exige novos produtos em prazos cada vez menores, tem levado projetistas a considerar a possibilidade de construir sistemas a partir da integração de componentes já existentes e previamente validados. Esses componentes podem ter sido desenvolvidos por diferentes equipes ou por terceiros e muitas vezes são projetados utilizando diferentes metodologias, linguagens e/ou níveis de abstração. Essa heterogeneidade torna complexo o processo de integração e validação de componentes, que normalmente é realizado através de simulação. O presente trabalho especifica mecanismos genéricos e extensíveis que oferecem suporte à cooperação entre componentes heterogêneos em um ambiente de simulação distribuída, sem impor padrões proprietários para formatos de dados e para a descrição do comportamento e interface dos componentes. Esses mecanismos são baseados na arquitetura DCB (Distributed Co-Simulation Backbone), voltada para co-simulação distribuída e heterogênea e inspirada nos conceitos de federado (componente de simulação) e federação (conjunto de componentes) que são definidos pelo HLA (High Level Architecture), um padrão de interoperabilidade para simulações distribuídas. Para dar suporte à co-simulação distribuída e heterogênea, esse trabalho descreve mecanismos que são responsáveis pelas tarefas de cooperação e distribuição, chamados de embaixadores, assim como o mecanismo gateway, que é responsável pela interoperabilidade entre linguagens e conversão de tipos de dados. Também é apresentada uma ferramenta de suporte à geração das interfaces de co-simulação, que são constituídas de dois embaixadores configuráveis e um gateway para cada federado, gerado a partir de templates pré-definidos.

Modelagem e simulação da câmara de reação de uma caldeira a carvão pulverizado

O presente trabalho é parte integrante de um projeto multidisciplinar cujo objetivo é construir uma ferramenta computacional capaz de simular a operação de uma Usina Termoelétrica tratando, especificamente, da câmara de reação (local da caldeira onde se desenvolve a queima do combustível). Foi desenvolvido um modelo matemático dinâmico, pseudo-homogêneo, unidimensional e, com propriedades físicas variáveis, com o objetivo de descrever a conversão do calor gerado pela combust em energia térmica, que está baseado em balanços de massa e energia. Este modelo foi implementado no simulador gPROMS e seus parâmetros foram estimados utilizando o módulo gEST deste software (método da máxima verossimilhança), com base em dados de operação da Usina e dados cinéticos de literatura. Para a simulação dinâmica da câmara de reação foi utilizado um método de integração implícita de passos múltiplos com ordem e passo variáveis. A câmara de combustão foi dividida em três zonas, sendo realizados, em cada zona, os balanços de massa e energia, adotando-se a hipótese de mistura perfeita. O modelo desenvolvido prevê o perfil de temperatura ao longo da câmara, as correntes de saída e a composição dos gases de combustão em cada zona; incluindo uma previsão aproximada da concentração dos poluentes formados, principalmente monóxido de nitrogênio e dióxido de enxofre. Os resultados apresentados mostraram que o modelo prediz satisfatoriamente a composição do gás de combustão na saída da câmara de reação e que a divisão da câmara em diferentes zonas de combustão é apropriada para avaliar a geração e o aproveitamento energético na caldeira.