Imperial College London
Portrait Picture

ProfessorColinCotter

Faculty of Natural Sciences, Department of Mathematics

Professor of Computational Mathematics
 
 
 
//

Contact

 

+44 (0)20 7594 3468colin.cotter

 
 
//

Location

 

755Huxley BuildingSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Melvin:2018:10.1002/qj.3268,
author = {Melvin, T and Benacchio, T and Thuburn, J and Cotter, C},
doi = {10.1002/qj.3268},
journal = {Quarterly Journal of the Royal Meteorological Society},
pages = {900--916},
title = {Choice of function spaces for thermodynamic variables in mixed finite-element methods},
url = {http://dx.doi.org/10.1002/qj.3268},
volume = {144},
year = {2018}
}
Download

RIS format (EndNote, RefMan)

TY  - JOUR
AB  - We study the dispersion properties of three choices for the buoyancy space in a mixed finiteelement discretization of geophysical fluid flow equations. The problem is analogous to that of the staggering of the buoyancy variable in finitedifference discretizations. Discrete dispersion relations of the twodimensional linear gravity wave equations are computed. By comparison with the analytical result, the best choice for the buoyancy space basis functions is found to be the horizontally discontinuous, vertically continuous option. This is also the space used for the vertical component of the velocity. At lowest polynomial order, this arrangement mirrors the Charney–Phillips vertical staggering known to have good dispersion properties in finitedifference models. A fully discontinuous space for the buoyancy corresponding to the Lorenz finitedifference staggering at lowest order gives zero phase velocity for high vertical wavenumber modes. A fully continuous space, the natural choice for scalar variables in a mixed finiteelement framework, with degrees of freedom of buoyancy and vertical velocity horizontally staggered at lowest order, is found to entail zero phase velocity modes at the large horizontal wavenumber end of the spectrum. Corroborating the theoretical insights, numerical results obtained on gravity wave propagation with fully continuous buoyancy highlight the presence of a computational mode in the poorly resolved part of the spectrum that fails to propagate horizontally. The spurious signal is not removed in test runs with higherorder polynomial basis functions. Runs at higher order also highlight additional oscillations, an issue that is shown to be mitigated by partial masslumping. In light of the findings and with a view to coupling the dynamical core to physical parametrizations that often force near the horizontal grid scale, the use of the fully continuous space should be avoided in favour of the horizontally discontinuous, vertically co
AU  - Melvin,T
AU  - Benacchio,T
AU  - Thuburn,J
AU  - Cotter,C
DO  - 10.1002/qj.3268
EP  - 916
PY  - 2018///
SN  - 0035-9009
SP  - 900
TI  - Choice of function spaces for thermodynamic variables in mixed finite-element methods
T2  - Quarterly Journal of the Royal Meteorological Society
UR  - http://dx.doi.org/10.1002/qj.3268
UR  - http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=PARTNER_APP&SrcAuth=LinksAMR&KeyUT=WOS:000443007800020&DestLinkType=FullRecord&DestApp=ALL_WOS&UsrCustomerID=1ba7043ffcc86c417c072aa74d649202
UR  - http://hdl.handle.net/10044/1/63111
VL  - 144
ER  -
Download