SEARCH · Engineering Papers

Results for “Source code instrumentation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS2011

An Icy Kuiper-Belt Around the Young Solar-Type Star HD 181327

HD 181327 is a young Main Sequence F5/F6 V star belonging to the Beta Pictoris moving group (age approx 12 Myr). It harbors an optically thin belt of circumstellar material at approx90 AU, presumed to result from collisions in a populat.ion of unseen planetesimals. Aims. We aim to study the dust properties in the belt in great details, and to constrain the gas-to-dust ratio. Methods. We obtained far-IR photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space Observatory, complemented by new 3.2 nun observations carried with the ATCA array. The geometry of the belt is constrained with newly reduced HST /NICMOS scattered light images that break the degeneracy between the disk geometry and the dust properties. We then use the radiative transfer code GRaTer to compute a large grid of dust models, and we apply a Bayesian inference method to identify the grain models that best reproduce the SED. We attempt to detect the oxygen and ionized carbon fine-structure lines with Herschel/PACS spectroscopy, providing observables to our photochemical code ProDiMo. Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an import.ant layer of ice for a total dust mass of approx 0.05 stellar Mass. We discuss evidences that the grains consists of fluffy aggregates. The upper limits on the gas atomic lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than approx 17 Stellar Mass Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Edgeworth-Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

Lebreton, J.↗

DOE OSTI2026

Machine learning at the Spallation Neutron Source accelerator and target

We describe the ongoing efforts to apply Machine Learning techniques to improve the performance of our accelerator and target. Specially, we are looking to minimize halo beam losses in the absence of a proper physics model, automatically detect and log anomalies in the target support systems such as cooling, and detect and prevent errant beam pulses in the linac. We also describe the infrastructure we use to acquire and stream data to the GPU cluster for training, our code development cycle, and edge computing for model inference. To minimize halo beam losses, we use a Reinforcement Learning technique tested on a virtual accelerator. The target anomaly detection is trained on archived data using incomplete physics models and is made part of the existing target reporting system. The errant beam prevention analyzes beam current and beam phase waveforms as well as accelerator configuration data to predict errant pulses. We also develop continual learning to adapt to changes in the accelerator.

Accelerator Physics↗

DOE OSTI2021

Expanded Capabilities and Opportunities for Virtual Beamline Code

In the 25 years since LLNL broke ground on the National Ignition Facility (NIF) - a cornerstone of the National Nuclear Security Administration’s Stockpile Stewardship Program - the NIF and Photon Science (NIF&PS) Directorate has steadily pushed the boundaries of laser physics, nonlinear optics, and photonics in service of inertial confinement fusion experiments, and advanced photon sources development. Beyond the 192 beamlines of the primary laser facility with its world-record energy and fusion yields, NIF&PS capabilities have expanded to include chirped-pulse amplification, kilojoule petawatt-class short-pulse systems to generate hard x-ray radiographic probes, high-average power lasers as secondary sources for generation or wake field particle acceleration along directed energy applications, and other scientific pursuits. The Computing Directorate partners with NIF&PS on numerous projects and technologies from diagnostic measurements and high performance control systems to data analysis, information technology infrastructure, and scientific simulation code development. “This strong collaboration enables cutting-edge science, particularly with regard to physics-based modeling and simulation,” says physicist Jean-Michel Di Nicola. One crucial capability for NIF&PS is the Virtual Beamline (VBL) laser simulation code. It can model all the major laser physics and technology involved in the design optimization, commissioning, and operations of advanced laser architectures from tabletop to NIF-scale. VBL provides researchers with high-fidelity models and high-resolution calculations of laser performance predictions—including for the entire NIF laser system, the Advanced Radiographic Capability, parts of the High Repetition-Rate Advanced Petawatt Laser System laser delivered to the Czech Republic, and the Optical Science Laser. After more than two decades of experimentally verified physics and computing enhancements, this unique workhorse code is wrapping up another major milestone: migration from Java to C++ with a flurry of user interface and optimization features, as well as laser physics and high-resolution enhancements thanks to parallel execution on Livermore Computing platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

