Adsorption of asphaltene molecules at the oil-water interface induces the formation of a complex microstructure, which stabilizes emulsions and impairs the efficiency of crude oil refining. In this work, we design a set of new shear rheology protocols to assess the effect of polar and non-polar solvents on indigenous Brazilian (BR) asphaltene adsorption. Moreover, the asphaltene morphology upon addition of solvents with distinct aromaticities is investigated by SEM microscopy. Our findings indicate that asphaltenes are a polycondensate aromatic island-type structure that forms reversible films when polar solvents are placed on top of the adsorbed film. The interfacial study also reveals that non-polar solvents induce precipitation/adsorption of asphaltenes at the oil-water interface. The asphaltene aggregates, upon the presence of weakly polar solvents, can consolidate into a more close-packed pattern, suggesting that network growth and asphaltene self-arrangement are directly related to the aromatic content. We explore the differences in asphaltene structuring and how it affects the extent of spontaneous emulsification. We find that the rate of emulsification is directly related to the chemical configuration of asphaltenes. Finally, we study the addition of stearic acid (SA) to asphaltene solutions in deionized water (DW) and synthetic water (SW) to better understand how surface and rheological properties are affected by competitive adsorption. We find that single SA are more prone to form liquid-like rather than solid-like films at the air-water interface. Furthermore, we show that the interfacial activity of our asphaltenes is enhanced in the presence of electrolytes and is dependent of the solvent aromaticity.

Link da defesa:

https://puc-rio.zoom.us/j/91264826112?pwd=UFdqQk0wQXBwcGNHUFNJbGZlVEdZZz09

Meeting ID: 912 6482 6112

Passcode: 175095

Reservoir characterization is an important tool for production/reservoir management. Well tests are commonly used in reservoir characterization and are the only source of dynamic data during the exploitation period. These tests typically measure the pressure and temperature responses at a well during controlled production, injection, or static conditions. Generally, only pressure data is post-processed in reservoir characterization. However, considering only pressure data can lead to misinterpretation associated with the neglected thermal effects, causing errors in reservoir properties estimation and consequently inefficient reservoir management. Besides that, pressure data have several noise sources that may compromise the accuracy of test results. Recent results have shown that temperature data can be used to improve reservoir parameter estimation. In this work, the ensemble smoother with multiple data assimilation method (ES-MDA) was applied in synthetic cases created by an in-house non-isothermal reservoir-well flow simulator that considers the Joule-Thomson heating and cooling, adiabatic fluid expansion/compression, conduction, and convection effects in the thermal energy balance equation. The synthetic measured data was obtained by adding gaussian and harmonics noises to the numerical predictions to simulate equipment and tidal effects, respectively. A sensitivity analysis of the effect of the CD matrix used for updating parameters of the ES-MDA method on the parameters estimations was carried out. The results show that adding temperature data to the observed data in the history matching improves the estimates of the reservoir parameters, especially for the skin region and reservoir porosity. For the analyses in which the pressure data had the addition of harmonic noise, the inclusion of temperature data also proved to be of great importance for an accurate characterization of the reservoir.

Link da defesa:

https://puc-rio.zoom.us/j/94290757371?pwd=UHh6WVB0dzUvL0pmM0RjT0xPQTU5dz09

Meeting ID: 942 9075 7371
Passcode: 000906
 

Suspensions correspond to a class of materials vastly used in a large set of applications and industries. Due to its extreme versatility, they have been the focus of numerous studies over the past decades. Suspensions are also very flexible and can display different rheological properties and macroscopic responses depending on the choice of parameters used as input in the system. More specifically, the rheological response of suspensions is intimately associated to the microstructural arrangement of the particles composing the medium and external factors, such as how strongly they are confined and particle rigidity. In the present study, the effect of particle rigidity, confinement, and flow rate on the microstructure of highly concentrated suspensions is studied using Core-Modified Dissipative Particle Dynamics. Preceding this main study, two other steps were necessary to guarantee a reliable and realistic simulation system, which consisted, essentially, on performing parametric studies to understand and estimate the appropriate values for wall-particle interaction parameters.
Chapter 2 and Chapter 3 of the present work address parametric studies performed to assist the input parameters choice to prevent particle penetration in a wall-bounded system. In Chapter 2 a simpler system, composed of solvent and walls, is built and the interaction parameters and wall densities are adjusted. In Chapter 3 the interactions are set for suspensions. In the latter case multiple parameters play a role in penetration and the traditional way to investigate these effects would be exhaustive and time consuming. Hence, we choose to use a Machine Learning approach to perform this study. Once the parameters were adjusted, the study of confinement could be carried out, in Chapter 4. The main goal of this Chapter was to understand how the microstructure of concentrated suspensions is affected by flow rate, particle rigidity and confinement. It was found that very soft particles always form a giant cluster regardless the confinement ratio; the difference being on how packed the particles are. In the rigid case, a stronger confinement leads the formation of larger clusters. Chapter 5 addresses a machine learning study carried out to predict the rheology of unconfined suspensions. The main contribution of this work is that it was possible to understand and adjust simulation parameters that affect unrealistic the physics of the system and develop a computational domain that enables to systematically study confinement effects on suspensions.

Link da defesa:

https://cwru.zoom.us/j/98796685273?pwd=ajBPa3N3dS9KOHdBOWJoeXlWRlVZQT09

Desde os primórdios da Revolução Industrial, caldeiras têm percorrido um longo caminho de desenvolvimento visando melhor segurança e maior eficiência. Em particular, o mercado de caldeiras flamotubulares tem exigido cada vez mais projetos personalizados dos fabricantes, abrindo caminho para o uso de modelos de caldeiras como eficientes ferramentas de projeto. Neste trabalho, foi desenvolvido um modelo termodinâmico para uma caldeira flamotubular, de três passes, operando em regime permanente com gás natural boliviano. Foi testada a sensibilidade no perfil de temperatura de várias correlações de troca de calor encontradas na literatura. Um planejamento composto central permitiu avaliar simultaneamente o impacto da variação de seis parâmetros (diâmetro dos tubos no primeiro, segundo e terceiro passes, comprimento da caldeira, número de tubos nos passes 2 e 3 e porcentagem de excesso de ar na combustão) e obter um modelo polinomial empírico a fim de otimizar as eficiências energética e exergética de caldeira. Um teste de malha foi realizado com temperaturas no final da chama, câmara de combustão, segundo e terceiro passes e apontou para uma divisão suficiente de 200 volumes de controle em cada passe. O modelo foi também validado com dados obtidos experimentalmente, observando-se que a diferença entre os valores experimentais e teóricos das temperaturas do primeiro e segundo passes é inferior a 22,9ºC e menor que 43,9ºC no terceiro passe. Com os modelos polinomiais robustos obtidos com o planejamento de experimentos, foi possível observar os parâmetros mais significativos, tanto para as eficiências energética e exergética, assim como para a destruição de exergia. São eles, por ordem decrescente de importância: o excesso de ar; as duas contribuições associadas ao comprimento da caldeira; o número de tubos nos segundo e terceiro passes e; o diâmetro dos tubos no segundo e terceiro passe e o diâmetro do tubo no primeiro passe. Além disso, foi possível quantificar o trade-off entre o aumento da superfície de troca de calor e a manutenção do nível de turbulência do escoamento.

Link da defesa:

https://puc-rio.zoom.us/j/93257958274?pwd=UmN3dnRFV2NFenJDTUZCSXlXR2dHUT09

Meeting ID: 932 5795 8274

Passcode: 976191

Flexible beams are becoming ubiquitous in several industrial applications, as new projects often aim for lighter and longer structures. This fact is directly related to the new challenging demands on structural performances, or it is a simple consequence of the engagement of industries in cost reduction programs (usage of less material). Flexible beams are usually modeled under the assumption of large displacements, finite rotations, but with small strains. Such hypotheses allow the equation of motion to be built using co-rotational finite elements. The co-rotational formulation decomposes the total motion of a flexible structure into two parts: a rigid body displacement and an elastic (small) deformation. This way, the geometric nonlinearity caused by the large displacements and rotations of the beam's cross sections can be efficiently computed. One of the novelties of this thesis is the direct usage of the equation of motion generated by a co-rotational finite element formulation in the computation of nonlinear normal modes (NNM). So far, most of the dynamic analyses with co-rotation finite element models were restricted to numerical integrations of the equation of motion. The knowledge of NNMs can be beneficial in the analysis of any nonlinear structure since it allows a thoroughly understanding of the vibratory response in the nonlinear regime. They can be used, for example, to predict a hardening/softening behavior, a localization of the responses, the interactions between modes, the existence of isolas, etc. The Rosenberg’s definition of NNM as periodic solutions (non-necessarily synchronous motion) is adopted here. The Harmonic Balance method and the Shooting methods are presented and used to compute periodic solutions of nonlinear systems. A numerical path continuation scheme is implemented to efficiently compute NNMs at different energy levels. Numerical examples show the capability of the proposed method when applied to co-rotational beam elements.

Link da defesa:
https://puc-rio.zoom.us/j/97235967592?pwd=Q205cUpEeU1xUDhpcTZ0ZVJ2Y09lUT09

Ultrasonic non-destructive evaluation (NDE) is of extreme importance inthe oil and gas industry (OGI), especially for assets and structures subjected to conditions that accelerate failure mechanisms. Despite being widely spread, ultrasonic non-destructive methods depend on a specialized workforce, thus being error-prone and time-consuming. In this context, pattern recognition methods, like Machine Learning (ML), fit conveniently to solve the challenges of the task. Hence, this work aims at applying ML techniques to address the interpretation of data acquired through ultrasonic NDE in the context of the OGI. For that purpose, this dissertation involves three case studies. Firstly, Ultrasonic Guided Wave (UGW) signals are used to classify defects present in welded thermoplastic composite joints. Results have shown that, when using features extracted with autoregressive models, the accuracy of the ML model improves by at least 72.5%. Secondly, ultrasonic image data is used to constructo an automatic weld diagnostic system. The proposed framework resulted in a lightweight model capable of performing classification with over 99% accuracy. Finally, simulation data was used to create a deep learning model for estimating the severity of corrosion-like defects in pipelines. R2 results superior to 0.99 were achieved.

Em diferentes momentos no seu ciclo de vida, o poço passa pela situação na qual ele está desconectado da superfície, mas é necessário monitorá-lo. As soluções atuais têm custo muito elevado e não estão alinhadas com as novas diretrizes ambientais da indústria. A Elétrons Livres está desenvolvendo um sistema de captação de energia submarina com o objetivo de fornecer energia para os sensores de baixo consumo que farão esse monitoramento. Este sistema consiste em uma miniturbina, um sistema de geração e um sistema de armazenamento.

 

The Nonlinear Hydraulic Diffusivity Equation (NHDE) models the single-phase flow of fluids in porous media considering the variation in the properties of the rock and the fluid present inside its pores. Often, the dimensionless linear solution for the flow of oil is computed using the Laplace and Fourier transform or Boltzmann transformation and provides the unsteady pressure field in Cartesian coordinates given by complementary error function erfc(xD, yD, tD) and in cylindrical coordinates described by the exponential integral function Ei(rD, tD).
This work develops a new analytical model based on an integro-differential solution to predict the formation mechanical damage caused by the permeability loss during the well-reservoir life-cycle for several oil flow problems. The appropriate Green’s function (GF) to solve NHDE for each well-reservoir setting approached in this thesis is used. The general solution is implemented in the Matlab® and the mathematical model calibration will be carried out by comparing the solution obtained in this work to the porous media finite difference oil flow simulator named IMEX®. The general solution of the NHDE is computed by the sum of the linear solution (constant permeability) and the first order term of the asymptotic series expansion, composed of the nonlinearities present in the hydraulic diffusivity function η(p). The second term of asymptotic series expansion is obtained by solving a Volterra’s second kind integro- differential equation in Matlab® and the permeability functions for some types of reservoir rock will be obtained through laboratory correlations, generated from field data. The instantaneous permeability loss effect is clearly noticed in the diagnostic and specialized plots.

Em reservatórios de gás retrógrado com pressões inferiores à pressão de orvalho, a produtividade dos poços pode ser comprometida devido ao aparecimento e acúmulo da fase líquida nas suas imediações. Este fenômeno é conhecido como bloqueio por condensação retrógrada e está associado à uma série de desafios para compreendê-lo. Um deles é a determinação da permeabilidade relativa das fases líquida e gasosa, que comumente advém de curvas de permeabilidade obtidas a partir da extrapolação de poucos dados experimentais. Dessa forma, elas tendem a não representar fielmente efeitos importantes para o escoamento, comprometendo a precisão da modelagem do fluxo no meio poroso. A fim de investigar o efeito das curvas de permeabilidade relativa sobre a formação de bancos de condensado, foi desenvolvido um modelo composicional em escala de reservatório para o estudo do escoamento de gás e condensado. Com o modelo, o uso de curvas de permeabilidade relativa obtidas através de simulação do escoamento de gás retrógrado na escala de poros e de correlações propostas na literatura foi avaliado. Considera-se: sistema isotérmico, escoamento bifásico e incorporação dos efeitos de forças capilares por meio do modelo de permeabilidade relativa. Equações de balanço molar e consistência de volume formam um sistema não linear resolvido pelo método de Newton que fornece pressão e número de mols de cada componente, em todos os volumes de controle do modelo, a cada passo de tempo. Para o equilíbrio de fases, a equação de Peng & Robinson foi implementada com uma rotina de flash a pressão e temperaturas constantes. O modelo foi validado contra a solução analítica para sistema monofásico e por fim, o simulador obteve a evolução temporal das curvas de pressão e saturação em função da distância do poço, a fim de comparar o efeito dos diferentes modelos de curvas de permeabilidade relativa na predição do bloqueio por condensado. Os resultados foram obtidos variando-se a permeabilidade absoluta do meio e a vazão de gás imposta no poço, e o impacto desses parâmetros no acúmulo de condensado foi avaliado.

Resumo: In different industries and in our daily lives, a wide range of soft materials present intriguing behaviors under applied stresses. Classical examples include: personal care products, foods, polymeric materials, muds, cement systems, biological tissues and blood. All these materials and products undergo different transport processes either in nature or during production and end use. To better understand, predict, design, and control these transport processes it is crucial to know (i) which methods are appropriate to measure the flow properties of these complex soft materials in different conditions and in a reliable and reproducible way, (ii) how these materials and products deform and possibly flow under applied stresses, and (iii) which mathematical models can be used to describe and predict the wide range of complex behaviors observed. In this talk, I discuss ongoing research on these topics with particular emphasis on industrial processes. Several examples are provided, including well cementing, waxy crude oil transportation, polymer extrusion and 3D printing.

Biografia: Flávio é engenheiro químico formado pela UFRJ e fez mestrado e doutorado no Departamento de Engenharia Mecânica da PUC-Rio. Em 2011, trabalhou no centro de pesquisas da Halliburton aqui no Rio. Lá, ele atuou em projetos de pesquisa e desenvolvimento como engenheiro de pesquisa sênior e líder da área de mecânica dos fluidos. Em 2016, ele retornou para a PUC-Rio para um pós-doutorado em reologia e, em 2018, foi para a Ghent University na Bélgica como pesquisador de pós-doutorado na área de impressão 3D. Atualmente ele é professor assistente do Departamento de Engenharia Química, Têxtil e de Materiais da Ghent University.

Resumo: Ascending thoracic aortic aneurysms (ATAAs) are bulging enlargements of the aorta in the section just after it exits the heart. They affect approximately 15,000 people in the US per year and can be fatal in 80% of the cases when ruptured. On the other hand, surgical repair of the aorta has its own risks and may not be necessary for many patients. Surgical intervention criteria are defined primarily in terms of the maximum diameter and growth rate, but it has been shown that about 60% of the patients present complications before the intervention criteria are reached, and about 40% of patients who meet the criteria for surgery do not, in fact, experience complications even when surgery is not performed. Clearly, a better way to predict the long-term outcome of an ATAA is needed.

Victor Barocas and his research group, including Ph.D. candidate Marisa Bazzi, use a combination of experimental and computational techniques to explore the mechanisms of ATAA formation, growth, and failure. For the talk, Victor will introduce the problem and discuss the various approaches his group is using, and he will provide context for Marisa's work in particular.

Marisa will then discuss her work exploring how macroscale biofluidic, geometrical and biomechanical factors affect the microscale growth and remodeling (G&R) process of the aortic tissue leading to the formation of ATAA. She will describe her combined experimental and computational approach using data from mice and humans to study the relevance of the main biomarkers on the ATAA outcome. She will then shift to her future work and discuss her planned creation of a fluid-structure-remodeling computational platform to study ATAA and other diseases that involve aberrant remodeling of the arterial tree.

Resumo: Em períodos de chuva, certa expectativa se cria no cenário brasileiro sobre os próximos eventos naturais que se consolidariam em uma catástrofe. Movimentos de massa são fenômenos naturais tipicamente brasileiros, ocasionados por desestabilização de taludes pelo efeito pluvial. Entretanto, quando tais eventos se constituem com volumes significativos, podem impactar consideravelmente infraestruturas e centros urbanos, podendo causar efeitos catastróficos. Monitorar e prever a gênese destes eventos e sobretudo quantificar sua evolução são tarefas de interesse científico, social e industrial. De um ponto de vista mecanicista, o fluxo de materiais granulares, lamosos e/ou misturas aquosas pode ser entendido como o escoamento de superfície livre de um fluido-tipo. Um fator agravante do comportamento destes escoamentos é a susceptibilidade ao aparecimento de instabilidades hidrodinâmica, as quais potencializam a imprevisibilidade do fenômeno. Prever a evolução destes escoamentos, em condições de permanência e uniformidade, é uma tarefa complexa. Quando ainda se considera tais instabilidades, torna-se ainda uma tarefa desafiadora. Nesta palestra, buscaremos entender, sobre um viés matemático, quais as complexidades da modelagem matemática do problema, como estas instabilidades surgem, e como podemos mensurar suas propriedades. Trata-se de uma tarefa em desenvolvimento pela comunidade científica, e que busca aplicar ainda tais modelos à determinação de áreas de risco e dimensionamento de estruturas de segurança e contenção como forma de aumentar a resiliência de centros urbanos e industriais à eventos naturais deste tipo.

