Orientador: Luiz Otávio Saraiva Ferreira / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-26T05:32:16Z (GMT). No. of bitstreams: 1 Horta_AgnusAzevedo_M.pdf: 4762311 bytes, checksum: 1bcc55cdb06dfc294f352fc1ac27c0ad (MD5) Previous issue date: 2014 / Resumo: Este trabalho tem como objetivo principal a implementação do motor de simulação de um framework de simulação com partículas, o motor de simulação utilizará o método DPD (Dissipative Particle Dynamics), baseando-se no paradigma de programação orientada a objeto (POO) e no uso de estruturas de dados otimizadas. O motor de simulação foi escrito em linguagem C++. A concepção do sistema foi realizada de forma a facilitar e promover a reutilização e manutenção do código. Buscou-se, também, a flexibilidade e generalização através do uso da linguagem Python na geração dos arquivos de entrada correspondentes a distribuição espacial das partículas, sendo utilizada a linguagem de marcação XML (eXtensible Markup Language) na estruturação dos arquivos resultantes da simulação. No final, o motor de simulação é avaliado aplicando o problema do fluxo de um fluído entre placas paralelas e o resultado comparado com os resultados obtidos no simulador Hoomd-Blue / Abstract: This work aims the implementation of the simulation of a simulation framework with particle engine, the simulation engine will use DPD (Dissipative Particle Dynamics) method, based on the paradigm of object oriented programming (OOP) and use optimized data structures. The simulation engine is written in C language ++. The system design was performed in order to facilitate and promote the reuse and maintainability of the code. Also, we sought flexibility and generalization through the use of Python in the generation of the input files corresponding to the spatial distribution of particles, (eXtensible Markup Language) XML markup being used in structuring the files resulting from the simulation. At the end, the simulation engine is the problem of applying rated flow of a fluid between parallel plates and the result compared with the results obtained in the simulation Hoomd-Blue / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-26T06:44:29Z (GMT). No. of bitstreams: 1 AlvarezCastro_HelverCrispiniano_D.pdf: 5280150 bytes, checksum: 96fb29e2849d11c06890186cb8af2f51 (MD5) Previous issue date: 2014 / Resumo: Conversão de matéria-prima e rendimento são estudados através de um modelo 3D que simula o principal reator do processo de craqueamento catalítico fluidizado (FCC). Fluidodinâmica computacional (CFD) é usada como técnica com abordagem Euleriana-Euleriana para prever o comportamento do riser da FCC. O modelo considera 12 grupos com desativação do catalisador dada por coque e envenenamento por nitritos alcalinas e adsorção aromático policíclico para estimar o comportamento cinético, o qual partindo de uma determinada matéria prima, se craquea a diversos produtos. Composições diferentes para matérias-primas foram consideradas. Para validação o modelo foi comparado com dados experimentais em condições de operação industrial. O modelo é capaz de representar com precisão a geração de produtos para as diferentes condições de operação consideradas. Todas as condições consideradas foram resolvidos usando o solver Ansys CFX 14.0. As diferentes variáveis do processo de operação e efeitos hidrodinâmicos do riser industrial de craqueamento catalítico fluidizado (FCC) foram avaliadas. As previsões do modelo são apresentados e comparadas; conversão e produtos são apresentados. Recomendações para criar condições para obter maior rendimento do produto no processo industrial foram dadas / Abstract: Feedstock conversion and yield products are studied through a 3D model simulating the main reactor of the fluid catalytic cracking (FCC) process. Computational Fluid Dynamic (CFD) is used with Eulerian-Eulerian approach to predict the fluid catalytic cracking behavior. The model considers 12 lumps with catalyst deactivation by coke and poisoning by alkaline nitrides and polycyclic aromatic adsorption to estimate the kinetic behavior which, starting from a given feedstock, produces several cracking products. Different feedstock compositions are considered. The model is compared with sampling data at industrial operation conditions. The simulation model is able to represent accurately the products behavior for the different operating conditions considered. All the conditions considered were solved using a solver Ansys CFX 14.0. The different operation process variables and hydrodynamic effects of the industrial riser of a fluid catalytic cracking (FCC) are evaluated. Predictions from the model are shown and comparison with experimental conversion and yields products are presented; recommendations are drawn to establish the conditions to obtain higher product yields in the industrial process / Doutorado / Engenharia de Processos / Doutor em Engenharia Química

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-21T06:39:27Z (GMT). No. of bitstreams: 1 GonzalezSilva_German_D.pdf: 4427269 bytes, checksum: bb45c122c3faad0613f060ad4cbd3b61 (MD5) Previous issue date: 2012 / Resumo: O objetivo principal desta tese é propor uma metodologia para aplicar simulação de grandes escalas (LES) em uma unidade de craqueamento catalítico industrial. Para atingir este objetivo, iniciou-se propondo uma metodologia para a construção da malha, sendo uma malha quase-uniforme. As malhas propostas foram implementadas para minimizar o esforço computacional e procedimento para a refinação uniforme no domínio do sistema. Inicialmente, foi estudada a fluidodinâmica de um leito fluidizado gás-sólido na escala de laboratório sem reação química, utilizando Simulação de Grandes Escalas. Com base nos resultados numéricos obtidos a partir de resultados da escala de laboratório foi aplicado o modelo cinético de 4 lump em uma simulação CFD tridimensional para um reator industrial FCC, utilizando LES para a fase gasosa, e considerando o catalisador como uma fase contínua (Euleriana). Os pacotes de simulação utilizados foram Ansys ICEM versão 13 para a construção da malha e Ansys CFX versão 13 para o pós-processamento dos resultados. No pós-processamento dos resultados foi proposta uma metodologia para determinar as médias azimutais das variáveis em planos perpendiculares ao escoamento e finalmente foram validadas as simulações com dados reportados na literatura. As principais conclusões do trabalho foram que foi possível diminuir o tamanho da malha e o tempo de processamento. Notou-se também que, ao usar malhas com discretização quase-uniforme não foi necessário fazer um refinamento de malha elevado, nem refinar perto da parede para o sistema gás-sólido / Abstract: The main objective of this thesis is to propose a methodology in how to apply Large Eddy Simulation (LES) on a unit of catalytic cracking. In order to achieve this, it was proposed an alternative way of constructing the computational mesh, by using quasi-uniform meshes. The proposed meshes were implemented to minimize the computational effort and procedure for refining them in the entire domain of the system. Initially it was studied the fluid dynamics of a lab scale gas-solid system without chemical reaction, using Large Eddy Simulation. Based on the numerical results obtained from lab scale results it was implemented the 4 lump kinetic model in a three dimensional CFD simulation of an FCC industrial reactor, using LES for the gas phase and considering the catalyst as a continuous phase (Eulerian). The simulation packages used were Ansys ICEM, version 13 for mesh construction and Ansys CFX 13 for computation and post-processing of the results. In the data post-processing it was proposed a methodology for calculating average values of fluctuating variables between two circular sections in the azimuthal direction. The results were compared with data reported in literature. The main conclusions of the results showed that it was possible to decrease the mesh size and the computational time. It was also noticed that by using quasi-uniform discretization it was not necessary to make a high mesh refinement near the wall for a gas-solid system / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientadores: José Roberto Nunhez, Everton Moraes Matos / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-21T22:56:39Z (GMT). No. of bitstreams: 1 Fontoura_DienerVolpinRibeiro_D.pdf: 11658919 bytes, checksum: 6a2be1fd48301240f5f4b18c84757db1 (MD5) Previous issue date: 2013 / Resumo: O escoamento de fluidos é, sem dúvida, de grande importância na indústria, especialmente no transporte de fluidos entre as unidades de operação. Sistemas multifásicos são bastante comuns, e em muitos casos, também ocorrem processos de transferência de massa e energia, como é o caso dos fornos de pré-aquecimento na indústria petroquímica. O escoamento gás-líquido nessas linhas é muitas vezes complexos, pois durante a vaporização da carga ocorre também o craqueamento térmico do petróleo. É comum nestas operações a formação de coque no interior do tubo, o que é muito indesejável. Para entender as condições nas quais esta formação de coque é minimizada, ou aumentada, foi proposto um modelo uidodinâmico computacional para simular o escoamento bifásico gás - liquido do petróleo no interior desses tubos. O modelo foi implementado na plataforma OpenFOAM, ao qual foram acrescidas novas rotinas para a estimativa das temperaturas e concentrações através das equações de conservação de massa e energia. O modelo k - epsilon foi utilizado para descrever a turbulência e dois modelos de vaporização foram utilizado para caracterizar a mudança de fase. Uma formulação não conservativa foi adotada na escrita das equações do modelo para maior estabilidade numérica em altas frações volumétricas. Foi também implementada uma rede cinética dependente da temperatura para descrever o craqueamento térmico. Através do modelo foi possível simular o escoamento tridimensional de petróleo no interior de um tubo com as dimensões da planta obtendo como resultado os perfis de velocidades, temperaturas e concentrações das fases líquida e gasosa / Abstract: Fluid flow is of great importance in the industry, especially for its transport between the operating units. Multiphase systems are very frequent. In many cases, mass and energy transfer processes, such as the case of preheating furnaces in petrochemical industry, are also observed. The gas-liquid flow in these lines is often complex, because during the vaporization of the charge it also occurs the petroleum thermal cracking. It is common in these operations the very undesirable formation of coke inside the tube. With the aim of understanding under which conditions the formation of coke is minimized or increased it was proposed a computational fluid dynamic model for simulating the petroleum gas-liquid two-phase flow inside of these tubes. The model was implemented on the OpenFOAM software, to which new routines for the calculation of temperatures and concentrations using the equations of conservation of mass and energy were added. The model k - epsilon was used to describe the turbulence and two vaporization models were used to characterize the phase change. A non-conservative formulation was adopted for describing the equations of the model in order to obtain numerical stability at high volumetric fractions. It was also implemented a temperature dependent kinetic net to describe the thermal cracking. With this model it was possible to simulate the three-dimensional flow of petroleum inside a tube with real industrial dimensions as a result of the gas and liquid phase's profiles of velocities, temperatures and concentrations / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientadores: Sérgio Tonini Button, José Divo Bressan / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-21T04:08:51Z (GMT). No. of bitstreams: 1 Martins_MarceloMatos_D.pdf: 11233728 bytes, checksum: 163e9bff5c80ef1c025e8c1083bed472 (MD5) Previous issue date: 2012 / Resumo: A simulação numérica computacional é nos dia de hoje frequentemente aplicada na elaboração de projetos ou análise dos processos de conformação plástica dos metais. A extrusão de metais é um dos principais processos de conformação plástica e largamente aplicado na fabricação de produtos e peças na indústria metal-mecânica. Tradicionalmente, essas análises são feitas utilizando o Método dos Elementos Finitos. Entretanto, há um aumento no interesse dos pesquisadores na utilização do Método dos Volumes Finitos para este fim. A literatura sugere que o escoamento na extrusão de metais pode ser analisado pela formulação do escoamento plástico (flow Formulation). No qual, pode-se assumir como o escoamento de um fluido incompressível e viscoso. Essa hipótese pode ser assumida já que o processo de extrusão é um processo isocórico. O método MacCormack é geralmente aplicado para simular os escoamentos de fluidos compressíveis pelo Método do Volumes Finitos. No escoamento de um fluido incompressível ou no escoamento de metal não existe uma equação para a evolução da variável pressão, sendo necessário a utilização de um método de acoplamento entre a pressão e a velocidade. Este trabalho trata da apresentação de um novo esquema numérico para a determinação de informações sobre o escoamento de um fluido incompressível e viscoso e sobre o escoamento de metal em um processo de extrusão direta, ambos em regime permanente. As equações governantes foram discretizadas pelo Método dos Volume Finitos através do Método de MacCormack explícito para uma malha estruturada e co-localizada. O acoplamento entre a pressão e a velocidade foi feita pelo método SIMPLE. O novo esquema numérico foi aplicado em escoamentos incompressíveis e viscosos para a glicerina e em escoamento de metais em processos de extrusão direta para o chumbo e uma liga de alumínio. O escoamento da glicerina foi avaliado para o caso entre placas paralelas e em dutos circulares sob condição axissimétrica e obtiveram boa concordância em relação ao resultados analíticos. Os campos de velocidades obtidos para a extrusão de metal alcançaram rápida convergência, em torno de 20000 iterações, essa quantidade de iterações foi inferior a quantidade que a glicerina necessitou. Para todos os materiais analisados os resultados numéricos tiveram boa concordância em comparação com resultados analíticos e experimentais obtidas da literatura. O método MacCormack produziu resultados coerentes para o escoamento da glicerina e dos metais sem a necessidade da adição de viscosidade artificial, como sugere a sua definição. Portanto, os resultados numéricos sugerem que o método MacCormack com o SIMPLE pode ser aplicado na resolução de escoamentos de fluidos incompressíveis e na conformação de metais além da sua tradicional aplicação na resolução de escoamentos compressíveis / Abstract: Computational numerical simulation is nowadays largely applied in the design and analysis of metal forming process. Extrusion of metals is one main forming process largely applied in the manufacturing of metallic products or parts. Historically, the Finite Element Method has been applied for decades in extrusion analysis. However, recently in the academy, there is a trend to use Finite Volume Method: literature suggests that metal flow by extrusion can be analysed by the flow formulation. Thus, metal flow can be modelled such us an incompressible viscous fluid. This hypothesis can be assumed because extrusion process is an isochoric process. The MacCormack Method is commonly used to simulate compressible fluid flow by the finite volume method. However, metal extrusion and incompressible fluid flow do not present state equations for the evolution of pressure, and therefore, a velocity-pressure coupling method is necessary to obtain a consistent velocity and pressure fields. Present work proposes a new numerical scheme to obtain information about both incompressible viscous fluid flow and metal flow in the extrusion process, in steady state. The governing equations were discretized by Finite Volume Method, using the Explicit MacCormack Method to structured and collocated mesh. The SIMPLE Method was applied to attain pressure-velocity coupling. These new numerical scheme was applied to incompressible viscous fluid flow of glycerine and forward extrusion process of lead and an aluminium alloy. The numerical results for glicerine fluid flow for parallel plates and axisymmetric flow in circular tube cases had quite good agreement in relation to the analytical solutions. The incompressible metal extrusion velocity fields achieved faster convergence than for liquid glycerine after 20.000 iterations and a good agreement with analytical and experimental results obtained from literature. The MacCormack Method applied for both glycerine and metals produced consistent results without the need of artificial viscosity as employed by the compressible flow simulation approaches. Hence, the present numerical results also suggest that MacCormack Method and SIMPLE can be applied in the solution of incompressible fluid flow and metal forming processes in adition to the traditional application for compressible fluid flow / Doutorado / Materiais e Processos de Fabricação / Doutor em Engenharia Mecânica

