Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDetailed Balance and Exact Results for Density Fluctuations in Supersonic Turbulence
84
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The probabilistic approach to turbulence is applied to investigate density fluctuations in supersonic turbulence. We derive kinetic equations for the probability distribution function (PDF) of the logarithm of the density field, $s$, in compressible turbulence in two forms: a first-order partial differential equation involving the average divergence conditioned on the flow density, $langle abla cdot {bs u} | srangle$, and a Fokker-Planck equation with the drift and diffusion coefficients equal to $-langle {bs u} cdot abla s | srangle$ and $langle {bs u} cdot abla s | srangle$, respectively. Assuming statistical homogeneity only, the detailed balance at steady state leads to two exact results, $langle abla cdot {bs u} | s rangle =0$, and $langle {bs u} cdot abla s | srangle=0$. The former indicates a balance of the flow divergence over all expanding and contracting regions at each given density. The exact results provide an objective criterion to judge the accuracy of numerical codes with respect to the density statistics in supersonic turbulence. We also present a method to estimate the effective numerical diffusion as a function of the flow density and discuss its effects on the shape of the density PDF.
rate research
Read More
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.
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.
Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.
This paper has been withdrawn
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
