Liquid Storage Tanks (LSTs) are essential infrastructure systems that are used in various municipal and industrial settings. They play a critical role in storing and transporting liquids such as drinking water, oil, and gas that are used in daily life. Failure of these structures due to their poor seismic performance can have devastating consequences including the release of the stored liquid and damage to the surrounding area with potentially irreversible environmental impacts. In addition, the damage caused to the tank structure can be extensive, resulting in significant financial losses. Furthermore, the disruption of services provided by the tank such as water supply, oil and gas storage can be considerable. It is therefore crucial to study the seismic behaviour of these structures and ensure their safety and reliability to minimize the potential damages. The aim of this study is to investigate the behavior of LSTs in terms of the contained liquid, the tank’s structure, and its supports when subjected to seismic excitations. To obtain accurate results, different numerical techniques are applied in different phases of the study, including the Finite Volume Method (FVM), Finite Elements Method (FEM), Volume of Fluids (VoF) method, and Smoothed Particles Hydrodynamic (SPH) method. These techniques allow for a detailed analysis of the behavior of the tank and its supports during seismic excitation, providing a comprehensive understanding of the performance of LSTs during earthquakes. In order to examine the reliability and accuracy of the numerical model, the first part of the study includes the validation of the developed numerical model through comparison of the model and experimental results. The validated numerical model is then used to obtain the hydrodynamic pressures at different locations on the roof of tanks subjected to base excitations. The effect of liquid impact and hydrodynamic pressures on the roof of LSTs can be significant, however, limited studies have been completed on this issue. In the second part of the study, Artificial Intelligence (AI) techniques such as Genetic Programming are used to formulate these pressure values so that the maximum pressure can be obtained using the tank characteristics such as size and fill depth by the relationship obtained based on the AI approach. In the third part of the study, the supporting structures of LSTs subjected to base excitations are analyzed, and their shear forces are extracted and compared with the National Building Code of Canada (NBC 2020) in order to evaluate the reliability of the code and discuss possible improvements. The results of this study can be used to evaluate and make improvements in standards and guidelines for the seismic design of LSTs, which can help ensure the safety and reliability of these crucial infrastructures during seismic events.

The present master's thesis deals with numerical modeling of solid-state micrometrical-sized polymeric dye lasers, commonly denoted as microcavities. It is part of a large research initiative carried out in the optics group, at the MAP (Microelectronics and Applied Physics) department in KTH (Kungliga Tekniska Högskolan - Royal Technical School) and targeted towards the design and manufacturing of micro- and nano-scaled polymeric components for nano-photonics, primarily lasers. The finite element method (FEM) in frequency domain is used as a primary modeling tool through the simulation software COMSOL Multiphysics. Models for spontaneous emission, optical losses and gain are developed and demonstrated. A specic layout is studied: the double hexagonal microcavity. While it was expected to be a good candidate for a laser, the design shows unexpected properties making it useful for sensing applications. Finally, the transposition of models to time domain is initiated : a replacement solution for the lacking perfectly matched layer (PML) in Comsol is developed and demonstrated. Methods for modeling materials parameters in time domain are investigated, together with the possible use of a more suitable algorithm : finite dierences in time domain (FDTD) or Yee's scheme.

