Este trabalho apresenta métodos para simulação de escoamentos a altos números de Mach utilizando diversos esquemas de discretização espacial. Dentre estes esquemas, considerou-se um esquema centrado e os esquemas de separação de vetores de fluxo de va Leer e de Liou. Os esquemas "upwind"foram implementados tanto em suas versões de primeira ordem quanto de segunda ordem de precisão espacial. As versões de segunda ordem de precisão estão baseadas numa abordagem "MUSCL", onde os estados à direita e à esquerda são reconstruídos linearmente e afetados por um limitador. O trabalho também discute uma técnica de geração de malhas não estruturadas bidimensionais, baseada no método de avanço de frente de geração. A malha é gerada através da introdução automática de pontos de domínio computacional. Utilizam-se termos fonte locais, que auxiliam na concentração de pontos em regiões específicas de domínio. Foram implementadas algumas técnicas de melhoramento de malha, tais como troca de diagonais e suavização. A metodologia de simulação de escoamento implementada utiliza ainda um procedimento de refinamento adaptativo, sendo que o processo de refinamento de triângulos usa um sensor baseado em propriedades físicas do escoamento. Os resultados foram obtidos através da simulação do escoamento em uma configuração típica de entrada de ar supersônica/hipersônica. Os diferentes esquemas de discretização espacial utilizados neste trabalho foram capazes de reproduzir adequadamente as características do escoamento.

A duplicação de código é um mau que deve ser combatido a qualquer custo. Existem situações em que parece não haver alternativa, o ambiente de desenvolvimento parece requerer que a informação seja replicada, ou então simplesmente não se percebe que ela está sendo duplicada e isso pode afetar negativamente o desenvolvimento de software de diversas maneiras. A duplicação de código resulta no crescimento do código-fonte e de sua complexidade, prejudicando a manutenção do sistema. Mesmo se fosse possível encontrar e alterar todas as réplicas de um trecho de código, é impossível garantir que todas as mudanças foram realizadas consistentemente - que as regiões comuns permanecem idênticas e que as diferenças foram mantidas - sem comparar manualmente cada trecho de código replicado, palavra por palavra, na esperança de que nenhum detalhe tenha sido deixado de lado. Independente da razão, há meios de se evitar a replicação de código. A geração de código-fonte é uma ferramenta conhecida para evitar que um trecho de código seja espalhado para diversas regiões da aplicação. Em diversas linguagens, estruturas podem ser criadas a partir de uma fonte comum de metadados utilizando um gerador de código-fonte sempre que o código é compilado. Particularmente, se tratando de interfaces com o usuário, outra abordagem seria utilizar uma fonte de metadados para gerar código-fonte em tempo de execução. Esse trabalho apresenta o MAGIU (Modelo Arquitetural de Geração de Interfaces com o Usuário), um modelo conceitual para criação de geradores de código-fonte para o desenvolvimento de interfaces com o usuário.

Sabe-se que um bom projeto de coletor de admissão é um fator importante para um bom desempenho de um motor de combustão interna. Modelos computacionais têm sido usados para simular escoamentos não-permanentes e compressíveis fornecendo assim suporte para melhor projeto do coletor. Neste trabalho, o estudo de escoamento compressível, intermitente e turbulento de fluido newtoniano no interior de cilindros de seção circular tem por objetivo entender os fenômenos que ocorrem dentro de um sistema complexo de dutos que conduzem o ar até os cilindros de um motor. Para tal é necessário primeiramente conhecer quais são estes fenômenos e identificar os parâmetros cruciais que influenciam no processo de solução do modelo matemático em uma geometria simplificada e assim possibilitar a avaliação do sistema em operação em uma geometria complexa. O modelo matemático usado para o escoamento compressível pulsátil dentro do duto circular reto é o URANS (Unsteady Reynolds Averaged Navier-Stokes). Para resolver as equações governantes do modelo matemático (continuidade, quantidade de movimento, energia e modelo de turbulência) é usado o software ANSYS Fluent CFD. Esta ferramenta é baseada na discretização de volumes finitos. Este trabalho é composto por três etapas. Na primeira etapa a abordagem baseada na pressão é usada resolvendo as equações governantes de forma segregada. Após um estudo de malha, os resultados obtidos mostraram que um duto excitado por um único pulso tem uma resposta periódica com frequência constante, e certo grau de amortecimento ao longo do comprimento do duto. Foi mostrado que o comprimento do duto é o parâmetro de maior influencia no comportamento dinâmico do sistema. Na segunda etapa a mesma geometria foi empregada e o mesmo tipo de abordagem foi usada. Foi observada a resposta dinâmica de uma condição de contorno de pressão oscilante na saída do duto. É possível verificar a importância da compressibilidade do ar que atenua as amplitudes das oscilações de pressão e vazão em massa ao longo do duto. E finalmente, a terceira etapa consistiu em uma geometria complexa, formada por um duto principal, um plenum e dois runners. Como condições de contorno na saída dos dois dutos secundários foram impostas perturbações de pressão de mesma amplitude e fases variáveis, e assim verificada a vazão em massa para cada defasagem pré-estabelecida.

