We investigate the nature of the Alfv\'enic turbulence cascade in two fluid MHD simulations in order to determine if turbulence is damped once the ion and neutral species become decoupled at a critical scale called the ambipolar diffusion scale (L$_{AD}$). Using mode decomposition to separate the three classical MHD modes, we study the second order structure functions of the Alfv\'en mode velocity field of both neutrals and ions in the reference frame of the local magnetic field. On scales greater than L$_{AD}$ we confirm that two fluid turbulence strongly resembles single fluid MHD turbulence. Our simulations show that the behavior of two fluid turbulence becomes more complex on scales less than L$_{AD}$. We find that Alfvenic turbulence can exist past L$_{AD}$ when the turbulence is globally super-Alfv\'enic, with the ions and neutrals forming separate cascades once decoupling has taken place. When turbulence is globally sub-Alfvenic and hence strongly anisotropic with a large separation between the parallel and perpendicular decoupling scales, turbulence is damped at L$_{AD}$. We also find that the power spectrum of the kinetic energy in the damped regime is consistent with a $k^{-4}$ scaling (in agreement with the predictions of Lazarian, Vishniac & Cho 2004).

We present high resolution simulations of two-fluid (ion-neutral) MHD turbulence with resolutions as large as 512^3. The simulations are supersonic and mildly sub-Alfvenic, in keeping with the conditions present in molecular clouds. Such turbulence is thought to influence star formation processes in molecular clouds because typical cores form on length scales that are comparable to the dissipation scales of this turbulence in the ions. The simulations are motivated by the fact that recent studies of isophotologue lines in molecular clouds have found significant differences in the linewidth-size relationship for neutral and ion species. The goals of this paper are to explain those observations using simulations and analytic theory, present a new set of density-based diagnostics by drawing on similar diagnostics that have been obtained by studying single-fluid turbulence, and show that our two-fluid simulations play a vital role in reconciling alternative models of star formation. The velocity-dependent diagnostics display a complementarity with the density-dependent diagnostics. We find that the linewidth-size relationships show a prominent difference between ions and neutrals when the line of sight is orthogonal to the mean field. This is because the MHD waves in the ions differ from hydrodynamic waves in the neutrals. We also find that the density probability distribution functions (PDFs) show prominent differences between the ions and neutrals when the line of sight is parallel to the mean field. This is because the velocity fluctuations in the ions tend to produce column density fluctuations along field lines. When the magnetic field makes an angle to the line of sight, both observable differences should be visible. These diagnostics should be easy for observers to test. This analysis assumes optically thin lines and a mean magnetic field that is uniform in direction within a cloud.

This report describes a new magnetohydrodynamic numerical model based on a hexagonal spherical geodesic grid. The model is designed to simulate astrophysical flows of partially ionized plasmas around a central compact object, such as a star or a planet with a magnetic field. The geodesic grid, produced by a recursive subdivision of a base platonic solid (an icosahedron), is free from control volume singularities inherent in spherical polar grids. Multiple populations of plasma and neutral particles, coupled via charge-exchange interactions, can be simulated simultaneously with this model. Our numerical scheme uses piecewise linear reconstruction on a surface of a sphere in a local two-dimensional "Cartesian" frame. The code employs HLL-type approximate Riemann solvers and includes facilities to control the divergence of magnetic field and maintain pressure positivity. Several test solutions are discussed, including a problem of an interaction between the solar wind and the local interstellar medium, and a simulation of Earth's magnetosphere.

We have studied saturated, MRI-driven turbulence using three-dimensional, isothermal simulations with resolutions that extend from 64 to 192 zones in each direction. The simulations were performed with several higher order Godunov algorithms. A variety of reconstruction strategies as well as a variety of Riemann solvers are tried. We show that the details of the isothermal MRI-driven turbulence depend principally on the Riemann solver and secondarily on the reconstruction strategy. Furthermore, we find that the effective viscosity parameter parameter tends to show progressively smaller decrements with increasing resolution when the best reconstruction strategy (WENO) and the best Riemann solver (linearized)are used. We attribute this result to the more sophisticated dissipation mechanisms that are used in higher-order Godunov schemes. Spectral analysis and transfer functions have been used to quantify the dissipative processes in these higher-order Godunov schemes.