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.

Simulação numérica da cavitação em turbomáquinas usando uma formulação Euler-Lagrange.

Turbomáquinas são máquinas operacionais que transferem energia mecânica entre um rotor e um fluido. Estas máquinas têm muitas aplicações industriais. Um dos componentes de uma turbomáquina responsável pela transferência da energia, ou receber a rotação do eixo e transformar em energia de fluido em caso de bomba ou transferir a energia do fluido para o eixo em caso de uma turbina, é o impelidor ou rotor. O fenómeno da cavitação envolve escoamento bifásico: o líquido a ser bombeado e as bolhas de vapor que são formadas durante o processo de bombeamento. O processo de formação dessas bolhas é complexo, mas ocorre principalmente devido a presença de regiões de pressões muito baixas. O colapso dessas bolhas pode muitas vezes levar a deterioração do material, dependendo da intensidade ou da velocidade de colapso das bolhas. O principal objetivo deste trabalho foi estudar o comportamento hidrodinâmico do escoamento nos canais do impelidor de uma turbomáquina do tipo radial usando recursos de fluidodinâmica computacional (CFD). Uma abordagem Euler-Lagrange acoplada com o modelo da equação de Langevin foi empregada para estimar a trajetória das bolhas. Resultados das simulações mostram as particularidades de um escoamento líquido-bolha de vapor passando em um canal de geometria curva, fornecendo assim informações que podem nos ajudar na prevenção da cavitação nessas máquinas.

Identificação da fase sigma em uma junta soldada pelo processo TIG manual autógeno de um aço hiperduplex SAF 2707HD (UNS S32707) por difração de raios-x.

O aço inoxidável hiperduplex possui alta resistência a corrosão por pite em ambientes contendo cloretos, quando comparado a outros aços inoxidáveis comercialmente conhecidos. Possui boas propriedades mecânicas, com limite de escoamento superior a 700MPa e limite de resistência a tração em torno de 1000MPa. Essas propriedades o tornam muito atrativos para aplicações em ambientes contendo cloretos, e por isso tem tido destaque na indústria de óleo e gás, refinarias, plataformas offshore, etc. A liga hiperduplex é composta por uma estrutura bifásica, contendo proporções aproximadamente iguais de ferrita e austenita. Esse material possui boa soldabilidade, mas por ser termodinamicamente metaestável, em altas temperaturas pode ocorrer a precipitação de fases intermetálicas não desejáveis, o que resulta em perda de propriedades mecânicas e diminuição da resistência a corrosão. A fase sigma tem sido fortemente estudada, pois é comum sua precipitação nos aços inoxidáveis da família duplex durante o procedimento de soldagem se este não for muito bem controlado. A fase sigma precipita preferencialmente na fase ferrítica, devido a maior concentração de Cr e Mo, que são os elementos formadores da fase. A resistência a corrosão é reduzida e as propriedades mecânicas do material são alteradas o tornando frágil devido a presença da fase sigma. É formada entre 600C e 1000C e possui uma estrutura tetragonal complexa. O objetivo do trabalho foi identificar a possível presença da fase sigma na junta soldada do aço inoxidável hiperduplex SAF 2707 HD (UNS S32707) pelo processo TIG autógeno manual através da difração de raios-x. Nessa pesquisa, foram analisadas uma junta soldada do material pelo processo TIG autógeno manual com arco pulsado. Complementando o estudo foram analisadas seis amostras do aço inoxidável superduplex, sendo que cinco amostras sofreram tratamento térmico para a proposital formação da fase sigma. O refinamento do resultado da difração das amostras foi feito utilizando o método de Rietveld no software Topas Academic versão 4.1. O resultado da amostra soldada de hiperduplex apresentou as fases austenita, ferrita e alguns prováveis óxidos. Os resultados das amostras de superduplex tratadas termicamente apresentaram a fase sigma, conforme esperado na pesquisa, e as fases austenita e ferrita.

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.

Efeitos da atividade física de baixa intensidade na estrutura das túnicas íntima e média da aorta em ratos espontaneamente hipertensos (SHR)

