skip to main content


Title: Shock–multicloud interactions in galactic outflows – I. Cloud layers with lognormal density distributions
ABSTRACT We report three-dimensional hydrodynamical simulations of shocks (${\cal M_{\rm shock}}\ge 4$) interacting with fractal multicloud layers. The evolution of shock–multicloud systems consists of four stages: a shock-splitting phase in which reflected and refracted shocks are generated, a compression phase in which the forward shock compresses cloud material, an expansion phase triggered by internal heating and shock re-acceleration, and a mixing phase in which shear instabilities generate turbulence. We compare multicloud layers with narrow ($\sigma _{\rho }=1.9\bar{\rho }$) and wide ($\sigma _{\rho }=5.9\bar{\rho }$) lognormal density distributions characteristic of Mach ≈ 5 supersonic turbulence driven by solenoidal and compressive modes. Our simulations show that outflowing cloud material contains imprints of the density structure of their native environments. The dynamics and disruption of multicloud systems depend on the porosity and the number of cloudlets in the layers. ‘Solenoidal’ layers mix less, generate less turbulence, accelerate faster, and form a more coherent mixed-gas shell than the more porous ‘compressive’ layers. Similarly, multicloud systems with more cloudlets quench mixing via a shielding effect and enhance momentum transfer. Mass loading of diffuse mixed gas is efficient in all models, but direct dense gas entrainment is highly inefficient. Dense gas only survives in compressive clouds, but has low speeds. If normalized with respect to the shock-passage time, the evolution shows invariance for shock Mach numbers ≥10 and different cloud-generating seeds, and slightly weaker scaling for lower Mach numbers and thinner cloud layers. Multicloud systems also have better convergence properties than single-cloud systems, with a resolution of eight cells per cloud radius being sufficient to capture their overall dynamics.  more » « less
Award ID(s):
1715876
NSF-PAR ID:
10244183
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
499
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2173 to 2195
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($\sigma _{\rho /\rho _0}$) to the turbulent sonic Mach number ($\mathcal {M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterized by the turbulence driving parameter: $b=\sigma _{\rho /\rho _0}/\mathcal {M}$. The turbulence driving parameter ranges from b = 1/3 (purely solenoidal) to b = 1 (purely compressive), with b = 0.38 characterizing the natural mixture (1/3 compressive, 2/3 solenoidal) of the two driving modes. Here, we present a new method for recovering $\sigma _{\rho /\rho _0}$, $\mathcal {M}$, and b, from observations on galactic scales, using a roving kernel to produce maps of these quantities from column density and centroid velocity maps. We apply our method to high-resolution ${\rm H}\,\rm{\small I}$ emission observations of the Small Magellanic Cloud (SMC) from the GASKAP-HI survey. We find that the turbulence driving parameter varies between b ∼ 0.3 and 1.0 within the main body of the SMC, but the median value converges to b ∼ 0.51, suggesting that the turbulence is overall driven more compressively (b > 0.38). We observe no correlation between the b parameter and ${\rm H}\,\rm{\small I}$ or H α intensity, indicating that compressive driving of ${\rm H}\,\rm{\small I}$ turbulence cannot be determined solely by observing ${\rm H}\,\rm{\small I}$ or H α emission density, and that velocity information must also be considered. Further investigation is required to link our findings to potential driving mechanisms such as star-formation feedback, gravitational collapse, or cloud–cloud collisions.

     
    more » « less
  2. ABSTRACT

    In star-forming clouds, high velocity flow gives rise to large fluctuations of density. In this work, we explore the correlation between velocity magnitude (speed) and density. We develop an analytic formula for the joint probability distribution function (PDF) of density and speed, and discuss its properties. In order to develop an accurate model for the joint PDF, we first develop improved models of the marginalized distributions of density and speed. We confront our results with a suite of 12 supersonic isothermal simulations with resolution of $1024^3$ cells in which the turbulence is driven by 3 different forcing modes (solenoidal, mixed, and compressive) and 4 rms Mach numbers (1, 2, 4, 8). We show, that for transsonic turbulence, density and speed are correlated to a considerable degree and the simple assumption of independence fails to accurately describe their statistics. In the supersonic regime, the correlations tend to weaken with growing Mach number. Our new model of the joint and marginalized PDFs are a factor of 3 better than uncorrelated, and provides insight into this important process.

     
    more » « less
  3. Abstract Stellar winds contain enough energy to easily disrupt the parent cloud surrounding a nascent star cluster, and for this reason they have long been considered candidates for regulating star formation. However, direct observations suggest most wind power is lost, and Lancaster et al. recently proposed that this is due to efficient mixing and cooling processes. Here we simulate star formation with wind feedback in turbulent, self-gravitating clouds, extending our previous work. Our simulations cover clouds with an initial surface density of 10 2 –10 4 M ⊙ pc −2 and show that star formation and residual gas dispersal are complete within two to eight initial cloud freefall times. The “efficiently cooled” model for stellar wind bubble evolution predicts that enough energy is lost for the bubbles to become momentum-driven; we find that this is satisfied in our simulations. We also find that wind energy losses from turbulent, radiative mixing layers dominate losses by “cloud leakage” over the timescales relevant for star formation. We show that the net star formation efficiency (SFE) in our simulations can be explained by theories that apply wind momentum to disperse cloud gas, allowing for highly inhomogeneous internal cloud structure. For very dense clouds, the SFE is similar to those observed in extreme star-forming environments. Finally, we find that, while self-pollution by wind material is insignificant in cloud conditions with moderate density (only ≲10 −4 of the stellar mass originated in winds), our simulations with conditions more typical of a super star cluster have star particles that form with as much as 1% of their mass in wind material. 
    more » « less
  4. ABSTRACT

    The probability distribution of density in isothermal, supersonic, turbulent gas is approximately lognormal. This behaviour can be traced back to the shock waves travelling through the medium, which randomly adjust the density by a random factor of the local sonic Mach number squared. Provided a certain parcel of gas experiences a large number of shocks, due to the central limit theorem, the resulting distribution for density is lognormal. We explore a model in which parcels of gas undergo finite number of shocks before relaxing to the ambient density, causing the distribution for density to deviate from a lognormal. We confront this model with numerical simulations with various rms Mach numbers ranging from subsonic as low as 0.1 to supersonic at 25. We find that the fits to the finite formula are an order of magnitude better than a lognormal. The model naturally extends even to subsonic flows, where no shocks exist.

     
    more » « less
  5. Abstract

    Foreshock transients such as foreshock bubbles (FBs), hot flow anomalies (HFAs), and spontaneous hot flow anomalies (SHFAs) display heated, tenuous cores and large flow deflections bounded by compressional boundaries. THEMIS and Cluster observations show that some cores contain local density enhancements which can be studied to better understand the evolution processes of foreshock transients. However, closer examinations of these substructures were not feasible until the availability of the higher resolution data from the Magnetospheric Multiscale mission (MMS). We identify 164 FB‐like, HFA‐like, and SHFA events from two MMS dayside phases for a statistical study to investigate their solar wind conditions, properties, and substructure properties. Occurrence rates of the three event types are higher for lower magnetic field strengths, higher solar wind speeds and Mach numbers, and quasi‐parallel bow shocks. Events usually span up to 3REalong the bow shock surface and extend up to 6REupstream from the bow shock. Though events with and without substructures exhibit similar solar wind conditions, events with substructures are more likely to have longer core durations and larger sizes. Substructure densities display a positive correlation with bulk flows and a negative correlation with temperatures. Substructure sizes vary between 4 and 24 ion inertial lengths, indicating multiple formation mechanisms. Substructures could be the boundary between two foreshock transient events that have merged into a single event, fast‐mode variations, generated by slow or mirror mode instabilities, or produced from instabilities due to parameter gradients at the compressional boundaries or shocks.

     
    more » « less