Clusters of galaxies, filled with hot magnetized plasma, are the largest bound objects in existence and an important touchstone in understanding the formation of structures in our Universe. In such clusters, thermal conduction follows field lines, so magnetic fields strongly shape the cluster's thermal history; that some have not since cooled and collapsed is a mystery. In a seemingly unrelated puzzle, recent observations of Virgo cluster spiral galaxies imply ridges of strong, coherent magnetic fields offset from their centre. Here we demonstrate, using three-dimensional magnetohydrodynamical simulations, that such ridges are easily explained by galaxies sweeping up field lines as they orbit inside the cluster. This magnetic drape is then lit up with cosmic rays from the galaxies' stars, generating coherent polarized emission at the galaxies' leading edges. This immediately presents a technique for probing local orientations and characteristic length scales of cluster magnetic fields. The first application of this technique, mapping the field of the Virgo cluster, gives a startling result: outside a central region, the magnetic field is preferentially oriented radially as predicted by the magnetothermal instability. Our results strongly suggest a mechanism for maintaining some clusters in a 'non-cooling-core' state.

In this paper we discuss the effect of recombinations to highly excited states (n > 100) in hydrogen during the cosmological recombination epoch. For this purpose, we developed a new ODE solver for the recombination problem, based on an implicit Gear's method. This solver allows us to include up to 350 l-resolved shells or ~61 000 separate levels in the hydrogen model and to solve the recombination problem for one cosmology in ~27 hours. This is a huge improvement in performance over our previous recombination code, for which a 100-shell computation (5050 separate states) already required ~150 hours on a single processor. We show that for 350 shells down to redshift z ~200 the results for the free electron fraction have practically converged. The final modification in the free electron fraction at z ~200 decreases from about \DeltaNe/Ne ~2.8% for 100 shells to \DeltaNe/Ne ~1.6% for 350 shells. However, the associated changes in the CMB power spectra at large multipoles l are rather small, so that for accurate computations in connection with the analysis of Planck data already ~100 shells are expected to be sufficient. Nevertheless, the total value of \tau could still be affected at a significant level. We also briefly investigate the effect of collisions on the recombination dynamics. With our current estimates for the collisional rates we find a correction of \DeltaNe/Ne ~ -0.088% at z ~ 700, which is mainly caused by l-changing collisions with protons. Furthermore, we present results on the cosmological recombination spectrum, showing that at low frequencies collisional processes are important. However, the current accuracy of collisional rates is insufficient for precise computations of templates for the recombination spectrum at \nu<~1 GHz, and also the effect of collisions on the recombination dynamics suffers from the uncertainty in these rates.

High-resolution X-ray observations have revealed cavities and `cold fronts' with sharp edges in temperature, density, and metallicity within galaxy clusters. Their presence poses a puzzle since these features are not expected to be hydrodynamically stable, or to remain sharp in the presence of diffusion. However, a moving core or bubble in even a very weakly magnetized plasma necessarily sweeps up enough magnetic field to build up a dynamically important sheath around the object; the layer's strength is set by a competition between `plowing up' of field and field lines slipping around the core. We show that a two-dimensional approach to the problem is quite generally not possible. In three dimensions, we show with analytic arguments and in numerical experiments, that this magnetic layer modifies the dynamics of a plunging core, greatly modifies the effects of hydrodynamic instabilities on the core, modifies the geometry of stripped material, and even slows the fall of the core through magnetic tension. We derive an expression for the maximum magnetic field strength, the thickness of the layer, and the opening angle of the magnetic wake. The morphology of the magnetic draping layer implies the suppression of thermal conduction across the layer, thus conserving strong temperature gradients over the contact surface. The intermittent amplification of the magnetic field as well as the injection of MHD turbulence in the wake of the core is identified to be due to vorticity generation within the magnetic draping layer. These results have important consequences for understanding the physical properties and the complex gasdynamical processes of the intra-cluster medium, and apply quite generally to motions through other magnetized environments, e.g., the ISM.

