Made available in DSpace on 2015-03-04T18:51:18Z (GMT). No. of bitstreams: 1 Apresentacao.pdf: 121883 bytes, checksum: da58d642d3c4e264571ed24d208c16d7 (MD5) Previous issue date: 2006-03-15 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / In this dissertation, we study the stationary heat equation in a heterogeneous tridimensional plate, using a "dimension reduction" techinique called hierarchical modelling and we generate model the original problem in a two-dimensional domain. To estimate the error modelling, we develop an asymptotic expansion for the original problem solution and for the aproximate solution. Comparing both solutions with their own asymptotic expansions, we obtain an estimative of the error modelling. We perform some computational experiments, using the Residual Free Bubbles (RFB) Method and the Multiscale Finite Element Method for the diffusion problem and for the diffusion-reaction problem in a two-dimensional domain, with small parameters. Finally, we extend the numeric solutions found the original tridimensional problem. / Neste trabalho, estudamos a equação do calor estacionária em uma placa heterogênea tridimensional. Para a modelagem deste problema, utilizamos uma técnica de redução de dimensão conhecida por Modelagem Hierárquica. Desta forma, geramos um modelo para o problema original em um domínio bidimensional. Com o objetivo de estimar o erro de modelagem, desenvolvemos a expansão assintótica da solução do problema original e da solução aproximada. Comparando as soluções com suas respectivas expansões assintóticas, obtemos uma estimativa para o erro de modelagem. Realizamos alguns experimentos computacionais, desenvolvendo o método Residual Free Bubbles (RFB) e o método de Elementos Finitos Multiescala (MEFM) para o problema de difusão e para o problema de difusão-reação em um domínio bidimensional, com parâmetros pequenos. Com base nestes experimentos, encontramos algumas soluções numéricas para o problema da placa tridimensional.

Made available in DSpace on 2017-07-24T19:22:12Z (GMT). No. of bitstreams: 1 Miguel Munoz Perez.pdf: 2838722 bytes, checksum: 47bde0e445ce247501fdda653ceefbb2 (MD5) Previous issue date: 2010-02-25 / The aim of this study was to evaluate with the finite element method (MEF), the stress distribution under simulated functional loading in different cavity preparations using the pulp chamber of an endodontically-treated second mandible left molar to observe the stability and retention of the Endocrown. Computer tomography data of the tooth was obtained to generate a 3D volumetric model. In this 3D digital tooth, different cavity preparations without (Wo) and with (W) cusp coverage were modeled and edited, as well as a restoration and luting layer. All cavity preparations with cusp coverage had a mesio-occlusal-distal extension. To obtain the Poisson ratio and the elastic modulus of the restorative material used in the MEF, ceramic specimens were subjected to flexural tests. The following groups were obtained according to the cavity preparation: intact tooth (control) [IT]; occlusal-distal extension [WoOD], mesioocclusal- distal extension [WoMOD]; with buccal cusp coverage [WB], with lingual cusps coverage [WL]; with coverage of all cusps [WT]; and remnant of 2 mm above of the cement-enamel junction [Rm]. All groups were subjected to a 200 N loading in the Ansys Workbench v12 software and the stress distribution was analyzed by Rankine principal stresses criteria. "None" experimental group simulated stress distribution as the "IT". "WoOD" showed lowest stress levels at the internal angles. Stress levels were higher at internal angles for "W" and "WoMOD". No difference was observed between the other groups when compared to "Wo" and "W" with "MOD" extension. Within the limits of this in vitro study, it was concluded that the presence of a proximal wall was an important biomechanical factor for the remaining dental tissue and that the cusps coverage did not show any advantages for the stress distribution when compared to the more conservative cavity preparations. / O objetivo deste trabalho foi analisar pelo método de elementos finitos (MEF) a distribuição de tensões em diferentes preparos cavitários que utilizaram a câmara pulpar disponível de um dente molar desvitalizado para a estabilidade e retenção de restauração “Endocrown”, submetida à simulação de forças oclusais funcionais. A tomografia computadorizada de uma mandíbula foi digitalmente reconstruída resultando num modelo 3D. Deste modelo foi extraída apenas a parte correspondente ao dente 37. Neste dente digital 3D, diferentes preparos cavitários sem e com envolvimento das cúspides foram modelados e editados, bem como uma restauração e camada de cimento. Todos os modelos com envolvimento das cúspides tiveram uma extensão mesio-ocluso-distal (MOD). Para a obtenção do coeficiente de Poisson e módulo de elasticidade do material restaurador utilizados no MEF, uma cerâmica foi submetida a teste de resistência à flexão. Foram obtidos os seguintes grupos: DH, dente hígido (controle); dois grupos sem envolvimento das cúspides, SeOD (extensão ocluso-distal) e SeMOD (extensão mésio-ocluso-distal); três grupos com envolvimento das cúspides, CeV (envolvimento das cúspides vestibulares), CeL (envolvimento das cúspides linguais) e CeT (envolvimento de todas as cúspides); e RM, remanescente de 2 mm acima do limite esmalte–cimento (LEC). Todos os grupos foram submetidos à simulação de uma carga de 200 N, no software Ansys Workbench v12, avaliando a distribuição de tensões principais de Rankine geradas nos preparos. As simulações demonstraram que os melhores resultados foram obtidos no grupo SeOD. Os resultados dos outros grupos apresentaram leve diferença, sendo muito próximos os resultados obtidos entre preparos cavitários sem e com envolvimento de cúspides, de extensão MOD. Com base nos resultados obtidos, pôde-se concluir que a presença de uma parede proximal é importante fator para beneficiar a biomecânica do remanescente dental, e que o envolvimento cuspídeo não forneceu um melhor prognóstico da situação quando comparando a preparos cavitários mais conservadores.