Bio: O Prof. Guilherme Henrique Fiorot é professor do Departamento de Engenharia Mecânica da Universidade Federal do Rio Grande do Sul (UFRGS) e pesquisador na área de Mecânica dos Fluidos, atuando no Programa de Pós-Graduação em Engenharia Mecânica (UFRGS). Possui graduação em Engenharia Elétrica pela Universidade Estadual Paulista FEIS/UNESP (2009) e mestrado em Engenharia Mecânica pela mesma Instituição (2012). O pesquisador é doutor pelo Institut National des Sciences Appliquées de Rennes (França) na área de Hidráulica de Superfície Livre (Engenharia Civil), e doutor pela Universidade Estadual Paulista FEIS/UNESP em Mecânica dos Fluidos (Engenharia Mecânica). Desenvolve pesquisa nas vertentes teórica, numérica e experimental em escoamentos multifásicos de diferentes naturezas, principalmente aqueles relacionados a escoamentos geofísicos. Entre os temas centrais de pesquisa de interesse se encontram os escoamentos naturais de superfície livre para fluidos não newtonianos, o desenvolvimento de instabilidades hidrodinâmicas naqueles escoamentos, além de interações entre escoamento e estruturas, e o desenvolvimento de experimentos empregando sistemas de medição fotométricos.

Foam is widely used in oil recovery operations to improve sweep efficiency, in gas storage and acidization operations, and to solve problems
caused by either a thief zone or gravity override. Foam, which can be preformed and injected into the reservoir or produced in-situ through the pore space, fills the high permeability areas known as thief zones and diverts the displacing fluid into the direction of trapped oil, reducing the relative permeability of gas and leading to a more stable displacement front. The efficiency of these processes largely depends on the generation and stability of the foam films (lamellae) residing in the pores. The mobility of the injected gas is reduced when foam is formed; this reduction is attributed to the reduction of the gas phase relative permeability. The lamellae formed create resistance against the gas flow, impeding its free
motion inside the porous media. The lamellae population that composes the foam is direct related to surfactant concentration, and their flow and
mobility are functions of the pore geometry and foam properties. However, the dynamics of foam formation in porous media is not fully understood due to its complexity. The goal of the first part of this research is to understand the dynamic process of gas invading a two-dimensional porous media glass model occupied by a surfactant solution and forming foam. The second part focus on foam formation and its impact on oil displacement during SAG (surfactant-alternating-gas) injection. A microfluidic setup composed of a glass micromodel, syringe pump, pressure transducer and microscope, was used to visualize the pore-scale displacement and correlate the evolution of lamellae formation during the injection process with pressure difference for different flow conditions through image processing. The dynamics of lamellae formation is reported and related to macroscopic flow behavior.

Este trabalho consiste no estudo de um sistema de geração isolada de energia. Um sistema isolado é aquele que não está ligado à rede. Seus objetivos consistem na construção de um simulador e na realização de testes experimentais. O sistema possui como componentes: painéis fotovoltaicos, grupo de baterias, inversor de corrente e um sistema integrado (reformador + célula combustível). Onde, o reformador utiliza o etanol como combustível para geração de gás reformado, e a célula combustível é do tipo LTPEM (Low Temperature Proton Exchange Membrane).  O gás reformado é composto, principalmente, por hidrogênio e óxidos de carbono, além do metano em baixas porcentagens. O sistema foi desenvolvido para que um consumidor fosse atendido por um SIGFI. Sendo seus diferenciais: a utilização de gás reformado, não possui compressor e nem armazenamento. Fez-se um volume de controle em cada componente do sistema, analisando alguns dados como: potência, vazão, consumo elétrico, entre outros. Com estes dados, foi possível o desenvolvimento de um simulador em Excel e VBA (Visual Basic for Applications). Após isso, testes em conjunto foram feitos, para que estes dados experimentais fossem comparados aos simulados, a fim de se validar o simulador. Ou seja, comparar os dados simulados com dados reais. Com o simulador validado, foram feitas diversas análises, com alguns cenários, pré-estabelecidos, a fim de se comprovar ou não a viabilidade técnica do sistema de geração isolada de energia elétrica.

Durante a produção de petróleo e gás em águas profundas, problemas de garantia de escoamento podem ocorrer. Dentre estes problemas, a formação de hidratos de gás natural pode interromper a produção e causar problemas para a operação. Neste seminário vamos conversar sobre a ocorrência deste fenômeno e sobre os desafios enfrentados pelas operadoras de petróleo para gerenciar este problema.

Bio: Adriana Teixeira trabalha no Centro de Pesquisas da Petrobras, CENPES, há 19 anos. Como Química de Petróleo Sênior e especialista em hidratos em Garantia de Escoamento, coordena projetos de pesquisa que buscam atender as necessidades operacionais da empresa. Graduada em Química pela USP-SP, obteve seu título de Mestre no Instituto de Químca da UFRJ. É aluna de Doutorado na Escola de Química da UFRJ, e o tema de sua tese está relacionado ao estudo de suspensões de hidratos.   

A deposição de parafina é um dos principais problemas quando se trata de garantia de escoamento. O fluido que sai quente dos reservatórios é transportado por tubulações que trocam calor com o ambiente frio. Quando o fluido atinge temperaturas abaixo da TIAC (Temperatura Inicial de Aparecimento de Cristais), partículas de sólido se precipitam. Isto causa não só aumento da viscosidade, mas também a formação de depósito nas paredes da tubulação, diminuindo a seção livre de escoamento e levando a um aumento da potência necessária para o deslocamento do fluido, ou até mesmo obstruindo totalmente a tubulação, resultando em um grande esforço de manutenção e causando grandes perdas de capital. Este problema é bastante recorrente, sendo assim, prever a deposição da parafina é crucial para projetos e operações de tubulações, sendo preciso desenvolver modelos numéricos que forneçam resultados de forma acurada e eficiente. Considerando que as tubulações são muito longas, o presente trabalho propõe um modelo hidrodinâmico unidimensional, com um modelo bidimensional de transferência de calor obtido através de um processo de marcha ao longo do duto. Todas as propriedades e fração volumétrica de sólido são determinadas em função da composição do fluido, pressão e temperatura, a partir de interpolações de tabelas criadas com um modelo termodinâmico em uma etapa de pré-processamento. O modelo considera a existência de depósito na parede da tubulação quando a fração de volume de sólido é igual ou maior do que 2 %. A modelagem proposta foi avaliada em diferentes situações, em escala laboratorial e de campo, apresentando dados de temperatura, pressão e espessura do depósito com concordância razoável com dados da literatura, mostrando que a modelagem implementada reproduz satisfatoriamente o comportamento físico do fenômeno de deposição de parafina.

Este trabalho faz uma análise da variação dos parâmetros que têm importância na propagação de fraturas hidráulicas e da influência desses parâmetros na contenção do fraturamento. Os experimentos numéricos foram feitos em um modelo 2D utilizando um simulador de elementos finitos com acoplamento sequencial hidrodinâmico, tendo como premissa o comportamento estacionário dos processos envolvidos. Inicialmente, foram feitos testes de validação das soluções numéricas empregadas neste trabalho a partir de casos cujas soluções são bem conhecidas. Então, efeitos de variações de poropressão, de estado de tensões, propriedades das rochas, intervalos de início da fratura hidráulica, efeitos térmicos e o dano à formação permoporosa foram utilizados para avaliar a contenção da fratura hidráulica. Primeiramente, os efeitos foram avaliados separadamente e, em seguida, foram combinados aos pares, por meio de sorteio, e então avaliados. Os estudos levaram à conclusão de que o fator de maior influência para o início da propagação da fratura hidráulica na rocha capeadora (primeiros metros) é o valor da tensão mínima de confinamento do reservatório e a continuidade da propagação vertical na rocha selante é dominada pelo contraste de tensões entre rochas reservatório e capeadora. Entretanto, os demais parâmetros exercem influência na contenção do fraturamento hidráulico e devem ser levados em consideração neste tipo de estudo, principalmente os que servirão de insumo para a tomada de decisões.

O presente trabalho apresenta evidências e visa confirmar as tendências gerais de que reparos feitos de materiais compósitos podem ser utilizados para restaurar completamente a integridade estrutural, desconsiderando efeitos de tempo, de segmentos de dutos que apresentam anomalias (defeitos) de perda de metal em sua superfície externa, causados por corrosão ou erosão. Sete espécimes tubulares de API 5L Gr. B e dez espécimes tubulares de API 5L X65 foram testados até a ruptura por meio de testes hidrostáticos monotônicos, com diâmetro externo nominal de 323,9 mm e espessura nominal de 8,38 mm. Duas amostras testadas não tinham defeitos. Seis espécimes apresentavam defeitos e não foram reparados. Nove espécimes apresentavam defeitos e foram reparados usando três sistemas de reparo de compósito diferentes. Os defeitos continham profundidades variando entre 40%, 55% e 70% da espessura do tubo. Os sistemas de reparos foram baseados em compósitos de epóxi e fibra de carbono, epóxi e fibra de vidro, e poliuretano e fibra de vidro. (FREIRE, VIEIRA, et al., 2019)
Para a constatação dos resultados experimentais foram implementadas metodologias analíticas e numéricas. O método analítico utiliza equações normativas semiempíricas (ISO 24817 e ASME PCC-2) e da teoria de reparos para dutos (Equação de Equilíbrio). A avaliação numérica consiste no método de elementos finitos para modelagem do problema através do software ANSYS®.
Os sistemas de reparo dos três fornecedores mostraram ser eficazes e que reparos de compósitos podem restaurar completamente a integridade estrutural de dutos com defeito externo, desconsiderando efeitos de tempo, conforme o Fator de Resistência Remanescente (RSF) calculado. Os resultados dos testes mostraram que métodos analíticos e numéricos podem ser usados para compreensão do comportamento mecânico do problema.

Este trabalho analisa a propagação de uma doença epidemiológica com uma abordagem estocástica. Na análise, o número de indivíduos que cada membro infectado da população pode infectar é modelado como uma variável aleatória e o número de indivíduos infectados ao longo do tempo é modelado como um processo estocástico de ramificação. O foco do trabalho é caracterizar a influência do modelo probabilístico da variável aleatória que modela o contágio entre indivíduos na disseminação da doença e na probabilidade de extinção, e analisar a influência de uma vacinação em massa no controle da propagação da doença. A comparação é feita com base em histogramas e estatísticas amostrais do número de indivíduos infectados ao longo do tempo, como média e variância. Os modelos estatísticos referentes à parte que trata de uma população não vacinada são calculados usando simulações de Monte Carlo para 3 diferentes famílias de variáveis aleatórias: binomial, geométrica-1 e geométrica-0. Para cada família, 21 distribuições diferentes foram selecionadas e, para cada distribuição, 4000 simulações do processo de ramificação foram computadas. Os modelos estatísticos referentes a uma população parcialmente vacinada foram calculados usando simulações de Monte Carlo para a família de variável aleatória binomial. Para essa família, 21 distribuições diferentes foram selecionadas e, para cada uma delas foram escolhidas 6 diferentes percentagens de população vacinada. Para cada percentagem, foram analisadas vacinas com 4 diferentes eficácias. No total, foram realizadas 2.2 milhões de simulações, caracterizando o problema como big data.

A termografia por infravermelho (TIV) tem sido usada como ferramenta de avaliação não destrutiva para detectar falhas em componentes estruturais, desempenhando um papel importante nos programas de inspeção, de fabricação e manutenção. A TIV também pode ser aplicada para avaliar limites de fadiga para materiais estruturais a partir de testes de tensão cíclica uniaxial ou de tração uniaxial quase-estáticos. Recentemente, o método quase-estático para medir o limite de fadiga via TIV, foi estendido para estados de tensão biaxial atuando na superfície de corpos de prova de tamanho real de dutos com anomalias geométricas. Os resultados mostraram que o método biaxial quase-estático previu satisfatoriamente a localização de pontos críticos para a iniciação da fadiga em espécimes tubulares com mossas.

Neste trabalho, as distribuições de tensões e deformações nos pontos críticos foram determinadas usando correlação de imagens digitais (DIC), aplicada a sete espécimes “dog bone”, padrão uniaxial, e oito espécimes de dutos de três metros de comprimento com tampas de cabeça plana fabricados a partir de tubos de aço API 5L Grau B. A determinação do estado de deformação superficial via DIC e a determinação das variações de temperatura correspondentes para os pontos do material de superfície observados ao longo dos testes, permitiram determinar não apenas o limite de resistência ao escoamento do material, mas também o que conveniou-se chamar de temperatura mínima característica (TMC). Os testes uniaxiais usando espécimes “dog bone” para o mesmo material indicam que as TMCs medidas estão intimamente relacionadas com o limite de escoamento do material (MYS).

Esta pesquisa destaca as observações experimentais relacionadas às medidas de tensão-deformação-temperatura no início do escoamento em pontos observados em ensaios uniaxiais e biaxiais. O objetivo deste trabalho é usar o monitoramento da temperatura de estruturas sob carga crescente quase-estática para prever a que “distância”  do início do MYS estão os estados de tensão dos pontos observados.

This work studies the nonlinear dynamic behavior of oil well drillstring, which is a long slender flexible structure responsible for the drilling. Its elements and functions are presented, and numerical analyses are performed later. The work develops a computational code using the software MATLAB® for the numerical simulation of the column’s dynamic behavior using the finite element method. The corotational formulation is used for the implementation of geometric nonlinearity. The structure's discretization uses a beam element with six degrees of freedom per node and employs Euler-Bernoulli’s beam formulation. The Newton-Raphson method is responsible for solving the nonlinear system of equations. In addition, the solution procedure users the Newmark method for the time integration of the problem’s movement equations. A linear setup spring model is proposed to represent the contact between the borehole wall and the column. The proposed methodology and computational code capabilities are evaluated by comparing some results to analytical or numerical results of examples available in the literature. These results give reliability to analyze drillstring problems, which present the displacements and forces time series of the whole column and the buckling modes generated. The results show that the column is very sensitive to any boundary condition changing, which corroborates the complexity of the problem. Hence, the work proposes a reasonable basis for further developments, allowing the entire coupled drillstring analysis.

A crescente demanda energética verificada ao redor do mundo e a conscientização pública acerca dos efeitos deletérios do excesso de gases estufa na atmosfera vem colaborando para a articulação de compromissos de grande alcance em nome da adaptação das matrizes energéticas a formas ambiental e economicamente sustentáveis. A adesão à energias renováveis (como solar e eólica) e a descentralização da matriz energética por meio de tecnologias de geração distribuída (visando a melhoria da eficiência do uso da energia) são alguns dos movimentos mais relevantes realizados para fazer frente a essas demandas. Neste interim, o presente trabalho é dedicado à simulação numérica mediante o conceito 4E (Energy, Exergy, Environmental and Economic) de um sistema híbrido CHP (Combined Heat and Power) on-grid para atendimento de pequenas demandas residenciais ou industriais, tendo gás natural e energia solar como vetores energéticos preferenciais. O sistema inclui um reformador de gás natural para produção de gás de síntese rico em hidrogênio, uma célula a combustível com membrana de troca de prótons (PEM), painéis fotovoltaicos, baterias conectadas à rede elétrica por um inversor bidirecional, trocadores de calor e componentes auxiliares como compressores e boilers. Os componentes do sistema foram modelados separadamente com base em equações de conservação e seus modelos devidamente validados. Uma análise energética e exergética do reformador de gás natural foi conduzida mediante a metodologia de planejamento de experimentos a fim de avaliar a necessidade de considerar uma formulação complexa do combustível em vez de um substituto (metano puro). Posteriormente, estes modelos foram inseridos como módulos de uma rotina mais ampla destinada a simular o desempenho econômico do sistema integrado num intervalo de tempo de até 20 anos. Tal rotina, implementada no MATLAB, permite a flexibilização de critérios operacionais importantes como número de consumidores, configuração do sistema híbrido (armazenamento e participação de painéis fotovoltaicos), diferentes tipos de tarifa (convencional ou branca) e o possível uso de rejeito térmico para cogeração, enriquecendo o escopo de resultados obtidos. Paybacks entre 7 e 20 anos de operação do sistema foram alcançados para diferentes combinações dos parâmetros examinados considerando-se a adesão no ano de 2020, onde consumidores residenciais obtiveram resultados predominantemente melhores do que os industriais em virtude da demanda menos exigente dos primeiros. Foram também previstas reduções de até 50% no custo cumulativo total para consumidores residenciais referente a adesão ao sistema proposto por 20 anos, levando-se em conta a queda prevista nos custos de aquisição dos componentes para as próximas décadas. A avaliação do sistema em termos ambientais foi feita através da quantidade equivalente de CO2 por unidade de energia. Concluiu-se que a configuração completa, mesmo auxiliada por cogeração, supera a média de emissões da matriz energética brasileira (devido à alta participação das fontes renováveis nessa matriz), permanecendo, ainda assim, como uma opção melhor do que a combustão pura do gás natural, especialmente no que diz respeito ao atendimento de demanda térmica.

Fuel consumption in companies in the rail transport sector represents one of the largest operating expenses and one of the biggest concerns in terms of pollutant emissions. The high fuel consumption also entails a high representation in the emissions scope matrix (more than 90% of railroad emissions come from fossil fuel consumption). Aiming to seek constant operational improvement, numerous studies have been carried out proposing new tools to reduce fuel consumption in the operation of a freight train. In this way, it is important to highlight the improvement of train driving parameters that can be calibrated to reduce fuel consumption. To accomplish this goal, the present work implements two machine learning models to predict the energy efficiency of a freight train: random forest and artificial neural networks. The random forest achieves the best performance against the models, with an accuracy of 91%. To calculate how much each parameter influences the prediction model, this work also uses the technique of accumulated local effects for each parameter related to energy efficiency. The final results show that, within the four analyzed calibration parameters, the traction per transported ton indicator presented greater representation in terms of absolute impact on the energy efficiency of a freight train.

I will discuss methods to compute and model dry and saturated granular flows of non-Brownian, non-cohesive, and non-spherical rigid bodies at the particle scale. In the case of rigid bodies suspended in a fluid, saturated means that the surrounding fluid is a single-phase fluid. I will leave most of the hard core scientific computing details aside and elaborate on the potential of these methods as discovery-enabling tools complementary to experiments and theory. While promising, these methods feature limitations and offer challenges, both conceptual and computational. One of the conceptual challenges is related to how the particle-scale models can be upscaled (i.e. coarse grained) to tractable macro-scale models or constitutive equations, while computational challenges are primarily associated to resolving finer length and time scales in a faster way. I will attempt to partly address these limitations and challenges and illustrate what can be achieved with such computational models on 3 examples: (i) buoyancy driven rising and settling on single non-spherical particles in a Newtonian fluid, (ii) coarse grained model of drag, lift and torque coefficients in inertial suspensions of spheres and (iii) dry granular slumping of non-convex and entangled rigid bodies.

 Biografia:

Anthony Wachs is currently a full professor at the University of British Columbia (UBC). He is affiliated with the department of Mathematics and the department of Chemical and Biological Engineering, and is also an associate member of the department of Computer Science. He received his PhD in Mechanical Engineering in 2000 from the University of Grenoble, France, and right after worked for 15 years for IFP Energies nouvelles, a national lab in France specialized in the field of energy. He joined UBC in 2015 and since that date has been an active member of the well established UBC Fluids group (Frigaard, Balmforth, Feng, Elfring, Wachs) specialized in multiphase and non-Newtonian fluid flows. His and his group's research interests are particle-laden flows, viscoplastic fluids, computational methods and high performance computing. 

In the last few years, the automotive industry has reinvented itself to meet the demands of the international market, which has been increasingly competitive in a context with environmental laws each year more severe. One alternative to lower harmful greenhouse gases emissions over the life of the vehicle is electric cars. However, the production and disposal of electric batteries is still a major problem to be solved. Therefore, companies are also searching for other potentialities to increase the internal combustion engine's efficiency and develop green technology, such as Homogeneous Charge Compression Ignition (HCCI) or Spark-Assisted Compression Ignition (SACI). A MATLAB routine was created to predict the performance of SACI multi-mode combustion of natural gas using a two-zone thermodynamic model. This work performs sensitivity analysis for five performance parameters: thermal efficiency (η_th), indicated mean effective pressure (IMEP), NOx emissions, mean in-cylinder temperature (T ̅_avg), and auto-ignition timing (AIT), with several variables such as engine speed, fuel-air equivalence ratio, spark timing, compression ratio, and intake pressure, using the design of experiments tools to assess the factors’ impact. The Central Composite Design indicates that the engine speed and the fuel-air equivalence ratio were the most important SACI factors since they influence all engine performance parameters. On the other hand, the intake pressure was significant in three performance parameters (η_th, IMEP and T ̅_avg), as was the spark timing (NOx, T ̅_avg and AIT). The compression ratio, however, was relevant in only one of them (AIT). Furthermore, a Univariate Analysis was done to compare Spark-Ignition (SI) and SACI combustion techniques. The results show that at the same regime as traditional SI engines, SACI engines tend to be around 9% more efficient, NOx emissions drop notably, more than 99%, IMEP presents an increase of 45.4%, and T ̅_avg decreases 24.8%.

Ondas acústicas, sônicas ou ultrassônicas, podem ser empregadas para a telemetria sem fio como alternativa a sistemas eletromagnéticos, transferindo dados e energia ao longo de um canal formado por uma ou mais camadas de sólidos elásticos ou fluidos acústicos. Um exemplo é a interrogação de sensores passivos através de uma parede metálica. Nesta configuração, pelo menos um transdutor acústico é fixado em um lado da parede (face externa), onde uma fonte de alimentação elétrica é disponível. No lado oposto (face interna), onde os sensores estão instalados, são fixados um ou mais transdutores. Na maioria das aplicações estes transdutores são cerâmicas piezelétricas que geram e recebem sinais ultrassônicos. Ondas acústicas ultrassônicas se propagam ao longo do sólido elástico, transferindo energia e dados entre as duas faces, possibilitando a alimentação e interrogação dos sensores. Este tipo de configuração pode ser empregado em aplicações onde o uso de penetradores elétricos ou ópticos não é recomendado. Entretanto, a resposta das piezocerâmicas pode sofrer influências de variações de temperatura e da própria deformação mecânica da parede metálica na qual são fixados.
O presente trabalho procurou quantificar a influência da deformação mecânica e da variação de temperatura na comunicação entre dois transdutores piezocerâmicos ultrassônicos, aderidos à uma placa metálica por meio de adesivo epóxi. No estudo, tomou-se como parâmetro quantitativo o sinal S21, que é o logaritmo da razão entre a potência recebida pela saída do sistema (face interna da parede) pela potência transmitida pela entrada (face externa da parede). O trabalho apresenta comparações entre resultados experimentais e simulados através de um modelo numérico de elementos finitos desenvolvido no COMSOL Multiphysics®.
Os ensaios experimentais foram realizados com pastilhas piezocerâmicas circulares, do tipo PZT4, com diâmetro e espessura de 25 e 2 mm, respectivamente. Os transdutores foram fixados, de forma concentricamente alinhada e por meio de um adesivo epóxi, nas duas superfícies de uma placa de aço inoxidável AISI 316 L com 6 mm de espessura. O trabalho apresenta tabelas e funções para a amplitude do sinal S21 na frequência onde a transferência de potência é maximizada. Para os casos estudados, observou-se que a frequência ideal muda muito pouco com a temperatura ou a deformação da placa sobre a qual os transdutores são fixados, permanecendo com valores entre 0,988 e 0,995 MHz em todas as condições avaliadas. Em função da deformação da placa metálica, a amplitude do sinal S21 também variou muito pouco, de -3,70 para -3,14 dB, desde a condição indeformada da placa até a máxima deformação aplicada, que foi de 1250 µm/m. Quanto à variação com a temperatura, na faixa de 30 a 100 ˚C, mais uma vez observou-se apenas um pequeno aumento de 0,8 dB na amplitude do sinal S21. Entretanto, para temperaturas acima de 100 ˚C, o sinal passa a cair rapidamente. Em nenhuma das condições estudadas neste trabalho foi observado prejuízo na transferência de potência entre os transdutores, indicando que este tipo de comunicação pode ser uma alternativa robusta ao uso de penetradores elétricos.

When wells reach the end of their life cycle, there are permanently plugged and abandoned.  Plug and abandonment can easily contribute to 25% of the total cost of drilling exploration wells offshore Norway. The cost of running a plug and abandonment operation on some offshore production wells may have a cost impact similar to the cost of the original drilling operation. Therefore, cost efficient plug and abandonment technology is a necessity without compromising the scope of the operation.

Biografia:

Mahmoud Khalifeh is Associate Professor at University of Stavanger (UiS).  He holds a PhD in plug and abandonment of wells with a focus on materials optimized for permanent P&A.  Dr. Khalifeh has developed the permanent P&A (Plug and Abandonment) subject (MSc level) at UiS.  He has published several articles in scientific journals and conferences.  Dr. Khalifeh represents UiS as observer at Norwegian Plug and Abandonment Forum (PAF) since 2012. 

Thin liquid films play a big role in many real-life applications and are of indisputable interest to scientific and industrial researchers. Evidence of thin films are observed in nature in large scales such as snow avalanches in the mountains, lava flows on volcanoes and landslides, and in small scales such as the pulmonary airways and the eye surface. They are also widespread in many industrial applications, ranging from high-resistance thin film resistors, atomization, soft-lithography methods and several coating techniques such as dip, roll, slot, spin and curtain coating. Understanding the physical mechanisms contributing to the stability of thin liquid films is a challenging problem, as thin films' flows present a fluid-fluid interface which is free to deform. The interface is bounded between two liquids or a liquid and a gas, typically having its own dynamic properties from which interfacial tension effects and complex interfacial rheological behavior arises. Instability is usually driven by long-range intermolecular forces, also known as van der Waals attractions, and may result in the rupture of the layer. Numerical investigation is often used to understand the breakup dynamics of thin liquid sheets by addressing the evolution of the film thickness using either asymptotic derivations of the lubrication theory or interface tracking techniques. In this work, a computational investigation of the breakup dynamics of a stationary thin liquid sheet bounded by a passive gas with a viscous interface is presented. The Arbitrary Lagrangian-Eulerian method (ALE) is used to track the interface position. The rheological behavior of the viscous interface is modeled by the Boussinesq-Scriven constitutive law, and the numerical solution is obtained through finite element approximation. The results show that thin liquid film stability is influenced both by surface rheology and disjoining effects and that the viscous character of the interface delays the sheet breakup, leading to more stable films.

Resumo: TideWise develops and operates its own unmanned marine vessels (USVs). Our vision is one of a marine industry where most tasks can be done either autonomously or from the shore, removing people from harm's way, and reducing environmental impact while doing so. This presentation will talk briefly about the company, and we will then dive a bit in the technology - how we believe that system integration of hardware and software components is the key to turning the algorithms into a functioning remotely operated - and in the future fully autonomous - system.

Na última década a participação dos campos do Pré-Sal brasileiro na produção nacional de petróleo aumentou substancialmente, tornando esses campos responsáveis por mais da metade da produção nacional e com perspectiva decrescimento para os próximos anos. Os reservatórios de petróleo encontrados nessa região são caracterizados por espessuras que podem variar de poucas dezenas a centenas de metros, rochas com boas qualidades permo-porosas e presença de óleo leve, com elevado teor de gás associado e com contaminantes como CO2. Por estas características, diversos sistemas de produção instalados nestes reservatórios foram preparados e deverão adotar o método de recuperação suplementar com injeção alternada de água e gás. No presente trabalho foi feita uma análise paramétrica sobre a influência que propriedades de reservatório como espessura porosa, permeabilidade horizontal, permeabilidade vertical e razão kv/kh, e variáveis operacionais como vazão de operação, razão WAG e tempo de ciclo podem gerar no fator de recuperação em um cenário de produção típico do Pré-Sal brasileiro com e sem o efeito de segregação gravitacional. Pela comparação dos resultados em diversos casos de simulação numérica, foi possível identificar as variáveis com maior impacto. Em seguida, a partir de uma análise de sensibilidade foi gerada uma equação para estimar o fator de recuperação em função das variáveis selecionadas. Números adimensionais propostos na literatura para avaliar a segregação gravitacional de fluido sem meio poroso foram calculados e utilizados para gerar outras equações para estimar o fator de recuperação. As estimativas do fator de recuperação feitas por cada função foram comparadas com os valores simulados para cada caso e foram identificadas as funções que apresentaram as estimativas mais próximas. Tais funções poderão ser utilizadas para estimar o fator de recuperação no cenário proposto e com aplicação em análises preliminares para projetos de desenvolvimento de campos de petróleo.

Resumo:

Embora comumente tratados como homogêneos, a maioria dos materiais policristalinos apresenta comportamento altamente heterogêneo na microescala, que pode ser controlado por fatores internos (como parâmetros microestruturais) ou fatores externos (como os parâmetros de carregamento). O estudo da relação entre o comportamento não-homogêneo presente na microescala e a resposta do material na macroescala requer o uso de técnicas de homogeneização. Tais técnicas são, na maioria das vezes, baseadas no conceito do Elemento Volumétrico Representativo (RVE – Representative Volume Element) – o menor volume de material necessário para representar seu comportamento médio global. Neste seminário serão discutidos métodos de medição do RVE, através de campos de deformação de altíssima resolução obtidos através da técnica de Correlação Digital de Imagens (DIC – Digital Image Correlation), aliados a dados microestruturais obtidos de escaneamentos de Difração de Elétrons Retroespalhados (EBSD – Electron Backscatter Diffraction). A medição dos tamanhos dos RVE's em materiais sujeitos a carregamentos de plasticidade e fluência também serão apresentados, demonstrando que o tamanho do elemento volumétrico representativo é de fato uma função do comportamento microestrutural em si e, portanto, deve ser estudado caso a caso e não encarado como uma propriedade do material. Além disso, uma investigação sobre as relações entre mecanismos microestruturais de deformação e o grau de heterogeneidade do comportamento na microescala de materiais também será apresentada, com o tamanho do RVE sendo proposto como um indicador de transições entre mecanismos dominantes.

Biografia:

Renato Vieira se formou em Engenharia Mecânica pela PUC-Rio em 2013 e obteve seu grau de Mestre em Engenharia Mecânica também pela PUC-Rio em 2016, sob orientação do prof. José Luiz Freire. Em 2021 recebeu seu grau de PhD em Engenharia Aeroespacial pela Universidade de Illinois, sob orientação do prof. John Lambros. Seus interesses de pesquisa incluem o comportamento microestrutural de materiais submetidos a carregamentos de fluência e fadiga em altas temperaturas, e as relações entre os parâmetros microestruturais e a resposta macroestrutural dos materiais.  

This work proposes to investigate normal stress measurements for yield stress materials. Using an ARES-G2 TA Instruments’ strain-controlled rheometer, three different elastoviscoplastic materials were studied. A simple but effec- tive experimental procedure was developed to obtain repetitive results for both steady-state and transient measurements. These tests were put to comparison, and interesting results were obtained. For the Carpobol dispersion and the com- mercial hair gel, transient and steady-state results were quite similar. Although, the same was not observed for the highly concentrated emulsion. Also, vari- ables such as N1, ψ1, τ, γ ̇, G′, G′′, η, τy, σs and others were explored to better understand the stress state of these materials. Furthermore, normal stresses were significantly higher than shear stresses, while all materials were submitted to low shear rates. For high shear rates, different from the other two fluids, the highly concentrated emulsion displayed normal stress values that were lower than shear stress ones. Finally, in contrast with most results found in literature, only nega- tive normal stress results were obtained for the materials studied.

Fuzzy logic systems are machine learning techniques that can model mathematically uncertainties. They are divided into type-1 fuzzy, and type-2 fuzzy logic systems. The type-1 fuzzy logic system has been widely applied to solve several problems related to machine learning, such as control, classification, clustering, prediction, among others. However, as it presents a better mathematical modeling of uncertainties, the type-2 fuzzy logic system has received much attention over the years. This improvement in modeling is also accompanied by an increase in mathematical and computational complexity. Aiming to reduce these complexities to solve classification problems, this work presents the development and comparison of two Gaussian membership functions for a type-2 interval fuzzy logic system using the upper and lower method. Gaussian membership functions with uncertainty in the mean and with uncertainty in the standard deviation are used. Both fuzzy models covered in this work are trained by algorithms based on first order information. Furthermore, this work proposes the extension of interval type-2 fuzzy models to present multiple outputs, significantly reducing the computational cost in solving multiclass classification problems. Finally, aiming to contextualize the use of these models in mechanical engineering applications, this work presents the solution of a problem of fault detection in aircraft gas turbines.

Vehicle acceleration performance is an important feature to be evaluated when internal combustion engines and new fuels are being developed. Predicting this parameter is also of great significance, once automotive track tests requires long periods of time to be done and high costs with equipment, rental of the track, hiring people and displacement of vehicles and fuels. In addition, the results are directly affected by track surface irregularities and variations in weather conditions such as ambient pressure, temperature, air humidity, wind speed and crosswinds. Therefore, this work aims to use the data collected in bench tests with an internal combustion engine, in order to computational modelling speed recovery tests of its respective car. Doing these experiments in laboratory has the advantage of a greater control of the room ambient conditions and the engine operating parameters. Also, it reduces the expenses related to on-track testing. The proposed methodology simulates the traction force on the wheels based on the measured torque in dynamometer or from the pressure curves inside the combustion chamber with the aid of friction models for spark ignition engines. In order to validate the proposed model, it became necessary to perform speed recovery tests with the car on a chassis dynamometer. Also, seven different mixtures of ethanol and gasoline were used, and it was concluded that pure anhydrous ethanol promoted a higher acceleration capacity in most of the experiments but it had higher fuel consumption. It was also seen that hydrated fuels reduced performance. The simulations demonstrated a high precision in relation to the experiment, with a recovery time difference average of 0.51 seconds and standard deviation of 0.078. Also, the acceleration performances had errors smaller than 5.25%.

Em um reservatório de petróleo, a eficiência do deslocamento do óleo é bastante reduzida pela presença de fraturas e camadas de alta permeabilidade. Quando uma fase aquosa é injetada, esta flui preferencialmente por caminhos de baixa resistência, deixando grandes volumes do reservatório não afetados pelo processo de injeção. Uma alternativa para minimizar esse problema é bloquear as fraturas com um sistema líquido que forma um gel após algum tempo, forçando a fase aquosa a fluir pela matriz porosa. Os géis de silicato de sódio são uma das formulações que podem ser utilizadas para esse fim. No entanto, a cinética do processo de formação do gel de silicato é difícil de controlar, pois depende de muitas variáveis. Essa incerteza do processo de gelificação pode levar à formação de uma fase de gel fora da posição desejada, levando a diversos problemas, como perda de injetividade. Uma formulação líquida utilizando microcápsulas é uma proposta com um método inovador que tem como objetivo controlar o tempo de liberação do agente ativador e, consequentemente, melhorar o controle sobre o início do processo de gelificação. O primeiro passo no desenvolvimento desta técnica foi estudar a evolução da reologia do gel formado a partir do silicato de sódio com ácido clorídrico para avaliar a taxa de formação do gel em função de diferentes parâmetros do processo, como a concentração de Na-Si e HCl. Os resultados mostram que o processo de gelificação é função do pH da solução e que o tempo de gelificação é menor quanto maior a concentração de ácido clorídrico. No método proposto, a solução de ácido clorídrico é encapsulada e a gelificação só se inicia após o rompimento do invólucro da cápsula, que é desencadeado pela imposição de um gradiente de pressão osmótica.

Thin liquid sheets are present in a variety of systems and applications. Here, we are interested in double-layered sheets, which are common in the curtain coating process. In curtain coating, the liquid falls from a die forming a thin curtain before wetting the moving substrate. One of the most important process limits is the curtain breakup, which sets a lower limit for the coating liquid flow rate. Consequently, this flow rate lower limit defines the feasible minimum deposited film thickness. Experimental evidence have shown that using a two-layer curtain, with a viscoelastic thin layer, may delay the curtain breakup to lower flow ratios. The breakup of two-layer liquid sheets, composed of a Newtonian and a viscoelastic liquid, is studied by solving the differential equations that describe the evolution of the liquid sheet configuration until breakup. The effect of different parameters on the breakup time is determined. The results show the same behavior observed experimentally, thin viscoelastic liquid layer delays the breakup, stabilizing the liquid sheet.

