Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

Stratospheric aircraft exhaust plume and wake chemistry studies

This report documents progress to date in an ongoing study to analyze and model emissions leaving a proposed High Speed Civil Transport (HSCT) from when the exhaust gases leave the engine until they are deposited at atmospheric scales in the stratosphere. Estimates are given for the emissions, summarizing relevant earlier work (CIAP) and reviewing current propulsion research efforts. The chemical evolution and the mixing and vortical motion of the exhaust are analyzed to track the exhaust and its speciation as the emissions are mixed to atmospheric scales. The species tracked include those that could be heterogeneously reactive on the surfaces of the condensed solid water (ice) particles and on exhaust soot particle surfaces. Dispersion and reaction of chemical constituents in the far wake are studied with a Lagrangian air parcel model, in conjunction with a radiation code to calculate the net heating/cooling. Laboratory measurements of heterogeneous chemistry of aqueous sulfuric acid and nitric acid hydrates are also described. Results include the solubility of HCl in sulfuric acid which is a key parameter for modeling stratospheric processing. We also report initial results for condensation of nitric acid trihydrate from gas phase H2O and HNO3.

Miake-Lye, R. C.↗

Spacecraft Orbital Debris Reentry: Aerothermal Analysis

In the past 40 years, thousands of objects have been placed in Earth orbit and are being tracked. Space hardware reentry survivability must be evaluated to assess risks to human life and property on the ground. The objective of this paper is to present results of a study to determine altitude of demise (burn-up) or survivability of reentering objects. Two NASA/JSC computer codes - Object Reentry Survival Analysis Tool (ORSAT) and Miniature ORSAT (MORSAT) were used to determine trajectories, aerodynamic aerothermal environment, and thermal response of selected spacecraft components. The methodology of the two codes is presented, along with results of a parametric study of reentering objects modeled as spheres and cylinders. Parameters varied included mass, diameter, wall thickness, ballistic coefficient, length, type of material, and mode of tumbling/spinning. Two fragments of a spent Delta second stage undergoing orbital decay, stainless steel cylindrical propellant tank and titanium pressurization sphere, were evaluated with ORSAT and found to survive entry, as did the actual objects. Also, orbital decay reentry predictions of the Japanese Advanced Earth Observing Satellite (ADEOS) aluminum and nickel box-type components and the Russian COSMOS 954 satellite beryllium cylinders were made with MORSAT. These objects were also shown to survive reentry.

Rochelle, Wm. C.↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

Accurate effective fluid approximation for ultralight axions

Ultralight axions are theoretically interesting and phenomenologically rich dark sector candidates, but they are difficult to track across cosmological timescales because of their fast oscillations. We resolve this problem by developing a novel method to evolve them efficiently and accurately. We first construct an exact effective fluid which at late times matches the axion but which evolves in a simple way. We then approximate this evolution with a carefully chosen equation of state and sound speed. With our scheme we find that we can obtain subpercent accuracy for the linear theory suppression of axion density fluctuations relative to that of cold dark matter without tracking even a single complete oscillation of the axion field. Furthermore, we use our technique to test other approximation schemes and to provide a fitting formula for the transfer function for the matter power spectrum in linear theory in axion models. Implementing our approach in existing cosmological axion codes is straightforward and will help unleash the potential of high-precision next-generation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NGEE Arctic Plant Traits: Plant Community Composition, Kougarok Road Mile Marker 64, Seward Peninsula, Alaska, 2016