There has been significant recent work which examines a situation where a thin magnetic layer is `draped' over a core merging into a larger cluster; the same process also appears to be at work at a bubble rising from the cluster centre. Such a thin magnetic layer could thermally isolate the core from the cluster medium, but only if the same shear process which generates the layer does not later disrupt it. On the other hand, if the magnetized layer can stabilize against the shear instabilities, then the magnetic layer can have the additional dynamical effect of reducing the shear-driven mixing of the core's material during the merger process. These arguments could equally well to underdense cluster bubbles, which would be even more prone to disruption. While it is well known that magnetic fields can suppress instabilities, it is less clear that a thin layer can suppress instabilities on scales significantly larger than its thickness. Here we consider the stability imparted by a thin magnetized layer. Such a layer can have a significant stabilizing effect even on modes with wavelengths much larger than the thickness of the layer l; to stabilize modes ten times larger requires only that the Alfv\'en speed in the magnetized layer is comparable to the relevant destabilizing velocity -- the shear velocity in the case of pure Kelvin-Helmholtz like instability, or a typical buoyancy velocity in the case of pure Rayleigh-Taylor. We confirm our calculations with two-dimensional numerical experiments using the Athena code.

We consider the competition of the different physical processes that can affect the evolution of a flame bubble in a Type Ia supernovae -- burning, turbulence and buoyancy. Even in the vigorously turbulent conditions of a convecting white dwarf, thermonuclear burning that begins at a point near the center (within 100 km) of the star is dominated by the spherical laminar expansion of the flame, until the burning region reaches kilometers in size. Consequently flames that ignite in the inner ~20 km promptly burn through the center, and flame bubbles anywhere must grow quite large--indeed, resolvable by large-scale simulations of the global system--for significant motion or deformation occur. As a result, any hot-spot that successfully ignites into a flame can burn a significant amount of white dwarf material. This potentially increases the stochastic nature of the explosion compared to a scenario where a simmering progenitor can have small early hot-spots float harmlessly away. Further, the size where the laminar flame speed dominates other relevant velocities sets a characteristic scale for fragmentation of larger flame structures, as nothing--by definition--can easily break the burning region into smaller volumes. This makes possible the development of semi-analytic descriptions of the earliest phase of the propagation of burning in a Type Ia supernovae, which we present here. Our analysis is supported by fully resolved numerical simulations of flame bubbles.

The details of ignition of Type Ia supernovae remain fuzzy, despite the importance of this input for any large-scale model of the final explosion. Here, we begin a process of understanding the ignition of these hotspots by examining the burning of one zone of material, and then investigate the ignition of a detonation due to rapid heating at single point. We numerically measure the ignition delay time for onset of burning in mixtures of degenerate material and provide fitting formula for conditions of relevance in the Type Ia problem. Using the neon abundance as a proxy for the white dwarf metallicity, we then find that ignition times can decrease by ~20% with addition of even 5% of neon by mass. When temperature fluctuations that successfully kindle a region are very rare, such a reduction in ignition time can increase the probability of ignition by orders of magnitude. If the neon comes largely at the expense of carbon, a similar increase in the ignition time can occur. We then consider the ignition of a detonation by an explosive energy input in one localized zone, eg a Sedov blast wave leading to a shock-ignited detonation. Building on previous work on curved detonations, we find that surprisingly large inputs of energy are required to successfully launch a detonation, leading to required matchheads of ~4500 detonation thicknesses - tens of centimeters to hundreds of meters - which is orders of magnitude larger than naive considerations might suggest. This is a very difficult constraint to meet for some pictures of a deflagration-to-detonation transition, such as a Zel'dovich gradient mechanism ignition in the distributed burning regime.

Supernovae of Type Ia are used as standard candles for cosmological observations despite the as yet incomplete understanding of their explosion mechanism. In one model, these events are thought to result from subsonic burning in the core of an accreting Carbon/Oxygen white dwarf that is accelerated through flame wrinkling and flame instabilities. Many such white dwarfs have significant magnetic fields. Here we derive the linear effects of such magnetic fields on one flame instability, the well-known Landau-Darrieus instability. When the magnetic field is strong enough that the flame is everywhere sub-Alfvenic, the instability can be greatly suppressed. Super-Alfvenic flames are much less affected by the field, with flames propagating parallel to the field somewh at destabilized, and flames propagating perpendicular to the field somewhat stabili zed. Trans-Alfvenic parallel flames, however, like trans-Alfvenic parallel shocks, are seen to be non-evolutionary; understanding the behavior of these flames will require careful numerical simulation.

Block-structured AMR meshes are often used in astrophysical fluid simulations, where the geometry of the domain is simple. We consider potential efficiency gains for time sub-cycling, or time refinement (TR), on Berger-Collela and oct-tree AMR meshes for explicit or local physics (such as explict hydrodynamics), where the work per block is roughly constant with level of refinement. We note that there are generally many more fine zones than there are coarse zones. We then quantify the natural result that any overall efficiency gains from reducing the amount of work on the relatively few coarse zones must necessarily be fairly small. Potential efficiency benefits from TR on these meshes are seen to be quite limited except in the case of refining a small number of points on a large mesh -- in this case, the benefit can be made arbitrarily large, albeit at the expense of spatial refinement efficiency.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

Recent Chandra and XMM-Newton observations of galaxy cluster cooling flows have revealed X-ray emission voids of up to 30 kpc in size that have been identified with buoyant, magnetized bubbles. Motivated by these observations, we have investigated the behavior of rising bubbles in stratified atmospheres using the FLASH adaptive-mesh simulation code. We present results from 2-D simulations with and without the effects of magnetic fields, and with varying bubble sizes and background stratifications. We find purely hydrodynamic bubbles to be unstable; a dynamically important magnetic field is required to maintain a bubble's integrity. This suggests that, even absent thermal conduction, for bubbles to be persistent enough to be regularly observed, they must be supported in large part by magnetic fields. Thermal conduction unmitigated by magnetic fields can dissipate the bubbles even faster. We also observe that the bubbles leave a tail as they rise; the structure of these tails can indicate the history of the dynamics of the rising bubble.

Critically understanding the `standard candle'-like behavior of Type Ia supernovae requires understanding their explosion mechanism. One family of models for Type Ia Supernovae begins with a deflagration in a Carbon-Oxygen white dwarf which greatly accelerates through wrinkling and flame instabilities. While the planar speed and behavior of astrophysically-relevant flames is increasingly well understood, more complex behavior, such as the flame's response to stretch and curvature, has not been extensively explored in the astrophysical literature; this behavior can greatly enhance or suppress instabilities and local flame-wrinkling, which in turn can increase or decrease the bulk burning rate. In this paper, we explore the effects of curvature on both nuclear flames and simpler model flames to understand the effect of curvature on the flame structure and speed.