O exercício contínuo de baixa intensidade é capaz de atenuar a hipertrofia da célula muscular lisa da parede da aorta de ratos espontaneamente hipertensos (SHR), e parece atuar sobre a distribuição de fibras oxitalânicas, elaunínicas e elásticas. As alterações funcionais das grandes artérias relacionadas com a idade, deposição de colágeno e elastina na parede arterial e a elasticidade, são também melhoradas com a atividade física. No presente trabalho objetivamos estudar os efeitos da atividade física aeróbica de baixa intensidade no remodelamento estrutural da artéria aorta em modelo de hipertensão genética de ratos SHR. Através da análise da distribuição das fibras oxitalânicas, elaunínicas e elásticas na aorta dos animais controles e SHR submetidos ou não a atividade física de baixa intensidade. Foram utilizados 32 ratos, sendo 16 ratos SHR machos e 16 ratos normotensos Wistar Kyoto (WKY) machos com 8 semanas de idade. Ratos machos foram alocados em 4 grupos: WKY sedentário (WKY-SED), WKY exercitado (EX-WKY), SHR sedentário (SED-SHR), e SHR exercitado (EX-SHR). Os ratos sedentários foram limitados à atividade na caixa, enquanto que os ratos exercitados foram submetidos a um exercício de 1 h / dia, 5 dias / semana. Esses grupos passaram pelo protocolo de atividade física de 20 semanas e a pressão arterial foi mensurada semanalmente (PA). As aortas foram colhidas e processadas para microscopia de luz, microscopia eletrônica e western blotting. Foram realizadas as colorações orcinol neo-fucsina e resorcina-fucsina de Weigert. No grupo hipertenso, o exercício mantém a PA em um nível relativamente semelhante ao inicio do protocolo, mostrando a capacidade de prevenir o aumento da PA ao longo das 20 semanas. No grupo dos animais hipertensos não tratados, a PA aumenta. A PA aumentou progressivamente nos ratos SED-SHRs atingindo 1894 mmHg, mas o exercício físico impediu este processo. Ao final do experimento a PA nos ratos EX-SHRs foi similar o dos ratos WKY (1184 vs. 1144 mmHg), respectivamente. Observou-se maior expressão de elastina e maior distribuição de fibras oxitalânicas, elaunínicas e elásticas em animais que foram submetidos ao protocolo de exercício físico. A porcentagem de fibras elásticas e oxitalânicas foi menor em SED-SHR comparado com SED-WKY, mas o exercício físico aumentou a porcentagem dessas fibras em ambos os grupos. Através da imuno-histoquímica ultra-estrutural para elastina e fibrilina, os grupos EX-WKY e EX-SHR apresentaram uma marcação mais intensa para elastina e fibrilina. Animais hipertensos que não sofreram o protocolo de exercício físico apresentam espessura da parede da aorta maior que a dos animais que sofreram o exercício. O número de lamelas elásticas, bem como as fibras oxitalânicas e elaunínicas, é maior no grupo EX-SHR, em relação ao SED-SHR. Os grupos exercitados tiveram maior expressão de eNOS que seus respectivos grupos sedentários, e as células endoteliais apresentaram características morfológicas preservadas. A associação da atividade física com modelos de hipertensão genética mostra que o exercício físico tem efeitos benéficos nessa situação, uma vez que atenua a hipertensão e o remodelamento adverso da parede da aorta.

O direito prestacional à saúde e sua implementação judicial limites e possibilidades

A presente dissertação centra-se no estudo do direito prestacional à saúde e de sua sindicação por intermédio de demandas individuais e coletivas. Para tanto, demonstrada a fundamentalidade formal e material do direito à saúde, analisa-se como os tribunais brasileiros se comportam diante das muitas objeções opostas à sua atuação. Constatada a elasticidade do tratamento conferido pela jurisprudência às demandas individuais, propõem-se critérios objetivos limitadores; diante da contrastante timidez dos julgados na seara coletiva, o foco passa a ser o escrutínio das possibilidades de controle judicial inexploradas.

Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

Uma nova descrição para a transferência de massa em meios porosos com transição saturado-insaturado.

Esse texto trata do problema de um fluido contaminado escoando por um meio poroso, tratando os componentes na mistura como meios contínuos. Na primeira parte, desenvolvemos a teoria de misturas de meios contínuos e discutimos equações da continuidade, momento linear e momento angular. A seguir, descrevemos o problema em detalhe e fazemos hipóteses para simplificar o escoamento. Aplicamos as equações encontradas anteriormente para encontrarmos um sistema de equações diferenciais parciais. Desse ponto em diante, o problema se torna quase puramente matemático. Discutimos o caso insaturado, e depois a saturação do meio poroso. Finalmente, adicionamos um contaminante à mistura e, em seguida, N contaminantes.