As atividades de exploração e produção de petróleo em estruturas Onshore e offshore estão associadas ao risco iminente de vazamento de gás ou óleo, que encontrando qualquer tipo de ignição, pode causar incêndios e explosões com o potencial de causar perdas financeiras substanciais e mortes. No presente trabalho, a metodologia em CFD foi utilizada para analisar o vazamento de gás natural em um ponto localizado no módulo de compressão de gás em uma plataforma offshore. A turbulência do escoamento foi modelada com a metodologia RANS, empregando o modelo de duas equações κ-ε. Considerando a situação crítica, o vazamento foi analisado em regime permanente, e modelado como uma fonte pontual. Foram consideradas duas intensidades de vento: 0,5 m/s representando a condição de calmaria e 6,5 m/s representando a velocidade máxima registrada. Para cada condição de vento, foram estudados quatro pontos de falha do sistema (pontos de liberação do gás), ou seja, vazamentos no lado norte, sul, leste e oeste do compressor. A partir dos resultados obtidos, pontos críticos foram identificados. O volume da nuvem foi determinado para cada cenário e, posteriormente, um cálculo de TNT equivalente e de sobrepressão (comparação do potencial da mistura explosiva com a massa de TNT necessária para produzir efeitos semelhantes) foram realizados para identificação dos potenciais riscos a estrutura offshore em caso de ignição do material gasoso disperso. Os resultados obtidos coincidiram com medidas experimentais de sobre pressão em situações semelhantes. Identificou-se como o pior cenário, vazamento no compressor na direção sul, correspondente a ventos fortes, com alta concentração de gás natural numa grande região da plataforma.

Com o crescimento da extração de petróleo offshore nos últimos anos, em particular oriundo do pré-sal, faz-se necessário adequar as plataformas a estas novas necessidades, particularmente em relação ao flare. Para que este processo seja feito com segurança, demandam-se cada vez mais modelos precisos de determinação da troca de calor radiativa e tecnologias de bloqueio de radiação, tais quais os heat shields, para que os equipamentos da plataforma operem em condições de segurança adequada. O presente trabalho introduz duas técnicas de medição da emissividade da superfície de três materias: cobertura de Thermal Sprayed Aluminum (TSA), cobertura de otatemp1000 e uma superfície crua de aço naval. Foi efetuado um levantamento das etodologias de medição de emissividade e dos conceitos básicos dos processos de radiação e convecção natural, a partir disso foi determinada uma modelagem do problema a ser observado para as estimativas de calor e emissividade. O primeiro experimento consistiu na medição da emissividade com o auxílio de um termógrafo, nesse foi realizada a comparação de medição sobre uma placa plana e um tubo revestido para o TSA e o Jotatemp1000. No segundo experimento, realizou-se a medição da emissividade através da determinação das transferências de calor por convecção natural e radiação em um espaço anular em cilindros longos. As medições neste experimento foram realizadas para três diferentes pressões, para que se quantifique o efeito da convecção natural. No primeiro experimento, houve uma discrepância significativa entre os valores medidos nos tubos e nas placas, para o TSA obteve-se um valor médio de ε = 0, 414 ± 0, 030 para a medição na placa e ε = 0, 325 ± 0, 030 para a medição no tubo, no caso do Jotatemp ε = 0, 952 ± 0, 030 para a placa e ε = 0, 874 ± 0, 030 para o tubo. Notou-se também que em ambos os experimentos a emissividade de todos os materiais permaneceu praticamente constante para os valores acima de 100 oC, isto se deve ao fato de que abaixo desta temperatura a convecção natural compõe parte considerável do calor total transferido. Para o Jotatemp por exemplo, obteve-se os seguintes valores para a pressão de 1 atm (maior influência da convecção natural), ε = 0, 405 ± 0, 045 (T = 100oC), ε = 0, 445 ± 0, 046 (T = 120oC) e ε = 0, 465 ± 0, 046 (T = 142oC). Apesar disso, para o Jotatemp1000 e o aço naval, observou-se uma disparidade considerável nos valores de emissividade encontradas no primeiro experimento em comparação com o segundo. Constata-se que a medição de emissividade com termógrafo possui uma precisão maior, por ser uma medida direta, e dentre as metodologias utilizadas para a medição no espaço anular o mais confiável para pequenos intervalos de temperatura e a parametrização pela superfície dada pelo matlab. Ainda que os resultados apresentem tais discrepâncias, pode-se observar que o TSA obteve os menores valores de emissividade em todos os experimentos e, portanto, é o material mais indicado para aplicação em heat shields a fim de mitigar os efeitos de radiação na plataforma.

16:00 Nicole Lopes M de B Junqueira; Modelos de ordem reduzida aplicados à chamas laminares não pré-misturada, 16:15 Rebeca Mendes Carrijo; Otimização Estrutural da Seção Mestra de um Navio Petroleiro. 16:30 Rodolfo Luiz Gaudereto de Freitas; Redução da demanda energética nacional e da emissão de gases ligados ao efeito estufa: a importância da eficiência energètica.16;45 Rodrigo Borges Rabelo; Determinação dos esforços axiais em tubulações devido aos transientes hidráulicos a partir da análise numérica dos gradientes de pressão

16:00 Anderson Raniere Sobral da Silva; Tensões na parede de um poço horizontal induzidas por flutuação de pressão no fluído de perfuração.16:15 Fabio Pinheiro Cardoso; Projeto, simulação dinâmica e controle de um multirotor omnidirecional. 16:30 Lara Cardido Alvim: Controle de força e posição aplicado a sistemas de contato. 16:45 Lara Cristina Pinheiro de Araujo; Previsão do acúmulo de deformações plásticas em grãos de metais policristalinos através de Aprendizado de Máquina.

16:00 Gustavo de Oliveira e Souza; Análise do potencial mercado consumidor interno brasileiro no contexto da inserção do hidrogênio na matriz energética e seu potencial no comércio internacional. 16:15 João Pedro Callado Biscaia;  Análise 4E de um sistema híbrido solar/ SOFC com armazenamento em baterias e chiller de absorção para geração de potência, calor e frio. 16:30 Mariana Mano da Silva Pessoa; Modelagem numérica de uma unidade de gaseificação de resíduos sólidos análise energética, exergetica, econômica e ambiental. 16:45 Lara Schimith Berghe; Efeito de solicitações mecânicas na estabilidade de emulsões de água em óleo (A/O) em uma célula Four-Roll MIll.

This work investigates the impact of functionalized graphene suspensions on the rheology of a Carbopol aqueous dispersion. The graphene derivatives used were graphene oxide (GO) and amino-functionalized graphene oxide (AFGO). The variation of the functionalized graphene concentration and oxidation was evaluated and related with the suspensions rheological behavior. GO nanosheets were produced from synthesis of graphite oxide (GrO) by modified Hummers method, and they were characterized by XRD, Raman,TGA, FTIR, XPS, TEM and AFM techniques. The GO was functionalized with triethylenetetramine by microwave assisted reaction to produce the AFGO, which was characterized by TGA and XPS techiniques. The rheological behavior of these suspensions was characterized by oscillatory and steady-state flow experiments and zeta potential characterized the suspensions colloidal aspects. The GO structure characterization shows that oxygenated functional groups were incorporated in its graphitical surface. GO oxidized for 96 hours (GO 96 h) showed higher interplanar distance and also presented fewer layers when compared with GO oxidized for 2 hours (GO 2 h). The AFGO characterization points out that aminated groups were covalently attached to the GO nanosheets and the GO with a higher oxidation degree produced an AFGO with higher amination level. All synthesized nanosheets are highly stable when suspended at neutral media. The functionalized graphene suspensions were well modeled by the Hershel-Bulkley equation. The increase of the nanosheets concentration in the suspensions impairs the level of fluid structure and leads to a decrease in viscosity, yield stress, and elasticity. The GO 96 h promoted a lower decrease in viscosity, yield stress and elasticity than the GO 2 h suspension. In the case of AFGO, the greater amination degree can lead to a more pronounced drop in the suspension rheological properties. For the suspension with a higher concentration of GO 96 h, it was observed the appearance of hysteresis at low shear rates. These results show that small changes in the graphene functionalized nanosheets surface can significantly influence the rheological responses of a non-Newtonian fluid.

Os sistemas dinâmicos desempenham um papel crucial no que diz respeito à compreensão de fenômenos inerentes a diversos campos da ciência. Desde a última década, todo aporte tecnológico alcançado ao longo de anos de investigação deram origem a uma estratégia suportada por uma poderosa entidade, permitindo a inferência de modelos capazes de representar sistemas dinâmicos. A entidade citada que, por sua vez, pode aparecer de várias formas, é amplamente conhecida como dados. Além disso, independentemente dos tipos de sensores adotados a fim de realizar o procedimento de aquisição de dados, é natural verificar a existência de uma certa corrupção ruidosa nos referidos dados. Genericamente, a tarefa de identificação é diretamente afetada pelo cenário ruidoso previamente descrito, implicando na falsa descoberta de um modelo generalizável. Em outras palavras, a corrupção ao ruído pode ser responsável pela geração de uma representação matemática infiel de um determinado sistema. Neste trabalho, sob as perspectivas de classificação e regressão, demonstramos como a hibridação de várias técnicas de aprendizado de máquinas melhora a robustez ao ruído no que diz respeito à tarefa de identificação. Especificamente, no trabalho atual, mostramos o sucesso da estratégia proposta a partir de exemplos numéricos, tais como o crescimento logístico, oscilador Duffing, modelo FitzHugh-Nagumo, atrator de Lorenz e uma modelagem Suscetível-Infeccioso-Recuperado do SARS-CoV-2.

Os risers são importantes componentes na produção e exploração de petróleo e derivados. São responsáveis pelo transporte do óleo encontrado no reservatório até a Unidade Estacionária de Produção. A crescente demanda por esse produto faz com que a exploração seja feita em regiões com condições cada vez mais adversas. Tipicamente, um projeto deste porte exige um número muito grande de análises numéricas de elementos finitos e exigem uma experiência grande por parte do projetista a fim de obter uma solução viável. Esse desafio leva engenheiros a buscarem ferramentas consistentes e seguras que auxiliem nas etapas iniciais do projeto das configurações de risers e que sejam capazes de diminuir o número de análises totais exigidas. Uma dessas ferramentas éa utilização de métodos de otimização para obter de maneira consistente e segura os parâmetros que definem uma configuração. Este trabalho apresenta o método de Otimização Bayesiana, um método baseado em técnicas de aprendizado de máquina capaz de resolver problemas de otimização do tipo caixa-preta, de maneira eficiente, explorando o uso de aproximações analíticas da função objetivo, que se deseja otimizar. O método é aplicado em diferentes estudos de casos visando validá-lo como capaz de resolver uma gama de problemas de configuração de riser de maneira eficiente e consistente. Dentre os problemas aplicados estão diferentes configurações, diferentes casos realistas, mono-objetivos e multi-objetivos.

The motivation of this work consists in the use of four-roll mill in order to increase the phase separation of emulsions water-in-oil (W/O) present in the primary process of oil industry. With the mass and momentum conservation, the continuous phase is modeled by an incompressible, bi-dimensional and isothermal flow. Numerical simulations employing the finite element method were implemented to reveal the influence of the several flow configurations in the material mechanical behavior. From the obtained results, the standard way of classifying the flow in the four-roll mill according to the literature was proved inefficient. This work suggests local flow classifications for each position depending if it is occupied by the continuous or dispersed phase. The effect of the dispersed phase was described by a post-processing scheme. Microelements in shape of vectors were inserted in the domain and their deformations and pathlines were investigated. Thus, the deformation of droplets and their respective influences in the emulsion instability were analyzed.

O setor de refrigeração tem um papel fundamental e crescente na vida moderna. Seus usos vão desde processos industriais até conservação de alimentos passando por aplicações de conforto térmico. Tamanho peso econômico também se reflete nas emissões de poluentes do segmento. Acordos internacionais têm sido implementados com o objetivo de substituir fluidos refrigerantes que são danosos à camada de ozônio ou com alto potencial de aquecimento global. Essas iniciativas buscam mitigar o impacto direto da liberação desses gases na atmosfera. Outras abordagens visam reduzir o nível de emissões por meio da diminuição do consumo energético da operação desses sistemas. Melhorar as estratégias de monitoramento e controle de equipamentos residenciais de condicionamento ar certamente concorre para esse fim, uma que a venda desse tipo de aparelho cresce a cada ano. Neste trabalho investiga-se a implementação de uma estratégia baseada em controle preditivo para sistemas de expansão direta. O objetivo foi verificar ser o controlador baseado nessa abordagem seria capaz de controlar a temperatura interna de um ambiente. Os resultados mostram que controlador proposto foi eficaz em controlar a temperatura interna de em condições externas constantes, porém ineficaz em contexto dinâmico. Em paralelo, também se estudou a aplicação do método de estimação estados em horizonte móvel. Objetivo foi discutir se era possível inferir parâmetros do ciclo de refrigeração sem efetivamente mensurá-los. Os resultados mostram que o modelo de refrigeração proposto é observável em condições externas constantes. Além disso, uma comparação de desempenho entre o sistema operando com fluidos diferentes foi realizada. O intuito foi mostrar as possibilidades que método de estimação no que tange o cálculo de indicadores de desempenho.

In structure life prediction analysis occurrence of crack defects are of paramount importance to be considered. Basic studies using the Linear Elastic Fracture Mechanics approach shows that Stress Intensity Factor (SIF) parameter controls Fatigue crack growth (FCG). However, overloads may induce material “memory effects” that delay, arrest or accelerate FCG rate, a behavior not described by considering a single elastic parameter. To account for service variable amplitude loadings, the use of a prescribed stress-strain distribution has been proposed, in recent research studies, as the driving force of FCG using the critical damage approach. Due to the assumptions considered in the analytical derivations, the solutions provided by the method violate important equilibrium and compatibility conditions specially as using an idealized singular stress-strain field at the crack front region. In this work a comprehensive review of the recently published analytical results for solutions under the elastic and elastoplastic material behavior regimens are presented: Williams and HRR with singular stress distribution field and Creager-Paris with a non-singular stress field, but with material elastic approach. These solutions are compared, throughout the study, to results obtained from the numerical analysis using a 3D finite element discretization with elastoplastic von Mises yielding criteria employed for the modeling of the crack tip blunt and some comprehensive conclusions are derived regarding application limits of the theoretical models as well as the required extent of the numerical model representation. In addition to monotonic increasing applied loads, unloading was also considered in the numerical analysis, although no reference to this condition is present in the literature considering the analytical approach. For the numerical simulations considered elastoplastic algorithms were implemented in a plugin based framework.

Liquid dropout and accumulation in gas-condensate reservoirs, especially in the near wellbore region, hinder gas flow and affect negatively the produced fluid composition. Yet, condensate banking forecasting is commonly inaccurate, as experiments seldom reproduce reservoir extreme conditions and complex fluid composition, while most pore-scale models oversimplify the physical phenomena associated with phase transitions between gas and condensate. To address this gap, a fully implicit isothermal compositional pore-network model for gas and condensate flow is presented. The proposed pore-networks consist of 3D structures of constricted circular capillaries. Condensation modes and flow patterns are attributed to the capillaries according to the medium's wettability, local saturations and influence of viscous and capillary forces. At the network nodes, pressure and molar contents are determined via the coupled solution of molar balance and volume consistency equations. Concomitantly, a PT-flash based on the Peng-Robinson equation of state is performed for each node, updating the local phases saturations and compositions. For the proposed model validation, flow analyses were carried out based on coreflooding experiments reported in the literature, with matching fluid composition and flow conditions, and approximated pore-space geometry. Predicted and measured relative permeability curves showed good quantitative agreement, for two values of interfacial tension and three values of gas flow velocity. Following the validation, the model was used to evaluate wettability alteration and gas injection as prospect enhanced recovery methods for gas-condensate reservoirs. Results exhibited similar trends observed in coreflooding experiments and conditions for optimal flow enhancement were identified.

Com as novas regulamentações ambientais proibindo o uso de fluidos de perfuração que apresentassem em sua composição elementos como óleo diesel, se fez necessário a busca por novos componentes que, além de serem eficazes em sua função, apresentassem característica de biodegradabilidade sendo, assim, menos danoso ao meio ambiente. A base do fluido de perfuração escolhida para se enquadrar as novas regras foi estudada neste trabalho com objetivo de avaliar a sua estabilidade e o seu comportamento reológico. O fluido trabalhado é uma emulsão inversa cuja fase contínua é a olefina e a fase dispersa, salmoura de NaCl, acrescido de um emulsificante primário, Cybermul, e um secundário, Cyberplus. Para retardar o processo de desestabilização, majoritariamente caracterizado pela sedimentação, em um segundo momento, foi adicionada a cal hidratada. Os testes realizados foram divididos em duas etapas: a primeira avaliando emulsões preparadas a 5000 rpm, 7500 rpm e 10000 rpm sem cal em sua composição e a segunda avaliando emulsões preparadas com a mesma velocidade de rotação, entretanto com cal hidratada. A principal variável adotada foi a velocidade de rotação com objetivo de ratificar o conhecimento de que quanto maior a velocidade de agitação durante a emulsificação, menores são as gotas geradas e, consequentemente, mais estável a emulsão. Os testes foram compostos por tensão interfacial, bottle test, distribuição do tamanho de gotas (microscopia e espalhamento dinâmico de luz) e testes de comportamento reológico, que incluem, tensão constante, curva de escoamento e varredura de tensão e tempo. Os resultados obtidos permitiram relacionar o tamanho das gotas geradas com a estabilidade da emulsão e com a sua viscosidade.

Well planning presents a major challenge in the oil & gas industry. Wellbore stability analysis is one of the most important steps during well planning and provides the technical basis for a safe drilling operation. Therefore, to be able to comprehend and to predict the physical response of wells is extremely valuable to drilling engineers. In order to achieve wellbore stability, drilling operators control, among other factors, the internal wellbore pressure using the perforation mud weight. The mud pressure is used to balance the stresses at the wellbore region and its value must remain inside a stable window to ensure a safe operation. In this work, we propose a novel simulator to perform wellbore stability analysis and to compute the optimal mud pressure window. The simulator is based on a plugin architecture, which provides a more flexible environment to develop and to extend the simulator or even perform different analyses by exchanging one or more plugins. The ideal internal pressure and mud weight window are computed by solving a root-finding problem based on different failure criteria for the well. The simulator uses the Finite Element Method to solve each geomechanical analysis during the solution procedure, assuming the fluid-mechanical coupling and elastoplastic behavior around the wellbore. The simulator was validated using various examples and wellbore stability analyses were performed on a set of case studies from the literature.

During the last decades, large investments were made towards the development of numerical models and methods to forecast and analyze the different aspects of oil and gas recovery. In this context, modern simulators must be able to incorporate a wide variety of options to answer questions related to reservoir management accurately and effectively. In this work, we present a reservoir simulator based on a plugin architecture, where different formulations, solvers, and models can be developed, extended, and enhanced. With this approach, we use the black-oil model to implement traditional and state-of-the-art techniques, including fully- and adaptive- implicit methods, heuristic and PID time-step controllers, Newton-Raphson and Inexact Newton, and C1-continuous and conventional phase-potential two-point flux approximations. Several plugin configurations were tested and validated with commercial simulators, and their performances were used to determine which are the most suitable to solve multiphase flow problems.

