O aterramento elétrico cumpre uma função essencial em um sistema de distribuição de energia elétrica, devendo atender, dentre vários requisitos, aqueles estabelecidos quanto aos valores de tensão de passo e de toque. Entretanto, em locais onde o solo apresenta valores elevados de resistividade, surgem dificuldades que nem sempre são superadas no sentido de obedecer às exigências citadas nas normas técnicas e especificações. Dentro deste contexto, este trabalho tem por objetivo avaliar a utilização de condutores e hastes envolvidos em concreto em sistemas de aterramento, visando à redução dos valores de resistência de terra. Inicialmente é abordada uma configuração convencional de malha de aterramento adotada por uma concessionária de energia elétrica, sendo realizados estudos de modo a verificar a influência do espaçamento e da quantidade das hastes nos valores de resistência de terra. São também consideradas modificações nas dimensões da malha e a utilização de hastes profundas e, na seqüência, são analisados os resultados de outras configurações. De modo a verificar a efetividade do uso do concreto, foram realizadas várias simulações computacionais considerando o sistema de aterramento com a utilização de condutores e hastes envolvidos em concreto. Diversas situações foram investigadas, sendo possível concluir que o uso do concreto resulta em menores valores de resistência de terra em relação às malhas convencionais. Esse resultado, de certo modo, vem ao encontro das expectativas, porém deve ser ressaltada a contribuição do trabalho, a qual possibilita a análise, através de cálculos teóricos, do comportamento das malhas de aterramento com o uso do concreto. A literatura pesquisada apresenta poucos resultados teóricos, sendo encontradas algumas investigações nesse tema onde são consideradas somente medições realizadas no campo. / The electrical grounding plays an essential role in a system of electric power distribution, and must follow, among many requirements, those ones established as step voltage values and touch voltage values. However, in those places where the soil has high resistivity values, some difficulties arise and may not always be overcomed in order to follow the requirements set forth in technical standards and specifications. In this context, the aim of this study is to evaluate the use of electrical conductors and rods encased in concrete inside grounding systems, in order to reduce the ground resistance values. Initially it was observed a conventional configuration of grounding grid that is used by an electric utility, and studies were made in order to verify the influence of the number of rods and the distance between them in the values of grounding resistance. There are also been taken into consideration the changes in the grounding grid dimensions and the use of deep rods and, right after, the results of other configurations are analyzed. In order to verify the effectiveness of the use of concrete, several computer simulations were performed considering the grounding system that used conductors and rods encased in concrete. Several situations were investigated, and it was possible to conclude that the use of concrete results in lower values of ground resistance compared to conventional grounding grid. This result, in a certain way, matches with the expectations, but it should be emphasized the contribution of the study which enables, through theoretical calculations, the analysis of the grounding grid behavior with the use of concrete. The literature shows few theoretical results and presents only some researches on this issue in which ones only measurements made in the field are considered.

Meshless methods provide a highly continuous approximation field, convenient for thin structures like shells. Nevertheless, the lack of Kronecker Delta property makes the formulation of essential boundary conditions not straightforward, as the trial and test fields cannot be tailored to boundary values. Similar problem arise when different approximation regions must be joined, in a multi-region problem, such as kinks, folds or joints. This work presents three approaches to impose both kinematic conditions: the well known Lagrange Multiplier method, used since the beginning of the Element Free Galerkin method; a pure penalty approach; and the recently rediscovered alternative of Nitsche\'s Method. We use the EFG discretization technique for thick Reissner-Mindlin shells and adapt the weak form as to separate displacement and rotational degrees of freedom and obtain suitable and separate stabilization parameters. This approach enables the modeling of discontinuous shells and local refinement on multi-region problems. / Métodos sem malha geram campos de aproximação com alta continuidade, convenientes para estruturas finas como cascas. No entanto, a ausência da propriedade de Delta de Kronecker dificulta a formulação de condições essenciais de contorno, já que os campos de aproximação e teste não podem ser moldados aos valores de contorno. Um problema similar aparece quando diferentes regiões de aproximação precisam ser juntadas em um problema multi-regiões como dobras, vincos ou junções. Este trabalho apresenta três métodos de imposição ambas condições cinemáticas: o já conhecido método dos multiplicadores de Lagrange, usado desde o começo do método de Galekin sem elementos (EFG); uma abordagem de penalidade pura; e o recentemente redescoberto método de Nitsche. Nós usamos a técnica de discretização com EFG para cascas espessas de Reissner-Mindlin e adaptamos a forma fraca de forma a separar graus de liberdade de deslocamento e rotação e obter coeficientes de estabilização diferentes e apropriados. Essa abordagem permite a modelagem de cascas discontínuas e o refinamento local em problemas multi-regiões.

