Dinâmica populacional e caracterização conquiliomorfológica de Diplodon ellipticus Spix in Wagner, 1827 (Bivalvia: Unionoida: Hyriidae) da lagoa dos Caiçaras (Piraí, Rio de Janeiro, Brasil)

No presente estudo monitoramos uma população de Diplodon ellipticus Spix in Wagner, 1827 da lagoa dos Caiçaras, Piraí, Rio de Janeiro, e avaliamos alguns aspectos conquiliomorfológicos da larva e do adulto, a relação peso-comprimento, os padrões populacionais tais como crescimento, mortalidade e expectativa de vida, e o ciclo larval. Amostragens mensais foram realizadas de novembro/2012 a novembro/2013 em três pontos da lagoa. Em cada ponto, definimos uma area de 15 m2, que foi subdividida em 15 quadrados de 1m2 cada. Os bivalves foram procurados por três coletores, usando mãos e pés, totalizando 45 minutos de coleta/área. Posteriormente, foram medidos com um paquímetro em relação ao comprimento total, marcados e devolvidos a lagoa. Quinze bivalves foram coletados e analisados em laboratório durante dois anos. Os fatores abióticos (condutividade, temperatura da água, temperatura do ambiente, umidade, pH e oxigênio dissolvido) foram mensurados e amostragens do sedimento e da água foram realizados em cada ponto. Os indivíduos foram agrupados em classes de comprimento com intervalo de 2,0 mm para análise da estrutura de comprimento da população. Os parâmetros de curva de crescimento foram estimados pela rotina ELEFAN (distribuição de frequência) e pelo método de Gulland-Holt (marcação e recaptura), ambos no programa FISAT. A mortalidade foi calculada pelo método da curva de captura convertida em comprimento e a expectativa de vida foi calculada pela equação invertida de von Bertalanffy. Ao todo, 3474 bivalves foram marcados e 1849 recapturados, alguns deles mais de uma vez. O menor bivalve recapturado mediu 23,58 mm e foi recapturado 127 dias depois com 25,50 mm e o maior bivalve recapturado mediu 62,14 mm e foi recapturado 31 dias depois com 62,20 mm. Pelo método indireto analizamos 6922 bivalves com comprimentos variando de 11,37 a 62,49 mm. A maior frequência de comprimento foi encontrada em tamanhos intermediários de 44-46 mm, como foi observado em outras populações de Diplodon. A análise do crescimento, mortalidade e expectativa de vida foram similares em ambos os métodos, por isso, utilizamos em conjunto para ajustar a curva de crescimento da população. A relação peso-comprimento foi alta (r = 0,7-0,8). A avaliação do ciclo reprodutivo indica uma continuidade no desenvolvimento larval ao longo dos meses, à exceção de janeiro/2014. Contabilizamos um total de 54.617 gloquídios, havendo mais indivíduos grávidos com gloquídios (n=5) nos meses de julho e setembro de 2014. Padrão similar foi encontrado em outras espécies de bivalves de água doce, inclusive em Diplodon ellipticus no estado do Paraná. Os gloquídios maduros foram encontrados em amplitudes de temperaturas de 20,1-25,27C e acreditamos que a temperatura pode influenciar na liberação dos mesmos. Os bivalves com o maior e o menor número de gloquídios foram encontrados em novembro/2013 (n= 4759) e setembro/2014 (n= 212), respectivamente. A correlação entre o comprimento do marsúpio e o tamanho da ninhada foi testada pela primeira vez nesse gênero e foi significativa (r= 0,303 , p<0,05). As informações deste estudo são relevantes ao conhecimento da espécie e às futuras ações de gestão que tenham como objetivo preservar as espécies de Diplodon

Determinação experimental e modelagem termodinâmica do equilíbrio de fases em sistemas com CO2, n-hexano, n-hexadecano e tetralina

