Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAlgorithmic comparisons of decaying, isothermal, supersonic turbulence
131
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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)
rate research
Read More
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.
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.
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.
In this study we investigate the formation and properties of prestellar and protostellar cores using hydrodynamic, self-gravitating Adaptive Mesh Refinement simulations, comparing the cases where turbulence is continually driven and where it is allowed to decay. We model observations of these cores in the C$^{18}$O$(2to 1)$, NH$_3(1,1)$, and N$_2$H$^+(1to 0)$ lines, and from the simulated observations we measure the linewidths of individual cores, the linewidths of the surrounding gas, and the motions of the cores relative to one another. Some of these distributions are significantly different in the driven and decaying runs, making them potential diagnostics for determining whether the turbulence in observed star-forming clouds is driven or decaying. Comparing our simulations with observed cores in the Perseus and $rho$ Ophiuchus clouds shows reasonably good agreement between the observed and simulated core-to-core velocity dispersions for both the driven and decaying cases. However, we find that the linewidths through protostellar cores in both simulations are too large compared to the observations. The disagreement is noticably worse for the decaying simulation, in which cores show highly supersonic infall signatures in their centers that decrease toward their edges, a pattern not seen in the observed regions. This result gives some support to the use of driven turbulence for modeling regions of star formation, but reaching a firm conclusion on the relative merits of driven or decaying turbulence will require more complete data on a larger sample of clouds as well as simulations that include magnetic fields, outflows, and thermal feedback from the protostars.
We use unpublished and published VLBI results to investigate the geometry and the statistical properties of the velocity field traced by H2O masers in five galactic regions of star formation -- Sgr B2(M), W49N, W51(MAIN), W51N, and W3(OH). In all sources the angular distribution of the H2O hot spots demonstrates approximate self-similarity (fractality) over almost four orders of magnitude in scale, with the calculated fractal dimension d between (approximately) 0.2 and 1.0. In all sources, the lower order structure functions for the line-of-sight component of the velocity field are satisfactorily approximated by power laws, with the exponents near their classic Kolmogorov values for the high-Reynolds-number incompressible turbulence. These two facts, as well as the observed significant excess of large deviations of the two-point velocity increments from their mean values, strongly suggest that the H2O masers in regions of star formation trace turbulence. We propose a new conceptual model of these masers in which maser hot spots originate at the sites of ultimate dissipation of highly supersonic turbulence produced in the ambient gas by the intensive gas outflow from a newly-born star. Due to the high brightness and small angular sizes of masing hot spots and the possibility of measuring their positions and velocities with high precision, they become a unique probe of supersonic turbulence.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
