Removal of CO2 directly from anthropogenic sources (capture) and its disposal in geological formations can take place for medium-term time periods (storage), or it can be permanent (sequestration), with the CO2 eventually becoming dissolved in the aqueous phase. The latter is the main subject of this dissertation. Carbon dioxide sequestration covers a wide range of strategies and alternatives. The main objective of CO2 sequestration alternatives is secure disposal of carbon in large amounts and for a lengthy time scale (typically 1000 years). Injection of CO2 into subsurface formations is generally considered as the main option for CO2 sequestration. Geological sequestration through injection covers a broad variety of target formations: disposal in depleted oil and gas reservoirs, trapping in oil reservoirs, replacing CH4 in coal bed methane recovery processes, trapping in deep aquifers, and salt cavern placement are the major CCS alternatives in geologic formations. In this thesis, hydrogeologic interaction between the injectant (CO2) and the host fluid (saline water) during injection is the main subject of the project. Because of the density and viscosity contrast of displacing and displaced fluids, the pattern of saturation progression is complicated. A set of semi-analytical solutions is developed for quick estimation of the position of isosats (contours of saturation) during primary injection in homogenous cases with simple geometry. All of the mathematical solutions are developed based on two assumptions; incompressible fluids and rocks and vertical equilibrium (capillary-gravity condition) for geometries with large aspect ratio (L >> H). First, a series of analytical solutions for primary drainage for a set of linear relative permeability functions is developed. The first analytical solution is based on the assumption of locally linearized Leverett-J functions, and by using the method of characteristics, a formulation for the isosats’ geometry is obtained. A semi-analytical solution is then proposed for calculation of the position of isosats with linearized relative permeability functions and arbitrary capillary-saturation correlation. The analytical solution is extended to incorporate a specific form of nonlinearity of the relative permeability function. Nonlinear relative permeability functions are also incorporated in another semi-analytical solution, and the positions of the isosats for any arbitrary Leverett-J function and relative permeability functions are developed. Sequential gas-saline injection is also modeled in that chapter. For approximate verification of the analytical solutions, a FEM numerical model is developed and the results of the analytical solutions are compared with the numerical solutions. These new analytical solutions provide powerful tools for prediction of saturation distribution during injection in vertical and horizontal wells, as well as for carrying out stochastic assessments (Monte Carlo simulations) and parametric weight assessment. The domain of applications of the new solutions go far beyond the limited question of CO2 sequestration: they can be used for injection of any less viscous fluid into a reservoir, whether the fluid is lighter or denser than the host fluid (gas injection, water-alternating gas injection, water injection into viscous oil reservoirs, solvent injection).

The majority of commercially available ultrasonic transducers used to excite and measure wave propagation in structures can be coupled only to a free surface. While convenient, this method is likely to excite multiple structural modes, making data interpretation difficult. Furthermore, the many modes excited make predicting the structure’s response a computationally intensive task. Here the dynamic radial displacement induced by a transient radial point load is calculated at more than 230,000 points on the outer surface of a virgin steel pipe to simulate a typical experiment. The radial component of the displacement field is calculated by convolving the Green’s functions of the pipe with the transient load. These functions are calculated on personal computers (in a distributed arrangement) by employing modal summation. The mode shapes are obtained from a Semi-Analytical Finite Element formulation used in conjunction with a separation of variables. The results are presented in a four dimensional animation, providing easier interpretations and insight into how to best select observation points for the detection of defects. The accuracy of the calculated displacements is verified experimentally. Agreement is good when magnitude and phase corrections are incorporated from the frequency response curves of the transducers used.