Optical Tomography is an imaging modality that enhances early diagnosis of disease through use of harmless Near-Infrared rays instead of conventional x-rays. The subsequent images are used to reconstruct the object. However Optical Tomography has not been effectively utilized due to the complicated photon scattering phenomenon and ill-posed nature of the corresponding image reconstruction scheme. The major method for reconstruction of the object is based on an iterative loop that constantly minimizes the difference between the predicted model of photon scattering with acquired images. Currently the most effective method of predicting the photon scattering pattern is the solution of the Radiative Transfer Equation (RTE) using the Finite Elements Method (FEM). However, the conventional FEM uses classical C0 interpolation functions, which have shortcomings in terms of continuity of the solution over the domain as well as proper representation of geometry. Hence higher discretization is necessary to maintain accuracy of gradient-based results which may significantly increase the computational cost in each iteration. This research implements the recently developed Isogeometric Approach (IGA) and particularly IGA-based FEM to address the aforementioned issues. The IGA-based FEM has the potential to enhance adaptivity and reduce the computational cost of discretization schemes. The research in this study applies the IGA method to solve the RTE with the diffusion approximation and studies its behavior in comparison to conventional FEM. The results show comparison of the IGA-based solution with analytical and conventional FEM solutions in terms of accuracy and efficiency. While both methods show high levels of accuracy in reference to the analytical solution, the IGA results clearly excel in accuracy. Furthermore, FE solutions tend to have shorter runtimes in low accuracy results. However, in higher accuracy solutions, where it matters the most, the IGA proves to be considerably faster. / Doctor of Philosophy / CT scans can save lives by allowing medical practitioners observe inside the patient's body without use of invasive surgery. However, they use high energy, potentially harmful x-rays to penetrate the organs. Due to limits of the mathematical algorithm used to reconstruct the 3D figure of the organs from the 2D x-ray images, many such images are required. Thus, a high level of x-ray exposure is necessary, which in periodic use can be harmful. Optical Tomography is a promising alternative which replaces x-rays with harmless Near-infrared (NIR) visible light. However, NIR photons have lower energy and tend to scatter before leaving the organs. Therefore, an additional algorithm is required to predict the distribution of light photons inside the body and their resulting 2D images. This is called the forward problem of Optical Tomography. Only then, like conventional CT scans, can another algorithm, called the inverse solution, reconstruct the 3D image by diminishing the difference between the predicted and registered images. Currently Optical Tomography cannot replace x-ray CT scans for most cases, due to shortcomings in the forward and inverse algorithms to handle real life usages. One obstacle stems from the fact that the forward problem must be solved numerous times for the inverse solution to reach the correct visualization. However, the current numerical method, Finite Element Method (FEM), has limitations in generating accurate solutions fast enough using economically viable computers. This limitation is mostly caused by the FEM's use of a simpler mathematical construct that requires more computations and is limited in accurately modelling the geometry and shape. This research implements the recently developed Isogeometric Analysis (IGA) and particularly IGA-based FEM to address this issue. The IGA-based FEM uses the same mathematical construct that is used to visualize the geometry for complicated applications such as some animations and computer games. They are also less complicated to apply due to much lower need for partitioning the domain. This study applies the IGA method to solve the forward problem of diffuse Optical Tomography and compare the accuracy and speed of IGA solution to the conventional FEM solution. The comparison reveals that while both methods can reach high accuracy, the IGA solutions are relatively more accurate. Also, while low accuracy FEM solutions have shorter runtimes, in solutions with required higher accuracy levels, the IGA proves to be considerably faster.

