Aims.
Inspired by the statistical mechanics of an ensemble of interacting particles (BBGKY hierarchy), we propose to account for small-scale inhomogeneities in self-gravitating astrophysical fluids ...by deriving a nonideal virial theorem and nonideal Navier-Stokes equations. These equations involve the pair radial distribution function (similar to the two-point correlation function used to characterize the large-scale structures of the Universe), similarly to the interaction energy and equation of state in liquids. Within this framework, small-scale correlations lead to a nonideal amplification of the gravitational interaction energy, whose omission leads to a missing mass problem, for instance, in galaxies and galaxy clusters.
Methods.
We propose to use a decomposition of the gravitational potential into a near- and far-field component in order to account for the gravitational force and correlations in the thermodynamics properties of the fluid. Based on the nonideal virial theorem, we also propose an extension of the Friedmann equations in the nonideal regime and use numerical simulations to constrain the contribution of these correlations to the expansion and acceleration of the Universe.
Results.
We estimate that the nonideal amplification factor of the gravitational interaction energy of the baryons lies between 5 and 20, potentially explaining the observed value of the Hubble parameter (since the uncorrelated energy accounts for ∼5%). Within this framework, the acceleration of the expansion emerges naturally because the number of substructures induced by gravitational collapse increases, which in turn increases their contribution to the total gravitational energy. A simple estimate predicts a nonideal deceleration parameter
q
ni
≃ −1; this is potentially the first determination of the observed value based on an intuitively physical argument. We also suggest that small-scale gravitational interactions in bound structures (spiral arms or local clustering) could yield a transition to a viscous regime that can lead to flat rotation curves. This transition could also explain the dichotomy between (Keplerian) low surface brightness elliptical galaxy and (nonkeplerian) spiral galaxy rotation profiles. Overall, our results demonstrate that nonideal effects induced by inhomogeneities must be taken into account, potentially with our formalism, in order to properly determine the gravitational dynamics of galaxies and the large-scale Universe.
Aims. We present a new radiation hydrodynamics code called ARK-RT which uses a two-moment model with the M1 closure relation for radiative transfer. This code was designed to be ready for ...high-performance computing on exascale architectures. Methods. The two-moment model is solved using a finite-volume scheme. The scheme is designed to be asymptotic preserving in order to accurately capture both optically thick and thin regimes. We also propose a well-balanced discretization of the radiative flux source term which allows users to capture constant flux steady states with discontinuities in opacity. We use the library Trilinos for linear algebra and the package Kokkos allows us to reach high-performance computing and portability across different architectures, such as multi-core, many-core, and GP-GPU. Results. ARK-RT is able to reproduce standard tests in both free-streaming and diffusive limits, including purely radiative tests and radiation hydrodynamics ones. Using a time-implicit solver is profitable as soon as the time-step given by the hydrodynamics is between 50 and 100 times larger than the explicit time-step for radiative transfer, depending on the preconditioner and the architecture. Nevertheless, more work is needed to ensure stability in all circumstances. Using ARK-RT, we study the propagation of an ionization front in convective dense cores. We show that the ionization front is strongly stable against perturbations even with destabilizing convective motions. As a result, the presence of instabilities should be interpreted with caution. Overall, ARK-RT is well-suited to studying many astrophysical problems involving convection and radiative transfer such as the dynamics of H II regions in massive pre-stellar dense cores and future applications could include planetary atmospheres.
By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth's oceans, fingering convection in ...stellar atmospheres, and moist convection in Earth's atmosphere are derived from the same general diabatic convective instability. We also show that "radiative convection" triggered by the CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analogous to moist and thermohaline convection. We derive a generalization of the mixing-length theory to include the effect of source terms in 1D codes. We show that CO/CH4 "radiative" convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth's atmosphere, thus possibly explaining the reddening in brown dwarf spectra. By using idealized 2D hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism, with other chemical transitions, could be present in many giant and Earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility of using brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.
Aims.
Clouds are expected to form in a broad range of conditions in the atmosphere of exoplanets given the variety of possible condensible species. This diversity, however, might lead to very ...different small-scale dynamics depending on radiative transfer in various thermal conditions. Here, we aim to provide some insight into these dynamical regimes.
Methods.
We performed an analytical linear stability analysis of a compositional discontinuity with a heating source term that depends on a given composition. We also performed idealized two-dimensional simulations of an opacity discontinuity in a stratified medium, using the
ARK
code. We used a two-stream gray model for radiative transfer and explored the brown-dwarf and Earth-like regimes.
Results.
We revealed the existence of a radiative Rayleigh-Taylor instability (RRTI, hereafter, which is a particular case of diabatic Rayleigh-Taylor instability) when an opacity discontinuity is present in a stratified medium. This instability is similar in nature to diabatic convection and relies only on buoyancy with radiative transfer heating and cooling. When the temperature is decreasing with height in the atmosphere, a lower-opacity medium on top of a higher-opacity medium is shown to be dynamically unstable, whereas a higher-opacity medium on top of a lower-opacity medium is stable. This stability-instability behavior is reversed if the temperature is increasing with height.
Conclusions.
The existence of a RRTI could have important implications for the stability of the cloud cover with regard to a wide range of planetary atmospheres. In our Solar System, it could help explain the formation of mammatus cloud in Earth atmospheres and the existence of the Venus cloud deck. Likewise, it suggests that stable and large-scale cloud covers could be ubiquitous in strongly irradiated exoplanets, but might be more patchy in low-irradiated or isolated objects such as brown dwarfs and directly imaged exoplanets.
Aims.
We present a new radiation hydrodynamics code called ARK-RT which uses a two-moment model with the M
1
closure relation for radiative transfer. This code was designed to be ready for ...high-performance computing on exascale architectures.
Methods.
The two-moment model is solved using a finite-volume scheme. The scheme is designed to be asymptotic preserving in order to accurately capture both optically thick and thin regimes. We also propose a well-balanced discretization of the radiative flux source term which allows users to capture constant flux steady states with discontinuities in opacity. We use the library Trilinos for linear algebra and the package Kokkos allows us to reach high-performance computing and portability across different architectures, such as multi-core, many-core, and GP-GPU.
Results.
ARK-RT is able to reproduce standard tests in both free-streaming and diffusive limits, including purely radiative tests and radiation hydrodynamics ones. Using a time-implicit solver is profitable as soon as the time-step given by the hydrodynamics is between 50 and 100 times larger than the explicit time-step for radiative transfer, depending on the preconditioner and the architecture. Nevertheless, more work is needed to ensure stability in all circumstances. Using ARK-RT, we study the propagation of an ionization front in convective dense cores. We show that the ionization front is strongly stable against perturbations even with destabilizing convective motions. As a result, the presence of instabilities should be interpreted with caution. Overall, ARK-RT is well-suited to studying many astrophysical problems involving convection and radiative transfer such as the dynamics of H
II
regions in massive pre-stellar dense cores and future applications could include planetary atmospheres.
Inspired by the statistical mechanics of an ensemble of interacting particles (BBGKY hierarchy), we propose to account for small-scale inhomogeneities in self-gravitating astrophysical fluids by ...deriving a non-ideal Virial theorem and non-ideal Navier-Stokes equations using a decomposition of the gravitational force into a near- and far-field component. These equations involve the pair radial distribution function (similar to the two-point correlation function), similarly to the interaction energy and equation of state in liquids. Small-scale correlations lead to a non-ideal amplification of the gravitational interaction energy, whose omission leads to a missing mass problem, e.g., in galaxies and galaxy clusters. We also propose an extension of the Friedmann equations in the non-ideal regime. We estimate the non-ideal amplification factor of the gravitational interaction energy of the baryons to lie between 5 and 20, potentially explaining the observed value of the Hubble parameter. Within this framework, the acceleration of the expansion emerges naturally because of the increasing number of sub-structures induced by gravitational collapse, which increases their contribution to the total gravitational energy. A simple estimate predicts a non-ideal deceleration parameter qni~-1;this is potentially the first determination of the observed value based on an intuitively physical argument. We suggest that correlations and gravitational interactions could produce a transition to a viscous regime that can lead to flat rotation curves. This transition could also explain the dichotomy between (Keplerian) LSB elliptical galaxy and (non-Keplerian) spiral galaxy rotation profiles. Overall, our results demonstrate that non-ideal effects induced by inhomogeneities must be taken into account in order to properly determine the gravitational dynamics of galaxies and the larger scale universe.
Aims. We present a new radiation hydrodynamics code, called "ARK-RT" which uses a two-moment model with the M1 closure relation for radiative transfer. This code aims at being ready for ...high-performance computing, on exascale architectures. Methods. The two-moment model is solved using a finite volume scheme. The scheme is asymptotic preserving to capture accurately both optically thick and thin regimes. We also propose a well-balanced discretization of the radiative flux source term able to capture constant flux steady states with discontinuities in opacity. We use the library Trilinos for linear algebra and the package Kokkos allows us to reach high-performance computing and portability across different architectures, such as multi-core, many-core, and GP-GPU. Results. ARK-RT is able to reproduce standard tests in both free-streaming and diffusive limits, including purely radiative tests and radiation hydrodynamics ones. Using a time-implicit solver is profitable as soon as the time step given by the hydrodynamics is 50-100 times larger than the explicit time step for radiative transfer, depending on the preconditioner and the architecture. Albeit more work is needed to ensure stability in all circumstances. Using ARK-RT, we study the propagation of an ionization front in convective dense cores. We show that the ionization front is strongly stable against perturbations even with destabilizing convective motions. As a result, the presence of instabilities should be interpreted with caution. Overall, ARK-RT is well-suited to study many astrophysical problems involving convection and radiative transfer such as the dynamics of H ii regions in massive pre-stellar dense cores and future applications could include planetary atmospheres.
Clouds are expected to form in a wide range of conditions in the atmosphere of exoplanets given the large range of possible condensible species. However this diversity might lead to very different ...small-scale dynamics depending on radiative transfer in various thermal conditions: we aim at providing some insights into these dynamical regimes. We perform an analytical linear stability analysis of a compositional discontinuity with a heating source term that depends on composition. We also perform idealized two-dimensional (2D) simulations of an opacity discontinuity in a stratified medium with the code ARK. We use a two-stream grey model for radiative transfer and explore the brown-dwarf and earth-like regimes. We reveal the existence of a Radiative Rayleigh-Taylor Instability (RRTI hereafter, a particular case of diabatic Rayleigh-Taylor instability) when an opacity discontinuity is present in a stratified medium. This instability is similar in nature to diabatic convection and relies only on buoyancy with radiative transfer heating and cooling. When the temperature is decreasing with height in the atmosphere, a lower-opacity medium on top of a higher-opacity medium is dynamically unstable while a higher-opacity medium on top of a lower-opacity medium is stable. This stability/instability behavior is reversed if the temperature is increasing with height. The existence of the RRTI could have important implications for the stability of the cloud cover of a wide range of planetary atmospheres. In our solar system, it could help explain the formation of mammatus cloud in Earth atmospheres and the existence of Venus cloud deck. Likewise, it suggests that stable and large scale cloud covers could be ubiquitous in strongly irradiated exoplanets but might be more patchy in low-irradiated or isolated objects like brown dwarfs and directly imaged exoplanets.
Convection is an important physical process in astrophysics well-studied using numerical simulations under the Boussinesq and/or anelastic approximations. However these approaches reach their limits ...when compressible effects are important in the high Mach flow regime, e.g. in stellar atmospheres or in the presence of accretion shocks. In order to tackle these issues, we propose a new high performance and portable code, called "ARK" with a numerical solver well-suited for the stratified compressible Navier-Stokes equations. We take a finite volume approach with machine precision conservation of mass, transverse momentum and total energy. Based on previous works in applied mathematics we propose the use of a low Mach correction to achieve a good precision in both low and high Mach regimes. The gravity source term is discretized using a well-balanced scheme in order to reach machine precision hydrostatic balance. This new solver is implemented using the Kokkos library in order to achieve high performance computing and portability across different architectures (e.g. multi-core, many-core, and GP-GPU). We show that the low-Mach correction allows to reach the low-Mach regime with a much better accuracy than a standard Godunov-type approach. The combined well-balanced property and the low-Mach correction allowed us to trigger Rayleigh-Bénard convective modes close to the critical Rayleigh number. Furthermore we present 3D turbulent Rayleigh-Bénard convection with low diffusion using the low-Mach correction leading to a higher kinetic energy power spectrum. These results are very promising for future studies of high Mach and highly stratified convective problems in astrophysics.
By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in ...stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analog to moist and thermohaline convection. We derive a generalization of the mixing length theory to include the effect of source terms in 1D codes. We show that CO/CH4 radiative convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth atmosphere thus possibly explaining the reddening in brown-dwarf spectra. By using idealized two-dimensional hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and could be seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism with other chemical transitions could be present in many giant and earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility to use brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.