Simulação do processo de transferência de calor não linear condução-radiação por meio de um esquema linear em diferenças finitas.

Neste trabalho o processo não linear de transmissão de calor condução-radiação é abordado num contexto bidimensional plano e simulado com o uso de um esquema linear em diferenças finitas. O problema original é tratado como o limite de uma sequencia de problemas lineares, do tipo condução-convecção. Este limite, cuja existência é comprovada, é facilmente obtido a partir de procedimentos básicos, accessíveis a qualquer estudante de engenharia, permitindo assim o emprego de hipóteses mais realistas, já que não se tem o limitante matemático para a abordagem numérica de uma equação diferencial parcial elíptica. Neste trabalho foi resolvido o problema de condução de calor em regime permanente em uma placa com condições de contorno convectivas e radioativas utilizando-se o software MatLab, vale ressaltar, que a mesma metodologia é aplicável para geometrias mais complexas.

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.

Análise comparativa entre simuladores de linha de fluxo e diferenças finitas para um reservatório de petróleo submetido a injeção de água

Water injection is the most widely used method for supplementary recovery in many oil fields due to various reasons, like the fact that water is an effective displacing agent of low viscosity oils, the water injection projects are relatively simple to establish and the water availability at a relatively low cost. For design of water injection projects is necessary to do reservoir studies in order to define the various parameters needed to increase the effectiveness of the method. For this kind of study can be used several mathematical models classified into two general categories: analytical or numerical. The present work aims to do a comparative analysis between the results presented by flow lines simulator and conventional finite differences simulator; both types of simulators are based on numerical methods designed to model light oil reservoirs subjected to water injection. Therefore, it was defined two reservoir models: the first one was a heterogeneous model whose petrophysical properties vary along the reservoir and the other one was created using average petrophysical properties obtained from the first model. Comparisons were done considering that the results of these two models were always in the same operational conditions. Then some rock and fluid parameters have been changed in both models and again the results were compared. From the factorial design, that was done to study the sensitivity analysis of reservoir parameters, a few cases were chosen to study the role of water injection rate and the vertical position of wells perforations in production forecast. It was observed that the results from the two simulators are quite similar in most of the cases; differences were found only in those cases where there was an increase in gas solubility ratio of the model. Thus, it was concluded that in flow simulation of reservoirs analogous of those now studied, mainly when the gas solubility ratio is low, the conventional finite differences simulator may be replaced by flow lines simulator the production forecast is compatible but the computational processing time is lower.

Análise de métodos numéricos de diferenças finitas para solução da equação de Poisson em domínios irregulares

Pós-graduação em Matematica Aplicada e Computacional - FCT

Modelagem sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Modelagem e imageamento 2.5D no domínio do tempo através de diferenças finitas

A modelagem 2.5D consiste em simular a propagação do campo de ondas em 3D em meios com simetria de translação em uma direção. Nesta tese esta abordagem é formulada para meios elásticos e anisotrópicos com classe de simetria arbitrária e a geometria de aquisição não precisa coincidir com um plano de simetria do meio. A migração por reversão no tempo do campo de ondas é formulada e implementada através de diferenças finitas 2.5D. Para reduzir os efeitos de retro-espalhamento e melhorar a recuperação da amplitude dos eventos migrados, propomos uma nova condição de imagem para migração reversa no tempo baseada na análise assintótica da condição de imagem clássica por correlação cruzada. Experimentos numéricos indicam que a migração reversa no tempo 2.5D com a nova condição de imagem proposta, melhora a resolução da imagem em relação à migração reversa no tempo 2D e reduz acentuadamente os ruídos causados por retro-espalhamento.

Modelagem acústica por diferenças finitas e elementos finitos em 2-D e 2,5-D

A modelagem acústica fornece dados úteis para avaliação de metodologias de processamento e imageamento sísmico, em modelos com estrutura geológica complexa. Esquemas de diferenças finitas (DF) e elementos finitos (EF) foram implementados e avaliados em modelos homogêneos e heterogêneos. O algoritmo de diferenças finitas foi estendido para o caso 2,5-D em modelos com densidade variável. Foi apresentada a modelagem de alvos geológicos de interesse exploratório existentes na Bacia Paleozóica do Solimões na Amazônia. Reflexões múltiplas de longo período produzidas entre a superfície livre e a discordância Cretáceo-Paleozóica, a baixa resolução da onda sísmica nas proximidades do reservatório e as fracas reflexões na interface entre as rochas reservatório e as rochas selantes são as principais características dos dados sintéticos obtidos, os quais representam um grande desafio ao imageamento sísmico.

Modelagem de perfis de indução no domínio do tempo por diferenças finitas