Σε αυτή τη Διπλωματική εργασία μελετάται η ταλαντωτική συμπεριφορά μιας κυκλικής πλάκας γραφενίου προκειμένου να διαπιστωθεί αν μπορεί να λειτουργήσει σαν νανομηχανική συσκευή ανίχνευσης μάζας. Για το λόγο αυτό δημιουργήθηκε και χρησιμοποιήθηκε ένα μοντέλο πεπερασμένων στοιχείων. Σε αυτό το μοντέλο οι αλληλεπιδράσεις μεταξύ των ατόμων και οι σχετικές κινήσεις αυτών προσομοιώθηκαν με κατάλληλα ελατήρια των οποίων η δυσκαμψία έχει υπολογιστεί χρησιμοποιώντας την ενέργεια δεσμού από τη θεωρία των μορίων. Ακόμη, τα άτομα του άνθρακα προσομοιώνονται με μάζες στη θέση των ατόμων άνθρακα. Το αναπτυχθέν μοντέλο αναπαριστά το δίσκο του γραφενίου σαν ορθότροπο υλικό. Προκείμενου να αποκτηθεί μια καλή εικόνα της συμπεριφοράς του κυκλικού γραφενίου γίνεται αρχικά διερεύνηση της ελεύθερης ταλάντωσης για διάφορες διαμέτρους αυτού. Επιπλέον, εξετάζεται η ταλαντωτική συμπεριφορά του γραφενίου με την προσκόλληση μάζας διαφόρων μεγεθών και σε διάφορες θέσεις πάνω στο γραφένιο κατά μήκος της ακτίνας του. Η επεξεργασία των μοντέλων έχει ως αποτέλεσμα τις φυσικές συχνότητες ταλάντωσης του γραφενίου και τις αντίστοιχες ιδιομορφές των φύλλων γραφενίου με ή χωρίς την προσκολλημένη μάζα παρέχοντας τα ταλαντωτικά τους χαρακτηριστικά. Επομένως, έχοντας δημιουργήσει μοντέλα για διάφορες διαμέτρους και μάζες μπορεί να εξεταστεί και να αποκτηθεί μια πλήρης εικόνα για την ενδεχόμενη ταλαντωτική συμπεριφορά του κυκλικού μονοστρωματικού δίσκου γραφενίου σαν αισθητήριο μάζας. / In this Diploma thesis is studied the vibrational behavior of a circular graphene sheet so as to investigate if this can operate as a nanomechanical system of mass sensing. For this purpose a finite elements’ model was created and being used. In this model the interatomic interactions and their relative positions are simulated by equivalent spring elements whose stiffness derived by corresponding potential energies provided by molecular theory and expressing the resistance in relative movements between carbon atoms under deformation. Moreover, the inertial effects of the system are simulated by using appropriate lumped masses in atomic positions. The model created represents the graphene plate as an orthotropic material. In order to obtain a full view of the behavior of circular graphene, firstly is been done investigation of its free vibration for different diameters of the plate. After that, is studied the Vibrational behavior with a mass lying on it. The solution gives the natural frequencies and corresponding mode shapes of vibration of individual graphene sheets and graphene-attached mass systems providing their vibrational characteristics. Subsequently, assuming the geometric characteristics of graphene as well as the weight and position of the attached mass as the global design parameters, a parametric study on mass sensing characteristics is presented in order to examine the potential behavior of a circular graphene monolayer as mass sensor.

Nesta pesquisa, estuda-se o fenômeno de localização de modos de vibração em estruturas moduladas quase periódicas de comportamento linear e não-linear. Em particular, contempla-se uma aplicação na Engenharia Civil, os painéis de placas periódicos com pequenas imperfeições fracamente acoplados entre si através de viga de grande rigidez, e principalmente submetidos à variação de forças de membrana introduzidas por meio da protensão, o que só pode ser levado em conta introduzindo a rigidez geométrica no modelo matemático. No caso de sistemas lineares, a presença de pequenas desordens nas características de rigidez ou massa de subsistemas fracamente acoplados pode causar confinamento espacial nas vibrações livres, conhecido como Localização de Modos, e pode inibir a propagação da resposta forçada. È o que tem sido mostrado na literatura técnica, em especial nas áreas de Engenharia Mecânica e Aeroespacial. Os efeitos de localização serão obtidos numa perspectiva modal. O programa de elementos finitos DYMPLATE implementado pelo autor para análise dinâmica não-linear de estruturas de placas, será utilizado para modelar estruturas periódicas (ordenadas) e quase periódicas (desordenadas). Os modelos são linearizados em torno de configurações deformadas de referência. O problema algébrico de autovalores é resolvido para obter as freqüências naturais e correspondentes modos de vibração. Estruturas planas constituídas por placas protendidas com módulos repetidos, pequenas imperfeições e diferentes condições de apoio e de carregamento, serão utilizadas na investigação numérica da influência de diversos fatores na Localização de Modos, em especial as forças de membrana. / In this research, the phenomenon of vibration modes localization in nearly periodic modular structures of linear and nonlinear behavior is studied. Of special interest is an application in Civil Engineering, lightly coupled periodic plate panels with small imperfections, mainly submitted to the variation of membrane forces introduced by prestress forces, which can be only considered by introducing geometric stiffness in the mathematical model. In the linear case, the presence of small disorders in the stiffness or mass characteristics of lightly coupled sub-systems can cause spatial confinement of free vibrations, known as Mode Localization, and can inhibit the propagation of the forced response. That is what has been shown in the literature, especially in the areas of Mechanical and Aerospace Engineering. The effects of localization are viewed from a modal perspective. DYMPLATE, a finite element software implemented by the author for nonlinear dynamic analysis of plates, will be utilized to model ordered and disordered plate periodic structures. The models are linearized about a deformed reference configuration. The algebraic eigenvalue problem is solved to obtain the natural frequencies and corresponding modes shapes. Plane structures constituted by prestressed plates with repetitive dynamic characteristics, small imperfections and different boundary conditions and loads, will be utilized in the numerical investigation of the influence of numerous factors in the Mode Localization.

