Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.

Adaptação via movimento de malhas em escoamentos compressíveis

Resumo não disponível.

Análise e modelagem do escoamento transitório de fluidos de perfuração

Drilling fluids present a thixotropic behavior and they usually gel when at rest. The sol-gel transition is fundamental to prevent the deposit of rock fragments, generated by drilling the well, over the drill bit during eventual stops. Under those conditions, high pressures are then required in order to break-up the gel when circulation is resumed. Moreover, very high pressures can damage the rock formation at the bottom of the well. Thus, a better understanding of thixotropy and the behavior of thixotropic materials becomes increasingly important for process control. The mechanisms that control thixotropy are not yet well defined and modeling is still a challenge. The objective of this work is to develop a mathematical model to study the pressure transmission in drilling fluids. This work presents a review of thixotropy and of different mathematical models found in the literature that are used to predict such characteristic. It also shows a review of transient flows of compressible fluids. The problem is modeled as the flow between the drillpipe and the annular region (space between the wall and the external part of the drillpipe). The equations that describe the problem (mass conservation, momentum balance, constitutive and state) are then discretized and numerically solved by using a computational algorithm in Fortran. The model is validated with experimental and numerical data obtained from the literature. Comparisons between experimental data obtained from Petrobras and calculated by three viscoplastic and one pseudoplastic models are conducted. The viscoplastic fluids, due to the yield stress, do not fully transmit the pressure to the outlet of the annular space. Sensibility analyses are then conducted in order to evaluate the thixotropic effect in pressure transmission.

Otimização de geometrias aerodinâmicas utilizando métodos inversos

O objetivo deste trabalho é a obtenção de uma técnica para a modelagem otimizada de corpos submetidos a fluxos de alta velocidade, como aerofólios em escoamentos transônicos e outras geometrias aerodinâmicas. A técnica é desenvolvida através de expansões em séries de Fourier para um conjunto de equações diferenciais com interrelação com as condições de contorno, sendo uma equação para a parte superior e outra para a parte inferior do aerofólio. O método de integração temporal empregado baseia-se no esquema explícito de Runge-Kutta de 5 estágios para as equações da quantidade de movimento e na relação de estado para a pressão. Para a aproximação espacial adota-se um esquema em volumes finitos no arranjo co-localizado em diferenças centrais. Utiliza-se dissipação artificial para amortecer as frequências de alta ordem do erro na solução das equações linearizadas. A obra apresenta a solução de escoamentos bi e tridimensionais de fluidos compressíveis transônicos em torno de perfis aerodinâmicos. Os testes num´ericos são realizados para as geometrias do NACA 0012 e 0009 e asas tridimensionais usando as equações de Euler, para número de Mach igual a 0.8 e ® = 0o. Os resultados encontrados comparam favoravelmente com os dados experimentais e numéricos disponíveis na literatura.

Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).

Análise das dosagens do fator de crescimento endotelial vascular no plasma e fluidos peritoneais de pacientes com câncer epitelial de ovário

A carcinogênese epitelial ovariana tem sido foco de estudos científicos em todo o mundo desenvolvido. A angiogênese tumoral ovariana é um processo multifatorial que resulta em vários produtos pró-angiogênicos. Entre eles, o fator de crescimento vascular endotelial (VEGF) é predominante. Os objetivos deste estudo foram relacionar as dosagens do VEGF dos fluidos peritoneais, do plasma periférico e do infundíbulo pélvico aos níveis de citorredução em pacientes operadas de adenocarcinoma epitelial de ovário (CEO); formular um modelo probabilístico de citorredução e utilizar estas dosagens para estimar a probabilidade do desfecho de citorredução. Além disto, foi criada uma nova variável chamada carga de VEGF. Pelo procedimento step-wise a citorredução foi melhor descrita pela carga de VEGF, mas faltou Normalidade aos resíduos, não sendo possível a adequação de um modelo matemático. Porém, a curva ROC, forneceu uma área sob a curva de 0,84, com sensibilidade de 71,4 % e especificidade variando de 69,5 a 73,9%. O ponto de corte ótimo foi 15,52 log de picograma de carga de VEGF. A odds-ratio (OR) calculada para citorredução ótima descrita pela carga de VEGF foi de 11 (IC= 2,59 ; 46,78). No grupo com estágio avançado (III e IV), a OR foi de 6 (IC= 1,15; 31,22). Apesar do pequeno número de casos, esta nova variável pode vir a ser um auxiliar na determinação de situações onde cirurgia citorredutora deixa de ser a pedra fundamental do tratamento primário do CEO e a indução quimioterápica passe a ter o principal papel na citorredução química antes da cirugia nestes casos.

Sobresaturación y nucleación primaria en fluidos biológicos respecto a la hidroxiapatita.

