skip to main content


Title: Cooling and instabilities in colliding flows
ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.  more » « less
Award ID(s):
1813298
NSF-PAR ID:
10356029
Author(s) / Creator(s):
; ; ; ; ; ;
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
  1. 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
  2. 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
  3. 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
  4. 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.

     
    more » « less
  5. 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