Análise de métodos numéricos para precificação de opções

Analisa e compara, sob as óticas de acurácia e custo, os principais métodos numéricos (simulação de Monte carlo, métodos de diferenças finitas, e modelos baseados em lattice para precificação de opções financeiras, exóticas e reais. Aponta as situações onde cada método deve ser empregado, e apresenta os algoritmos necessários para implementação dos métodos numéricos na prática.

Aplicação de sensoriamento remoto orbital no mapeamento de unidades vulcano-sedimentar no Platô da Ramada, Vila Nova do Sul, RS

Os produtos de sensoriamento remoto gerados pelos novos sensores orbitais, aliados ao desenvolvimento de sistemas computacionais e de técnicas de processamento de imagens possibilitam a manipulação numérica das imagens, visando o realce de informações de interesse para a interpretação visual ou automática. O uso dessas técnicas traz boas perspectivas no mapeamento geológico básico, pois possibilita uma visão sinóptica de amplas áreas e a integração de dados litológicos e estruturais de várias fontes. Este trabalho selecionou as imagens do sensor ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) devido as características de suas resoluções espectral e espacial, como as mais promissoras para identificação dos diferentes litotipos da área do Platô da Ramada, em especial as unidades vulcânicas e hipabissais relacionadas e sedimentares subjacentes. O Platô da Ramada, região de Vila Nova do Sul (RS) constitui uma feição geomorfológica de destaque na área, facilmente identificável nas imagens orbitais, sendo formado por uma seqüência vulcânica caracterizada por depósitos efusivos e piroclásticos, de composição dominantemente ácida. Representa uma fração significativa do magmatismo alcalino sódico neoproterozóico no Escudo Sul-rio-grandense e marca um dos ciclos vulcânicos do período pós-colisional do Ciclo Brasiliano Este estudo testou e avaliou diversos processamentos das imagens multiespectrais para diferenciação litológica, a identificação de alvos e a definição de morfoestruturas da área do Platô da Ramada. A integração de dados geológicos existentes com os produtos obtidos neste estudo possibilitou novas informações anexadas ao mapa geológico utilizado como verdade terrestre. O processamento utilizando a técnica de Transformação por Componentes Principais Seletivas proporcionou os melhores resultados realçando diferenças espectrais existentes entre as rochas vulcânicas e hipabissais e as rochas sedimentares. A partir dessa técnica, selecionou-se as imagens CP2 dos pares das bandas 4-5 (R), 3-4 (G) e 8-9 (B) na geração de uma composição colorida. Esta imagem permitiu a diferenciação espectral entre as rochas vulcânicas do Platô da Ramada e as rochas sedimentares do Grupo Maricá, bem como a individualização, no Grupo Maricá, de duas subunidades, levando-se em conta a densidade de intrusões de diques riolíticos ao norte da área. Identificou-se, ainda, um corpo diorítico com forma elíptica na borda SW do Platô da Ramada, não registrado anteriormente.

Solução semi-analítica da equação de Langevin assintótica para o deslocamento aleatório pelo método Picard

Neste trabalho é desenvolvida uma solução semi-analítica para a Equação de Langevin assintótica (Equação de Deslocamento Aleatório) aplicada à dispersão de poluentes na Camada Limite Convectiva (CLC). A solução tem como ponto de partida uma equação diferencial de primeira ordem para o deslocamento aleatório, sobre a qual é aplicado o Método Iterativo de Picard. O novo modelo é parametrizado por um coeficiente de difusão obtido a partir da Teoria de Difusão Estatística de Taylor e de um modelo para o espectro de turbulência, assumindo a supersposição linear dos efeitos de turbulência térmica e mecânica. A avaliação do modelo é realizada através da comparação com dados de concentração medidos durante o experimento de dispersão de Copenhagen e com resultados obtidos por outros quatro modelos: modelo de partículas estocástico para velocidade aleatória (Modelo de Langevin), solução analítica da equação difusão-advecção, solução numérica da equação difusão-advecção e modelo Gaussiano. Uma análise estatística revela que o modelo proposto simula satisfatoriamente os valores de concentração observados e apresenta boa concordância com os resultados dos outros modelos de dispersão. Além disso, a solução através do Método Iterativo de Picard pode apresentar algumas vantagem em relação ao método clássico de solução.

Uso do método das interfaces coesivas na análise da simulação do comportamento de fratura em materiais frágeis