HarDWR - Harmonized Water Rights Records

Economics↗

DOE OSTI2015

Simulation of nonlinear propagation of biomedical ultrasound using PZFlex and the KZK Texas code

In biomedical ultrasound nonlinear acoustics can be important in both diagnostic and therapeutic applications and robust simulations tools are needed in the design process but also for day-to-day use such as treatment planning. For most biomedical application the ultrasound sources generate focused sound beams of finite amplitude. The KZK equation is a common model as it accounts for nonlinearity, absorption and paraxial diffraction and there are a number of solvers available, primarily developed by research groups. We compare the predictions of the KZK Texas code (a finite-difference time-domain algorithm) to an FEM-based commercial software, PZFlex. PZFlex solves the continuity equation and momentum conservation equation with a correction for nonlinearity in the equation of state incorporated using an incrementally linear, 2nd order accurate, explicit algorithm in time domain. Nonlinear ultrasound beams from two transducers driven at 1 MHz and 3.3 MHz respectively were simulated by both the KZK Texas code and PZFlex, and the pressure field was also measured by a fibre-optic hydrophone to validate the models. Further simulations were carried out a wide range of frequencies. The comparisons showed good agreement for the fundamental frequency for PZFlex, the KZK Texas code and the experiments. For the harmonic components, the KZK Texas code was in good agreement with measurements but PZFlex underestimated the amplitude: 32% for the 2nd harmonic and 66% for the 3rd harmonic. The underestimation of harmonics by PZFlex was more significant when the fundamental frequency increased. Furthermore non-physical oscillations in the axial profile of harmonics occurred in the PZFlex results when the amplitudes were relatively low. These results suggest that careful benchmarking of nonlinear simulations is important.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2018

Description of the code ANVIL (ANisotropic Vhf Impulse Location)