Nas últimas três décadas os phase locked loops (PLLs) totalmente digitais têm recebido muita atenção devido, principalmente, às vantagens que eles oferecem em comparação aos PLLs analógicos. Essas vantagens incluem melhor desempenho, maior velocidade e confiabilidade, tamanho reduzido e menor custo. Os PLLs também são amplamente utilizados em sistemas de comunicações e em outras aplicações digitais. A presente dissertação é uma contribuição no campo dos PLLs digitais adaptativos e otimizados para a sua implementação em hardware. É feito uma análise de suas características dinâmicas e proposta uma nova estrutura de PLL digital capaz de melhorar a resposta da malha em termos de tempo de aquisição e largura de banda. A Malha Síncrona Digital \"Tanlock\" com Estimação de Frequência e Ganho Adaptativo para Convergência Rápida, como é chamada, foi desenvolvida a partir da malha digital \"Tanlock\", utilizando-se teoremas de ponto fixo e mapas contrativos para determinar as condições de ganho que garantam convergência rápida e melhor utilização da largura de banda. Resultados das simulações são comparados com os obtidos teoricamente para avaliar o desempenho da malha proposta. / In the last three decades, fully-digital Phase-Locked-Loops (PLLs) systems have received a lot of attention due to its advantages in comparison with analog PLLs. These advantages include improved transient response, reliability and also reduced size and cost. The PLLs are widely used in communications systems and many other digital applications. This dissertation is a contribution to the field of digital adaptive PLLs optimized to hardware implementation. Here, a new PLL structure is presented; the Frequency Sensing Adaptive TDTL is an improvement to the classic Time-Delay Tanlock structure, alowing fast convergence to the synchronous states, using fixed-point theorems and contractive maps to determine the gain conditions which ensure the rapid convergence and also providing wider bandwidth. The results of simulations are compared with those obtained theoretically in order to assess the loop performance.

O método numérico para fluidos incompressíveis desenvolvido no presente estudo é o Moving Particle Semi-Implicit Method (MPS) que enxerga o domínio discretizado em partículas, é baseado em representação lagrangeana e não tem a necessidade de utilização de malhas. O método MPS tem como equações governantes uma forma particular da equação de Navier-Stokes e a equação da continuidade para fluidos incompressíveis e não viscosos. Os métodos de simulação de fluidos mais comumente utilizados são baseados em representação euleriana e utilizam malhas para descrever a geometria do domínio a ser simulado. Devido a essas diferenças, uma das grandes virtudes do método de partículas é a facilidade de investigação de fenômenos altamente não-lineares como o de superfície livre com quebra de ondas, de líquidos no interior de uma embarcação em movimento, de ondas batendo na parte externa do casco de um navio, etc. Em artigos já publicados, resultados de experimentos físicos mostraram boa aderência aos resultados numéricos de simulações realizadas com o método MPS. No presente trabalho, resultados das forças de excitação das simulações com ondas regulares foram comparados com os resultados do programa Wave Analysis MIT (WAMIT) que é um programa consagrado no meio científico. Houve uma boa concordância de resultados entre os dois programas. A otimização do cálculo de vizinhança forneceu uma grande economia de tempo computacional. A maior contribuição deste estudo foi a otimização da função que resolve o sistema linear implementando no código desenvolvido um código paralelizado de uso público existente chamado Portable, Extensible Toolkit for Scientific Computation (PETSc) que proporcionou um bom ganho de desempenho. / A numerical method called Moving Particle Semi-implicit (MPS) method was developed in this study to analyze incompressible fluids. It is a particle method using a lagrangean representation without any grid. The governing equations are the Navier-Stokes equation and continuity equation for incompressible and non-viscous flow. Most of the computational fluid dynamics (CFD) methods are based on eulerian representation and use grids to describe the geometry of the simulated domain. These differences make the MPS method easier to analyze highly nonlinear phenomena as free surface with wave breaking, sloshing, slamming, etc. In previously published articles, results of physical experiments had shown good agreement with the numerical results obtained with MPS method. In the present work, results of exciting forces were compared with the results obtained with a validated program called Wave Analysis MIT (WAMIT). It had a good agreement of results between these two programs. The optimization of the neighborhood calculation function got a good economy of computational time. The greatest contribution of this study was the optimization of the linear system solver. It was made implementing in the developed code a parallelized public code called Portable, Extensible Toolkit for Scientific Computation (PETSc) that provided a good performance profit.