Orientadores: José Roberto Nunhez, Nicolas Spogis / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-22T07:47:05Z (GMT). No. of bitstreams: 1 Vergel_JoseLuisGomez_M.pdf: 8156840 bytes, checksum: f0b5d7b4b0875d4716285f8ec3ce5fa2 (MD5) Previous issue date: 2013 / Resumo: O uso de malhas não estruturadas recentemente tem recebido uma considerável atenção para discretização do domínio computacional em simulações CFD (Vakili & Esfahany, 2009; Cabrales, et al., 2011). A malha não estruturada (comumente tetraédricas) permite uma adequada adaptação em geometrias irregulares e uma geração automática na maioria dos softwares comerciais em CFD. Na literatura constantemente tem sido sugerido que as malhas hexaédricas são preferidas as tetraédricas, porém esta última pode ser usada, desde que alguns cuidados sejam tomados para proporcionar uma boa representação do fenômeno físico do problema. Alguns trabalhos recentemente publicados mostram que as malhas tetraédricas quando utilizadas com cuidado, podem obter resultados satisfatórios (Spogis & Nunhez, 2009). A maioria dos trabalhos desenvolvidos com uso da simulação CFD para processos de mistura utilizam malha hexaédrica, enquanto os resultados obtidos com malhas tetraédricas não possuem a mesma aprovação (Joaquim Junior et al., 2007), parte deste trabalho esta focado em mostrar que elas podem fornecer bons resultados em simulações de processos de mistura. Outro objetivo desse trabalho é mostrar o efeito do refinamento em malhas tetraédricas. Investigou-se a influência dos elementos prismáticos perto da parede sobre parâmetros globais, tais como o número de potência e número de bombeamento. Foi também realizada uma análise inicial da sensibilidade de alguns modelos de turbulência e esquemas de discretização sobre o campo de fluxo produzido. Além disso, simulações com diferentes refinamentos de malha foram feitas usando a abordagem de múltiplos sistemas de referência (MFR); enquanto em outros casos foi proporcionada uma abordagem de malhas deslizantes (SG) no impelidor. Os esquemas de discretização foram limitados aos esquemas Upwind, High resolution e blend fator, uma vez que esquemas diferentes não funcionam adequadamente com elementos tetraédricos. O software comercial de CFD CFX 14.0 foi utilizado para simular os resultados. Um impelidor de pás inclinadas (PBT 45°, com bombeamento para baixo) foi simulado em regime turbulento. Os resultados obtidos pelas simulações CFD para o número de potência e bombeamento, e perfil da componente axial da velocidade são discutidos e comparados com dados experimentais (Machado et al. et al., 2011). A velocidade foi medida usando a técnica PIV (Particle Image Velocimetry). Uma significativa influencia foi observada em alguns parâmetros na determinação de importantes variáveis nos processos de mistura. De forma geral, a predição do perfil da velocidade é influenciada, de alguma forma, pela resolução da malha, pelo modelo de turbulência e pelo esquema de discretização / Abstract: The use of non-structured meshes has received recently considerable attention for discretization of the computational domain in CFD simulations (Vakili & Esfahany, 2009; Cabrales et al., 2011). The non-structured meshes (the most common are tetrahedral) allow adequate adaptation of irregular geometries and easy use of automated algorithms for mesh generation in most commercial CFD software. It has been consistently suggested in the literature that hexahedral meshes are preferred over tetrahedral. However, they can be used, provided that some care is taken to provide a good representation of the physical phenomena of the problem. Recently published works have been proving that if they are used with care, satisfactory results are obtained (Spogis & Nunhez, 2009). Most of works developed with the use of CFD simulation for mixing processes use hexahedral and to date tetrahedral do not have the same acceptance of hexahedral (Joaquim Junior et al., 2007) and part of this work is aimed at showing that they can provide good results when simulating mixing processes. Another objective of this work is to show the effect of mesh refinement in tetrahedral meshes. It is investigated the influence of the prismatic elements near wall on global parameters such as the Power and Flow numbers. An initial analysis of the sensitivity of some turbulence models (Standard k-? model, the shear stress transport (SST) and BSL-Reynolds stress model) and the discretization scheme on the flow field produced were also carried out in this work. In addition, simulations with different mesh resolutions were made using the multiple reference frame (MFR) approach; whereas other cases a sliding mesh modeling was provided for the impeller. The discretization schemes were limited to upwind, high resolution and blend factor, since other schemes do not work well with tetrahedral. The commercial CFD software CFX 14.0 was used to simulate the results. A pitched blade turbine (PBT 45°, down-pumping), was simulated in turbulent flow. The results obtained for the Power and Flow numbers, the axial velocity profile component of the model in CFD are discussed and compared with experimental data (Machado et al. et al., 2011). The velocity was measured using the PIV technique (Particle Image Velocimetry). A significant influence of some parameters is observed on the determination of some important variable mixing. In general, the predictions of the velocity profile are influenced in some way by the resolution of the mesh, turbulence model and the discretization scheme / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