In this paper, we have developed a general, three-dimensional method to locate sources of earth-directed radiation which takes into account the flight paths of reflected signals. While time-of-arrival algorithms exist for locating radio sources using line-of-sight propagation paths, radio sources originating above the surface of the earth will not necessarily emit strong power along direct paths to satellites. A combination of direct and ground-reflected pulses or only ground-reflected signals from such sources could be received by satellite-borne sensors. The work presented here applies to satellite detection of sub-ionospheric radio sources in a vacuum environment with ideal reflection off the earth's surface. Because satellites are not static receivers, their configuration with respect to a radio source is not always optimal. Therefore, a statistical study is performed using 1000 randomized configurations of 24 satellites in middle earth orbit. For each configuration, the radio source latitude and longitude is fixed, and its altitude is varied from 1 to 97 km. An analytic direct-path algorithm using the five satellites nearest the radio source provides an initial guess for the radio source latitude and longitude. We find that, using this approximation, the mean error in the calculated nadir location has a maximum value of 134 m for the 97-km source altitude. The coordinates of the initial guess are used to define the origin of a grid within which an all-points search refinement is performed. Finally, using this procedure, the overall maximum mean error in radio source position is found to be on the order of the computational grid size.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a series of various files organized by state sub-directories. The first two characters of each file name is the abbreviation for the state the in which the file contains data for. After the abbreviation is the text which describes the contents of the file. Here is each file type described in detail: XXFullHarmonizedRights.csv: A file of the combined groundwater and surface water records for each state. Essentially, this file is the merging of XXGroundwaterHarmonizedRights.csv and XXSurfaceWaterHarmonizedRights.csv by state. The column headers for each of this type of file are: state - The name of the state the data comes from. FIPS - The two-digit numeric state ID code. waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. priorityDate - The priority date associated with the right. origWaterUse - The original stated water use(s) from the state. waterUse - The water use category under the unified use categories established here. source - Whether the right is for surface water or groundwater. basinNum - The alpha-numeric identifier of the WMA the record belongs to. CFS - The maximum flow of the allocation in cubic feet per second (ft3s-1). Arizona is unique among the states, as its surface and groundwater resources are managed with two different sets of boundaries. So, for Arizona, the basinNum column is missing and instead there are two columns: surBasinNum - The alpha-numeric identifier of the surface water WMA the record belongs to. grdBasinNum - The alpha-numeric identifier of the groundwater WMA the record belongs to. XXStatePOD.shp: A shapefile which identifies the location of the Points of Diversion for the state's water rights. It should be noted that not all water right records in XXFullHarmonizedRights.csv have coordinates, and therefore may be missing from this file. XXStatePOU.shp: A shapefile which contains the area(s) in which each water right is claimed to be used. Currently, only Idaho and Washington provided valid data to include within this file. XXGroundwaterHarmonizedRights.csv: A file which contains only harmonized groundwater rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. XXSurfaceWaterHarmonizedRights.csv: A file which contains only harmonized surface water rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. Additionally, one file, stateWMALabels.csv, is not stored within a sub-directory. While we have referred to the spatial boundaries that each state uses to manage its water resources as WMAs, this term is not shared across all states. This file lists the proper name for each boundary set, by state. For those whom may be interested in exploring our code more in depth, we are also making available an internal data file for convenience. The file is in .RData format and contains everything described above as well as some minor additional objects used within the code calculating the cumulative curves. For completeness, here is a detailed description of the various objects which can be found within the .RData file: states: A character vector containing the state names for those states in which data was collected for. More importantly, the index of the state name is also the index in which that state's data can be found in the various following list objects. For example, if California is the third index in this object, the data for California will also be in the third index for each accompanying list. rightsByState_ground: A list of data frames with the cleaned ground water rights collected from each state. This object holds the the data that is exported to created the xxGroundwaterHarmonizedRights.csv files. rightsByState_surface: A list of data frames with the cleaned surface water rights collected from each state. This object holds the the data that is exported to created the xxSurfaceWaterHarmonizedRights.csv files. fullRightsRecs: A list of the combined groundwater and surface water records for each state. This object holds the the data that is exported to created the xxFullHarmonizedRights.csv files. projProj: The spatial projection used for map creation in the beginning of the project. Specifically, the World Geodetic System (WGS84) as a coordinate reference system (CRS) string in PROJ.4 format. wmaStateLabel: The name and/or abbreviation for what each state legally calls their WMAs. h2oUseByState: A list of spatial polygon data frames which contain the area(s) in which each water right is claimed to be used. It should be noted that not all water right records have a listed area(s) of use in this object. Currently, only Idaho and Washington provided valid data to be included in this object. h2oDivByState: A list of spatial points data frames which identifies the location of the Point of Diversion for the state's water rights. It should be noted that not all water right records have a listed Point of Diversion in this object. spatialWMAByState: A list of spatial polygon data frames which contain the spatial WMA boundaries for each state. The only data contained within the table are identifiers for each polygon. It is worth reiterating that Arizona is the only state in which the surface and groundwater WMA boundaries are not the same. wmaIDByState: A list which contains the unique ID values of the WMAs for each state. plottingDim: A character vector used to inform mapping functions for internal map making. Each state is classified as either "tall" or "wide", to maximize space on a typical 8x11 page. The code related to the creation of this dataset can be viewed within HarDWR GitHub Repository/dataHarmonization.

Economics↗

NASA NTRS2012

An Icy Kuiper Belt Around the Young Solar-type Star HD 181327