Este trabalho apresenta uma abordagem numérica da fratura baseada no Método dos Elementos Finitos na qual são considerados materiais frágeis como objeto de estudo. As simulações realizadas são possíveis através de um algoritmo computacional que consiste na implementação de molas entre as faces dos elementos finitos. Esta composição constitui o método denominado Interfaces Coesivas. As molas fornecem o efeito de coesão e é representado por leis constitutivas do tipo bi-linear. Sua descrição é apresentada, bem como sua implementação na análise. Tanto o modo I como o modo II de ruptura são considerados. Com o objetivo de verificar o presente modelo com outros resultados e soluções, são analisados exemplos numéricos e experimentais extraídos da literatura. Tais exemplos tratam da nucleação, propagação e eventuais paradas de trincas em peças estruturais submetidas aos mais diversos carregamentos e condições de contorno. É ainda discutido como os parâmetros de fratura do material interferem no processo de ruptura. Uma tentativa é feita de determinar estes parâmetros, em especial no caso de concreto.

Estudo de recapeamentos asfálticos de pavimentos através de ensaios acelerados

O objetivo deste estudo é avaliar o desempenho de recapeamentos asfálticos, aplicados sobre pavimentos severamente trincados, através do acompanhamento da degradação ocorrida em função do tráfego. Foram testados dois tipos de recapeamentos: uma camada de 5,6 cm de espessura, em concreto asfáltico convencional e um tratamento superficial duplo, modificado por polímero SBS. Os recapeamentos asfálticos foram aplicados sobre uma pista experimental já trafegada e realizou-se, artificialmente, um trincamento padronizado buscando reproduzir efeitos similares àqueles provocados pelo fenômeno de fadiga. Utilizou-se a técnica de ensaios acelerados de pavimentos através da utilização de um simulador linear de tráfego, com rodado duplo. A carga de semi-eixo de 50 kN adotada proporcionou a degradação do pavimento três meses de ensaio para cada recapeamento. O monitoramento do pavimento, em termos estruturais e funcionais, fez-se pelo levantamento periódico de deflexões, bacias deflectométricas, afundamentos de trilha de roda e o controle do surgimento de trincas na superfície. Para melhor quantificar e compreender os fenômenos de degradação estrutural, instrumentou-se a pista experimental com sensores (células de tensão total e strain gages) na interface entre pavimento trincado e recapeamentos, obtendo-se tensões e deformações nestes locais, considerados críticos no desempenho de pavimentos recapeados. Durante os períodos de ensaio, foram monitoradas as condições pluviométricas e temperaturas do ar e pavimento. A determinação dos módulos elásticos das camadas das estruturas testadas foi efetuada pela aproximação das bacias obtidas em campo com os dados gerados em uma simulação numérica, através do programa computacional FLAPS. Percebeu-se que, além do tráfego, o trincamento no recapeamento em tratamento superficial foi devido à condição estrutural inicial do pavimento existente, severidade do trincamento artificial incompatível com a espessura desse tipo de restauração e baixas temperaturas. Quanto ao recapeamento em concreto asfáltico, percebeu-se que o maior trincamento surgiu em áreas sem trincas artificiais subjacentes, partindo do fundo da camada asfáltica do pavimento existente, em função da pouca rigidez da camada granular.

O custo de capital em operações de securitização de recebíveis de empresas não-financeiras por meio da emissão de quotas de Fundos de Investimento em Direitos Creditórios - FIDCs

