In this study we investigate the flow of a Boussinesq fluid contained in a rotating, differentially heated spherical shell. Previous work, on the spherical shell of Boussinesq fluid, differentially heated the shell by prescribing temperature on the inner boundary of the shell, setting the temperature deviation from the reference temperature to vary proportionally with -cos 20, from the equator to the pole. We change the model to include an energy balance equation at the earth's surface, which incorporates latitudinal solar radiation distribution and ice-albedo feedback mechanism with moving ice boundary. For the fluid velocity, on the inner boundary, two conditions are considered: stress-free and no-slip. However, the model under consideration contains only simple representations of a small number of climate variables and thus is not a climate model per se but rather a tool to aid in understanding how changes in these variables may affect our planet's climate. The solution of the model is followed as the differential heating is changed, using the pseudo arc-length continuation method, which is a reliable method that can successfully follow a solution curve even at a turning point. Our main result is in regards to hysteresis phenomenon that is associated with transition from one to multiple convective cells, in a dfferentially heated, co-rotating spherical shell. In particular, we find that hysteresis can be observed without transition from one to multiple convective cells. Another important observation is that the transition to multiple convective cells is significantly suppressed altogether, in the case of stress-free boundary conditions on the fluid velocity. Also, the results of this study will be related to our present-day climate. / UOIT

Modern geodynamic modeling is more complex than ever, and has been used to answer questions about Earth pertaining to the dynamics of the convecting mantle and core, layers humans have never directly interacted with. While the insights gleaned from these models cannot be argued, it is important to ensure calculations are understood and behaving correctly according to known math and physics. Here I perform several thermal 3-D spherical shell tests using the geodynamic code ASPECT, and compare the results against the legacy code CitcomS. I find that these two codes match to within 1.0% using a number of parameters. The application of geodynamic modeling is also traditionally to expand our understanding of Earth; however, even with a scarcity of data modern methods can provide insight into other planetary bodies. I use machine learning to show that coronae, circular features on the surface of the planet Venus, are not randomly distributed. I suggest the idea of coronae being fed by secondary mantle plumes in connected clusters. The entirety of the Venusian surface is poorly understood as well, with a large percentage being topographically smooth and much younger than the planet's hypothesized age. I use modeling to test the hypothesis of a large impact being responsible for a major resurfacing event in Venus's history, and find three distinct scenarios following impact: relatively little change, some localized change evolving into resurfacing through geologic time, or large-scale overturn and injection of heat deep into the Venusian mantle. / Doctor of Philosophy / Modern geodynamic modeling has been used to answer questions about Earth in wide-ranging fields. Despite technological improvements, it is important to ensure the calculations are understood and behaving correctly. Here I perform several tests using a code called ASPECT and compare the results against another code, CitcomS. I find that the two codes are in good agreement. Application of these techniques is also traditionally done for Earth, but modern methods can provide insight into other planets or moons as well. Coronae are circular features on the surface of Venus that are poorly understood. I use machine learning to show that these are not randomly distributed, and suggest a mechanism for the formation of clusters of coronae. The surface of Venus is also strange: it is both too flat and too young based on current ideas in planetary science. I use modeling to test whether a large impact could cause the details of Venus's surface we see today.

Sabe-se que o Método dos Elementos Finitos em sua forma convencional é uma ferramenta poderosa no cálculo estrutural moderno. Porém, se o problema apresenta singularidades, como os efeitos de borda tipicamente introduzidos pelos vínculos nas estruturas em casca, a análise pode exigir alto refinamento da malha. Procurando resolver mais eficientemente esse tipo de problema, especificamente em estruturas com simetria de revolução como os tubos cilíndricos e as cascas esféricas, apresentam-se neste trabalho alternativas não convencionais para o emprego do Método dos Elementos Finitos. Dadas as simetrias de forma e carregamento, a abordagem pode ser feita em campo unidimensional. São apresentadas as respostas analíticas, em termos de deslocamentos e esforços, para as estruturas citadas, partindo-se de suas equações diferenciais governantes. Em seguida, as formas fracas correspondentes são resolvidas pelo Método dos Elementos Finitos, incorporando-se alguns tipos de enriquecimento que aproveitam a estrutura deste método. Por fim, são comparados os resultados aproximados entre si e com relação aos analíticos, comprovando o grande potencial das alternativas sugeridas. / It is known that the Finite Element Method in its conventional form is a powerful technique in the modern structural calculus. However, if the problem has singularities, as boundary effects typically introduced by the support into shell structures, the analysis may demand high refinement of the mesh. In order to resolve this type of problem more efficiently, particularly in structures with symmetry of revolution as cylindrical tubes and spherical shells, this dissertation presents non-conventional alternatives for the employment of the Finite Element Method. Due to the form and load symmetry, the treatment can be made in one-dimensional system. The analytical responses, in terms of displacements and efforts, are recovered to the mentioned structures, from their governing differential equations. Next, corresponding weak forms are resolved by the Finite Element Method and some types of enrichment that utilize the structure of this method are incorporated. Finally, the approach results are compared among themselves and with analytical results, proving the great potential of the suggested alternatives