Context. HD 181327 is a young main sequence F5/F6 V star belonging to the Beta Pictoris moving group (age approx.. 12 Myr). It harbors an optically thin belt of circumstellar material at radius approx.. 90 AU, presumed to result from collisions in a population of unseen planetesimals. Aims. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. Methods. We obtained far-infrared photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space Observatory, complemented by new 3.2 mm observations carried with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS scattered light images that allow the degeneracy between the disk geometry and the dust properties to be broken. We then use the radiative transfer code GRaTeR to compute a large grid of models, and we identify the grain models that best reproduce the spectral energy distribution (SED) through a Bayesian analysis. We attempt to detect the oxygen and ionized carbon fine-structure lines with Herschel/PACS spectroscopy, providing observables to our photochemical code ProDiMo. Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes.We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of approx.. 0.05 Solar Mass (in grains up to 1 mm). We discuss evidences that the grains consists of fluffy aggregates. The upper limits on the gas atomic lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than approx. 17 Solar Mass. Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Edgeworth-Kuiper belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

HD 181327↗

DOE OSTI2017

Electron Temperature Fluctuation Measurements and Transport Model Validation at Alcator C-Mod

The tokamak is a type of toroidal device used to confine a fusion plasma using large magnetic fields. Tokamaks and stellarators the leading devices for confining plasmas for fusion, and the capability to predict performance in these magnetically confined plasmas is essential for developing a sustainable fusion energy source. The magnetic configuration of tokamaks and stellarators does not exist in Nature, yet, the fundamental processes governing transport in fusion plasmas are universal – turbulence and instabilities, driven by inhomogeneity and asymmetry in the plasma, conspire to transport heat and particles across magnetic field lines and can play critical roles in impurity confinement and generation of intrinsic rotation. Turbulence exists in all plasmas, and in neutral fluids as well. The study of turbulence is essential to developing a fundamental understanding of the nature of the fourth state of matter, plasmas. Experimental studies of turbulence in tokamaks date back to early scattering observations from the late 1970s. Since that time, great advances in turbulence diagnostics have been made, all of which have significantly enhanced our knowledge and understanding of turbulence in tokamaks. Through comparisons with advanced gyrokinetic theory and turbulent-transport models a great deal of evidence exists to implicate turbulent-driven transport as an important mechanism determining transport in all channels: heat, particle and momentum However, prediction and control of turbulent-driven transport remains elusive. Key to development of predictive transport models for magnetically confined fusion plasmas is validation of the nonlinear gyrokinetic transport model, which describes transport due to turbulence. Validation of gyrokinetic codes must include detailed and quantitative comparisons with measured turbulence characteristics, in addition to comparisons with inferred transport levels and equilibrium profiles. For this reason, advanced plasma diagnostics for studying core turbulence are needed in order to assess the accuracy of gyrokinetic models for turbulent-driven particle, heat and momentum transport. New core turbulence diagnostics at the world-class tokamaks Alcator C-Mod at MIT and ASDEX Upgrade at the Max Planck Institute for Plasma Physics have been designed, developed, and operated over the course of this project. These new instruments are capable of measuring electron temperature fluctuations and the phase angle between density and temperature fluctuations locally and quantitatively. These new data sets from Alcator C-Mod and ASDEX Upgrade are being used to fill key gaps in our understanding of turbulent transport in tokamaks. In particular, this project has results in new results on the topics of the Transport Shortfall, the role of ETG turbulence in tokamak plasmas, profile stiffness, the LOC/SOC transition, and intrinsic rotation reversals. These data are used in a rigorous process of “Transport model validation”, and this group is a world-leader on using turbulence models to design new hardware and new experiments at tokamaks. A correlation electron cyclotron emission (CECE) diagnostic is an instrument used to measure micro-scale fluctuations (mm-scale, compared to the machine size of meters) of electron temperature in magnetically confined fusion plasmas, such as those in tokamaks and stellarators. These micro-scale fluctuations are associated with drift-wave type turbulence, which leads to enhanced cooling and mixing of particles in fusion plasmas and limits achieving the required temperatures and densities for self-sustained fusion reactions. A CECE system can also be coupled with a reflectometer system that measured micro-scale density fluctuations, and from these simultaneous measurements, one can extract the phase between the density (n) and temperature (T) fluctuations, creating an nT phase diagnostic. Measurements of the fluctuations and the phase angle between them are extremely useful for testing and validating predictive models for the transport of heat and particles in fusion plasmas due to turbulence. Once validated, the models are used to predict performance in ITER and other burning plasmas, such as the MIT ARC design. Most recently, data from the newly developed, so-called “CECE diagnostic” [Cima 1995, White 2008] and “nT phase angle measurements” [Haese 1999, White 2010] ]will be combined with data from density fluctuation diagnostics at ASDEX Upgrade to support a long-term program of physics research in turbulence and transport that will allow for more stringent testing and validation of gyrokinetic turbulent-transport codes. This work directly impacts the development of predictive transport models in the U.S. FES program, such as TGLF, developed by General Atomics, which are used to predict performance in ITER and other burning plasma devices as part of advancing the development of fusion energy sciences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2006