Esta dissertação investiga, em caráter exploratório, se há evidência de que o custo de capital proporcionado a empresas brasileiras não-financeiras em operações de securitização por meio da estruturação de Fundos de Investimento em Direitos Creditórios – FIDCs (espécie de fundo de investimento regulamentada em 2001 pela Comissão de Valores Mobiliários – CVM) pode ser inferior ao custo de capital associado a operações de crédito tradicionais. A investigação baseia-se em comparação das notas de rating atribuídas a FIDCs estruturados para adquirir direitos creditórios de empresas não-financeiras e ao endividamento de longo prazo dessas mesmas empresas. A comparação entre notas de rating é estruturada a partir da premissa de que há forte correlação entre rating e spread (principal componente do custo de capital de terceiros), conforme documentado em diversos estudos empíricos precedentes. Foram selecionadas 23 operações de emissão de cotas de FIDCs, desde sua regulamentação até 2007, que atendem aos seguintes critérios: terem sido objeto de distribuição pública registrada na Comissão de Valores Mobiliários; terem mais de 5 cotistas seniores; e terem sido constituídos para adquirir recebíveis de uma única empresa não-financeira. A análise das 23 ocorrências fornece evidência de que a estruturação de FIDCs pode resultar na emissão de títulos securitizados com risco de crédito inferior àquele associado ao endividamento de longo prazo da empresa cedente. A redução do risco de crédito é observada por meio do diferencial de notas de rating em 8 das 23 operações analisadas. A ocorrência de 5 operações em que não há diferencial sugere, porém, que a redução do custo de capital não é a única razão pela qual os FIDCs vêm sendo estruturados. As limitações do rating como índice do custo de capital e as limitações numérica e temporal da amostra inviabilizam afirmações conclusivas acerca do impacto potencial do FIDC sobre o custo de capital. Não obstante, a análise sugere novas linhas de pesquisa que podem lançar mais luz sobre o assunto.

Dissipação de energia em estruturas com utilização de polímeros termoplásticos

Esta dissertação tem por objetivo contribuir à compreensão dos mecanismos e efeitos de dissipação de energia em estruturas mediante a utilização de materiais viscoelásticos. São obtidas experimentalmente as propriedades destes polímeros termoplásticos enfatizando seu uso como meio de redução das amplitudes de vibrações em estruturas submetidas a ações dinâmicas. Inicialmente, apresenta-se uma breve resenha histórica sobre o uso de materiais viscoelásticos com finalidade de amortecimento estrutural. São descritos modelos reológicos referentes ao comportamento mecânico, químico e térmico e também modelos matemáticos relacionados ao fenômeno do amortecimento necessários para o desenvolvimento do projeto dos amortecedores. A seguir, descreve-se detalhadamente o projeto, construção e análise de dissipadores de energia confeccionados com polímeros termoplásticos, mediante aplicações de carregamentos harmônicos cíclicos, é análisado seu comportamento. Finalmente, faz-se uma análise numérica da influência da incorporação de dissipadores numa estrutura real com problemas de vibrações excessivas. Conclusões sobre a aplicação deste tipo de amortecedor em estruturas em geral são apresentadas.

Aplicação de um modelo de dispersão atmosférico na região de Candiota-RS e sua relação com as condições meteorológicas

Este trabalho tem por objetivo, estudar a dispersão de poluentes na região da Usina Termelétrica Presidente Médici (UTPM), localizada no município de Candiota-RS, utilizando a modelagem numérica (mais especificamente o modelo ISCST - Industrial Surce Complex Term Short Term) e relacionar o seu comportamento com as condições meteorológicas. Esse modelo de dispersão utilizado é recomendado pela EPA (Environmental Protection Agency) para tratamento de dispersão de poluentes emitidos por fontes industriais como refinarias, termelétricas, siderúrgicas, etc. O período de estudo foi de 19/01/03 à 28/01/03, quando se realizou a primeira campanha experimental, pertencente ao projeto "Estudo Ambiental Aplicando a Modelagem Numérica na Região de Candiota-RS". Os resultados simulados pelo modelo foram comparados com os dados das concentrações de Material Particulado (MP) e Dióxido de Enxofre (SO 2 ), medidos através dos pontos receptores em duas estações de qualidade do ar próximas à Usina Termoelétrica. Uma correlação entre as concentrações de poluentes modelados e os dados coletados foi realizada, assim como, uma comparação com os padrões de qualidade do ar estipulados CONAMA (Conselho Nacional do Meio Ambiente). Os resultados denotam que o modelo ISCST subestima as concentrações de MP e SO 2 observadas nas estações de qualidade do ar, o que provavelmente ocorre, devido à presença de outras fontes de emissão (como indústrias de cimento e calcário) nas proximidades das estações. A situação com a melhor correlação nos dez dias, entre modelo e amostrador, apresentou-se com 90%. Em apenas um dos dez dias, foi excedido o padrão secundário (100 m g/m 3 ) com o poluente SO 2 , caracterizado a qualidade do ar da região como boa para esse período.

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.

Estudo crítico da metodologia de cálculo para silos prismáticos multicelulares

