skip to main content


Title: Finite shock model of density in supersonic turbulence
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
Award ID(s):
1616026
NSF-PAR ID:
10400978
Author(s) / Creator(s):
;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society: Letters
Volume:
521
Issue:
1
ISSN:
1745-3925
Format(s):
Medium: X Size: p. L64-L69
Size(s):
["p. L64-L69"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. ABSTRACT

    Owing to the complexity of turbulent magnetic fields, modelling the diffusion of cosmic rays is challenging. Based on the current understanding of anisotropic magnetohydrodynamic (MHD) turbulence, we use test particles to examine the cosmic rays’ superdiffusion in the direction perpendicular to the mean magnetic field. By changing Alfvén Mach number MA and sonic Mach number MS of compressible MHD simulations, our study covers a wide range of astrophysical conditions including subsonic warm gas phase and supersonic cold molecular gas. We show that freely streaming cosmic rays’ perpendicular displacement increases as 3/2 to the power of the time travelled along local magnetic field lines. This power-law index changes to 3/4 if the parallel propagation is diffusive. We find that the cosmic rays’ parallel mean free path decreases in a power-law relation of $M_\mathrm{ A}^{-2}$ in supersonic turbulence. We investigate the energy fraction of slow, fast, and Alfvénic modes and confirm the dominance of Alfvénic modes in the perpendicular superdiffusion. In particular, the energy fraction of fast mode, which is the main agent for pitch-angle scattering, increases with MA, but is insensitive to MS ≥ 2. Accordingly, our results suggest that the suppressed diffusion in supersonic molecular clouds arises primarily due to the variations of MA instead of MS.

     
    more » « less
  3. Abstract

    We study the impact of compressibility on two-dimensional turbulent flows, such as those modeling astrophysical disks. We demonstrate that the direction of cascade undergoes continuous transition as the Mach numberMaincreases, from inverse atMa = 0, to direct atMa=. Thus, atMa1comparable amounts of energy flow from the pumping scale to large and small scales, in accord with previous data. For supersonic turbulence withMa1the cascade is direct, as in three dimensions, which results in multifractal density field. For that regime (Ma1) we derive a Kolmogorov-type law for potential forcing and obtain an explicit expression for the third order correlation tensor of the velocity. We further show that all third order structure functions are zero up to first order in the inertial range scales, which is in sharp contrast with incompressible turbulence where the third order structure function, that describes the energy flux associated with the energy cascade is non-zero. The properties of compressible turbulence have significant implications on the amplification of magnetic fields in conducting fluids. We thus demonstrate that imposing external magnetic field on compressible flows of conducting fluids allows to manipulate the flow producing possibly large changes even at small Mach numbers. Thus Zeldovich’s antidynamo theorem, by which atMa = 0 the magnetic field is zero in the steady state, must be used with caution. Real flows have finiteMaand, however small it is, for large enough values ofI, the magnetic flux through the disk, the magnetic field changes the flow appreciably, or rearranges it completely. This renders the limitMa → 0 singular for non-zero values ofI. Of particular interest is the effect of the density multifractality, atMa1which is relevant for astrophysical disks. We demonstrate that in that regime, in the presence of non-zeroIthe magnetic field energy is enhanced by a large factor as compared to its estimates based on the mean field. Finally, based on the insights described above, we propose a novel two-dimensional Burgers’ turbulence, whose three-dimensional counterpart is used for studies of the large-scale structure of the Universe, as a model for supersonic two-dimensional magnetohydrodynamic flows.

     
    more » « less
  4. Abstract Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities, and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits an increased average sonic Mach number, and maintains a slower star formation rate (SFR) over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power law, and then to the stars, becomes significantly slower, resulting in slower SFRs and longer depletion times. 
    more » « less
  5. 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