Apresentamos aqui uma metodologia alternativa para modelagem de ferramentas de indução diretamente no domínio do tempo. Este trabalho consiste na solução da equação de difusão do campo eletromagnético através do método de diferenças finitas. O nosso modelo consiste de um meio estratificado horizontalmente, através do qual simulamos um deslocamento da ferramenta na direção perpendicular às interfaces. A fonte consiste de uma bobina excitada por uma função degrau de corrente e o registro do campo induzido no meio é feito através de uma bobina receptora localizada acima da bobina transmissora. Na solução da equação de difusão determinamos o campo primário e o campo secundário separadamente. O campo primário é obtido analiticamente e o campo secundário é determinado utilizando-se o método de Direção Alternada Implícita, resultando num sistema tri-diagonal que é resolvido através do método recursivo proposto por Claerbout. Finalmente, determina-se o valor máximo do campo elétrico secundário em cada posição da ferramenta ao longo da formação, obtendo-se assim uma perfilagem no domínio do tempo. Os resultados obtidos mostram que este método é bastante eficiente na determinação do contato entre camadas, inclusive para camadas de pequena espessura.

Modelagem sísmica por diferenças finitas em meios bidimensionais com difratores

Este trabalho tem por objetivo a modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações. Como parte deste estudo foi feita a análise numérica visando as condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas. Para a validação do programa de diferenças finitas foi feita a comparação cinemática com a Teoria do Raio para um modelo simples. O motivo deste estudo é ter uma melhor compreensão e controle sobre os problemas de modelagem, visando contribuir para a solução de problemas na interpretação de dados sísmicos. Segundo vários autores na literatura geológica, Derby (1877), Evans (1906), Paiva (1929) e Moura (1938). A Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o recente, atingindo espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançando espessuras de centenas de metros, perfazendo um volume total de 90.000 quilômetros cúbicos. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica, impossibilitando uma melhor interpretação dos horizontes refletores que se encontram abaixo das soleiras. Para representar situações geológicas desse tipo foi usado um modelo acústico de velocidades. Para o cálculo dos sismogramas foi utilizado um programa de diferenças finitas com aproximação de quarta ordem da equação da onda acústica no espaço e no tempo. As aplicações dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas têm melhorado a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observado a influência significativa das múltiplas devido a camada de alta velocidade, o que faz com que haja maior perda de energia dificultando a interpretação dos alvos. Por esta razão recomendo a integração de dados de superfície com os dados de poço, com o objetivo de se obter uma melhor imagem dos alvos abaixo das soleiras de diabásio.

Estudo da deposição de energia das radiofrequências através do método das diferenças finitas dependentes do tempo Monografia

Radiotherapy in veterinary practice is already known and widely distributed in large specialized centers of developed countries. In early 2000, there were about 30 radiotherapy equipment specifically designed for the veterinary clinic in the United States. In Brazil, the veterinary radiotherapy is still confined to research in universities, where most of the procedures is radiation therapy performed with superficial x-ray machines, with a voltage between 50 and 150 kVp, focus-distance surface (DFS) between 20,0cm and 40,0cm. As that occurs in human medicine, new research strengthens the development and prospects for the use of radiotherapy as a safe option for treating cancer in animals. This paper presents a methodology for calculating the exposure time for superficial radiotherapy procedures in veterinary medicine for small animals (dogs and cats). The dosimetric parameters of X-rays are determined using a spreadsheet tool for Microsoft Office Excel, developed in this paper for a device Dermopan 2, Veterinary Hospital of UNESP in Araçatuba. Using the worksheet helps the veterinarian to determine the time of exposure to radiation determined for each clinical case, optimize the workflow for professionals in veterinary radiotherapy procedures, which often lack the medical physics in team and at the time of radiotherapy. The correct use of spreadsheet decreases the chances of errors in dose rates of radiation, providing a higher quality of care

Modelos dinâmicos simplificados, para verificação do comportamento de rede de satélites amarrados por cabos flexíveis, usando-se o método explícito das diferenças finitas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise comparativa entre simuladores de linha de fluxo e diferenças finitas para um reservatório de petróleo submetido a injeção de água

