Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-04-30T17:11:57Z No. of bitstreams: 1 MORGANA DE VASCONCELLOS ARAÚJO - DISSERTAÇÃO PPGEQ 2014..pdf: 6485490 bytes, checksum: ac573f517948e0ae4fa9b0adcf3f0f41 (MD5) / Made available in DSpace on 2018-04-30T17:11:57Z (GMT). No. of bitstreams: 1 MORGANA DE VASCONCELLOS ARAÚJO - DISSERTAÇÃO PPGEQ 2014..pdf: 6485490 bytes, checksum: ac573f517948e0ae4fa9b0adcf3f0f41 (MD5) Previous issue date: 2014-07-24 / CNPq / Capes / Diante do número de acidentes com dutos na indústria petrolífera, diversas pesquisas vêm buscando o aperfeiçoamento das tecnologias de detecção de vazamento, fazendo-se necessário o conhecimento do comportamento do(s) fluido(s) no interior da tubulação. O trabalho tem como objetivo geral o estudo termofluidodinâmico do escoamento água-óleo em uma conexão tê na presença de vazamento. Para isto, foi construído um domínio tridimensional e foi simulado o escoamento não isotérmico de água testando-se os modelos de turbulência k-Ɛ padrão, SST e RNG k-Ɛ e utilizando as equações da conservação de massa, momento e energia. Foram simulados os escoamentos bifásicos de água-óleo em uma conexão tê com vazamento, variando-se a fração volumétrica de óleo e as temperaturas das misturas nas seções de entrada. Comparando os perfis de velocidade e temperatura com dados experimentais observou-se que o modelo k-Ɛ padrão representou satisfatoriamente o comportamento do escoamento. Notou-se que quanto menor a fração volumétrica de óleo no escoamento, maior a pressão a montante do vazamento e que a região de encontro entre o duto principal e o duto secundário apresentou os maiores gradientes de pressão. / From the large number of accidents involving pipelines in the oil industry, many researcher have the objective of improving the tecnologies of leak detection, where the knowledge of the behavior (s) fluid (s) inside the pipe is necessary. The work aims to describe the thermofluid dynamics study of oil-water flow in a pipe with tee junction and with the presence of leakage. For this, a three-dimensional domain was built and it was simulated the nonisothermal flow of water by testing models of turbulence k-Ɛ standard, SST e RNG k-Ɛ, using the equations of conservation of mass, momentum and energy. Comparing the profiles of velocity and temperature with experimental data, it were observed that the k-Ɛ standard model well represented the flow pattern behavior satisfactorily. The water-oil two-phase flow in a tee junction with leakage was simulated, varying the volume fraction of oil mixtures and temperatures of the input sections. It was noted that the smaller the volume fraction of oil in the flow, the greater the pressure upstream of the leak and the region where the main duct is linked to the secondary duct showed higher pressure gradients.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The present work is dedicated to perform the mathematical modeling and DNS and LES simulations of a three-dimensional, temporally evolving incompressible plane wake are performed, seeking to evidence differences in stability, transition and onset of both coherent and small scale structures, when the flow is subjected to random perturbations of different amplitudes. The perturbations are generated using the Random-Flow-Generator (RFG) technique, being imposed to the flow as initial conditions. The Navier-Stokes equations are solved in a prismatic domain, with periodic boundary conditions in all directions, using Fourier pseudospectral method. The invariants of the velocity gradient tensor, Q and R, are analyzed for random perturbations with magnitudes 10−3, 10−4 and 10−5, showing the onset of their characteristic teardrop correlation map. Moreover, maps of the second and third invariants of the rate-of-strain tensor, QS and RS, are shown, in order to evidence the differences in local flow strain and topological characteristics of the dissipation of kinetic energy. Isosurface plots of Q and QW, as well as vorticity contours are shown, allowing visual identification of the coherent structures and confirming patterns predicted by the invariant maps. / O presente trabalho se dedica a modelagem matemática e a simulações numéricas DNS e LES de uma esteira tridimensional, incompressível, em desenvolvimento temporal, buscando evidenciar diferenças na estabilidade, transição e no desenvolvimento de estruturas coerentes e de pequena escala, quando o escoamento é submetido a perturbações randômicas de diferentes amplitudes. As perturbações são geradas utilizando-se a técnica Random Flow Generator (RFG), sendo sobrepostas à condição inicial do escoamento. As equações de Navier-Stokes são resolvidas em um domínio prismático, com condições de contorno periódicas em todas as direções, utilizando-se o método pseudoespectral de Fourier. Os invariantes do tensor gradiente de velocidade, Q e R, são analisados para perturbações de magnitude 10−3, 10−4 and 10−5, mostrando a formação de uma correlação no formato de gota, característica da resolução das equações de Navier-Stokes. Além disso, são apresentados mapas do segundo e terceiro invariante do tensor taxa de deformação, QS e RS, a fim de evidenciar as diferenças locais no escoamento e as características topológicas na taxa de dissipação de energia cinética. Isosuperfícies de Q e QW, bem como contornos de vorticidade são apresentados, possibilitando a identificação visual das estruturas coerentes, e confirmando os padrões de estruturas previstos pelos mapas de invariância. / Mestre em Engenharia Mecânica