This thesis presents a theoretical study of thermoelastic guided waves in cylinders in the context of Lord-Shulman generalized theory of thermoelasticity. Two different methods were formulated to study dispersion relations in infinite cylinders. One of them is a Semi Analytical Finite Element (SAFE) method and the other is an analytical method. In the SAFE method, the dispersion equation has been formulated as a generalized eigenvalue problem by treating radial displacement and temperature with a one dimensional finite element model through the thickness of the cylinder. In the analytical method, displacement potentials are introduced to obtain the dispersion relations of guided wave modes. This method is applicable to isotropic cylinders and has been developed primarily to cross check the SAFE formulation. Frequency spectra obtained by both methods for an isotropic cylinder have shown excellent agreement with each other. Since the SAFE method can be used for an anisotropic composite cylinder, guided wave modes for anisotropic and composite cylinders are presented. Transient analysis of ultrasonic guided waves generated by concentrated heating of the outer surface of an infinite anisotropic cylinder has also been studied. The SAFE method is employed to model the response of a cylinder due to a pulsed laser focused on its surface. Green’s functions were constructed numerically by superposition of guided wave modes in frequency and wave number domains. Time histories of the propagating modes are then calculated by applying an inverse Fourier transformation in the time domain. Transient radial displacements of longitudinal and flexural modes of a silicon nitride cylinder are presented. Propagation of thermoelastic waves in finite length circular cylinders have also been investigated. The SAFE method is used to simulate the guided wave modes in the cylinder. Frequency spectra obtained by the SAFE formulation, for a finite length transversely isotropic cylinder, are validated by comparing the numerical results with relevant publications. Frequency spectra for axisymmetric and asymmetric modes in a silicon nitride finite cylinder with both ends insulated and restrained by frictionless rigid walls are presented. The plain strain problem of circumferential guided waves is also studied and the results are validated for an isothermal case.

The majority of commercially available ultrasonic transducers used to excite and measure wave propagation in structures can be coupled only to a free surface. While convenient, this method is likely to excite multiple structural modes, making data interpretation difficult. Furthermore, the many modes excited make predicting the structure’s response a computationally intensive task. Here the dynamic radial displacement induced by a transient radial point load is calculated at more than 230,000 points on the outer surface of a virgin steel pipe to simulate a typical experiment. The radial component of the displacement field is calculated by convolving the Green’s functions of the pipe with the transient load. These functions are calculated on personal computers (in a distributed arrangement) by employing modal summation. The mode shapes are obtained from a Semi-Analytical Finite Element formulation used in conjunction with a separation of variables. The results are presented in a four dimensional animation, providing easier interpretations and insight into how to best select observation points for the detection of defects. The accuracy of the calculated displacements is verified experimentally. Agreement is good when magnitude and phase corrections are incorporated from the frequency response curves of the transducers used.

The advent of additive manufacturing has opened up new frontiers in developing metal structures that can have complex geometries, composite structures made of dissimilar metals, and metal structures with embedded sensing and actuation capabilities. These types of structures are possible with ultrasonic additive manufacturing (UAM); a novel manufacturing technology that combines additive manufacturing through the ultrasonic welding of thin metal foils with computer numerical control (CNC) milling. However, the process suffers from a critical limitation, i.e., a range of build heights within which bonding between a foil and the substrate cannot be originated. <br>This work has two research objectives, the first is a fundamental understanding of the complex dynamic interaction between the substrate and ultrasonic horn, or sonotrode. Specifically, it focuses on the effects that specific modes of vibration have on the dynamic response of the substrate. The second objective is to utilize the UAM process to create metal structures with an embedded sensor that can detect contact or impact. In addressing the first objective, a semi-analytical model was developed to determine the response to three forcing descriptions that approximate the interfacial friction between the foil and substrate induced by sonotrode compression and excitation. Several observations can be seen in the results: as the height increases the dominant modes of vibration change, the modes of vibration excited also change during a single weld cycle as the sonotrode travels across the length of the substrate, and finally the three forcing models do not have a significant impact on the substrate response trends with height and during the weld cycle. <br>In addressing the second objective, three prototypes were created by embedding a triboelectric nanogenerator (TENG) sensor within an AL3003 metal-matrix. TENGs utilize contact electrification between surfaces of dissimilar materials, typically polymers, combined with electrostatic induction to generate electrical energy from a mechanical excitation. The sensors demonstrate a discernible response over a 1-5 Hz frequency range. In addition, the sensors have a linear relationship between output voltage and a mechanically applied load, and have the ability to sense contact through both touch and due to an impacting object.