[pt] Meios que apresentam descontinuidades como fraturas e falhas em um maciço rochoso ou reservatórios de petróleo impõem algumas dificuldades na simulação numérica pelo Método dos Elementos Finitos. Uma dessas dificuldades é a necessidade de geração de malhas muito refinadas, principalmente na região próxima à descontinuidade, até a obtenção de uma resposta confiável do modelo, o que pode consumir um tempo significativo. Ao mesmo tempo, a discretização da descontinuidade com elementos muito pequenos, quando comparados ao restante do modelo, pode conduzir a um aumento no tempo de simulação. Neste trabalho é apresentada a formulação de um elemento finito cortado por uma descontinuidade para aplicação em problemas com acoplamento fluido-mecânico em meios saturados com um único fluido. A inserção da descontinuidade no elemento é obtida pela adição de novos termos à função de interpolação, dispensando a sua discretização. Esses termos adicionais conseguem reproduzir a mudança no campo de deslocamento e poro-pressão no elemento devido à presença da descontinuidade. A resposta do elemento é verificada através da comparação com uma solução analítica unidimensional e com exemplos simples simulados em um programa comercial. / [en] Media that present discontinuities as fractures and faults in a rock mass or oil reservoirs impose some difficulties in numerical simulation using standard Finite Element Method. One of these difficulties is the need for very refined mesh generation, especially in the region near the discontinuities, to obtain a reliable answer of the model, which can consume significant time. At the same time, the discretization of the discontinuity with very small elements compared to the rest of the model may lead to an increase in simulation time. This work presents the formulation of an element that is crossed by a discontinuity for use in coupling fluid-mechanical problems in single fluid saturated mediums. The insertion of the discontinuity in the element is obtained by adding new terms in the interpolation function, which eliminates its discretization. These additional terms can reproduce the change in the displacement and pore pressure field in the element due to the presence of discontinuity. The response of the element is validated by comparing it with one dimensional analytical solution and simple examples simulated in a commercial program.