CoordenaÃÃo de AperfeiÃoamento de Pessoal de NÃvel Superior / A animaÃÃo de tecido vem sendo estudada abundantemente nas Ãltimas dÃcadas por haver demanda na indÃstria do entretenimento bem como no comÃrcio eletrÃnico de roupas. Esse tipo de animaÃÃo, na maioria das vezes, à feita com base em simulaÃÃo fÃsica, havendo muito gasto computacional. Esse trabalho tenta usufruir de vÃrios modelos para diminuir esse gasto. à feito uma modelagem do tecido com uma malha triangular e usa-se um modelo massa-mola-amortecedor para simular as forÃas entre as partÃculas, que sÃo os vÃrtices dessa malha. Dependendo da disposiÃÃo do modelo do tecido e dos objetos da cena, sÃo detectadas colisÃes entre eles. A malha à discretizada ou simplificada, levando em consideraÃÃo a compressÃo, a colisÃo e a curvatura das molas, para que seja usada uma malha boa em cada passo ao longo da animaÃÃo. Portanto, o objetivo principal desse trabalho à estudar o comportamento do tecido utilizando o modelo de remalhamento para adaptar essa malha. / In the last decades, cloth animation has been the focus of much research, due to demands from the entertainment industry and from e-commerce. That type of animation is most often the result of a physics-based simulation and has a great computational cost. This work investigates how to reduce the computational cost of the simulation, by refining the mesh only in regions that need a fine level of detail. The fabric model consists of a triangular mesh and uses a spring-mass-damper system to compute the forces among the particles, which are located at the meshâs vertices. The collision detection depends on the arrangement of the cloth model and the objects in the scene. The mesh is refined or simplified, taking into account the spring compression, collision and curvature, so the simulation uses a better mesh every time step. Therefore, this workâs main objective is to study the dynamic behavior of cloth, using a remeshing procedure in order to adapt the mesh.