Sabe-se que o Método dos Elementos Finitos em sua forma convencional é uma ferramenta poderosa no cálculo estrutural moderno. Porém, se o problema apresenta singularidades, como os efeitos de borda tipicamente introduzidos pelos vínculos nas estruturas em casca, a análise pode exigir alto refinamento da malha. Procurando resolver mais eficientemente esse tipo de problema, especificamente em estruturas com simetria de revolução como os tubos cilíndricos e as cascas esféricas, apresentam-se neste trabalho alternativas não convencionais para o emprego do Método dos Elementos Finitos. Dadas as simetrias de forma e carregamento, a abordagem pode ser feita em campo unidimensional. São apresentadas as respostas analíticas, em termos de deslocamentos e esforços, para as estruturas citadas, partindo-se de suas equações diferenciais governantes. Em seguida, as formas fracas correspondentes são resolvidas pelo Método dos Elementos Finitos, incorporando-se alguns tipos de enriquecimento que aproveitam a estrutura deste método. Por fim, são comparados os resultados aproximados entre si e com relação aos analíticos, comprovando o grande potencial das alternativas sugeridas. / It is known that the Finite Element Method in its conventional form is a powerful technique in the modern structural calculus. However, if the problem has singularities, as boundary effects typically introduced by the support into shell structures, the analysis may demand high refinement of the mesh. In order to resolve this type of problem more efficiently, particularly in structures with symmetry of revolution as cylindrical tubes and spherical shells, this dissertation presents non-conventional alternatives for the employment of the Finite Element Method. Due to the form and load symmetry, the treatment can be made in one-dimensional system. The analytical responses, in terms of displacements and efforts, are recovered to the mentioned structures, from their governing differential equations. Next, corresponding weak forms are resolved by the Finite Element Method and some types of enrichment that utilize the structure of this method are incorporated. Finally, the approach results are compared among themselves and with analytical results, proving the great potential of the suggested alternatives

La Décomposition Modale Empirique (EMD) est un outil de traitement de signal piloté par les données et dédié aux signaux non-stationnaires issus ou non de systèmes linéaires. L'idée de base de l'EMD est l'interpolation des extrema par des splines pour extraire de composantes oscillantes appelées modes empiriques intrinsèques (IMFs) et un résidu. Dans cette thèse, un nouvel algorithme de l'EMD est introduit où au lieu d'une interpolation rigide, un lissage est utilisé pour la construction des enveloppes supérieures et inférieures du signal à décomposer. Ce nouvel algorithme est plus robuste au bruit que l'EMD conventionnelle et réduit le nombre d'IMFs "artificielles" (sur-décomposition). En combinant le nouvel algorithme et la méthode de séparation d'énergie (ESA) basée sur l'Opérateur d'Energie de Teager-Kaiser (OETK), un nouveau schéma de démodulation des signaux AM-FM multi-composante appelé EMD-ESA est introduit. Différentes versions de l'EMD-ESA sont analysées en terme de performance. Pour l'analyse Temps-Fréquence (TF), une nouvelle formulation de la carte TF de l'EMD-ESA appelée Transformation de Teager-Huang (THT) est présentée. Cette nouvelle Représentation TF (RTF) ne présentant pas de termes d'interférences est comparée aux RTF classiques telles que le spectrogramme, le scalogramme, la distribution de Wigner-Ville Distribution (WVD), la Pseudo-WVD et la réallocation de la Pseudo-WVD. En combinant la nouvelle formulation de la THT et la transformée de Hough, une nouvelle méthode de détection des signaux multi-composante à modulation linéaire de fréquence dans le plan TF est présentée. Cette méthode de détection est appelée transformation de Teager-Huang-Hough (THHT). Les résultats de la THHT sont comparés à ceux de la transformée WVD-Hough. Finalement, l'analyse TF par THT et par des RTF classiques (WVD, spectrogramme, etc.) de signaux réels de rétrodiffusion par des coques cylindriques de dimensions et de caractéristiques physiques différentes est présentée. Les résultats obtenus montrent l'apport de la THT comme un outil TF.