Numerical approach for high precision 3D relativistic star models (original) (raw)
Related papers
A new numerical approach to the oscillation modes of relativistic stars
Arxiv preprint gr-qc/9503014, 1995
The oscillation modes of a simple polytropic stellar model are studied. Using a new numerical approach (based on integration for complex coordinates) to the problem for the stellar exterior we have computed the eigenfrequencies of the highly damped w-modes. The results obtained agree well with recent ones of Leins, Nollert and Soffel [Phys. Rev. D 48 3467 (1993)]. Specifically, we are able to explain why several modes in this regime of the complex frequency plane could not be identified within the WKB approach of Kokkotas and Schutz [Mon. Not. R. astr. Soc. 255 119 ]. Furthermore, we have established that the "kink" that was a prominent feature of the spectra of Kokkotas and Schutz, but did not appear in the results of Leins et al., was a numerical artefact. Using our new numerical code we are also able to compute, for the first time, several of the slowly damped (p) modes for the considered stellar models. For very compact stars we find, somewhat surprisingly, that the damping of these modes does not decrease monotonically as one proceeds to higher oscillation frequencies. The existence of low-order modes that damp away much faster than anticipated may have implications for questions regarding stellar stability and the lifetime of gravitational-wave sources. The present results illustrate the accuracy and reliability of the complex-coordinate method and indicate that the method could prove to be of great use also in problems involving rotating stars. There is no apparent reason why the complex-coordinate approach should not extend to rotating stars, whereas it is accepted that all previous methods will fail to do so.
Physical Review D, 2002
This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.
Phys Rev D, 2002
This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The following aspects of our code have been validated: the generation of initial data representing perturbed general relativistic polytropic models (both rotating and nonrotating), the long-term evolution of relativistic stellar models, and the coupling of our evolution code to analysis modules providing, for instance, the detection of apparent horizons or the extraction of gravitational waveforms. The tests involve single nonrotating stars in stable equilibrium, nonrotating stars undergoing radial and quadrupolar oscillations, nonrotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, nonrotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasiradial oscillations. We have carried out evolutions in full general relativity and compared the results to those obtained either with perturbation techniques, or with lower dimensional numerical codes, or in the Cowling approximation (in which all the perturbations of the spacetime are neglected). In all cases an excellent agreement has been found. The numerical evolutions have been carried out using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. In particular, we have found the ``monotonized central differencing'' limiter to be particularly effective in evolving the relativistic stellar models considered. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.
Modeling rotating stars in two dimensions
EAS Publications Series, 2013
In this lecture I present the way stars can be modeled in two dimensions and especially the fluid flows that are driven by rotation. I discuss some of the various ways of taking into account turbulence and conclude this contribution by a short presentation of some of the first results obtained with the ESTER code on the modeling of interferometrically observed fast rotating early-type stars.
2D modelling of pulsating stars with rapid rotation
Frontiers in Astronomy and Space Sciences
Rapid stellar rotation is an important phenomenon in stellar physics, particularly for massive and intermediate mass main-sequence stars. This affects all aspects of the star’s physics including its structure, evolution, and pulsations, and makes it necessary to use 2D numerical approaches rather than the 1D approaches typically used. In this contribution, we will review 2D numerical methods for modelling and interpreting pulsation modes in rapidly rotating stars. We will start by deriving the pulsation equations, both in an adiabatic and non-adiabatic setting, then provide a description of the 2D numerical implementation. We will then explain approximate implementations of the effects of rotation, namely first, second, and third order perturbative approaches, as well as the traditional approximation. This will then be followed by a description on how to calculate disk-integrated mode visibilities in various photometric bands, and how to apply this to mode identification in rapid ro...
Self-consistent 2D models of fast-rotating early-type stars
Astronomy & Astrophysics, 2013
Aims. This work aims at presenting the first two-dimensional models of an isolated rapidly rotating star that include the derivation of the differential rotation and meridional circulation in a self-consistent way. Methods. We use spectral methods in multidomains, together with a Newton algorithm to determine the steady state solutions including differential rotation and meridional circulation for an isolated non-magnetic, rapidly rotating early-type star. In particular we devise an asymptotic method for small Ekman numbers (small viscosities) that removes the Ekman boundary layer and lifts the degeneracy of the inviscid baroclinic solutions. Results. For the first time, realistic two-dimensional models of fast-rotating stars are computed with the actual baroclinic flows that predict the differential rotation and the meridional circulation for intermediate-mass and massive stars. These models nicely compare with available data of some nearby fast-rotating early-type stars like Ras Alhague (α Oph), Regulus (α Leo), and Vega (α Lyr). It is shown that baroclinicity drives a differential rotation with a slow pole, a fast equator, a fast core, and a slow envelope. The differential rotation is found to increase with mass, with evolution (here measured by the hydrogen mass fraction in the core), and with metallicity. The core-envelope interface is found to be a place of strong shear where mixing will be efficient. Conclusions. Two-dimensional models offer a new view of fast-rotating stars, especially of their differential rotation, which turns out to be strong at the core-envelope interface. They also offer more accurate models for interpreting the interferometric and spectroscopic data of early-type stars.
An algorithm for computing the 2D structure of fast rotating stars: Part I steady solutions
Journal of Computational Physics, 2016
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.
2D Computations of g-modes in Fast Rotating Stars
2012
We present complete 2D computations of g modes in distorted polytropic models of stars performed with the Two-dimensional Oscillation Program (TOP). We computed low-degree modes (ℓ = 1 modes with radial order n = −1 . . . −14, and ℓ = 2, 3 modes with n = −1 . . . −5 and −16 . . . −20) of a nonrotating model and followed them by slowly increasing the rotation rate up to 70 % of the Keplerian break-up velocity. We use these computations to determine the domain of validity of perturbative methods up to the 3rd order. We study the evolution of the regularities of the spectrum and show quantitative agreement with the traditional approximation for not too large values of the ratio of the rotation rate to the pulsation frequency. We also show the appearance of new types of modes, called "rosette" modes due to their spatial structure. Thanks to the ray theory for gravito-inertial waves that we developed, we can associate these modes with stable periodic rays.
Numerical calculation of linear modes in stellar disks
We present a method for solving the twodimensional linearized collisionless Boltzmann equation using Fourier expansion along the orbits. It resembles very much solutions present in the literature, but it differs by the fact that everything is performed in coordinate space instead of using action-angle variables. We show that this approach, though less elegant, is both feasible and straightforward. This approach is then incorporated in a matrix method in order to calculate self-consistent modes, using a set of potential-density pairs which is obtained numerically. We investigated the stability of some unperturbed disks having an almost flat rotation curve, an exponential disk and a non-zero velocity dispersion. The influence of the velocity dispersion, halo mass and anisotropy on the stability is further discussed.