Search NASA⌕ Search

Engineering topics

Baraffe, I.

Publications and source records attributed to Baraffe, I..

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of stratification on overshooting and waves atop the convective core of M ⊙ main-sequence stars

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a $5\,\mathrm{ M}_{\odot }$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2017

Cloudless Atmospheres for Young Low-Gravity Substellar Objects

Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains challenging. The presence of very thick clouds is a possible source of this challenge, because of their extremely red near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with the evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provide a very good model to match the observed VL-G NIR spectra. The sequence of extremely red colors in the NIR for atmospheres with effective temperatures from approx. 2000 K down to approx. 1200 K is very well reproduced with predicted radii typical of young low-gravity objects. Future observations with NIRSPEC and MIRI on the James Webb Space Telescope (JWST) will provide more constraints in the mid-infrared, helping to confirm or refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence or absence of clouds will be directly determined by the silicate absorption features that can be observed with MIRI. JWST will therefore be able to better characterize the atmosphere of these hot young brown dwarfs and their low-gravity exoplanet analogs.

Tremblin, P.↗
NASA NTRS2017

Cloudless Atmospheres for Young Low-Gravity Substellar Objects

Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains challenging. The presence of very thick clouds is a possible source of this challenge, because of their extremely red near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with the evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provide a very good model to match the observed VL-G NIR spectra. The sequence of extremely red colors in the NIR for atmospheres with effective temperatures from ∼2000 K down to ∼1200 K is very well reproduced with predicted radii typical of young low-gravity objects. Future observations with NIRSPEC and MIRI on the James Webb Space Telescope (JWST) will provide more constraints in the mid-infrared, helping to confirm or refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence or absence of clouds will be directly determined by the silicate absorption features that can be observed with MIRI. JWST will therefore be able to better characterize the atmosphere of these hot young brown dwarfs and their low-gravity exoplanet analogs.

Tremblin, P.↗
NASA NTRS2017

Advection of Potential Temperature in the Atmosphere of Irradiated Exoplanets: A Robust Mechanism to Explain Radius Inflation

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a two-dimensional steady-state atmospheric circulation model, the validity of which is assessed by comparison to three-dimensional calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occurring in the Earth's atmosphere or oceans. In the deep atmosphere, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD 209458b, but also reproduce the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non-uniform atmospheric heating thus provides a robust mechanism to explain the inflated radii of irradiated hot Jupiters.

tmospheric effects; methods: numerical; planets an↗
NASA NTRS2016

Cloudless Atmospheres for L/T Dwarfs and Extrasolar Giant Planets

The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BDs) since their first detection 20 years ago has always been the key role played by micron-size condensates, called "dust" or "clouds," in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this Letter, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/ CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near-infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.

Brown dwarf stars↗
NASA NTRS2016

Near-Infrared Spectroscopy of the Y0 WISEP J173835.52+273258.9 and the Y1 WISE J035000.32-565830.2: the Importance of Non-Equilibrium Chemistry

We present new near-infrared spectra, obtained at Gemini Observatory, for two Y dwarfs: WISE J035000.32−565830.2 (W0350) and WISEP J173835.52+273258.9 (W1738). A FLAMINGOS-2 R = 540 spectrum was obtained for W0350, covering 1.0 < λ micrometer < 1.7, and a cross-dispersed GNIRS R = 2800 spectrum was obtained for W1738, covering 0.993 - 1.087 micrometer, 1.191 - 1.305 micrometer, 1.589 - 1.631 micrometer, and 1.985 - 2.175 micrometer, in four orders. We also present revised Y JH photometry for W1738, using new NIRI Y and J imaging, and a re-analysis of the previously published NIRI H band images. We compare these data, together with previously published data for late-T and Y dwarfs, to cloud-free models of solar metallicity, calculated both in chemical equilibrium and with disequilibrium driven by vertical transport. We find that for the Y dwarfs the non-equilibrium models reproduce the near-infrared data better than the equilibrium models. The remaining discrepancies suggest that fine-tuning the CH4/CO and NH3/N2 balance is needed. Improved trigonometric parallaxes would improve the analysis. Despite the uncertainties and discrepancies, the models reproduce the observed near-infrared spectra well. We find that for the Y0, W1738, T(sub eff) = 425 +/- 25K and log g = 4.0 +/- 0.25, and for the Y1, W0350, T(sub eff) = 350 +/- 25K and log g = 4.0 +/- 0.25. W1738 may be metal-rich. Based on evolutionary models, these temperatures and gravities correspond to a mass range for both Y dwarfs of 3 - 9 Jupiter masses, with W0350 being a cooler, slightly older, version of W1738; the age of W0350 is 0.3 - 3 Gyr, and the age of W1738 is 0.15 - 1 Gyr.

two Y dwarfs↗