Este trabalho consiste na análise do comportamento do método WAMG (Wavelet-Based Algebraic Multigrid), método numérico de resolução de sistemas de equações lineares desenvolvido no LMAG-Laboratório de Eletromagnetismo Aplicado, com relação a diversos suavizadores. O fato dos vetores que compõem os operadores matriciais Pronlongamento e Restrição do método WAMG serem ortonormais viabiliza uma série de análises teóricas e de dados experimentais, permitindo visualizar características não permitidas nos outros métodos Multigrid (MG), englobando o Multigrid Geométrico (GMG) e o Multigrid Algébrico (AMG). O método WAMG V-Cycle com Filtro Haar é testado em uma variedade de sistemas de equações lineares variando o suavizador, o coeficiente de relaxação nos suavizadores Damped Jacobi e Sobre Relaxação Sucessiva (SOR), e a configuração de pré e pós-suavização. Entre os suavizadores testados, estão os métodos iterativos estacionários Damped Jacobi, SOR, Esparsa Aproximada a Inversa tipo Diagonal (SPAI-0) e métodos propostos com a característica de suavização para-otimizada. A título de comparação, métodos iterativos não estacionários são testados também como suavizadores como Gradientes Conjugados, Gradientes Bi-Conjugados e ICCG. Os resultados dos testes são apresentados e comentados. / This work is comprised of WAMG (Wavelet-Based Algebraic Multigrid) method behavioral analysis based on variety of smoothers, numerical method based on linear equation systems resolution developed at LMAG (Applied Electromagnetism Laboratory). Based on the fact that the vectors represented by WAMG Prolongation and Restriction matrix operators are orthonormals allows the use of a variety of theoretical and practical analysis, and therefore gain visibility of characteristics not feasible through others Multigrid (MG) methods, such as Geometric Multigrid (GMG) and Algebraic Multigrid (AMG). WAMG V-Cycle method with Haar Filter is tested under a variety of linear equation systems, by varying smoothers, relaxation coefficient at Damped Jacobi and Successive Over Relaxation (SOR) smoothers, and pre and post smoothers configurations. The tested smoothers are stationary iterative methods such as Damped Jacobi, SOR, Diagonal type-Sparse Approximate Inverse (SPAI-0) and suggested ones with optimized smoothing characteristic. For comparison purposes, the Conjugate Gradients, Bi-Conjugate Gradient and ICCG non-stationary iterative methods are also tested as smoothers. The testing results are formally presented and commented.

Esta pesquisa apresenta uma metodologia, baseada em duas técnicas não invasivas para identificação e diagnóstico de pontos quentes em transformadores de potência durante plena operação ou em fase de projeto. A primeira técnica é baseada na medição de radiação infravermelha, emitida pelo equipamento em funcionamento e registrada por meio de análise termográfica. A segunda técnica é possível a partir do conhecimento prévio das características construtivas do transformador, fazendo uso do Método dos Elementos Finitos (MEF). A segunda técnica pode ser validada a partir das medições realizadas utilizando a primeira técnica. A formação de gases no interior dos transformadores de potência, devido à elevação da temperatura do óleo isolante em função dos pontos quentes, é discutida em detalhes com base nas normas técnicas estabelecidas pelo IEEE e IEC. As técnicas e procedimentos abordados ao longo dessa pesquisa foram obtidos a partir de um transformador com potência nominal de 120 MVA e relação de tensão 13,8/230 kV, projetado para integrar uma fazenda eólica ao Sistema Interligado Nacional (SIN). / This research presents a methodology based on two noninvasive techniques for identification and diagnostic of hot spots in power transformers during operation or project development. The first is based on measurements of infrared radiation from the equipment during operation and recording by thermography. The second technique is possible from the previous knowing of the constructive characteristics of the power transformer, by using the Finite Element Method (FEM). The second technique can be validated from measurements obtained using the first technique. The gas formation into the power transformers, because of the high temperatures in the insulating oil due to the hot spots, is discussed in details based on normative recommendations well established by the IEEE and IEC. All techniques and procedures to be approached in this research were obtained using a 120-MVA power transformer with voltage relationship of 13.8/230 kV that was projected to interconnect a wind farm to the Interconnected Brazilian System (Sistema Interligado Nacional - SIN).

