Issue |
A&A
Volume 555, July 2013
|
|
---|---|---|
Article Number | A81 | |
Number of page(s) | 12 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/201220725 | |
Published online | 04 July 2013 |
Comparison of different nonlinear solvers for 2D time-implicit stellar hydrodynamics
1
Max-Planck-Institut für Astrophysik, Karl Schwarzschild Strasse 1, 85741
Garching, Germany
e-mail: mviallet@mpa-garching.mpg.de
2
Physics and Astronomy, University of Exeter,
Stocker Road, Exeter, EX4
4QL, UK
3
École Normale Supérieure, Lyon, CRAL (UMR CNRS 5574), Université
de Lyon 1, 69007
Lyon,
France
Received:
12
November
2012
Accepted:
11
May
2013
Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, at a cost: the solution of a system of nonlinear equations is required at each time step. In this work, the nonlinear system results from the discretization of the hydrodynamical equations with the Crank-Nicholson scheme. We compare the cost of different methods, based on Newton-Raphson iterations, to solve this nonlinear system, and benchmark their performances against time-explicit schemes. Since our general scientific objective is to model stellar interiors, we use as test cases two realistic models for the convective envelope of a red giant and a young Sun. Focusing on 2D simulations, we show that the best performances are obtained with the quasi-Newton method proposed by Broyden. Another important concern is the accuracy of implicit calculations. Based on the study of an idealized problem, namely the advection of a single vortex by a uniform flow, we show that there are two aspects: i) the nonlinear solver has to be accurate enough to resolve the truncation error of the numerical discretization; and ii) the time step has be small enough to resolve the advection of eddies. We show that with these two conditions fulfilled, our implicit methods exhibit similar accuracy to time-explicit schemes, which have lower values for the time step and higher computational costs. Finally, we discuss in the conclusion the applicability of these methods to fully implicit 3D calculations.
Key words: hydrodynamics / methods: numerical / stars: interiors
© 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.