Non Radiation Hardened Microprocessors in Spaced Based Remote Sensing Systems

The CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) mission is a comprehensive suite of active and passive sensors including a 20Hz 230mj Nd:YAG lidar, a visible wavelength Earth-looking camera and an imaging infrared radiometer. CALIPSO flies in formation with the Earth Observing System Post-Meridian (EOS PM) train, provides continuous, near-simultaneous measurements and is a planned 3 year mission. CALIPSO was launched into a 98 degree sun synchronous Earth orbit in April of 2006 to study clouds and aerosols and acquires over 5 gigabytes of data every 24 hours. The ground track of one CALIPSO orbit as well as high and low intensity South Atlantic Anomaly outlines is shown. CALIPSO passes through the SAA several times each day. Spaced based remote sensing systems that include multiple instruments and/or instruments such as lidar generate large volumes of data and require robust real-time hardware and software mechanisms and high throughput processors. Due to onboard storage restrictions and telemetry downlink limitations these systems must pre-process and reduce the data before sending it to the ground. This onboard processing and realtime requirement load may mean that newer more powerful processors are needed even though acceptable radiation-hardened versions have not yet been released. CALIPSO's single board computer payload controller processor is actually a set of four (4) voting non-radiation hardened COTS Power PC 603r's built on a single width VME card by General Dynamics Advanced Information Systems (GDAIS). Significant radiation concerns for CALIPSO and other Low Earth Orbit (LEO) satellites include the South Atlantic Anomaly (SAA), the north and south poles and strong solar events. Over much of South America and extending into the South Atlantic Ocean the Van Allen radiation belts dip to just 200-800km and spacecraft entering this area are subjected to high energy protons and experience higher than normal Single Event Upset (SEU) and Single Event Latch-up (SEL) rates. Although less significant, spacecraft flying in the area around the poles experience similar upsets. Finally, powerful solar proton events in the range of 10MeV/10pfu to 100MeV/1pfu as are forecasted and tracked by NOAA's Space Environment Center in Colorado can result in Single Event Upset (SEU), Single Event Latch-up (SEL) and permanent failures such as Single Event Gate Rupture (SEGR) in some technologies. (Galactic Cosmic Rays (GCRs) are another source, especially for gate rupture) CALIPSO mitigates common radiation concerns in its data handling through the use of redundant processors, radiation-hardened Application Specific Integrated Circuits (ASIC), hardware-based Error Detection and Correction (EDAC), processor and memory scrubbing, redundant boot code and mirrored files. After presenting a system overview this paper will expand on each of these strategies. Where applicable, related on-orbit data collected since the CALIPSO initial boot on May 4, 2006 will be noted.

Decoursey, Robert J.↗

NASA NTRS2000

Optical Performance Modeling of FUSE Telescope Mirror

