Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of radiative turbulent mixing layers
340
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Radiative turbulent mixing layers should be ubiquitous in multi-phase gas with shear flow. They are a potentially attractive explanation for the high ions such as OVI seen in high velocity clouds and the circumgalactic medium (CGM) of galaxies. We perform 3D MHD simulations with non-equilibrium (NEI) and photoionization modeling, with an eye towards testing simple analytic models. Even purely hydrodynamic collisional ionization equilibrium (CIE) calculations have column densities much lower than observations. Characteristic inflow and turbulent velocities are much less than the shear velocity, and the layer width $h propto t_mathrm{cool}^{1/2}$ rather than $h propto t_mathrm{cool}$. Column densities are not independent of density or metallicity as analytic scalings predict, and show surprisingly weak dependence on shear velocity and density contrast. Radiative cooling, rather than Kelvin-Helmholtz instability, appears paramount in determining the saturated state. Low pressure due to fast cooling both seeds turbulence and sets the entrainment rate of hot gas, whose enthalpy flux, along with turbulent dissipation, energizes the layer. Regardless of initial geometry, magnetic fields are amplified and stabilize the mixing layer via magnetic tension, producing almost laminar flow and depressing column densities. NEI effects can boost column densities by factors of a few. Suppression of cooling by NEI or photoionization can in principle also increase OVI column densities, but in practice is unimportant for CGM conditions. To explain observations, sightlines must pierce hundreds or thousands of mixing layers, which may be plausible if the CGM exists as a `fog of tiny cloudlets.
rate research
Read More
We present models of turbulent mixing at the boundaries between hot (T~10^{6-7} K) and warm material (T~10^4 K) in the interstellar medium, using a three-dimensional magnetohydrodynamical code, with radiative cooling. The source of turbulence in our simulations is a Kelvin-Helmholtz instability, produced by shear between the two media. We found, that because the growth rate of the large scale modes in the instability is rather slow, it takes a significant amount of time (~1 Myr) for turbulence to produce effective mixing. We find that the total column densities of the highly ionized species (C IV, N V, and O VI) per interface (assuming ionization equilibrium) are similar to previous steady-state non-equilibrium ionization models, but grow slowly from log N ~10^{11} to a few 10^{12} cm^{-2} as the interface evolves. However, the column density ratios can differ significantly from previous estimates, with an order of magnitude variation in N(C IV)/N(O VI) as the mixing develops.
The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relationship between large-scale gas kinematics and magnetic fields through non-radiative magnetohydrodynamical simulations of the creation, evolution and disruption of AGN jet-inflated lobes in an isolated Perseus-like galaxy cluster, with and without pre-existing turbulence. In particular, we connect cluster velocity measurements with mock RM maps to highlight their underlying physical connection, which opens up the possibility of comparing turbulence levels in two different observables. For single jet outbursts, we find only a local impact on the velocity field, i.e. the associated increase in velocity dispersion is not volume-filling. Furthermore, in a setup with pre-existing turbulence, this increase in velocity dispersion is largely hidden. We use mock X-ray observations to show that at arcmin resolution, the velocity dispersion is therefore dominated by existing large-scale turbulence and is only minimally altered by the presence of a jet. For the velocity structure of central gas uplifted by buoyantly rising lobes, we find fast, coherent outflows with low velocity dispersion. Our results highlight that projected velocity distributions show complex structures which pose challenges for the interpretation of observations.
A crucial aspect of 3D Monte Carlo radiative transfer is the choice of the spatial grid used to partition the dusty medium. We critically investigate the use of octree grids in Monte Carlo dust radiative transfer, with two different octree construction algorithms (regular and barycentric subdivision) and three different octree traversal algorithms (top-down, neighbour list, and the bookkeeping method). In general, regular octree grids need higher levels of subdivision compared to the barycentric grids for a fixed maximum cell mass threshold criterion. The total number of grid cells, however, depends on the geometry of the model. Surprisingly, regular octree grid simulations turn out to be 10 to 20% more efficient in run time than the barycentric grid simulations, even for those cases where the latter contain fewer grid cells than the former. Furthermore, we find that storing neighbour lists for each cell in an octree, ordered according to decreasing overlap area, is worth the additional memory and implementation overhead: using neighbour lists can cut down the grid traversal by 20% compared to the traditional top-down method. In conclusion, the combination of a regular node subdivision and the neighbour list method results in the most efficient octree structure for Monte Carlo radiative transfer simulations.
We present results from multifrequency radiative hydrodynamical chemistry simulations addressing primordial star formation and related stellar feedback from various populations of stars, stellar energy distributions (SEDs) and initial mass functions. Spectra for massive stars, intermediate-mass stars and regular solar-like stars are adopted over a grid of 150 frequency bins and consistently coupled with hydrodynamics, heavy-element pollution and non-equilibrium species calculations. Powerful massive population III stars are found to be able to largely ionize H and, subsequently, He and He$^+$, causing an inversion of the equation of state and a boost of the Jeans masses in the early intergalactic medium. Radiative effects on star formation rates are between a factor of a few and 1 dex, depending on the SED. Radiative processes are responsible for gas heating and photoevaporation, although emission from soft SEDs has minor impacts. These findings have implications for cosmic gas preheating, primordial direct-collapse black holes, the build-up of cosmic fossils such as low-mass dwarf galaxies, the role of AGNi during reionization, the early formation of extended disks and angular-momentum catastrophe.
Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in neon novae, yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
