Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: I. Turbulence Statistics

Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.Comment: 37 pages, 11 figures, 4 table

Effect of freeze-dryer design on heat transfer variability investigated using a 3D mathematical model

International audienceIn the freeze-drying process, vials located at the border of the shelf usually present higher heat flow rates which in turn result in higher product temperatures than central vials. This phenomenon, named edge vial effect, can result in product quality variability within the same batch of vials and between batches at different scales. Our objective was to investigate the effect of various freeze-dryer design features on the heat transfer variability. A 3D mathematical model previously developed in COMSOL Multiphysics and experimentally validated was used to simulate heat transfer of a set of vials located at the edge and in the centre of the shelf. The design features considered were the loading configurations of the vials, the thermal characteristics of the rail, the walls and the shelves and some relevant dimensions of the drying chamber geometry. The presence of the rail in the loading configuration and the value of the shelf emissivity strongly impacted on the heat flow rates received by the vials. Conversely, the heat transfer was not significantly influenced by modifications of the thermal conductivity of the rail, the emissivity of the walls and by the geometry of the drying chamber. The developed model revealed to be a powerful tool to predict the heat transfer variability between edge and central vials for the cycle development and scale-up and to compare various freeze-dryer design features

Gravitational Collapse in Turbulent Molecular Clouds. I. Gasdynamical Turbulence

Observed molecular clouds often appear to have very low star formation efficiencies and lifetimes an order of magnitude longer than their free-fall times. Their support is attributed to the random supersonic motions observed in them. We study the support of molecular clouds against gravitational collapse by supersonic, gas dynamical turbulence using direct numerical simulation. Computations with two different algorithms are compared: a particle-based, Lagrangian method (SPH), and a grid-based, Eulerian, second-order method (ZEUS). The effects of both algorithm and resolution can be studied with this method. We find that, under typical molecular cloud conditions, global collapse can indeed be prevented, but density enhancements caused by strong shocks nevertheless become gravitationally unstable and collapse into dense cores and, presumably, stars. The occurance and efficiency of local collapse decreases as the driving wave length decreases and the driving strength increases. It appears that local collapse can only be prevented entirely with unrealistically short wave length driving, but observed core formation rates can be reproduced with more realistic driving. At high collapse rates, cores are formed on short time scales in coherent structures with high efficiency, while at low collapse rates they are scattered randomly throughout the region and exhibit considerable age spread. We suggest that this naturally explains the observed distinction between isolated and clustered star formation.Comment: Minor revisions in response to referee, thirteen figures, accepted to Astrophys.

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

Nonlinear theory of mirror instability near threshold

An asymptotic model based on a reductive perturbative expansion of the drift kinetic and the Maxwell equations is used to demonstrate that, near the instability threshold, the nonlinear dynamics of mirror modes in a magnetized plasma with anisotropic ion temperatures involves a subcritical bifurcation,leading to the formation of small-scale structures with amplitudes comparable with the ambient magnetic field

Turbulent Control of the Star Formation Efficiency

