- Award ID(s):
- 1813298
- NSF-PAR ID:
- 10356029
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 508
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2266 to 2278
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
null (Ed.)ABSTRACT Zahn’s widely used model for turbulent mixing induced by rotational shear has recently been validated (with some caveats) in non-rotating shear flows. It is not clear, however, whether his model remains valid in the presence of rotation, even though this was its original purpose. Furthermore, new instabilities arise in rotating fluids, such as the Goldreich–Schubert–Fricke (GSF) instability. Which instability dominates when more than one can be excited, and how they influence each other, were open questions that this paper answers. To do so, we use direct numerical simulations of diffusive stratified shear flows in a rotating triply periodic Cartesian domain located at the equator of a star. We find that either the GSF instability or the shear instability tends to take over the other in controlling the system, suggesting that stellar evolution models only need to have a mixing prescription for each individual instability, together with a criterion to determine which one dominates. However, we also find that it is not always easy to predict which instability ‘wins’ for given input parameters, because the diffusive shear instability is subcritical, and only takes place if there is a finite-amplitude turbulence ‘primer’ to seed it. Interestingly, we find that the GSF instability can in some cases play the role of this primer, thereby providing a pathway to excite the subcritical shear instability. This can also drive relaxation oscillations, which may be observable. We conclude by proposing a new model for mixing in the equatorial regions of stellar radiative zones due to differential rotation.more » « less
-
ABSTRACT Without additional heating, radiative cooling of the halo gas of massive galaxies (Milky Way-mass and above) produces cold gas or stars exceeding that observed. Heating from active galactic nucleus (AGN) jets is likely required, but the jet properties remain unclear. This is particularly challenging for galaxy simulations, where the resolution is orders-of-magnitude insufficient to resolve jet formation and evolution. On such scales, the uncertain parameters include the jet energy form [kinetic, thermal, cosmic ray (CR)]; energy, momentum, and mass flux; magnetic fields; opening angle; precession; and duty cycle. We investigate these parameters in a $10^{14}\, {\rm M}_{\odot }$ halo using high-resolution non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios qualitatively meet observational constraints on the halo gas and show that CR-dominated jets most efficiently quench the galaxy by providing CR pressure support and modifying the thermal instability. Mildly relativistic (∼MeV or ∼1010K) thermal plasma jets work but require ∼10 times larger energy input. For fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo mass, kinetic jets are inefficient at quenching unless they have wide opening or precession angles. Magnetic fields also matter less except when the magnetic energy flux reaches ≳ 1044 erg s−1 in a kinetic jet model, which significantly widens the jet cocoon. The criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with a long enough cooling time at the cooling radius.more » « less
-
We investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ( $M_S \sim 7$ ), super-Alfvénic ( $M_A > 1$ ) magnetized flows with frozen-in magnetic flux ( $R_M \gg 1$ ). These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.more » « less
-
ABSTRACT We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma β ≡ Pg/PB to fall, the CGM becomes globally unstable. At the cool disc-to-hot−halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower β haloes, CR heating drives a hot, single-phase diffuse wind with velocities v ∝ (theat/tff)−1, which exceeds the escape velocity when theat/tff ≲ 0.4. In higher β haloes, where the Alfven Mach number is higher, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to ‘trapping’ of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.
-
null (Ed.)Abstract The objective of this work is to investigate the aerodynamics and thermal interactions between a spreading flame and the surrounding walls as well as their effects on fire behaviors. A three-dimensional transient computational fluid dynamics (CFD) combustion model is used to simulate concurrent-flow flame spread over a thin solid sample in a narrow flow duct. The height of the flow duct is the main parameter. The numerical results predict a quenching height for the flow duct below which the flame fails to spread. For duct heights sufficiently larger than the quenching height, the flame reaches a steady spreading state before the sample is fully consumed. The flame spread rate and the pyrolysis length at steady-state first increase and then decrease when the flow duct height decreases. The detailed gas and solid profiles show that flow confinement has multiple effects on the flame spread process. On one hand, it accelerates flow during thermal expansion from combustion, intensifying the flame. On the other hand, increasing flow confinement reduces the oxygen supply to the flame and increases conductive heat loss to the walls, both of which weaken the flame. These competing effects result in the aforementioned nonmonotonic trend of flame spread rate as duct height varies. Near the quenching duct height, the transient model reveals that the flame exhibits oscillation in length, flame temperature, and flame structure. This phenomenon is suspected to be due to thermodiffusive instability.more » « less