Moldagem por injeção é um dos mais importantes processos industriais para produção de produtos plásticos finos. Esse processo é dividido essencialmente em quatro estágios: plastificação, preenchimento, empacotamento e resfriamento. O escoamento de um fluido caracterizado por alta viscosidade em uma cavidade estreita é um problema tipicamente encontrado em processos de moldagem por injeção.Neste caso, o escoamento pode ser descrito por uma formulação conhecida como aproximação de Hele-Shaw. Tal formulação pode ser derivada das equações de conservação tridimensionais usando um número de suposições a respeito do polímero injetado e da geometria da cavidade do molde, juntamente com a integração e o acoplamento das equações da conservação da quantidade de movimento e da continuidade. Essa formulação, referindo às limitações da geometria do molde como sendo canais estreitos e quase sem curvatura, é comumente denominada formulação 2 1/2D. Neste trabalho, é apresentada uma técnica para a simulação da fase de preenchimento de um processo de moldagem por injeção, usando essa formulação 2 1/2D, com um método de volumes finitos e malhas não estruturadas. O modelo de Cross modificado com dependência da temperatura de Arrhenius é empregado para descrever a viscosidade do polímero fundido. O campo de distribuição de temperatura é tridimensional e é resolvido usando um esquema semi-Lagrangeano baseado em volumes finitos. As malhas não estruturadas utilizadas são geradas por triangulação de Delaunay e o método numérico implementado usa a estrutura de dados topológica SHE - Singular Handle Edge, que é capaz de lidar com condições de contorno e singularidades, aspectos comumente encontrados em simulações numéricas de escoamento de fluidos. / Injection molding is one of the most important industrial processes for the manufacturing of thin plastic products. This process can be divided into four stages: plastic melting, filling, packing and cooling phases. The flow of a fluid characterized by high viscosity in a narrow gap is a problem typically found in injection molding processes. In this case, the flow can be described by a formulation known as Hele-Shaw approach. Such formulation can be btained from the three-dimensional conservation equation using a number of assumptions regarding the injected polymer and the geometry of the mold, together with the integration and the coupling of the momentum and continuity equations. This approach, referring to limitations of the mould geometry to narrow, weakly curved channels, is usually called 2 1/2D approach. In this work a technique for the simulation of the filling stage of the injection molding process, using this 2 1/2D approach, with a finite volume method and unstructured meshes, is presented. The modified-Cross model with Arrhenius temperature dependence is employed to describe the viscosity of the melt. The temperature field is 3D and it is solved using a semi-Lagrangian scheme based on the finite volume method. The employed unstructured meshes are generated by Delaunay triangulation and the implemented numerical method uses the topological data structure SHE - Singular Handle Edge, capable to deal with boundary conditions and singularities, aspects commonly found in numerical simulation of fluid flow.