Desenvolvimento e implementação de algoritmo para análise de dutos fabricados ou reparados com material compósito carregados com pressão interna utilizando o método dos elementos finitos.

A integridade de dutos ganha importância à medida em que o desenvolvimento da indústria expande a malha para transportar líquidos e gases num contexto de aumento das exigências ambientais. Um importante aliado para manutenção da integridade de dutos são reparos de materiais compósitos. Estes materiais apresentam baixa densidade, capacidade de direcionar resistência de acordo com as direções das solicitações, execução de reparo a frio sem necessidade de interromper produção ou grande maquinário. Este trabalho tem como objetivo desenvolver e implementar um algoritmo de elementos finitos que permita avaliar os esforços e a resistência das paredes de um tubos fabricados ou reparados com laminados de material compósito carregados com pressão interna. Entre as vantagens de desenvolver um programa tem-se: agilidade de avaliação, menor custo com licença, menores exigências computacionais, possibilidade de desenvolver o programa e o melhor entendimento da modelagem dos fenômenos. Utiliza-se como entrada do programa o diâmetro do duto, pressão interna e parâmetros do laminado. A modelagem em elementos finitos é realizada a partir da teoria clássica de laminados. Aplicando o carregamento resultante da pressão interna, determina-se os deslocamentos e são calculadas as tensões e aplicado o critério de falha de Tsai-Hill em cada camada. Estudos experimentais e numéricos encontrados na literatura foram simulados com o programa gerado e os resultados para propriedades do laminado, tensões nos dutos e pressão de ruptura apresentam concordância com os resultados da literatura.O programa ainda tem sua estrutura modificada para encontrar a pressão de falha a partir dos dados do laminado. O programa implementado permite uma avaliação rápida de resistência do reparo e possibilita avaliar rapidamente a resposta a mudanças nos parâmetros de projeto do laminado.

Simulação de problemas de transferência de calor em regime permanente com uma relação entre condutividade térmica e temperatura constante por partes.

Este trabalho estuda a transferência de calor por condução considerando a condutividade térmica como uma função constante por partes da temperatura. Esta relação, embora fisicamente mais realista que supor a condutividade térmica constante, permite obter uma forma explícita bem simples para a inversa da Transformada de Kirchhoff (empregada para tratar a não linearidade do problema). Como exemplo, apresenta-se uma solução exata para um problema com simetria esférica. Em seguida, propôe-se uma formulação variacional (com unicidade demonstrada) que introduz um funcional cuja minimização é equivalente à solução do problema na forma forte. Finalmente compara-se uma solução exata obtida pela inversa da Transformada de Kirchhoff com a solução obtida via formulação variacional.

Comportamento mecânico de cerâmica produzida a partir de reciclagem de para-brisa de automóvel

É crescente o interesse nos materiais cerâmicos, devido as suas características como baixa massa específica e maior resistência a ambientes agressivos do que a maioria das ligas metálicas. Este trabalho tem o objetivo de produzir a partir de material descartado, cerâmicas em diferentes temperaturas de sinterização e avaliar a sua tensão de ruptura em ensaio de flexão de três pontos e a confiança desta medida. Devido ao custo de produzir pó de vidro, tanto pelo alto gasto de energia para fundir a matéria-prima como pelo consumo de minerais industriais, foi proposto utilizar vidro de para-brisa obsoleto reduzindo despesas operacionais e definindo um destino econômico e ambiental viável para estes rejeitos. A metodologia consistiu-se na obtenção do pó de vidro com características adequadas para ser conformado e sinterizado. Foram usadas duas composições e quatro tratamentos térmicos para obter oito materiais. Uma composição com apenas o pó oriundo da moagem de para-brisa e outra com este pó mais 4% de óxido de nióbio. A resistência à flexão dos produtos obtidos foi avaliada. Utilizou-se a estatística de Weibull para caracterizar estes resultados. Os resultados obtidos indicam que o material de composição pó de vidro e temperatura final de sinterização de 650C obteve a maior resistência mecânica entre os materiais sintetizados. A adição do óxido de nióbio provoca um decréscimo na resistência mecânica se comparada com o material sem a adição deste óxido. Entretanto, comparando as duas composições na mesma temperatura final de sinterização, a adição de óxido de nióbio provocou um aumento no módulo de Weibull, excetuando-se dois de oito materiais obtidos. As diferentes composições e temperaturas de sinterização afetaram as propriedades mecânicas dos materiais obtidos.

