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

Estudo transversal sobre a exposição às tecnologias não-invasivas de cuidado de enfermagem obstétrica na atenção ao parto e sua repercussão na vitalidade do recém-nascido

Estudo cujo objeto foi o emprego de tecnologias não-invasivas de cuidado de enfermagem obstétrica (TNICEO) por enfermeiras obstétricas durante o acompanhamento do trabalho de parto e parto, e suas repercussões sobre a vitalidade do recém-nascido. A tese é: a disponibilização de TNICEO pelas enfermeiras obstétricas na atenção ao trabalho de parto, parto e nascimento está associada a recém-natos com índice de Apgar (IA) >8, quando comparado com o índice de Apgar de recém-nascidos cujas mães não puderam optar pelo uso destas tecnologias. Objetivou: a) descrever as características obstétricas das mulheres e de seu trabalho de parto e parto acompanhados por enfermeiras obstétricas; b) descrever as TNICEO disponibilizadas pelas enfermeiras obstétricas no cuidado à parturiente; c) medir e comparar a associação entre o IA de primeiro e quinto minutos de vida dos recém-nascidos cujas mães fizeram uso das TNICEO com: o IA do primeiro e quinto minutos de vida dos recém-nascidos cujas mães foram submetidas ao tratamento tradicional (intervenções); o IA do primeiro e quinto minutos de vida dos recém-nascidos cujas mães utilizaram tanto TNICEO e assistência tradicional (AT); o IA do primeiro e quinto minutos de vida dos recém-nascidos cujas mães e não foram submetidas a nenhum tipo de assistência (TNICEO ou AT) durante o trabalho de parto. Tratou-se de estudo observacional descritivo, transversal, retrospectivo. Realizado em Hospital Maternidade Municipal localizado na zona norte do Município do Rio de Janeiro. A amostra foi 6.790 parturientes, que tiveram parto vaginal acompanhado por enfermeiras obstétricas, entre setembro/2004 e dezembro/2011. A fonte de informações foi o Livro de Registros de Partos (LRP) da maternidade. Os resultados evidenciaram que: 91,9% utilizaram um ou mais recurso relacionado às TNICEO e 60,2% foram submetidas a uma ou mais intervenção da AT. Quanto ao IA no 1 e no 5 minuto de vida, com relação às variáveis relacionadas aos tipos de assistência que utilizaram, constatou-se que os neonatos cujas mães durante o trabalho de parto e parto utilizaram algum tipo de TNICEO apresentaram percentuais mais elevados de IA > 8, tanto no 1 minuto (93,4%) como no 5 minuto de vida (99,0%). Em contrapartida os neonatos cujas mães foram submetidas a algum procedimento relacionado à AT apresentaram os menores percentuais de IA (82,8%) no 1 minuto e (94,7%) no 5 minuto de vida. Confirmando a tese proposta, conclui-se que a razão de chance do IA ser > 8 no primeiro minuto de vida é aumentada (OR 4,564; IC95%: 1,887 11,038; p valor 0,0008) a favor do grupo de mulheres que utilizaram apenas as TNICEO durante TP e/ou P, comparado ao grupo de mulheres submetidas às intervenções da AT. Referente ao quinto minuto de vida, a razão de chance de um recém-nascido cuja mãe que teve seu TP e/ou P acompanhados pelas enfermeiras obstétricas e utilizaram as somente as TNICEO na assistência foi maior (OR = 4,927; IC95% 2,349 10,334 ; p valor 0,00020), quando comparado com o grupo de parturientes da AT.

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

Inversão de funções do plano no plano aplicada ao cálculo de azeótropos

Azeotropia é um fenômeno termodinâmico onde um líquido em ebulição produz um vapor com composição idêntica. Esta situação é um desafio para a Engenharia de Separação, já que os processos de destilação exploram as diferenças entre as volatilidades relativas e, portanto, um azeótropo pode ser uma barreira para a separação. Em misturas binárias, o cálculo da azeotropia é caracterizado por um sistema não-linear do tipo 2 × 2. Um interessante e raro caso é o denominado azeotropia dupla, que pode ser verificado quando este sistema não-linear tem duas soluções, correspondendo a dois azeótropos distintos. Diferentes métodos tem sido utilizados na resolução de problemas desta natureza, como métodos estocásticos de otimização e as técnicas intervalares (do tipo Newton intervalar/bisseção generalizada). Nesta tese apresentamos a formulação do problema de azeotropia dupla e uma nova e robusta abordagem para a resolução dos sistemas não-lineares do tipo 2 × 2, que é a inversão de funções do plano no plano (MALTA; SALDANHA; TOMEI, 1996). No método proposto, as soluções são obtidas através de um conjunto de ações: obtenção de curvas críticas e de pré-imagens de pontos arbritários, inversão da função e por fim, as soluções esperadas para o problema de azeotropia. Esta metodologia foi desenvolvida para resolver sistemas não-lineares do tipo 2 × 2, tendo como objetivo dar uma visão global da função que modela o fenômeno em questão, além, é claro, de gerar as soluções esperadas. Serão apresentados resultados numéricos para o cálculo dos azeótropos no sistema benzeno + hexafluorobenzeno a baixas pressões por este método de inversão. Como ferramentas auxiliares, serão também apresentados aspectos numéricos usando aproximações clássicas, tais como métodos de Newton com técnicas de globalização e o algorítmo de otimização não-linear C-GRASP, para efeito de comparação.

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.

