Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.

Resistência ao cisalhamento de solos não saturados : análise experimental e teórica

Um equipamento triaxial convencional originalmente utilizado para ensaios em amostras de solos saturados saturada foi modificado para que fosse possível a realização de ensaios triaxiais em solos não saturados. As principais modificações foram a instalação de uma pedra cerâmica de alto valor de entrada de ar (300 kPa) na nova base da câmara triaxial e uma nova linha de aplicação de pressão para aplicação de pressão de ar na amostra. A técnica de translação de eixos é utilizada para evitar cavitação no sistema de aplicação de pressão. O solo estudado nesta pesquisa consiste de um colúvio de arenito localizado nas encostas da Formação Serra Geral, entre os municípios de Timbé do Sul (SC) e São José dos Ausentes (RS). Ensaios triaxiais convencionais foram realizados nas condições drenado e não drenado em amostras indeformadas para tensões de confinamento iniciais variando de 50 até 500 kPa; ensaios triaxiais convencionais foram realizados na condição drenada em amostras remoldadas para as tensões de confinamento iniciais de 50, 100 e 200kPa. Ensaios triaxiais com sucção controlada foram realizados na condição drenado em amostras indeformadas para as tensões normais líquidas de 50, 100, 150 e 200 kPa em níveis de sucção variando de 25 até 150 kPa. A partir dos resultados dos ensaios triaxiais determinou-se os parâmetros de resistência ao cisalhamento: c’ – intercepto coesivo efetivo; f´ - ângulo de atrito interno efetivo; fb – ângulo de atrito interno que quantifica a contribuição da sucção na resistência ao cisalhamento. O comportamento encontrado foi de uma envoltória de resistência ao cisalhamento não saturada bilinear com acréscimo na resistência ao cisalhamento para baixos níveis de sucção até aproximadamente 75 kPa, e após um decréscimo na resistência ao cisalhamento. A tensão normal líquida não influencia no valor de fb e a sucção não provoca alteração no valor de f´. O valor de fb apresenta valores superiores ao valor de f´, quando se utiliza o critério de ruptura clássico de resistência ao cisalhamento. Alguns critérios de determinação dos valores máximos de ruptura foram aplicados para um melhor entendimento do comportamento resistente deste material.