Volume 631, November 2019
|Number of page(s)||15|
|Section||Numerical methods and codes|
|Published online||16 October 2019|
Eigenvectors, Circulation, and Linear Instabilities for Planetary Science in 3 Dimensions (ECLIPS3D)
Ecole Normale Supérieure de Lyon, CRAL, UMR CNRS 5574, 69364 Lyon Cedex 07, France
2 School of Physics and Astronomy, University of Exeter, Exeter EX4 4QL, UK
3 IRAP, Université de Toulouse, CNRS, UPS, Toulouse, France
4 College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter EX4 4QF, UK
Accepted: 3 September 2019
Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using global circulation models (GCMs). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results.
Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states indicated to be essential in existing theories of the acceleration of hot Jupiter superrotation.
Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode.
Results. Our code is tested using classical wave and instability test cases in multiple geometries (2D, 3D, two-layer equivalent depth). The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes.
Conclusions. We developed and tested a code for the study of linear processes in planetary atmospheres with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably used to further study the initial spin up of superrotation of GCM simulations of hot Jupiters.
Key words: methods: numerical / hydrodynamics / waves / instabilities / planets and satellites: atmospheres
© F. Debras et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.