Numerical simulations of flow–topography interaction using unstructured grids
Author(s)
McVicar, Alistair J.
Type
Thesis
Abstract
Eddies and jets are important components of global ocean momentum and heat budgets
but are typically unresolved in low resolution global climate models. Herein, they
are evaluated with an idealised model set–up that incorporates barotropic flow, past a
cylinder on a β –plane. The flow dynamics are a function of two non–dimensional numbers:
the Reynolds number and the [Symbol appears here. To view, please open pdf attachment] –parameter.
The model used, Fluidity–ICOM, utilises unstructured meshes and a new stable mixed
discontinuous/continuous finite element pair (P1DGP2). Unstructured meshes decrease
the computational cost; the simulations using a non–uniform unstructured mesh had approximately
40% fewer nodes and ran at twice the speed of a uniform structured mesh for
a comparable drag coefficient (Cd). The validation of Fluidity–ICOM was performed for
a range of Reynolds numbers (0:0 < Re [Mathematical symbol appears here. To view, please open pdf attachment] 3 x 10[to the power of six]) and the percentage difference between
published and Fluidity–ICOM values of Cd was found to be less than 10% for the regimes
where the dynamics are essentially two–dimensional. The validation highlighted two important
considerations: the position of the lateral domain boundary and the boundary
mesh resolution.
The wake structure for a moderate Reynolds number (1000) and [Symbol appears here. To view, please open pdf attachment] –parameter (75)
changed significantly between coarse and fine boundary resolutions. The former was
comprised of a double jet structure and the latter a single jet in the lee of the cylinder.
This study demonstrated that resolving the frictional boundary layer dynamics is crucially
important, as they substantially impact on the downstream flow. Evaluation of the single
jet structure for a large parameter space [Mathematical formula appears here. To view, please open pdf attachment] revealed the presence
of interfacial Rossby waves with both eastward and westward propagation with respect
to the mean flow. The Rossby wave occurred due to the presence of a strong staircase
gradient in absolute vorticity. As the Reynolds number increased for a fixed [Symbol appears here. To view, please open pdf attachment] –parameter,
the presence of a stronger shear resulted in a faster phase speed of the Rossby wave and
a stronger mean–flow. This parameter–space also showed a large dependence on drag to
the [Symbol appears here. To view, please open pdf attachment] –parameter. Overall, this study has implications for the Gulf Stream separation
and for understanding the interaction of the Antarctic Circumpolar Current (ACC) with
topography.
but are typically unresolved in low resolution global climate models. Herein, they
are evaluated with an idealised model set–up that incorporates barotropic flow, past a
cylinder on a β –plane. The flow dynamics are a function of two non–dimensional numbers:
the Reynolds number and the [Symbol appears here. To view, please open pdf attachment] –parameter.
The model used, Fluidity–ICOM, utilises unstructured meshes and a new stable mixed
discontinuous/continuous finite element pair (P1DGP2). Unstructured meshes decrease
the computational cost; the simulations using a non–uniform unstructured mesh had approximately
40% fewer nodes and ran at twice the speed of a uniform structured mesh for
a comparable drag coefficient (Cd). The validation of Fluidity–ICOM was performed for
a range of Reynolds numbers (0:0 < Re [Mathematical symbol appears here. To view, please open pdf attachment] 3 x 10[to the power of six]) and the percentage difference between
published and Fluidity–ICOM values of Cd was found to be less than 10% for the regimes
where the dynamics are essentially two–dimensional. The validation highlighted two important
considerations: the position of the lateral domain boundary and the boundary
mesh resolution.
The wake structure for a moderate Reynolds number (1000) and [Symbol appears here. To view, please open pdf attachment] –parameter (75)
changed significantly between coarse and fine boundary resolutions. The former was
comprised of a double jet structure and the latter a single jet in the lee of the cylinder.
This study demonstrated that resolving the frictional boundary layer dynamics is crucially
important, as they substantially impact on the downstream flow. Evaluation of the single
jet structure for a large parameter space [Mathematical formula appears here. To view, please open pdf attachment] revealed the presence
of interfacial Rossby waves with both eastward and westward propagation with respect
to the mean flow. The Rossby wave occurred due to the presence of a strong staircase
gradient in absolute vorticity. As the Reynolds number increased for a fixed [Symbol appears here. To view, please open pdf attachment] –parameter,
the presence of a stronger shear resulted in a faster phase speed of the Rossby wave and
a stronger mean–flow. This parameter–space also showed a large dependence on drag to
the [Symbol appears here. To view, please open pdf attachment] –parameter. Overall, this study has implications for the Gulf Stream separation
and for understanding the interaction of the Antarctic Circumpolar Current (ACC) with
topography.
Date Issued
2012-03
Date Awarded
2012-08
Copyright Statement
Attribution NoDerivatives 4.0 International Licence (CC BY-ND)
Advisor
Allison, Peter
Piggott, Matthew
Czaja, Arnaud
Sponsor
Imperial College London
Publisher Department
Earth Science and Engineering and the Grantham Institute for Climate Change
Publisher Institution
Imperial College London
Qualification Level
Doctoral
Qualification Name
Doctor of Philosophy (PhD)