A survey of relevant literature on the topic of wave propagation and scattering in pipes is given first. This review is followed by a theoretical framework which is pertinent to wave propagation in homogeneous, isotropic, pipes. Emphasis is placed on approximate solutions stemming from a computer based, Semi-Analytical Finite Element (SAFE) formulation. A modal analysis of the dynamic response of homogeneous, isotropic pipes, when subjected to a transient ultrasonic excitation, demonstrates that dominant features, i.e., singularities in an unblemished pipe’s displacement Frequency Response Function (FRF) coincide with its cutoff frequencies. This behaviour is confirmed experimentally. A novel technique is developed to deduce such a pipe’s wall thickness and elastic properties from three cutoff frequencies. The resulting procedure is simulated numerically and verified experimentally. Agreement between the new ultrasonic procedure and traditional destructive tests is within experimental uncertainty. Then a hybrid-SAFE technique is used to simulate waves scattered by various open rectangular notches. The simulations show, for the first time, that singularities distinct from the unblemished pipe’s cutoff frequencies arise in a displacement FRF when an axisymmetric notch is introduced. They also suggest that the new singularities depend on the properties of the parent pipe and the finite element region but effects are local to a notch. It is demonstrated further that the difference between the frequency at which a singularity introduced by a notch occurs and the nearest corresponding unblemished pipe’s cutoff frequency is a function of the notch’s dimensions. By plotting contours of constant frequency differences, it is shown that it is usually possible to characterize the notch’s dimensions by using two modes. However, the frequency difference for a third mode may be also needed occasionally. The more general case of nonaxisymmetric notches is shown to be a straightforward extension of the axisymmetric case.

A survey of relevant literature on the topic of wave propagation and scattering in pipes is given first. This review is followed by a theoretical framework which is pertinent to wave propagation in homogeneous, isotropic, pipes. Emphasis is placed on approximate solutions stemming from a computer based, Semi-Analytical Finite Element (SAFE) formulation. A modal analysis of the dynamic response of homogeneous, isotropic pipes, when subjected to a transient ultrasonic excitation, demonstrates that dominant features, i.e., singularities in an unblemished pipe’s displacement Frequency Response Function (FRF) coincide with its cutoff frequencies. This behaviour is confirmed experimentally. A novel technique is developed to deduce such a pipe’s wall thickness and elastic properties from three cutoff frequencies. The resulting procedure is simulated numerically and verified experimentally. Agreement between the new ultrasonic procedure and traditional destructive tests is within experimental uncertainty. Then a hybrid-SAFE technique is used to simulate waves scattered by various open rectangular notches. The simulations show, for the first time, that singularities distinct from the unblemished pipe’s cutoff frequencies arise in a displacement FRF when an axisymmetric notch is introduced. They also suggest that the new singularities depend on the properties of the parent pipe and the finite element region but effects are local to a notch. It is demonstrated further that the difference between the frequency at which a singularity introduced by a notch occurs and the nearest corresponding unblemished pipe’s cutoff frequency is a function of the notch’s dimensions. By plotting contours of constant frequency differences, it is shown that it is usually possible to characterize the notch’s dimensions by using two modes. However, the frequency difference for a third mode may be also needed occasionally. The more general case of nonaxisymmetric notches is shown to be a straightforward extension of the axisymmetric case.

