The first ocean general circulation models developed in the late sixties were based on finite differences schemes on structured grids. Many improvements in the fields of engineering have been achieved since three decades with the developments of new numerical methods based on unstructured meshes. Some components of the first models may now seem out of date and new second generation models are therefore under study, with the aim of taking advantage of the potential of modern numerical techniques such as finite elements. In particular, unstructured meshes are believed to be more efficient to resolve the large range of time and space scales present in the ocean.
Besides the classical continuous finite element or finite volume methods, another popular new trend in engineering applications is the Discontinuous Galerkin (DG) method, i.e. discontinuous finite elements presenting many interesting numerical properties in terms of dispersion and dissipation, errors convergence rates, advection schemes, mesh adaptation, etc. The method is especially efficient at high polynomial orders. The motivation for this PhD research is therefore to investigate the use of the high-order DG method for geophysical flow modeling.
A first part of the thesis is devoted to the mesh adaptation using the DG method. The inter-element jumps of the fields are used as error estimators. New mesh size fields or polynomial orders are then derived and local h- or p-adaptation is performed. The technique is applied to standard benchmarks and computations in more realistic domains as the Gulf of Mexico.
A second part deals with the use of the high order DG method with high-order representation of geometrical features. On one hand, a method is proposed to deal with complex representations of the coastlines. Computations are performed using high-order mappings around the Rattray island, located in the Great Barier Reef. Numerical results are then compared to in-situ measurements. On the other hand, a new method is proposed to deal with curved manifolds in order to represents oceanic or atmospheric flows on the sphere. The approach is based on the use of a local high-order non-orthogonal basis, and is equivalent to the use of vectorial shape and test functions to represent the vectorial conservation laws on the manifold's surface.
A method is finally proposed to analyze the dispersion and dissipation properties of any numerical scheme on any kind of grid, possibly unstructured. The DG method is then compared to other techniques as the mixed non-conforming linear elements, and the impact of unstructured meshes is studied.
| Identifer | oai:union.ndltd.org:BICfB/oai:ucl.ac.be:ETDUCL:BelnUcetd-10282008-142923 |
| Date | 14 November 2008 |
| Creators | Bernard, Paul-Emile |
| Publisher | Universite catholique de Louvain |
| Source Sets | Bibliothèque interuniversitaire de la Communauté française de Belgique |
| Language | English |
| Detected Language | English |
| Type | text |
| Format | application/pdf |
| Source | http://edoc.bib.ucl.ac.be:81/ETD-db/collection/available/BelnUcetd-10282008-142923/ |
| Rights | unrestricted, J'accepte que le texte de la thèse (ci-après l'oeuvre), sous réserve des parties couvertes par la confidentialité, soit publié dans le recueil électronique des thèses UCL. A cette fin, je donne licence à l'UCL : - le droit de fixer et de reproduire l'oeuvre sur support électronique : logiciel ETD/db - le droit de communiquer l'oeuvre au public Cette licence, gratuite et non exclusive, est valable pour toute la durée de la propriété littéraire et artistique, y compris ses éventuelles prolongations, et pour le monde entier. Je conserve tous les autres droits pour la reproduction et la communication de la thèse, ainsi que le droit de l'utiliser dans de futurs travaux. Je certifie avoir obtenu, conformément à la législation sur le droit d'auteur et aux exigences du droit à l'image, toutes les autorisations nécessaires à la reproduction dans ma thèse d'images, de textes, et/ou de toute oeuvre protégés par le droit d'auteur, et avoir obtenu les autorisations nécessaires à leur communication à des tiers. Au cas où un tiers est titulaire d'un droit de propriété intellectuelle sur tout ou partie de ma thèse, je certifie avoir obtenu son autorisation écrite pour l'exercice des droits mentionnés ci-dessus. |
Page generated in 0.0046 seconds