Issue |
A&A
Volume 555, July 2013
|
|
---|---|---|
Article Number | A148 | |
Number of page(s) | 12 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/201321725 | |
Published online | 17 July 2013 |
Stable higher order finite-difference schemes for stellar pulsation calculations
1
Institut d’Astrophysique et Géophysique de l’Université de
Liège,
Allée du 6 Août 17,
4000
Liège,
Belgium
e-mail: daniel.reese@ulg.ac.be
2
Observatoire de Paris, LESIA, CNRS UMR 8109,
92195
Meudon,
France
Received:
17
April
2013
Accepted:
5
June
2013
Context. Calculating stellar pulsations requires high enough accuracy to match the quality of the observations. Many current pulsation codes apply a second-order finite-difference scheme, combined with Richardson extrapolation to reach fourth-order accuracy on eigenfunctions. Although this is a simple and robust approach, a number of drawbacks exist that make fourth-order schemes desirable. A robust and simple finite-difference scheme that can easily be implemented in either 1D or 2D stellar pulsation codes, is therefore required.
Aims. One of the difficulties in setting up higher order finite-difference schemes for stellar pulsations is the so-called mesh-drift instability. Current ways of dealing with this defect include introducing artificial viscosity or applying a staggered grid approach; however, these remedies are not well-suited to eigenvalue problems, especially those involving non-dissipative systems, because they unduly change the spectrum of the operator, introduce supplementary free parameters, or lead to complications when applying boundary conditions.
Methods. We propose here a new method, inspired from the staggered grid strategy, which removes this instability while bypassing the above difficulties. Furthermore, this approach lends itself to superconvergence, a process in which the accuracy of the finite differences is boosted by one order.
Results. This new approach is successfully applied to stellar pulsation calculations, and is shown to be accurate, flexible with respect to the underlying grid, and able to remove mesh drift.
Conclusions. Although specifically designed for stellar pulsation calculations, this method can easily be applied to many other physical or mathematical problems.
Key words: methods: numerical / methods: analytical / stars: oscillations / stars: rotation
© ESO, 2013
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.