The numerical simulation of fluid flow involving complex geometries is greatly limited by the required spatial grid resolution. These flows often contain small regions with complex motions, while the remaining flow is relatively smooth. Adaptive mesh refinement (AMR) enables the spatial grid to be refined in local regions that require finer grids to resolve the flow. This work describes an approach to parallelization of a structured adaptive mesh refinement (SAMR) algorithm. This type of methodology is based on locally refined grids superimposed on coarser grids to achieve the desired resolution in numerical simulations. Parallel implementations of SAMR methods offer the potential for accurate simulations of high complexity fluid flows. However, they present interesting challenges in dynamic resource allocation, data-distribution and load-balancing. The overall efficiency of parallel SAMR applications is limited by the ability to partition the underlying grid hierarchies at run-time to expose all inherent parallelism, minimize communication and synchronization overheads, and balance load. The methodology is based on a message passing interface model (MPI) using the recursive coordinate bisection (RCB) for domain partition. For this work, a semi-implicit projection method has been implemented to solve the incompressible Navier Stokes equations. All numerical implementations are an extension of a sequential Fortran 90 code, called "AMR3D", developed in the work of Nós (2007) .The efficiency and robustness of the applied methodology are verified via convergence analysis using the method of manufactured solutions. Validations were performed by simulating an incompressible jet flow and a lid driven cavity flow. / A simulação numérica de escoamentos envolvendo geometrias complexas é fortemente limitada pela resolução da malha espacial. Na grande maioria dos escoamentos, há pequenas regiões do domínio onde o fluido se movimenta de forma complexa gerando gradientes elevados, enquanto que no restante do domínio o escoamento é relativamente calmo". O Refinamento Adaptativo de Malhas (Adaptive Mesh Refinement - AMR), possibilita que o refinamento da malha espacial seja mais apurado em regiões especificas, enquanto que nas demais regiões o refinamento pode ser mais grosseiro. O presente trabalho consiste no desenvolvimento de uma metodologia de paralelização para a solução das equações de Navier-Stokes em malhas adaptativas bloco-estruturadas (Structured Adaptive Mesh Refinement - SAMR) utilizando a interface MPI (Message Passing Interface) e o método de bisseção por coordenadas RCB (Recursive Coordinate Bisection) para o balanço de carga. Implementações de métodos SAMR em processamento paralelo oferecem a possibilidade de simulações precisas de escoamentos de elevada complexidade. No entanto, apresentam desafios interessantes quanto à dinamicidade na alocação e distribuição dos dados e no balanceamento de carga. Cabe ressaltar que a é ciência total das aplicações envolvendo métodos SAMR em processamento paralelo é fortemente dependente da qualidade do particionamento dinâmico de domínio, efetuado em tempo de execução, para que se garanta os menores custos de comunicação e sincronização possíveis, além de uma boa distribuição da carga computacional. Neste trabalho, utilizou-se o esquema semi-implícito proposto por Ceniceros et al. (2010) para avanço temporal. Todas as implementações foram efetuadas como uma extensão do código AMR3D", proposto por Nós (2007). A é ciência e a robustez do método proposto são verificadas por meio do método das soluções manufaturadas. As validações foram feitas por meio da simulação do escoamento em uma cavidade com tampa deslizante e de um jato incompressível. / Doutor em Engenharia Mecânica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / As técnicas de injeção de traçadores têm sido amplamente utilizadas na investigação de escoamentos em meios porosos, principalmente em problemas envolvendo a simulação numérica de escoamentos miscíveis em reservatórios de petróleo e o transporte de contaminantes em aquíferos. Reservatórios subterrâneos são em geral heterogêneos e podem apresentar variações significativas das suas propriedades em várias escalas de comprimento. Estas variações espaciais são incorporadas às equações que governam o escoamento no interior do meio poroso por meio de campos aleatórios. Estes campos podem prover uma descrição das heterogeneidades da formação subterrânea nos casos onde o conhecimento geológico não fornece o detalhamento necessário para a predição determinística do escoamento através do meio poroso. Nesta tese é empregado um modelo lognormal para o campo de permeabilidades a fim de reproduzir-se a distribuição de permeabilidades do meio real, e a geração numérica destes campos aleatórios é feita pelo método da Soma Sucessiva de Campos Gaussianos Independentes (SSCGI). O objetivo principal deste trabalho é o estudo da quantificação de incertezas para o problema inverso do transporte de um traçador em um meio poroso heterogêneo empregando uma abordagem Bayesiana para a atualização dos campos de permeabilidades, baseada na medição dos valores da concentração espacial do traçador em tempos específicos. Um método do tipo Markov Chain Monte Carlo a dois estágios é utilizado na amostragem da distribuição de probabilidade a posteriori e a cadeia de Markov é construída a partir da reconstrução aleatória dos campos de permeabilidades. Na resolução do problema de pressão-velocidade que governa o escoamento empregase um método do tipo Elementos Finitos Mistos adequado para o cálculo acurado dos fluxos em campos de permeabilidades heterogêneos e uma abordagem Lagrangiana, o método Forward Integral Tracking (FIT), é utilizada na simulação numérica do problema do transporte do traçador. Resultados numéricos são obtidos e apresentados para um conjunto de realizações amostrais dos campos de permeabilidades. / Tracer injection techniques have been widely used to investigate flows in heterogeneous porous media, especially in problems related to numerical simulation of miscible flows in oil reservoirs and to contaminant transport in aquifers. Oil reservoirs are generally heterogeneous and may possess spatially significant variations in their properties on several length scales. These spatial variations are incorporated into the governing equations for flow problems in porous media on the basis of random fields. Random fields provide a natural description of rock heterogeneities in the typical case in which the geological knowledge of rock is much less detailed than is necessary to predict flow properties through it deterministically. In this thesis we adopt a scalar log-normal permeability field k(x) to reproduce the statistical distribution of the permeability values of a real medium, and the numerical generation of these random fields is based on a Successive Sum of Independent Gaussian Fields defined on multiple length scales. The aim of this work is to study the uncertainty quantification in inverse problems for tracer transport in heterogeneous porous media in a Bayesian framework and propose the permeability update based on observed measurements of spatially sparse tracer concentration at certain times. A two-stage Markov chain Monte Carlo (MCMC) method is used to sample posterior probability distribution with hierarchical priors and the Markov chain is constructed from random reconstruction of the permeability fields. To solve the Darcys law we use a mixed finite elements method which are suitable to compute accurately the relevant fluxes in heterogeneous permeability fields and a Lagrangian strategy, the Forward Integral Tracking (FIT) method, for the numerical simulation of tracer transport problem. Numerical results are presented for a set of sampled realizations of the permeability fields.