Avaliação de propriedades mecânicas e térmicas de compósito à base de polietileno de alta densidade e hidroxiapatita deficiente de cálcio

No presente trabalho, foram processados compósitos de polietileno de alta densidade (PEAD) com hidroxiapatita deficiente de cálcio (HA), com o objetivo de obter materiais com melhores propriedades mecânicas e bioatividade. A adição da HA deficiente de cálcio proporcionou um aumento no módulo de elasticidade (maior rigidez), menor resistência ao impacto e decréscimo do grau de cristalinidade do PEAD, proporcionando uma maior bioatividade ao material. A análise térmica exploratória (sistema não isotérmico) foi realizada por meio da técnica de calorimetria exploratória diferencial (DSC) e foram avaliados os teores de fosfato de cálcio e a velocidade de rotação da rosca no processamento dos materiais. No estudo da cristalização não-isotérmica observou-se uma diminuição da temperatura de cristalização com o aumento da taxa de resfriamento para todos os materiais sintetizados. A energia de ativação (Ea) da cristalização dos materiais foi avaliada por meio dos métodos Kissinger e Ozawa. A amostra com 5% de HA deficiente de cálcio e velocidade de processamento de 200 rpm foi a que apresentou menor valor de energia de ativação, 262 kJ/mol, menor desvio da linearidade e a que mais se assemelhou à matriz de PEAD sem HA. O teor de hidroxiapatita deficiente de cálcio não favorece o processo de cristalização devido à alta energia de ativação determinada pelos métodos descritos. Provavelmente, a velocidade de rotação, favorece a dispersão da carga na matriz de PEAD, dificultando o processo de cristalização. Na aplicação do método de Osawa-Avrami, os coeficientes de correlação indicaram perda na correlação linear. Estas perdas podem estar associadas a uma pequena percentagem de cristalização secundária e/ou à escolha das temperaturas utilizadas para determinar a velocidade de cristalização. Na determinação dos parâmetros pelo método de Mo, as menores percentagens de cristalização apresentaram um grande desvio da linearidade, com coeficiente de correlação bem menor que 1 e com o aumento da percentagem de cristalização, o desvio da linearidade diminui, ficando próximo de 1. Os resultados obtidos mostraram que o modelo de Mo e de Osawa-Avrami não foram capazes de definir o comportamento cinético dos materiais produzidos neste trabalho.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Métodos sem malha: aplicações do Método de Galerkin sem elementos e do Método de Interpolação de Ponto em casos estruturais.

Apesar de serem intensamente estudados em muitos países que caminham na vanguarda do conhecimento, os métodos sem malha ainda são pouco explorados pelas universidades brasileiras. De modo a gerar uma maior difusão ou, para a maioria, fazer sua introdução, esta dissertação objetiva efetuar o entendimento dos métodos sem malha baseando-se em aplicações atinentes à mecânica dos sólidos. Para tanto, são apresentados os conceitos primários dos métodos sem malha e o seu desenvolvimento histórico desde sua origem no método smooth particle hydrodynamic até o método da partição da unidade, sua forma mais abrangente. Dentro deste contexto, foi investigada detalhadamente a forma mais tradicional dos métodos sem malha: o método de Galerkin sem elementos, e também um método diferenciado: o método de interpolação de ponto. Assim, por meio de aplicações em análises de barras e chapas em estado plano de tensão, são apresentadas as características, virtudes e deficiências desses métodos em comparação aos métodos tradicionais, como o método dos elementos finitos. É realizado ainda um estudo em uma importante área de aplicação dos métodos sem malha, a mecânica da fratura, buscando compreender como é efetuada a representação computacional da trinca, com especialidade, por meio dos critérios de visibilidade e de difração. Utilizando-se esses critérios e os conceitos da mecânica da fratura, é calculado o fator de intensidade de tensão através do conceito da integral J.

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.