A esfera é utilizada como domínio computacional na modelagem de diversos fenômenos físicos, como em previsão numérica do tempo. Sua discretização pode ser feita de diversas formas, sendo comum o uso de malha regulares em latitude/longitude. Recentemente, também para melhor uso de computação paralela, há uma tendência ao uso de malhas mais isotrópicas, dentre as quais a icosaédrica. Apesar de já existirem modelos atmosféricos que usam malhas icosaédricas, não há consenso sobre as metodologias mais adequadas a esse tipo de malha. Nos propusemos, portanto, a estudar em detalhe diversos fatores envolvidos no desenvolvimento de modelos atmosféricos globais usando malhas geodésicas icosaédricas. A discretização usual por volumes finitos para divergente de um campo vetorial utiliza como base o Teorema da Divergência e a regra do ponto médio nas arestas das células computacionais. A distribuição do erro obtida com esse método apresenta uma forte relação com características geométricas da malha. Definimos o conceito geométrico de alinhamento de células computacionais e desenvolvemos uma teoria que serve de base para explicar interferências de malha na discretização usual do divergente. Destacamos os impactos de certas relações de alinhamento das células na ordem da discretização do método. A teoria desenvolvida se aplica a qualquer malha geodésica e também pode ser usada para os operadores rotacional e laplaciano. Investigamos diversos métodos de interpolação na esfera adequados a malhas icosaédricas, e abordamos o problema de interpolação e reconstrução vetorial na esfera em malhas deslocadas. Usamos métodos alternativos de reconstrução vetorial aos usados na literatura, em particular, desenvolvemos um método híbrido de baixo custo e boa precisão. Por fim, utilizamos as técnicas de discretização, interpolação e reconstrução vetorial analisadas em um método semi-lagrangiano para o transporte na esfera em malhas geodésicas icosaédricas. Realizamos experimentos computacionais de transporte, incluindo testes de deformações na distribuição do campo transportado, que mostraram a adequação da metodologia para uso em modelos atmosféricos globais. A plataforma computacional desenvolvida nesta tese, incluindo geração de malhas, interpolações, reconstruções vetoriais e um modelo de transporte, fornece uma base para o futuro desenvolvimento de um modelo atmosférico global em malhas icosaédricas. / Spherical domains are used to model many physical phenomena, as, for instance, global numerical weather prediction. The sphere can be discretized in several ways, as for example a regular latitude-longitude grid. Recently, also motivated by a better use of parallel computers, more isotropic grids have been adopted in atmospheric global circulation models. Among those, the icosahedral grids are promising. Which kind of discretization methods and interpolation schemes are the best to use on those grids are still a research subject. Discretization of the sphere may be done in many ways and, recently, to make better use of computational resources, researchers are adopting more isotropic grids, such as the icosahedral one. In this thesis, we investigate in detail the numerical methodology to be used in atmospheric models on icosahedral grids. The usual finite volume method of discretization of the divergence of a vector field is based on the divergence theorem and makes use of the midpoint rule for integration on the edges of computational cells. The error distribution obtained with this method usually presents a strong correlation with some characteristics of the icosahedral grid. We introduced the concept of cell alignment and developed a theory which explains the grid imprinting patterns observed with the usual divergence discretization. We show how grid alignment impacts in the order of the divergence discretization. The theory developed applies to any geodesic grid and can also be used for other operators such as curl and Laplacian. Several interpolation schemes suitable for icosahedral grids were analysed, including the vector interpolation and reconstruction problem on staggered grids. We considered alternative vector reconstruction methods, in particular, we developed a hybrid low cost and good precision method. Finally, employing the discretizations and interpolations previously analysed, we developed a semi-Lagrangian transport method for geodesic icosahedral grids. Several tests were carried out, including deformational test cases, which demonstrated that the methodology is suitable to use in global atmospheric models. The computational platform developed in this thesis, including mesh generation, interpolation, vector reconstruction and the transport model, provides a basis for future development of global atmospheric models on icosahedral grids.