Les lanceurs spatiaux sont des structures complexes associant une multitude de composants. L’assemblage de ces éléments doit répondre à un niveau de performance élevé. Le collage structural demeure un bon candidat en raison des nombreux avantages qu’il présente. Cependant, cette technologie montre des inconvénients. En raison des changements brusques de géométrie et de propriétés matériaux, des concentrations de contraintes apparaissent aux extrémités du joint de colle. Ce phénomène appelé effets de bords est néfaste pour la tenue mécanique de l’assemblage collé. La présence des effets de bords exclut l’utilisation de critères en contrainte utilisés classiquement. Le dimensionnement d’assemblages collés requiert des outils fiables prenant en compte ces effets de bords. Dans cette étude, un modèle de ruine incrémentale, associant une approche en contrainte et en énergie, est utilisé. L’utilisation de cet outil dans un cadre industriel, impose de répondre aux besoins d’un Bureau d’Études, notamment en termes de coût de calculs. Afin de le diminuer, une implémentation semi-analytique, est tout d’abord développée. Puis, une seconde méthode d’implémentation, basée sur la méthode des Eléments Finis, permet une prévision plus précise de la ruine d’un assemblage. La pertinence de ces deux approches a été vérifiée pour plusieurs configurations de joints collés. Des campagnes d’essais, destinées à confronter les résultats expérimentaux aux prévisions numériques, ont été réalisées. Dans le cadre de ce travail, un montage de collage et d’essai pour assemblages tubulaires a en particulier été développé. L’objectif du pré-dimensionnement est d’identifier une zone d’intérêt dans l’ensemble du domaine d’étude. Aussi, une étude paramétrique peut être requise afin de réaliser cette tâche. Afin de réduire le coût de calcul, une méthode d’interpolation spatiale appelée Krigeage a été mise en œuvre, et permet la construction efficace d’une surface de réponse. / Space Launchers are complex structures composed of a large number of elements. The assembling of these components must show a high level of reliability. The use of adhesive bonding technology is an interesting solution since it presentsseveral assets compared to “classical” joint techniques (such as riveting, bolting and welding), mainly because it can help to construct lighter and less energy consuming systems However„ the implementation of adhesives also has somedrawbacks. Due to the strong variations of geometrical and material properties, stress concentrations appear at the extremities of the joint. This phenomenon; called edge effects; has a great influence on the failure of the bond. As a result, the simple use of a classical stress or energetic criteria is not appropriate to predict the fracture of such structures. Therefore, it is obvious that the design of bonded assemblies requires reliable tools to take the edge effects into account. In this work an incremental failure model, which combines the stress and energetic criteria, is used. In order to decrease the computational cost, a semi-analytical application of this model is proposed. This is intended to make the approach more interesting to be implemented in an industrial environment. The accuracy of the prediction of the failure load is enhanced by means of the Finite Element method. The reliability of both the semi-analytical and Finite Element approaches is verified by comparing the model predictions with experimental data issued from double-notched Arcan and tubular specimen geometries. The aim of the pre-design phase is to identify the critical area in the whole range of the application of the studied geometry. Therefore, the realization of a parametric study is required in order to build a response surface. In the present study, this has been achieved by means of spatial interpolation using the Kriging model.