A interação entre fluidos e estruturas caracteriza um problema multi-físico não linear e está presente numa grande variedade de áreas da engenharia. Este trabalho apresenta o desenvolvi mento de ferramentas computacionais com base no Método dos Elementos Finitos (MEF) para a análise de interação fluido-estrutura (IFE) considerando escoamentos com baixas velocidades. Dada a interdisciplinaridade do tema, se faz necessário o estudo em três diferentes assuntos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional, e o problema de acoplamento. No caso da dinâmica das estruturas empregar-se um elemento finito que seja adequado para a simulação de problemas de IFE, que claramente demandam uma análise não linear geométrica, optando-se pelo emprego de uma formulação descrita em posições, a qual evita problemas relativos à aproximação de rotações finitas. Quanto à dinâmica dos fluidos computacional, é empregado um método estável e ao mesmo tempo sensível à movimentação da estrutura, utilizando a descrição Lagrangeana-Euleriana Arbitrária (ALE). Os casos considerados neste trabalho, assim como muitos dos problemas de engenharia, ocorrem com escoamentos em baixas velocidades, implicando na incompressibilidade do fluido, o que demanda, para um método estável, a utilização de elementos que atendam à condição de Ladyzhenskaya-Babuska-Brezzi (LBB). Além disso, é necessário também o emprego de métodos que consigam neutralizar as variações espúrias decorrentes da não-linearidade de possíveis escoamentos com convecção dominante e que surgem com a aplicação do processo clássico de Galerkin. Para superar esse problema, é aplicado o método Streamline-Upwind/Petrov-Galerkin (SUPG), que adiciona difusividade artificial na direção do escoamento, controlando a amplitude dos termos convectivos. No que se refere ao acoplamento fluido-casca, buscam-se modularidade e versatilidade adotando-se o modelo particionado. O modelo de acoplamento implementado garante ainda a utilização de malhas do fluido e da estrutura sem a necessidade de coincidência de nós. / Interaction between fluids and structures characterizes a nonlinear multi-physics problem presente in a wide range of engineering fields. This works presets the development of computational tools based on finite element method (FEM) for fluid-structure interaction (FSI) analysis considering low speed flows (incompressible), as a great part of the engineering problems. Given the topic multidisciplinary nature, it is necessary to study three different subjects: the computational structural dynamics, the computational fluid mechanics and the coupling problem. Regarding structural mechanics, we seek to employ a finite element adequate to FSI simulation, what clearly demands a geometric nonlinear analysis. We chose to employ shell elements with formulation in terms of positions, which avoids problems related to finite rotations approximations. Concerning computational fluid dynamics, we employ a stable method, at same time sensible o structural movements, which is written in the arbitrary Lagrangian-Eulerian (ALE) description. The flow incompressibility demands, for a stable method, the use of elements according to the Ladyzhenskaya-Bbuska-Brezzi (LBB) condition. It is also necessary to employ methods able to neutralize the spurious variations that appears from convection dominated flows when applying the standard Galerking method. In order to overcome this problem, we apply the Streamline-Upwind/Petrov-Galerkin (SUPG) method, which adds artificial diffusivity to the streamline direction, controlling spurious variations. Considering the fluid-shell coupling, we seek modularity and versatility, adopting the partitioned model. The developed coupling model ensure the use of fluid and structure meshes with no need for matching nodes.