Supersonic turbulence plays a dual role in molecular clouds: On one hand, it contributes to the global support of the clouds, while on the other it promotes the formation of small-scale density fluctuations, identifiable with clumps and cores. Within these, the local Jeans length \Ljc is reduced, and collapse ensues if \Ljc becomes smaller than the clump size and the magnetic support is insufficient (i.e., the core is ``magnetically supercritical''); otherwise, the clumps do not collapse and are expected to re-expand and disperse on a few free-fall times. This case may correspond to a fraction of the observed starless cores. The star formation efficiency (SFE, the fraction of the cloud's mass that ends up in collapsed objects) is smaller than unity because the mass contained in collapsing clumps is smaller than the total cloud mass. However, in non-magnetic numerical simulations with realistic Mach numbers and turbulence driving scales, the SFE is still larger than observational estimates. The presence of a magnetic field, even if magnetically supercritical, appears to further reduce the SFE, but by reducing the probability of core formation rather than by delaying the collapse of individual cores, as was formerly thought. Precise quantification of these effects as a function of global cloud parameters is still needed.Comment: Invited review for the conference "IMF@50: the Initial Mass Function 50 Years Later", to be published by Kluwer Academic Publishers, eds. E. Corbelli, F. Palla, and H. Zinnecke

Statistics of Dissipation and Enstrophy Induced by a Set of Burgers Vortices

Dissipation and enstropy statistics are calculated for an ensemble of modified Burgers vortices in equilibrium under uniform straining. Different best-fit, finite-range scaling exponents are found for locally-averaged dissipation and enstrophy, in agreement with existing numerical simulations and experiments. However, the ratios of dissipation and enstropy moments supported by axisymmetric vortices of any profile are finite. Therefore the asymptotic scaling exponents for dissipation and enstrophy induced by such vortices are equal in the limit of infinite Reynolds number.Comment: Revtex (4 pages) with 4 postscript figures included via psfi

Turbulent dissipation in the ISM: the coexistence of forced and decaying regimes and implications for galaxy formation and evolution

We discuss the dissipation of turbulent kinetic energy Ek in the global ISM by means of 2-D, MHD, non-isothermal simulations in the presence of model radiative heating and cooling. We argue that dissipation in 2D is representative of that in three dimensions as long as it is dominated by shocks rather than by a turbulent cascade. Energy is injected at a few isolated sites in space, over relatively small scales, and over short time periods. This leads to the coexistence of forced and decaying regimes in the same flow. We find that the ISM-like flow dissipates its turbulent energy rapidly. In simulations with forcing, the input parameters are the radius l_f of the forcing region, the total kinetic energy e_k each source deposits into the flow, and the rate of formation of those regions, sfr_OB. The global dissipation time t_d depends mainly on l_f. In terms of measurable properties of the ISM, t_d >= Sigma_g u_rms^2/(e_k sfr_OB), where Sigma_g is the average gas surface density and u_rms is the rms velocity dispersion. For the solar neighborhood, t_d >= 1.5x10^7 yr. The global dissipation time is consistently smaller than the crossing time of the largest energy-containing scales. In decaying simulations, Ek decreases with time as t^-n, where n~0.8-0.9. This suggests a decay with distance d as Ek\propto d^{-2n/(2-n)} in the mixed forced+decaying case. If applicable to the vertical direction, our results support models of galaxy evolution in which stellar energy injection provides significant support for the gas disk thickness, but not models of galaxy formation in which this energy injection is supposed to reheat an intra-halo medium at distances of up to 10-20 times the optical galaxy size, as the dissipation occurs on distances comparable to the disk height.Comment: 23 pages, including figures. To appear in ApJ. Abstract abridge

Diffusion in supersonic, turbulent, compressible flows

We investigate diffusion in supersonic, turbulent, compressible flows. Supersonic turbulence can be characterized as network of interacting shocks. We consider flows with different rms Mach numbers and where energy necessary to maintain dynamical equilibrium is inserted at different spatial scales. We find that turbulent transport exhibits super-diffusive behavior due to induced bulk motions. In a comoving reference frame, however, diffusion behaves normal and can be described by mixing length theory extended into the supersonic regime.Comment: 11 pages, incl. 5 figures, accepted for publication in Physical Review E (a high-resolution version is available at http://www.aip.de./~ralf/Publications/p21.abstract.html

The Fractal Density Structure in Supersonic Isothermal Turbulence: Solenoidal versus Compressive Energy Injection

In a systematic study, we compare the density statistics in high resolution numerical experiments of supersonic isothermal turbulence, driven by the usually adopted solenoidal (divergence-free) forcing and by compressive (curl-free) forcing. We find that for the same rms Mach number, compressive forcing produces much stronger density enhancements and larger voids compared to solenoidal forcing. Consequently, the Fourier spectra of density fluctuations are significantly steeper. This result is confirmed using the Delta-variance analysis, which yields power-law exponents beta~3.4 for compressive forcing and beta~2.8 for solenoidal forcing. We obtain fractal dimension estimates from the density spectra and Delta-variance scaling, and by using the box counting, mass size and perimeter area methods applied to the volumetric data, projections and slices of our turbulent density fields. Our results suggest that compressive forcing yields fractal dimensions significantly smaller compared to solenoidal forcing. However, the actual values depend sensitively on the adopted method, with the most reliable estimates based on the Delta-variance, or equivalently, on Fourier spectra. Using these methods, we obtain D~2.3 for compressive and D~2.6 for solenoidal forcing, which is within the range of fractal dimension estimates inferred from observations (D~2.0-2.7). The velocity dispersion to size relations for both solenoidal and compressive forcing obtained from velocity spectra follow a power law with exponents in the range 0.4-0.5, in good agreement with previous studies.Comment: 17 pages, 11 figures, ApJ in press, minor changes to language, simulation movies available at http://www.ita.uni-heidelberg.de/~chfeder/videos.shtml?lang=e