-
Effect of Rotation on Wave Mixing in Intermediate Mass Stars
Authors:
A. Varghese,
R. P. Ratnasingam,
R. Vanon,
P. V. F. Edelmann,
S. Mathis,
T. M. Rogers
Abstract:
Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$,…
▽ More
Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$, $2\times 10^{-5}$, $3\times 10^{-5}$, $4\times 10^{-5}$ and $1\times 10^{-4}$ rad.s$^{-1}$) and observe that the mixing decreases with decreasing Rossby number. This can be attributed to the effect of rotation on convection which influences the amplitude with which the waves are excited near the convective-radiative interface.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
Three-Dimensional Simulations of Massive Stars: II. Age Dependence
Authors:
R. Vanon,
P. V. F. Edelmann,
R. P. Ratnasingam,
A. Varghese,
T. M. Rogers
Abstract:
We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three $7M_\odot$ stars of different ages (ZAMS, midMS and TAMS). A comparison with several theoretical prescriptions shows the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly abo…
▽ More
We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three $7M_\odot$ stars of different ages (ZAMS, midMS and TAMS). A comparison with several theoretical prescriptions shows the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly above the convective boundary (CB); overshooting to the second, deeper layer becomes increasingly more infrequent with stellar age. Internal gravity wave (IGW) propagation is significantly impacted in the midMS and TAMS models as a result of some IGWs getting trapped within their Brunt-Väisälä frequency spikes. A fundamental change in the wave structure across radius is also observed, driven by the effect of density stratification on IGW propagation causing waves to become evanescent within the radiative zone, with older stars being affected more strongly. We find that the steepness of the frequency spectrum at the surface increases from ZAMS to the older models, with older stars also showing more modes in their spectra.
△ Less
Submitted 27 July, 2023;
originally announced July 2023.
-
Internal Gravity Waves in Massive Stars II: Frequency Analysis Across Stellar Mass
Authors:
R. P. Ratnasingam,
T. M. Rogers,
S. Chowdhury,
G. Handler,
R. Vanon,
A. Varghese,
P. V. F. Edelmann
Abstract:
Stars with masses above 1.6 solar masses generally possess convective cores and radiative envelopes, which allows the propagation of outward-travelling internal gravity waves. We have studied the generation and propagation of IGWs in such stars using two-dimensional, fully nonlinear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code,…
▽ More
Stars with masses above 1.6 solar masses generally possess convective cores and radiative envelopes, which allows the propagation of outward-travelling internal gravity waves. We have studied the generation and propagation of IGWs in such stars using two-dimensional, fully nonlinear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for 5 different stellar masses at the middle of main sequence: 3 - 13 solar masses. From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compare our simulation results with stellar photometric data from a recent survey and find that for frequency intervals above 8 microHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars, indicating that the brightness variations are caused by core-generated IGWs.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Chemical Mixing Induced by Internal Gravity Waves in Intermediate Mass Stars
Authors:
A. Varghese,
R. P. Ratnasingam,
R. Vanon,
P. V. F. Edelmann,
T. M. Rogers
Abstract:
Internal gravity waves (IGWs) can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into two - dimensional (2D) hydrodynamic simulations. Following the work of Rogers & McElwaine (2017), arXiv:1709.04920, we extend our study to different masses (3 M$_{\odot}$, 7 M$_{\odot}$ and 20 M$_{\odot}$) and ages (ZAMS, midMS and TAMS). The diffusion profi…
▽ More
Internal gravity waves (IGWs) can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into two - dimensional (2D) hydrodynamic simulations. Following the work of Rogers & McElwaine (2017), arXiv:1709.04920, we extend our study to different masses (3 M$_{\odot}$, 7 M$_{\odot}$ and 20 M$_{\odot}$) and ages (ZAMS, midMS and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt-Väisälä frequency, density, thermal damping, the geometric effect and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases towards the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in one-dimensional (1D) stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.
△ Less
Submitted 11 November, 2022;
originally announced November 2022.
-
The role of zonal flows in disc gravito-turbulence
Authors:
Riccardo Vanon
Abstract:
The work presented here focuses on the role of zonal flows in the self-sustenance of gravito-turbulence in accretion discs. The numerical analysis is conducted using a bespoke pseudo-spectral code in fully compressible, non-linear conditions. The disc in question, which is modelled using the shearing sheet approximation, is assumed to be self-gravitating, viscous, and thermally diffusive; a consta…
▽ More
The work presented here focuses on the role of zonal flows in the self-sustenance of gravito-turbulence in accretion discs. The numerical analysis is conducted using a bespoke pseudo-spectral code in fully compressible, non-linear conditions. The disc in question, which is modelled using the shearing sheet approximation, is assumed to be self-gravitating, viscous, and thermally diffusive; a constant cooling timescale is also considered. Zonal flows are found to emerge at the onset of gravito-turbulence and they remain closely linked to the turbulent state. A cycle of zonal flow formation and destruction is established, mediated by a slow mode instability (which allows zonal flows to grow) and a non-axisymmetric instability (which disrupts the zonal flow), which is found to repeat numerous times. It is in fact the disruptive action of the non-axisymmetric instability to form new leading and trailing shearing waves, allowing energy to be extracted from the background flow and ensuring the self-sustenance of the gravito-turbulent regime.
△ Less
Submitted 14 April, 2018;
originally announced April 2018.
-
Zonal flow evolution and overstability in accretion discs
Authors:
Riccardo Vanon,
Gordon Ogilvie
Abstract:
This work presents a linear analytical calculation on the stability and evolution of a compressible, viscous self-gravitating (SG) Keplerian disc with both horizontal thermal diffusion and a constant cooling timescale when an axisymmetric structure is present and freely evolving. The calculation makes use of the shearing sheet model and is carried out for a range of cooling times. Although the sol…
▽ More
This work presents a linear analytical calculation on the stability and evolution of a compressible, viscous self-gravitating (SG) Keplerian disc with both horizontal thermal diffusion and a constant cooling timescale when an axisymmetric structure is present and freely evolving. The calculation makes use of the shearing sheet model and is carried out for a range of cooling times. Although the solutions to the inviscid problem with no cooling or diffusion are well known, it is non-trivial to predict the effect caused by the introduction of cooling and of small diffusivities; this work focuses on perturbations of intermediate wavelengths, therefore representing an extension to the classical stability analysis on thermal and viscous instabilities. For density wave modes the analysis can be simplified by means of a regular perturbation analysis; considering both shear and thermal diffusivities, the system is found to be overstable for intermediate and long wavelengths for values of the Toomre parameter $Q \lesssim 2$; a non-SG instability is also detected for wavelengths $\gtrsim 18H$, where $H$ is the disc scale height, as long as $γ\lesssim 1.305$. The regular perturbation analysis does not however hold for the entropy and potential vorticity slow modes as their ideal growth rates are degenerate. To understand their evolution, equations for the axisymmetric structure's amplitudes in these two quantities are analytically derived and their instability regions obtained. The instability appears boosted by increasing the value of the adiabatic index and of the Prandtl number, while it is quenched by efficient cooling.
△ Less
Submitted 13 December, 2016;
originally announced December 2016.
-
Non-axisymmetric instabilities in discs with imposed zonal flows
Authors:
R. Vanon,
G. I. Ogilvie
Abstract:
We conduct a linear stability calculation of an ideal Keplerian flow on which a sinusoidal zonal flow is imposed. The analysis uses the shearing sheet model and is carried out both in isothermal and adiabatic conditions, with and without self-gravity (SG). In the non-SG regime a structure in the potential vorticity (PV) leads to a non-axisymmetric Kelvin-Helmholtz (KH) instability; in the short-wa…
▽ More
We conduct a linear stability calculation of an ideal Keplerian flow on which a sinusoidal zonal flow is imposed. The analysis uses the shearing sheet model and is carried out both in isothermal and adiabatic conditions, with and without self-gravity (SG). In the non-SG regime a structure in the potential vorticity (PV) leads to a non-axisymmetric Kelvin-Helmholtz (KH) instability; in the short-wavelength limit its growth rate agrees with the incompressible calculation by Lithwick (2007), which only considers perturbations elongated in the streamwise direction. The instability's strength is analysed as a function of the structure's properties, and zonal flows are found to be stable if their wavelength is $\gtrsim 8H$, where $H$ is the disc's scale height, regardless of the value of the adiabatic index $γ$. The non-axisymmetric KH instability can operate in Rayleigh-stable conditions, and it therefore represents the limiting factor to the structure's properties. Introducing SG triggers a second non-axisymmetric instability, which is found to be located around a PV maximum, while the KH instability is linked to a PV minimum, as expected. In the adiabatic regime, the same gravitational instability is detected even when the structure is present only in the entropy (not in the PV) and the instability spreads to weaker SG conditions as the entropy structure's amplitude is increased. This eventually yields a non-axisymmetric instability in the non-SG regime, albeit of weak strength, localised around an entropy maximum.
△ Less
Submitted 4 September, 2017; v1 submitted 8 September, 2016;
originally announced September 2016.