Neste trabalho estudamos analiticamente e numericamente a topologia do campo magnético em tokamaks, nos focando, basicamente, no tokamak ITER (International Thermonuclear Experimental Reactor) em construção na França. Para as análises numéricas utilizamos como ferramenta o conjunto integrado de códigos CRONOS, que realiza simulações de descargas em tokamaks, e acopla uma série de códigos computacionais, tais como o módulo ';Helena';, usado para o estudo do equilíbrio MHD, entre outros. Já a parte analítica foi realizada pela simulação de superfícies magnéticas através de equações diferenciais ordinárias. Simulações numéricas e analíticas são apresentadas para densidades de corrente toroidal com buraco (oca), relatando também o aparecimento de cadeias de ilhas próximas ao eixo magnético, causadas pela inversão de corrente. Simulações das superfícies magnéticas do tokamak ITER, através de campos magnéticos gerados por fios, também são analisadas, assim como alterações na topologia das superfícies causadas pela adição de um ruído que simula colisões de partículas dentro da coluna de plasma. Identificamos neste último modelo estruturas que aprisionam linhas de campo magnético por muitas voltas toroidais, alterando os padrões de escape para as placas do divertor através do efeito stickiness.

Este trabalho apresenta um estudo comparativo entre três ferramentas computacionais diferentes na determinação do coeficiente de pressão em um caso típico de aerofólio supercrítico, em particular o Boeing A4, com sua conseqüente comparação com os dados experimentais obtidos em túnel de vento para este aerofólio. Esta comparação contribui para o processo de validação de uma ferramenta de mecânica dos fluidos computacional (CFD) que vem sendo desenvolvida no presente grupo de trabalho. Os códigos de CFD empregados resolvem uma formulação de Navier-Stokes com média de Reynolds (RANS) ou uma formulação de Euler acoplada a soluções de camada limite. Neste último caso, tem-se a opção de fixar a transição sobre a superfície do aerofólio ou deixá-la livre. Todos os códigos utilizam modelos de turbulência do tipo viscosidade de vórtice para levar em consideração os efeitos de transporte de turbulência. Em particular, o modelo de uma equação de Spalart-Allmaras e o modelo de duas equações k-e RT são empregados nos códigos com formulação RANS. O código de camada limite utiliza o modelo algébrico de Cebeci-Smith. O trabalho apresenta os aspectos da integração no tempo e a integração espacial, bem como os procedimentos de aceleração de convergência usados no código que está sendo desenvolvido, o que possibilita numa análise posterior estabelecer relações de tendências de resultados de acordo com os modelamentos e equações que estão sendo utilizados. Foram utilizados onze casos teste diferentes, todos utilizando a configuração do aerofólio Boeing A4. Os resultados obtidos permitiram analisar a qualidade e a confiabilidade dos resultados dos códigos no que se diz respeito à correta identificação de possíveis ondas de choque, valores de coeficientes de pressão, posicionamento da transição de escoamento laminar para turbulento e tendências que os códigos possam apresentar de oscilações numéricas na solução. Estes resultados e análises realizadas dão subsídios importantes no processo de validação do código em desenvolvimento.

