Um método de matriz resposta com esquema iterativo de inversão parcial por região para problemas unidimensionais de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas

Um método de matriz resposta (RM) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas de transporte de nêutrons monoenergéticos e com fonte fixa, em geometria unidimensional na formulação de ordenadas discretas (SN). O método RM com esquema iterativo de inversão parcial por região (RBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações SN, afora os erros de arredondamento da aritmética finita computacional. Desenvolvemos um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em qualquer ponto do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Aplicações do formalismo de Keldysh ao transporte e ao bombeamento de calor em nanoestruturas

Num regime balísstico e a baixas temperaturas, a fórmula de Landauer dá uma boa descrição do transporte de calor para nano-junções conectadas a dois fios acoplados a banhos térmicos a temperaturas diferentes. Partindo de um modelo microscópico e utilizando o método de funções de Green fora do equilíbrio, é possível obter uma expressão para a condutância térmica na nano-junção equivalente a fórmula de Landauer. Esta depende dos valores das constantes de acoplamento entre os modos de fônons da região central e dos fios, além do gradiente térmico. A expressão para a condutância térmica é muito semelhante aquela obtida para a condutância elétrica. Neste trabalho nós apresentamos o método para o cálculo de grandezas relacionadas ao transporte térmico em um regime onde não há um gradiente de temperatura entre os reservatórios mas o sistema sofre uma perturbação dependente do tempo. Ou seja, com uma escolha conveniente da parametrização temporal dos termos de acoplamento entre a nano-junção e os fios é possível produzir uma corrente de calor na ausência de diferença de temperaturas entre os banhos térmicos aos quais os fios estão conectados. Esse fenômeno caracteriza o bombeamento de calor. Desenvolvemos uma teoria de transporte dependente do tempo para descrever o bombeamento. A teoria é geral, dependendo da densidade de fônons, da intensidade e dependência temporal do acoplamento. Aplicamos o formalismo em um modelo simples demonstrando que, em princípio, é possível bombear calor através de uma cadeia linear de átomos sem gradiente térmico.

Análise espectral das equações de transporte de nêutrons na formulação de ordenadas discretas em meios multiplicativos

É presentada nesta dissertação uma análise espectral das equações de transporte de nêutrons, independente do tempo, em geometria unidimensional e bidimensional, na formulação de ordenadas discretas (SN), utilizando o modelo de uma velocidade e multigrupo, considerando meios onde ocorrem o fenômeno da fissão nuclear. Esta análise espectral constitui-se na resolução de problemas de autovalores e respectivos autovetores, e reproduz a expressão para a solução geral analítica local das equações SN (para geometria unidimensional) ou das equações nodais integradas transversalmente (geometria retangular bidimensional) dentro de cada região homogeneizada do domínio espacial. Com a solução geral local determinada, métodos numéricos, tais como os métodos de matriz de resposta SN, podem ser derivados. Os resultados numéricos são gerados por programas de computadores implementados em MatLab, versão 2012, a fim de verificar a natureza dos autovalores e autovetores correspondentes no espaço real ou complexo.

Transporte eletrônico e quiralidade molecular: um estudo de dispositivos orgânicos em sistemas de dois terminais

No presente trabalho, investigamos o transporte eletrônico molecular em dois compostos orgânicos, o Ponceau SS (PSS) e o Oligo-(para)fenileno-vinileno (PPV) através de cálculos ab initio e função de Green de não equilíbrio (FGNE). Estes métodos demonstraram equivalência para a descrição destes dispositivos moleculares. Fizemos cálculos quânticos para o Hamiltoniano derivado de Hartree-Fock (HF) e obtivemos as propriedades de corrente-voltagem (I-V) para as duas estruturas moleculares. Com o método FGNE conseguimos modelar o transporte através de um sistema de multiníveis eletrônicos obtendo a corrente descrevendo as regiões de ressonância e a assimetria do sistema. Como resposta o PSS demonstrou assimetria para polarizações direta e reversa e a ressonância é alcançada mostrando que o dispositivo opere como um transistor molecular bi-direcional. Para o PPV investigamos também as propriedades geométricas através da conexão entre transporte eletrônico e o grau de quiralidade molecular que foi calculado usando o índice quiral que depende apenas das posições atômicas. Obtivemos que moléculas quirais e propriedades estruturais podem induzir uma assimetria no transporte eletrônico, resultando num processo de retificação. Também obtivemos que a resposta elétrica (I-V) e momento de dipolo elétrico são proporcionais ao grau de quiralidade molecular. Estes resultados sugerem que o transporte eletrônico neste sistema pode ser explorado na avaliação do seu grau de quiralidade.

Estudo da quantificação de incertezas para o problema de contaminação de meios porosos heterogêneos

As técnicas de injeção de traçadores têm sido amplamente utilizadas na investigação de escoamentos em meios porosos, principalmente em problemas envolvendo a simulação numérica de escoamentos miscíveis em reservatórios de petróleo e o transporte de contaminantes em aquíferos. Reservatórios subterrâneos são em geral heterogêneos e podem apresentar variações significativas das suas propriedades em várias escalas de comprimento. Estas variações espaciais são incorporadas às equações que governam o escoamento no interior do meio poroso por meio de campos aleatórios. Estes campos podem prover uma descrição das heterogeneidades da formação subterrânea nos casos onde o conhecimento geológico não fornece o detalhamento necessário para a predição determinística do escoamento através do meio poroso. Nesta tese é empregado um modelo lognormal para o campo de permeabilidades a fim de reproduzir-se a distribuição de permeabilidades do meio real, e a geração numérica destes campos aleatórios é feita pelo método da Soma Sucessiva de Campos Gaussianos Independentes (SSCGI). O objetivo principal deste trabalho é o estudo da quantificação de incertezas para o problema inverso do transporte de um traçador em um meio poroso heterogêneo empregando uma abordagem Bayesiana para a atualização dos campos de permeabilidades, baseada na medição dos valores da concentração espacial do traçador em tempos específicos. Um método do tipo Markov Chain Monte Carlo a dois estágios é utilizado na amostragem da distribuição de probabilidade a posteriori e a cadeia de Markov é construída a partir da reconstrução aleatória dos campos de permeabilidades. Na resolução do problema de pressão-velocidade que governa o escoamento empregase um método do tipo Elementos Finitos Mistos adequado para o cálculo acurado dos fluxos em campos de permeabilidades heterogêneos e uma abordagem Lagrangiana, o método Forward Integral Tracking (FIT), é utilizada na simulação numérica do problema do transporte do traçador. Resultados numéricos são obtidos e apresentados para um conjunto de realizações amostrais dos campos de permeabilidades.

Modelagem computacional da estabilização de sistemas subcríticos segundo o modelo unidimensional de difusão de nêutrons monoenergéticos

Um método numérico espectronodal (END) livre de erros de truncamento espacial é desenvolvido para problemas unidimensionais de difusão de nêutrons monoenergéticos em duas versões. Na versão de problemas de autovalor, o método gera soluções numéricas para o perfil do fluxo escalar e para o fator de multiplicação efetivo (k), que coincidem com a solução analítica dominante, afora os erros da aritmética finita computacional. Na versão de fonte fxa, o método também gera soluções numéricas analíticas para o problema de fonte fixa correspondente, onde a fonte de fissão, com dependência espacial, é obtida analiticamente, a partir da reconstrução espacial do fluxo escalar gerado pelo método END para problemas de autovalor. Alguns experimentos numéricos são apresentados para dois problemas modelos a fim de ilustrar a precisão do método.

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.

A review of particle transport theory in a binary stochastic medium

A finalidade deste trabalho ´e apresentar uma revis˜ao da teoria do transporte de part´ıculas em meios compostos por uma mistura aleat´oria bin´aria. Para atingir este objetivo n´os apresentamos brevemente alguns conceitos b´asicos de teoria do transporte, e ent˜ao discutimos em detalhes a deriva¸c˜ao de duas abordagens desenvolvidas para a solu¸c˜ao de tais problemas: os modelos de mistura atˆomica e de Levermore-Pomraning. Providenciamos ainda, com o uso da formula¸c˜ao LTSN, compara ¸c˜oes num´ericas destes modelos com resultados de benchmark gerados atrav´es de um processo de Monte Carlo.

Solução analítica do problema bidimensional transiente de dispersão de poluentes atmosféricos pelo método GITT dupla

Neste trabalho é apresentada uma solução analítica de um problema bidimensional e transiente de dispersão de poluentes atmosféricos. O modelamento utilizado é conhecido na literatura como modelo Kzz para dispersão de poluentes atmosféricos e é representado por uma equação difusivo-advectiva com coeficientes de difusão e advecção variáveis. São utilizados três diferentes coeficientes de difusão nas simulações, bem como as componentes horizontal e vertical do vento são tomadas como variáveis. A solução analítica é gerada através da aplicação da técnica GITT (Generalized Integral Transform Technique) dupla com problema transformado resolvido por Transformada de Laplace e diagonalização de matrizes. Filtros matemáticos são usados para homogenizar as condições de contorno viabilizando o uso da técnica citada. Além disso, o tipo de filtro matemático utilizado permite a sensível diminuição do custo computacional. Resultados numéricos são obtidos e comparados com dados experimentais e outras soluções da literatura.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Análise do processo de negociação de tarifas de transporte através da utilização da teoria dos jogos

O objetivo deste trabalho é apresentar uma visão diferenciada sobre um dos processos que compõe a fase de distribuição da cadeia de suprimentos: a negociação de tarifas de transportes. É possível analisar o desenvolvimento de um processo de negociação de tarifas de transporte expresso através da linguagem da teoria dos jogos. Na seqüência, um caso real de negociação e precificação de tarifas demonstra a aplicação prática do modelo utilizado. A aplicação da teoria dos jogos leva a resultados que demonstram os fatores que mais influenciam na negociação e formação dos preços das tarifas. O setor de logística tem crescido no Brasil e sua importância é cada vez maior dentro das organizações. Deste modo, prover ferramentas que tornem os agentes envolvidos nos processos de negociação mais racionais, com o uso da teoria econômica, pode ser um grande diferencial na agenda das organizações. As conclusões do trabalho alinham-se com as tendências observadas nas cadeias de suprimento, mais especificamente na logística de distribuição das organizações. A evolução do grau de preparo e profissionalismo dos operadores logísticos vai exigir que o processo de contratação e gestão de relacionamento com esses fornecedores seja mais sofisticada e efetivo por parte das empresas contratantes. Além disso, os processos tradicionais de compra e contratação de serviços logísticos podem não ser suficientes para maximizar os resultados das organizações através de redução de custos operacionais e elevação do nível de serviço ao cliente final (ambos diretamente relacionados à logistica de distribuição).

Desenvolvimento de um método espectronodal livre de erros de truncamento espacial para problemas adjuntos de transporte de partículas neutras monoenergéticas na formulação de ordenadas discretas em geometria unidimensional

Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto.