Avaliação de métodos numéricos de análise linear de estabilidade para perfis de aço formados a frio.

Para o projeto de estruturas com perfis de aço formados a frio, é fundamental a compreensão dos fenômenos da instabilidade local e global, uma vez que estes apresentam alta esbeltez e baixa rigidez à torção. A determinação do carregamento crítico e a identificação do modo de instabilidade contribuem para o entendimento do comportamento dessas estruturas. Este trabalho avalia três metodologias para a análise linear de estabilidade de perfis de aço formados a frio isolados, com o objetivo de determinar os carregamentos críticos elásticos de bifurcação e os modos de instabilidade associados. Estritamente, analisa-se perfis de seção U enrijecido e Z enrijecido isolados, de diversos comprimentos e diferentes condições de vinculação e carregamento. Determinam-se os carregamentos críticos elásticos de bifurcação e os modos de instabilidade globais e locais por meio de: (i) análise com o Método das Faixas Finitas (MFF), através do uso do programa computacional CUFSM; (ii) análise com elementos finitos de barra baseados na Teoria Generalizada de Vigas (MEF-GBT), via uso do programa GBTUL; e (iii) análise com elementos finitos de casca (MEF-cascas) por meio do uso do programa ABAQUS. Algumas restrições e ressalvas com relação ao uso do MFF são apresentadas, assim como limitações da Teoria Generalizada de Viga e precauções a serem tomadas nos modelos de cascas. Analisa-se também a influência do grau de discretização da seção transversal. No entanto, não é feita avaliação em relação aos procedimentos normativos e tampouco análises não lineares, considerando as imperfeições geométricas iniciais, tensões residuais e o comportamento elastoplástico do material.

Elastografia em imagens de ultrassom utilizando elementos de contorno.

Este trabalho apresenta uma nova metodologia para elastografia virtual em imagens simuladas de ultrassom utilizando métodos numéricos e métodos de visão computacional. O objetivo é estimar o módulo de elasticidade de diferentes tecidos tendo como entrada duas imagens da mesma seção transversal obtidas em instantes de tempo e pressões aplicadas diferentes. Esta metodologia consiste em calcular um campo de deslocamento das imagens com um método de fluxo óptico e aplicar um método iterativo para estimar os módulos de elasticidade (análise inversa) utilizando métodos numéricos. Para o cálculo dos deslocamentos, duas formulações são utilizadas para fluxo óptico: Lucas-Kanade e Brox. A análise inversa é realizada utilizando duas técnicas numéricas distintas: o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC), sendo ambos implementados em Unidades de Processamento Gráfico de uso geral, GpGPUs ( \"General Purpose Graphics Units\" ). Considerando uma quantidade qualquer de materiais a serem determinados, para a implementação do Método dos Elementos de Contorno é empregada a técnica de sub-regiões para acoplar as matrizes de diferentes estruturas identificadas na imagem. O processo de otimização utilizado para determinar as constantes elásticas é realizado de forma semi-analítica utilizando cálculo por variáveis complexas. A metodologia é testada em três etapas distintas, com simulações sem ruído, simulações com adição de ruído branco gaussiano e phantoms matemáticos utilizando rastreamento de ruído speckle. Os resultados das simulações apontam o uso do MEF como mais preciso, porém computacionalmente mais caro, enquanto o MEC apresenta erros toleráveis e maior velocidade no tempo de processamento.

Modelagem do absorvedor e do gerador de ciclos de refrigeração por absorção de calor com o par amônia/água baseados na tecnologia de filme descendente sobre placas inclinadas.

Esse trabalho constitui o desenvolvimento da modelagem térmica e simulação por métodos numéricos de dois componentes fundamentais do ciclo de refrigeração por absorção de calor com o par amônia/água: o absorvedor e o gerador. A função do absorvedor é produzir mistura líquida com alta fração mássica de amônia a partir de mistura líquida com baixa fração mássica de amônia e mistura vapor mediante retirada de calor. A função do gerador é produzir mistura líquido/vapor a partir de mistura líquida mediante o fornecimento de calor. É proposto o uso da tecnologia de filmes descendentes sobre placas inclinadas e o método de diferenças finitas para dividir o comprimento da placa em volumes de controle discretos e realizar os balanços de massa, espécie de amônia e energia juntamente com as equações de transferência de calor e massa para o filme descendente. O objetivo desse trabalho é obter um modelo matemático simplificado para ser utilizado em controle e otimização. Esse modelo foi utilizado para calcular as trocas de calor e massa no absorvedor e gerador para diversas condições a partir de dados operacionais, tais como: dimensões desses componentes, ângulo de inclinação da placa, temperatura de superfície e condições de entrada da fase líquida e vapor. Esses resultados foram utilizados para estabelecer relações de causa e efeito entre as variáveis e parâmetros do problema. Os resultados mostraram que o ângulo de inclinação da placa ótimo tanto para o absorvedor como para o gerador é a posição vertical, ou 90°. A posição vertical proporciona o menor comprimento de equilíbrio (0,85 m para o absorvedor e 1,27 m para o gerador com as condições testadas) e se mostrou estável, pois até 75° não foram verificadas variações no funcionamento do absorvedor e gerador. Dentre as condições testadas para uma placa de 0,5 m verificou-se que as maiores efetividades térmicas no absorvedor e gerador foram respectivamente 0,9 e 0,7 e as maiores efetividades mássicas no absorvedor e gerador foram respectivamente 0,6 e 0,5. É esperado que os dados obtidos sejam utilizados em trabalhos futuros para a construção de um protótipo laboratorial e na validação do modelo.

Models and methods for the direct simulation of rough and micropatterned surfaces

Friction in hydrodynamic bearings are a major source of losses in car engines ([69]). The extreme loading conditions in those bearings lead to contact between the matching surfaces. In such conditions not only the overall geometry of the bearing is relevant, but also the small-scale topography of the surface determines the bearing performance. The possibility of shaping the surface of lubricated bearings down to the micrometer ([57]) opened the question of whether friction can be reduced by mean of micro-textures, with mixed results. This work focuses in the development of efficient numerical methods to solve thin film (lubrication) problems down to the roughness scale of measured surfaces. Due to the high velocities and the convergent-divergent geometries of hydrodynamic bearings, cavitation takes place. To treat cavitation in the lubrication problem the Elrod- Adams model is used, a mass-conserving model which has proven in careful numerical ([12]) and experimental ([119]) tests to be essential to obtain physically meaningful results. Another relevant aspect of the modeling is that the bearing inertial effects are considered, which is necessary to correctly simulate moving textures. As an application, the effects of micro-texturing the moving surface of the bearing were studied. Realistic values are assumed for the physical parameters defining the problems. Extensive fundamental studies were carried out in the hydrodynamic lubrication regime. Mesh-converged simulations considering the topography of real measured surfaces were also run, and the validity of the lubrication approximation was assessed for such rough surfaces.

Estudo fenomenológico e numérico do escoamento estratificado óleo-água ondulado e com mistura na interface

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura.

Calibração de redes de distribuição de água pelo método inverso aplicado a escoamento transiente

Desenvolve-se um método para estimar os parâmetros de uma rede hidráulica a partir de dados observados de cargas hidráulicas transientes. Os parâmetros físicos da rede como fatores de atrito, rugosidades absolutas, diâmetros e a identificação e quantificação de vazamentos são as grandezas desconhecidas. O problema transiente inverso é resolvido utilizando uma abordagem indireta que compara os dados disponíveis de carga hidráulica transiente observados com os calculados através de um método matemático. O Método Transiente Inverso (MTI) com um Algoritmo Genético (AG) emprega o Método das Características (MOC) na solução das equações do movimento para escoamento transiente em redes de tubos. As condições de regime permanente são desconhecidas. Para avaliar a confiabilidade do MTI-AG desenvolvido aqui, uma rede-exemplo é usada para os vários problemas de calibração propostos. O comportamento transiente é imposto por duas manobras distintas de uma válvula de controle localizada em um dos nós da rede. Analisam-se, ainda, o desempenho do método proposto mediante a variabilidade do tamanho do registro transiente e de possíveis erros de leitura nas cargas hidráulicas. Ensaios numéricos realizados mostram que o método é viável e aplicável à solução de problema inverso em redes hidráulicas, sobretudo recorrendo-se a poucos dados observados e ao desconhecimento das condições iniciais de estado permanente. Nos diversos problemas de identificação, as informações transientes obtidas da manobra mais brusca produziu estimações mais eficientes.

Investigação de camadas aqüíferas por métodos geoelétricos no Estado de São Paulo

No Estado de São Paulo, a maioria dos poços profundos perfurados localiza-se em terrenos sedimentares. Assim sendo foi realizado este trabalho a fim de verificar a aplicabilidade de métodos geofísicos à procura de água subterrânea no Estado. A água da chuva atingindo a superfície do terreno toma caminhos a saber: evaporação, escoamento e infiltração. Em face da porosidade, permeabilidade e força da gravidade, a água de infiltração toma o caminho descendente até: - atingir uma zona permeável (saturada ou não) e seu fluxo. - atingir uma rocha sedimentar impermeável ou rocha ígnea sã. - saturar fissuras existentes em rochas ígneas sãs. Basicamente que um poço seja produtor de água, é necessária a extração de água acumulada nas seguintes condições: 1- água acumulada em rochas permeáveis saturadas. 2- água acumulada em fissuras de rochas ígneas sãs. 3- água acumulada em depressões existentes ou na superfície de rochas ígneas ou do embasamento cristalino, coberto por rochas porosas e permeáveis. Baseado nestas formas de ocorrência de água subterrânea, as pesquisas de campo foram orientadas no sentido de: - identificar pela interpretação das curvas ou mesmo uma curva de sondagem elétrica, aquela camada aqüífera, saturada de água e já conhecida \"a priori\" como aqüífero, bem como estabelecer a sua profundidade, espessura e comportamento geológico. - verificar a profundidade do nível freático. - mapear a superfície do embasamento cristalino, procurando achar nela, depressões causa das por erosões pretéritas. Procurou-se verificar também: até que ponto, em um dado local uma estratigrafia geológica conhecida, corresponde a uma \"estratigrafia geofísica\" obtida por medidas na superfície. - determinar a resistividade, ou resistividades, das diferentes formações litológicas. - procurar correlacionar variações laterais de resistividade de uma certa camada geológica com as variações na sua granulometria e litologia ou sais ) dissolvidos na solução que a satura. De um modo geral a eletrorresistividade fornece bons resultados quando aplicada à procura de água subterrânea no Estado de São Paulo. Nos sedimentos do Grupo Tubarão (Permo-Carbonífero) verifica-se uma correspondência entre altos valores de resistividade e as altas vazões dos poços (vazão específica). A \"estratigrafia geofísica\" obtida corresponde muito bem à estratigrafia geológica. Espessuras dos derrames de basalto (Eocretáceo), intrusões de diabásio, arenito da Série Bauru (Neocretáceo), arenito da Formação Botucatu (Eocretáceo), sedimentos do Grupo Tubarão (Permo-Carbonífero) e camadas de sedimentos recentes são facilmente determinadas. Na região da Praia Grande existem lençóis de água doce sobrejazendo à água salobra e salgada.

Estudo numérico do efeito de blocagem hidrodinâmica em um túnel de cavitação sobre o escoamento ao redor de um submarino.

O projeto de submarinos é um tema que vem se desenvolvendo há muitos anos. No Brasil, vivemos um momento de pleno avanço na área, sendo encabeçado pela Marinha do Brasil com o desenvolvimento do primeiro submarino de propulsão nuclear brasileiro. No decorrer desse projeto, diversos estudos devem ser desenvolvidos, dentre eles, ensaios com modelos em escala, que serão realizados em um túnel de cavitação. Durante estes experimentos, espera-se que surjam diversas fontes de incerteza, que podem comprometer os resultados obtidos, destacando-se o efeito de blocagem hidrodinâmica. Nesse contexto, esta dissertação apresenta uma investigação do fenômeno de blocagem no modelo de submarino DARPA SUBOFF na seção de testes de um túnel de cavitação. Para isso, simulações tridimensionais de média de Reynolds (RANS) foram realizadas para o modelo de submarino posicionado no túnel (sob blocagem) e sob condição de águas profundas (referência), para diferentes ângulos de ataque e número de Reynolds (Re) constante. Coeficientes de força e momento hidrodinâmicos foram obtidos e os métodos de correção de blocagem clássico e WICS foram empregados para obter fatores que pudessem corrigir esses valores. A eficiência de cada método foi avaliada comparando os valores corrigidos com os de referência. Os resultados mostraram que ambos métodos de correção são eficientes na recuperação de coeficientes de força e momento hidrodinâmicos para razões de blocagem inferiores a 10%, mas ressalvas são feitas quanto à utilização de cada um desses métodos, sendo o método clássico indicado para o estudo de geometrias simples e o método WICS para modelos mais complexos. Paralelamente, a topologia do escoamento ao redor do submarino foi estudada através do método dos pontos críticos, para avaliar se a blocagem hidrodinâmica pode causar alterações nas estruturas presentes no escoamento. A análise topológica permitiu avaliar a evolução nas estruturas do escoamento para os ângulos de ataque e Re simulados. Foram capturadas as linhas de separação no corpo e avaliada sua relação com a tensão de cisalhamento ao longo do casco. A blocagem hidrodinâmica não causou alterações drásticas na topologia do escoamento.