Submitted by Fabio Sobreira Campos da Costa (fabio.sobreira@ufpe.br) on 2016-06-14T12:31:20Z No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Dissertacao Saulo Menezes.pdf: 3142876 bytes, checksum: 87a4e508de7bbbb2df2129de180c6f4d (MD5) / Made available in DSpace on 2016-06-14T12:31:20Z (GMT). No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Dissertacao Saulo Menezes.pdf: 3142876 bytes, checksum: 87a4e508de7bbbb2df2129de180c6f4d (MD5) Previous issue date: 2014-08-07 / CAPEs / Neste trabalho, abordamos o processo de adaptação de malhas via “remeshing”, a adap-tação de malhas não estruturadas é usada aqui como uma ferramenta de alto desempe-nho que visa simultaneamente garantir a qualidade do resultado e reduzir o tempo com-putacional envolvido em toda a simulação. Usamos um simulador numérico de escoa-mento óleo-água em reservatório de petróleo em duas dimensões. Os fluidos e a rocha reservatório são tidos como incompressíveis e não se considera variações térmicas. Um método de volumes finitos centrado nos nós e baseado em uma estrutura de dados por arestas é empregado na discretização das equações de pressão e de saturação. Uma for-mulação agregada IMPES (Implicit Pressure Explicit Saturation) é utilizada para resol-ver as variáveis de saturação e de pressão. Aproximações de baixa e alta ordem na eq. de saturação e seus efeitos na propagação dos erros é calculada. A cada avanço no tem-po, novos campos de saturações e de pressões são avaliados e, em seguida, a qualidade destes resultados é estimada. Isto é realizado por meio de um estimador de erros basea-do na hessiana dos campos de pressão e de saturação, os quais são usados para calcular os erros local e global referentes a cada campo. Estes erros são comparados a uma tole-rância que irá definir se a malha deve ou não ser adaptada, onde e em que grau de refi-namento. O método de adaptação desenvolvido é o do “remeshing” ou redefinição de malhas, que se baseia na reconstrução total ou parcial da malha. A redefinição da malha é feita com o gerador open-source Gmsh que possibilita a criação da malha com contro-le do grau de refinamento local através da malha de “background”. Após a adaptação, os campos de pressão e de saturações são interpolados da malha anterior para a nova malha. Métodos lineares, quadráticos e adaptativos de interpolação são explorados e avaliados. Para um mesmo instante de tempo, o processo de adaptação se repete até que a qualidade dos resultados atinja a tolerância exigida. Para a integração e conformidade de todas estas etapas foi desenvolvido um software na linguagem C++ usando a biblio-teca para gerenciamento de malhas FMDB (Flexible Distributed Mesh Data Base). Foi efetuado um estudo da qualidade dos resultados obtidos e da eficiência da simulação quando da utilização do procedimento de adaptação de malhas desenvolvido em simula-ções de escoamento monofásico e bifásico óleo-água em meios porosos usando malhas não estruturadas 2D triangulares, considerando meios homogêneos e heterogêneos, iso-trópicos e anisotrópicos, comprovando a robustez da metodologia implementada nas simulações efetuadas. / In this work, we address the mesh adaptation process using "remeshing" tehcnique. The adaptation of unstructured meshes is used here as a high performance tool intended both to ensure the quality of the result and to reduce the computational time involved in the whole simulation. The numerical simulation are performed using a two dimensional two phase flow numerical simulation of oil and water in porous media. Fluids and reservoir rock are considered incompressible and does not consider thermal variations are consid-ered. A cell centered finite volume method using an edge based data structure used for of the error the discretization of the pressure and saturation equations. The IMPES for-mulation (Implicit Pressure Explicit Saturation) composes the structure of the simula-tor.The hyperbollic saturation equation is discretized using either a first or higher order aproximation and the effect of such choice is analysed. After advancing the solution in time the obtained result is analysed using an a-posteriori erro indicator. At each time step a new field of saturation and pressure is evaluated, and then, the quality of results is assessed. This is accomplished by means of an estimator based on the Hessian of the fields of pressure and saturation, which are used to calculate both the local and global errors for each field. These errors are compared to a tolerance that will define whether or not the mesh must be adapted where the mesh density required. The adaptive method developed is refered to as "remeshing" and considers either total or partial reconstruc-tion of the mesh. The remeshing is done with the open-source mesh generator named Gmsh that enables to build mesh containing the level of refinement through the mesh refinement "background." mesh. After adaptation, the fields of pressure and saturations are interpolated from the previous mesh to the new mesh. Linear and quadratic interpo-lation methods are explored and evaluated. For the same level of time, the adaptation process is repeated until the quality of results reached the required tolerance. For the integration and accomplishement of all these steps a software was developed in C + + language and using the library for unstructured mesh management FMDB (Flexible Distributed Mesh Data Base) mesh. A study of the quality of the analysed results and of the efficiency of the simulation was performed when using the adaptive mesh procedure developed through the analysis of single-phase flow and two-phase oil-water in porous media using 2D models with triangular unstructured meshes. The analysed examples considers homogeneous and heterogeneous media, isotropic and anisotropic, demon-strating the robustness and efficiency of the methodology implemented for the cases analysed in the simulations performed.