L'objectif général de la thèse est la simulation numérique des mouvements forts du sol dûs aux séismes. Les déformations importantes du sol engendrent des comportements nonlinéaires dans les couches superficielles. L'apport principal de la thèse est la prise en compte de la nonlinéarité des milieux dans un contexte éléments finis Galerkin discontinus. Différentes lois de comportement sont implémentées et analysées. Le cas particulier du modèle élastoplastique de Masing-Prandtl-Ishlinskii-Iwan (MPII) est approfondi. Cette étude est divisée en deux parties. Une première qui vise à poser la structure du problème en présentant les équations et modèles utilisés pour décrire les mouvements du sol. Dans cette partie nous présentons également la méthode d’approximation spatiale Galerkin Discontinue ainsi que les différents schémas temporels que nous avons considérés. Une attention particulière est portée sur la complexité algorithmique du modèle nonlinéaire élastoplastique MPII en vue de réduire le temps de calcul des simulations. La deuxième partie est dédiée aux applications numériques. Ces applications sont réparties en trois catégories distinctes. 1) Nous nous intéressons toutd’abord à la configuration unidimensionnelle où une seule onde de cisaillement est propagée. Dans ce contexte, un flux numérique décentré est établi et des applications aux cas nonlinéaire élastique et nonlinéaire élastoplastique sont étudiées. Une solution analytique concernant le cas nonlinéaire élastique est proposée, ce qui permet de réaliser une étude numérique de convergence. 2) Le problème unidimensionnel étendu aux trois composantes du mouvement est étudié et utilisé comme un premier pas vers le 3D compte tenu du couplage entre les ondes de cisaillement et de compression. Nous nous intéressons ici à des signaux synthétiques et réels. L’application d’une méthode permettant de réduire significativement le temps de calcul du modèle élastoplastique est détaillée. 3) Une configuration tridimensionnelle est examinée. Après différentes applications de vérification en milieu linéaire, deux cas d’étude élastoplastique sont analysés. Une première sur un mode propre d’un cube puis une seconde sur un milieu plus réaliste composé d’un bassin hémisphérique à couches sédimentaires ayant un comportement élastoplastique / The general objective of this thesis is the numerical simulation of strong ground motions due to earthquakes. Significant deformations of the soil generate nonlinear behaviors in the superficial layers. The main contribution of this work is to take into account the nonlinearity of the media in a discontinuous Galerkin finite elements context. Different constitutive laws are implemented and analyzed. The particular case of theMasing-Prandtl-Ishlinskii-Iwan (MPII) elastoplastic model is looked at in-depth. This study is divided into two parts. A first one that aims at defining the framework of the problem by presenting the equations and models used to describe the soil motion. In this part we also present the Galerkin Discontinuous spatial approximation method as well as the different temporal schemes that we considered. Particular attention is paid to the algorithmic complexity of the nonlinear elastoplastic MPII model in order to reduce the computation time of simulations. The second part is dedicated to numerical applications. These applications are divided into three distinct categories. 1) We are first interested in the one-dimensional configuration where a single shear wave is propagated. In this context, an upwind numerical flux is established and applications to nonlinear elastic and nonlinear elastoplastic cases are studied. Ananalytical solution concerning the nonlinear elastic case is proposed, which makes it possible to carry out a numerical study of convergence. 2) The one-dimensional problem extended to the three components of the motion is studied and used as a first step towards 3D applications considering the coupling between the shear and compression waves. We are interested here in synthetic and real input signals. The application of a method that significantly reduces the calculation time of the elastoplastic model