Este trabalho tem como objetivo desenvolver um método numérico para simular escoamentos incompressíveis, isotérmicos, confinados ou com superfícies livres, de fuidos viscoelásticos governados pelos modelos integrais de Maxwell e K-BKZ (Kaye-Bernstein, Kearsley e Zapas). A técnica numérica apresentada é uma extensão do método GENSMAC (Tomé McKee - J. Comp. Phys., (110), pp 171--186, 1994 ) para a solução das equações de conservação, juntamente com as equações constitutivas integrais de Maxwell e K-BKZ. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada. O tensor de Finger, B_t\'(t) é calculado com base nas idéias do método de campos de deformação (Peters et al. - J. Non-Newtonian Fluid Mech. (89), de maneira que não há a necessidade de seguir a trajetória da partícula de fuido para descrever a história de deformação da partícula. Uma abordagem diferente para a discretização do instante passado é utilizada e o tensor de Finger e o tensor das tensões são calculados utilizando um método de segunda ordem. A validação do método numérico descrito nesse trabalho foi feita utilizando o escoamento em um canal bidimensional e a solução numérica obtida para a velocidade e para as componentes de tensão com o modelo de Maxwell foram comparadas com as respectivas soluções analíticas no estado estacionário, mostrando excelente concordância. Os resultados numéricos para a simulação do escoamento em uma contração planar 4 : 1 mostraram bons resultados, tanto qualitativos quanto quantitativos, quando comparados com os resultados experimentais de Quinzani et al. ( J. Non-Newtonian Fluid Mech. (52), pp 1?36, 1994 ). Além disso, utilizando os modelos Maxwel e K-BKZ, o escoamento em uma contração planar 4 : 1 foi simulado para vários números de Weissenberg e os resultados obtidos estão de acordo com os encontrados na literatura. Resultados numéricos de escoamentos com superfícies livres modelados pelas equações integrais de Maxwell e K-BKZ são apresentados. Em particular, a simulação numérica do jato oscilante para diferentes números de Weissenberg e diferentes números de Reynolds é apresentada. / The aim of this work is to develop a numerical technique for simulating incompressible, isothermal, free surface (also con¯ned) viscoelastic flows of fuids governed by the integral models of Maxwell and K-BKZ (Kaye-Bernstein, Kearsley and Zapas). The numerical technique described herein is an extension of the GENSMAC method (Tome and McKee, J. Comput. Phys., 110, pp. 171-186, 1994) to the solution of the momentuum and mass conservation equations together with the integral constitutive Maxwell and K-BKZ equations. The governing equations are solved by the finite difference method on a staggered grid using a Marker-and-Cell approach. The fluid is represented by marker particles on the fluid surface only. This provides the visualization and location of the fluid free surface so that the free surface stress conditions can be applied. The Finger tensor Bt0(t) is computed using the ideias of the deformation fields method (Peters et al. J. Non-Newtonian Fluid Mech., 89, pp. 209-228, 2001) so that it is not necessary to track a fluid particle in order to calculate its deformation history. However, in this work modifcations to the deformation fields method are introduced: the past time is discretized using a different formula, the Finger tensor Bt0(x; t) is obtained by a second order method and the stress tensor ? (x; t) is computed by a second order quadrature formula. The numerical method presented in this work is validated by simulating the flow of a Maxwell fluid in a two-dimensional channel and the numerical solutions of the velocity and the stress components are compared with the respective analytic solutions providing a good agreement. Further, the flow through a 4:1 planar contraction of a specific fuid studied experimentally by Quinzani et al. (J. Non-Newtonian Fluid Mech., 52, pp. 1-36, 1994) was simulated and the numerical results were compared qualitatively and quantitatively with the experimental results and very good agreement was obtained. The Maxwell and the K-BKZ models were applied to simulate the 4:1 planar contraction problem using various Weissenberg numbers and the numerical results were in agreement with those published in the literature. Finally, numerical results of free surface flows using the Maxwell and K-BKZ integral constitutive equations are presented. In particular, the numerical simulation of jet buckling using several Weissenberg numbers and various Reynolds numbers are presented

