ترغب بنشر مسار تعليمي؟ اضغط هنا

An algorithm for computing the 2D structure of fast rotating stars

80   0   0.0 ( 0 )
 نشر من قبل Michel L. E. Rieutord
 تاريخ النشر 2016
  مجال البحث فيزياء
والبحث باللغة English
 تأليف M. Rieutord




اسأل ChatGPT حول البحث

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.



قيم البحث

اقرأ أيضاً

Many early-type stars have been measured with high angular velocities. In such stars, mode identification is difficult as the effects of fast and differential rotation are not well known. Using fundamental parameters measured by interferometry, the E STER structure code and the TOP oscillation code, we investigate the oscillation spectrum of Rasalhague (alpha Ophiuchi), for which observations by the MOST satellite found 57 oscillations frequencies. Results do not show a clear identification of the modes and highlight the difficulties of asteroseismology for such stars with a very complex oscillation spectrum.
Early-type stars generally tend to be fast rotators. In these stars, mode identification is very challenging as the effects of rotation are not well known. We consider here the example of $alpha$ Ophiuchi, for which dozens of oscillation frequencies have been measured. We model the star using the two-dimensional structure code ESTER, and we compute both adiabatic and non-adiabatic oscillations using the TOP code. Both calculations yield very complex spectra, and we used various diagnostic tools to try and identify the observed pulsations. While we have not reached a satisfactory mode-to-mode identification, this paper presents promising early results.
In our previous study of low mass stars using TESS, we found a handful which show a periodic modulation on a period <1 d but also displayed no flaring activity. Here we present the results of a systematic search for Ultra Fast Rotators (UFRs) in the southern ecliptic hemisphere which were observed in 2 min cadence with TESS. Using data from Gaia DR2, we obtain a sample of over 13,000 stars close to the lower main sequence. Of these, we identify 609 stars which lie on the lower main sequence and have a periodic modulation <1 d. The fraction of stars which show flares appears to drop significantly at periods <0.2 d. If the periods are a signature of the rotation rate, this would be a surprise, since faster rotators would be expected to have a stronger magnetic field and, therefore, produce more flares. We explore possible reasons for our finding: the flare inactive stars are members of binaries, in which case the stars rotation rate could have increased as the binary orbital separation reduced due to angular momentum loss over time, or that enhanced emission occurs at blue wavelengths beyond the pass band of TESS. Follow-up spectroscopy and flare monitoring at blue/ultraviolet wavelengths of these flare inactive stars are required to resolve this question.
We present a new fast algorithm for computing the Boys function using nonlinear approximation of the integrand via exponentials. The resulting algorithms evaluate the Boys function with real and complex valued arguments and are competitive with previously developed algorithms for the same purpose.
A recent photometric survey in the NGC~3766 cluster led to the detection of stars presenting an unexpected variability. They lie in a region of the Hertzsprung-Russell (HR) diagram where no pulsation are theoretically expected, in between the $delta$ Scuti and slowly pulsating B (SPB) star instability domains. Their variability periods, between $sim$0.1--0.7~d, are outside the expected domains of these well-known pulsators. The NCG~3766 cluster is known to host fast rotating stars. Rotation can significantly affect the pulsation properties of stars and alter their apparent luminosity through gravity darkening. Therefore we inspect if the new variable stars could correspond to fast rotating SPB stars. We carry out instability and visibility analysis of SPB pulsation modes within the frame of the traditional approximation. The effects of gravity darkening on typical SPB models are next studied. We find that at the red border of the SPB instability strip, prograde sectoral (PS) modes are preferentially excited, with periods shifted in the 0.2--0.5~d range due to the Coriolis effect. These modes are best seen when the star is seen equator-on. For such inclinations, low-mass SPB models can appear fainter due to gravity darkening and as if they were located between the $delta$~Scuti and SPB instability strips.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا