Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.

High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternative, we explore the injection of SLRs via dust grains as a way to overcome the mixing barrier. We numerically model the interaction of a supernova remnant containing SLR-rich dust grains with a nearby molecular cloud. The dust grains are subject to drag forces and both thermal and non-thermal sputtering. We confirm that the expanding gas shell stalls upon impact with the dense cloud and that gas-phase SLR injection occurs slowly due to hydrodynamical instabilities at the cloud surface. In contrast, dust grains of sufficient size (> 1 micron) decouple from the gas and penetrate into the cloud within 0.1 Myr. Once inside the cloud, the dust grains are destroyed by sputtering, releasing SLRs and rapidly enriching the dense (potentially star-forming) regions. Our results suggest that SLR transport on dust grains is a viable mechanism to explain SLR enrichment.