A formação de gelo em aeronaves é um tópico de grande interesse para a indústria aeronáutica. Vários acidentes fatais já ocorreram na aviação comercial em decorrência da formação de gelo em partes sensíveis de aeronaves. Por conta disso, os órgãos certificadores elaboraram uma série de requisitos em relação aos mecanismos de proteção e de eliminação de gelo, assim como na avaliação da qualidade de vôo e desempenho de uma aeronave que opere em condições de gelo acumulado nas suas superfícies de sustentação e de controle. Ao longo dos anos, muitos estudos têm sido conduzidos de maneira a se alcançar uma rápida e precisa maneira de avaliar esse fenômeno e com isso reduzir o ciclo de projeto. Nesse trabalho, é apresentado um método simples de avaliação do mecanismo de acúmulo de gelo em aerofólios, baseado principalmente no cálculo da eficiência de coleta, um dos principais parâmetros envolvidos no processo. Partindo-se de um código numérico desenvolvido por Lozowski e Oleskiw [1], foi realizada uma reestruturação do mesmo. Algumas subrotinas foram atualizadas, assim como novos e melhores procedimentos de cálculo foram adicionados ao código original. A versão atual do código foi validada através da comparação com a versão original e, posteriormente, foram conduzidas análises da influência de diversos parâmetros na eficiência de coleta resultante de um determinado corpo submetido a condições conhecidas. Além disso, foram avaliados alguns casos relativos à formação de gelo e os resultados são comentados aqui.

Simulações de grandes escalas (LES) atualmente têm desempenhado um papel importante no estudo de turbulência, aeroacústica e combustão. Em LES, as maiores escalas são resolvidas e as escalas menores são modeladas por um modelo de submalha. A grande vantagem deste tipo de simulação para o estudo de turbulência está no fato de que a produção de turbulência acontece primordialmente nas grandes escalas. Além disso, o comportamento das grandes escalas é extremamente afetado pela topologia do escoamento analisado. Neste contexto, o objetivo primordial deste trabalho é estudar escoamentos próximos a paredes, em regiões com curvatura, para avaliar o efeito de gradientes de pressão adversos nas estatísticas do escoamento, analisando as consequências destes gradientes de pressão na produção, transporte e dissipação de energia turbulenta. A formulação empregada nas simulações é uma formulação LES, e estão sendo analisados números de Reynolds moderados, para os quais, em geral, não se têm dados de simulações diretas na literatura. Desta forma, pretende-se entender um pouco mais sobre as estruturas turbulentas geradas nas regiões de gradiente de pressão adverso e os mecanismos de produção de turbulência em tais regiões.

Nas últimas décadas, as melhorias na tecnologia de extit{hardware} principalmente no que diz respeito a unidade de processamento central, memória e capacidade de armazenamento, diversas ferramentas numéricas aplicadas na engenharia foram desenvolvidas com o intuito de se obter respostas rápidas e mais baratas do que ensaios, experimentos e testes laboratoriais. A Dinâmica dos Fluidos Computacional (DFC) faz uma aproximação numérica como um meio de simular escoamentos. Um dos grandes desafios em soluções numéricas na DFC é a garantia da consistência, estabilidade e convergência garantindo com isso, a confiabilidade do código numérico. As turbomáquinas operam em uma ampla faixa de velocidades no campo do escoamento, podendo ter regiões com baixos números de Mach. Nesse ponto, surge um problema, pois códigos da DFC para o cálculo de escoamento compressível geralmente são programados seguindo formulações matemáticas baseadas no vetor de variáveis conservadas, podendo ter problemas numéricos e apresentar soluções inconsistentes devido a pequena variação da densidade do fluido. Chorin através de alguns estudos constatou que os códigos de DFC possuem uma diminuição de sua convergência e estabilidade quando o número de Mach do escoamento torna-se menor do que o limite de compressibilidade. Com o objetivo de solucionar este problema, Chorin criou a técnica de compressibilidade artificial. Choi e Merkle; Turkel e outros pesquisadores modificaram a técnica de compressibilidade artificial e desenvolveram a técnica de pré-condicionamento para diversas aplicações para garantir a convergência, estabilidade e melhorias na soluções numéricas dos códigos DFC para uma ampla faixa de velocidades do escoamento. No presente trabalho, foi realizada a implementação do método de pré-condicionamento de Turkel, Vatsa e Radaspiel em um programa de DFC desenvolvido por Tomita para calcular escoamentos compressíveis, internos e externos. A fim de melhorar a robustez do código numérico o pré-condicionador foi modificado seguindo o proposto no trabalho de Damorfal e Siu. O pré-condicionamento foi implementado primeiramente para as Equações de Euler e depois para escoamentos viscosos e turbulentos. Diversas experiências numéricas foram realizadas e resultados interessantes utilizando o código de DFC compressível e o código de DFC compressível com o pré-condicionador original e modificado foram obtidos. Todos estes casos estão apresentados e comentados no presente trabalho. A implementação do pré-condicionador em um código de DFC desenvolvido para o cálculo de escoamento tipicamente compressível resultou em uma ferramenta numérica com maior potencial para o cálculo do escoamento numa ampla faixa de número de Mach.

