Solução analítica da equação de difusão-advenção pelo método GILTT aplicada à dispersão de poluentes atmosféricos

O objetivo deste trabalho é obter uma nova solução analítica para a equação de advecção-difusão. Para tanto, considera-se um problema bidimensional difusivo-advectivo estacionário com coeficiente de difusão turbulenta vertical variável que modela a dispersão de poluentes na atmosfera. São utilizados três coeficientes difusivos válidos na camada limite convectiva e que dependem da altura, da distância da fonte e do perfil de velocidade. A abordagem utilizada para a resolução do problema é a técnica da Transformada Integral Generalizada, na qual a equação transformada do problema difusivo-advectivo é resolvida pela técnica da Transformada de Laplace com inversão analítica. Nenhuma aproximação é feita durante a derivação da solução, sendo assim, esta é exata exceto pelo erro de truncamento. O modelo ´e avaliado em condições moderadamente instáveis usando o experimento de Copenhagen. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com dados experimentais e com os resultados da literatura. O modelo proposto mostrou-se satisfatório em relação aos dados dos experimentos difusivos considerados.

A equação de Euler e a IS: estimações para o Brasil

Neste artigo estudamos a relação entre a taxa de juros e o hiato do produto no Brasil através da estimação de modelos Novo-Keynesianos. Para tanto, estimamos os modelos por três métodos: (1) método generalizados dos momentos, (2) máxima verossimilhança utilizando dados de expectativa divulgados pelo Banco Central do Brasil e (3) máxima verossimilhança utilizando variáveis de expectativa estimadas por um modelo VAR. As conclusões são altamente dependentes do método de estimação. Ao utilizar (1), os resultados indicam uma relação espúria entre a taxa de juros real e o hiato. Entretanto, as conclusões originadas de (2) indicam que somente o hiato defasado seria uma variável relevante para a nossa especificação. Ao estimar o modelo por (3), as estimativas corroboram os resultados obtidos com o método generalizados dos momentos.

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

Solução LTSn da equação adjunta de transporte de nêutrons com fonte arbitrária para elevada ordem de quadratura numa placa homogênea

O objetivo deste trabalho consiste em estender o método LTSN à solução do problema adjunto de transporte de nêutrons. A solução adjunta é interpretada fisicamente como uma função importância que designa a capacidade de contribuição de cada cela do espaço de fase para um funcional resposta. A derivação desta interpretação, através do princípio variacional, está sucintamente apresentada. Surgida da necessidade de generalização da fonte adjunta, também propõe-se uma nova formulação LTSN capaz de resolver problemas de transporte, tanto direto quanto adjunto, com fonte arbitrária, para elevada ordem de quadratura em geometria de placa. Esta nova formulção inspira-se na propriedade de invariância de projeção dos meios isotrópicos mas também é válida para os meios anisotrópicos. Todos os resultados apresentados pelas simulações numéricas de problemas adjuntos são calculados pela nova formulação LTSN e são comparados ou com a definição de função importância ou pelas relações de reciprocidade ou pelo código ANISN.

Solução da equação advectivo-difusiva utilizando regras para manipulação de exponenciais de operadores diferenciais e simetrias de Lie : aplicações em engenharia ambiental

Vários métodos analíticos, numéricos e híbridos podem ser utilizados na solução de problemas de difusão e difusão-advecção. O objetivo deste trabalho é apresentar dois métodos analíticos para obtenção de soluções em forma fechada da equação advectivo-difusiva em coordenadas cartesianas que descreve problemas de dispersão de poluentes na água e na atmosfera. Um deles é baseado em regras de manipulação de exponenciais de operadores diferenciais, e o outro consiste na aplicação de simetrias de Lie admitidas por uma equação diferencial parcial linear. Desenvolvem-se regras para manipulação de exponenciais de operadores diferenciais de segunda ordem com coeficientes constantes e para operadores advectivo-difusivos. Nos casos em que essas regras não podem ser aplicadas utiliza-se uma formulação para a obtenção de simetrias de Lie, admitidas por uma equação diferencial, via mapeamento. Define-se um operador diferencial com a propriedade de transformar soluções analíticas de uma dada equação diferencial em novas soluções analíticas da mesma equação. Nas aplicações referentes à dispersão de poluentes na água, resolve-se a equação advectivo-difusiva bidimensional com coeficientes variáveis, realizando uma mudança de variáveis de modo a reescrevê-la em termos do potencial velocidade e da função corrente correspondentes ao respectivo escoamento potencial, estendendo a solução para domínios de contornos arbitrários Na aplicação referente ao problema de dispersão de poluentes na atmosfera, realiza-se uma mudança de variáveis de modo a obter uma equação diferencial parcial com coeficientes constantes na qual se possam aplicar as regras de manipulação de exponenciais de operadores diferenciais. Os resultados numéricos obtidos são comparados com dados disponíveis na literatura. Diversas vantagens da aplicação das formulações apresentadas podem ser citadas, a saber, o aumento da velocidade de processamento, permitindo a obtenção de solução em tempo real; a redução da quantidade de memória requerida na realização de operações necessárias para a obtenção da solução analítica; a possibilidade de dispensar a discretização do domínio em algumas situações.