This work presents the result achieved through investigation of RH-2 gas-lift valve dynamic performance. It is an unbalanced gas-lift valve injection pressure operated, for high injection pressure applications. First, a CFD model was developed in order to obtain the valve orifice flow pattern performance curve. The CFD curve was compared with the curve experimentally obtained, conceiving very promising results. A test bench, and the necessary test devices, have been projected and built looking forward to determine the bellows load rate. Second, a CFD mesh deformation model was developed, which allowed the determination of dynamic performance curves on orifice, transition and throttling flow patterns for a given bellows pressure. The accuracy of the mesh deformation model was verified comparing the obtained curve with experimental results. The results of numerous CFD model simulations allowed to correlate data and to find a transition pressure, the production pressure for the maximum gas throughput and the product of the discharge coefficient by the expansion factor. The mathematical model developed allows predicting the gas throughput capacity of the gas-lift valve using injection pressure, production pressure and bellows pressure, for a defined load rate and a defined port diameter. / Nos primeiros estágios de produção, um poço de petróleo apresenta um fluxo natural, sendo a energia necessária para elevação dos fluidos até a superfície fornecida pelo próprio reservatório. Neste caso, diz-se que o poço é surgente. Entretanto, na medida em que as reservas são produzidas, a pressão do reservatório diminui e surge a necessidade de complementar artificialmente a energia requerida para a elevação. Diz-se, então, que o poço produz por elevação artificial. De todos os métodos de elevação artificial, o gas lift é o mais amplamente utilizado, e para a sua implementação são requeridas as chamadas válvulas de gas lift. O conhecimento do comportamento dinâmico da válvula de gas lift é fundamental na etapa de projeto e no estabelecimento de estratégias operacionais eficientes. Este trabalho apresenta os resultados obtidos no estudo do desempenho dinâmico da válvula de gas-lift RH-2, a qual consiste em uma válvula de gas lift não balanceada operada por pressão e de alta pressão de carregamento no fole, fabricada pela Weatherford, e que não dispõe de um modelo para o levantamento de suas curvas de desempenho dinâmico. Inicialmente, foi desenvolvido um modelo fluidodinâmico utilizando-se Dinâmica de Fluidos Computacional (CFD) para a obtenção da curva de desempenho no regime de orifício e este foi comparado com dados obtidos experimentalmente, produzindo resultados muito promissores. Uma bancada de testes com os dispositivos necessários foi projetada e construída para o levantamento da taxa de carga do fole. A seguir, foi desenvolvido um modelo fluidodinâmico da válvula no aplicativo CFX, utilizando uma malha deformável, que permitiu a geração das curvas de desempenho nos regimes de orifício, transição e estrangulamento para uma dada pressão de carregamento do fole. A precisão da malha deformável foi verificada comparando a curva obtida para diversas posições da haste com resultados experimentais. Os resultados das várias simulações do modelo fluidodinâmico permitiram correlacionar dados e calcular a pressão de transição, a pressão de produção onde ocorre a vazão máxima e o produto do coeficiente de descarga pelo fator de expansão. O modelo matemático desenvolvido permite predizer a vazão de gás que escoa pela válvula de gas lift a partir da pressão de injeção, da pressão de produção e da pressão de carregamento no fole, para uma taxa de carga do fole e um diâmetro de porta definidos.