[pt] Soluções poliméricas diluídas podem ser usadas na indústria de Óleo e Gás para injeção em reservatórios de petróleo como agentes de controle de mobilidade em técnicas de recuperação avançada de óleo, com o objetivo de aumentar o fator de recuperação do campo e postergar a indesejada alta taxa de produção de água. No entanto, o fenômeno macroscópico esperado, aumento da eficiência de varrido, é influenciado diretamente por fatores microscópicos que ocorrem na escala de poro. Devido às suas grandes cadeias moleculares, ao avançar através do meio poroso, soluções poliméricas de alto peso molecular dificultam o avanço da água alterando a razão de mobilidades água/óleo. Porém, o estudo de escoamento de soluções poliméricas ainda é um desafio devido ao comportamento não Newtoniano da relação vazão x diferença de pressão e à sensibilidade a diversos fatores como: tipo de polímero utilizado, concentração, razão de viscosidades, temperatura, salinidade e taxa de cisalhamento. O presente trabalho analisa a implementação de dois modelos constitutivos, Oldroyd-B e FENE-CR, aplicados ao método de solução de elementos finitos, EF, utilizando a formulação dos resíduos ponderados de Galerkin para modelar o escoamento de soluções poliméricas diluídas através de capilares com constrição e entender os fenômenos viscosos e elásticos envolvidos. Os resultados obtidos fornecem uma descrição mais detalhada da aplicação dos modelos testados e do escoamento de soluções poliméricas diluídas em gargantas de poros. / [en] Diluted polymeric solutions can be used in Oil & Gas industry for reservoir injection as mobility control agents in enhanced oil recovery technics, which goal is to increase the oil field recovery factor and postpone the high water production rates. However, the macroscopic expected phenomena, increase of the areal sweep, is directly influenced by microscopic factors that occurs at the pore scale. Due to their huge molecular chains, by advancing through the porous media, high molecular weight polymer solutions hinder the water flow changing the water/oil mobility ratio. Nevertheless, the comprehension of polymeric solutions flow is still a great challenge due to the Non-Newtonian behavior of the flow rate x pressure gradient relation and the sensibility to several variables such as: type of polymer used, concentration, viscosity ratio, temperature, salinity and shear rate. The present work analyzes the implementation of two constitutive models, Oldroyd-B and FENE-CR, coupled with the Finite Element method using the Galerkin weighted residual formulation to model the flow of diluted polymeric solutions through constricted capillaries and understand the viscous and elastic phenomena involved. The obtained results provide a more detailed portrait of the constitutive models and the flow of diluted polymeric solutions through pore throats.

[pt] Neste trabalho apresentam-se três tipos de técnicas de análise do colapso estrutural através do método dos elementos finitos: análise linearizada da carga crítica, análise incremental da carga crítica e análise não linear completa. Na análise linearizada da carga crítica formulou-se um problema de autovalor empregando matrizes de rigidez baseadas na configuração indeformada da estrutura e materiais com comportamento linear elástico. No caso da análise incremental da carga crítica, o problema de autovalor foi formulado empregando matrizes de rigidez incrementais para levar em consideração os grandes deslocamentos e propriedades não lineares do material. Finalmente, na análise não linear completa a configuração deformada da estrutura e propriedades não lineares do material são atualizadas durante todo o processo incremental-iterativo até atingir a carga crítica. Desenvolveu-se uma implementação computacional para estudar as três técnicas de análise em estruturas planas como vigas, colunas, pórticos e arcos, empregando elementos isoparamétricos bidimensionais para estado plano de tensões. A configuração deformada da estrutura, devido aos grandes deslocamentos e rotações dos elementos, foi considerada através de uma formulação Lagrangeana Total, enquanto o comportamento inelástico do material foi modelado empregando um modelo elastoplástico de Von Mises (J2) com encruamento isotrópico. Nos exemplos apresentados mostrou-se a influência da não linearidade geométrica e física na estimativa de cargas críticas e no comportamento pós-crítico, podendo ocorrer bifurcações ao longo da trajetória de equilíbrio fundamental definida no espaço carga-deslocamentos. / [en] This work presents three kinds of techniques for collapse analysis using the finite element method: linear buckling analysis, nonlinear buckling analysis and full nonlinear analysis. The linear buckling analysis requires the definition of an eigenvalue problem using a stiffness matrix formulation based on the initial configuration of the structure and under the assumption of a linear elastic material behavior. In the case of nonlinear buckling analysis, the eigenvalue problem was formulated employing an incremental stiffness matrix in order to consider the effects of large displacements and nonlinear material properties in the critical load estimation. Finally, the full nonlinear analysis takes into account the deformed configuration and the nonlinear material properties of the structure, updating both of them through all the incremental-iterative process up to reaching the critical load. A Finite Element computational program, using plane stress isoperimetric bidimensional elements, was developed to study the three analysis techniques applied to plane structures such as beams, columns, frames and arches. The deformed configuration of the structure, due to large displacements and rotations, was considered through the Total Lagrangian formulation, whereas the inelastic material behavior was modeled using the Von Mises plasticity model with isotropic hardening. The examples presented in this article show the influence of geometric and material nonlinearity in the critical load estimation and the postcritical behavior, being this the reason for the potential occurrence of bifurcation points over the fundamental equilibrium path defined in the load-displacement space.