We describe the Metrology Data Processor (METDAT), the Optical Surface Analysis Code (OSAC), and their application to the image evaluation of the Far Ultraviolet Spectroscopic Explorer (FUSE) mirrors. The FUSE instrument - designed and developed by the Johns Hopkins University and launched in June 1999 is an astrophysics satellite which provides high resolution spectra (lambda/Delta(lambda) = 20,000 - 25,000) in the wavelength region from 90.5 to 118.7 nm The FUSE instrument is comprised of four co-aligned, normal incidence, off-axis parabolic mirrors, four Rowland circle spectrograph channels with holographic gratings, and delay line microchannel plate detectors. The OSAC code provides a comprehensive analysis of optical system performance, including the effects of optical surface misalignments, low spatial frequency deformations described by discrete polynomial terms, mid- and high-spatial frequency deformations (surface roughness), and diffraction due to the finite size of the aperture. Both normal incidence (traditionally infrared, visible, and near ultraviolet mirror systems) and grazing incidence (x-ray mirror systems) systems can be analyzed. The code also properly accounts for reflectance losses on the mirror surfaces. Low frequency surface errors are described in OSAC by using Zernike polynomials for normal incidence mirrors and Legendre-Fourier polynomials for grazing incidence mirrors. The scatter analysis of the mirror is based on scalar scatter theory. The program accepts simple autocovariance (ACV) function models or power spectral density (PSD) models derived from mirror surface metrology data as input to the scatter calculation. The end product of the program is a user-defined pixel array containing the system Point Spread Function (PSF). The METDAT routine is used in conjunction with the OSAC program. This code reads in laboratory metrology data in a normalized format. The code then fits the data using Zernike polynomials for normal incidence systems or Legendre-Fourier polynomials for grazing incidence systems. It removes low order terms from the metrology data, calculates statistical ACV or PSD functions, and fits these data to OSAC models for the scatter analysis. In this paper we briefly describe the laboratory image testing of FUSE spare mirror performed in the near and vacuum ultraviolet at John Hopkins University and OSAC modeling of the test setup performed at NASA/GSFC. The test setup is a double-pass configuration consisting of a Hg discharge source, the FUSE off-axis parabolic mirror under test, an autocollimating flat mirror, and a tomographic imaging detector. Two additional, small fold flats are used in the optical train to accommodate the light source and the detector. The modeling is based on Zernike fitting and PSD analysis of surface metrology data measured by both the mirror vendor (Tinsley) and JHU. The results of our models agree well with the laboratory imaging data, thus validating our theoretical model. Finally, we predict the imaging performance of FUSE mirrors in their flight configuration at far-ultraviolet wavelengths.

Saha, Timo T.↗

DOE OSTI2024

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

A Roadmap for a Lightning Modeling Grand Challenge

This document is a roadmap for building an interconnected model of the physical processes that produce a lightning discharge, and its observable optical and radio signals. We call this a Lightning Modeling Grand Challenge, recognizing that significant effort and coordination of human and financial resources is required to realize the capability. The roadmap serves to outline the coordination of resources necessary to enable stitching together existing knowledge and model components to make a lightning prediction, and to test these predictions with observations. Such a capability does not currently exist. The roadmap is motivated not only by a spirit of scientific inquiry, but by practical challenges faced by US Federal and societal stakeholders. Advancements in lightning observations have outpaced our tests of integrated understanding, leaving many stakeholders unsure how to design their missions to properly detect and discriminate lightning, and unsure how to apply the sometimes-disagreeing lightning signals from diverse instruments. The time is right to connect existing theories and models to support stakeholders in understanding the signals they observe, for needs as diverse as climate monitoring, national security, weather forecasting, public safety, and protection of natural and built environments. The roadmap’s two main technical sections describe the components of a linked physical model, followed by a description of models of lightning signals and sensors that are driven by outputs from the physical model. The goal is to predict the time-varying physical properties of lightning that are self-consistent with the thunderstorm’s structure and dynamics. These lightning signals then propagate through the storm, with realistic dispersion and attenuation, to receivers on the ground or in space. At a high level, the model begins with weather (cloud) model output, including explicit prediction of the electrification of cloud particles. The cloud’s electrical structure drives a model of lightning physics, from initiation, through channel development, and discharges along those channels. Key lightning parameters, such as the temperature and currents in the channel, and their space and time distribution, are then used to produce optical and electromagnetic signal sources that propagate to modeled receivers. This architecture therefore generates a dataset suitable for comparison to existing and envisioned observing systems. The need for additional measurements and field campaigns to support model development is described. In each model sub-component, inputs, outputs, uncertainties, evaluation methods, and next steps are summarized, interleaved with references to the scientific literature. Identifying boundaries between the model sub-components aids in segmenting an integrated, complex model into practical work packages and system sub-components, allowing a diverse team to contribute and maintain the system. We estimate that at least five years of effort and a $\$$10M initial investment is necessary to make a significant step forward. Mechanisms to facilitate community coordination, including annual workshops and open-source code repositories, are described.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

