Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDense Regions in Supersonic Isothermal Turbulence
61
0
0.0
(
0
)
Authors
Brant Robertson
Ask ChatGPT about the research
No Arabic abstract
The properties of supersonic isothermal turbulence influence a variety of astrophysical phenomena, including the structure and evolution of star forming clouds. This work presents a simple model for the structure of dense regions in turbulence in which the density distribution behind isothermal shocks originates from rough hydrostatic balance between the pressure gradient behind the shock and its deceleration from ram pressure applied by the background fluid. Using simulations of supersonic isothermal turbulence and idealized waves moving through a background medium, we show that the structural properties of dense, shocked regions broadly agree with our analytical model. Our work provides a new conceptual picture for describing the dense regions, which complements theoretical efforts to understand the bulk statistical properties of turbulence and attempts to model the more complex features of star forming clouds like magnetic fields, self-gravity, or radiative properties.
rate research
Read More
Supersonic isothermal turbulence establishes a network of transient dense shocks that sweep up material and have a density profile described by balance between ram pressure of the background fluid versus the magnetic and gas pressure gradient behind the shock. These rare, densest regions of a turbulent environment can become Jeans unstable and collapse to form pre-stellar cores. Using numerical simulations of magneto-gravo-turbulence, we describe the structural properties of dense shocks, which are the seeds of gravitational collapse, as a function of magnetic field strength. In the regime of a weak magnetic field, the collapse is isotropic. Strong magnetic field strengths lead to significant anisotropy in the shocked distribution and collapse occurs preferentially parallel to the field lines. Our work provides insight into analyzing the magnetic field topology and density structures of young protostellar collapse, which the theory presented here predicts are associated with large-scale strong shocks that persist for at least a free-fall time.
The tight correlation between turbulence and luminosity in Giant HII Regions is not well understood. While the luminosity is due to the UV radiation from the massive stars in the ionizing clusters, it is not clear what powers the turbulence. Observations of the two prototypical Giant HII Regions in the local Universe, 30 Doradus and NGC604, show that part of the kinetic energy of the nebular gas comes from the combined stellar winds of the most massive stars - the cluster winds, but not all. We present a study of the kinematics of 30 Doradus based on archival VLT FLAMES/GIRAFFE data and new high resolution observations with HARPS. We find that the nebular structure and kinematics are shaped by a hot cluster wind and not by the stellar winds of individual stars. The cluster wind powers most of the turbulence of the nebular gas, with a small but significant contribution from the combined gravitational potential of stars and gas. We estimate the total mass of 30 Doradus and we argue that the region does not contain significant amounts of neutral (HI) gas, and that the giant molecular cloud 30Dor-10 that is close to the center of the nebula in projection is in fact an inflating cloud tens of parsecs away from R136, the core of the ionizing cluster. We rule out a Kolmogorov-like turbulent kinetic energy cascade as the source of supersonic turbulence in Giant HII Regions.
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.
Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