A descoberta de petróleo na camada de Pré-Sal possibilita a geração de ganhos em relação à dependência energética do país, mas também grandes desafios econômicos e tecnológicos. Os custos de extração são maiores devido a vários fatores como a exigência de equipamentos de exploração que suportem elevadas pressões, altas temperaturas e grandes concentrações de gases ácidos, tais quais, dióxido de carbono (CO2) e sulfeto de hidrogênio (H2S). Uma das principais preocupações com o CO2 é evitar liberá-lo para a atmosfera durante a produção. Com a modelagem termodinâmica de dados de equilíbrio de sistemas envolvendo CO2 supercrítico e hidrocarbonetos é possível projetar equipamentos utilizados em processos de separação. A principal motivação do trabalho é o levantamento de dados de equilíbrio de fases de sistemas compostos de CO2 e hidrocarbonetos, possibilitando assim prever o comportamento dessas misturas. Os objetivos específicos são a avaliação do procedimento experimental, a estimação e predição dos parâmetros de interação binários para assim prever o comportamento de fases dos sistemas ternários envolvendo CO2 e hidrocarbonetos. Duas metodologias foram utilizadas para obtenção dos dados de equilíbrio: método estático sintético (visual) e método dinâmico analítico (recirculação das fases). Os sistemas avaliados foram: CO2 + n-hexano, CO2 + tetralina, CO2 + n-hexadecano, CO2 + n-hexano + tetralina e CO2 + tetralina + n-hexadecano à alta pressão; tetralina + n-hexadecano à baixa pressão. Para o tratamento dos dados foi utilizada equação de estado cúbica de Peng-Robinson e a regra de mistura clássica

Modelagem matemática determinística do fluxo do tráfego veicular

Nesta dissertação consideramos duas abordagens para o tráfego de veículos: a macroscópica e a microscópica. O tráfego é descrito macroscopicamente por três grandezas físicas interligadas entre si, a saber, a velocidade, a densidade e o fluxo, descrevendo leis de conservação do número de veículos. Há vários modelos para o tráfego macroscópico de veículos. A maioria deles trata o tráfego de veículos como um fluido compressível, traduzindo a lei de conservação de massa para os veículos e requer uma lei de estado para o par velocidade-densidade, estabelecendo uma relação entre eles. Já o modelo descrito pela abordagem microscópica considera os veículos como partículas individuais. Consideramos os modelos da classe "car - following". Estes modelos baseiam-se no princípio de que o (n - 1)-ésimo veículo (denominado de "following-car") acelera em função do estímulo que recebe do n-ésimo veículo. Analisamos a equação de conservação do número de veículos em modelos macroscópicos para fluxo de tráfego. Posteriormente resolvemos esta equação através da linearização do modelo, estudando suas retas características e apresentamos a resolução do problema não linear em domínios limitados utilizando o método das características

Emissividade de neutrinos pelo processo Urca direto na presença de um campo magnético

As estrelas de nêutrons nascem com altas temperaturas (~ 1011 K) e durante alguns segundos sofrem um rápido resfriamento por emissão de neutrinos. O processo Urca direto é o principal mecanismo para explicar essa perda de energia. O problema do resfriamento das estrelas de nêutrons é um problema de grande interesse porque seu entendimento pode fornecer informações importantes sobre a constituição do interior da estrela. Na literatura existente até o momento, a emissividade de neutrinos é calculada considerando os núcleons como partículas não relativísticas quando considerados todos os níveis de Landau das partículas carregadas. Por outro lado, a emissividade de neutrinos para núcleons relativísticos é calculada considerando somente o primeiro nível de Landau (para campo magnético forte). Para campos magnéticos fracos, onde mais de um nível de Landau é ocupado, é usada a emissividade correspondente à do campo nulo. Neste trabalho aplicamos a teoria de Weinberg-Salan para interações fracas no cálculo da emissividade de neutrinos com e sem campo magnético presente, num cálculo totalmente relativístico para os núcleons e considerando todos os níveis de Landau. Esta é a contribuição original do trabalho. Para descrever a matéria a altas densidades, utilizamos uma teoria relativística de campo médio a temperatura zero que inclui apenas o octeto bariônico e os léptons mais leves. São apresentados os resultados para a emissividade de neutrinos, onde é evidente a ocupação dos diferentes níveis de Landau como função do campo magnético.

Uma Simulação computacional do passeio aleatório simples

