Thermal Instability-Induced Interstellar Turbulence

We study the dynamics of phase transitions in the interstellar medium by means of three-dimensional hydrodynamic numerical simulations. We use a realistic cooling function and generic nonequilibrium initial conditions to follow the formation history of a multiphase medium in detail in the absence of gravity. We outline a number of qualitatively distinct stages of this process, including a linear isobaric evolution, transition to an isochoric regime, formation of filaments and voids (also known as "thermal" pancakes), the development and decay of supersonic turbulence, an approach to pressure equilibrium, and final relaxation of the multiphase medium. We find that 1%-2% of the initial thermal energy is converted into gas motions in one cooling time. The velocity field then randomizes into turbulence that decays on a dynamical timescale E_k ~ t^-n, 1 < n < 2. While not all initial conditions yield a stable two-phase medium, we examine such a case in detail. We find that the two phases are well mixed with the cold clouds possessing a fine-grained structure near our numerical resolution limit. The amount of gas in the intermediate unstable phase roughly tracks the rms turbulent Mach number, peaking at 25% when M_rms ~ 8, decreasing to 11% when M_rms ~ 0.4.Comment: To appear in the ApJ Letters, April 2002; 5 pages, 3 color figures, mpeg animations available at http://akpc.ucsd.edu/ThermalLetter/thermal.htm

Simulating Supersonic Turbulence in Magnetized Molecular Clouds

We present results of large-scale three-dimensional simulations of weakly magnetized supersonic turbulence at grid resolutions up to 1024^3 cells. Our numerical experiments are carried out with the Piecewise Parabolic Method on a Local Stencil and assume an isothermal equation of state. The turbulence is driven by a large-scale isotropic solenoidal force in a periodic computational domain and fully develops in a few flow crossing times. We then evolve the flow for a number of flow crossing times and analyze various statistical properties of the saturated turbulent state. We show that the energy transfer rate in the inertial range of scales is surprisingly close to a constant, indicating that Kolmogorov's phenomenology for incompressible turbulence can be extended to magnetized supersonic flows. We also discuss numerical dissipation effects and convergence of different turbulence diagnostics as grid resolution refines from 256^3 to 1024^3 cells.Comment: 10 pages, 3 figures, to appear in the proceedings of the DOE/SciDAC 2009 conferenc

Is the Scaling of Supersonic Turbulence Universal?

The statistical properties of turbulence are considered to be universal at sufficiently small length scales, i. e., independent of boundary conditions and large-scale forces acting on the fluid. Analyzing data from numerical simulations of supersonic turbulent flow driven by external forcing, we demonstrate that this is not generally true for the two-point velocity statistics of compressible turbulence. However, a reformulation of the refined similarity hypothesis in terms of the mass-weighted velocity rho^(1/3)v yields scaling laws that are almost insensitive to the forcing. The results imply that the most intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence.Comment: 4 pages, 6 figures, accepted for publication by Phys. Rev. Let

Dissipative Structures in Supersonic Turbulence

We show that density-weighted moments of the dissipation rate, $\epsilon_l$, averaged over a scale $l$, in supersonic turbulence can be successfully explained by the She and L\'ev\^eque model [Phys. Rev. Lett. {\bf 72}, 336 (1994)]. A general method is developed to measure the two parameters of the model, $\gamma$ and $d$, based directly on their physical interpretations as the scaling exponent of the dissipation rate in the most intermittent structures ($\gamma$) and the dimension of the structures ($d$). We find that the best-fit parameters ($\gamma=0.71$ and $d=1.90$) derived from the $\epsilon_l$ scalings in a simulation of supersonic turbulence at Mach 6 agree with their direct measurements, confirming the validity of the model in supersonic turbulence.Comment: 4 pages, 3 figures, accepted by Phys. Rev. Let

Numerical dissipation control in high order shock-capturing schemes for LES of low speed flows