O Método dos Elementos Finitos Generalizados, MEFG, compartilha importantes características dos métodos sem malha. As funções de aproximação do MEFG, atreladas aos pontos nodais, são enriquecidas de modo análogo ao refinamento p realizado no Método das Nuvens hp. Por outro lado, por empregar uma malha de elementos para construir as funções partição da unidade, ele também pode ser entendido como uma forma não convencional do Método dos Elementos Finitos. Neste trabalho, ambas as interpretações são consideradas. Os métodos sem malha, particularmente o Método de Galerkin Livre de Elementos e o Método das Nuvens hp, são introduzidos com o propósito de estabelecer os conceitos fundamentais para a descrição do MEFG. Na seqüência, apresentam-se aplicações numéricas em análise linear e evidenciam-se características que tornam o MEFG interessante para a simulação da propagação de descontinuidades. Após discutir os modelos de dano adotados para representar o comportamento não-linear do material, são introduzidos exemplos de aplicação, inicialmente do Método das Nuvens hp e depois do MEFG, na análise de estruturas de concreto. Os resultados obtidos servem de argumento para a implementação de um procedimento p-adaptativo, particularmente com o MEFG. Propõe-se, então a adaptação do Método dos Resíduos em Elementos Equilibrados à formulação do MEFG. Com vistas ao seu emprego em problemas não-lineares, algumas modificações são introduzidas à formulação do estimador. Mostra-se que a medida obtida para representar o erro, apesar de fundamentada em diversas hipóteses nem sempre possíveis de serem satisfeitas, ainda assim viabiliza a análise não-linear p-adaptativa. Ao final, são enumeradas propostas para a aplicação do MEFG em problemas caracterizados pela propagação de defeitos / The Generalized Finite Element Method, GFEM, shares several features with the so called meshless methods. The approximation functions used in the GFEM are associated with nodal points like in meshless methods. In addition, the enrichment of the approximation spaces can be done in the same fashion as in the meshless hp-Cloud method. On the other hand, the partition of unity used in the GFEM is provided by Lagrangian finite element shape functions. Therefore, this method can also be understood as a variation of the Finite Element Method. Indeed, both interpretations of the GFEM are valid and give unique insights into the method. The meshless character of the GFEM justified the investigation of meshless methods in this work. Among them, the Element Free Galerkin Method and the hp-Cloud Method are described aiming to introduce key concepts of the GFEM formulation. Following that, several linear problems are solved using these three methods. Such linear analysis demonstrates several features of the GFEM and its suitability to simulate propagating discontinuities. Next, damage models employed to model the nonlinear behavior of concrete structures are discussed and numerical analysis using the hp-Cloud Method and the GFEM are presented. The results motivate the implementation of a p-adaptive procedure tailored to the GFEM. The technique adopted is the Equilibrated Element Residual Method. The estimator is modified to take into account nonlinear peculiarities of the problems considered. The hypotheses assumed in the definition of the error measure are sometimes violated. Nonetheless, it is shown that the proposed error indicator is effective for the class of p-adaptive nonlinear analysis investigated. Finally, several suggestions are enumerated considering future applications of the GFEM, specially for the simulation of damage and crack propagation