This dataset reports the characteristics of the dominant vegetation communities at the NGEE Arctic Kougarok field site, Kougarok Road Mile Marker 64 on the Seward Peninsula. Selected plots were surveyed from 18-23 July 2016. Environmental data (e.g., elevation, slope, soil moisture regime, disturbance type and degree, mean canopy height, etc.) were also recorded for each plot. Elevation measurements were updated 2020-08-18. All plant species (vascular plants, bryophytes and lichens) were recorded along with their percent cover in the plots as determined visually by the lead author. Plots were chosen subjectively in areas of homogeneous and representative vegetation and varied in size from 1-25 m2 depending on canopy structure and height. The site is an east-facing hillslope with vegetation that varies from the summit to the toeslope. The 30 total plots were selected as five replicate plots sampled from each of the six identified habitat types: 1) non-acidic mountain complex at the hillcrest, 2) dwarf-shrub lichen tundra on the shoulder, 3) alder shrublands predominantly on an elevational band across the backslope, 4) willow-birch tundra on the upper backslope, 5) tussock tundra on the footslope in inter-water tracks, and 6) tussock tundra mixed with willow-birch tundra or alder savanna in poor developed water tracks on the footslope. This dataset is comprised of two comma-separated (*.csv) files containing species and environmental data for the plant community composition plots. It also contains seven zipped folders of plot photographs, one map showing plot locations (*.pdf), one blank datasheet with the key to codes used in the field (*.pdf), select photos organized in a (*.pdf), and one User file (*.pdf). The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Using Microwave and Infrared Radiances from Off-Nadir Pixels: Application of Radiative Transfer to Slanted Line-of-Sight and Comparisons with NASA EOS Aqua Data

The passive infrared and microwave nadir sounders such as (A)TOVS observe the atmosphere from a polar orbit by directing their scan pointed at the ground up to about 49 degrees from nadir. Except for the pixels located right on the satellite ground track, the radiance measurements collected by these instruments characterize hence atmospheric emission paths which are slanted with respect to the zenithal direction at the ground. At the outer swath edges, the deviations from nadir reach about 60 degrees in terms of Satellite Zenith Angle (SZA). The radiative transfer codes used in operational Numerical Weather Prediction applications make the appropriate corrections to account for the extra path induced by the non-zero SZA. However, no corrections are made to account for the fact that the atmospheric profiles along the slanted line-of-sight (LOS) are different from the vertical because of horizontal gradients in the atmosphere. Using NASA EOS Aqua satellite's orbits, zenith and azimuth angles, as well as three-dimensional fields of temperature, water vapor, and ozone produced by the NASA Global Modeling and Assimilation Office, we extracted slanted atmospheric profiles for actual soundings performed by the AIRS and AMSU-A instruments onboard EOS Aqua. We will present the results of our study comparing the calculated brightness temperatures along slanted LOS and vertical LOS with AIRS and AMSU-A observations.

Poli, Paul↗

Analysis of Air Traffic Track Data with the AutoBayes Synthesis System

The Next Generation Air Traffic System (NGATS) is aiming to provide substantial computer support for the air traffic controllers. Algorithms for the accurate prediction of aircraft movements are of central importance for such software systems but trajectory prediction has to work reliably in the presence of unknown parameters and uncertainties. We are using the AutoBayes program synthesis system to generate customized data analysis algorithms that process large sets of aircraft radar track data in order to estimate parameters and uncertainties. In this paper, we present, how the tasks of finding structure in track data, estimation of important parameters in climb trajectories, and the detection of continuous descent approaches can be accomplished with compact task-specific AutoBayes specifications. We present an overview of the AutoBayes architecture and describe, how its schema-based approach generates customized analysis algorithms, documented C/C++ code, and detailed mathematical derivations. Results of experiments with actual air traffic control data are discussed.

Schumann, Johann Martin Philip↗

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗

The Rodinia Benchmark Suite in SYCL

As opposed to the Open Computing Language (OpenCL) programming model in which host and device codes are generally written in different programming languages, SYCL can combine host and device codes for an application in a type-safe way to improve development productivity and performance portability. Rodinia is a widely used benchmark suite for heterogeneous computing. Hence, we port the OpenCL implementations of the benchmark suite to SYCL manually. The SYCL benchmark suite is an opensource project for tracking the development of the mainstream SYCL compilers, and for developers and researchers interested in programming productivity, performance analysis, and portability across different computing platforms. We organize the remainder of the report as follows. Section II introduces the SYCL programming model, shows the major differences between an OpenCL program and a SYCL program, and gives a summary of the benchmark suite. In Section III, we describe our SYCL implementations in more details. We evaluate the SYCL benchmarks on the Intel® central processing units (CPUs) and graphics processing units (GPUs) in Section IV. Section V concludes the report.

97 MATHEMATICS AND COMPUTING↗

A Green's function method for heavy ion beam transport