We look at the spreading of localized temperature perturbations in the accreted fuel layer of a non-rotating neutron star. The pressure at the base of the accreted fuel layer is large and the material is only partially degenerate. Any temperature perturbations and resulting pressure gradients will lead to enormous accelerations (both laterally and vertically) of the material in the fuel layer. If the burning of this fuel cannot proceed more rapidly than the spreading of this perturbation, then localized burning cannot take place, and it is likely that the ignition would have to proceed simultaneously throughout the envelope. We present some multidimensional simulations of the spreading of temperature perturbations in a helium atmosphere on a neutron star.

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.

Observed novae abundances and explosion energies estimated from observations indicate that there must be significant mixing of the heavier material of the white dwarf (C+O) into the lighter accreted material (H+He). Accordingly, nova models must incorporate a mechanism that will dredge up the heavier white dwarf material, and fluid motions from an early convection phase is one proposed mechanism. We present results from two-dimensional simulations of classical nova precursor models that demonstrate the beginning of a convective phase during the `simmering' of a Nova precursor. We use a new hydrostatic equilibrium hydrodynamics module recently developed for the adaptive-mesh code FLASH. The two-dimensional models are based on the one-dimensional models of Ami Glasner, and were evolved with FLASH from a pre-convective state to the onset of convection.

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.

We present direct numerical simulations of mixing at Rayleigh-Taylor unstable interfaces performed with the FLASH code, developed at the ASCI/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. We present initial results of single-mode studies in two and three dimensions. Our results indicate that three-dimensional instabilities grow significantly faster than two-dimensional instabilities and that grid resolution can have a significant effect on instability growth rates. We also find that unphysical diffusive mixing occurs at the fluid interface, particularly in poorly resolved simulations.

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.

Recent observations of Cepheids in the Virgo cluster have bolstered the evidence that supports a Hubble constant in 70-90 km/s/Mpc range. This evidence, by and large, probes the expansion of the Universe within 100 Mpc. We investigate the possibility that the expansion rate within this region is systematically higher than the true expansion rate due to the presence of a local, large underdense region or void. We begin by calculating the expected deviations between the locally measured Hubble constant and the true Hubble constant for a variety of models. We also discuss the expected correlations between these deviations and mass fluctuation for the sample volume. We find that the fluctuations are small for the standard cold dark matter as well as mixed dark matter models but can be substantial in a number of interesting and viable nonstandard scenarios. However, deviations in the Hubble flow for a region of radius 200 Mpc are small for virtually all reasonable models. Therefore, methods based on supernovae or the Sunyaev-Zel'dovich effect, which can probe 200 Mpc scales, will be essential in determining the true Hubble constant. We discuss, in detail, the fluctuations induced in the cosmic background radiation by voids at the last scattering surface. In addition, we discuss the dipole and quadrupole fluctuations one would expect if the void enclosing us is aspherical or if we lie off-center.