A armazenagem de grãos é atualmente um importante tópico a ser explorado na região sul do Brasil, e vem se desenvolvendo consideravelmente nos últimos anos. Atualmente a armazenagem é feita através de silos cilíndricos que podem ser de grandes diâmetros ou elevados. Na área da armazenagem é importante implantar novas tecnologias como o caso de silos multicelulares prismáticos que são bastante difundidos nos Estados Unidos. Para estudar estes silos propõe-se analisar um caso real sob a ótica morfológica e dimensional. O estudo baseia-se na comparação entre métodos de análise estrutural analítica convencional com análise numérica através do cálculo por elementos finitos, além de testar as análises de cálculo mais avançadas. Para estudar um caso real buscaram-se informações detalhadas através de desenhos fornecidos por empresas que fabricam estes tipos de silos. Adotando um caso prático como padrão, na seqüência fez-se a análise da interação entre o número de módulos do silo com o objetivo de verificar a influência do número de módulos no dimensionamento dos componentes do silo. Analisou-se o referido projeto utilizando procedimentos clássicos de análise de dimensionamento (métodos analíticos) e compararam-se os mesmos com a análise utilizando o método dos elementos finitos. Adicionalmente executaram-se análises mais avançadas utilizando o método dos elementos finitos como não linearidade geométrica e física, flambagem e análise de freqüências e modos.Após as referidas análises pôde-se concluir que analisando dimensionalmente somente um módulo do silo, pode-se obter resultados que representam todo o conjunto de módulos. Além disso, verificou-se que a análise por elementos finitos é mais apropriada, pois permite menos simplificações de cálculo, apresentando assim resultados mais reais. Quanto às análises adicionais, verificou-se que a carga crítica de flambagem considerando a análise não linear geométrica por elementos finitos é maior do que a carga crítica obtida por análises convencionais.

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.

Coordenação de projeto de obra de edificação : proposta de ferramenta computacional para programação e controle do fluxo de informações com uso de sistema colaborativo

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.

Alterações morfológicas na área pré-óptica medial, núcleo periventricular anteroventral e amígdala medial póstero-dorsal induzidas pela manipulação neonatal

A manipulação neonatal imprime alterações no desenvolvimento neuroendócrino, morfológico e comportamental de ratos. A área pré-óptica medial (MPOA) contém neurônios produtores do hormônio liberador do hormônio luteinizante (LHRH) que estão intimamente relacionados com a reprodução. Os núcleos periventricular anteroventral (AVPV) e a amígdala medial póstero-dorsal (MePD) são áreas do sistema nervoso central (SNC) que apresentam receptores de estrógeno, portanto poderiam estar atuando no feedback das gonadotrofinas. O objetivo deste trabalho foi verificar os efeitos da manipulação neonatal sobre o volume do núcleo, a densidade de células, o diâmetro do corpo celular de neurônios e o número total estimado de neurônios da MePD e AVPV dos lados direito e esquerdo, em ratas aos 11 e 90 dias de idade. Na MPOA, foram avaliadas a densidade de células e o diâmetro do corpo celular de neurônios. A MPOA apresentou uma redução da densidade numérica de 50% em ratas manipuladas de 11 dias e 50% em ratas manipuladas aos 90 dias de idade. A manipulação neonatal provocou alterações nos parâmetros analisados nas três áreas estudadas, quando comparadas aos grupos não manipulados. Na MePD a redução foi 63% e de 47% respectivamente aos 11 e aos 90 dias de idade. No AVPV a redução foi de 44% e de 54% aos 11 e 90 dias de idade. As dimensões lineares do soma dos neurônios da MePD e MPOA apresentaram uma redução de tal forma que o diâmetro dos neurônios no grupo manipulado foi significativamente menor do que o grupo não manipulado. Mas as dimensões lineares dos neurônios do AVPV não apresentaram alterações significativas. A manipulação neonatal por uma causa ainda desconhecida induz a alterações estáveis na MPOA, MePD e AVPV. A diminuição do número de neurônios nessas estruturas pode explicar a redução da atividade do eixo hipotálamo-hipófise-gonadal (HPG) em ratas adultas, manipuladas no período neonatal.

Simulações da dinâmica de fitoplâncton no sistema hidrológico do Taim (RS)