Orientador: Luiz Otávio Saraiva Ferreira / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-27T02:13:59Z (GMT). No. of bitstreams: 1 Oliveira_FabiolaMartinsCamposde_M.pdf: 13513306 bytes, checksum: fd54d4a28df68e3217dbd8316b6aeadd (MD5) Previous issue date: 2015 / Resumo: Recentes avanços na tecnologia de processadores multinúcleos e vários-núcleos popularizaram a computação paralela, acelerando a execução de programas e possibilitando a simulação de domínios maiores. Dentre os problemas complexos que requerem alta velocidade de processamento, os problemas de fluidodinâmica computacional se destacam, já que suas simulações tendem a ter um alto custo computacional e exigem grandes domínios de simulação. O método baseado nas Redes de Boltzmann é um método de fluidodinâmica computacional apropriado para o uso de paralelismo que vem ganhando destaque na comunidade científica. Embora haja trabalhos que explorem o paralelismo em GPU nesse método, um simulador eficiente na execução e visualização interativa ainda não foi explorado adequadamente. Assim, a proposta deste trabalho é implementar em GPU um simulador interativo de fluidodinâmica com o método das Redes de Boltzmann. Inicialmente, o simulador foi desenvolvido em linguagem C e foi paralelizado em CPU usando MPI. Em seguida, foi paralelizado em GPU usando CUDA e convertido para linguagem orientada a objetos em C++. Depois, a visualização interativa foi acrescentada utilizando técnicas como interoperabilidade entre CUDA e OpenGL, texturização 3D, fluxo programável da GPU, além de funções de interação com o usuário. O simulador foi validado para casos 2D e 3D em fluxos monocomponentes monofásicos. Além disso, para demonstrar o ganho de desempenho em velocidade de processamento de problemas paralelizados em relação a execuções sequenciais, um conjunto de testes com tamanhos crescentes de domínio foi desenvolvido. O resultado dessa comparação indicou que o simulador implementado em GPU com visualização interativa teve desempenho 71.3 vezes maior em relação à versão sequencial em CPU sem visualização interativa. Dessa forma, a abordagem de paralelização em GPU com visualização interativa mostrou-se muito adequada à execução de simulações fluidodinâmicas, sendo uma ferramenta útil no estudo de escoamentos fluídicos, capaz de executar inúmeros cálculos e lidar com a grande quantidade de memória exigida por simulações fluídicas / Abstract: Recent advances on multicore and many-core processor technology have popularized the parallel computing, accelerating program execution and enabling the simulation of larger domains. Within the complex problems that require a high processing speed, the computational fluid dynamics problems stand out, since their simulations tend to have high computational cost and demand large simulation domains. The method based on the Lattice Boltzmann is an appropriate computational fluid dynamics algorithm to explore parallelism that has been noteworthy in scientific community. Although there are several works that approach GPU parallelism in this method, an efficient simulator implementation and interactive visualization have not been explored adequately. Thus, the purpose of this work is to implement in GPU an interactive fluid dynamics simulator based on the Lattice Boltzmann method. First, the simulator was developed in C language and was parallelized in CPU using MPI. Next, it was parallelized in GPU using CUDA and converted into C++ object-oriented language. Then, the interactive visualization was added using techniques such as CUDA-OpenGL interoperability, 3D texturing, GPU programmable pipeline, and interaction features. The simulator was validated for 2D and 3D cases in single component, single phase flows. Besides that, to show the performance gain in processing velocity of parallelized problems in relation to sequential executions, a test set with increasing domain sizes was developed. This comparison result indicated the GPU-implemented interactive simulator was 71.3 times faster in relation to sequential CPU version without interactive visualization. Thereby, the GPU parallelization approach with interactive visualization showed to be very adequate to fluid dynamics simulations, being a useful tool in fluid flow study, capable of simulating numerous calculations and dealing with the large amount of memory required in fluidic simulations / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestra em Engenharia Mecânica

