No Arabic abstract
The spectrum of oscillation modes of a star provides information not only about its material properties (e.g. mean density), but also its symmetries. Spherical symmetry can be broken by rotation and/or magnetic fields. It has been postulated that strong magnetic fields in the cores of some red giants are responsible for their anomalously weak dipole mode amplitudes (the dipole dichotomy problem), but a detailed understanding of how gravity waves interact with strong fields is thus far lacking. In this work, we attack the problem through a variety of analytical and numerical techniques, applied to a localised region centred on a null line of a confined axisymmetric magnetic field which is approximated as being cylindrically symmetric. We uncover a rich variety of phenomena that manifest when the field strength exceeds a critical value, beyond which the symmetry is drastically broken by the Lorentz force. When this threshold is reached, the spatial structure of the g-modes becomes heavily altered. The dynamics of wave packet propagation transitions from regular to chaotic, which is expected to fundamentally change the organisation of the mode spectrum. In addition, depending on their frequency and the orientation of field lines with respect to the stratification, waves impinging on different parts of the magnetised region are found to undergo either reflection or trapping. Trapping regions provide an avenue for energy loss through Alfven wave phase mixing. Our results may find application in various astrophysical contexts, including the dipole dichotomy problem, the solar interior, and compact star oscillations.
Recent photometric observations of massive stars show ubiquitous low-frequency red-noise variability, which has been interpreted as internal gravity waves (IGWs). Simulations of IGWs generated by convection show smooth surface wave spectra, qualitatively matching the observed red-noise. On the other hand, theoretical calculations by Shiode et al (2013) and Lecoanet et al (2019) predict IGWs should manifest at the surface as regularly-spaced peaks associated with standing g-modes. In this work, we compare these theoretical approaches to simplified 2D numerical simulations. The simulations show g-mode peaks at their surface, and are in good agreement with Lecoanet et al (2019). The amplitude estimates of Shiode et al (2013) did not take into account the finite width of the g-mode peaks; after correcting for this finite width, we find good agreement with simulations. However, simulations need to be run for hundreds of convection turnover times for the peaks to become visible; this is a long time to run a simulation, but a short time in the life of a star. The final spectrum can be predicted by calculating the wave energy flux spectrum in much shorter simulations, and then either applying the theory of Shiode et al (2013) or Lecoanet et al (2019).
Detailed modeling of stellar evolution requires a better understanding of the (magneto-)hydrodynamic processes which mix chemical elements and transport angular momentum. Understanding these pro- cesses is crucial if we are to accurately interpret observations of chemical abundance anomalies, surface rotation measurements and asteroseismic data. Here, we use two-dimensional hydrodynamic simula- tions of the generation and propagation of internal gravity waves (IGW) in an intermediate mass star to measure the chemical mixing induced by these waves. We show that such mixing can generally be treated as a diffusive process. We then show that the local diffusion coefficient does not depend on the local fluid velocity, but rather on the wave amplitude. We then use these findings to provide a simple parametrization for this diffusion which can be incorporated into stellar evolution codes and tested against observations.
Typical flows in stellar interiors are much slower than the speed of sound. To follow the slow evolution of subsonic motions, various sound-proof equations are in wide use, particularly in stellar astrophysical fluid dynamics. These low-Mach number equations include the anelastic equations. Generally, these equations are valid in nearly adiabatically stratified regions like stellar convection zones, but may not be valid in the sub-adiabatic, stably stratified stellar radiative interiors. Understanding the coupling between the convection zone and the radiative interior is a problem of crucial interest and may have strong implications for solar and stellar dynamo theories as the interface between the two, called the tachocline in the Sun, plays a crucial role in many solar dynamo theories. Here we study the properties of gravity waves in stably-stratified atmospheres. In particular, we explore how gravity waves are handled in various sound-proof equations. We find that some anelastic treatments fail to conserve energy in stably-stratified atmospheres, instead conserving pseudo-energies that depend on the stratification, and we demonstrate this numerically. One anelastic equation set does conserve energy in all atmospheres and we provide recommendations for converting low-Mach number anelastic codes to this set of equations.
Previous numerical studies have identified phase mixing of low-frequency Alfven waves as a mean of parallel electric field amplification and acceleration of electrons in a collisionless plasma. Theoretical explanations are given of how this produces an amplification of the parallel electric field, and as a consequence, also leads to enhanced collisionless damping of the wave by energy transfer to the electrons. Our results are based on the properties of the Alfven waves in a warm plasma which are obtained from drift-kinetic theory, in particular, the rate of their electron Landau damping. Phase mixing in a collisionless low-$beta$ plasma proceeds in a manner very similar to the visco-resistive case, except for the fact that electron Landau damping is the primary energy dissipation channel. The time and length scales involved are evaluated. We also focus on the evolution of the parallel electric field and calculate its maximum value in the course of its amplification.
Zonal flows in rotating systems have been previously shown to be suppressed by the imposition of a background magnetic field aligned with the direction of rotation. Understanding the physics behind the suppression may be important in systems found in astrophysical fluid dynamics, such as stellar interiors. However, the mechanism of suppression has not yet been explained. In the idealized setting of a magnetized beta plane, we provide a theoretical explanation that shows how magnetic fluctuations directly counteract the growth of weak zonal flows. Two distinct calculations yield consistent conclusions. The first, which is simpler and more physically transparent, extends the Kelvin-Orr shearing wave to include magnetic fields and shows that weak, long-wavelength shear flow organizes magnetic fluctuations to absorb energy from the mean flow. The second calculation, based on the quasilinear, statistical CE2 framework, is valid for arbitrary wavelength zonal flow and predicts a self-consistent growth rate of the zonal flow. We find that a background magnetic field suppresses zonal flow if the bare Alfven frequency is comparable to or larger than the bare Rossby frequency. However, suppression can occur for even smaller magnetic fields if the resistivity is sufficiently small enough to allow sizable magnetic fluctuations. Our calculations reproduce the $eta/B_0^2 = text{const.}$ scaling that describes the boundary of zonation, as found in previous work, and we explicitly link this scaling to the amplitude of magnetic fluctuations.