Aunque en las últimas décadas se han producido considerables progresos en el diagnóstico y el tratamiento de la litiasis renal, la compresión de su patogenia no ha ido en paralelo con estos avances. De hecho, a pesar de los grandes esfuerzos en investigación y de un considerable número de teorías, aún se desconoce la secuencia exacta de acontecimientos que llevan a la formación de los cálculos renales. En este trabajo se realiza un estudio de los principales tipos de cálculos renales y sus posibles causas. Se estudia, concretamente, la formación de la hidroxiapatita, un cálculo renal inorgánico relativamente frecuente. Se utiliza, para tal fin, el software comercial Medusa.

Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Inflação em sistemas com fluidos de radiação com viscosidade e dissipação e análise de estabilidade

Este trabalho apresenta um estudo da estabilidade das equações da inflação morna com um fluido de radiação viscoso. A viscosidade do fluido é proveniente do constante decaimento de partículas neste, devido à dissipação do campo escalar da inflação, o ínflaton.Esta viscosidade, que pode ser volumar ou laminar, é tratada em termos de teorias termodinâmicas fora do equilíbrio. Este estudo se limita às equações de fundo da inflação morna, de modo que somente a viscosidade volumar tem um efeito significativo, sendo a viscosidade laminar importante somente no contexto de perturbações cosmológicas. A descrição da viscosidade em termos de uma termodinâmica fora do equilíbrio, porém, não pode ser realizada univocamente, pois a única informação que temos sobre processos irreversíveis é a segunda lei da termodinâmica. Portanto, parte-se em busca de teorias que estejam de acordo com esta lei e que, por argumentos plausíveis, sejam capazes de descrever o comportamento dos fluxos dissipativos próximo ao equilíbrio. O objetivo deste trabalho é estudar a estabilidade da inflação morna viscosa para teorias causais e não causais para o fluido de radiação com viscosidade, de forma que se possa observar o impacto da viscosidade no regime inflacionário e a relevância de se passar a considerar a causalidade. Para o fluido de radiação, as teorias consideradas são a teoria não causal de Eckart e as teorias causais de Israel-Stewart e de Denicol et al (hidrodinâmica dissipativa causal não linear). Obtém-se que as teorias causais, como era de se esperar, além de serem, por definição, consistentes no tocante à finitude da velocidade de propagação dos fluxos dissipativos, tornam o sistema dinâmico estável para valores de viscosidade mais distantes do equilíbrio. Observa-se também, nitidamente, que a teoria de Denicol et al é a mais robusta nesse sentido. Este trabalho, portanto, visa dar continuidade ao estudo dos efeitos não-isentrópicos na inflação, já que, além da dissipação do ínflaton na inflação morna, o impacto da viscosidade tem despertado bastante interesse.

Modelo simple para la permeabilidad dinámica de fluidos viscoelásticos

In the present work a simple form to obtain analytical expression for the dynamic permeability of Maxwellian fluids is presented. This expression gives the frequency dependent form of this dynamic permeability. In particular case, the analytic expression for the sinusoidal pressure pump fluid is illustrated in the configuration space. As an example of the feasibility of this expression the flow of human blood in a tube is presented finding that the human heart frequency has the same order that the frequencies where the dynamic permeability shows resonances. In order to make clear the above aspect of the dynamic permeability a model of pulsing pressure drops (gaussian like) are analyzed.

Sondagens mecânicas e prospecção geofísica na caracterização de fluidos

Esta dissertação aborda a utilização de métodos geofísicos de superfície e de profundidade na caracterização de fluidos, compreendendo várias abordagens distintas, consoante o tipo de temática em análise. Assim, o assunto proposto para a presente dissertação relaciona-se com a utilização conjunta de sondagens mecânicas e prospecção geofísica na caracterização de fluidos através da exploração de três temáticas distintas, designadamente a Hidrogeologia, a Geotermia e o Ambiente. O presente trabalho de dissertação baseado na experiência profissional da autora, é feito de acordo com o estabelecimento de uma linha condutora, fundamentada e estruturada, para uma visão global do modo como a interpretação integrada das sondagens mecânicas e dos métodos geofísicos dá um forte contributo para a melhoria do conhecimento de uma região, seja ele geológico, hidrogeológico ou ambiental. Os temas a abordar foram seleccionados não só pela sua indiscutível importância, devido ao interesse estratégico que actualmente suscitam, como também pelos diferentes modos de processar e interpretar a informação geofísica adquirida consoante os objectivos da mesma, de forma a salientar particularidades na respectiva interpretação e integração de dados.

Desenvolvimento de padrão gravimétrico de medição de caudal de fluidos entre 10 mL/min e 100 nL/min

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

Desarrollo de un método analítico para la determinación de aflatoxinas M1 y B1 en fluidos biológicos empleando cromatografía de líquidos de alta resolución

Tesis (Maestría en Ciencias con Especialidad en Química Analítica Biomédica) UANL