O objetivo deste trabalho é o desenvolvimento de um pós-processador para um software baseado no método dos elementos finitos destinado a problemas de eletromagnetismo. Três métodos de determinação de torque de origem eletromagnética são abordados e comparados do ponto de vista da precisão e implementação, quais sejam: variação da energia magnética, integração do tensor de Maxwell e princípio do trabalho virtual (derivada do jacobiano local). Dois métodos são propostos para a determinação de reatâncias de eixo direto e em quadratura de máquinas síncronas. A metodologia é aplicada na análise de um motor síncrono de relutância de 15 hp, 380v, 4 pólos, cujos resultados são apresentados no capítulo IV. / The aim of this work is to develop a post-processor for use with an electromagnetism-oriented software based upon the Finite Elements Method. Three methods for determining electromagnetic torque are presented and compared in terms of precision and implementation: (i) Variation of Magnetic Energy, (ii) Maxwell Stress Method and (iii) Virtual-work Principle (Local Jacobian Derivative). Two methods for calculating direct-axis and quadrature-axis reactances of synchronous machines are proposed. The methodology is applied in the analysis of a 15 HP, 380 v, 4-pole reluctance synchronous motor and the results are presented in Chapter IV.

Os hidrocarbonetos são encontrados sob alta pressão em rochas porosas, denominadas rochas reservatório. A camada pré-sal apresenta rochas reservatório cobertas por uma camada impermeável de sal de grande espessura. Ocorrências de estruturas salinas são favoráveis para o aprisionamento dos hidrocarbonetos e aumentam a probabilidade de sucesso na prospecção de óleo e gás, pois são excelentes rochas capeadoras, de porosidade e permeabilidade praticamente nulas. As rochas salinas apresentam deformação lenta e contínua quando submetidas a tensões constantes, fenômeno conhecido como fluência. Durante a perfuração de poços através de espessas camadas de sal podem ocorrer problemas operacionais, como o aprisionamento da coluna de perfuração e o colapso do poço. Ainda, a fluência pode levar ao colapso os revestimentos de um poço de petróleo revestido. Estes contratempos geram grandes desafios e criam oportunidades de evolução na indústria do petróleo. Neste contexto, este trabalho consiste na análise numérica da escavação de poços em rochas salinas para exploração de petróleo, com o estudo do comportamento geomecânico do sal utilizando o programa de elementos finitos Abaqus®. O efeito da fluência do sal durante e após a perfuração e cimentação dos poços foi verificado utilizando diferentes pesos de fluido de perfuração, simulado como não penetrante. Análises dos deslocamentos, deformações e tensões na parede do poço e em sua vizinhança foram realizadas por meio de análises de deformação plana e análises axissimétricas. Após a instalação do revestimento e cimentação do poço que engloba tanto o processo de endurecimento do cimento, quanto a resposta de fluência da camada de sal, pôde-se analisar os deslocamentos, deformações e o comportamento das tensões na fronteira sal-cimento, evitando possíveis intervenções em poços que acarretam perdas econômicas. Sendo assim, este estudo auxilia no monitoramento e controle do fechamento de poços de petróleo em evaporito, após a escavação e cimentação do poço, evitando os diversos problemas decorrentes do comportamento de fluência do sal. / Hydrocarbons are found under high pressure in porous rocks, called reservoir rocks. The presalt layer shows reservoir rock covered by an impermeable salt thick layer. Occurrences of salt structures are favorable for trapping of hydrocarbons and increase the probability of success in oil and gas prospecting. They are excellent cap rocks with porosity and permeability practically nil. Salt rock present creep when subjected to continuous and constant stress. During well drilling through thick salt layers operational problems may occur like the imprisonment of the drill string and the collapse of the well. Also, creep may cause the collapse of the well casing. These setbacks create great challenges and opportunities for the evolution of the oil industry. In this context, this work proposes the numerical analysis of well excavation in salt rock for oil exploration with the study of the geomechanical behavior of salt using a finite element method (FEM) software Abaqus®. The effect of salt creep during and after drilling and cementing of wells was verified using different weights of the drilling fluid assumed to be non-penetrable with respect to the wellbore formation. Analysis of displacements, strains and stresses on the face of the wellbore and into the salt formation were performed with plane strain and axisymmetric techniques. After installation of the well casing, cementing is simulated encompassing cement hardening as well as salt creep. It was possible to analyze displacements, strains and the behavior of the stress interactions between the existing boundary cement-salt formation, avoiding possible unnecessary workover operations that cause economical losses. Thus, this study assists in the control and monitoring of closing oil wells in evaporite, after excavation and cementing the well, avoiding many problems stemming from the behavior creep of salt.