Les structures mécaniques utilisées de nos jours ne cessent d’évoluer en utilisant des matériaux composites ou à gradient fonctionnel afin de répondre aux enjeux de résistance accrue, allégement de la structure et amélioration des performances. Ceux-ci nécessitent un contrôle adéquat de leur état de santé afin de s’assurer de l’intégrité de la structure. L’utilisation des ondes guidées ultrasonores fournit un moyen efficace et rapide d’inspection sur de longues distances. Néanmoins, ces ondes présentent certaines caractéristiques complexes qui rendent la tâche très difficile. L’utilisation d’outils d’analyse tels que les modèles numériques constitue un grand atout pour ce type d’application. Dans ce contexte, l’objectif de cette de thèse est le développement d’un outil de modélisation performant, permettant d’étudier la propagation des ondes guidées ultrasonores avec une grande précision et une faible consommation de ressources et de temps de calculs. De ce fait, l’intérêt est porté sur des méthodes numériques d’ordres élevés dont les propriétés de convergence sont beaucoup améliorées que les méthodes classiques. En particulier, la méthode semi-analytique éléments finis de Galerkin discontinue pour la détermination des courbes de dispersion des ondes guidées est développée. La méthode est applicable aux structures planes et cylindriques fabriquées de matériaux isotropes, anisotropes et hétérogènes (à gradient fonctionnel de propriétés). Une étude comparative sur l’analyse des performances de ces méthodes est effectuée. Celle-ci a démontré la capacité de la méthode à modéliser la propagation des ondes guidées ultrasonores dans des guides d’ondes à section arbitraire avec des performances prometteuses par rapport à la méthode des éléments finis classique. / The mechanical structures used today are constantly evolving using composite or functionally gradient materials to meet the challenges of increased strength, lightening the structure and improving performance. These require adequate control of their state of health to ensure the integrity of the structure. The use of Ultrasound Guided Waves (UGW) provides an efficient and fast way of inspection over long distances. Nevertheless, these waves have some complex features that make the task very difficult. The use of analysis tools such as numerical models is a great asset for this type of application. In this context, the objective of this thesis is the development of a powerful modeling tool, allowing to study the propagation of UGWs with a great precision, less computational time and consumption of resources. Accordingly, we are interested in higher order numerical methods whose convergence properties are much improved than the classical methods. In particular, a semi-analytical discontinuous Galerkin finite element method (SADG-FE) is developped for the determination of the dispersion properties of guided waves in arbitrary cross-section waveguides. The method is applicable to plates and cylindrical structures made of isotropic, anisotropic heterogeneous (functionally graded) materials. The performance analysis of these methods and their comparisons are performed with respect to the models based on the classical finite element method. The results demonstrated the ability of the proposed method to model the propagation of ultrasounic guided waves in arbitrary section waveguides with promising performance over the conventional finite element method.

In this thesis, I study numerical and semi-analytical models of self-gravitating protostellar discs, with the aim of furthering our understanding of the role of disc-self gravity in planet formation. At the time of writing, the ALMA era of observational astronomy is upon us. Therefore, I place my research into this context with synthetic images of both numerical and semi-analytical models. I begin with an examination into the apparent lack of convergence, with increasing resolution, of the fragmentation boundary in Smoothed Particle Hydrodynamics (SPH) simulations of a protostellar disc. I run a suite of SPH with different numerical implementations, and find that even very similar implementations can fundamentally change the final answer. I analyse a suite of SPH simulations that fragment to form gravitationally bound objects, with the motivation of informing future population synthesis model development. I find that fragment-fragment and fragment-disc interaction dominates the orbital evolution of the system even at very early times, and any attempt to produce a population of objects from the gravitational instability process must include these interactions. Before a disc fragments, it will go through a self-gravitating phase. If the disc cools globally on a timescale such that it is balanced by heating due to gravitational stresses, the disc will be in a state of quasi-equilibrium. So long as the disc mass is sufficiently low, and spirals are sufficiently tightly wound, then angular momentum transport can be described by the local approximation, for which there is an analytical description. Using this analytical description, I develop an existing 1D model into 3D, and examine a wide range of parameter space for which disc self-gravity produces significant non-axisymmetry. Using radiative transfer calculations coupled with synthetic observations, I determine that there is a very narrow range of parameter space in which a disc will have sufficiently large gravitational stresses so as to produce detectable spirals, but the stresses not be so large as to cause the disc to fragment. By developing a simple analytical prescription for dust, I show that this region of parameter space can be broadened considerably. However, it requires grains that are large enough to become trapped by pressure maxima in the disc, so I conclude that if self-gravitating spiral arms are detected in the continuum, it is likely that at least some grain growth has taken place.