Solução da equação de difusão unidimensional transiente para o estudo da dispersão de poluentes na camada limite planetária

Neste trabalho apresenta-se uma solução analítica para a dispersão vertical turbulenta em uma Camada Limite Convectiva e em uma Camada Limite Estável. A equação analisada considera a difusão com velocidades finitas, o que representa o transporte turbulento fisicamente correto. Considerando o caráter não-local, adicionam-se na equação que representa uma fonte área instantânea, termos como: o tempo de relaxação, a assimetria, a escala de tempo Lagrangeana e a velocidade turbulenta vertical. A solução é obtida utilizando-se a técnica da Transformada de Laplace. Os parâmetros que encerram a turbulência são derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade. Foram utilizados coeficientes de difusão especáficos para cada uma das camadas. A transformada inversa é obtida através do esquema numérico de quadratura Gaussiana. São apresentadas várias simulações para diferentes alturas de fonte área e obtém-se o valor da concentração para alturas próximas ao solo e próximas ao topo da Camada Limite Planetária. A inserção do termo de contra-gradiente na equação resultou em uma pequena influência na concentração de poluentes, observada de forma mais expressiva na Camada Limite Convectiva.

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.

Solução analítica da equação da energia estacionária e bidimensional para simulação de escoamento plenamente desenvolvido em placa paralela pelo método da giltt

Neste trabalho e apresentado um avanço na tecnica GILTT(Generalized Integral and Laplace Transform Technique) solucionando analiticamente um sistema de EDO's(Equações Diferenciais Ordinarias) de segunda ordem resultante da transformação pela GITT(Generalized Integral Transform Technique). Este tipo de problema usualmente aparece quando esta tecnica é aplicada na solução de problemas bidimensionais estacionários. A principal idéia consiste na redução de ordem do problema transformado em outro sistema de EDO's lineares de primeira ordem e a solução analítica deste problema, pela técnica da transformada de Laplace. Como exemplo de aplicação é resolvida a equação da energia linear bidimensional e estacionária. São apresentadas simulações numéricas e comparações com resultados disponíveis na literatura.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Equação de burgers propriedades e comportamento assintótico

No presente trabalho, obtemos e analisamos diversas propriedades das soluções u(·, t) da equação de difusão linear (equação do calor em meios unidimensionais homogêneos) ut = μuxx x 2 R, t > 0 correspondentes a estados iniciais u(x, 0) = u0(x), com u0 2 Lp(R), para algum 1 p < 1; bem como da equação de Burgers ut + cuux = μuxx x 2 R, t > 0 onde c, μ são constantes dadas, sendo c 6= 0 e μ > 0 e ainda assumindo u(x, 0) = u0(x) com u0 2 Lp(R) para 1 p < 1, e limitado. Estudamos também a equação mais geral da forma ut + f(u)x = μuxx x 2 R, t > 0 discutindo várias propriedades importantes das soluções, associadas a estados iniciais u0 2 Lp(R) \ L1(R) para algum 1 p < 1. Em particular, examinamos o comportamento de ku(·, t)kLr(R), p r 1, para t >> 1, e diversas propriedades relacionadas.

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.

Alguns resultados para equação do calor e equações de advecção-difusão

Neste trabalho, são obtidas diversas propriedades (em especial, referentes ao comportamento ao t -+ +00) das soluções u(', t) da equação linear do calor, Ut = div(AV'u), x E JRn, t > O onde A E JRnxné uma matriz constante simétrica e positiva definida, correspondentes a estados iniciais p-somáveis, i.e., u(x, O) = uo(x), Uo E LP(JRn), onde 1 :::;p < 00. Em particular, é examinado o comportamento de Ilu(., t)IILP(lRn) ao t -+ +00, mostrando-se que Ilu(., t)IILl(lRn)-+ Ikn u(x, O)dXI quando p = 1, e Ilu(-' t)IILP(lRn)-+ O quando p > 1. São analisadas, também, as taxas de decaimento e o comportamento assintótico das soluções u(', t) de equações de advecção-difusão da forma Ut + divf(u) = div(A(u)V'u), x E JRn, t > O correspondentes a estados iniciais p-somáveis e limitados, i.e., u(x, O)= uo(x), u(', O) E LP(JRn) n LOO(JRn), onde 1 :::;p :::; 2. Novamente, é examinado o comportamento de Ilu(" t)IILP(lRn)ao t -+ +00, mostrando-se que Ilu(., t)IILl(lRn)-+ Ikn u(x, O)dxl quando p = 1, e Ilu(" t)IILP(lRn)-+ O quando p > 1. Várias outras propriedades importantes são também discutidas, seguindo principalmente [Silva, 2003], [Crandall e Tartar, 1980], [Hagstrom et al., 2003], [Zingano, 1999], [Zingano, 2004a], [Zingano, 2004b].