Submitted by Franciele Moreira (francielemoreyra@gmail.com) on 2017-09-06T17:43:44Z No. of bitstreams: 2 Dissertação - Rafael Ribeiro de Carvalho Vaz - 2017.pdf: 6167014 bytes, checksum: 8f2f0dbd609aac9c2bd59458c6d0c243 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2017-09-15T15:30:56Z (GMT) No. of bitstreams: 2 Dissertação - Rafael Ribeiro de Carvalho Vaz - 2017.pdf: 6167014 bytes, checksum: 8f2f0dbd609aac9c2bd59458c6d0c243 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2017-09-15T15:30:56Z (GMT). No. of bitstreams: 2 Dissertação - Rafael Ribeiro de Carvalho Vaz - 2017.pdf: 6167014 bytes, checksum: 8f2f0dbd609aac9c2bd59458c6d0c243 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2017-08-24 / Fundação de Amparo à Pesquisa do Estado de Goiás - FAPEG / This work presents an end-to-end methodology for the deployment of Self-Healing systems in the distribution grid in order to decrease DEC in feeders. The methodology developed consists in choosing feeders to receive Self-Healing systems based on their performance in relation to the DEC. It is determined load transfer units considering power flow by the current injection method and protection coordination for these feeders. Once the units are determined, a feeders rank is set to order the system implementation considering hypothesis testing and binary linear programming. For this rank, besides the DEC, the number of customers as well as critical loads are considered. The communication system is based on optical mesh networks where single and multi objective algorithms are applied to optimize the network in order to design a reliable and low cost communication system. / Este trabalho apresenta uma metodologia fim-a-fim para implantação de sistemas Self-Healing na rede de distribuição com o objetivo de redução da duração equivalente de interrupção por unidade consumidora em alimentadores. A metodologia desenvolvida consiste na definição dos alimentadores candidatos a receberem os sistemas de Self-Healing baseado no desempenho dos mesmos em relação a duração equivalente de interrupção por unidade consumidora. Nesses alimentadores são determinadas as unidades de transferência de carga considerando fluxo de potência pelo método de injeção de correntes e a coordenação da proteção. Uma vez escolhidas essas unidades, é realizada uma classificação dos alimentadores com o intuito de estabelecer uma ordem de prioridade para implantação dos sistemas utilizando teste de hipótese e programação linear binária. Para essa determinação de prioridade, além do DEC, são considerados a quantidade de unidades consumidoras bem como o de cargas críticas. Neste trabalho, o sistema de comunicação considerado para interligar os religadores é baseado em redes de fibra óptica em malha. Com intuito de dimensionar o sistema de comunicação de forma que este seja confiável e de baixo custo, são aplicados algoritmos mono e multiobjetivo para otimização da rede, sendo proposta uma metodologia envolvendo cálculo de confiabilidade da rede.

Filtered two-fluid models used to perform large scale simulations of gas-solid fluidized flows of industrial risers require closures for filtered parameters such as filtered and residual stresses, and interphase interaction forces mainly effective drag. Closure models for those filtered parameters may be derived by averaging over results of highly resolved simulations with microscopic two-fluid modeling. This work is a contribution in that context. Recent models for filtered parameters have been written as functions of filter size, filtered solid volume fraction, and filtered slip velocity. A recent study showed that macro-scale variables like domain average solid volume fraction and gas Reynolds number also significantly affect the filtered parameters. In the current work, in addition to these filtered and macro-scale variables, the effects of two new variables over the filtered parameters are investigated: filtered solid kinetic energy and sub-grid gas turbulence. It is shown that the filtered solid kinetic energy should be accounted for in the concerning correlations, thereby improving accuracy. Regarding gas turbulence, literature shows it has no significant effects on the motion of high Stokes particles. Extending on literature, this work investigates the sub-grid gas turbulence effects on meso-scale structures formed of high Stokes particles. Results showed that sub-grid gas turbulence has no significant effects on the meso-scale structures and corresponding filtered parameters. The open source code MFIX was used for all simulations. Ranges of dilute concentration of solid and gas Reynolds number typical of riser flow regimes were considered. A modified two-fluid model with microscopic formulation was used. The sub-grid gas turbulence was generated by means of a forcing function procedure which was implemented in the physical space, over the gravitational term of the momentum conservation equation of the gas phase. First, numerical simulations of the gas phase alone were performed, accounting for literature available data, in order to set a turbulent gas field and calibrate the turbulence intensity. Then, the forcing function was introduced in to the two-fluid model and various gas-solid flows were simulated. While the current results show the necessity of accounting for additional variables in the filtered parameter correlation, they also make it clear the necessity of further developments that are required in the search for better accuracy. / Modelos filtrados de dois-fluidos usados em simulações de grandes escalas de escoamentos fluidizados de gás-sólido de risers industriais exigem fechamentos para parâmetros filtrados tais como as tensões filtradas e residuais, e forças interativas interfases, principalmente arrasto efetivo. Modelos de fechamento para estes parâmetros filtrados podem ser gerados a partir de procedimentos de media aplicados sobre resultados de simulações altamente resolvidas com modelagem microscópica de dois-fluidos. Este trabalho é uma contribuição neste contexto. Modelos de fechamento recentes para parâmetros filtrados tem sido formulados em função de tamanho de filtro, fração volumétrica de sólido filtrada, e velocidade de deslizamento filtrada. Estudo recente mostrou que variáveis de macro-escalas como fração volumétrica de sólido e número de Reynolds de gás médios no domínio também afetam significativamente os parâmetros filtrados. No presente trabalho, além dessas variáveis filtradas e de macro-escala, os efeitos de duas novas variáveis sobre os parâmetros filtrados são investigados: energia cinética filtrada do sólido e turbulência submalha do gás. Em relação à energia cinética filtrada do sólido, mostra-se que a sua consideração refina as correlações em questão, contribuindo assim para melhor acuracidade. Com relação à turbulência do gás, a literatura mostra que não tem efeitos significativos no movimento de particulados de elevados números de Stokes. Acrescentando à literatura, este trabalho investiga os efeitos da turbulência sub-malha do gás sobre estruturas de meso-escala formados de particulados de elevados números de Stokes. Os resultados mostraram que a turbulência sub-malha do gás não tem efeitos significativos sobre estruturas de meso-escalas e parâmetros filtrados correspondentes. O código aberto MFIX foi usado para todas as simulações. Faixas de concentração diluída de sólido e número de Reynolds típicos de escoamentos em risers foram considerados. Um modelo modificado de dois fluidos com formulação microscópica foi utilizado. A turbulência sub-malha do gás foi gerada por meio de um procedimento de \'forcing function\' que foi implementado no espaço físico, sobre o termo fonte gravitacional da equação de momentum da fase gás. Primeiramente, simulações numéricas da fase gás foram realizadas separadamente, levando-se em conta dados disponíveis na literatura, a fim de gerar um campo de gás turbulento e calibrar a intensidade de turbulência. Posteriormente, a \'forcing function\' foi introduzida no modelo de dois-fluidos e vários escoamentos de gás-sólido foram simulados. Enquanto os resultados obtidos mostram a necessidade de consideração de variáveis adicionais para correlação de parâmetros filtrados, também deixam claro a necessidade de desenvolvimentos mais aprofundados na busca de melhor acuracidade.

