Formulação teórica dos fundamentos da otimização global topográfica com análise de desempenho e aplicações à estabilidade de fases de misturas termodinâmicas

Métodos de otimização que utilizam condições de otimalidade de primeira e/ou segunda ordem são conhecidos por serem eficientes. Comumente, esses métodos iterativos são desenvolvidos e analisados à luz da análise matemática do espaço euclidiano n-dimensional, cuja natureza é de caráter local. Consequentemente, esses métodos levam a algoritmos iterativos que executam apenas as buscas locais. Assim, a aplicação de tais algoritmos para o cálculo de minimizadores globais de uma função não linear,especialmente não-convexas e multimodais, depende fortemente da localização dos pontos de partida. O método de Otimização Global Topográfico é um algoritmo de agrupamento, que utiliza uma abordagem baseada em conceitos elementares da teoria dos grafos, a fim de gerar bons pontos de partida para os métodos de busca local, a partir de pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem dois objetivos. O primeiro é realizar uma nova abordagem sobre método de Otimização Global Topográfica, onde, pela primeira vez, seus fundamentos são formalmente descritos e suas propriedades básicas são matematicamente comprovadas. Neste contexto, propõe-se uma fórmula semi-empírica para calcular o parâmetro chave deste algoritmo de agrupamento, e, usando um método robusto e eficiente de direções viáveis por pontos-interiores, estendemos o uso do método de Otimização Global Topográfica a problemas com restrições de desigualdade. O segundo objetivo é a aplicação deste método para a análise de estabilidade de fase em misturas termodinâmicas,o qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. A solução deste problema de otimização global é necessária para o cálculo do equilíbrio de fases, que é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Além disso, afim de ter uma avaliação inicial do potencial dessa técnica, primeiro vamos resolver 70 problemas testes, e então comparar o desempenho do método proposto aqui com o solver MIDACO, um poderoso software recentemente introduzido no campo da otimização global.

Otimização global topográfica aplicada a problemas de otimização com restrições de igualdade e desigualdade

Os métodos de otimização que adotam condições de otimalidade de primeira e/ou segunda ordem são eficientes e normalmente esses métodos iterativos são desenvolvidos e analisados através da análise matemática do espaço euclidiano n-dimensional, o qual tem caráter local. Esses métodos levam a algoritmos iterativos que são usados para o cálculo de minimizadores globais de uma função não linear, principalmente não-convexas e multimodais, dependendo da posição dos pontos de partida. Método de Otimização Global Topográfico é um algoritmo de agrupamento, o qual é fundamentado nos conceitos elementares da teoria dos grafos, com a finalidade de gerar bons pontos de partida para os métodos de busca local, com base nos pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem como objetivo a aplicação do método de Otimização Global Topográfica junto com um método robusto e eficaz de direções viáveis por pontos-interiores a problemas de otimização que tem restrições de igualdade e/ou desigualdade lineares e/ou não lineares, que constituem conjuntos viáveis com interiores não vazios. Para cada um destes problemas, é representado também um hiper-retângulo compreendendo cada conjunto viável, onde os pontos amostrais são gerados.

Um método SN híbrido direto para cálculos de sistemas combustível-moderador em geometria unidimensional

Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto.

A mobilidade eletroforética e o perfil de potencial da membrana celular

O objetivo do presente trabalho foi estudar o comportamento dos potenciais superficiais e do perfil de potencial atraves da membrana de eritr ocito em func ao da forca i onica e das cargas superficiais, usando um modelo que leva em conta as cargas el etricas do glicoc alix e das proteınas citoplasm aticas, al em das cargas superficiais da bicamada lipıdica e os efeitos dos eletr olitos divalentes. Programas especıficos em linguagem C foram elaborados para o c alculo desses potenciais, tomando como dados num ericos resultados experimentais de medidas de mobilidade eletrofor etica de eritr ocitos para diferentes valores de forca i onica. Neste c alculo, o metodo para tratamento dos dados eletrofor eticos indicado por Hsu et al.[57] foi incluıdo em nosso modelo. A equac ao de Poisson-Boltzmann nao linear foi resolvida por computac ao num erica, usando o metodo de Runge-Kutta de quarta ordem, obtendo-se os perfis de potencial. Os resultados mostraram que a estimativa da densidade de carga el etrica na superfıcie de c elulas usando a equac ao cl assica de Helmholtz-Smoluchowski conduz a valores que nao conseguem refletir as forcas que regem o comportamento eletrofor etico das mesmas. O presente modelo gerou valores de potenciais superficiais e perfis de potencial para a membrana do eritr ocito bem distintos daqueles obtidos anteriormente para um modelo descrito por uma equac ao de Poisson-Boltzmann linear. Nossos resultados confirmam que a avaliac ao de parametros el etricos superficiais da membrana de eritr ocito, envolvendo dados oriundos de eletroforese, deve incluir c alculos hidrodin amicos al em de eletroest aticos, como sugerido por Hsu et al. [57].

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.