Nesta tese é apresentado um esquema novo de alta resolução upwind (denominado TDPUS-C3) para reconstrução de fluxos numéricos para leis de conservação não lineares e problemas relacionados em DFC. O esquema é baseado nos critérios de estabilidade CBC e TVD e desenvolvido utilizando condições de diferenciabilidade \'C POT. 3\'. Além disso, é realiozada a implementação da associação do esquema TDPLUS-C3 com a modelagem de turbulência RNG \'\\kappa - \\epsilon\'. O propósito é obter soluções numéricas de sistemas hiperbólicos de leis de conservação para dinâmica dos gases e equações de Navier-Stokes para escoamento incompreensível de fluidos newtonianos e não newtonianos (viscoelásticos). Fazendo o uso do esquema TDPUS-C3, a precisão global dos métodos numéricos é verificada acessando o erro em problemas teste (benchmark) 1D e 2D. Um estudo comparativo entre os resultados do esquema TDPUS-C3 e os esquemas upwind convencionais para leis de conservação hiperbólicas complexas é também realizado. A Associação das modelagens numéricas (upwinding mais RNG \'\\kappa - \\epsilon\') é , então, examinada na simulação de escoamentos turbulentos de fluidos newtonianos envolvendo superfícies livres móveis, usando a metodologia URANS. No geral, em termos do comportamento global, concordância satisfatória é observada / In this thesis, a new high-resolution upwind scheme (named TDPUS-C3) for reconstruction of numerical fluxes for nonlinear conservation laws and related CFD problems in presented. The scheme is based on CBC and TVD stability criteria and developed by employing differentiability condictions (\'C POT. 3\'). In additon, the implementation of an association of the TDPUS-C3 scheme with the RNG \'\\kappa - \\epsilon\' turbulence modelling is also performed. The purpose is to obtain numerical solutions of systems of hyperbolic conservation laws for gas dynamics and Navier-Stokes equations for incompressible flow of Newtonian and non-Newtonian (viscoelstic) fluids. By using the TDPUS-C3 scheme, the global accuracy of the numerical methods is verified by assessing the error on 1D and 2D benchmark test cases. A comparative study between the TDPUS-C3 scheme and convectional upwind schemes to solve standard and complex hyperbolic conservation laws is also accomplished. The association of the numerical modelling (upwinding plus RNG \'\\kappa - epsilon\') is then examined in the simulation of turbulent Newtonian fluid flows involving moving free surfaces, by using URANS methodology. Overall, satisfactory agreement is found in terms of the overall behaviour

