Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA multi-shock model for the density variance of anisotropic, highly-magnetised, supersonic turbulence
84
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $sigma^2_{rho/rho_0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parameterises the star formation (SF) rate of a cloud. However, models for $sigma^2_{rho/rho_0}$ all share a common feature -- the variance is assumed to be isotropic. This assumption does not hold when a trans/sub-Alfvenic mean magnetic field, $vec{B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $sigma_{rho/rho_0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $sigma_{rho/rho_0}^2$ in the high Mach number $(mathcal{M})$ limit as small-scale density fluctuations become suppressed by the strong $vec{B}_0$. The model reduces to the isotropic $sigma_{rho/rho_0}^2-mathcal{M}$ relation in the hydrodynamical limit. To validate our model, we calculate $sigma_{rho/rho_0}^2$ from 12~high-resolution, three-dimensional, supersonic, sub-Alfvenic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented over-densities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvenic $vec{B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.
rate research
Read More
It is widely accepted that supersonic, magnetised turbulence plays a fundamental role for star formation in molecular clouds. It produces the initial dense gas seeds out of which new stars can form. However, the exact relation between gas compression, turbulent Mach number, and magnetic field strength is still poorly understood. Here, we introduce and test an analytical prediction for the relation between the density variance and the root-mean-square Mach number in supersonic, isothermal, magnetised turbulent flows. We approximate the density and velocity structure of the interstellar medium as a superposition of shock waves. We obtain the density contrast considering the momentum equation for a single magnetised shock and extrapolate this result to the entire cloud. Depending on the field geometry, we then make three different assumptions based on observational and theoretical constraints: B independent of density, B proportional to the root square of the density and B proportional to the density. We test the analytically derived density variance--Mach number relation with numerical simulations, and find that for B proportional to the root square of the density, the variance in the logarithmic density contrast, $sigma_{ln rho/rho_0}^2=ln[1+b^2mathscr{M}^2beta_0/(beta_0+1)]$, fits very well to simulated data with turbulent forcing parameter b=0.4, when the gas is super-Alfvenic. However, this result breaks down when the turbulence becomes trans-Alfvenic or sub-Alfvenic, because in this regime the turbulence becomes highly anisotropic. Our density variance--Mach number relations simplify to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure $beta_0$ approaches infinite.
We investigate the clustering and dynamics of nano-sized particles (nano-dust) in high-resolution ($1024^3$) simulations of compressible isothermal hydrodynamic turbulence. It is well-established that large grains will decouple from a turbulent gas flow, while small grains will tend to trace the motion of the gas. We demonstrate that nano-sized grains may cluster in a turbulent flow (fractal small-scale clustering), which increases the local grain density by at least a factor of a few. In combination with the fact that nano-dust grains may be abundant in general, and the increased interaction rate due to turbulent motions, aggregation involving nano dust may have a rather high probability. Small-scale clustering will also affect extinction properties. As an example we present an extinction model based on silicates, graphite and metallic iron, assuming strong clustering of grain sizes in the nanometre range, could explain the extreme and rapidly varying ultraviolet extinction in the host of GRB 140506A.
The rich structure that we observe in molecular clouds is due to the interplay between strong magnetic fields and supersonic (turbulent) velocity fluctuations. The velocity fluctuations interact with the magnetic field, causing it too to fluctuate. Using numerical simulations, we explore the nature of such magnetic field fluctuations, $vec{delta B}$, over a wide range of turbulent Mach numbers, $mathcal{M} = 2 - 20$ (i.e., from weak to strong compressibility), and Alfven Mach numbers, $mathcal{M}_{text{A}0} = 0.1 - 100$ (i.e., from strong to weak magnetic mean fields, $B_0$). We derive a compressible quasi-static fluctuation model from the magnetohydrodynamical (MHD) equations and show that velocity gradients parallel to the mean magnetic field give rise to compressible modes in sub-Alfvenic flows, which prevents the flow from becoming two-dimensional, as is the case in incompressible MHD turbulence. We then generalise an analytical model for the magnitude of the magnetic fluctuations to include $mathcal{M}$, and find $|vec{delta B}| = delta B = c_ssqrt{pirho_0}mathcal{M}mathcal{M}_{text{A}0}$, where $c_s$ is the sound speed and $rho_0$ is the mean density of gas. This new relation fits well in the strong $B$-field regime. We go on to study the anisotropy between the perpendicular ($ B_{perp}$) and parallel ($ B_{parallel}$) fluctuations and the mean-normalised fluctuations, which we find follow universal scaling relations, invariant of $mathcal{M}$. We provide a detailed analysis of the morphology for the $delta B_{perp}$ and $delta B_{parallel}$ probability density functions and find that eddies aligned with $B_0$ cause parallel fluctuations that reduce $B_{parallel}$ in the most anisotropic simulations. We discuss broadly the implications of our fluctuation models for magnetised gases in the interstellar medium.
We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary, and we introduce a new diagnostic, the dimensionless star formation efficiency versus density ($epsilon_{rm ff}(s)$) curve, which provides a sensitive diagnostic of the PDF shape and dynamics. We show that the PDF follows a universal functional form consisting of a log-normal at low density with two distinct power law tails at higher density; the first of these represents the onset of self-gravitation, and the second reflects the onset of rotational support. Once the star formation efficiency reaches a few percent, the PDF becomes statistically steady, with no evidence for secular time-evolution at star formation efficiencies from about five to 20 percent. We show that both the Mach number and the virial parameter influence the characteristic densities at which the log-normal gives way to the first power-law, and the first to the second, and we extend (for the former) and develop (for the latter) simple theoretical models for the relationship between these density thresholds and the global properties of the turbulent medium.
We study density fluctuations in supersonic turbulence using both theoretical methods and numerical simulations. A theoretical formulation is developed for the probability distribution function (PDF) of the density at steady state, connecting it to the conditional statistics of the velocity divergence. Two sets of numerical simulations are carried out, using either a Riemann solver to evolve the Euler equations or a finite-difference method to evolve the Navier-Stokes (N-S) equations. After confirming the validity of our theoretical formulation with the N-S simulations, we examine the effects of dynamical processes on the PDF, showing that the nonlinear term in the divergence equation amplifies the right tail of the PDF and reduces the left one, the pressure term reduces both the right and left tails, and the viscosity term, counter-intuitively, broadens the right tail of the PDF. Despite the inaccuracy of the velocity divergence from the Riemann runs, as found in our previous work, we show that the density PDF from the Riemann runs is consistent with that from the N-S runs. Taking advantage of their much higher effective resolution, we then use the Riemann runs to study the dependence of the PDF on the Mach number, $mathcal{M}$, up to $mathcal{M}sim30$. The PDF width, $sigma_{s}$, follows the relation $sigma_{s}^2 = ln (1+b^2 {mathcal M}^2)$, with $bapprox0.38$. However, the PDF exhibits a negative skewness that increases with increasing $mathcal{M}$, so much of the growth of $sigma_{s}$ is accounted for by the growth of the left PDF tail, while the growth of the right tail tends to saturate. Thus, the usual prescription that combines a lognormal shape with the standard variance-Mach number relation greatly overestimates the right PDF tail at large $mathcal{M}$, which may have a significant impact on theoretical models of star formation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
