Source author record

David C. Collins

David C. Collins appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal(rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multi-phase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(σ^{2 Γ}_{ρ/\langle ρ\rangle}-1)^{1/2} (1-σ_{ρ\langle ρ\rangle}^{-2})^{-1/2}\mathcal{M}^{-1}Γ^{-1/2}$ with the density dispersion $σ_{ρ/ \langle ρ\rangle}$, the turbulent Mach number $\mathcal{M}$, and the polytropic exponent $Γ$. Purely solenoidal and purely compressive driving correspond to $b \sim 1/3$ and $b \sim 1$, respectively. Using simulations in which a shock is driven into a multi-phase medium with structures of different sizes and $Γ< 1$, we find $b \sim 1$ for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

preprint2022arXiv

Turbulence Generation by Shock Interaction with a Highly Non-Uniform Medium

An initially planar shock wave propagating into a medium of non-uniform density will be perturbed, leading to the generation of post-shock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly-non-uniform density (order-unity normalized variance, $σ_ρ/\overlineρ \sim 1$) and strong shocks (shock Mach numbers $\overline{M}_s \gtrsim 10$). This leads to a highly disrupted shock and a turbulent post-shock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the post-shock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the pre-shock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The post-shock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

preprint2021arXiv

Physics of Thermonuclear Explosions: Magnetic Field Effects on Deflagration Fronts and Observable Consequences

We present a study of the influence of magnetic field strength and morphology in Type Ia Supernovae and their late-time light curves and spectra. In order to both capture self-consistent magnetic field topologies as well evolve our models to late times, a two stage approach is taken. We study the early deflagration phase (1s) using a variety of magnetic field strengths, and find that the topology of the field is set by the burning, independent of the initial strength. We study late time (~1000 days) light curves and spectra with a variety of magnetic field topologies, and infer magnetic field strengths from observed supernovae. Lower limits are found to be 106G. This is determined by the escape, or lack thereof, of positrons that are tied to the magnetic field. The first stage employs 3d MHD and a local burning approximation, and uses the code Enzo. The second stage employs a hybrid approach, with 3D radiation and positron transport, and spherical hydrodynamics. The second stage uses the code HYDRA. In our models, magnetic field amplification remains small during the early deflagration phase. Late-time spectra bear the imprint of both magnetic field strength and morphology. Implications for alternative explosion scenarios are discussed.

preprint2019arXiv

The Power Spectra of Polarized, Dusty Filaments

We develop an analytic model for the power spectra of polarized filamentary structures as a way to study the Galactic polarization foreground to the Cosmic Microwave Background. Our approach is akin to the cosmological halo-model framework, and reproduces the main features of the Planck 353 GHz power spectra. We model the foreground as randomly-oriented, three-dimensional, spheroidal filaments, accounting for their projection onto the sky. The main tunable parameters are the distribution of filament sizes, the filament physical aspect ratio, and the dispersion of the filament axis around the local magnetic field direction. The abundance and properties of filaments as a function of size determine the slopes of the foreground power spectra, as we show via scaling arguments. The filament aspect ratio determines the ratio of $B$-mode power to $E$-mode power, and specifically reproduces the Planck-observed dust ratio of one-half when the short axis is roughly one-fourth the length of the long axis. Filament misalignment to the local magnetic field determines the $TE$ cross-correlation, and to reproduce Planck measurements, we need a (three-dimensional) misalignment angle with a root mean squared dispersion of about 50 degrees. These parameters are not sensitive to the particular filament density profile. By artificially skewing the distribution of the misalignment angle, this model can reproduce the Planck-observed (and parity-violating) $TB$ correlation. The skewing of the misalignment angle necessary to explain $TB$ will cause a yet-unobserved, positive $EB$ dust correlation, a possible target for future experiments.

preprint2015arXiv

Observational Diagnostics of Self-Gravitating MHD Turbulence in Giant Molecular Clouds

We study the observable signatures of self-gravitating MHD turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term "early," "intermediate," and "advanced." At early times, i.e. $t<0.15t_{ff}$, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e. $0.15t_{ff}< t \leq 0.35t_{ff}$, there exists signatures of the prestellar cores in the shallower PDF and power spectrum power law slopes. The column density PDF power law tails at these times have line of sight averaged slopes ranging from -2.5 to -1.5 with shallower values belonging to simulations with lower magnetic field strength. The density power spectrum slope becomes shallow and can be characterized by $P(k)=A_1k^{β_2}e^{-k/k_c}$, where $A_1$ describes the amplitude, $k^{β_2}$ describes the classical power law behavior and the scale $k_c$ characterizes the turn over from turbulence dominated to self-gravity dominated. At advanced stages of collapse, i.e. $\approx t>0.35t_{ff}$, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.

preprint2015arXiv

Self-generated turbulence in magnetic reconnection

Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number ($\Lund$). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional instability in the plane of a current sheet between reversing field lines without a guide field. The instability leads to rapid reconnection of magnetic field lines at a rate independent of $\Lund$ over at least the range $3.2\times 10^3 \lesssim \Lund \lesssim 3.2 \times 10^5$ resolved by the simulations. We find that this instability occurs even for cases with $\Lund \lesssim 10^4$ that in our models appear stable to the recently described, two-dimensional, plasmoid instability. Our results suggest that three-dimensional, MHD processes alone produce fast (resistivity independent) reconnection without recourse to kinetic effects or external turbulence. The unstable reconnection layers provide a self-consistent environment in which the extensively studied turbulent reconnection process can occur.

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

preprint2012arXiv

Comparisons of Cosmological MHD Galaxy Cluster Simulations to Radio Observations

Radio observations of galaxy clusters show that there are $μ$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($σ_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $σ_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

preprint2012arXiv

Cosmological MHD Simulations of Galaxy Cluster Radio Relics: Insights and Warnings for Observations

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, ρ~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 μG and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

preprint2011arXiv

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters II. The Effects of Cluster Size and Dynamical State

Theory and simulations suggest that magnetic fields from radio jets and lobes powered by their central super massive black holes can be an important source of magnetic fields in the galaxy clusters. This is paper II in a series of studies where we present self-consistent high-resolution adaptive mesh refinement cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial masses ranging from 1 $\times$ 10$^{14}$ to 2 $\times$ 10$^{15}$ M$_{\odot}$. In this work we examine the effects of the mass and merger history on the final magnetic properties. We find that the evolution of magnetic fields is qualitatively similar to those of previous studies. In most clusters, the injected magnetic fields can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process with hierarchical mergers, while the amplification history and the magnetic field distribution depend on the cluster formation and magnetism history. This can be very different for different clusters. The total magnetic energies in these clusters are between 4 $\times$ 10$^{57}$ and $10^{61}$ erg, which is mainly decided by the cluster mass, scaling approximately with the square of the total mass. Dynamically older relaxed clusters usually have more magnetic fields in their ICM. The dynamically very young clusters may be magnetized weakly since there is not enough time for magnetic fields to be amplified.

preprint2010arXiv

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters. I. The Effect of Injection Energy and Redshift

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

preprint2010arXiv

Mass and Magnetic distributions in Self Gravitating Super Alfvenic Turbulence with AMR

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(ρ) \propto ρ^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto ρ^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.

preprint2010arXiv

The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships

Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.

preprint2009arXiv

Cosmological AMR MHD with Enzo

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

preprint2009arXiv

Turbulence and Dynamo in Galaxy Cluster Medium: Implications on the Origin of Cluster Magnetic Fields

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($μ$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.