The use of Green's function has played a fundamental role in transport calculations for high-charge high-energy (HZE) ions. Two recent developments have greatly advanced the practical aspects of implementation of these methods. The first was the formulation of a closed-form solution as a multiple fragmentation perturbation series. The second was the effective summation of the closed-form solution through nonperturbative techniques. The nonperturbative methods have been recently extended to an inhomogeneous, two-layer transport media to simulate the lead scattering foil present in the Lawrence Berkeley Laboratories (LBL) biomedical beam line used for cancer therapy. Such inhomogeneous codes are necessary for astronaut shielding in space. The transport codes utilize the Langley Research Center atomic and nuclear database. Transport code and database evaluation are performed by comparison with experiments performed at the LBL Bevalac facility using 670 A MeV 20Ne and 600 A MeV 56Fe ion beams. The comparison with a time-of-flight and delta E detector measurement for the 20Ne beam and the plastic nuclear track detectors for 56Fe show agreement up to 35%-40% in water and aluminium targets, respectively.

NASA Discipline Radiation Health↗

A Monte-Carlo Model for the Formation of Radiation-induced Chromosomal Aberrations

Purpose: To simulate radiation-induced chromosome aberrations in mammalian cells (e.g., rings, translocations, and dicentrics) and to calculate their frequency distributions following exposure to DNA double strand breaks (DSBs) produced by high-LET ions. Methods: The interphase genome was assumed to be comprised of a collection of 2 kbp rigid-block monomers following the random-walk geometry. Additional details for the modeling of chromosomal structure, such as chromosomal domains and chromosomal loops, were included. A radial energy profile for heavy ion tracks was used to simulate the high-LET pattern of induced DSBs. The induced DSB pattern depended on the ion charge and kinetic energy, but always corresponded to the DSB yield of 25 DSBs/cell/Gy. The sum of all energy contributions from Poisson-distributed particle tracks was taken to account for all possible one-track and multi-track effects. The relevant output of the model was DNA fragments produced by DSBs. The DSBs, or breakpoints, were defined by (x, y, z, l) positions, where x, y, z were the Euclidian coordinates of a DSB, and where l was the relative position along the genome. Results: The code was used to carry out Monte Carlo simulations for DSB rejoinings at low doses. The resulting fragments were analyzed to estimate the frequencies of specific types of chromosomal aberrations. Histograms for relative frequencies of chromosomal aberrations and P.D.F.s (probability density functions) of a given aberration type were produced. The relative frequency of dicentrics to rings was compared to empirical data to calibrate rejoining probabilities. Of particular interest was the predicted distribution of ring sizes, irrespective of their frequencies relative to other aberrations. Simulated ring sizes were . 4 kbp, which are far too small to be observed experimentally (i.e., by microscopy) but which, nevertheless, are conjectured to exist. Other aberrations, for example, inversions, translocations, as well as multi-centrics were also recorded. Conclusion: High-LET DNA damage affects the frequencies of chromosomal aberrations. The ratio of rings to dicentrics is correct for the genomic size cut-offs corresponding to available experimental data. The present work predicts a relative abundance of small rings following irradiation by heavy ions.

Ponomarev, Artem L.↗

NASA Tech Briefs, July 2012