A variação temporal e espacial da ecologia de lagos, estuários e reservatórios está relacionada com a variação na composição da população de fitoplâncton, já que esta comunidade é a base da cadeia alimentar marinha e lacustre. Esta pesquisa desenvolveu um modelo para estimativa da dinâmica de biomassa de fitoplâncton, considerando os mecanismos de transporte no meio aquático, crescimento e perdas de sua biomassa. O modelo possui três módulos: (a) hidrodinâmico, que trata dos fluxos quantitativos do meio, associado a um algoritmo de secagem/inundação; (b) de transporte, que trata dos mecanismos de transporte das substâncias no meio; e (c) biológico, que retrata os mecanismos biológicos relacionados ao fitoplâncton. O módulo hidrodinâmico foi resolvido por um esquema semi-implícito de diferenças finitas, com uma abordagem mista Euleriana-Lagrangiana para os termos advectivos. Para os módulos restantes foi utilizado um esquema mais simples, de diferenças centrais. O modelo foi aplicado no sistema formado pela Lagoa Mangueira e Banhado do Taim localizado na parte sul costeira do estado do Rio Grande do Sul. Este sistema está inserido no Sistema Hidrológico do Taim como unidade de conservação federal, onde se desenvolve o projeto PELD (Pesquisas Ecológicas de Longa Duração). A hidrodinâmica do sistema foi ajustada e verificada apresentando bons resultados. Testes de secagem e inundação mostraram que as tendências esperadas foram obtidas. O aprimoramento da solução numérica do modelo foi comprovado através de testes de preservação de volume e massa em sistemas conservativos. A análise de sensibilidade dos parâmetros do módulo biológico mostrou que os parâmetros para os quais o modelo é mais sensível a alterações são aqueles relacionados aos efeitos de temperatura nas algas e as perdas por respiração e consumo por zooplâncton. Os cenários de análise propostos proporcionaram uma maior compreensão do transporte do material particulado suspenso e do material escalar genérico, assim como a dinâmica e as florações de algas no sistema em condições físicas e meteorológicas reais.

Sincronização idêntica em redes de osciladores caóticos

A presente dissertação versa sobre a sincronização idêntica em redes de osciladores caóticos. Uma perspectiva razoavelmente histórica sobre a literatura da área é apresentada . O conceito de caos é introduzido junto com outras idéias da dinâmica não-linear: sistemas dinâmicos, exemplos de sistemas, atratores, expoentes de Liapunov, etc. A integração numérica de equações diferenciais é largamente utilizada, principalmente, para o cálculo de expoentes e o desenho do diagrama de fases. A sincronização idêntica é definida, inicialmente, em redes que não passam de um par de osciladores. A variedade de sincronização (conjunto de pontos no espaço de fases no qual a solução do sistema é encontrada se há sincronização) é determinada. Diferentes variantes de acoplamentos lineares são enfocadas: acoplamento interno, externo, do tipo mestre-escravo e birecional, entre outras. Para detectar sincronização, usa-se o conceito de expoente de Liapunov transversal, uma extensão do conceito clássico de expoente de Liapunov que caracteriza a sincronização como a existência de um atrator na variedade de sincronização. A exposição é completada com exemplos e atinge relativo detalhe sobre o assunto, sem deixar de ser sintética com relação à ampla literatura existente. Um caso de sincronização em antifase que usa a mesma análise é incluído. A sincronização idêntica também é estudada em redes de osciladores idênticos com mais de dois osciladores. As possibilidades de sincronização completa e parcial são explanadas. As técnicas usadas para um par de osciladores são expandidas para cobrir este novo tipo de redes. A existência de variedades de sincronização invariantes é considerada como fator determinante para a sincronização. A sincronização parcial gera estruturas espaciais, analisadas sob a denominação de padrões. Algumas relações importantes entre as sincronizações são explicitadas, principalmente as degenerescências e a relação entre a sincronização parcial e a sincronização completa do respectivo estado sincronizado para alguns tipos de acoplamento. Ainda são objetos de interesse as redes formadas por grupos de osciladores idênticos que são diferentes dos osciladores dos outros grupos. A sincronização parcial na qual todos os grupos de osciladores têm seus elementos sincronizados é chamada de sincronização primária. A sincronização secundária é qualquer outro tipo de sincronização parcial. Ambas são exemplificadas e analisadas por meio dos expoentes transversais e novamente por meio da existência de invariantes de sincronização. Obtém-se, então, uma caracterização suficientemente ampla, completada por casos específicos.