Uma família de esquemas upwind denominada FUS-RF (Family of Upwind Scheme via Rational Functions), que é derivada via funções racionais e dependentes de parâmetros, é proposta para o cálculo de soluções aproximadas de equações de conservação. A fim de ilustrar a capacidade dos novos esquemas, vários resultados computacionais para sistemas hiperbólicos de leis de conservação são apresentados. Esses testes mostram a inflluência dos parâmetros escolidos sobre a qualidade dos resultados numéricos. Fazendo o uso de alguns testes de padrões, comparação dos novos limitadores de fluxo correspondentes com o esquema bem estabelecido van Albada e esquema atual EPUS (Eight-degree Polynomial Upwind Scheme) é também realizada. Os testes numéricos realizados em transporte de escalares e problemas de dinâmica dos gases confirmam que alguns esquemas da família FUS-RF são não oscilatórios e fornecem resultados confiáveis quando perfis descontínuos são transportados. Um esquema particular dessa nova família de esquemas upwind é então selecionado e utilizado para resolver escoamentos complexos com superfícies livres móveis / A family of upwind schemes named as FUS-RF (Family of Upwind Scheme via Rational Functions), which is derived via rational functions and dependent of parameters, is proposed for computing approximated solutions of conservation equations. In order to illustrate the capability of the new schemes, several computational results for system of hyperbolic conservation laws are presented. These results clarify the influence of the chosen parameters on the quality of the numerical calculations. Using some standard test cases, comparison of the new corresponding limiters with the well established van Albada and the recently introduced EPUS (Eight-degree Polynomial Upwind Scheme) limiters is also done. Numerical tests on both scalar and gas dynamics problems confirm that some schemes of the FUS-RF family are non-oscillatory and yield sharp results when solving profiles with discontinuities. A particular upwind scheme of this new family is then slected and used for solving complex incompressible moving free surface flows