[pt] A Mecânica da Fratura é uma área da mecânica clássica que estuda os processos que podem resultar no surgimento e propagação de fraturas e trincas nos materiais. Seus conceitos podem ser empregados em diversas áreas, como, por exemplo, no desenvolvimento de novos materiais, na área industrial (auxiliando na definição dos parâmetros de criticidade de defeitos) e em projetos estruturais, onde a presença de descontinuidades pode resultar na diminuição da resistência do material, levando-o ao colapso estrutural. As fraturas podem surgir como defeitos básicos nos materiais constituintes dos elementos, podem ser induzidas nos processos de construção ou, ainda, podem ser desenvolvidas durante a vida útil da estrutura. O estudo do comportamento de uma fratura é fundamental para a verificação coerente de tensões e deformações nos projetos estruturais. Alguns autores apresentaram diferentes metodologias para o cálculo de parâmetros importantes associados à propagação de fraturas. Tais parâmetros podem variar de acordo com o comportamento do material, características geométricas, carregamentos, condições de contorno e configuração das trincas. Essas análises também podem ser validadas através de metodologias numéricas. O presente projeto tem por objetivo avaliar o Fator de Intensidade de Tensões a partir de análises numéricas de três casos de fratura com base nos conceitos da Mecânica da Fratura Linear Elástica (MFLE). Para tal, aplicaram-se duas técnicas de modelagem em Elementos Finitos: elementos quarter-point 2D e 3D e o Método dos Elementos Finitos Estendidos (XFEM) 3D. Por fim, os resultados obtidos são comparados com os resultados publicados na literatura. Apesar das vantagens associadas ao uso do XFEM para modelagem de fraturas, o cálculo pela integral de domínio para esta técnica apresenta oscilações nos valores fornecidos para as diferentes solicitações de contornos. A técnica de modelagem com Elementos Quarter-Points 2D e 3D apresenta resultados mais estáveis e próximos das soluções analíticas. / [en] Fracture Mechanics is an area of classical mechanics that studies processes that can result in the creation and propagation of fractures and cracks in materials. Its concepts can be utilized in many areas, such as, for example, the development of new materials, in the industrial area (to assist in the definition of defect criticality parameters) and in structural projects, where the presence of discontinuities can result in decreased material resistance, leading to its structural collapse. The fractures can emerge as basic defects in materials that constitute the structural elements, can be inducted in construction processes or could even be developed during the lifespan of the structure. Studying the behavior of a fracture is fundamental to verifying coherently stress and deformations on structural projects. Some authors presented different methodologies to calculate important parameters associated to the propagation of cracks. These parameters could vary according to the material behavior, geometrical characteristics, loads, boundary conditions and cracking patterns. Those analyses can also be validated through numerical methodologies. This present project aims at evaluating the Stress Intensity Factor from numerical analyses of three cracks cases based on the concepts of Linear Elastic Fracture Mechanics (LEFM). Two techniques of Finite Element modeling were considered: quarter-point elements 2D and 3D and Extended Finite Element Method (XFEM) 3D. In the end, the results obtained are compared with results already published in the available literature. Despite the advantages associated with the use of XFEM for fracture modeling, the calculation of the stress intensity factor by the domain integral for this technique presents oscillations in the values provided for the different contour requests. The modeling technique with 2D and 3D quarter-points elements presents results that are more stable and closer to the analytical solutions.

[pt] Este trabalho analisa excitações verticais dinâmicas de máquinas para fundações com embutimento. Combinações de fundações superficiais e profundas, para simulação de bloco de fundações, são analisadas com uso de diferentes metodologias. Analisam-se soluções de Novak e Beredugo em combinação com Novak e El-Sharnouby, com uso de porcentagens de recalque adicional retiradas do trabalho de Poulos e Davis. Também são analisadas soluções simplificadas de fundações superficiais com embutimento de Wolf em combinação com fundações profundas de Novak e El-sharnouby, com uso de porcentagens de recalque adicional dinâmicas,obtidas com elementos finitos neste trabalho. São comparados os valores de deslocamentos de regime permanente com os obtidos por elementos finitos tridimensionais com uso do programa Plaxis 3D, comprovando a eficácia dos valores de taxa de recalque adimensional dinâmico e uso de modelo simplificado de fundações superficiais de Wolf. / [en] This paper analyzes dynamic vertical loadings due to machine foundations on embedded foundations. Combinations of shallow and deep foundations, to pile raft foundations, are analyzed using diferent approaches. First, solutions of Novak and Beredugo combined with Novak and El- Sharnouby, using additional settlement ratios from Poulos and Davis. Then it is compared to simplified embedded shallow foundations solutions by Wolf combined with Novak and El-Sharnouby, using dynamic additional settlement ratio obtained by finite elements on this paper. Settlements obtained on the steady state regime are compared to the values obtained by tridimensional finite elements using the program Plaxis 3D, proving the efficiency of dynamic additional settlement ratio and the simplified model of shallow foundations by Wolf.

