An algorithm for computing the 2D structure of fast rotating stars
From MaRDI portal
Publication:726841
DOI10.1016/J.JCP.2016.05.011zbMATH Open1349.85004arXiv1605.02359OpenAlexW2349617672MaRDI QIDQ726841FDOQ726841
Authors: M. Rieutord, Francisco Espinosa Lara, Bertrand Putigny 
Publication date: 5 December 2016
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.
Full work available at URL: https://arxiv.org/abs/1605.02359
Recommendations
- 2-D models of rapidly rotating stars
- 2-dimensional models of rapidly rotating stars I. Uniformly rotating zero age main sequence stars
- A multiple shooting approach for the numerical treatment of stellar structure and evolution
- scientific article; zbMATH DE number 3900383
- Highly accurate calculation of rotating neutron stars
Gas dynamics (general theory) (76N15) General theory of rotating fluids (76U05) Magnetohydrodynamics and electrohydrodynamics (76W05) Galactic and stellar structure (85A15)
Cites Work
- CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear systems
- Stellar structure and evolution.
- Spectral Methods
- Title not available (Why is that?)
- Jacobian-free Newton-Krylov methods: a survey of approaches and applications.
- Title not available (Why is that?)
- Implementation of the Jacobian-free Newton-Krylov method for solving the first-order ice sheet momentum balance
- The dynamics of a fully radiative rapidly rotating star enclosed within a spherical box
- 2-dimensional models of rapidly rotating stars I. Uniformly rotating zero age main sequence stars
- The dynamics of the radiative envelope of rapidly rotating stars
- CESAM
- Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method
- Stellar structure and evolution.
- Do eddington–sweet circulations exist?
- A finite element solutionof three-dimensional mixed convection gas flows in horizontal chnnels using preconditioned iterative metrix methods
Cited In (12)
- Stress-driven spin-down of a viscous fluid within a spherical shell
- A multiple shooting approach for the numerical treatment of stellar structure and evolution
- 2D solar modeling
- 2-D models of rapidly rotating stars
- Projected rotation velocities of several Algol-like irregular variable stars with H\(_{\alpha}\)-emission (ALIVARS)
- A finite volume code for meridional circulation in stars.
- From axiomatics of quantum probability to modelling geological uncertainty and management of intelligent hydrocarbon reservoirs with the theory of open quantum systems
- Title not available (Why is that?)
- Algorithmic Representation of Astrophysical Structures
- Title not available (Why is that?)
- Compressible fluid modes in rigid ellipsoids: towards modal acoustic velocimetry
- 2-dimensional models of rapidly rotating stars I. Uniformly rotating zero age main sequence stars
Uses Software
This page was built for publication: An algorithm for computing the 2D structure of fast rotating stars
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q726841)
