Engineering topics

Turk, Matthew J.

Publications and source records attributed to Turk, Matthew J..

8 records

DOE OSTI2020

The challenges of theory-software translation

Background: Software is now ubiquitous within research. In addition to the general challenges common to all software development projects, research software must also represent, manipulate, and provide data for complex theoretical constructs. Ensuring this process of theory-software translation is robust is essential to maintaining the integrity of the science resulting from it, and yet there has been little formal recognition or exploration of the challenges associated with it. Methods: We thematically analyse the outputs of the discussion sessions at the Theory-Software Translation Workshop 2019, where academic researchers and research software engineers from a variety of domains, and with particular expertise in high performance computing, explored the process of translating between scientific theory and software. Results: We identify a wide range of challenges to implementing scientific theory in research software and using the resulting data and models for the advancement of knowledge. We categorise these within the emergent themes of design, infrastructure, and culture, and map them to associated research questions. Conclusions: Systematically investigating how software is constructed and its outputs used within science has the potential to improve the robustness of research software and accelerate progress in its development. We propose that this issue be examined within a new research area of theory-software translation, which would aim to significantly advance both knowledge and scientific practice.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Computing environments for reproducibility: Capturing the “Whole Tale”

The act of sharing scientific knowledge is rapidly evolving away from traditional articles and presentations to the delivery of executable objects that integrate the data and computational details (e.g., scripts and workflows) upon which the findings rely. This envisioned coupling of data and process is essential to advancing science but faces technical and institutional barriers. The Whole Tale project aims to address these barriers by connecting computational, data-intensive research efforts with the larger research process transforming the knowledge discovery and dissemination process into one where data products are united with research articles to create "living publications" or tales. The Whole Tale focuses on the full spectrum of science, empowering users in the long tail of science, and power users with demands for access to big data and compute resources. We report here on the design, architecture, and implementation of the Whole Tale environment.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

THE AGORA HIGH-RESOLUTION GALAXY SIMULATIONS COMPARISON PROJECT. II. ISOLATED DISK TEST

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Building a Community-Focused Data Platform

The National Data Service (NDS) presents a broad-reaching vision of integrating discovery, publishing, and sharing of data. Metcalfe's Law states that the value of a network increases with the number of connections in the network, not the number of nodes. The NDS must develop and steward connections not only between data, as is in its vision statement, but between individuals as well. We propose that this be accomplished by explicitly developing it as a platform on which data exploration can be constructed by individual researchers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

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 one, two, and three 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Dwarf galaxies with ionizing radiation feedback. II. Spatially resolved star formation relation

We investigate the spatially resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate (SFR). Because we have self-consistently calculated the location of ionized gas, we are able to make simulated, spatially resolved observations of star formation tracers, such as Hα emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3 × 10{sup 11} M {sub ☉}, we find that the correlation between SFR density (estimated from mock Hα emission) and H{sub 2} density shows large scatter, especially at high resolutions of ≲75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution and because Hα traces hot gas around star-forming regions and is displaced from the H{sub 2} peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces and molecular clouds being dispersed via stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

DWARF GALAXIES WITH IONIZING RADIATION FEEDBACK. I. ESCAPE OF IONIZING PHOTONS

We describe a new method for simulating ionizing radiation and supernova feedback in the analogs of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3 × 10{sup 11} M{sub ☉}, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ∼20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60° from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f{sub esc}(i), and how it evolves as the particle ages. We discover that the average escape fraction f{sub esc} is dominated by a small number of SFMC particles with high f{sub esc}(i). On average, the escape fraction from an SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myr. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump and from a galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

POPULATION III STAR FORMATION IN LARGE COSMOLOGICAL VOLUMES. I. HALO TEMPORAL AND PHYSICAL ENVIRONMENT

We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z {approx} 10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z = 10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z = 30 to z = 10, and ranging in mass from 2.3 Multiplication-Sign 10{sup 5} M{sub Sun} to 1.2 Multiplication-Sign 10{sup 10} M{sub Sun }. At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z {approx} 10.

79 ASTRONOMY AND ASTROPHYSICS↗