Spelling suggestions: "subject:"volumes finitos."" "subject:"volumes tinitos.""
31 |
Simulação numérica dos campos de temperatura e velocidade em armazenadores térmicosBandini, Marcos Augusto January 1997 (has links)
No presente trabalho avalia-se o comportamento da difusão e da convecção natural em sistemas térmicos que possuem armazenadores de calor, como os que envolvem coletores solares ligados a reservatórios de água quente. Através da solução das equações diferenciais que regem o problema, obtém-se os campos de velocidade e temperatura no interior destes armazenadores. Este último, especialmente, é importante para a determinação da performance térmica de todo o sistema. A intensidade da estratificação térmica depende do fluxo de entrada e saída de água no reservatório, da sua razão de aspecto, da posição da resistência elétrica de apoio, se ela existir, da espessura e difusividade térmica da parede e da espessura e material do isolamento térmico. A perda de calor para o meio ambiente causa um fluxo de líquido descendente, junto às paredes do reservatório, incorrendo em camadas limites hidrodinânlica e térmica. Além do papel destas camadas limites, também investiga-se a importância relativa da difusão de calor através da parede metálica do armazenador, na sua direção axial (altura). Não são investigados neste trabalho os efeitos de entrada e saída do líquido do armazenador, e também as correntes convectivas originadas do suprimento de calor por resistência elétrica. As equações da energia e do movimento são resolvidas de forma acoplada, em coordenadas cilíndricas axissimétricas, em regime transiente, empregando o Método dos Volumes Finitos. / In the present work the diffusion behavior as well as the natural convection in thermal systems that possess heat reservoirs, like the ones that involve solar collectors connected to hot water reservoirs, have been evaluated. Through the solution of the differential equations that describe the problem, the velocity fields and the internal temperature of these tanks have been obtained. The last one is specially important for the determination of the thermal performance of the whole system. Intensity of the thermal stratification depends on the entrance and exit water flow in the reservoir, on the tanks aspect ratio, on the position support of the electrical resistance, if it exists, on the thickness and thermal diffusion of the wall and on the material of thermal insulation. The loss of heat for the environment brings about a descending liquid flow, next to the walls of the reservoir, ending up by creating hidrodynamic and thermal boundary layers. Besides the role of these boundary layers, the relative importance o f the heat diffusion through the metallic wall of the reservoir in its axial direction (height), also has been investigated. The effects of the entrance and exit of liquid in the reservoir, as well as the convective currents originated from the heat supplied by electric resistance have not been investigated. The energy and moment equations are solved in a coupled way, in axissimetrical cilyndrical coordinates, in transient regime, employing a Finite Volume Method.
|
32 |
Simulação numerica do escoamento livre em um canal helicoidal de seção retangularMelgarejo Morales, Rigoberto Eleazar 27 October 2000 (has links)
Orientador: Eugenio Spano Rosa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-28T23:56:11Z (GMT). No. of bitstreams: 1
MelgarejoMorales_RigobertoEleazar_D.pdf: 8900310 bytes, checksum: 9f57b0d7b80e8430b0d0c8941dda51e1 (MD5)
Previous issue date: 2000 / Resumo: Neste trabalho é realizado uma simulação numérica do escoamento com superficie livre completamente desenvolvido em um canal helicoidal de seção retangular, em regime laminar e turbulento. As equações de conservação foram escritas para um sistema de coordenadas ortogonais locais. O sistema de equações foi discretizado e resolvido pelo método de volumes fmitos com o esquema ln'brido de interpolação, implementados no PHOENICS-CFD. A interface livre é calculada utilizando o algoritmo HOL (Height-Of-Liquid). As simulações numéricas foram realizadas para fluidos com viscosidade de 1 a 100cP e Reynolds de 352 a 45500. A turbulência é modelada com o modelo de duas equações 1(-E para baixo números de Reynolds com uma correção devido ao efeito da curvatura. As simulações mostram que o método utilizado, quando comparado com dados experimentais, calcula satisfatoriamente a localização da interface. Apresenta-se os campos da velocidade média, do escoamento secundário, da energia cinética turbulenta e da dissipação e produção da energia cinética turbulenta. A partir dos resultados numéricos calcula-se o fator de atrito local e médio / Abstract: A numerical simulation of the fully developed flow in an open helical channel with rectangular cross section is accomplished, for the laminar and turbulent regime. The conservation equations are written in a local orthogonal coordinate system. The conservation equations are solved by the finite volume method with the hybrid scheme of interpolation, implemented in PHOENICS-CFD. The ftee surface is calculated using the HOL algorithm (Height-Of-Liquid). The numerical simulations were accomplished for fluids with viscosity varying ftom 1 to 100cP and Reynolds ftom 352 to 45500. The turbulence is modeled by two equations k-E model for low Reynolds numbers with a correction due to the effect of the curvature. The numerical results for the tree surface position agrees within :t4 mm against the experimental data. The averaged axial velocity and the secondary flow fields are presented, as well as the turbulent kinetic energyand the dissipation of the turbulent kinetic energy. Finally the local and the averaged ftiction factors are presented. / Doutorado / Doutor em Engenharia Mecânica
|
33 |
Simulação de escoamento em coordenadas generalizadas utilizando o metodo dos volumes finitos com esquema UNIFAES de discretizaçãoVilela, Carlos Alberto de Almeida 31 July 2001 (has links)
Orientador: Jose Ricardo Figueiredo / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T18:04:25Z (GMT). No. of bitstreams: 1
Vilela_CarlosAlbertodeAlmeida_D.pdf: 6593008 bytes, checksum: 706d78c356008d4b1f84aa2b735170e6 (MD5)
Previous issue date: 2001 / Resumo: O esquema de discretização UNIFAES, é aqui avaliado e comparado com alguns dos mais tradicionais e utilizados esquemas de discretização para o método dos volumes finitos da atualidade. Diversos casos testes lineares do problema de transporte convectivo-difusivo em meio fluido e das equações de Navier-Stokes foram utilizados para comparação. Os esquemas utilizados para comparação foram o Central, Exponencial, LOADS e QUICK em duas versões diferentes de formulação. Os casos testes realizados para avaliação do comportamento numérico realizados em coordenadas generalizadas com malhas regulares e irregulares foram: Transporte convectivo-difusivo em, a) escoamento paralelo b) escoamento rotacional de fluido como sólido, c) escoamento em canto, d) escoamento em canto com perfil de entrada parabólico, escoamento em cavidade fechada inclinada e escoamento em canal divergente gradual. Os resultados obtidos nos testes demonstraram que o esquema UNIFAES, apresentou comportamento bastante promissor comparado aos comportamentos dos demais esquemas. Na grande maioria dos testes os resultados obtidos pelo esquema UNIFAES foram os melhores em questão de acuidade, estando próximo aos obtidos pelo esquema QUICK, que foi o segundo esquema que melhores resultados apresentou / Abstract: The UNIFAES discretization scheme is evaluate and compare with some classical and most used schemes for finite volume method in present time. A lot of linear convection-diffusion problems in fluid flow and the solution of Navier-Stokes equations were used for comparison. The compared used schemes were Central, Exponential, LOADS, and QUICK with two different formulations. The chosen test cases for comparisons, using regular and irregular grid arrangements in generalized coordinates formulations were: a) Convective-diffusive transport in parallel fluid flow, in rotational fluid flow, in fluid flow near a wall, in fluid flow near a wall with parabolic inlet profile of the transported property , b) Inclined side walls lid driven cavity, and c) Gradually expanding channel. The obtained results shown that UNIFAES scheme has a very promise behavior compared to the other schemes. In most studied cases, UNIFAES scheme demonstrated better accurate results, elose to QUICK, the second better one / Doutorado / Doutor em Engenharia Mecânica
|
34 |
Introdução a modelagem em 3D para reatores air-lift empregando o "metodo dos volumes finitos" na resolução da fluidodinamica (CFD)Barreira, Mardonny Nazareno 03 August 2018 (has links)
Orientador: Rubens Maciel Filho / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T15:03:24Z (GMT). No. of bitstreams: 1
Barreira_MardonnyNazareno_M.pdf: 5649911 bytes, checksum: 74e112ba543ced12ed31508868adfdde (MD5)
Previous issue date: 2003 / Mestrado
|
35 |
Analise parametrica de um trocador de calor atraves da simulação numericaSaito, Marcelo Batista 13 September 2002 (has links)
Orientador: Eugênio Spanó Rosa, Marcelo Moreira Ganzarolli / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T15:19:17Z (GMT). No. of bitstreams: 1
Saito_MarceloBatista_M.pdf: 6699986 bytes, checksum: 845433a3bcb3aa43a33251e4350ebf40 (MD5)
Previous issue date: 2002 / Resumo: O presente trabalho tem por objetivo determinar numericamente a queda de pressão e o fluxo de calor de um trocador de calor ar-ar utilizado para resfriamento de um ambiente selado que não pode ser contaminado pelo ar externo. Esta análise tem por finalidade otimizar a geometria que maximiza o calor transferido para uma determinada queda de pressão. O trocador é formado por dutos paralelos de seções retangulares com escoamento em contra-corrente no regime turbulento. As duas correntes são separadas por paredes sólidas e lisas. A análise paramétrica é efetuada resolvendo-se equações de conservação da massa, quantidade de movimento e energia, uti1izando o modelo de turbulência K-E padrão e considerando propriedades variáveis. A condução de calor nas paredes sólidas também foi incluída nas simulações. A solução numérica tridimensional é obtida através da discretização das equações de conservação em volumes finitos, sendo empregado um pacote numérico comercial (pHOENICS 3.3). As simulações são realizadas para número de Reynolds entre 8000 e 160000. Os resultados obtidos permitem avaliar o desempenho do trocador de calor para valores especificados da queda de pressão e determinar a geometria que maximiza o calor transferido / Abstract: The pressure drop and heat transfer flux on an air-air heat exchanger employed as a cooling device for a closed cavity with heat dissipation are determined numerically. The analysis is further extended to optimize the cross section dimensions to render the maximum heat transfer rate at a given pressure drop. The heat exchanger is of rectangular cross section ducts placed side by side in counter flow turbulent regime. The duct walls are of smooth and flat aluminum. The conservation equations of mass, momentum and energy are solved with the K-B turbulence model using the finite volume code embodied on the PHOENICS 3.3. The heat conduction along the solid walls is considered employing conjugate heat transfer treatment. The simulations are performed for Reynolds number between 8000 and 160000. The pressure drop and heat transfer rate are determined and expressed in terms of dimensionless parameters. The analysis is further advanced determining cross section aspect ratio which gives the maximum heat transfer rate for a given pressure dro / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica
|
36 |
Modelagem matematica e simulação numerica do escoamento liquido-vapor num prato de destilaçãoNoriler, Dirceu, 1978- 04 January 2003 (has links)
Orientadores: Maria Regina Wolf Maciel, Henry França Meier, Antonio Andre Chivanga Barros / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T15:31:22Z (GMT). No. of bitstreams: 1
Noriler_Dirceu_M.pdf: 4837242 bytes, checksum: ed71f25bb9219c76441c2042ff4f6a07 (MD5)
Previous issue date: 2003 / Resumo: Os modelos para representação de colunas de destilação atuais baseiam-se no conceito de estágios de equilíbrio e de não equilíbrio, ambos considerando o equilíbrio hidrodinâmico em uma abordagem macroscópica. No entanto, em meados da década de 90, iniciou-se a aplicação de modelos microscópicos para subsidiar e incrementar a análise e otimização dos processos de refino, dando-se início à era dos códigos comerciais de Fluidodinâmica Computacional. O principal objetivo deste trabalho é desenvolver um modelo microscópico para descrever a fluidodinâmica do escoamento líquido-vapor em um prato perfurado de uma coluna de destilação, com base na conservação de quantidade de movimento sob condições de turbulência e aplicar e implementar o modelo matemático no código comercial de fluido dinâmica computacional CFX, da ABA Technology, para a realização de um estudo de experimentação numérica de modelos e métodos. Foram propostos três modelos distintos para representação do escoamento em um prato perfurado de destilação: um modelo bidimensional e monofásico, 2-DSM (Liu et ai. (1999)) , um modelo tridimensional bifásico e homogêneo, 3-DHomM, e por fim um modelo tridimensional bifásico e heterogêneo, 3-DHetM. Desta forma, os três modelos foram implementados no código comercial CFX 4.4 da ABA technology, o qual utiliza o método dos volumes finitos com variáveis co-localizadas resolvendo o sistema em coordenadas generalizadas. As malhas numéricas foram construídas por intermédio da técnica dos multiblocos, utilizando o pré-processador "BUILD" presente no código. Os principais resultados mostram que o modelo 3-DTHomM representa, qualitativamente, de forma satisfatória a hidrodinâmica de um prato perfurado de destilação. A metodologia proposta neste trabalho se mostrou adequada e apta a reproduzir tal comportamento fluidodinâmico / Abstract: The current models for predicting distillation column behavior are based on equilibrium and non-equilibrium stage concepts, but both use the hydro dynamical equilibrium in a macroscopic point of view. The objectives of this work are: developing the microscopic model to describe the gas-liquid flow on a distillation sieve tray, based on momentum equations under turbulence conditions, inter-phase momentum transfer and other modeling features; and applying the mathematical modeling in a commercial CFD code, for numerical experiments in order to verify the software and to validate the models. Three models were proposed: the first model named Two Dimension Single Model (2-DSM model), with two dimension and single phase assumptions; the second one, named Three Dimension Homogeneous Model (3-DHomM model), with two phases and three dimension with homogeneous flow, i. e., the liquid and the gas velocities are the same; and the last model is a model with two phases and three dimension assumptions, named Three Dimension Heterogeneous Model (3-DHetM model). These three models were applied in CFX 4.4 CFD commercial code from ABA Technology, by using the finite-volume methods with co-located variable in a generalized coordinate system, to obtain numerical solutions. The pre-processor BUILD were used to build up the numerical grid with the multi-block technique. The main results show that the 3-DHomM model, despite of its simplicity when compared with 3-DhetM, describes the fluid dynamics on a sieve distillation tray satisfactorily. The methodology proposed in this work also showed to be adequate, and with capability to describe the gas-liquid flow on a sieve distillation tray / Mestrado / Engenharia de Processos / Mestre em Engenharia Química
|
37 |
Métodos de volumes finitos robustos para a simulação de escoamentos bifásicos de água e óleo em reservatórios de petróleoCONTRERAS, Fernando Raul Licapa 10 January 2017 (has links)
Submitted by Pedro Barros (pedro.silvabarros@ufpe.br) on 2018-09-21T19:43:58Z
No. of bitstreams: 2
license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5)
TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5) / Approved for entry into archive by Alice Araujo (alice.caraujo@ufpe.br) on 2018-09-24T23:13:56Z (GMT) No. of bitstreams: 2
license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5)
TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5) / Made available in DSpace on 2018-09-24T23:13:56Z (GMT). No. of bitstreams: 2
license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5)
TESE Fernando Raul Licapa Contreras.pdf: 6260685 bytes, checksum: d6d66a2bde6c2b142b3b9f666d333f4e (MD5)
Previous issue date: 2017-01-10 / FACEPE / No presente trabalho apresentamos formulações numéricas para a discretização de um modelo matemático que descreve o escoamento bifásico de água e óleo em reservatórios de petróleo heterogêneos e anisotrópicos. Estas equações são resolvidas através da metodologia IMPES (Implicit Pressure Explicit Saturation) que é, em particular, um método segregado para tratar escoamento bifásico (água-óleo) em que, o campo de pressão é obtida implicitamente e o campo de saturação é obtida explicitamente. Para resolver a equação de pressão propormos formulações numéricas baseadas no método de volumes finitos lineares (a saber: MPFA-HD e MPFA-H) e no método de volumes finitos não lineares (a saber: NLFV-PP e NLFV-DMP), estas formulações em geral lidam com meios altamente heterogêneos e anisotrópicos sobre malhas poligonais em geral. De outra parte, neste trabalho são estudadas as varias propriedades numéricas que satisfazem as formulações numéricas propostas. Além disso, no caso do método não linear que satisfaz o principio do máximo discreto (DMP), mostramos que este possui um estêncil pequeno desde que os pontos de interpolação sejam calculados a partir dos pontos harmônicos. Para resolver a equação de saturação, utilizamos um método multidimensional nomeada MOOD (Multi-dimensional Optimal Order Detection), esta método é baseada numa discretização espacial de alta ordem com procedimento de limitação a posteriori das variáveis de estado (i. e. saturação) não físicos. Em cada volume de controle o grau polinomial ótimo é determinado para construir uma aproximação das soluções respeitando o principio do máximo discreto. Mediante um processo iterativo o grau polinomial é reduzido quando o volume de controle em questão não satisfaz as condições de estabilidade. Neste método, uma estratégia baseada na Técnica de Mínimos Quadrados é usada para a reconstrução polinomial. Todos métodos propostos são localmente conservativos, e todas as incógnitas são completamente centradas nos volumes de controle, e possuem a habilidade de representar reservatórios com propriedades geológicas complexas utilizando malhas poligonais quaisquer. O desempenho dessas foi avaliado a partir da solução de problemas relevantes encontrados na literatura. / In the present work we present numerical formulations for the discretization of a mathematical models that describes the two-phase flow of water and oil in heterogenous and anisotropic oil reservoirs. These equations are solved through Implicit Pressure Explicit Saturation (IMPES) methodology which is in particular a segregated method to treat two-phase (water-oil) flow in which the pressure field is implicitly obtained and the saturation field is obtained explicitly. In order to solve the pressure equation we propose numerical formulations based on the finite linear volume method (MPFA-HD and MPFA-H) and non-linear finite volume method (NLFV-PP and NLFV-DMP). Formulations generally deal with highly heterogeneous and anisotropic media over polygonal meshes in general. On the other hand, in this work the numerical properties that satisfy the proposed numerical formulations are studied. In addition, in the case of the non-linear method that satisfies the discrete maximum principle (DMP), we show that it has a small stencil since the interpolation points are calculated from the harmonic points. To solve the saturation equation, we use a multidimensional method named MOOD (Multidimensional Optimal Order Detection), this method is based on a high order spatial discretization with a posteriori limitation procedure of the non-physical state variables (saturation). In each control volume the optimum polynomial degree is determined to construct an approximation of the solutions respecting the discrete maximum principle. By means of an iterative process the polynomial degree is reduced when the control volume in question does not satisfy the stability conditions. In this method, a strategy based on the Least Squares technique is used for polynomial reconstruction. All proposed methods are locally conservative, and all unknowns are completely centered on control volumes, and have the ability to represent reservoirs with complex geological properties using any polygonal meshes. The performance of these was evaluated from the solution of relevant problems found in the literature.
|
38 |
[en] TURBULENT COMPLEX FLOW SIMULATION WITH CLASSICAL MODELING AND LARGE EDDY SIMULATION / [pt] SIMULAÇÃO DE ESCOAMENTO TURBULENTO COMPLEXO COM MODELAGEM CLÁSSICA E DE GRANDES ESCALASJOSE DINIZ MESQUITA ABRUNHOSA 06 January 2004 (has links)
[pt] Uma investigação da capacidade de previsão de modelos de turbulência baseados na modelagem estatística clássica e de grandes escalas é apresentada. A modelagem estatística clássica de turbulência (média de Reynolds) foi analisada, através da solução de escoamentos complexos, como, por exemplo, o escoamento turbulento em degrau (backstep). Especial atenção foi dada aos modelos kapa-epsilon de baixo Reynolds
e as variantes renormalizadas (RNG). O comportamento dos vários termos da equação da energia cinética turbulenta na região da parede foram analisados em detalhes, especialmente o termo de difusão de pressão. Avaliou-se a importância da correta modelagem do termo de difusão de pressão sobre as predições dos modelos de baixo número de Reynolds, nas regiões de recirculação. Alguns modelos,
propostos na literatura para o termo de difusão de pressão, foram também avaliados teórica e numericamente. A capacidade de previsão da metodologia de simulação de grandes escalas (LES por Large Eddy Simulation) também foi realizada. O desempenho do modelo de Smagorinsky para prever escoamentos limitados por fronteiras sólidas foi avaliado do ponto de vista computacional. Utilizou-se o método de volumes finitos para integrar tanto as equações médias de Reynolds quanto as equações LES. O escoamento turbulento em canal foi resolvido de modo bidimensional e tridimensional. Já o escoamento em degrau (backstep) foi resolvido exclusivamente de modo bidimensional, enquanto o escoamento em um duto de seção quadrada foi simulado de modo tridimensional. Os resultados foram comparados com
aqueles obtidos pelos modelos de baixo Reynolds, analisando-se a relação custo-benefício. / [en] An investigation of turbulence models prediction capacity
based on classical statistical modeling and large eddy
simulation (LES) is presented. The classical statistical
modeling (average of Reynolds) was analyzed, by
investigating the solution of complex flows, as, for
example, the turbulent flow past a backwardfacing-
step (backstep). Special attention was given to low
Reynolds number k-e models and models derived by
renormalization group theory (RNG). The behavior of the
different terms in the turbulent kinetic energy equation in
the near wall region was examined in details, specially the
pressure diffusion term. It was evaluated the importance of
the correct modeling of the pressure diffusion term on
the predictions of the low Reynolds number models, in
recirculating flows. A few models, proposed in the
literature for the pressure diffusion term, were also
evaluated theoretically and numerically. The prediction
capacity of large eddy simulation (LES) technique was also
investigated. The ability of Smagorinsky model to predict
complex limited wall flows was analyzed from a computational
standpoint. The finite-volume method was employed to
integrate both the Reynolds average and LES equations. The
fully developed turbulent channel flow was solved in two-
dimensional and three-dimensional numerical simulations. The
turbulent flow over a backward-facing-step was computed
exclusively in a twodimensional manner, while the fully
developed turbulent flow in a straight square duct was
simulated in a three-dimensional manner. The results were
compared with those obtained by the low Reynolds models,
analyzing the cost-benefit relation.
|
39 |
Eulerian-lagrangian simulation of a turbulent evaporating sprayRodrigo Badia Piccinini 25 November 2011 (has links)
This work consists in a numerical simulation of a spray jet: a two-phase flow composed by a gaseous and a liquid phase. The gaseous phase is treated as a continuous medium, and the liquid phase is treated as a dispersed phase. An asymptotic approximation of zero Mach number was applied to the gaseous phase in order to account for density variations due to temperature gradients without dealing with extra complexities of the fully compressible flow formulation. The effects of turbulence on the gas flow were modeled using the concept of turbulent viscosity determined by a system of two partial differential equations, the so called k-epsilon model. The liquid dispersed phase was modeled as point-particles, to whom thermodynamic properties and 0-d models for momentum, heat and mass transfers were assigned. The gaseous phase was discretized and numerically solved using the finite volume method. The numerical results were compared to measurements and a reasonable prediction was found. It was verified, though, underestimations of droplet velocity and evaporation rate.
|
40 |
Estruturas de dados topológicas aplicadas em simulações de escoamentos compressíveis utilizando volumes finitos e métodos de alta ordem / Topologic data structures applied on compressible flows simulations using finite volume and high-order methodsBarbosa, Fernanda Paula 18 December 2012 (has links)
A representação de malhas por meio de estrutura de dados e operadores topológicos e um dos focos principais da modelagem geométrica, onde permite uma implementação robusta e eficiente de mecanismos de refinamento adaptativo, alinhamento de células e acesso as relações de incidência e adjacência entre os elementos da malha, o que é de grande importância na maioria das aplicações em mecânica dos fluidos. No caso de malhas não estruturadas, a não uniformidade da decomposição celular e melhor representada por uma estrategia mais sofisticada, que são as estruturas de dados topológicas. As estruturas de dados topológicas indexam elementos de uma malha representando relações de incidência e adjacência entre elementos, garantindo acesso rápido às informações. Um dos aspectos mais comuns aos problemas tratados pela mecânica dos fluidos computacional é a complexidade da geometria do domínio onde ocorre o escoamento. O uso de estruturas de dados para manipular malhas computacionais e de grande importância pois realiza de modo eficiente as consultas às informações da malha e centraliza todas as operações sobre a malha em um único módulo, possibilitando sua extensão e adaptação em diversas situações. Este trabalho visou explorar o acoplamento de uma estrutura de dados topológica, a Mate Face, em um módulo simulador existente, de modo a gerenciar todos os acessos à malha e dispor operações e iteradores para pesquisas complexas nas vizinhanças de cada elemento na malha. O módulo simulador resolve as equações governantes da mecânica dos fluidos através da técnica de volumes finitos. Foi utilizada uma formulação que atribui os valores das propriedades aos centroides dos volumes de controle, utiliza métodos de alta ordem, os esquemas ENO e WENO, que tem a finalidade de capturar com eficiência descontinuidades presentes em problemas governados por equações diferenciais parciais hiperbólicas. As equações de Euler em duas dimensões representam os escoamentos de interesse no presente trabalho. O acoplamento da estrutura de dados Mate Face ao simulador foi realizada através da criação de uma biblioteca desenvolvida que atua como uma interface de comunicação entre os dois módulos, a estrutura de dados e o simulador, que foram implementados em diferentes linguagens de programação. Deste modo, todas as funcionalidades existentes na Mate Face tornaram-se acessíveis ao simulador na forma de procedimentos. Um estudo sobre malhas dinâmicas foi realizado envolvendo o método das molas para movimentação de malhas simulando-se operações de arfagem. A idéia foi verificar a aplicabilidade deste método para auxiliar simulações de escoamentos não estacionarios. Uma outra vertente do trabalho foi estender a estrutura Mate Face de forma a representar elementos não suportados a priori, de modo a flexibilizar o seu uso em simulações de escoamentos baseados no método de volumes finitos espectrais. O método dos volumes espectrais e utilizado para se obter alta resolução espacial do domínio computacional, que também atribui valores das propriedades aos centroides dos volumes de controle, porém, os volumes de controle são particionados em volumes menores de variadas topologias. Assim, uma extensão da Mate Face foi desenvolvida para representar a nova malha para a aplicação do método, representando-se cada particionamento localmente em cada volume espectral. Para todas as etapas deste trabalho, realizaram-se experimentos que validaram a utilizaação da estrutura de dados Mate Face junto a métodos numéricos. Desta forma, a estrutura pode auxiliar as ferramentas de simulações de escoamentos de fluidos no gerenciamento e acesso à malha computacional / The storage and access of grid files by data structures and topologic operators is one of the most important goals of geometric modeling research field, which allows an efficient and stable implementation of adaptive refinement mechanisms, cells alignment and access to incidence and adjacency properties from grid elements, representing great concernment in the majority of applications from fluid mechanics. In the case of non-structured grids, the cellular decomposition if non-uniform and is better suited by a more sophisticated strategy - the topologic data structs. The topologic data structs index grid elements representing incidence and adjacency properties from grid elements, ensuring quick access to information. One of most common aspects from problems solved by computational fluid mechanic is the complexity of the domain geometry where the fluid ows. The usage of data structures to manipulate computational grids is of great importance because it performs efficiently queries on grid information and centers all operations to the grid on a unique module, allowing its extension and flexible usage on many problems. This work aims at exploring the coupling of a topologic data structure, the Mate Face, on a solver module, by controlling all grid access providing operators and iterators that perform complex neighbor queries at each grid element. The solver module solves the governing equations from fluid mechanics though the finite volume technique with a formulation that sets the property values to the control volume centroids, using high order methods - the ENO and WENO schemes, which have the purpose of efficiently capture the discontinuities appearing in problems governed by hyperbolic conservation laws. The two dimensional Euler equations are considered to represent the flows of interest. The coupling of the Mate Face data structure to the solver module was achieved by a creation of a library that acts as an interface layer between both modules, the Mate Face and the solver, which had been implemented using different programming languages. Therefore, all Mate Face class methods are available to the solver module though the interface library in the form of procedures. A study of dynamic grids was made by using spring methods for the moving grid under pitch movement case. The goal was to analyze the applicability of such method to aid non stationary simulations. Another contribution of this work was to show how the Mate Face can be extended in order to deal with non-supported types of elements, allowing it to aid numeric simulations using the spectral finite volume method. The spectral nite volume method is used to obtain high spatial resolution, also by setting the property values to the control volume centroids, but here the control volumes are partitioned into smaller volumes of different types, from triangles to hexagons. Then, an extension of the Mate Face was developed in order to hold the new generated grid by the partitioning specfied by the spectral finite volume method. The extension of Mate Face represents all partitioned elements locally for each original control volume. For all implementations and proposals from this work, experiments were performed to validate the usage of the Mate Face along with numeric methods. Finally, the data structure can aid the fluid flow simulation tools by managing the grid file and providing efficient query operators
|
Page generated in 0.0725 seconds