Uma das principais preocupações consideradas na Engenharia de Software consiste na redução do custo de manutenção durante seu ciclo de vida, que pode compreender até 67% de seu custo total. A qualidade interna do código-fonte de um software representa um dos aspectos relacionados ao seu grau de manutenibilidade. Características de seu design como sua estruturação e facilidade de compreensão podem influenciar o esforço associado à sua manutenção. Especialistas em projetos de sistemas Orientados a Objetos (OO) defendem que a manutenção de baixo acoplamento entre os módulos resultam em produtos de software mais confiáveis e manuteníveis. Ao longo dos anos, autores formularam heurísticas para auxiliar os programadores a realizarem a tarefa de controlar o acoplamento entre módulos de um software. Uma destas heurísticas, a Lei de Demeter (Law of Demeter - LoD), foi proposta com a motivação de favorecer a modularização de um software através de sua aplicação. A Análise Estática de Código (AES) consiste em uma dentre as possíveis abordagens para automação da verificação da aderência ou violação da LoD no código-fonte de um software. Entretanto, sua aplicaçãoo pode apre- sentar alguns desafios como longos tempos de execução de análise, ruído excessivo e taxas inaceitáveis de falsos-positivos. Neste trabalho de pesquisa, se investiga, concebe, implementa e verifica um framework, denominado Law Of Demeter Reduced (LoDRed), para identificação automatizada de violações da LoD utilizando AES e redução do volume de ruídos resultantes. Nele, propõe-se ainda uma extensão para LoD, a Extensão da Lei de Demeter (Law of Demeter Extension - LoD-E) com a finalidade de investigar, definir e suportar, sob a ótica da LoD, elementos para análise de algumas construções particulares de linguagens OO como generics e lambda expressions, além de considerar a aplicação de filtros sobre agrupamentos lógicos. O framework LoDRed consiste na principal contribuição desta pesquisa. Ele define um processo, composto por quatro métodos que, em conjunto, têm a finalidade de identificar e reportar, de forma automatizada, violações da LoD e da LoD-E nos códigos-fonte de softwares. A sua implementação resultou em três produtos de software de código aberto e foi verificada nesta pesquisa por meio da execução de um experimento que apresentou redução no volume de ruídos da ordem de até 95%.

A previsão acurada de fenômenos aerodinâmicos não estacionários tem cada vez mais dependido do emprego de técnicas de Mecânica dos Fluidos Computacional (CFD). A aplicação de tais técnicas em problemas não estacionários depende crucialmente de tecnologias de geração de malha que, usualmente, são processos bastante lentos. Para contornar esta limitação, várias metodologias têm sido desenvolvidas a fim de calcular a dinâmica geométrica do problema físico sem gerar novamente uma malha computacional por completo. Neste sentido, este trabalho apresenta um estudo de técnicas eficientes de geração e movimentação de malhas computacionais para simulação não estacionária em CFD, com aplicações em Aeroelasticidade. Para este propósito, é implementada uma abordagem de movimentação de malhas baseada em Funções de Base Radial (RBF). Nesta metodologia, os deslocamentos dos pontos de malha na superfície da estrutura (móvel) são propagados para os pontos volumétricos (interiores), sem alterar a topologia da malha. Este trabalho conduz uma análise comparativa do emprego de diversos tipos de RBF na movimentação de malhas de níveis de refinamento diferentes, bidimensionais e tridimensionais, no que diz respeito à qualidade da malha deformada resultante. Além disso, aborda estratégias de movimentação e sua influência na qualidade de malha gerada. Como aplicação, é conduzida simulação aerodinâmica computacional com código BREXP3D. Este código utiliza formulação de mecânica dos fluidos não estacionária compressível e não viscosa (equações de Euler), implementada numericamente via método de diferenças finitas, o que exige malhas hexaédricas estruturadas. Este trabalho implementa condições de contorno para simulação de asa tridimensional em malhas monobloco e conduz simulação aerodinâmica estacionária em regime transônico, comparando com resultados experimentais. Problemas de malha e convergência não permitiram executar com sucesso a simulação aerodinâmica computacional não estacionária.