Técnica para dedução do modelo algébrico PTT, aplicações e análises numéricas em escoamentos bidimensionais

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Transporte hidráulico de misturas bentonita-água em condutos forçados

A argila bentonítica é amplamente utilizada em transporte de sólidos produzidos durante perfuração de poços. Teve por objetivo estudar o escoamento de misturas bentonita-água e determinar suas propriedades reológicas e parâmetros hidráulicos úteis nos projetos de instalações de recalque de misturas sólido-líquido. Foi montado um circuito fechado de tubulações para estudar dados de perda de carga e perfis de velocidade. Realizaram-se ensaios com misturas bentonita-água sob varias concentrações, algumas transportando areia. Observaram-se que a reologia da mistura bentonita-água é melhor descrita pela formulação de Herschell-Bulkley para fluidos não-Newtonianos. O coeficiente de atrito para descrever a perda de carga da mistura bentonita-água observada em tubulações no laboratório coloca-se entre as previsões de Tomita (1959) e Szilas et al (1981). A variação da velocidade da mistura na seção transversal do tubo é melhor aproximada pela equação de Bogue-Metzner (1963).

LAS ERUPCIONES DEL VOLCAN ARENAL, COSTA RICA

En este articulo hemos considerado el Volcán Arenal, como un ejemplo típico en America Central del vulcanismo explosivo basalticoandesitico. Presentamos en la primera parte una teoría general de los procesos que se encuentran  en el volcanismo explosivo, distinguiendo en primer lugar los productos de las erupciones, y en segundo lugar la combinación de los hechos eruptivos con las leyes de la mecánica de los fluidos. Esta combinación es muy eficiente para explicar las estructuras físicas de las erupciones y al mismo tiempo revela una gran complejidad de los problemas geofísicos, geoquímicos y fluido-mecánicos, que se son el fundamento del volcanismo. En la segunda parte parte enfocamos la actividad del Arenal en el contexto de la tectónica regional, haciendo al mismo tiempo un análisis detallado de la evolución de los cráteres y migración de la actividad. La estructura y dinámica de las varias unidades de flujo que constituyen el campo lávico  fueron explicadas de acuerdo a los principios aplicables a los fluidos no newtonianos andesitico y, en particular, de este volcán.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Simulação de escoamentos incompressíveis não newtonianos em dutos com expansão brusca

O objetivo deste trabalho é a simulação numérica de escoamentos incompressíveis bidimensionais em dutos com expansão brusca, considerando o raio de expansão de 3 : 1. As equações governantes implementadas são as de Navier, que junto com relações constitutivas para a tensão visam representar comportamentos não newtonianos. A integração temporal é feita usando o esquema explícito de Runge-Kutta com três estágios e de segunda ordem; as derivadas espaciais são aproximadas pelo método de diferenças finitas centrais. Escoamentos em expansões bruscas para fluidos newtonianos apresentam um número de Reynolds crítico, dependente do raio de expansão, na qual três soluções passam a ser encontradas: uma solução sim étrica instável e duas soluções assimétricas rebatidas estáveis. Aumentando o número de Reynolds, a solução passa a ser tridimensional e dependente do tempo. Dessa forma, o objetivo é encontrar as diferenças que ocorrem no comportamento do fluxo quando o fluido utilizado possui características não newtonianas. As relações constitutivas empregadas pertencem à classe de fluidos newtonianos generalizados: power-law, Bingham e Herschel-Bulkley. Esses modelos prevêem comportamentos pseudoplásticos e dilatantes, plásticos e viscoplásticos, respectivamente. Os resultados numéricos mostram diferenças entre as soluções newtonianas e não newtonianas para Reynolds variando de 30 a 300. Os valores de Reynolds críticos para o modelo power-law não apresentaram grandes diferenças em comparação com os da solução newtoniana. Algumas variações foram percebidas nos perfis de velocidade. Entretanto, os resultados obtidos com os modelos de Bingham e Herschel-Bulkley apresentaram diferenças significativas quando comparados com os newtonianos com o aumento do parâmetro adimensional Bingham; à medida que Bingham é aumentado, o tamanho dos vórtices diminui. Além disso, os perfis de velocidade apresentam diferenças relevantes, uma vez que o fluxo possui regiões onde o fluido se comporta como sólido.

Influência da faixa de taxas de deformação na precisão de caracterizações reológicas de fluidos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Coeficientes de perda de carga em válvulas e acessórios no escoamento de fluídos não-newtonianos puramente viscosos e viscoelásticos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Supertareas newtonianas : la relación del indeterminismo con la pérdida de la energía en sistemas newtonianos anómalos

[ES]La tesis tiene como objetivo profundizar en el conocimiento de la relación que existe entre el indeterminismo y el fallo del principio de conservación de la energía en la mecánica newtoniana bajo procesos que consisten en una sucesión infinita de colisiones perfectamente elásticas, llamados supertareas newtonianas. Se propone una supertarea, que nombramos STMNC, que muestra no ser conservativa pero que, bajo los mecanismos hasta ahora conocidos, no es manifiestamente indeterminista. Se demuestra, con un planteamiento general que llamamos G¿1, que tal proceso es de hecho indeterminista antes de que la supertarea llegue a ejecutarse; no obstante, se señalan las dificultades para que el proceso sea indeterminista una vez que la supertarea concluye, es decir, una vez que la energía se pierde. Con esto, y tras el análisis de las fuentes del indeterminismo y la pérdida de la energía, se sugiere que el único aspecto en común entre ambas anomalías, cuando emergen bajo estos procesos, es la infinidad de colisiones que involucran. Para obtener estos resultados se sigue la siguiente estrategia. En la primera parte, se analizan los argumentos en torno a las paradojas de Zenón, unos a favor y otros en contra de las ideas de espacio y tiempo que permiten la ejecución de las supertareas. En la segunda parte, se analizan los argumentos que amenazan a los procesos que tratamos de ser inconsistentes y de no ser newtonianos. En conjunto, con estas dos partes garantizamos que las supertareas newtonianas son procesos legítimos, brindando así sentido a nuestro resultado principal. En la tercera parte, se analizan las diversas fuentes del indeterminismo en la mecánica newtoniana bajo supertareas y la relación que dichas fuentes guardan con el indeterminismo.

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).

Desenvolvimento de um simulador numérico empregando o método Smoothed Particle Hydrodynamics para a resolução de escoamentos incompressíveis. Implementação computacional em paralelo (CUDA)

Neste trabalho, foi desenvolvido um simulador numérico baseado no método livre de malhas Smoothed Particle Hydrodynamics (SPH) para a resolução de escoamentos de fluidos newtonianos incompressíveis. Diferentemente da maioria das versões existentes deste método, o código numérico faz uso de uma técnica iterativa na determinação do campo de pressões. Este procedimento emprega a forma diferencial de uma equação de estado para um fluido compressível e a equação da continuidade a fim de que a correção da pressão seja determinada. Uma versão paralelizada do simulador numérico foi implementada usando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA) da NVIDIA Corporation. Foram simulados três problemas, o problema unidimensional do escoamento de Couette e os problemas bidimensionais do escoamento no interior de uma Cavidade (Shear Driven Cavity Problem) e da Quebra de Barragem (Dambreak).

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

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.