Actuators based on piezoelectric materials are portable, manage small displacements while maintaining high resolution, frequency and stiffness. Those properties make them ideal for application such as acoustic transmission and micromanipulation. But inherent nonlinearities to those actuators, such as hysteresis and creep, greatly increases the challenge to control such devices. Furthermore, the increasing need for more precise and faster actuators, allied with frequent changes in the environmental conditions and complex dynamics related to the mechanical displacement induced by the excitation signal fur-ther worsens the problem. Analytical models are application-specific, mean-ing that they are not easily and efficiently scalable to all systems. Also, with increased complexity, the understating of underlying phenomena is not fully documented, making it difficult to develop such models. This work investigates those challenges from the perspective of the system identification methodology and data-driven models for piezoelectric actuators. Those are powerful tools that vary considerably according to the modeling criteria and the black-box approach is tested with experimental data acquired in a laboratory setting for a micromanipulator and acoustic transmission test benches. The results show good capabilities to compensate the nonlinearities for the micromanipulator, by predicting the hysteresis at different input frequencies. The acoustic trans-mission system was successfully modeled and although the results show that are still room for improvements, it provides insights in possible optimizations for the setup as the models here devised are useful for short prediction win-dows. Furthermore, this encourages further research in the pursuit of better data-driven abstractions for the piezoacoustic transmission and micromanipu-lation applications.

Due to their physical-chemical behavior, water-soluble polymers are used extensively in various phases of drilling, completion, workover, and production of oil and gas wells. Therefore, it is fundamental to predict and to control in-situ porous medium behavior in order to understand polymer performance. Experiments were conducted to study the degradation of a semi diluted (2000 ppm) aqueous solution of PEO, using two capillaries with different entrance radii (50 mm and 100 mm) both with 25 mm radius constriction, creating Fast-Transient Flows in their center. Different injection rates were imposed in order to observe different shear and extensional rates in the system. The effluent of the flow was collected, and reinjected, and rheological properties of the fluids were used as proxies for the degradation of the solution. We observed that for the more abrupt contraction, the minimum flow rate needed for degrading the polymer solution is lower. This result, when analyzed purely under shear rate perspective, is not reasonable, since the shear rates to which the polymer is subjected are equal at both capillaries.  Therefore, we inferred that the abruptness of the contraction plays a role in the degradation, which means elongational rate may be responsible for the lower critical flow rate. It was also observed a pattern for how the degradation occurs with subsequent injections.  We could infer that subsequent injections cause incremental degradation before approaching a stabilization plateau and that higher flow rates generated lower degradation plateaus.

O sucesso da operação de perfuração de poços de petróleo provém da incessante pesquisa científica que busca soluções e melhorias às diversas etapas deste processo. A correta formulação dos fluidos de perfuração e suas medições reológicas são parte fundamental neste cenário. O objetivo deste estudo é investigar o funcionamento de um viscosímetro Fann 35A e a metodologia utilizada para caracterização reológica de fluidos de perfuração, mediante comparação com um reômetro rotacional. Por conseguinte, diferentes fluidos, tais como soluções poliméricas e fluidos de perfuração, foram analisados com a finalidade de propor novas recomendações e aprimorar o uso do equipamento, comumente utilizado pelas companhias de petróleo. Curvas de escoamento e testes de força gel foram executados e, os resultados mostraram que as equações API são responsáveis pela determinação errônea dos parâmetros reológicos destes fluidos. Além disso, foi realizado um estudo teórico da relação entre a pressão e vazão em escoamento em tubo e espaço anular. Observou-se que a escolha precisa da função de viscosidade é imprescindível para correto dimensionamento de bomba. Ademais, uma geometria ranhurada foi projetada e desenvolvida a fim de evitar os efeitos do deslizamento aparente em baixas taxas de cisalhamento. A força gel, medida importante para avaliar o reinício do escoamento em poços de petróleo, também foi avaliada e apresentou resultados bastante divergentes dos obtidos no reômetro. O gel adicional exibiu respostas significativas e assim, sua utilização é indicada. Por fim, as melhorias e metodologias propostas mostraram-se promissoras, embora em alguns casos não possua valia pelo custo benefício. A incorporação das recomendações obtidas no presente estudo garante aquisição com maior acurácia da reologia de fluidos de perfuração e, consequentemente, desempenho correto de algumas funções atribuídas ao mesmo, evitando assim ocorrência de sérios problemas operacionais, ambientais e financeiros.

Reynolds Average Navier Stokes (RANS) models are among the most employed models to solve turbulent flows, due to their low computational cost. The majority of RANS models use the Boussinesq approximation, based on a linear relation between the deviatoric part of Reynolds stress tensor and the rate of strain tensor, with the turbulent viscosity as the positive proportionality parameter. However, these models fail in several situations, and a great deal of effort has been made by the scientific community aiming to improve model prediction through the development of non-linear models. Analysis of higher-order models employing objective orthogonal tensors has shown that these are very promising to improve the prediction of the normal components of the Reynolds stress. At the present work, turbulent characteristic velocity and length based on the intensity of the rate of strain tensor were examined, and new wall damping functions were developed presenting very good results. Further, a new one-equation turbulent model based only on the turbulent kinetic energy transport equation was proposed coupled with a new algebraic closure equation to model the turbulent kinetic energy dissipation. Very good results were obtained with the new developed models, when compared with DNS data and with the predictions of previous works.

O estudo reológico de hidratos vem se tornando cada vez mais importante graças à constante expansão da indústria de óleo e gás, principalmente em águas ultra profundas. O processo de formação de hidratos é uma grande preocupação, principalmente porque, em muitos casos, leva ao bloqueio total dos dutos de produção, causando interrupção na produção, além de perda de tempo de dinheiro. Hidratos são compostos cristalinos formados por água e pequenas moléculas de gás, em condições termodinâmicas de alta pressão e baixa temperatura. A fim de estudar este fenômeno, por analogia, compostos químicos que facilitam esta formação, à pressão atmosférica, como ciclipentano (CP) e tetrahidrofurano (THF) são utilizados. No presente trabalho, emulsões compostas por CP, óleo Primol, água deionizada e Span 80 (agente estabilizador) foram utilizadas e analisadas no reômetro Physica MCR301. A partir de uma perturbação térmica, pôde-se caracterizar vários parâmetros que influenciam a formação dos hidratos, tais como: taxa de cisalhamento, temperatura de indução, fração volumétrica da água, taxa de resfriamento etc. Além disso, foram realizados testes que avaliam a habilidade de reconstrução dos hidrato com o tempo e a existência de uma tensão limite de escoamento, a partir de testes oscilatórios.

O fenômeno de retorno de chama em tubos é conhecido e estudado há várias décadas. Sua análise clássica é baseada na determinação do gradiente de velocidade crítico, Gc, que o delimita como função das propriedades das misturas combustíveis. Entretanto, não é conhecido o efeito da diluição por CO2, importante para a previsão da segurança das instalações do pré-sal. Por isto são aqui desenvolvidos estudos específicos do retorno de chamas pré-misturadas em escoamentos laminares.
O objetivo geral deste trabalho é determinar experimentalmente a influência da diluição por CO2 sobre o retorno de chamas (flashback) em misturas de hidrocarbonetos (CH4 e propano) e de ar. O levantamento do estado da arte permitiu especificar as características da instalação experimental para o estudo deste fenômeno e, também, identificar as principais questões a serem abordadas. Foi projetado e construído um aparato experimental para o estudo do flashback em escoamentos laminares.
Os resultados originais obtidos mostram como a propensão ao retorno de chama é influenciada pela natureza do combustível, pela estequiometria da mistura e pela diluição. Misturas de propano possuem maior propensão ao flashback e maiores valores de Gc do que as de metano. Também foi mostrado que há uma redução da propensão ao flashback com o aumento da diluição. Esta propensão foi relacionada aos números adimensionais que caracterizam a combustão, isto é, os números de Lewis, Péclet, Karlovitz e Zel'dovich. Para este último, uma proposta original visando sua determinação é apresentada, que envolve uma expressão da taxa de liberação de calor da reação química global controlada por uma variável de progresso.

Microcapsules are applied in several sectors of industry when a physical barrier between the core material and the external environment is required. They protect their cargo and ultimately release it in a controlled way. In the present work, microcapsules with hydrogel-based shells are produced. Monodispersed microcapsules are formed by ionotropic gelation of gellan gum from monodispersed oil-in-water-in-oil (O/W/O) double emulsion templates obtained using glass-capillary microfluidic devices. An oil extraction step was added after the shell gelation process to enable the dispersion of the microcapsules in an aqueous medium. We report the operability window for the production of monodispersed microcapsules as a function of the flow rate of each fluid phase and the dimensions of the device. Microcapsules with mean diameters ranging from 95 to 260 μm and a maximum coefficient of variation of 5% were formed. The results show how to independently control the capsule diameter and shell thickness by varying the outer and middle phase flow rates. After that, we experimentally investigate the flow of monodispersed gellan gum microcapsules through a constricted capillary tube by measuring the evolution of the pressure difference and flow visualization. The maximum pressure difference and capsule deformation is obtained for capsules with different diameter and shell thickness. We map the conditions at which the capsule membrane ruptures during the flow, releasing its internal phase. Then, the gastro-resistance of gellan microcapsules is verified through an in vitro test that mimics the gastric and intestinal phases of digestion. Confocal fluorescence microscopy is used to track microcapsules integrity and we show that microcapsules cargo is released in the intestine mostly due to its pH. Finally, we demonstrate that it is possible to produce magnetic microcapsules with well controlled magnetic response by adding different amounts of ferrofluid to their core or shell. The microcapsules produced have great potential for different applications in food, biomedical, pharmaceutical and oil and gas industries.

Corrosion is a natural phenomenon that affects virtually every material, consumes billions of dollars annually worldwide and may result in catastrophic events. Coating is the most utilized technique to prevent corrosion. In this sense, epoxies are widely applied in all sectors of industry. Despite being high performance polymers, epoxies present issues, such as the use of toxic crosslinkers and long curing time. Benzoxazines are an emerging class of high performance self-curing thermosets. While typical monomers are based on mono- and difunctional derivatives, higher functional benzoxazines, prepared without using solvents, but purified with nontoxic solvents, are more desirable towards a greener synthesis of PBZ for high performance applications. Herein, we describe an environment-friendly approach to synthesizing a trifunctional benzoxazine from melamine, phenol, and paraformaldehyde. The chemical structure of the synthesized benzoxazine was confirmed by spectroscopic analyses, while the polymerization was monitored by thermal analysis. The synthesized benzoxazine monomer is, then, copolymerized with a commercial difunctional epoxy resin. We investigate the synergistic effect on improved thermal properties using differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). Results showed that a glass transition temperature of up to 268 °C was obtained. A higher thermal stability was also achieved with an onset degradation of nearly 400 °C and char yield of 22 wt% at 800 °C. Later, the copolymer was combined with silica to produce anticorrosion coatings. Results revealed superhydrophobic surfaces with great scratch resistance and adhesion to steel substrate. Electrochemical test proved high effectiveness of the coatings as anti-corrosion barrier, with increasing performance with the incorporation of silica.

The mathematical abstraction of a physical process is essential in engineering problems. There are many reasons for developing mathematical models of existing systems: i) for security and/or economic reasons, as it may be impractical or impossible to carry out experiments on the real system, ii) mathematical models are more flexible than physical prototypes, allowing a quick refinement of system designs to optimize various performance measures. The applications of the models can be divided into four parts, namely: design, estimation, control and monitoring. Some specific applications are i) simulations, ii) soft sensors, iii) performance evaluation, iv) statistical quality control and, v) fault detection and diagnosis. This work aims to: i) develop different classes of models capable of accurately simulating the output variable of a system, ii) evaluate the efficiency of optimization algorithms used in the parameter estimation task, iii) assess which friction model is the most appropriate to describe this phenomenon in a positioning system. The results show that combining models is an effective alternative to obtain more accurate simulation results. Nonlinear friction models were more adequate to describe this phenomenon, which had an asymmetric character, in the positioning system. In addition, the simulations performed with parameters estimated by the evolutionary algorithms showed better results. The decision tree-based optimizer, used in the second case study, proved to be equally effective compared to evolutionary algorithms.

The flow of capsules suspended in a liquid phase through small channels and capillaries poses a complex problem presented in different applications, from red blood cells on hemodynamics to flow in porous media. In porous media applications, the understanding of micro-scale dynamics is fundamental to assess the macroscopic flow behavior. Constricted channels and capillaries can be used to model a pore-throat connecting two adjacent pore-bodies. The flow of a suspended capsule through such models was analyzed to evaluate the flow characteristics, including the maximum pressure difference required to push a capsule through the constriction as a function of capsule radius, initial membrane tension, membrane material, channel and capillary geometries, as well as flow conditions. Inner and outer liquid phases were described by the Navier-Stokes equations and capsule membrane dynamics was modeled by a 1-D spring-like flexible structure. The fluid-structure interaction problem was solved using the finite element method coupled with the immersed boundary method. Results showed the mobility reduction of the continuous phase due to the presence of a capsule as it flows through the constriction. Such results can be used to design microcapsules to block preferential water flow paths in oil displacement process in porous media.

Wax deposition in petroleum production and transportation lines is one of the most relevant problems faced by the industry in order to assure the flow of oil and gas at the designed economical rates. Significant losses occur due to decreased production, line replacements and maintenance costs associated with cleaning operations. Wax deposit formation on the inner wall of the pipes might occur when the warm oil from the well loses heat to the cold environment, typical of deep water production operations, and its temperature reaches a critical value at which wax crystal formation occurs. This critical temperature is the WAT, wax appearance temperature. The prediction of wax deposit formation by simulation models is of fundamental importance for the proper design and operation of petroleum lines. Several wax deposition models have been developed and employed over the years, incorporating different wax deposition mechanisms. In the present work, the wax deposition phenomenon was studied employing a model fluid flowing through an annular test section, built to offer simple and well-defined boundary and initial conditions for the deposition process. Optical access to the interior of the test section allowed for the registration of images of the wax deposit formation. A miniature temperature probe was designed and employed to obtain original information on the temperature profiles within the deposit as it was formed. Also, the probe registered the transient evolution of the deposit interface temperature for different flow rates and cooling rates. The transient formation of a cloud of wax crystals over the deposit interface and carried by the flow was registered by a high frame rate camera. The temperature within this cloud was measured by the temperature probe. The results have shown that the deposit interface temperature evolves from a value equal to the WAT for the fluid measured by microscopy, rapidly reaching a constant value which is intermediate between the WAT and the solution thermodynamic phase change temperature. This information contradicts one of the key assumptions included in the molecular diffusion deposition mechanism, and widely employed in academic and industrial simulation models. The presence of wax crystals in the solution at temperatures above the WAT indicates that, locally, the WAT might differ from that obtained by microscopy, and that the local prevailing cooling rates could be lower than those imposed at the cold wall, allowing crystal formation above the deposit interface. The data obtained in the present work are of good quality and can be used as a reference to test wax deposition simulation models.

O escoamento trifásico na indústria do petróleo é caracterizado pela presença das fases gás, óleo e água. A presença da terceira fase (água) traz complexidade a esse processo, visto que pode provocar a formação de diversos novos padrões de escoamento, além dos já conhecidos para escoamento bifásico. Adicionalmente, a presença de uma fase líquida dispersa na outra pode formar uma emulsão, alterando significativamente a viscosidade e, assim, influenciando diretamente na perda de carga. O foco do presente trabalho é na previsão do padrão de golfadas com dispersões de água e óleo utilizando um modelo transiente unidimensional de Dois Fluidos. A presença da água é modelada através da solução da equação de conservação de massa para a fase água. Visando prever com precisão a queda de pressão, assim como a distribuição das frações volumétricas de cada fase ao longo do domínio, desenvolveu-se um modelo de fechamento algébrico para avaliar o escorregamento entre as fases líquidas. Com o modelo proposto, os resultados obtidos para a velocidade de deslizamento no escoamento água/óleo foram comparados com dados experimentais e de outros modelos disponíveis na literaturaa, apresentando um excelente desempenho. O modelo foi então utilizado para analisar o escoamento trifásico no padrão de golfadas. As previsões para a queda de pressão e características das golfadas (comprimento, frequência e velocidade de translação) foram comparados com dados experimentais da literatura e os resultados são promissores.

Recently ultrasonic guided waves have shown great potential for nondestructive testing (NDT) and structural health monitoring (SHM) in a damage evaluation scenario. Measurements utilizing elastic waves are particularly useful due to their capability to propagate in different materials such as solid and fluid bounded media, and, also, the ability to cover broad areas. When enough guided waves measurements are available and advanced data-driven techniques such as machine learning can be applied to the problem, the damage evaluation procedure becomes then even more powerful and robust. Based on these circumstances, the present work deals with the application of machine learning models to provide fault evaluation inferences based on ultrasonic guided waves information. Two main case studies are tackled in the mentioned subject. Firstly, a carbon fiber reinforced polymer (CFRP) plate is assessed using open data of Lamb guided wave signals in the detection of dot type defects. Results demonstrated that a baseline dependent approach can obtain excellent results when using system identification feature extraction. Secondly, corrosion-like defects in an aluminium plate are classified according to its severity. The methodology is assisted by a mode separation scheme of SH guided waves signals of pre-acquired data. The achievement on the matter has proved that the adoption of mode separation can in fact improve the machine learning results.