As complexidades das configurações, da geometria e da topologia dos escoamentos encontrados em sistemas tais como câmaras de combustão de turbo-máquinas requerem métodos numéricos eficientes para modelar os escoamentos turbulentos reativos existentes. Este trabalho de modelagem numérica está relacionado com o desenvolvimento de ferramentas numéricas eficientes necessárias ao estudo dos escoamentos turbulentos reativos nessas aplicações. Na primeira parte desta tese, após relembrar os elementos de base que permitem a descrição de um escoamento turbulento reativo prée-misturado, são apresentadas as diversas etapas que presidem a elaboração do programa computacional. Inicialmente, introduz-se a técnica de compressibilidade artificial, que consiste em modificar as equações governantes acrescentando termos não-estacionários no pseudo-tempo, destinados a introduzir uma velocidade do som finita no domínio de cálculo. Em seguida é exposta a discretização por volumes finitos das equações governantes em malhas não-estruturadas bidimensionais híbridas, que contêm elementos triangulares e/ou quadriláteros. A adição de termos de dissipação artificial bem com o uso de um método do tipo Runge-Kutta para integrar os termos temporais no pseudo-tempo das equações governantes, que completam o método, são comentados em seguir. O programa de cálculo, desenvolvido na linguagem C++, foi validado considerando-se o problema de uma chama plana turbulenta, propagando-se livremente em uma turbulência "congelada". Este fenômeno foi modelado com um termo de produção química médio do tipo flamelet, no qual um corte representando extinção foi introduzido com a finalidade de selecionar uma só velocidade possível de propagação. Os resultados obtidos para uma chama lenta (0,5 m/s) ou rápida (10 m/s) foram comparados com aqueles obtidos na literatura. A precisão dos resultados é bastante satisfatória e os gradientes elevados de pressão estática, característicos do termo fonte empregado, são capturados corretamente, tanto no que se refere a sua intensidade quanto à sua localização, dentro da frente de chama média. Uma análise de sensibilidade e precisão da malha foi igualmente conduzida, assim como um estudo das características de convergência em função do valor do coeficiente de compressibilidade artificial. Tendo sido validado, o programa de cálculo foi, em seguida, utilizado com a finalidade de estudar esteiras inertes ou reativas encontradas à jusante de um obstáculo de seção transversal triangular, colocado no meio de um escoamento turbulento em um canal. Um modelo de turbulência do tipo k-e, acoplado a dois tipos de leis de parede foi utilizado. O desprendimento turbilhonar, típico dos escoamentos de esteiras inertes considerados, foi previsto corretamente com, em particular, uma freqüência de Strouhal de desprendimento com boa concordância com os resultados experimentais disponíveis. As esteiras reativas foram calculadas com um modelo de flamelet, onde os transportes turbulentos da variável de avanço da combustão foram calculados por aproximação do tipo gradiente. Três formas do termo de produção química médio, existentes na literatura, foram utilizados. Aparentemente, a estrutura da esteira reativa prevista é fortemente dependente da formulação adotada para o termo fonte médio. Embora os resultados obtidos sejam comparados àqueles apresentados por outro autores, em uma configuração de escoamento similar, a comparação com resultados experimentais mostra que é necessário uma melhoria na modelagem física de forma a prever este tipo de escoamento reativo de maneira satisfatória.

O presente trabalho consiste no estudo, implementação e análise de métodos não-oscilatórios de alta ordem de precisão em malhas não-estruturadas através da técnica de volumes finitos. A motivação deste trabalho surge da necessidade de se ter simulações acuradas para os escoamentos com altos números de Mach e, portanto, com fortes descontinuidades, que se pretende analisar. Nos últimos anos, o grupo de dinâmica dos fluidos computacional do Instituto de Aeronáutica e Espaço vem desenvolvendo ferramentas computacionais capazes de capturar com precisão as descontinuidades, ondas de choque e superfícies de contato, presentes nos escoamentos aerodinâmicos de interesse. O objetivo principal do presente trabalho de pesquisa consiste em estudar os esquemas essencialmente não-oscilatórios, do inglês essentially non-oscillatory (ENO), e os esquemas essencialmente não-oscilatórios ponderados, do inglês weighted essentially non-oscillatory (WENO). Os esquemas ENO e WENO foram desenvolvidos com o propósito de capturar com acurácia descontinuidades presentes em problemas governados por equações diferenciais parciais hiperbólicas. No desenvolvimento deste trabalho, vai-se utilizar uma ferramenta numérica de simulação de escoamentos compressíveis em duas dimensões já desenvolvida pelo atual grupo de pesquisa. Esta ferramenta resolve as equações governantes da mecânica dos fluidos através da técnica de volumes finitos com uma formulação que atribui os valores das propriedades conservadas aos centróides dos volumes de controle. O código numérico trata de malhas com qualquer geometria poligonal, como triângulos e quadriláteros, além de malhas híbridas compostas por ambos os elementos. As equações de Euler em duas dimensões representam os escoamentos de interesse no presente trabalho. O processo de reconstrução dos esquemas ENO e WENO é descrito em detalhes para polinômios lineares, quadráticos e cúbicos e, portanto, para esquemas numéricos de segunda, terceira e quarta ordem de precisão. Estes polinômios interpolam variáveis primitivas em pontos de quadraturas de Gauss utilizando moléculas de cálculo determinadas por vizinhanças compostas por volumes de controle adjacentes a um volume principal. Métodos upwind são utilizados para se aproximar os valores dos fluxos nas faces dos volumes de controle. As discretizações temporais são efetuadas por esquemas do tipo Runge-Kutta de alta ordem. Uma ferramenta multigrid é utilizada para acelerar a convergência dos esquemas para soluções de estado estacionário e uma técnica de refinamento adaptativo de malha é utilizada para melhorar a qualidade da malha em regiões com altos gradientes das propriedades do escoamento. Aplicações em escoamentos aerodinâmicos de interesse são realizadas e os resultados são comparados com a literatura e com os métodos previamente implementados pelo grupo de pesquisa.