A presente pesquisa tem como objetivo realizar um estudo teórico e desenvolver simulações computacionais envolvendo a dinâmica de sistemas gás-sólido. O foco principal do trabalho é a modelagem do particulado através da análise das forças de contato entre partículas de materiais granulares utilizando modelos contínuos baseados na mecânica dos solos e na teoria cinética dos escoamentos granulares (sistemas grandes com muitas partículas, formulação Euleriana - Volumes Finitos) e modelos discretos baseados nas características físicas dos materiais (sistemas intermediários e número limitado de partículas, formulação Lagrangeana - Método dos Elementos Discretos). Investigam-se os modelos existentes na literatura com intuito de melhorar os modelos contínuos e discretos baseados na interação entre as partículas que caracterizam a dinâmica do particulado em sistemas gás-sólido. Propõe-se uma nova abordagem para a determinação do coeficiente de rigidez da mola baseada em uma equivalência entre os modelos lineares e não-lineares. Utiliza-se o código fonte MFIX para realizar simulações computacionais da dinâmica de sistemas gás-sólido, analisando o processo de fluidização, mistura e segregação de partículas, influência das correlações de arrasto, e análise das forças de contato entre as partículas através do novo método para a determinação do coeficiente de rigidez da mola . Os resultados obtidos são comparados com dados numéricos e experimentais da literatura. / The purpose of the present study is to perform a theoretical study and develop numerical simulations involving dynamic in gas-solid systems. The focus of the work is the modeling of particulate matter using continuous models based on soil mechanics and the kinetic theory of granular flows (large systems with many particles, Eulerian formulation - Finite Volume) and discrete models based on physical characteristics of the particles (intermediate systems and limited number of particles, Lagrangian formulation - Discrete Element Method). It is proposed a new approach to determine the normal spring stiffness coefficient of the linear model through the numerical solution for the overlap between particles in non-linear models. The linear spring stiffness is determined using an equivalence between the linear and the non-linear models. It is used the MFIX computational code to perform numerical simulations of the dynamics of gas-solid systems. It is analyzed the processes of fluidization, mixing and particle segregation and the influence of drag correlations. The proposed approach for normal spring stiffness coefficient is applied in the numerical simulations of two problems: single freely falling particle and bubbling fluidized bed. The results were compared with numerical and experimental data from literature.