[pt] O crescente uso de antenas para a transmissão de dados tem demandado a instalação de um número crescente de torres, com destaque para as torres estaiadas. Essas estruturas são compostas por um mastro central e diversos níveis de estais. As torres são estruturas esbeltas e apresentam um comportamento altamente não linear sob cargas estáticas e dinâmicas, devido à interação mastro-cabos e altura elevada da torre. O objetivo deste trabalho é estudar a influência da não linearidade no comportamento estático e dinâmico de uma torre estaiada. Para isto, desenvolvem-se modelos de torre estaiada em elementos finitos, com dois níveis de cabos pré-tensionados. Consideram-se dois tipos de geometria para a disposição dos cabos: em leque e em paralelo. Em particular, investiga-se o efeito da ruptura de um ou mais cabos e do comportamento elastoplástico dos materiais que compõem a torre estaiada na estabilidade e vibrações não lineares da estrutura. Inicialmente, avalia-se o comportamento linear da estrutura, obtendo-se as frequências naturais e cargas críticas. A seguir, obtêm-se os caminhos não lineares de equilíbrio de cada modelo, aplicando o conceito de imperfeição modal. Finalmente, realiza-se a análise da vibração livres e forçada amortecida. Para isto utilizam-se as respostas no tempo e plano de fase, transformadas de Fourier e espectrográficos. Em todos os casos analisa-se em detalhe o efeito da perda de cabos e da plastificação dos elementos na capacidade de carga, estabilidade e deslocamentos da estrutura. Os resultados demonstram que a perda de cabos e plastificação dos elementos têm grande influência nas frequências naturais e cargas críticas. Finalmente, nessas condições observa-se um comportamento não linear com perda de capacidade de carga e aumento significativo de deslocamentos, velocidades e acelerações, evidenciando a necessidade de investigação do comportamento não linear desse tipo de estrutura. / [en] The growing use of antennas for data transmission calls for the installation of growing number of high towers, with the cable-stayed towers standing out. These structures are composed of a central mast and several levels of stays. These towers are slender structures and exhibit a highly non-linear behavior under static and dynamic loads, due to the mast-cable interaction and the tower height. The objective of this work is to study the influence of the nonlinearity on the static and dynamic behavior of a cable-stayed tower. For this, finite element models of cable-stayed towers with two levels of pre-tensioned cables are developed. Two cable disposition geometries are considered: fan and parallel. In particular, the effect of the rupture of one or more cables and the elastoplastic behavior of the cable-stayed tower materials on the stability and non-linear vibrations of the structure is investigated. Initially, the linear behavior of the structure is evaluated through the natural frequencies and critical loads. Next, the non-linear equilibrium paths of each model are obtained, applying the concept of modal imperfection. Finally, free and forced damped vibration analysis is performed. For this purpose, responses in time and phase plane, Fourier transforms and spectrographs are used. In all cases, the effect of the loss of cables and the plastification of the elements on the load capacity, stability and displacements of the structure is analyzed in detail. The results demonstrate that the loss of cables and plastification of elements have a great influence on natural frequencies and critical loads. Under such conditions, non-linear behavior with loss of load capacity and significant increase in displacements, velocities and accelerations stand out, and highlight the need for investigation of the nonlinear behavior of this type of structure.