The Yee & Sjögreen adaptive numerical dissipation control in high order scheme (High Order Filter Methods for Wide Range of Compressible Flow Speeds, ICOSAHOM 09, 2009) is further improved for DNS and LES of shock-free turbulence and low speed turbulence with shocklets. There are vastly different requirements in the minimization of numerical dissipation for accurate turbulence simulations of different compressible flow types and flow speeds. Traditionally, the method of choice for shock-free turbulence and low speed turbulence are by spectral, high order central or high order compact schemes with high order linear filters. With a proper control of a local flow sensor, appropriate amount of numerical dissipation in high order shock-capturing schemes can have spectral-like accuracy for compressible low speed turbulent flows. The development of the method includes an adaptive flow sensor with automatic selection on the amount of numerical dissipation needed at each flow location for more accurate DNS and LES simulations with less tuning of parameters for flows with a wide range of flow speed regime during the time-accurate evolution, e.g., time varying random forcing. An automatic selection of the different flow sensors catered to the different flow types is constructed. A Mach curve and high-frequency oscillation indicators are used to reduce the tuning of parameters in controlling the amount of shock-capturing numerical dissipation to be employed for shock-free turbulence, low speed turbulence and turbulence with strong shocks. In Kotov etal. (High Order Numerical Methods for LES of Turbulent Flows with Shocks, ICCFD8, Chengdu, Sichuan, China, July 14–18, 2014) the LES of a turbulent flow with a strong shock by the Yee & Sjögreen scheme indicated a good agreement with the filtered DNS data. A work in progress for the application of the adaptive flow sensor for compressible turbulence with time-varying random forcing is forthcoming. The present study examines the versatility of the Yee & Sjögreen scheme for DNS and LES of traditional low speed flows without forcing. Special attention is focused on the accuracy performance of this scheme using the Smagorinsky and the Germano–Lilly SGS models

Binding Energy and the Fundamental Plane of Globular Clusters

A physical description of the fundamental plane of Galactic globular clusters is developed which explains all empirical trends and correlations in a large number of cluster observables and provides a small but complete set of truly independent constraints on theories of cluster formation and evolution in the Milky Way. Within the theoretical framework of single-mass, isotropic King models, it is shown that (1) 39 regular (non--core-collapsed) globulars with measured core velocity dispersions share a common V-band mass-to-light ratio of 1.45 +/- 0.10, and (2) a complete sample of 109 regular globulars reveals a very strong correlation between cluster binding energy and total luminosity, regulated by Galactocentric position: E_b \propto (L^{2.05} r_{\rm gc}^{-0.4}). The observational scatter about either of these two constraints can be attributed fully to random measurement errors, making them the defining equations of a fundamental plane for globular clusters. A third, weaker correlation, between total luminosity and the King-model concentration parameter, c, is then related to the (non-random) distribution of globulars on the plane. The equations of the FP are used to derive expressions for any cluster observable in terms of only L, r_{\rm gc}, and c. Results are obtained for generic King models and applied specifically to the globular cluster system of the Milky Way.Comment: 60 pages with 19 figures, submitted to Ap

Turbulent Mixing in the Interstellar Medium -- an application for Lagrangian Tracer Particles

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.Comment: 10 pages, 5 figures, conference proceedings "Turbulent Mixing and Beyond 2007

The Statistics of Supersonic Isothermal Turbulence

We present results of large-scale three-dimensional simulations of supersonic Euler turbulence with the piecewise parabolic method and multiple grid resolutions up to 2048^3 points. Our numerical experiments describe non-magnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We discuss numerical resolution issues and demonstrate convergence, in a statistical sense, of the inertial range dynamics in simulations on grids larger than 512^3 points. The simulations allowed us to measure the absolute velocity scaling exponents for the first time. The inertial range velocity scaling in this strongly compressible regime deviates substantially from the incompressible Kolmogorov laws. The slope of the velocity power spectrum, for instance, is -1.95 compared to -5/3 in the incompressible case. The exponent of the third-order velocity structure function is 1.28, while in incompressible turbulence it is known to be unity. We propose a natural extension of Kolmogorov's phenomenology that takes into account compressibility by mixing the velocity and density statistics and preserves the Kolmogorov scaling of the power spectrum and structure functions of the density-weighted velocity v=\rho^{1/3}u. The low-order statistics of v appear to be invariant with respect to changes in the Mach number. For instance, at Mach 6 the slope of the power spectrum of v is -1.69, and the exponent of the third-order structure function of v is unity. We also directly measure the mass dimension of the "fractal" density distribution in the inertial subrange, D_m = 2.4, which is similar to the observed fractal dimension of molecular clouds and agrees well with the cascade phenomenology.Comment: 15 pages, 19 figures, ApJ v665, n2, 200