Devido aos efeitos nocivos dos combustíveis de origem fóssil, cresceu o interesse no uso de fontes sustentáveis de energia, em particular os derivados da biomassa. A gaseificação é um processo termoquímico de conversão de qualquer material orgânico carbonoso, occorendo a alta temperatura (geralmente entre 600 e 1000 °C) e com quantidade limitada de agente oxidante (vapor de água, ar ou oxigênio), que produz , carvões, bioleos e uma mistura gasosa, chamado gás de síntese ou syngas. O presente trabalho consiste no desenvolvimento de uma ferramenta de modelagem computacional baseado no equilíbrio termodinâmico. Assim, são avaliadas a flexibilidade, a confiabilidade e a estabilidade computacional das abordagens de equilíbrio estequiométrico e de minimização da energia livre de Gibbs para determinação das composições do syngas formado no processo. São investigados a influência da umidade relativa da biomassa (0, 20 e 40% em base mássica), da temperatura de gaseificação (800, 1100 e 1400 K), da razão ar combustível em relação a estoquiométria (0,2, 0,3 e 0,4) e da quantidade de oxigênio no agente de gaseificação (oxigênio puro, ar ou mistura equimolar dos dois) na composição do syngas, bem como nas eficiências energética e exergética, para biomassas compostas de carbono, oxigênio e hidrogênio. O melhor modelo é usado para otimizar as condições de gaseificação para biomassas com 20% de umidade relativa no ar, determinando as condições de concentração de agente oxidante e de temperatura que maximizem a eficiência energética. Esses dados, estão usados para desenvolver modelos matemáticos preditivos usando o conceito de análise de misturas. A robusteza desses modelos é avaliada por análise de variância (ANOVA) eeles estão usados para prever as condições de operações otimizadas da co-gaseificação de resíduos sólidos urbanos com bagaço de cana de açucar. Esses resultados estão comparados com as simulações realizadas em um modelo téorico considerando a cinética das reações. Observou-se que o metodo usando a minimização da energia livre de Gibbs apresentou os melhores resultados, mas os valores absolutos das eficiências estão fortemente dependente da aproximação feita nas correlações determinando o poder calorifico.

This thesis aims to investigate the cementing process of an oil well in presence of a filtrate loss zone. An experimental setup was designed, built and operated to simulate an annular well field-like geometry where cement would cure under controlled conditions. This well simulator consisted of an 8m concentric annular column with a section of semi-permeable external wall, equipped with pressure and temperature sensors. It allows the investigation of the cementing cure mechanisms, as well as the impact of the fluid loss zone in the pressure drop behavior. In order to predict this behavior on real oil wells, a 2D transient numerical model is proposed. A finite element model was implemented in Python with the aid of an open source library named FEniCs. Mass and momentum conservation equations are solved to obtain the pressure and velocity fields, and the cement mixture is considered an incompressible single-phase mixture composed of two chemical species: free-water and dissolved cement. The mass transport is modeled with an advection-diffusion equation and dissolved cement species is modeled as a viscoplastic fluid with shrinkage. Finally, the simulation results were confronted with the experimental data, and a good agreement was observed. Further investigation of numerical model parameters was performed, and a sensitivity analysis evaluated individual influence of those parameters in the pressure drop. The results indicate that pressure profile evolution has a strong dependency on thickening time, fluid loss flow rate magnitude and the Newtonian viscosity plateau for the evaluated shear rate regimen.

Polymers are used in a wide variety of applications and are present in all aspects of people’s daily lives; one of the main drawbacks related to their use, however, is their flammability. Major advances have been made in developing flame retardants (FRs) to be used in polymers, but most of the commercial FR additives in use today, mainly halogenated FRs, are extremely toxic both to humans and to the environment and are thus being increasingly regulated or banned throughout the world. Research into the use of renewable products as FRs for polymers has grown exponentially in the past two decades as a result, leading to the creation of many bio-based, environmentally-friendly FR systems. Acrylonitrile-butadiene-styrene (ABS) is one of the most widely used polymers in the world, but it is also one of the most flammable and one of the most challenging to flame retard due to its complex terpolymer structure. It is currently used commercially with toxic halogenated FRs, which face global restrictions or bans. Very little work has been performed on reducing ABS’s flammability using nature-derived or low-toxicity flame retardants. There is therefore an urgent need for the development of low-toxicity, bio-based FR systems to substitute halogenated FRs for commercial use in ABS.

The present work aims to contribute to the development of a bio-based flame-retardant solution for ABS. The first objective, a practical goal, is to develop a flame-retardant grade of ABS using mainly nature-derived additives, while maintaining acceptable mechanical properties. The second objective, having a scientific intent, is to comprehend the mechanisms of flame retardation of the bio-based FRs. The project consists of a two-phase study: (1) a Screening Phase, during which the effects of 8 different FRs and combinations thereof on ABS’s flammability are evaluated through microscale combustion calorimetry (MCC); and (2) a Detailed Analysis and Mechanistic Study Phase, in which the most promising candidates from Phase 1 are further analyzed through a variety of techniques with the objectives of better understanding their flammability and mechanical performances and of comprehending their mechanisms of flame retardation. Chapter 1 presents a brief introduction to the work followed by a literature review. Chapter 2 contains the experimental methodology used throughout the research. The Screening Phase is discussed in Chapter 3, which provides details about the experimental planning, preliminary tests, and the results of the screening experiments; the most promising samples and scientifically motivating synergies are identified. Chapter 4 contains the results of Phase 2, presenting in-depth discussions and hypotheses to explain the flame-retardation mechanisms of the bio-based FRs in ABS. Chapter 5 provides the main conclusions of the research and recommendations for future work.

Link da defesa:

https://cwru.zoom.us/j/99122416853?pwd=eGRDVzdKRXlvTi95aks3bDdRczVNdz09

Given the constant advances in technology, electronic devices have been going through a process of miniaturization while sustaining an increase in power. This trend proves to be a challenge for thermal management since commonly electronic cooling systems are air-based, so that the low heat transfer coefficient of air limits its capacity to keep up with the thermal needs of today’s industry. In this respect, two-phase cooling has been regarded as a promising solution to provide adequate cooling for electronic devices.
Two-phase thermosyphon loops combine the technology of two-phase cooling with its inherent passive nature, as the system does not require a pump to provide circulation for its working fluid, thanks to gravity and buoyancy forces. A micro-channel heat sink located right on top of the electronic device dissipates the heat generated. This makes for an energy and cost-efficient solution. Moreover, having a thermosyphon loop operating with a low GWP refrigerant such as R-1234yf results in low impact for the environment since it is an environmentally friendly refrigerant, and the system has low to none energy consumption.
This work provides a detailed numerical model for the simulation of a two-phase thermosyphon loop operating under steady-state conditions. The loop comprises an evaporator (chip and micro-fin heat sink), a riser, a tube-in-tube water-cooled condenser and a downcomer. Fundamental and constitutive equations were established for each component. A finite-difference method, 1-D for the flow throughout the thermoysphon’s components and 2-D for the heat conduction in the evaporator and chip, was employed. The model was validated against experimental data for refrigerant R134a, with good agreement. Additional simulations were performed to present a performance comparison between R134a and its environmentally friendly substitute, R1234yf.

A perfuração de poços marítimos em águas profundas e ultra-profundas é feita por sondas de perfuração flutuantes, que se conectam ao poço no leito marinho por meio do riser de perfuração. Essa conexão transmite esforços dinâmicos para a cabeça de poço e revestimentos, provocando fadiga nessa estrutura.
Levantou-se a hipótese de que as características dinâmicas do sistema mudam quando há aplicação de um carregamento elevado devido à deriva da sonda após um blackout, alterando o dano acumulado por fadiga.
A primeira etapa do trabalho envolveu a simulação estática deste carregamento extremo em um modelo solo-estaca em elementos finitos 3D. Com base nos resultados obtidos, realizou-se uma análise dinâmica para dois tipos de solo diferente, e para duas lâminas d’água diferentes.
Concluiu-se que o solo menos rígido sofre pouca influência sobre o dano acumulado, enquanto que para o solo mais rígido aumenta-se em cerca de 50% o dano percebido. A fadiga acumulada em lâmina d’água maior foi ordens de grandeza menor que para a lâmina d’água mais rasa.

O setor da refrigeração possui um papel essencial e crescente na economia global, com um aumento na quantidade de sistemas operantes. A necessidade de desenvolver novos refrigerantes tem sido cada vez mais frequente, a fim de atender a legislações ambientais cada vez mais rigorosas. Igualmente, medidas envolvendo a introdução de novos materiais, como os nanofluidos, tem sido uma constante. Neste trabalho, um sistema de refrigeração usando uma mistura nanolubrificante POE-diamante e refrigerante R410A foi simulado. Dados experimentais cedidos pela Universidade Federal de Uberlândia (UFU) foram usados para a elaboração e validação do modelo. O simulador utiliza a equação de Peng-Robinson para o cálculo das propriedades termodinâmicas e o método de fronteira móvel para a modelagem dos trocadores. O impacto das nanopartículas em relação aos parâmetros críticos foi avaliado a partir do princípio do isomorfismo e da natureza de ambos os materiais: fluido base e nanopartículas. A convergência da simulação do ciclo de refrigeração foi obtida com o método do simplex modificado mostrou-se adequado para tal aplicação, apresentando convergência satisfatória em todos os casos. As temperaturas de evaporação, condensação e de descarga do compressor são obtidas a partir das condições de operação do compressor, dos dois fluidos de transferência de calor, do grau de superaquecimento no evaporador e também do grau de subresfriamento no condensador. Superfícies de resposta foram criadas a fim de avaliar o efeito de cada uma das variáveis (temperatura de evaporação, frequência do compressor e concentração de nanopartículas) utilizadas no estudo s do coeficiente de performance (COP), da capacidade frigorífica e da potência do compressor. A temperatura de evaporação possui um impacto significativo sobre a capacidades frigorífica e o COP, enquanto que a potência é mais afetada pela frequência do compressor. A concentração de nanopartícula, apesar de possuir um efeito mais atenuado, não deve ser desprezada, devido à modificação que causa sobre as propriedades termofísicas da mistura.

Tradicionalmente, os testes de formação em poços de petróleo buscam caracterizar o campo de permeabilidades a partir da interpretação dos transientes de pressão (PTA) nos períodos de fluxo e estática baseados em modelos isotérmicos de escoamento em meios porosos. Com o avanço da instrumentação dos testes, registros mais precisos de temperatura passaram a estar disponíveis e fomentaram a pesquisa baseada em modelos não isotérmicos que possibilitaram a análise a partir dos transientes de temperatura (TTA). Além da caracterização de parâmetros do reservatório como permeabilidade e porosidade com a interpretação dos transientes de temperatura, os dados de pressão obtidos a partir de um modelo não isotérmico representa de forma mais fidedigna o fenômeno físico sobretudo quando os testes são submetidos a maiores diferenciais de pressão. Este trabalho consiste no desenvolvimento de um simulador para teste de formação que considera a modelagem não isotérmica de reservatório unidimensional radial acoplado a um poço produtor e na utilização deste simulador, associado a métodos de otimização multivariável, para resolução do problema inverso da caracterização de parâmetros do reservatório. Alguns métodos de otimização foram testados e o algoritmo do Simplex de Nelder-Mead apresentou melhor eficácia. Foram estabelecidos três tipos de problemas e utilizados em três casos hipotéticos considerando inclusive a imposição artificial de ruídos nos sinais de pressão e temperatura utilizados para resolução do problema inverso.

The optimized design of structures considering nonlinearities has been widely reasearched in the recent decades. The finite element analysis applied to topology optimization is jeopardized from excessive deformation of low-density elements under high compression, which hinders the process of finding an optimal solution. Two methods, Energy Interpolation and Additive Hyperelasticity, are implemented to overcome this difficulty in the minimum compliance problem, and hyperelastic material models are used to investigate their influence on the optimized topology. The Method of Moving Asymptotes is used to update the design variables whose sensitivities were calculated from the adjoint method. The state equation is solved through the Newton-Raphson method with an adjusting load step to reduce computational cost. Results for two benchmark problems are compared with those already estabilished in the literature. The use of different hyperelastic models presented little influence on the final design of the structure. The Energy Interpolation method was able to converge for much higher loads than the default method, while Additive Hyperelasticity did not show any significant improvements on plane strain.

O deslocamento de um líquido por outro em espaços anulares é comumente encontrado na indústria do petróleo, e a maioria deles envolve materiais não newtonianos. O espaço anular muitas vezes apresenta irregularidades causadas pela erosão, onde quantidades consideráveis de fluido de perfuração podem ser deixadas para trás durante o processo de deslocamento, comprometendo a qualidade da operação de cimentação. Motivados por esse processo industrial, testes de deslocamento entre líquidos a vazão constante foram realizados em espaços anulares cuja parede externa possui, em uma determinada posição axial, um aumento repentino de diâmetro seguido de uma diminuição repentina de diâmetro mais a jusante. O objetivo dos experimentos era determinar a eficiência do deslocamento em função da vazão, reologia dos fluidos e geometria da cavidade. Os resultados revelaram forte influência desses parâmetros na eficiência de deslocamento.
Ao mesmo tempo, um estudo numérico foi desenvolvido. Simulações numéricas das equações de Navier-Stokes em geometria axissimétrica para fluidos incompressíveis foram acopladas ao método Level-Set para captura da interface. Fluidos com viscosidade constante e o modelo newtoniano generalizado de Carreau-Yasuda foram utilizados. Isso permitiu simular deslocamentos entre dois fluidos newtonianos e entre um fluido newtoniano e outro não-newtoniano. Este foi utilizado tanto como fluido deslocador quanto como deslocado. Foram realizadas simulações para várias razões de diâmetros, viscosidades, tempos de relaxação, e números de capilaridade e de Reynolds. Identificamos quando a aproximação do espaço anular por duas placas paralelas pode ser aplicada e calculamos como a forma da interface depende dos parâmetros investigados.

A estimativa de produção de petróleo é um dos parâmetros essenciais para mensurar a economicidade de um campo e, para tanto, existem várias técnicas convencionais na área da engenharia de petróleo para predizer esse cálculo. Essas técnicas abrangem desde modelos analíticos simplificados até simulações numéricas mais complexas. Este trabalho utiliza Redes Neurais Artificias (RNA) para prever uma curva de produção de óleo que mais se aproxime da obtida por um simulador numérico. A metodologia consiste na utilização da rede neural do tipo feedforward para a previsão da vazão inicial e da curva de produção ao longo de dez anos para um poço produtor de óleo. Essa metodologia tem aplicação prática na área da exploração, visto que, nessa fase, ainda há muita incerteza sobre a acumulação de petróleo e, portanto, os modelos de reservatório tendem a não ser complexos. Os resultados foram obtidos a partir do treinamento de RNAs com dados coletados do simulador numérico IMEX, cujas saídas foram posteriormente comparadas com os dados originais da simulação numérica. Foi possível obter uma precisão de 97% na estimativa da vazão inicial do poço produtor de óleo. A previsão da curva de produção apresentou um erro percentual médio absoluto inferior a 10% nos dois primeiros anos. Apesar dos valores de erro terem crescido ao longo dos últimos anos, eles são menores quando comparados com a metodologia de declínio exponencial.

Two-phase intermittent pipeline flows are frequently found in industry, such as in the transportation of oil and gas. Their dynamics is commonly modeled through a one-dimensional approach, since long pipelines are employed, making 2D or 3D simulations unfeasible. The models always require an assessment of wall and interfacial friction factors, which are usually obtained by empirical correlations. In this work, from full 3D velocity fields across statistically developed slug unit cells, friction factors and momentum flux parameters were evaluated and new correlations for them are proposed. For that, several three-dimensional numerical simulations were performed by the software OpenFOAM. The Volume of Fluid (VOF) method and the k-ω Shear Stress Transport (SST) turbulence model were applied within a computational domain corresponding to a single slug unit cell, with periodic boundary conditions under imposed pressure gradient and overall void fraction. Liquid holdup, superficial gas and liquid velocities, slug parameters, turbulence quantities, liquid and gas velocity profiles, pressure distribution and wall shear stresses were also evaluated. Comparisons with experimental data and empirical correlations were performed, showing a quite good agreement. In summary, this work presents a methodology for assessing mean parameters of intermittent flows from three-dimensional fields, which may be further used to improve one-dimensional models.

A crescente demanda global por energia e a finitude dos recursos fósseis despertaram grande interesse pelo uso de energias renováveis e tecnologias menos poluentes. Neste contexto, este trabalho introduz uma simulação numérica de um sistema híbrido on-grid de uma unidade combinada de calor e potência (CHP) para microaplicações residenciais e industriais. O sistema é composto por uma célula a combustível tipo membrana polimérica (PEMFC) acoplada a um reformador de gás natural, painéis fotovoltaicos e baterias conectados à rede elétrica através de um inversor híbrido bidirecional. Uma análise energética foi desenvolvida para validar a rotina computacional e determinar a vazão de gás natural. Em seguida, foi realizada uma análise econômica baseada na evolução do fluxo de caixa dos usuários e no custo cumulativo total do sistema no horizonte de 2020 a 2040, de forma a investigar a influência das taxas de incremento das tarifas de energia elétrica e gás natural, diferentes configurações do sistema, do número de consumidores e do fator de aproveitamento de créditos na rede. Diferentes tipos de tarifa (convencional e branca) e a possibilidade de cogeração com o rejeito térmico da PEMFC também foram avaliados. Ao final, uma análise ambiental foi desenvolvida para avaliar a contribuição para o potencial de aquecimento global do micro-CHP. Paybacks entre 6 e 20 anos de operação do sistema foram alcançados para diferentes combinações dos parâmetros examinados considerando-se a adesão no ano de 2020. Adicionalmente, fortes reduções no custo cumulativo total foram obtidas levando-se em conta a queda prevista nos custos de aquisição dos componentes para as próximas décadas. Finalmente, emissões equivalentes até 30% inferiores às da eletricidade da matriz energética nacional e do fornecimento de calor por queima de gás natural foram calculadas com uso da cogeração no atendimento da demanda térmica dos usuários.

Com a necessidade de uma maior diversificação da matriz energética mundial e a inevitabilidade de investimento em fontes de energia limpa, a energia solar fotovoltaica se mostrou sem dúvidas a fonte energética mais promissora. Seu crescimento hoje é exponencial e se revela cada vez mais presente na matriz energética da maioria dos países. No Brasil, já se vê um mercado fotovoltaico sólido com geração de milhares de empregos e investimentos de bilhões de reais. Sua intermitência e imprevisibilidade é razão de grande pesquisa, no intuito de obter dados confiáveis que estimam a geração de um sistema solar fotovoltaico ao longo de seu ciclo de vida. Sendo assim, este trabalho conduz uma extensa revisão bibliográfica, passando pelos fundamentos básicos, até o desenvolvimento de um algorítmico robusto, onde é possível utilizar diversos modelos matemáticos presentes na literatura. Os dados de radiação solar, velocidade do vento e temperatura são obtidos de mapas gerados por 39 estações solarimétricas dentro e próximas das fronteiras do estado do Rio de Janeiro. O simulador também conta com uma extensa base de dados atual de módulos solares e inversores. Dessa forma, é possível simular e dimensionar sistemas solares fotovoltaicos para o estado do Rio de Janeiro, calculando indicadores de viabilidade econômico-financeira e ambiental e fluxo de caixa detalhado que refletem a viabilidade de um determinado projeto. Além disso, o simulador tem a capacidade de gerar mapas do Rio de Janeiro com dados destes indicadores em toda sua extensão.

Throughout the wax gelation, the network of paraffin crystals modifies the behavior of waxy crude oil. It changes from a low viscosity Newtonian to a high viscosity time-dependent material with yield strength. Now, it is totally challenging to find out the restart pressure for gelled crude oil flow with such a complex microstructure. Through my Master`s dissertation, we investigated two viscoplastic materials, namely a hair gel with a negligible thixotropy and a quite thixotropic 2% aqueous suspension of Laponite to mimic the startup flow of waxy crude oils. For both materials, the minimum axial pressure gradient required for the onset of flow was measured, and the measured values were in good agreement with the prediction of a conventional force balance. On the other hand, industry cases have exhibited that the just mentioned force balance leads to an overestimation of the minimum startup pressure gradient. In some studies, an elicited explanation is the thixotropic behavior of the gelled crude, but our above-mentioned results served to falsify it. Over the first part of my PhD thesis, we aimed to verify in the laboratory that why the force balance does not hold for gelled crude oil and then we sought a physically proper explanation for this discrepancy and also a reliable way to predict the minimum startup pressure gradient. Along these lines, we show the role of gelled crude oil's shrinkage in the discordance between static yield strength and required minimum pressure gradient to onset the flow, through rheometry and fluid flow in a tube. Then, we introduce a modified force balance equation with the role of shrinkage included to best estimate the minimum restart pressure gradient. Another essential element through the restart flow of gelled waxy crude is to find out a reliable strategy to mathematically model the material`s rheology. In most models that aim at predicting the rheological behavior of gelled waxy crude oil, the microstructure changes during flow are assumed to be thixotropic; reversible time dependent. But, we observed in our experiments with well-controlled flow and thermal histories that the irreversible character of time dependence is quite evident. Thus, in the second part of thesis we propose a model based on previous developments by Souza Mendes and co-workers that accounts for the irreversible time dependence observed experimentally in a waxy crude oil. The predictive capability of the proposed model is then assessed via comparisons with experimental data.

A camada geológica da costa brasileira conhecida como Pré-sal possui diversos reservatórios de petróleo cuja explotação encontra-se no início. Com o fim de explorar essas jazidas de forma mais segura e rentável, é essencial representá-las corretamente.Entretanto, não é uma tarefa simples. As rochas carbonáticas, típicas do Pré-sal, recorrentemente apresentam fraturas e carstes e o escoamento através dessas estruturas se comporta de forma diferente à da modelagem tradicional. Para solucionar o problema, desenvolveram-se técnicas que representam o fluxo através das demais mídias da rocha. No entanto, essas modelagens inserem diversas complexidades novas para a correta compreensão das formações e da previsão da produção. Nesse trabalho são analisados os impactos que as variáveis de reservatório têm no resultado das simulações com dupla permeabilidade.

A robótica móvel tem desempenhado um papel cada vez maior na sociedade moderna. Para cumprir suas tarefas, usualmente é necessário que os robôs móveis sejam capazes de navegar autonomamente pelo ambiente, partindo de um determinado local e chegando a um objetivo, evitando colisões. Uma etapa importante desse processo envolve a decisão dos caminhos a serem percorridos pelo robô ao longo de um mapa conhecido, que é executada por um programa chamado planejador de trajetórias. Muitos planejadores estão disponíveis na literatura, cada um com suas particularidades e aplicações. Esta dissertação apresenta um estudo comparativo entre cinco planejadores de trajetórias: Algoritmo de Dijkstra, Dynamic Programming, Probabilistic Roadmap (PRM), Rapidly-exploring Random Tree (RRT) e Hybrid A*, no qual são avaliados os tempos de processamento, os números de passos e os comprimentos das trajetórias resultantes. O estudo é executado por meio de simulações no ambiente MATLAB em três mapas diferentes, com objetivos diversos: o primeiro avalia a capacidade de se planejar uma trajetória em um ambiente complexo com múltiplos caminhos; o segundo testa a capacidade do planejador de evitar maiores obstáculos e lidar com passagens estreitas; e o último mapa, composto por longas retas e curvas abertas, tem o objetivo de avaliar o comportamento dos planejadores em trajetórias mais extensas. Cada planejador é submetido a séries de simulações em cada mapa, sendo que em cada uma alteram-se também suas características, para uma análise de sensitividade. Os melhores resultados de cada planejador para cada critério permitem identificar as principais vantagens e desvantagens de cada um.

Macroscopic rock properties such as porosity, permeability, and wave velocity describe petrophysical characteristics of rocks. In the particular case of heterogeneous rocks like carbonates, the estimation of equivalent petrophysical properties is challenging because of the presence of distinct rock structures in the same domain of interest. Also, the presence of voids with sizes of disparate orders of magnitude imposes extra difficulty to predict equivalent properties in carbonates. Equivalent properties can be determined through experimental measurements in the laboratory, which can be relatively slow and expensive, or via numerical simulation. In this study, we determined equivalent permeability of vuggy porous media using a Digital Rock Physics workflow, which is comprised of the microCT scan of the rock sample, the segmentation of the digital image, the mesh generation on the porous matrix and vuggy regions, and the solution of the Brinkman model using the finite element method. Digital Image Analysis (DIA) was used to determine the characteristics of the macropores (volume, surface area, pattern, among others), which were use to curve-fit the equivalent permeability resulting from solutions of the Brinkman model for many different samples. The goal was to be able to evaluate the equivalent permeability simply based on the image characteristics of the sample, without the need of solving the Brinkman model  

Wheeled-legged robots are an attractive solution for versatile locomotion in challenging terrain. They combine the speed and efficiency of wheels with the ability of legs to traverse challenging terrain. In general, the challenges with wheeled-legged locomotion involve trajectory generation and motion control for trajectory tracking. This thesis focuses in particular on the trajectory optimization task for wheeled-legged robots navigating in challenging terrain. For this, a motion planning framework is proposed that optimizes over the robot's base position and orientation, and the wheels' positions and contact forces in a single planning problem, taking into account the terrain information and the robot dynamics. The robot is modeled as a single rigid-body, which allows to plan complex motions for long time horizons and still keep a low computational complexity to solve the optimization quickly. The knowledge of the terrain map allows the optimizer to generate feasible motions for obstacle negotiation in a dynamic manner, at higher speeds. Such motions cannot be discovered without taking into account the terrain information. Two different formulations allow for either purely driving motions, where obstacle negotiation is enabled by the legs, or hybrid driving-walking motions, which are able to overcome discontinuities in the terrain profile. The optimization is formulated as a Nonlinear Programming Problem (NLP) and the reference motions are tracked by a hierarchical whole-body controller that computes the torque actuation commands for the robot. The trajectories are verified on the quadrupedal robot ANYmal equipped with non-steerable torque-controlled wheels in simulations and experimental tests. The proposed trajectory optimization framework enables wheeled-legged robots to navigate over challenging terrain, e.g., steps, slopes, stairs, while negotiating these obstacles with dynamic motions. 

Well formation tests are usually performed to determine the geological properties of a reservoir and the obtained data has often been interpreted based on an assumption that the reservoir is homogeneous in the vertical direction and described by a 1-D model. However, many reservoirs are found to be composed of different number of layers that have different characteristics. Production wells in such reservoirs may receive oil from more than one layer. In stratified reservoir system, the pressure and temperature behavior are not necessarily the same as in single layered system, and rarely reveals the same average properties of the entire system. The prediction of the characteristics of the individual layers is important to describe properly the reservoir and improve production management. This work presents a numerical transient-thermal model for a coupled wellbore/2D-reservoir considering Joule-Thompson heating/cooling, adiabatic fluid expansion, conduction and convection effects for both wellbore and reservoir for one phase fluid flow. The two-dimensional reservoir model allows the analysis of stratified zones and barriers. The model allows cross flow between the adjacent layers with different rock properties. Wellbore temperature and pressure at a certain gauge depth are evaluated numerically. Results show how pressure transient analysis (PTA) and temperature transient analysis (TTA) can be used to characterize different configuration of stratified reservoirs.

Recent advances in additive manufacturing techniques have increased their flexibility in making complex parts on a smaller scale. In this context, the design of porous microstructures has been standing out in the scientific community due to the ability to optimize the cell topology to meet the stress design requirements. However, there are several challenges that inhibit the fabrication of optimized parts obtained by the multi-scale topology optimization method, such as the connectivity of microstructures. The multiscale topological optimization consists of the optimization of both the macro-scale, global structure, and the micro-scale, microstructure of the material. The main objective of this work is to develop an effective scheme to guarantee compatibility in the transition between the different material microstructures obtained in multiscale optimization. The multiscale methodologies for simultaneous topological optimization of both scales and the homogenization procedures are described. The main numerical and computational aspects of these methods are presented, as well as representative examples to illustrate the capabilities of the proposed scheme. 

Neste estudo é apresentada a análise da interação da luz laser de alta potência com rochas nos processos de perfuração por fragmentação térmica, cujo objetivo principal do trabalho é estabelecer, através de um estudo experimental a possibilidade de utilizar a tecnologia a laser na perfuração de rochas duras, tais como o granito, quartzo entre outros. Encontrando sua motivação na aplicação de novas tecnologias para desenvolver ferramentas que melhorem a eficiência no processo de perfuração de poços (maior taxa de penetração), e conseguir retirar a maior quantidade de material fragilizado pela irradiação do laser. Para isto, foi necessário fazer uma montagem opto – mecânica envolvendo um sistema de limpeza que ajudasse a remover o material fragilizado pela ação do laser, usando-se gás nitrogênio a alta pressão. Foram analisados diferentes intervalos de tempo de exposição da radiação do laser nas rochas, para avaliar a taxa de volume retirado, a quantidade de energia especifica requerida para perfurar diferentes materiais, em função da potência e do tempo de irradiação do laser. Para este fim, foi necessário montar um sistema de limpeza que permitisse a remoção do material fragilizado pela ação do laser. Seguindo a revisão bibliográfica na seleção do material foram escolhidos três tipos de materiais (granito branco, granito cinza e travertino) conhecidos comercialmente no Brasil e que tem similitude com as rochas encontradas nos reservatórios do Pré-sal. A partir desta seleção foi indispensável conhecer a composição química dos materiais através de fluorescência, raios X, microscopia de varredura (MEV). Além disso, estudou-se o comportamento termomecânico das rochas carbonáticas a traves de analises térmicos (Termogravimetria e Analise Térmico Diferencial), para identificar e compreender os fenômenos envolvidos no processo de perfuração. Os resultados obtidos são analisados para parametrizar as variáveis em consideração, melhorando as condições do processo de perfuração por fragmentação térmica dependendo do material estudado.

O uso de juntas cardânicas ativas é restrito pela capacidade de torque de pequenos motorredutores e, atualmente, os dispositivos embarcados são obrigatórios para as aplicações robóticas. Nesta tese, um novo conceito de manipulador robótico acionado por juntas cardânicas com três graus de liberdade foi introduzido. O controle dinâmico é essencial para estudar as limitações desse dispositivo, portanto, o objetivo deste estudo foi controlar a junta cardânica ativa de três graus de liberdade usando simulações numéricas e experimentais. O manipulador foi projetado com apenas uma junta cardânica para que a sua cinemática e dinâmica sejam exploradas; por esse motivo, a junta foi construída com sensores de carga na base e sensor de unidade de movimento inercial na parte superior do efetuador do manipulador. Além disso, foram fabricadas três placas de controle: a primeira foi projetada para controlar os três acionamentos dos motores de passo; a segunda, para ler o sensor da unidade de movimento inercial; e a última, para ler os sensores de carga. Quatro problemas foram descritos para testar os limites deste dispositivo, analisando, além da cinemática e dinâmica, o atrito do rolamento, a identificação da folga e o torque do impacto. O primeiro problema mantém a posição do efetuador do manipulador constante enquanto transmite rotação entre os eixos. O segundo problema, o efetuador recebe um caminho planejado, por exemplo, um círculo, mas não transmite rotação entre os eixos. O terceiro problema é a combinação dos movimentos anteriores, em que o efetuador transmite rotação entre os eixos, enquanto segue por um caminho planejado. Para o quarto problema: uma nova abordagem é aplicada para mover o efetuador de um ponto para outro usando rotação cônica. Uma nova notação dinâmica é introduzida para esclarecer e simplificar a matemática usada para descrever a dinâmica de um corpo rígido 3D para resolver a cinemática e a dinâmica inversa, sendo os resultados numéricos do modelo físico comparados com os experimentais para definir os limites de aplicação da junta.

Em poços de petróleo, uma boa cimentação é fundamental para garantir o suporte estrutural do revestimento e isolar o poço de águas subterrâneas. A alteração das propriedades e o surgimento de defeitos no cimento podem levar a vazamentos e danos ambientais. Durante a desativação e o abandono do poço em operações conhecidas como P&A (Plug & Abandonment), é necessário avaliar a condição do cimento assim como a identificação de possíveis defeitos.
Desde 1950, ferramentas de perfilagem acústica capazes de analisar a condição do cimento vem sendo desenvolvidas e aprimoradas. As ferramentas mais famosas foram a “Cement Bond logging Tools” (CBT), baseada na técnica de operação CBL (Cement Bond Logging) e VDL (Variable Density Logging), e a “Cement Evaluation Tools” (CET). Contudo, essas ferramentas foram projetadas para operar somente em configurações com um único revestimento e sem a coluna de produção. Remover a coluna de produção é um processo caro, por isso, um desafio na indústria que permanece até os dias atuais é a criação de uma ferramenta capaz de operar em poços sem a necessidade de se remover a coluna de produção visando a redução dos custos do processo, ou em poços que possuem múltiplos revestimentos.
Nos últimos anos muitas patentes e artigos tem sido publicados abordando esse mesmo tema. Este trabalho possui como motivação analisar a eficiência da utilização de ondas ultrassônicas guiadas em poços com multicamadas para detecção de defeitos presentes na camada de cimento. Para isso, primeiramente se realizou uma vasta pesquisa bibliográfica buscando identificar os tipos mais frequentes de defeitos existentes, o princípio de funcionamento das principais ferramentas comerciais existentes e os avanços das pesquisas na fronteira do conhecimento desse assunto.
Modelou-se a partir da teoria de ondas e de suas equações analíticas o poço estudado, encontrando as equações de dispersão que regem o problema. Desenvolveu-se um código analítico-numérico capaz de resolver essas equações e de obter as curvas de dispersão.
Além disso, analisou-se como os defeitos e suas variações afetam as curvas de dispersão, identificando os comportamentos característicos para cada caso. Posteriormente, a distribuição da energia, obtida a partir do vetor de Poynting no poço, foi usado afim de se investigar os modos que concentram a energia na camada de cimento. Para o caso sem a coluna de produção, concluiu-se que os melhores modos de excitação estão em alta frequência, pois a energia se concentra no cimento. Já para o caso com a coluna de produção, encontrou-se que na frequência de 30,13 kHz a uma vagarosidade de 726,5 us/m a energia está praticamente toda concentrada no cimento o que é interessante para se analisar defeitos apenas no cimento e que em 39,17 kHz e 507 us/m a energia está distribuída entre o cimento e a coluna de produção, sendo interessante para analisar defeitos que envolvam ambos.
 

Após uma jazida de petróleo ser encontrada, a produção de óleo ou gás é feita através de um poço produtor que é perfurado até atingir as camadas de rocha onde os hidrocarbonetos estão alojados. Com a constante produção, a pressão de reservatório decresce até atingir um nível que é  insuficiente para o aproveitamento econômico. Geralmente utiliza-se a injeção de água para manter o nível de pressão do reservatório. Nos estudos de um reservatório de petróleo é fundamental o conhecimento de propriedades básicas da rocha e dos fluidos nela contidos. São essas propriedades que determinam as quantidades de fluidos existentes no meio poroso, a sua distribuição, a capacidade desses fluidos se moverem e,mais importante, a quantidade de fluidos que pode ser extraída. Através do método convencional de injeção de água objetiva-se a manutenção da pressão do reservatório e o deslocamento de óleo em direção aos poços produtores. A água (fluido deslocante) tende a ocupar gradualmente o espaço antes ocupado pelo óleo (fluido deslocado), contudo, por efeitos capilares, uma parcela do óleo não é retirada do meio poroso configurando o que chamamos óleo residual. Em função da razão de mobilidade da água e do óleo, a frente de deslocamento não é uniforme, e um grande volume do reservatório não é atingido pela água de injeção. A adição de polímero à água de injeção visa o aumento da viscosidade da água, e assim, melhorar a razão de mobilidade água-óleo, aumentando a eficiência de varrido uma vez que uniformiza a frente de avanço, reduzindo a formação de caminhos preferenciais no reservatório. Além de diminuir a razão de mobilidade, soluções poliméricas podem contribuir para um melhor deslocamento de óleo em escala de poro, a partir de seu efeito elástico, reduzindo, portanto, a saturação de óleo residual. Contudo, tal mecanismo em micro-escala, ou seja, em escala de poro não é totalmente compreendido. O presente trabalho preocupa-se principalmente em analisar o fator de recuperação do óleo e saturação de óleo residual após processo de deslocamento de óleo por água salgada, solução polimérica de poliacrilamida parcialmente hidrolisada (HPAM) e solução de glicerina em testemunhos de Arenito Bentheimer. Um porta-testemunho especial foi utilizado para a realização dos testes de deslocamento, sendo monitoradas a variação de pressão ao longo da amostra, além dos volumes de injeção e produção de fluidos em função do tempo.

 

O crescente uso de modelos integrados de produção (MIP) na indústria de petróleo como solução para representar o potencial de uma jazida se justifica pelos atuais cenários encontrados para o desenvolvimento de novos projetos, caracterizado por baixos preços de venda do petróleo e conceitos de produção de alta complexidade. Esta abordagem, onde diferentes partes de um sistema de produção são integradas, permite ao usuário um entendimento detalhado das interações entre reservatório, poços e rede de escoamento, e facilita a detecção de gargalos e consequentemente a otimização do plano de explotação.
Neste contexto, é fundamental obter uma modelagem satisfatória do fluido em todo o sistema de produção. O modelo deve honrar tanto o escoamento no meio poroso, isotérmico, quanto o escoamento nos poços e dutos, que precisa ser caracterizado em várias temperaturas. Além disso, o modelo deve ter tempo de simulação adequado.
Uma maneira criteriosa de modelar as propriedades de um fluido é através do ajuste de uma equação de estado (EOS). Uma EOS detalhada com 24 componentes determinada por cromatografia gasosa e EOS simplificadas com 14, 9, 7, 6, 5 e 4 pseudocomponentes foram geradas para avaliar este problema. As EOS foram usadas para representar as propriedades PVT em um MIP e ao final foram comparadas as respostas das EOS simplificadas e detalhadas, a fim de estabelecer resultados adequados com um tempo computacional adequado.
Os resultados obtidos mostram que o uso de EOS excessivamente simplificadas, apesar da melhoria no tempo computacional, podem gerar resultados insatisfatórios em modelos integrados de produção.