Water injection is the most widely used method for supplementary recovery in many oil fields due to various reasons, like the fact that water is an effective displacing agent of low viscosity oils, the water injection projects are relatively simple to establish and the water availability at a relatively low cost. For design of water injection projects is necessary to do reservoir studies in order to define the various parameters needed to increase the effectiveness of the method. For this kind of study can be used several mathematical models classified into two general categories: analytical or numerical. The present work aims to do a comparative analysis between the results presented by flow lines simulator and conventional finite differences simulator; both types of simulators are based on numerical methods designed to model light oil reservoirs subjected to water injection. Therefore, it was defined two reservoir models: the first one was a heterogeneous model whose petrophysical properties vary along the reservoir and the other one was created using average petrophysical properties obtained from the first model. Comparisons were done considering that the results of these two models were always in the same operational conditions. Then some rock and fluid parameters have been changed in both models and again the results were compared. From the factorial design, that was done to study the sensitivity analysis of reservoir parameters, a few cases were chosen to study the role of water injection rate and the vertical position of wells perforations in production forecast. It was observed that the results from the two simulators are quite similar in most of the cases; differences were found only in those cases where there was an increase in gas solubility ratio of the model. Thus, it was concluded that in flow simulation of reservoirs analogous of those now studied, mainly when the gas solubility ratio is low, the conventional finite differences simulator may be replaced by flow lines simulator the production forecast is compatible but the computational processing time is lower.

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.

Determinação da resposta dinâmica de feixes de condutores em linhas aéreas de transmissão

Até os dias de hoje, a abordagem estática tem sido a usualmente empregada para a avaliação da resposta de cabos de Linhas de Transmissão (LTs), apesar da resposta dinâmica ser, em muitos casos, reconhecidamente importante na avaliação do desempenho dos cabos. Para uma análise completa, não se pode desconhecer a natureza dinâmica da maioria dos fenômenos a que as LTs estão submetidas, sendo exemplos típicos as excitações mecânicas causadas pela ação do vento (carregamentos de baixa e alta freqüência) e a ruptura de cabos. Cabe salientar que mesmo a análise estática empregada tem sido simplificada, considerando normalmente apenas casos contemplados com soluções analíticas. A justificativa para emprego da análise estática baseia-se no menor esforço numérico exigido. Entretanto, essa justificativa já não se sustenta dado os grandes avanços na área computacional, possibilitando o estudo do desempenho de LTs na ocorrência de fenômenos que provocam carregamentos dinâmicos. Adicionalmente, há uma crescente demanda para que se obtenha um melhor entendimento de muitas questões relativas ao comportamento dinâmico das LTs, que possam explicar desempenhos observados e/ou sustentar o desenvolvimento de novas alternativas mais arrojadas de projeto e construção Entre os diversos métodos utilizados na engenharia estrutural, a integração direta das equações do movimento, através de métodos numéricos como as diferenças finitas centrais, se constitui numa poderosa ferramenta de cálculo. Esta ferramenta possibilita o tratamento de problemas envolvendo não linearidades geométricas e do material, como são os casos onde se avalia a resposta de cabos suspensos submetidos a carregamentos variáveis no tempo. Esta pesquisa objetiva a aplicação do método da integração direta das equações do movimento na análise de feixes de cabos de LT, quando submetidos à ação de carregamentos mecânicos variáveis no tempo, principalmente à excitação de ventos oriundos de fenômenos com natureza complexa (tormentas elétricas, por exemplo). É apresentado um método para determinar a resposta dinâmica de feixes que considera a interação entre o vento incidente e o movimento do condutor. A solução é obtida por integração numérica, no domínio do tempo, das equações de movimento de um modelo tridimensional não-linear discreto de feixe, as quais definem as forças nos espaçadores e cabos através de coeficientes aerodinâmicos obtidos experimentalmente. Também são apresentados modelos dos fenômenos meteorológicos mais comuns em nosso país (Brasil): tormenta extratropical (EPS) e tormenta elétrica (TS) Como ilustração, são apresentados exemplos de modelagem de vãos de LTs com condutor singelo e com feixes de condutores. Os exemplos demonstram a capacidade de avaliação de Estados Limites relacionados à distância relativa entre subcondutores, à estabilidade do feixe, à representação das suas propriedades e ao comportamento dinâmico, bem como aos carregamentos transmitidos às estruturas. A análise emprega conhecidas relações constitutivas para representar o comportamento tensão-deformação dos cabos. O enfoque utilizado possibilita a avaliação mais precisa de casos reais que ainda não podiam ser convenientemente tratados, além de permitir a extensão para estudos bem mais complexos, tais como feixes com disposições assimétricas.

Um estudo da otimização da geometria de um pára-quedas simplificado

Neste trabalho desenvolve-se um estudo numérico do fluxo de ar em torno da geometria de um pára-quedas tradicional simplificado, para alguns valores de Reynolds. O método baseia-se na solução das equações incompressíveis de Navier- Stokes discretizadas pelo método de diferenças finitas e integradas pelo método de Runge-Kutta. Utiliza-se o método dos contornos virtuais para representar a geometria numa malha cartesiana e o método de otimização não-linear dos poliedros flexíveis para otimização do coeficiente de arraste calculado através do código de dinâmica de fluidos computacional; esteé um método de busca multivariável, onde o pior vértice de um poliedro com n + 1 vérticesé substituído por um novo.