Dada a complexidade hidrodinâmica de escoamentos gás-sólidos em leitos fluidizados muitas pesquisas têm sido realizadas abordando questões relativas a natureza altamente instável desse processo de fluidização. O conhecimento dessa hidrodinâmica permite estabelecer os parâmetros corretos para reação e transporte de massa em um reator de leito fluidizado, permitindo tomada de decisões à respeito do desempenho do reator. Modelos Eulerianos do contínuo, onde as fases coexistentes são tratadas como meios contínuos interpenetrantes, representam a aproximação mais prática para simular a hidrodinâmica de escoamentos gás-sólido em reatores de leito fluidizado. Aplicando-se procedimentos de médias de Euler obtêm-se os modelos hidrodinâmicos A e B desenvolvidos por GIDASPOW [27] do IIT/ANL, representando as equações conservativas e constitutivas de modelos de duas fases separadas para escoamentos gás-sólido. Neste trabalho, aplica-se um modelo Euleriano de duas fases separadas, desenvolvendo simulações numéricas para diferentes dimensões de colunas ascendentes de leito fluidizado circulante. Para resolução das equações discretas no domínio bidimensional, aplicam-se malhas refinadas. Analisam-se os efeitos nos resultados de simulação causados por essas malhas. Apresenta-se uma discussão geral sobre o modelo aplicado em torno dos resultados obtidos. / Due to the hydrodynamic complexity of gas-solid flows in fluidized beds, several researches have been conducted addressing issues concerning the highly unstable nature of this fluidization process. The knownledge of its hydrodynamics allows establishing the correct parameters for the reaction and mass transport in a fluidized bed reactor, improving its performance. Continuous Eulerian models, where the phases are treated as a continuum, represent the most practical approach to simulate the hydrodynamics of gas-solid flows in fluidized bed reactors. By applying Euler average procedures, it is possible to obtain the hydrodynamic models A and B, developed by GIDASPOW [27] at IIT/ANL. These models represent the conservative and constitutive equations for the models of two separate phases for gas-solid flows. In this work, an Eulerian two-fluids model is applied, developing numerical simulations for different risers of circulating fluidized bed. For the resolution of the discrete equations in the two-dimensional domain, refined meshes are used. An analysis of their effects on the results of the simulation is performed. A general discussion about the results obtained by the model applied is also presented.