SU-E-T-150: Brachytherapy QA Employing a High Resolution Liquid Filled Ionisation Chamber Array: Initial Experience and Limitations

Purpose: Verifying a complex 3D brachytherapy dose distribution by measurement is non-trivial. Ideally a photon detector array should be independent of energy and angle, have high spatial resolution and be robust for routine clinical use. Methods: An iridium-192 source was used. A PMMA jig was constructed comprising an outer slab and a central insert with eight milled channels for 1.33mm (outer diameter) steel needles, see figure. All calculations were performed using an empty CT study reconstructing eight virtual needles (QA-CT), using the v2 source model (Elekta AG, Sweden). A high resolution liquid filled ionisation chamber array SRS1000, together with Verisoft software v6.0 (PTW Freiburg, Germany), was used to perform measurements of plans of increasing complexity to evaluate its suitability for device- and patient-specific QA. The dimension of backscatter material was investigated. The patient plan dwell time distribution was entered manually into the QA-CT and the dose distribution was calculated. Results: Our measurements indicate that the array is independent of energy and angle. The resulting measured dose values are linearly interpolated to 2025 values. Shifts of 1mm of the entire needle are readily detectable. Individual dwell position shifts (2.5mm) are also readily measurable. Moreover a dwell time increase of 1 second both in the edge and central region are detectable. Conclusion: The high resolution SRS1000 array is a powerful instrument for brachytherapy QA enabling 977 simultaneous measurements to be performed. Our measurements suggest 60mm of RW3 backscatter material upstream and downstream are sufficient. Local percentage difference analysis is useful for device based QA, normalized relative percentage difference is arguably better for patient specific QA. Automated transfer of patient plan dwell time distribution to the QA plan is required to enable a comprehensive patient QA study to be performed. Moreover the described measurement technique can be used for commissioning new brachytherapy applicators.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-IeP4-06: Feasibility of External Beam Treatment Field Verification Using Cherenkov Imaging

