Magnetic collimation of protostellar winds into bipolar outflows

Researchers describe self-consistent 2-D magnetohydrodynamic (MHD) simulations of the collimation of an isotropic protostellar wind into bipolar outflows by magnetic stresses in the ambient medium. A variety of ambient field strengths, wind luminosities, and density profiles were studied. Collimation occurs when the energy of the magnetic field swept up by the expanding bubble approaches the bubble thermal energy. Measured axial and radial expansion rates are in good agreement with the analytical predictions of Konigl (1982)

Cosmological Radiation Hydrodynamics with ENZO

We describe an extension of the cosmological hydrodynamics code ENZO to include the self-consistent transport of ionizing radiation modeled in the flux-limited diffusion approximation. A novel feature of our algorithm is a coupled implicit solution of radiation transport, ionization kinetics, and gas photoheating, making the timestepping for this portion of the calculation resolution independent. The implicit system is coupled to the explicit cosmological hydrodynamics through operator splitting and solved with scalable multigrid methods. We summarize the numerical method, present a verification test on cosmological Stromgren spheres, and then apply it to the problem of cosmological hydrogen reionization.Comment: 14 pages, 3 figures, to appear in Recent Directions in Astrophysical Quantitative Spectroscopy and Radiation Hydrodynamics, Ed. I. Hubeny, American Institute of Physics (2009

Simulating Radiating and Magnetized Flows in Multi-Dimensions with ZEUS-MP

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.Comment: Accepted for publication in the ApJ Supplement. 42 pages with 29 inlined figures; uses emulateapj.sty. Discussions in sections 2 - 4 improved per referee comments; several figures modified to illustrate grid resolution. ZEUS-MP source code and documentation available from the Laboratory for Computational Astrophysics at http://lca.ucsd.edu/codes/currentcodes/zeusmp2

Proceedings of the 24th annual Central Plains irrigation conference

Presented at Proceedings of the 24th annual Central Plains irrigation conference held on February 21-22 in Colby, Kansas.Includes bibliographical references

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

Pre-peak ram pressure stripping in the Virgo cluster spiral galaxy NGC 4501

VIVA HI observations of the Virgo spiral galaxy NGC 4501 are presented. The HI disk is sharply truncated to the southwest, well within the stellar disk. A region of low surface-density gas, which is more extended than the main HI disk, is discovered northeast of the galaxy center. These data are compared to existing 6cm polarized radio continuum emission, Halpha, and optical broad band images. We observe a coincidence between the western HI and polarized emission edges, on the one hand, and a faint Halpha emission ridge, on the other. The polarized emission maxima are located within the gaps between the spiral arms and the faint Halpha ridge. Based on the comparison of these observations with a sample of dynamical simulations with different values for maximum ram pressure and different inclination angles between the disk and the orbital plane,we conclude that ram pressure stripping can account for the main observed characteristics. NGC 4501 is stripped nearly edge-on, is heading southwest, and is ~200-300 Myr before peak ram pressure, i.e. its closest approach to M87. The southwestern ridge of enhanced gas surface density and enhanced polarized radio-continuum emission is due to ram pressure compression. It is argued that the faint western Halpha emission ridge is induced by nearly edge-on ram pressure stripping. NGC 4501 represents an especially clear example of early stage ram pressure stripping of a large cluster-spiral galaxy.Comment: 22 pages, 25 figures, accepted for publication in A&

Three-dimensional magnetohydrodynamic simulations of the evolution of magnetic fields in Fanaroff-Riley class II radio sources

Radio observations of Fanaroff-Riley class II sources often show correlations between the synchrotron emission and the linear-polarimetric distributions. Magnetic position vectors seem to align with the projected emission of both the radio jets and the sources' edges. Using statistics we study such relation as well as its unknown time evolution via synthetic polarisation maps of model FR II sources formed in 3D-MHD numerical simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is initially random with a Kolmogorov power spectrum, everywhere. We investigate the structure and evolution of magnetic fields in the sources as a function of the power of jets and the observational viewing angle. Our synthetic polarisation maps agree with observations, showing B-field vectors which are predominantly aligned with the jet axis, and show that magnetic fields inside sources are shaped by the jets' backflow. Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than inside more spherical ones. We see jets increase the magnetic energy in cocoons in proportion to the jet velocity and the cocoon width. Filaments in the synthetic emission maps suggest turbulence develops in evolved sources.Comment: Accepted for publication in the MNRAS. 21 pages, 11 figure

Local Star formation triggered by SN shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a SNR with a neutral magnetized cloud we derived analytically a set of conditions which are favorable for driving gravitational instability in the cloud and thus star formation. We have built diagrams of the SNR radius, versus the cloud density, that constrain a domain in the parameter space where star formation is allowed. The diagrams are also tested with fully 3-D MHD simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity ~1 $mu$G results only a small shrinking of the star formation triggering zone in the diagrams, a larger magnetic field (~10 $mu$G) causes a significant shrinking, as expected. Applications of the diagrams to a few regions of our own galaxy have revealed that star formation in those sites could have been triggered by shock waves from SNRs. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is smaller than the observed values in our own Galaxy (sfe ~0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star forming galaxies, not even when the magnetic field in the neutral clouds is neglected. (abridged)Comment: 19 pages, 13 figures, accepted for pubblication in MNRA