Submitted by JOSIANE SANTOS DE OLIVEIRA (josianeso) on 2018-12-21T13:22:18Z No. of bitstreams: 1 José Henrique Nazzi Ehms_.pdf: 4808489 bytes, checksum: b590c723912ffe1d5d90006182f8cff2 (MD5) / Made available in DSpace on 2018-12-21T13:22:18Z (GMT). No. of bitstreams: 1 José Henrique Nazzi Ehms_.pdf: 4808489 bytes, checksum: b590c723912ffe1d5d90006182f8cff2 (MD5) Previous issue date: 2018-10-31 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Armazenamento térmico apresenta grande potencial de utilização em diversas aplicações, como energia solar, climatização, conservação de alimentos ou aproveitamento de calor residual em processos industriais. O armazenamento térmico de calor latente é realizado com materiais de mudança de fase (PCM), através dos processos de fusão e solidificação. A representação destes processos através de simulação numérica é realizada com o acréscimo ao modelo matemático básico, composto pelas equações da conservação da massa, quantidade de movimento e energia, modelos para descrever o calor latente e a transição na velocidade entre as fases, tais como: Darcy STM (source term method), VVM (variable viscosity method) e SOM (switch-off method). No entanto, a grande maioria das pesquisas nesta área utiliza o primeiro método. Além disso, são poucos os estudos comparativos de diferentes métodos para descrever processos de mudança de fase. Assim, o objetivo do presente estudo é comparar os métodos Darcy STM, VVM e misto (utilizando elementos dos dois anteriores) na simulação numérica de processos de fusão e de solidificação do PCM RT27 no interior de uma esfera. O estudo foi realizado utilizando-se fluidodinâmica computacional, através do método dos volumes finitos. O Modelo numérico foi validado com resultados experimentais da literatura. Os resultados quantitativos e qualitativos de fração líquida mostram que o método Darcy STM é mais adequado ao processo de solidificação, enquanto o método VVM produz resultados mais próximos aos experimentais no processo de fusão. O custo computacional foi menor para o método Darcy STM e maior para o método VVM, enquanto o método misto apresenta custo computacional pouco inferior ao do método VVM. Na análise da camada líquida na fusão de contato, foram analisadas a velocidade descendente do sólido, a espessura da camada e a vazão na camada. Os resultados referentes a camada líquida indicaram significativa influência das configurações do método Darcy STM. No entanto, são pouco influenciados pelas configurações do método VVM. / Thermal energy storage presents great potential of utilization in several applications, such as solar energy, HVAC systems, food conservation or waste heat recovery in industrial processes. Latent heat thermal energy storage is realized with phase change materials (PCM), through solidification and melting processes. Representation of such processes through numerical simulation is performed with the addition to the basic numerical model, composed of the conservation equations of mass, momentum and energy, models to account the latent heat and the velocity transition between the phases, such as: Darcy STM (source term method), VVM (variable viscosity method) and SOM (switch-off method). However, the large majority of the research on such area employ the first method. Besides that, there are few comparative studies of different methods to describe phase change processes. Thus, the objective of the present work is compare Darcy STM, VVM and mixed method (using elements of the two prior) in the simulation of melting and solidification processes of PCM RT27 inside a sphere. The study was realized using computational fluid dynamics, with the finite volume method. The numerical model was validated with experimental results from literature. Quantitative and qualitative results of liquid fraction show that Darcy STM is most suitable to solidification process, while VVM produces results closer to experimental in the melting process. Computational cost was smaller for Darcy STM and greater for VVM, while mixed method presents computational cost slightly lower than the one of VVM. In the analysis of the liquid layer in close contact melting, were analyzed descending velocity of the solid, liquid layer thickness and the flow in the liquid layer. The results regarding the liquid layer denote significant influence of the configurations for Darcy STM. However, such results are little influenced by the configurations of VVM.