O estudo de bolhas de separação laminar em aerofólios que operam a baixos números de Reynolds é importante para diversas aplicações tais como: turbinas a gás, geradores eólicos, veículos aéreos não tripulados, dentre outros. A perda de performance desses equipamentos, normalmente, está associada a presença e a ruptura de bolhas de separação. Neste trabalho buscou-se investigar, experimentalmente, os efeitos relacionados a variação súbita no nível de perturbações do escoamento à montante da região de formação da bolha. A ideia é avaliar como o escoamento se comporta em uma situação onde o nível de turbulência do escoamento livre pode variar. No caso de turbinas, essa condição pode ser observada quando as esteiras dos aerofólios de um estágio influenciam àquelas do estágio seguinte. Apesar da relevância prática, o regime transiente foi pouco investigado na literatura. Para o presente estudo foi utilizado um canal de água do laboratório de Engenharia de Fluidos. O gradiente de pressão, necessário para a formação da bolha, foi gerado com uma placa convergente-divergente e o nível de perturbações na camada limite foi controlado através de um gerador do tipo fita vibrante. Os campos de velocidades na região da bolha foram medidos com a técnica de Velocimetria por Imagem de Partícula com alta resolução temporal. A geração de perturbações e as medições de velocidade foram sincronizadas, permitindo assim o uso de técnicas de extração de médias de eventos. Caracterizou-se a bolha em um regime estacionário, onde não houve excitação de perturbações controladas a montante da bolha. Os resultados foram comparados com referências da literatura, apresentando boa concordância. No regime transiente, analisou-se o escoamento desde o instante em que o gerador de perturbações foi subitamente desligado até a recuperação completa da bolha. Nessas condições, notou-se que bolha cresceu até atingir um comprimento maior do que o observado no regime estacionário. Somente após o desprendimento de um grande vórtice, a bolha voltou a exibir características similares àquelas do regime estacionário. Esse comportamento se assemelha àquele reportado na literatura para ruptura de bolhas (“bursting”) mas que ainda não é bem compreendido.

A produção de gás para os próximos anos continuará em ascensão devido ao aumento da demanda energética ao redor do mundo. A tendência global é a busca por energias mais limpas que o óleo e o carvão. O gás aparece como uma ótima opção pelo seu baixo custo e por ser menos poluente.
Gás em reservas mais profundas podem estar presentes em condição de gás retrógrado, a uma temperatura entre a temperatura crítica e cricondentérmica. O gás retrógrado em alta pressão e alta temperatura apresenta apenas a fase gasosa. Entretanto durante a produção, a pressão é reduzida na região mais próxima ao poço. Caso a redução da pressão atinja a pressão de orvalho da mistura, ocorrerá a formação de condensado. O acúmulo excessivo deste condensado ao redor do poço pode reduzir a permeabilidade efetiva do gás prejudicando a sua produção.
Para estudar a formação de condensado e como ela implica na redução da permeabilidade do gás e do líquido, foram realizados experimentos utilizando o Método do Regime Pseudo-Permanente em duas amostras de rochas em condições de alta pressão e alta temperatura.
Ao todo foram realizados 4 experimentos, sendo três utilizando uma amostra de arenito e o último uma amostra de carbonato. O fluido utilizado para representar o gás retrógrado foi uma mistura de etano/n-pentano em três níveis de formação de condensado. Os resultados obtidos mostraram que a permeabilidade relativa do gás e do líquido foram influenciadas pela quantidade de líquido produzido, vazão de gás e características da amostra.

A migração, induzida por cisalhamento, de partículas em suspensão tem sido o objeto de diversos estudos no passado, e alguns modelos já foram propostos a fim de prevê-la. A migração de gotas em emulsões, contudo, apresenta uma série de especificidades, especialmente nos casos em que fenômenos de floculação e coalescência estão presentes. Emulsões de água-em-óleo com surfactantes não-iônicos, tal qual o Span 80, tendem a formar agregados de gotas quando se encontram sob cisalhamento, e esses agregados podem afetar o deslocamento de partículas em um campo cisalhante não-homogêneo. Testes em diferentes taxas de cisalhamento e intervalos de tempo foram realizados em uma geometria transparente do tipo placa-placa na intenção de avaliar macro- e microscopicamente os processos de floculação e migração de gotas de água em uma fase contínua de óleo. Estes experimentos apresentaram uma rápida e intensa migração radial de gotas rumo ao centro da placa em altas taxas de cisalhamento, e uma formação significativa de flocos em taxas de cisalhamento mais baixas. O tamanho médio e concentração de gotas foram determinados para diferentes posições da placa ao longo do tempo.

The majority real-world engineering design, management and decision making problems, may be seen as a multi-objective optimization problem, while others may comprehend classification, prediction or regression tasks. The use of artificial intelligence (AI) techniques to address these problems is attractive since they i) are most recommended for MO problems when compared with classical mathematical programming methods, which requires significant computational cost and informations of the problem domain; and ii) present better results with a simpler structure, adaptability and the possibility to extract information from the problems, when compared with other methods. Therefore, this work considers the application of AI techniques in anti-lock braking system (ABS) improvement, heat exchanger design optimization and theft-sensitive leak detection system (LDS) problems to i) develop and improve AI techniques; ii) constitute an AI application guide for control and design MO problems, as well as for the feature extraction process and machine learning classifier training through supervised learning approach; and mainly iii) demonstrate AI methods potential in solving realworld mechanical engineering problems. The results showed that the novel AI optimization algorithm versions proposed performed well considering the Pareto front dominance, spacing, Euclidian distance and hypervolume for the ABS problem and two IGD-based metrics compared through wilcoxon rank sum test for the HE design problem. Also, the decision tree-based machine learning classifier, trained through supervised learning and principal component analysis feature extraction technique, detcted all fuel theft situations with no false alarms.

Ao longo do tempo de vida de um poço de petróleo, a produção de óleo pode sofrer uma variação significativa. Sendo assim, esta que, anteriormente, era composta em sua maioria de petróleo pode acabar se tornando em uma mistura de água e óleo, gerando a ocorrência de emulsões no sistema. Contudo, a presença dessas emulsões torna-se um problema recorrente para a indústria de petróleo, uma vez que acarretam em um aumento de custos de transporte e processamento do óleo.
Sendo assim, a verificação da eficiência de separação primária de emulsões ocorre, normalmente, com a utilização de emulsões reais retiradas dos campos. Todavia, tais emulsões são fluidos que podem apresentar comportamentos tóxicos, voláteis e inflamáveis [1]. Somado a isso, a falta da confiabilidade necessária para a aferição da eficiência dessas tecnologias de separação uma vez que há uma extrema dificuldade de repetibilidade das emulsões sintetizadas a partir do óleo cru, induz à necessidade de utilização de fluidos modelos que possuam comportamento similares ao petróleo.
O objetivo geral deste trabalho é, sobretudo, aprofundar o estudo das emulsões de petróleo e desenvolver um método para a síntese de fluidos de trabalho com características reológicas, de formação e de estabilidade semelhantes às de óleos e emulsões reais encontrados no campo. Assim, para o desenvolvimento deste trabalho, óleos modelo com características macroscópicas similares a dois diferentes petróleos (A e B) foram sintetizados. Verificou-se que as emulsões modelo obtiveram um comportamento mais similar às emulsões formadas a partir o petróleo A. Contudo, tais sistemas apresentaram um comportamento de shear-thinning no regime semi diluído, devido, sobretudo, a diferenças qualitativas na estrutura dos agentes tensoativos das emulsões de petróleo e os surfactantes comerciais utilizados em emulsões modelo.

This research focused on studying the rheology of W/O emulsions formed by different crude oils, with special attention to their behavior when subjected to Brownian and hydrodynamic forces. At a microscale level, particle-particle and particle-medium interactions, both of which are involved in the phenomenon of aggregation, define the emulsion’s rheological behavior due to the complex structures created by the droplets as their concentration rises. An experimental study was performed in order to visualize the emulsions characteristics according to its disperse phase concentration and the shear rate applied, verifying the existence of both coalescence and flocculation phenomena during shearing. The modeling of interactions between droplets allowed the prediction of the emulsion’s behavior by colloidal thermodynamics. With it, a methodology to fit the experimental data that included both the shear rate and disperse phase concentration parameters was proposed. This work emphasizes the importance of studying emulsion systems at a micro-scale level in order to obtain a better comprehension of the formation and breakage processes of complex and random aggregate structures. This allows the prediction of the emulsion’s rheological behavior, and the proposition of a phenomenological model to best describe it. risisualize the emulsions.

This work presents selected applications of topology optimization for non-Newtonian fluid flow problems using the virtual element method (VEM) in arbitrary two-dimensional domains. The objective is to design an optimal layout into a fluid flow domain to minimize dissipative energy governed by the Navier–Stokes–Brinkman and non-Newtonian Carreau–Yasuda model equations. The porosity approach proposed by (Borrvall and Petersson, 2003) [1] is used in the topology optimization formulation. To solve this problem numerically, the recently proposed VEM method is used. The key feature that distinguishes VEM from the standard finite element method (FEM) is that the interpolation functions in the interior of the elements do not need to be computed explicitly. This is because the integration is on lower-order polynomial and basis functions, and there is great flexibility by using a nonconvex element. Therefore, the computation of the main element matrices and vectors are reduced to the evaluation of geometric quantities on the boundary of the elements. Finally, several numerical examples are provided to demonstrate the efficiency of the VEM compared to FEM and the applicability of the topology optimization to fluid flow problems.

Infrared thermography has been used as a nondestructive evaluation (NDE) technique to detect flaws in structural components, playing an important role in manufacturing inspection, in-service inspection and maintenance programs.
An investigation program was launched with the objective of presenting combinations of analytical, experimental and numerical methods to predict and monitor fatigue initiation and fatigue damage progression in equipment such as pressure vessels, tanks, piping and pipelines with dents or complex-shaped anomalies.
The monitoring of fatigue initiation and propagation in the actual specimens used nondestructive inspection techniques such as thermoelastic stress analysis (TSA), three-dimensional digital image correlation (3D-DIC) and fiber optic Bragg strain gages (FBSG) to determine strains at fatigue hot spots locations. Strain fields determined from the experimental measurements and from the finite element method (FEM) were combined with the fatigue Coffin-Manson strain-life equation and the Miner’s fatigue damage rule to predict fatigue life (N).
Results from tested 3m long tubular (with nominal dimensions: 324mm external diameter and 6.35mm wall thickness) specimens containing complex-shaped dents were reported and fully analyzed.
This work confirmed that infrared rapid fatigue assessment methods are practical and efficient tools that can provide a reliable, non-destructive and faster results about the fatigue behavior of materials.
Good agreement among fatigue life estimations and actual fatigue lives of complex dent shapes in pipeline specimens can only be achieved if accurate measurements or numerical determinations of the circumferential strains actuating at the dent hot-spots were coupled with suitable fatigue strain-life curves. Coupling an experimental technique for accurately determining dent shapes to a numerical strain analysis technique will lead to good hot-spot strain estimations to be combined with the strain-life fatigue curves. The present conclusions can be applied to other structures that may present dents such as tanks and pressure vessels.

A industrialização e a crescente preocupação com o meio ambiente geram, cada vez mais, a busca por fontes de energia que emitam menos gases efeito estufa. A biomassa, devido a sua grande ocorrência ao redor do mundo e a sua diversidade, é uma forte alternativa aos combustíveis fósseis. Sua gaseificação gera um combustível gasoso chamado syngas. A problemática no manejo dos resíduos sólidos urbanos (RSU) e a grande disponibilidade do bagaço de cana-de-açúcar no Brasil, resíduo sólido agrícola, fizeram deles tipos de biomassa de interesse para este trabalho. Objetivou-se modelar no MATLAB ® a gaseificação de biomassa a partir de uma abordagem cinética para um gaseificador de leito fixo cocorrente e para o ar como agente gaseificador. O modelo foi validado com dados experimentais e numéricos da literatura e aplicado à simulação da co-gaseificação de RSU e bagaço de cana-de-açúcar, na qual a razão de co-gaseificação (RCG) representou a percentagem de RSU na biomassa de entrada. Um planejamento de experimentos composto central com 3 fatores e 3 níveis foi realizado, resultando em 27 ensaios variando os fatores RCG, umidade da biomassa e razão de equivalência. Foram criados modelos polinomiais para a composição do syngas obtido, para o PCI do syngas, para a eficiência energética do processo e para a soma das frações molares de CO e H2 em base úmida. Os modelos foram considerados robustos, com valores de R2 e R2 ajustado variando de 0,96082 a 0,99345 e 0,94007 a 0,98998, respectivamente. O impacto dos fatores escolhidos nas respostas foi analisado, e os modelos de eficiência energética e soma das frações molares de CO e H2 foram maximizados. O caso otimizado, com RCG 7,98%, umidade 5,00% e razão de equivalência 0,18, resultou em um syngas de composição 3,72% H2O, 29,68% CO, 7,87% CO2, 19,07% H2 e 0,80% de CH4 em mol, correspondendo a um PCI de 6,56 MJ/Nm3 e uma eficiência energética de 37,66%, Por fim, o processo demonstrou bom potencial para geração de um gás rico em CO e H2.

O objetivo dessa dissertação é apresentar a simulação termodinâmica de uma planta de etanol derivado de cana-de-açúcar, abrangendo o sistema de cogeração e hidrólise enzimática aplicáveis ao bagaço, e associada a uma planta de cultivo de microalgas para produção de biodiesel de terceira geração, por rota de transesterificação. Adicionalmente, foram avaliados estatisticamente, por planejamento composto central 32, os efeitos de duas variáveis – designação de bagaço para hidrólise enzimática (x1) e de etanol para desidratação (x2) – sobre parâmetros de performance energéticos e ambientais. Em um cenário base similar à realidade brasileira – coordenadas (0, 50%) –, observou-se que a integração pode produzir uma quantidade de biodiesel (26,6 kg.t-1 de cana) equiparável, em termos energéticos, à oferta de etanol (62,2 kg.t-1 de cana) na saída da planta. Nesse cenário, a demanda energética da seção de microalgas-biodiesel demonstrou-se semelhante à da planta de etanol tradicional. Os modelos produzidos demonstram-se robustos e precisos (R2 superiores a 0,998). As avaliações estatísticas indicaram que x1 afeta quadraticamente as respostas, e x2 promove efeitos nulos ou brandos. Destaca-se valores de NER entre 0,763 e 0,904, próximos a 1, e valores de FER entre 4,33 e 10,88, acima do limite mínimo de 3 para se determinar o caráter renovável do sistema. As otimizações das respostas para identificação de coordenadas ótimas, segundo cada parâmetro de performance definido, foram realizadas, encontrando-se um valor ótimo para eficiência exergética global nas coordenadas (92 %, 0). Por fim, a associação das plantas de etanol e de biodiesel forneceu resultados promissores, com perspectivas interessantes para eficiências energéticas, conversão de energia fóssil em renovável, e incremento de eficiência exergética original da planta de etanol.

Desde os primórdios da Revolução Industrial, caldeiras têm percorrido um longo caminho de desenvolvimento visando melhor segurança e maior eficiência. Em particular, o mercado de caldeiras flamotubulares tem exigido cada vez mais projetos personalizados dos fabricantes, abrindo caminho para o uso de modelos de caldeiras como eficientes ferramentas de projeto. Neste trabalho, foi desenvolvido um modelo termodinâmico para uma caldeira flamotubular, de três passes, operando em regime permanente com gás natural boliviano. Foi testada a sensibilidade no perfil de temperatura de várias correlações de troca de calor encontradas na literatura. Um planejamento composto central permitiu avaliar simultaneamente o impacto da variação de seis parâmetros (diâmetro dos tubos no primeiro, segundo e terceiro passes, comprimento da caldeira, número de tubos nos passes 2 e 3 e porcentagem de excesso de ar na combustão) e obter um modelo polinomial empírico a fim de otimizar as eficiências energética e exergética de caldeira. Um teste de malha foi realizado com temperaturas no final da chama, câmara de combustão, segundo e terceiro passes e apontou para uma divisão suficiente de 200 volumes de controle em cada passe. O modelo foi também validado com dados obtidos experimentalmente, observando-se que a diferença entre os valores experimentais e teóricos das temperaturas do primeiro e segundo passes é inferior a 22,9ºC e menor que 43,9ºC no terceiro passe. Com os modelos polinomiais robustos obtidos com o planejamento de experimentos, foi possível observar os parâmetros mais significativos, tanto para as eficiências energética e exergética, assim como para a destruição de exergia. São eles, por ordem decrescente de importância: o excesso de ar; as duas contribuições associadas ao comprimento da caldeira; o número de tubos nos segundo e terceiro passes e; o diâmetro dos tubos no segundo e terceiro passe e o diâmetro do tubo no primeiro passe. Além disso, foi possível quantificar o trade-off entre o aumento da superfície de troca de calor e a manutenção do nível de turbulência do escoamento.

LInk da defesa:

Join Zoom Meeting: https://puc-rio.zoom.us/j/93257958274?pwd=UmN3dnRFV2NFenJDTUZCSXlXR2dHUT09

Meeting ID: 932 5795 8274

Passcode: 976191

The present work focuses on a numerical investigation of the flow through porous media with macropores, such as carbonates. The presence of vugs and fractures have a strong effect on the flow characteristics. The flow through the porous matrix is usually described by the Darcy equation and the flow through the macropores by Stokes equation. The coupling between these two distinct approaches brings great complexity to the modeling of such flows. In this work, we use Brinkman formulation that is able to describe the flow, both in the porous matrix and macropores with a single differential equation. We solved the set of differential equations using the finite element model and implemented the code in the FEniCS platform. We first solved the 1-D flow through parallel plates with one of the walls being a porous material. The goal was to compare the predictions obtained with the Brinkman formulation to that obtained by using the Beavers-Joseph boundary condition. Then, we solved a 2-D flow through a porous medium with macropores. The geometry of the pore structure was obtained from 2D slices of tomographic images of carbonates. The goal of this analysis was to evaluate an equivalent permeability as a function of macropores area and structures.