[pt] A tecnologia de poços inteligentes já possui muitos anos de experiência de campo. Inúmeras publicações tem descrito como o controle de fluxo remoto e os sistemas de monitoração podem diminuir o número de intervenções, o número de poços e aumentar a eficiência do gerenciamento de reservatórios. Apesar da maturidade dos equipamentos de completação o conceito de poço inteligente integrado como um elemento chave do Digital Oil Field ainda não está completmente desenvolvido. Sistemas permanentes de monitoração nesse contexto tem um papel fundamental como fonte da informação a respeito do sistema de produção real visando calibração de modelos e minimização de incerteza. Entretanto, cada sensor adicional representa aumento de complexidade e de risco operacional. Um entendimento fundamentado do que realmente é necessário, dos tipos de sensores aplicáveis e quais técnicas de análises estão disponíveis para extrair as informações necessárias são pontos chave para o sucesso do projeto de um poço inteligente. Este trabalho propõe uma nova forma de tratar os dados em tempo real de poços inteligentes através da centralização do pré-processamento dos dados. Um modelo poço inteligente numérico para temperatura em regime transiente foi desenvolvido, testado e validado com a intenção de gerar dados sintéticos. A aplicação foi escolhida sem perda de generalidade como um exemplo representativo para validação dos algorítmos de limpeza e extração de características desenvolvidos. Os resultados mostraram aumento da eficiência quando comparados com o estado da arte e um potencial para capturar a influência mútua entre os processos de produção. / [en] Intelligent Well (IW) technology has built-up several years production experience. Numerous publications have described how remote flow control and monitoring capabilities can lead to fewer interventions, a reduced well count and improved reservoir management. Despite the maturity of IW equipment, the concept of the integrated IW as a key element in the Digital Oil Field still not fully developed. Permanent monitoring systems in this framework play an important role as source of the necessary information about actual production system aiming model calibration and uncertainty minimization. However, each extra permanently installed sensor increases the well s installation complexity and operational risk. A well-founded understanding of what data is actually needed and what analysis techniques are available to extract the required information are key factors for the success of the IW project. This work proposes a new framework to real-time data analysis through centralizing pre-processing. A numeric IW transient temperature model is developed, tested and validated to generate synthetic data. It was chosen without loss off generality as a representative application to test and validate the cleansing and feature extraction algorithms developed. The results achieved are compared with the state of the art ones showing advantages regarding efficiency and potential to capture mutual influence among processes.

[pt] Ao longo dos anos, a produção de energia tem dependido de recursos não sustentáveis, como os combustíveis fosséis. No entanto, com o aquecimento global e a crise energética urge-se investir em recursos de energia renovável, como o hidrogênio. O gás deve ser armazenado em um ambiente seguro para evitar vazamentos. Portanto, este trabalho foca no armazenamento de hidrogênio em cavernas de sal, uma vez que essas rochas possuem propriedades relevantes, como a baixa permeabilidade. Um fluxo de trabalho para análise de integridade de cavernas desde a construção até a operação é proposto, implementado e aplicado para o estudo de casos sintéticos e reais. O armazenamento de hidrogênio provoca variações de temperatura e pressões dentro da caverna. A termodinâmica do gás segue uma solução diabática, atualizando a pressão e a temperatura do gás a cada instante para representar cenários de campo. A formulação termomecânica é implementada no simulador GeMA, que acopla diferentes físicas. Casos sintéticos consideram modelos homogêneos e diferentes geometrias de caverna. Os resultados demonstraram a importância dos efeitos térmicos, pois as amplitudes térmicas podem comprometer a integridade da rocha, por exemplo, induzindo tensões de tração e afetando a permeabilidade. Um estudo hidráulico demonstrou risco mínimo de migração de gás para o exterior. Por último, dois casos reais foram investigados, litologia heterogênea e uma caverna irregular baseada em dados de sonar. Os resultados evidenciaram alguns desafios na operação de cavernas. / [en] Over the years, energy has been highly dependent on non-sustainable resources. However, global warming and the energy crisis urge for investments in renewable energy resources, such as hydrogen. The gas must be stored in a secure medium to avoid migration to the external environment. Thus, this work focuses on hydrogen storage in salt caverns, as these rocks present relevant properties for a storage site, such as low permeability. A workflow for cavern analysis from construction to operation is proposed, implemented, and applied to synthetic and actual field cases. Hydrogen storage provokes variations in temperature and pressure inside the cavern. The gas thermodynamics follows a diabatic solution, which updates gas pressure and temperature to represent the field conditions. The thermomechanical formulation is implemented into an in-house framework GeMA, which couples different physics. Synthetic case studies include homogeneous deposits and different cavern geometries. The results demonstrate the importance of thermal effects, as temperature amplitudes may compromise rock integrity, inducing tensile stresses and affecting its permeability. A hydraulic study demonstrated minimal hydrogen migration risk. Finally, two real field conditions were investigated, considering heterogeneous salt stratifications and a sonar-based cavern geometry. The results highlight some integrity challenges to be faced during cavern operation.