The processes studied by nuclear engineers generally include coupled physics phenomena (Thermal-Hydraulics, Neutronics, Material Mechanics, etc.) and modeling such multiphysics processes numerically can be computationally intensive. A way to reduce the computational burden is to use spatial meshes that are optimally suited for a specific solution; such meshes are obtained through a process known as Adaptive Mesh Refinement (AMR). AMR can be especially useful for modeling multiphysics phenomena by allowing each solution component to be computed on an independent mesh (Multimesh AMR). Using AMR on time dependent problems requires the spatial mesh to change in time as the solution changes in time. Current algorithms presented in the literature address this concern by adapting the spatial mesh at every time step, which can be inefficient. This Thesis proposes an algorithm for saving computational resources by using a spatially adapted mesh for multiple time steps, and only adapting the spatial mesh when the solution has changed significantly. This Thesis explores the mechanisms used to determine when and where to spatially adapt for time dependent, coupled physics problems. The algorithm is implemented using the Deal.ii fiinite element library [1, 2], in 2D and 3D, and is tested on a coupled neutronics and heat conduction problem in 2D. The algorithm is shown to perform better than a uniformly refined static mesh and, in some cases, a mesh that is spatially adapted at every time step.

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.