Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe.Stars born with masses larger than approximately eight times the mass of the Sun play a significant role in the evolution of galaxies. They are the chemical factories that produce and expel heavy elements through their wind and when they end their lives as supernovae and form a black hole or neutron star 1-3 . However, the chemical yields that enrich the interstellar medium and the remnant mass strongly depend on the progenitor star's interior properties 4 . The detectable progenitors of supernovae include blue supergiant stars, which are hot massive stars in a shell-hydrogen or core-helium burning stage of stellar evolution. Stellar evolution models of these post-main sequence stars contain by far the largest uncertainties in stellar astrophysics, as observational constraints on interior mixing, rotation and angular momentum transport are missing. These phenomena are further compounded when coupled with mass loss, binarity and magnetic fields 1-3 . Across astrophysics, from star formation to galactic evolution, it is imperative to calibrate theoretical models of massive stars using observations because they determine the evolution of the cosmos.A unique methodology for probing stellar interiors is asteroseismology 5 , which -similarly to seismology of earthquakes -uses oscillations to derive constraints on the structure of stars.The study of stellar interiors of low-mass stars like the Sun has undergone a revolution in

We present the first three-dimensional (3D), hydrodynamic simulations of the core convection zone (CZ) and extended radiative zone spanning from 1% to 90% of the stellar radius of an intermediate mass (3 M ) star. This allows us to selfconsistently follow the generation of internal gravity waves (IGWs) at the convective boundary and their propagation to the surface. We find that convection in the core is dominated by plumes. The frequency spectrum in the CZ and that of IGW generation is a double power law as seen in previous two-dimensional (2D) simulations. The spectrum is significantly flatter than theoretical predictions using excitation through Reynolds stresses induced by convective eddies alone. It is compatible with excitation through plume penetration. An empirically determined distribution of plume frequencies generally matches the one necessary to explain a large part of the observed spectrum. We observe waves propagating in the radiation zone and excited standing modes, which can be identified as gravity and fundamental modes. They show similar frequencies and node patterns to those predicted by the stellar oscillation code GYRE. The continuous part of the spectrum fulfills the IGW dispersion relation. A spectrum of tangential velocity and temperature fluctuations close to the surface is extracted, which are directly related to observable brightness variations in stars. Unlike 2D simulations we do not see the high frequencies associated with wave breaking, likely because these 3D simulations are more heavily damped.

Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

Internal gravity waves (IGW) propagate in the radiation zones of all stars. During propagation, their amplitudes are affected by two main features: radiative diffusion and density stratification. We have studied the implications of these two features on waves traveling within the radiative zones of non-rotating stars with stellar parameters obtained from the one dimensional stellar evolution code, MESA. As a simple measure of induced wave dynamics, we define a criterion to see if waves can become nonlinear and if so, under what conditions. This was done to understand the role IGW may play in angular momentum transport and mixing within stellar interiors. We find that the IGW generation spectrum, convective velocities and the strength of density stratification all play major roles in whether waves become nonlinear. With increasing stellar mass, there is an increasing trend in nonlinear wave energies. The trends with different metallicities and ages depend on the generation spectrum.

Context. Large-scale vortices in protoplanetary disks are thought to form and survive for long periods of time. Hence, they can significantly change the global disk evolution and particularly the distribution of the solid particles embedded in the gas, possibly explaining asymmetries and dust concentrations recently observed at submillimeter and millimeter wavelengths. Aims. We investigate the spatial distribution of dust grains using a simple model of protoplanetary disk hosted by a giant gaseous vortex. We explore the dependence of the results on grain size and deduce possible consequences and predictions for observations of the dust thermal emission at submillimeter and millimeter wavelengths. Methods. Global 2D simulations with a bi-fluid code are used to follow the evolution of a single population of solid particles aerodynamically coupled to the gas. Possible observational signatures of the dust thermal emission are obtained using simulators of ALMA and Nest Generation Very Large Array (ngVLA) observations. Results. We find that a giant vortex not only captures dust grains with Stokes number S t < 1 but can also affect the distribution of larger grains (with S t ∼ 1) carving a gap associated with a ring composed of incompletely trapped particles. The results are presented for different particle sizes and associated with their possible signatures in disk observations. Conclusions. Gap clearing in the dust spatial distribution could be due to the interaction with a giant gaseous vortex and their associated spiral waves without the gravitational assistance of a planet. Hence, strong dust concentrations at short sub-mm wavelengths associated with a gap and an irregular ring at longer mm and cm wavelengths could indicate the presence of an unseen gaseous vortex.