Foram desenvolvidos estudos de modelagem e simulação numérica de escoamentos bifásicos gás-sólido na coluna ascendente de leitos fluidizados circulantes utilizando um modelo Euleriano de duas fases separadas. O sistema de equações diferenciais parciais conservativas governantes foi obtido através de um procedimento tradicional. Ambas as fases foram assumidas como meio contínuo. Aplicou-se o procedimento de médias estatísticas de Euler, enfatizando a obtenção dos modelos hidrodinâmicos A e B desenvolvidos no IIT/ANL. Realizou-se análise comparativa de correlações para transferência de quantidade de movimento na interface. Discutiu-se a formulação de condições de contorno apropriadas. As equações diferenciais parciais médias foram discretizadas em volumes de controle Eulerianos. As equações de continuidade foram resolvidas implicitamente. As equações de quantidade de movimento foram resolvidas através de um procedimento explícito-implícito. Foram desenvolvidas simulações numéricas para uma coluna ascendente típica de leitos fluidizados circulantes. Desenvolveu-se análise paramétrica da influência de vários aspectos físicos e matemáticos sobre o escoamento. Avaliou-se resultados de simulação através de metodologia de identificação e caracterização de estruturas coerentes. Estudou-se o efeito da função de arrasto na interface sobre os processos dinâmicos que caracterizam estas estruturas coerentes. Foram realizados estudos numéricos de turbulência a partir de resultados de simulação direta. Várias conclusões e recomendações para futuros trabalhos foram propostas com base nas análises realizadas. Foram apresentadas algumas considerações gerais relativas a aspectos críticos na modelagem e simulação com modelo das duas fases separadas. / Studies were carried out on modeling and numerical simulation of gas-solid two-phase flows in the riser of circulating fluidized beds using an Eulerian two-fluids model. The system of conservative partial differential governing equations was derived through a traditional procedure. Both phases were assumed as a continuum. The Euler averaging procedure was applied emphasizing the derivation of the so called hydrodynamic models A and B developed at IIT/ANL. A comparative analysis was performed among correlations for momentum transfer at the interface. The formulation of suitable boundary conditions was discussed. The average partial differential conservative equations were discretized on Eulerian control volumes. The continuity equations were solved implicitly. The momentum equations were solved through an explicit-implicit procedure. Numerical simulation was performed for a typical circulating fluidized bed riser. A parametric analysis was carried out regarding the influence on the flow of various physical and mathematical aspects. Results of simulation were evaluated through a methodology of identification and characterization of coherent structures. The effect of the interface drag function on dynamic features of those coherent structures was addressed. Numerical studies on turbulence were performed from results of direct simulation. Several conclusions and recommendations for future work were put forward on the basis of the performed analyses. Some general considerations were presented regarding critical features of modeling and simulation through Eulerian two-fluids models.

Este trabalho é dedicado à análise e implementação de esquemas "upwind" de alta ordem modernos e o modelo de turbulência k-epsilon padrão no Freeflow-2D; um ambiente integrado para simulação numérica em diferenças finitas de problemas de escoamentos incompressíveis com superfícies livres. O propósito do estudo é a simulação de escoamentos de fluidos newtonianos incompressíveis, bidimensionais, confinados e/ou com superfícies livres e a altos valores do número de Reynolds. O desempenho do código Freeflow-2D atual é avaliada na simulação do escoamento numa expansão brusca e de um jato livre incidindo perpendicularmente sobre uma superfície rígida impermeável. O código é então aplicado na simulação de um jato planar turbulento em uma porção de fluido com superfície livre e estacionário. Os resultados numéricos obtidos são comparados com dados experimentais, soluções analíticas e soluções numéricas de outros trabalhos. / This work is devoted to the analysis and implementation of modern high-order upwind schemes and the standard k-epsilon turbulence model into the Freeflow-2D; a finite difference integrated environment for the numerical simulation of incompressible free surface flow problems. The purpose of this study is the two-dimensional simulation of high-Reynolds incompressible newtonian confined and/or free surface flows. The performance of the current Freeflow-2D code is assessed by applying it to the simulation of flow over a backward facing step and of an impinging free jet onto an impermeable rigid surface. The code is then applied to a turbulent planar jet into a pool. The numerical results are compared with experimental data, analytical solution, and numerical simulations of other works.