Em 1828 foi observado um fenômeno no microscópio em que se visualizava minúsculos grãos de pólen mergulhados em um líquido em repouso que mexiam-se de forma aleatória, desenhando um movimento desordenado. A questão era compreender este movimento. Após cerca de 80 anos, Einstein (1905) desenvolveu uma formulação matemática para explicar este fenômeno, tratado por movimento Browniano, teoria cada vez mais desenvolvida em muitas das áreas do conhecimento, inclusive recentemente em modelagem computacional. Objetiva-se pontuar os pressupostos básicos inerentes ao passeio aleatório simples considerando experimentos com e sem problema de valor de contorno para melhor compreensão ao no uso de algoritmos aplicados a problemas computacionais. Foram explicitadas as ferramentas necessárias para aplicação de modelos de simulação do passeio aleatório simples nas três primeiras dimensões do espaço. O interesse foi direcionado tanto para o passeio aleatório simples como para possíveis aplicações para o problema da ruína do jogador e a disseminação de vírus em rede de computadores. Foram desenvolvidos algoritmos do passeio aleatório simples unidimensional sem e com o problema do valor de contorno na plataforma R. Similarmente, implementados para os espaços bidimensionais e tridimensionais,possibilitando futuras aplicações para o problema da disseminação de vírus em rede de computadores e como motivação ao estudo da Equação do Calor, embora necessita um maior embasamento em conceitos da Física e Probabilidade para dar continuidade a tal aplicação.

Avaliação de modelos de permeabilidade em meios porosos não consolidados

As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

Um estudo de caso envolvendo a retroanálise de tensões residuais de cravação.

As tensões residuais de cravação ocorrem quando a estaca e o solo não são totalmente descomprimidos após dissipação completa da energia transferida à estaca durante sua instalação. Este assunto tem merecido a atenção de vários pesquisadores, uma vez que a presença das tensões residuais na estaca causa uma alteração no seu comportamento quando carregada. A presente dissertação reuniu estudos anteriores que estabeleceram os principais fatores que influenciam as tensões residuais de cravação. Procurou-se executar uma série de simulações objetivando melhor entender o desenvolvimento das tensões residuais com a porcentagem de resistência de ponta, em continuidade aos estudos de Costa (1994). Costa (1994) observou que a porcentagem de carga na ponta é um fator relevante na avaliação das tensões residuais de cravação. Uma análise paramétrica efetuada confirmou estudos anteriores, verificando que a razão entre a carga residual na ponta em relação à capacidade de carga global cresce, a medida que a porcentagem de carga na ponta aumenta, chegando a um valor máximo para, em seguida, diminuir. Este comportamento é similar ao da curva de compactação do solo, quando, ao aumentar a umidade, o peso específico seco aumenta, até um valor máximo, correspondente à umidade ótima, para em seguida reduzir. Verificouse, ainda, semelhantemente ao aumento da energia de compactação, que vai transladando a curva para cima e para à esquerda do gráfico umidade x peso específico seco, que o aumento do comprimento da estaca apresenta um comportamento similar. O aumento do comprimento leva a curva para cima e para a esquerda. Finalmente, selecionou-se um caso de obra com condições propícias para o desenvolvimento de altas tensões residuais de cravação. Trata-se de um caso de estacas metálicas longas, embutidas em solo residual jovem, de elevada resistência, bem documentado com extensa instrumentação, por ocasião da instalação. A retroanálise de cinco estacas do banco de dados mostrou que as cargas residuais previstas, num programa de simulação de cravação, se aproximaram muito dos valores experimentais. Foi observado também que a profundidade do ponto neutro previsto e medido apresentou uma excelente concordância. O resultado mais relevante desta pesquisa foi quando os valores de porcentagem de ponta foram introduzidos no eixo das abscissas e os valores da razão entre a carga residual na ponta e a capacidade de carga global, no eixo das ordenadas, e se observou o aspecto da curva semelhante à de compactação. A dissertação ilustra ser possível e simples a previsão das tensões residuais de cravação através de uma análise pela equação da onda.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Estimativa de fontes externas uniformes de nêutrons em regiões combustíveis para uma distribuição prescrita de potência

