Numerical simulations of atmospheric and oceanic flows are fundamentally limited by a lack of model resolution. This thesis describes the application of unstructured mesh finite element methods to geophysical fluid dynamics simulations. These methods permit the mesh resolution to be concentrated in regions of relatively increased dynamical importance. Dynamic mesh adaptivity can further be used to maintain an optimised mesh even as the flow develops. Hence unstructured dynamic mesh adaptive methods have the potential to enable efficient simulations of high Reynolds number flows in complex geometries. In this thesis, the thermally driven rotating annulus is used to test these numerical methods. This system is a classic laboratory scale analogue for large scale geophysical flows. The thermally driven rotating annulus has a long history of experimental and numerical research, and hence it is ideally suited for the validation of new numerical methods. For geophysical systems there is a leading order balance between the Coriolis and buoyancy accelerations and the pressure gradient acceleration: geostrophic and hydrostatic balance. It is essential that any numerical model for these systems is able to represent these balances accurately. In this thesis a balanced pressure decomposition method is described, whereby the pressure is decomposed into a ``balanced'' component associated with the Coriolis and buoyancy accelerations, and a ``residual'' component associated with other forcings and that enforces incompressibility. It is demonstrated that this method can be used to enable a more accurate representation of geostrophic and hydrostatic balance in finite element modelling. Furthermore, when applying dynamic mesh adaptivity, there is a further potential for imbalance injection by the mesh optimisation procedure. This issue is tested in the context of shallow-water ocean modelling. For the linearised system on an $f$-plane, and with a steady balance permitting numerical discretisation, an interpolant is formulated that guarantees that a steady and balanced state remains steady and in balance after interpolation onto an arbitrary target mesh. The application of unstructured dynamic mesh adaptive methods to the thermally driven rotating annulus is presented. Fixed structured mesh finite element simulations are conducted, and compared against a finite difference model and against experiment. Further dynamic mesh adaptive simulations are then conducted, and compared against the structured mesh simulations. These tests are used to identify weaknesses in the application of dynamic mesh adaptivity to geophysical systems. The simulations are extended to a more challenging system: the thermally driven rotating annulus at high Taylor number and with sloping base and lid topography. Analysis of the high Taylor number simulations reveals a direct energy transfer from the eddies to the mean flow, confirming the results of previous experimental work.
| Identifer | oai:union.ndltd.org:bl.uk/oai:ethos.bl.uk:547516 |
| Date | January 2011 |
| Creators | Maddison, James R. |
| Contributors | Read, Peter L. ; Marshall, David P. ; Piggott, Matthew D. |
| Publisher | University of Oxford |
| Source Sets | Ethos UK |
| Detected Language | English |
| Type | Electronic Thesis or Dissertation |
| Source | http://ora.ox.ac.uk/objects/uuid:4b95031b-4517-4aaf-9bb2-4d6d4a145499 |
Page generated in 0.0017 seconds