A esfera é utilizada como domínio computacional na modelagem de diversos fenômenos físicos, como em previsão numérica do tempo. Sua discretização pode ser feita de diversas formas, sendo comum o uso de malha regulares em latitude/longitude. Recentemente, também para melhor uso de computação paralela, há uma tendência ao uso de malhas mais isotrópicas, dentre as quais a icosaédrica. Apesar de já existirem modelos atmosféricos que usam malhas icosaédricas, não há consenso sobre as metodologias mais adequadas a esse tipo de malha. Nos propusemos, portanto, a estudar em detalhe diversos fatores envolvidos no desenvolvimento de modelos atmosféricos globais usando malhas geodésicas icosaédricas. A discretização usual por volumes finitos para divergente de um campo vetorial utiliza como base o Teorema da Divergência e a regra do ponto médio nas arestas das células computacionais. A distribuição do erro obtida com esse método apresenta uma forte relação com características geométricas da malha. Definimos o conceito geométrico de alinhamento de células computacionais e desenvolvemos uma teoria que serve de base para explicar interferências de malha na discretização usual do divergente. Destacamos os impactos de certas relações de alinhamento das células na ordem da discretização do método. A teoria desenvolvida se aplica a qualquer malha geodésica e também pode ser usada para os operadores rotacional e laplaciano. Investigamos diversos métodos de interpolação na esfera adequados a malhas icosaédricas, e abordamos o problema de interpolação e reconstrução vetorial na esfera em malhas deslocadas. Usamos métodos alternativos de reconstrução vetorial aos usados na literatura, em particular, desenvolvemos um método híbrido de baixo custo e boa precisão. Por fim, utilizamos as técnicas de discretização, interpolação e reconstrução vetorial analisadas em um método semi-lagrangiano para o transporte na esfera em malhas geodésicas icosaédricas. Realizamos experimentos computacionais de transporte, incluindo testes de deformações na distribuição do campo transportado, que mostraram a adequação da metodologia para uso em modelos atmosféricos globais. A plataforma computacional desenvolvida nesta tese, incluindo geração de malhas, interpolações, reconstruções vetoriais e um modelo de transporte, fornece uma base para o futuro desenvolvimento de um modelo atmosférico global em malhas icosaédricas. / Spherical domains are used to model many physical phenomena, as, for instance, global numerical weather prediction. The sphere can be discretized in several ways, as for example a regular latitude-longitude grid. Recently, also motivated by a better use of parallel computers, more isotropic grids have been adopted in atmospheric global circulation models. Among those, the icosahedral grids are promising. Which kind of discretization methods and interpolation schemes are the best to use on those grids are still a research subject. Discretization of the sphere may be done in many ways and, recently, to make better use of computational resources, researchers are adopting more isotropic grids, such as the icosahedral one. In this thesis, we investigate in detail the numerical methodology to be used in atmospheric models on icosahedral grids. The usual finite volume method of discretization of the divergence of a vector field is based on the divergence theorem and makes use of the midpoint rule for integration on the edges of computational cells. The error distribution obtained with this method usually presents a strong correlation with some characteristics of the icosahedral grid. We introduced the concept of cell alignment and developed a theory which explains the grid imprinting patterns observed with the usual divergence discretization. We show how grid alignment impacts in the order of the divergence discretization. The theory developed applies to any geodesic grid and can also be used for other operators such as curl and Laplacian. Several interpolation schemes suitable for icosahedral grids were analysed, including the vector interpolation and reconstruction problem on staggered grids. We considered alternative vector reconstruction methods, in particular, we developed a hybrid low cost and good precision method. Finally, employing the discretizations and interpolations previously analysed, we developed a semi-Lagrangian transport method for geodesic icosahedral grids. Several tests were carried out, including deformational test cases, which demonstrated that the methodology is suitable to use in global atmospheric models. The computational platform developed in this thesis, including mesh generation, interpolation, vector reconstruction and the transport model, provides a basis for future development of global atmospheric models on icosahedral grids.