Topics covered include: Instrument Suite for Vertical Characterization of the Ionosphere-Thermosphere System; Terahertz Radiation Heterodyne Detector Using Two-Dimensional Electron Gas in a GaN Heterostructure; Pattern Recognition Algorithm for High-Sensitivity Odorant Detection in Unknown Environments; Determining Performance Acceptability of Electrochemical Oxygen Sensors; Versatile Controller for Infrared Lamp and Heater Arrays; High-Speed Scanning Interferometer Using CMOS Image Sensor and FPGA Based on Multifrequency Phase-Tracking Detection; Ultra-Low-Power MEMS Selective Gas Sensors; Compact Receiver Front Ends for Submillimeter-Wave Applications; Dynamically Reconfigurable Systolic Array Accelerator; Blocking Losses With a Photon Counter; Motion-Capture-Enabled Software for Gestural Control of 3D Mod; Orbit Software Suite; CoNNeCT Baseband Processor Module Boot Code SoftWare (BCSW); Trajectory Software With Upper Atmosphere Model; ALSSAT Version 6.0; Employing a Grinding Technology to Assess the Microbial Density for Encapsulated Organisms; Demonstration of Minimally Machined Honeycomb Silicon Carbide Mirrors; Polyimide Aerogel Thin Films; Nanoengineered Thermal Materials Based on Carbon Nanotube Array Composites; Composite Laminate With Coefficient of Thermal Expansion Matching D263 Glass; Robust Tensioned Kevlar Suspension Design; Focal Plane Alignment Utilizing Optical CMM; Purifying, Separating, and Concentrating Cells From a Sample Low in Biomass; Virtual Ultrasound Guidance for Inexperienced Operators; Beat-to-Beat Blood Pressure Monitor; Non-Contact Conductivity Measurement for Automated Sample Processing Systems; An MSK Radar Waveform; Telescope Alignment From Sparsely Sampled Wavefront Measurements Over Pupil Subapertures; Method to Remove Particulate Matter from Dusty Gases at Low Pressures; Terahertz Quantum Cascade Laser With Efficient Coupling and Beam Profile; Measurement Via Optical Near-Nulling and Subaperture Stitching; 885-nm Pumped Ceramic Nd:YAG Master Oscillator Power Amplifier Laser System; Airborne Hyperspectral Imaging System; Heat Shield Employing Cured Thermal Protection Material Blocks Bonded in a Large-Cell Honeycomb Matrix; and Asymmetric Supercapacitor for Long-Duration Power Storage.

Source record↗

The Accelerator Neutrino Neutron InteractionExperiment: Design, Construction and Preparationfor First Physics

This thesis covers work performed for the Accelerator Neutrino Neutron Interaction Experiment (ANNIE) between 2015 and 2020. ANNIE is a beamline neutrino experiment located at the Fermi National Accelerator Laboratory in Chicago. Its primary physics goals are the measurement of neutrino interaction cross-sections on oxygen, and a measurement of the multiplicity of final state neutrons from these interactions, both as functions of final state kinematics. In order to achieve these goals ANNIE will be pioneering the use of Large Area Picosecond Photodetectors and the use of gadolinium doping to enhance neutron visibility in a water Cherenkov detector. Over the time period of this thesis the ANNIE detector was constructed, an initial 'Phase I' background run was completed, the detector was decommissioned and then substantially upgraded in preparation for the main 'Phase II' physics run. At the time of writing ANNIE is being re-commissioned, with Phase II data taking anticipated to begin in 2021. During this time I spent 21/2 years on-site contributing to almost all aspects of hardware work, while performing software development for the simulation and analysis frameworks. As my first software task I developed Monte Carlo simulations of the detector response, based on the WCSim framework. While thismore » framework provided a starting point for the tank simulation significant modifications were required to implement the complete ANNIE detector system. This included a rewrite of the geometry, along with modifications to the photosensor definitions, hit collections, digitization, triggering, input and output interfaces, etc. To take advantage of improved modelling of neutrons in Geant4.10 the physics lists were also reimplemented in the new version format. As part of this work I determined performance metrics based on the output files, calculating kinematic acceptance, detector efficiency and expected event rates. My contributions on reconstruction focused primarily on the MRD, for which I implemented two track reconstruction algorithms, specialised for high and low energy muons. The high energy algorithm achieves >90% efficiency for tracks with at least 8 hits, and is robust against dead paddles and unrelated hits. The low energy algorithm performs a search in combination with tank information to provide sensitivity down to a single MRD hit. Early data handling work for the Phase I background run included writing the code for handling TDC hits and performing studies on neutron detection without the requirement of multi-hit coincidence. As the collaboration moved to a new framework for the Phase II analysis I performed a substantial amount of work to develop the codebase, implementing tools to perform tasks such as file IO, data parsing, waveform processing, event building, timestamp alignment, time clustering, hit residual calculation, event visualisation etc. I helped develop the python API for integrating multi-language tools, and ported the framework into a Singularity container to provide a self-contained environment that…« less

43 PARTICLE ACCELERATORS↗

BatAnalysis - A Comprehensive Python Pipeline for Swift BAT Survey Analysis

The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80- channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, opensource pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ~18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