O demônio familiar: pensamento social e psiquismo em Os sertões

Esta tese se dedica a fazer uma análise da obra de Euclides da Cunha, publicada em 1902, Os sertões. Para tanto, toma como campos de investigação a história do pensamento social brasileiro, a história dos discursos e os saberes da psicologia de fins do século XIX, com o objetivo de examinar a semântica do conceito de psique, fundamental na construção do argumento do livro. Embora não caracterizado por sentido específico, o sistema semântico da psique pode ser compreendido no livro de Euclides através de uma regra de semelhança que se configura como regra de reflexão, isto é, tradução entre o organismo biológico e o organismo social. Na linguagem de Os sertões, como argumentamos, o psíquico se realiza como metáfora que sustenta a tradução e a regra de semelhança entre fisiologia e sociologia. Neste sentido, em reflexão com os ensaios de Hans Blumenberg e com a teoria dos sistemas de Niklas Luhmann, analisamos o livro de Euclides de modo a compreender o horizonte semântico ao redor do qual o seu psiquismo aparece sistematizado. Em Canudos, quando a ordem social republicana foi atacada, Euclides depura o argumento de que os crimes da nacionalidade derivam os seus motivos da inconsciência generalizada dos habitantes do litoral sobre as populações rurais caracterizando Canudos como um crime de consciência de onde o autor reclama para si a tarefa de vingar, isto é, tornar conhecidas as populações historicamente ignoradas, socialmente esquecidas pela civilização. A semântica da psique, nesse sentido, assume para Euclides uma técnica de observação, mas também uma hipótese política sobre as condições de sobrevivência da sociedade brasileira.

O que há de novo nas toxicomanias?

O presente trabalho de pesquisa envolve os casos de toxicomanias, como um novo sintoma, e a relação com a subjetividade de nossa época, provocando um questionamento sobre a clínica e a direção do tratamento. A mudança operada na cultura destaca o discurso do mestre atual como aquele que leva em si a marca do declínio dos ideais e das identificações não mais estabelecidos por um significante mestre e, sim, pelo mercado, provocando um imperativo feroz de Um gozo que promete a satisfação plena através do objeto, no caso, a droga. Vimos nos casos das toxicomanias e do alcoolismo, a ilusória completude capturada por um gozo auto- erótico que ludibria o sintoma com a droga e silencia o sujeito. O resultado do entrelaçamento do discurso capitalista com o discurso da ciência recai sobre a droga inserindo-a como mercadoria, cuja produção cada vez mais desenfreada, empurra a uma intoxicação generalizada, de difícil desintoxicação. Com a orientação psicanalítica é possível estabelecer um dispositivo que promova essa desintoxicação, reinventado caso a caso. É ofertar através do discurso do analista, um vazio de significações, de uma não resposta, em que, através da fala do sujeito, seja possível pinçar a função da droga na economia psíquica e os significantes que marcam sua história. O debate sobre a política pública e a função da psicanálise nos serviços públicos deve-se orientar pela singularidade e a dimensão ética do gozo no que concerne a proposta universalizante da política pública.

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.

Análise dos problemas direto e inverso em um modelo de evolução populacional através de transformações integrais e inferência bayesiana

Modelos de evolução populacional são há muito tempo assunto de grande relevância, principalmente quando a população de estudo é composta por vetores de doenças. Tal importância se deve ao fato de existirem milhares de doenças que são propagadas por espécies específicas e conhecer como tais populações se comportam é vital quando pretende-se criar políticas públicas para controlar a sua proliferação. Este trabalho descreve um problema de evolução populacional difusivo com armadilhas locais e tempo de reprodução atrasado, o problema direto descreve a densidade de uma população uma vez conhecidos os parâmetros do modelo onde sua solução é obtida por meio da técnica de transformada integral generalizada, uma técnica numérico-analítica. Porém a solução do problema direto, por si só, não permite a simulação computacional de uma população em uma aplicação prática, uma vez que os parâmetros do modelo variam de população para população e precisam, portanto, ter seus valores conhecidos. Com o objetivo de possibilitar esta caracterização, o presente trabalho propõe a formulação e solução do problema inverso, estimando os parâmetros do modelo a partir de dados da população utilizando para tal tarefa dois métodos Bayesianos.

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.