John A. ZuHone

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid-structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Evolution of FLASH, a multi-physics scientific simulation code for high-performance computing

The FLASH code has evolved into a modular and extensible scientific simulation software system over the decade of its existence. During this time it has been cumulatively used by over a thousand researchers to investigate problems in astrophysics, cosmology, and in some areas of basic physics, such as turbulence. Recently, many new capabilities have been added to the code to enable it to simulate problems in high-energy density physics. Enhancements to these capabilities continue, along with enhancements enabling simulations of problems in fluid-structure interactions. The code started its life as an amalgamation of already existing software packages and sections of codes developed independently by various participating members of the team for other purposes. The code has evolved through a mixture of incremental and deep infrastructural changes. In the process, it has undergone four major revisions, three of which involved a significant architectural advancement. Along the way, a software process evolved that addresses the issues of code verification, maintainability, and support for the expanding user base. The software process also resolves the conflicts arising out of being in development and production simultaneously with multiple research projects, and between performance and portability. This paper describes the process of code evolution with emphasis on the design decisions and software management policies that have been instrumental in the success of the code. The paper also makes the case for a symbiotic relationship between scientific research and good software engineering of the simulation software.

TESTING SECONDARY MODELS FOR THE ORIGIN OF RADIO MINI-HALOS IN GALAXY CLUSTERS

We present an MHD simulation of the emergence of a radio minihalo in a galaxy cluster core in a ?secondary? model, where the source of the synchrotron-emitting electrons is hadronic interactions between cosmic-ray protons with the thermal intracluster gas, an alternative to the ?reacceleration model? where the cosmic ray electrons are reaccelerated by turbulence induced by core sloshing, which we discussed in an earlier work. We follow the evolution of cosmic-ray electron spectra and their radio emission using passive tracer particles, taking into account the time-dependent injection of electrons from hadronic interactions and their energy losses. We find that secondary electrons in a sloshing cluster core can generate diffuse synchrotron emission with luminosity and extent similar to observed radio minihalos. However, we also find important differences with our previous work. We find that the drop in radio emission at cold fronts is less prominent than that in our reacceleration-based simulations, indicating that in this flavor of the secondary model the emission is more spatially extended than in some observed minihalos. We also explore the effect of rapid changes in the magnetic field on the radio spectrum. While the resulting spectra in some regions are steeper than expected from stationary conditions, the change is marginal, with differences in the synchrotron spectral index of 0.15?0.25, depending on the frequency band. This is a much narrower range than claimed in the best-observed minihalos and produced in the reacceleration model. Our results provide important suggestions to constrain these models with future observations.