We study the probability distribution function (PDF) of the mass density in simulations of supersonic turbulence with properties appropriate for molecular clouds. For this study we use Athena, a new higher-order Godunov code. We find there are surprisingly similar relationships between the mean of the time-averaged PDF and the turbulent Mach number for driven hydrodynamic and strong-field MHD turbulence. There is, however, a large scatter about these relations, indicating a high level of temporal and spatial variability in the PDF. Thus, the PDF of the mass density is unlikely to be a good measure of magnetic field strength. We also find the PDF of decaying MHD turbulence deviates from the mean-Mach relation found in the driven case. This implies that the instantaneous Mach number alone is not enough to determine the statistical properties of turbulence that is out of equilibrium. The scatter about the mean-Mach relation for driven turbulence, along with the large departure of decaying turbulence PDFs from those of driven turbulence, may illuminate one factor contributing to the large observed cloud-to-cloud variation in the star formation rate per solar mass.

We study energy dissipation and heating by supersonic MHD turbulence in molecular clouds using Athena, a new higher-order Godunov code. We analyze the dependence of the saturation amplitude, energy dissipation characteristics, power spectra, sonic scaling, and indicators of intermittency in the turbulence on factors such as the magnetic field strength, driving scale, energy injection rate, and numerical resolution. While convergence in the energies is reached at moderate resolutions, we find that the power spectra require much higher resolutions that are difficult to obtain. In a 1024^3 hydro run, we find a power law relationship between the velocity dispersion and the spatial scale on which it is measured, while for an MHD run at the same resolution we find no such power law. The time-variability and temperature intermittency in the turbulence both show a dependence on the driving scale, indicating that numerically driving turbulence by an arbitrary mechanism may not allow a realistic representation of these properties. We also note similar features in the power spectrum of the compressive component of velocity for supersonic MHD turbulence as in the velocity spectrum of an initially-spherical MHD blast wave, implying that the power law form does not rule out shocks, rather than a turbulent cascade, playing a significant role in the regulation of energy transfer between spatial scales.

The electron mean free path in the intracluster medium (ICM) of galaxy clusters is much larger than the gyroradius; thus, heat is transported anisotropically along magnetic field lines. We show that the intracluster medium is unstable to the magnetothermal instability (MTI) using MHD simulations with anisotropic thermal conduction. As a result of the MTI, we find that the temperature profile of the ICM can be substantially modified on timescales of several billion years while the magnetic field is amplified by dynamo action up to more than fifty times the original energy. We also show that the instability drives field lines to become preferentially radial leading to conduction that is a highly efficient fraction of the Spitzer conductivity. As such, we present the first self-consistent calculation of the effective thermal conductivity in the ICM.

Using fully three-dimensional hydrodynamic simulations, we investigate the effect of the Coriolis force on the hydrodynamic and observable properties of colliding wind binary systems. To make the calculations tractable, we assume adiabatic, constant velocity winds. The neglect of radiative driving, gravitational deceleration, and cooling limit the application of our models to real systems. However, these assumptions allow us to isolate the effect of the Coriolis force, and by simplifying the calculations, allow us to use a higher resolution (up to 640^3) and to conduct a larger survey of parameter space. We study the dynamics of collidng winds with equal mass loss rates and velocities emanating from equal mass stars on circular orbits, with a range of values for the ratio of the wind to orbital velocity. We also study the dynamics of winds from stars on elliptical orbits and with unequal strength winds. Orbital motion of the stars sweeps the shocked wind gas into an Archimedean spiral, with asymmetric shock strengths and therefore unequal postshock temperatures and densities in the leading and trailing edges of the spiral. We observe the Kelvin-Helmholtz instability at the contact surface between the shocked winds in systems with orbital motion even when the winds are identical. The change in shock strengths caused by orbital motion increases the volume of X-ray emitting post-shock gas with T > 0.59 keV by 63% for a typical system as the ratio of wind velocity to orbital velocity decreases to V_w/V_o = 2.5. This causes increased free-free emission from systems with shorter orbital periods and an altered time-dependence of the wind attenuation. We comment on the importance of the effects of orbital motion on the observable properties of colliding wind binaries.