Projetos de reatores nucleares foram classificados em quatro gerações (Gen) pelo Departamento de Energia dos Estados Unidos da América (DOE), quando o DOE introduziu o conceito de reatores de geração IV (Gen IV). Reatores Gen IV são um conjunto de projetos de reator nuclear, em sua maioria teóricos, atualmente sendo pesquisados. Entre os projetos Gen IV, incluem-se os projetos dos ADS (Accelerator Driven Systems), que são sistemas subcríticos estabilizados por fontes externas estacionárias de nêutrons. Estas fontes externas de nêutrons são normalmente geradas a partir da colisão de prótons com alta energia contra os núcleos de metais pesados presentes no núcleo do reator, fenômeno que é conhecido na literatura como spallation, e os prótons são acelerados num acelerador de partículas que é alimentado com parte da energia gerada pelo reator. A criticalidade de um sistema mantido por reações de fissão em cadeia depende do balanço entre a produção de nêutrons por fissão e a remoção por fuga pelos contornos e absorção de nêutrons. Um sistema está subcrítico quando a remoção por fuga e absorção ultrapassa a produção por fissão e, portanto, tende ao desligamento. Entretanto, qualquer sistema subcrítico pode ser estabilizado pela inclusão de fontes estacionárias de nêutrons em seu interior. O objetivo central deste trabalho é determinar as intensidades dessas fontes uniformes e isotrópicas de nêutrons, que se deve inserir em todas as regiões combustíveis do sistema, para que o mesmo estabilize-se gerando uma distribuição prescrita de potência elétrica. Diante do exposto, foi desenvolvido neste trabalho um aplicativo computacional em linguagem Java que estima as intensidades dessas fontes estacionárias de nêutrons, que devem ser inseridas em cada região combustível para que estabilizem o sistema subcrítico com uma dada distribuição de potência definida pelo usuário. Para atingir este objetivo, o modelo matemático adotado foi a equação unidimensional de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas (SN) e o convencional método de malha fina diamond difference (DD) foi utilizado para resolver numericamente os problemas SN físicos e adjuntos. Resultados numéricos para dois problemas-modelos típicos são apresentados para ilustrar a acurácia e eficiência da metodologia proposta.

Simulação numérica da dispersão advecção de pesticidas no solo sob efeito da temperatura.

É apresentado um modelo de dispersão-advecção de evolução unidimensional que simula a lixiviação de pesticidas em lisímetros ou colunas de solo sob efeito da temperatura média diária do perfil do solo. O modelo matemático e todo o seu conjunto de hipóteses será denominado DAPESTE. Nas simulações numéricas do modelo DAPESTE serão utilizados o método dos elementos finitos para a semi-discretização da variável espacial e o método de Eüler regressivo para a discretização da variável temporal. Serão utilizados métodos de elementos finitos apropriados para problemas de dispersão-advecção nos quais o transporte advectivo predomina sobre o transporte dispersivo. O modelo DAPESTE supõe que as difusividades do pesticida nas fases gasosa e aquosa do solo dependem da temperatura média diária do solo a qual varia periodicamente com a profundidade e com o tempo. O coeficiente de dispersão hidrodinâmico do modelo DAPESTE dependerá da temperatura do solo. O coeficiente de partição água-ar, variando com a temperatura, será determinado pela equação de Clausius-Clapeyron. A equação de van?t Ho® será usada para determinar o coeficiente de sorção do pesticida no solo como função da temperatura. Com a equação de Arrhenius será estimado o efeito da temperatura na taxa de degradação do pesticida. Estas relações de dependência entre os parâmetros do modelo e a temperatura do solo auxiliam sobremaneira na compreensão do destino de pesticidas no solo sob diferentes cenários de temperaturas médias diárias, especialmente a meia vida do pesticida e a concentração lixiviada no perfil do solo.

Novo método para estimativa da evapotranspiração de referência com base na temperatura do ar.

A região de Jales, localizada no noroeste de São Paulo, é uma importante área produtora de uvas de mesa do Estado, onde predominam pequenas propriedades com mão-de-obra familiar. Em todos os vinhedos da região são empregados sistemas de irrigação, para suprir a deficiência hídrica que ocorre durante o ciclo produtivo da cultura. Para um manejo adequado da irrigação, são necessárias informações sobre a necessidade hídrica da cultura. A evapotranspiração de referência (ETo) é um parâmetro fundamental para a determinação dessa necessidade, sendo desejável que se tenha um método que estime a ETo com boa precisão e a partir de dados meteorológicos fáceis de serem obtidos (VILLA NOVA; PEREIRA, 2006). O método de Penman-Monteith é considerado, atualmente, como padrão para a estimativa de ETo. Para seu uso, entretanto, são necessárias variáveis meteorológicas nem sempre disponíveis, principalmente aos pequenos produtores. Por essa razão, métodos que empregam um menor número de variáveis são também utilizados para estimar a ETo. No presente trabalho foram calculados valores diários de ETo pelo método de Hargreaves original; e por um novo método, em que se emprega a equação de Bristow-Campbell para a estimativa de Rs. Esses valores foram comparados aos obtidos pelo método padrão de Penman-Monteith, nas condições do noroeste paulista.