Purpose: Cherenkov light emission has been shown to correlate with ionizing radiation (IR) dose delivery in solid tissue. In order to properly correlate Cherenkov light images with real time dose delivery in a patient, we must account for geometric and intensity distortions arising from observation angle, as well as the effect of monitor units (MU) and field size on Cherenkov light emission. To test the feasibility of treatment field verification, we first focused on Cherenkov light emission efficiency based on MU and known field size (FS). Methods: Cherenkov light emission was captured using a PI-MAX4 intensified charge coupled device(ICCD) system (Princeton Instruments), positioned at a fixed angle of 40° relative to the beam central axis. A Varian TrueBeam linear accelerator (linac) was operated at 6MV and 600MU/min to deliver an Anterior-Posterior beam to a 5cm thick block phantom positioned at 100cm Source-to-Surface-Distance(SSD). FS of 10×10, 5×5, and 2×2cm{sup 2} were used. Before beam delivery projected light field images were acquired, ensuring that geometric distortions were consistent when measuring Cherenkov field discrepancies. Cherenkov image acquisition was triggered by linac target current. 500 frames were acquired for each FS. Composite images were created through summation of frames and background subtraction. MU per image was calculated based on linac pulse delay of 2.8ms. Cherenkov and projected light FS were evaluated using ImageJ software. Results: Mean Cherenkov FS discrepancies compared to light field were <0.5cm for 5.6, 2.8, and 8.6 MU for 10×10, 5×5, and 2×2cm{sup 2} FS, respectably. Discrepancies were reduced with increasing field size and MU. We predict a minimum of 100 frames is needed for reliable confirmation of delivered FS. Conclusion: Current discrepancies in Cherenkov field sizes are within a usable range to confirm treatment delivery in standard and respiratory gated clinical scenarios at MU levels appropriate to standard MLC position segments.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

SU-E-T-47: A Monte Carlo Model of a Spot Scanning Proton Beam Based On a Synchrotron Proton Therapy Accelerator

Purpose: To build the model of a spot scanning proton beam for the dose calculation of a synchrotron proton therapy accelerator, which is capable of accelerating protons from 50 up to 221 MeV. Methods: The spot scanning beam nozzle is modeled using TOPAS code, a simulation tool based on Geant4.9.6. The model contained a beam pipe vacuum window, a beam profile monitor, a drift chamber, two plane-parallel ionization chambers, and a spot-position monitor consisted of a multiwire ionization chamber. A water phantom is located with its upstream surface at the isocenter plane. The initial proton beam energy and anglar deflection are modeled using a Gaussian distribution with FWHM (Full Widths at Half Maximum) deponding on its beam energy. The phase space file (PSF) on a virtual surface located at the center between the two magnets is recorded. PSF is used to analyze the pencil beam features and offset the pencil beam position. The source model parameters are verificated by fitting the simulated Result to the measurement. Results: The simulated percentage depth dose (PDD) and lateral profiles of scanning pencil beams of various incident proton energies are verificated to the measurement. Generally the distance to agreement (DTA) of Bragg peaks is less than 0.2cm. The FWHM of Gaussian anglar distribution was adjusted to fit the lateral profile difference between the simulation and the measurement to less than 2∼3cm. Conclusion: A Monte Carlo model of a spot scanning proton beam was bullt based on a synchrotron proton therapy accelerator. This scanning pencil beam model will be as a block to build the broad proton beam as a proton TPS dose verification tool.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Development of a true single line of sight 3D hot-spot imaging for the National Ignition Facility

High resolution 3D self-emission x-ray imaging during inertial confinement fusion capsule implosion enables the measurement of the shape of the hotspot. While current 3D imaging capabilities use multiple lines of sight to perform image reconstruction, it would be highly desirable to use only one line of sight, as this would significantly reduce the number of windows in the target hohlraum and decrease their impacts on implosion symmetry. Such a goal is achievable using the zone plate coded imaging technique developed by N. M. Ceglio [Proc. SPIE 0106, 55-62 (1977)]. It consists of fielding a visible light Fresnel zone plate on shot to record a shadowgraph. Here, the image can then be reconstructed either by printing the shadowgraph on a transparent film and shining a suitable wavelength light through it or by numerical reconstruction. A new approach using numerical reconstruction is presented, and it relaxes the constraint by an order of magnitude on the optic design, thus enabling an easier fabrication process, as it allows a scaling-up of the optic dimensions. The design, fabrication process, and testing with an x-ray source of a prototype is presented. The reconstruction of an ~14.5 × 17 μm 2 broadband x-ray source was successful and shows that the performances are in line with expectation with an at least 5 mm axial resolution.

47 OTHER INSTRUMENTATION↗