Vulnerabilidade dos solos à erosão: estimativa da perda de solos na bacia hidrográfica do Rio Guapi-Macacu - RJ.

Este trabalho apresenta metodologias para estimativa de perda de solos na Bacia Hidrográfica do Rio Guapi-Macacu - RJ. Para a estimativa anual de perdas de solo foi utilizada a ferramenta InVest (Integrated Valuation of Environmental Services and Tradeoffs) através do módulo que desenvolve retenção de sedimentos. Esse módulo permite calcular a perda de solo média anual de cada parcela de terra, além de determinar o quanto de solo pode chegar a um determinado ponto de interesse. Para identificar a perda potencial de solo, o modelo emprega a Equação Universal de Perda de Solo (USLE) na escala do pixel, integrando informações sobre relevo, precipitação, padrões de uso da terra e propriedades do solo. Seus resultados são dados em toneladas por sub-bacias. Os resultados obtidos neste trabalho permitiram concluir que embora haja limitações no uso da Equação Universal de Perda de Solo, o modelo possibilitou a espacialização de classes de perdas de solos com indicações de áreas consideradas mais ou menos vulneráveis aos processos erosivos, considerando os dados disponíveis e suas escalas. O uso do InVest para calcular a USLE apresentou como principal vantagem a possibilidade de integração dos dados necessários em um único ambiente, reduzindo a possibilidade de erros na conversão de dados. Contudo, a maior limitação encontrada a sua aplicação está na dificuldade de obtenção de dados de entrada necessários ao modelo.

Estudo do processo de difusão de gases em carvões com base em isotérmicas de Langmuir

Tese de Doutoramento apresentada à Universidade Fernando Pessoa como parte dos requisitos para obtenção do grau de Doutror em Ciências da Terra.

Estimação robusta em modelos lineares de equações simultâneas

Os Modelos de Equações Simultâneas (SEM) são modelos estatísticos com muita tradição em estudos de Econometria, uma vez que permitem representar e estudar uma vasta gama de processos económicos. Os estimadores mais usados em SEM resultam da aplicação do Método dos Mínimos Quadrados ou do Método da Máxima Verosimilhança, os quais não são robustos. Em Maronna e Yohai (1997), os autores propõem formas de “robustificar” esses estimadores. Um outro método de estimação com interesse nestes modelos é o Método dos Momentos Generalizado (GMM), o qual também conduz a estimadores não robustos. Estimadores que sofrem de falta de robustez são muito inconvenientes uma vez que podem conduzir a resultados enganadores quando são violadas as hipóteses subjacentes ao modelo assumido. Os estimadores robustos são de grande valor, em particular quando os modelos em estudo são complexos, como é o caso dos SEM. O principal objectivo desta investigação foi o de procurar tais estimadores tendo-se construído um estimador robusto a que se deu o nome de GMMOGK. Trata-se de uma versão robusta do estimador GMM. Para avaliar o desempenho do novo estimador foi feito um adequado estudo de simulação e foi também feita a aplicação do estimador a um conjunto de dados reais. O estimador robusto tem um bom desempenho nos modelos heterocedásticos considerados e, nessas condições, comporta-se melhor do que os estimadores não robustos usados no estudo. Contudo, quando a análise é feita em cada equação separadamente, a especificidade de cada equação individual e a estrutura de dependência do sistema são dois aspectos que influenciam o desempenho do estimador, tal como acontece com os estimadores usuais. Para enquadrar a investigação, o texto inclui uma revisão de aspectos essenciais dos SEM, o seu papel em Econometria, os principais métodos de estimação, com particular ênfase no GMM, e uma curta introdução à estimação robusta.