Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source 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 1,639 records · Page 91

Nuclear and in-source laser spectroscopy with the ISAC yield station

A new decay station has been built for the ISAC facility at TRIUMF for the rapid and reliable characterization of radioactive ion beam (RIB) compositions and intensities with the capability of simultaneously collecting α, β, and γ decay data from RIB with intensities between a few and ≈10{sup 11} ions per second. It features user-friendly control, data acquisition, and analysis software. The analysis of individual decay time structures allows the unambiguous assignment of α and γ lines even with substantial isobaric contamination present. The capability for accurate half-life measurements is demonstrated with the example of {sup 46}K. The coupling of the yield station to the laser ion source, TRILIS, allows the correlation of radiometric data with automated laser frequency scans. First results of in-source laser spectroscopy measurements on astatine are discussed.

43 PARTICLE ACCELERATORS↗

eqtools: Modular, extensible, open-source, cross-machine Python tools for working with magnetic equilibria

As plasma physics research for fusion energy transitions to an increasing emphasis on cross-machine collaboration and numerical simulation, it becomes increasingly important that portable tools be developed to enable data from diverse sources to be analyzed in a consistent manner. This paper presents eqtools, a modular, extensible, open-source toolkit implemented in the Python programming language for handling magnetic equilibria and associated data from tokamaks. eqtools provides a single interface for working with magnetic equilibrium data, both for handling derived quantities and mapping between coordinate systems, extensible to function with data from different experiments, data formats, and magnetic reconstruction codes, replacing the diverse, non-portable solutions currently in use. Moreover, while the open-source Python programming language offers a number of advantages as a scripting language for research purposes, the lack of basic tokamak-specific functionality has impeded the adoption of the language for regular use. Implementing equilibrium-mapping tools in Python removes a substantial barrier to new development in and porting legacy code into Python. In this paper, we introduce the design of the eqtools package and detail the workflow for usage and expansion to additional devices. Additionally, the implementation of a novel three-dimensional spline solution (in two spatial dimensions and in time) is also detailed. Finally, verification and benchmarking for accuracy and speed against existing tools are detailed. Wider deployment of these tools will enable efficient sharing of data and software between institutions and machines as well as self-consistent analysis of the shared data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fermi-Hubbard Transport Solver v. 1.0

SAND2021-1511 O. The code takes a description of a system comprised of two leads that are sources and sinks of electrons, a set of abstract sites that those electrons can occupy, and their connectivity. It then predicts the current passing between the two leads through the set of abstract sites. There are numerical parameters that define the rate at which electrons tunnel between sites and another parameter that determines the energy of having two electrons on the same site. When the ratio of these energies and biases are set just right, they can be used to predict the existence of Coulomb blockade in which no current flows because the energetic penalty would be too high. That is, the code predicts a current and the absence of a current which is predictive of Coulomb blockade. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Brickson, MitchellI↗

ARPA-E Grid Optimization (GO) Competition Challenge 1

The ARPA-E Grid Optimization (GO) Competition Challenge 1, from 2018 to 2019, focused on the basic Security Constrained AC Optimal Power Flow problem (SCOPF) for a single time period. The Challenge utilized sets of unique datasets generated by the ARPA-E GRID DATA program. Each dataset consisted of a collection of power system network models of different sizes with associated operating scenarios (snapshots in time defining instantaneous power demand, renewable generation, generator and line availability, etc.). The datasets were of two types: Real-Time, which included starting-point information, and Online, which did not. Week-Ahead data is also provided for some cases but was not used in the Competition. Although most datasets were synthetic and generated by GRIDDATA, a few came from industry and were only used in the Final Event. All synthetic Input Data and Team Results for the GO Competition Challenge 1 for the Sandbox, Trial Events 1 to 3, and the Final Event along with problem, format, scoring and rules descriptions are available here. Data for industry scenarios will not be made public. Challenge 1, a minimization problem, required two computational steps. Solver 1 or Code 1 solved the base SCOPF problem under a strict wall clock time limit, as would be the case in industry, and reported the base case operating point as output, which was used to compute the Objective Function value that was used as the scenario score. The feasibility of the solution was provided by the Solver 2 or Code 2, which solves the power flow problem for all contingencies based on the results from Solver 1. This is not normally done in industry, so the time limits were relaxed. In fact, there were no time limits for Trial Event 1. This proved to be a mistake, with some codes running for more than 90 hours, and a time limit of 2 seconds per contingency was imposed for all other events. Entrants were free to use their own Solver 2 or use an open-source version provided by the Competition. Containers, such as Docker, were considered to improve the portability of codes, but none that could reliably support a multi-node parallel computing environment, e.g., MPI, could be found. For more information on the competition and challenge see the "GO Competition Challenge 1 Information" and "GO Competition Challenge 1 Additional Information" resources below.

ACOPF↗

Coded-mask-based wavefront sensing technique for APS nanofocusing beamline diagnostics

Here, we extend our recently developed coded-mask wavefront sensing technique to enable single-shot measurements of nanofocused x-ray beams. This method accurately reconstructs the focal beam profile by backpropagating the wavefront measured downstream of the beam focus. To validate its performance, we benchmarked it against the conventional fluorescence wire scan method, successfully measuring ∼120 nm focal spots at the 28-ID-B beamline of the Advanced Photon Source using a polymeric compound refractive lens. The results highlight the effectiveness of coded-mask wavefront sensing for high-precision beam profiling and its application as a real-time wavefront monitoring tool.

Shi, Xianbo [Argonne National Laboratory (ANL), Ar↗

SU-C-201-03: Coded Aperture Gamma-Ray Imaging Using Pixelated Semiconductor Detectors

Purpose: Improved localization of gamma-ray emissions from radiotracers is essential to the progress of nuclear medicine. Polaris is a portable, room-temperature operated gamma-ray imaging spectrometer composed of two 3×3 arrays of thick CdZnTe (CZT) detectors, which detect gammas between 30keV and 3MeV with energy resolution of <1% FWHM at 662keV. Compton imaging is used to map out source distributions in 4-pi space; however, is only effective above 300keV where Compton scatter is dominant. This work extends imaging to photoelectric energies (<300keV) using coded aperture imaging (CAI), which is essential for localization of Tc-99m (140keV). Methods: CAI, similar to the pinhole camera, relies on an attenuating mask, with open/closed elements, placed between the source and position-sensitive detectors. Partial attenuation of the source results in a “shadow” or count distribution that closely matches a portion of the mask pattern. Ideally, each source direction corresponds to a unique count distribution. Using backprojection reconstruction, the source direction is determined within the field of view. The knowledge of 3D position of interaction results in improved image quality. Results: Using a single array of detectors, a coded aperture mask, and multiple Co-57 (122keV) point sources, image reconstruction is performed in real-time, on an event-by-event basis, resulting in images with an angular resolution of ∼6 degrees. Although material nonuniformities contribute to image degradation, the superposition of images from individual detectors results in improved SNR. CAI was integrated with Compton imaging for a seamless transition between energy regimes. Conclusion: For the first time, CAI has been applied to thick, 3D position sensitive CZT detectors. Real-time, combined CAI and Compton imaging is performed using two 3×3 detector arrays, resulting in a source distribution in space. This system has been commercialized by H3D, Inc. and is being acquired for various applications worldwide, including proton therapy imaging R&D.

60 APPLIED LIFE SCIENCES↗

Radionuclide Inventory and Decay Heat Quantification Methodology for Severe Accident Simulations.

The MELCOR and MACCS computer codes require inputs for radionuclide inventory and decay heat in order to simulate severe accident phenomenology, source term, and consequences. Therefore, Sandia National Laboratories ( SNL ) has developed accurate and automated methods using SCALE 6 in conjunction with automation scripts to directly generate MELCOR and MACCS input records for radionuclide inventories and decay heating. Consistent information between the two codes is essential for realistic modeling of severe nuclear accidents it is the radionuclide inventory and decay power that are the principal safety problems of any severe accident that assumes successful reactor shutdown. SNL plans to conduct best-estimate calculations and uncertainty analyses using MELCOR and MACCS for several reactors, including units 1 - 3 of Fukushima Daiichi, which entail accurate and technically-scrutable information for the input models. The SNL methodology for generating radionuclide and decay power inputs directly from SCALE6 calculations is mostly presented in the context of supporting Fukushima modeling efforts. Furthermore, certain historical assumptions used by MELCOR and MACCS in the abstraction of radionuclide classes and their properties are reviewed, and these assumptions are compared to modern code predictions by SCALE6.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Considerations for Developing a Regulatory Roadmap for Distributed Energy Resource (DER) Integration in Colombia

The USAID-NREL Partnership in Colombia, in conjunction with SURE and USEA, selected four action plan teams from the Young Professionals Leadership Program to receive tailored NREL technical assistance to support action plan implementation. The four plans selected represent various aspects associated with planning for the efficient integration of DERs, including electric mobility, residential and commercial energy applications, distributed generation modeling, and regulatory considerations. In response to national grid modernization, decentralization, digitalization, and electrification trends, CREG is looking to update Colombia's electric distribution code, with an emphasis on the integration of electric vehicles and charging stations, but also energy storage systems, non-conventional sources of renewable energy, and other distributed energy resources (DERs). CREG aims to establish the general steps that must be carried out to update the distribution code for the integration of DERs, with a view to the design and development of new markets. The objective of the TA provided to CREG was to support them in identifying and outlining the general processes and approaches to consider when updating Colombia's distribution code for the integration of electric vehicles, distributed renewable energy systems, and energy storage systems, with a view to the design and development of new markets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SU-E-T-29: A Web Application for GPU-Based Monte Carlo IMRT/VMAT QA with Delivered Dose Verification

Purpose: To enable an existing web application for GPU-based Monte Carlo (MC) 3D dosimetry quality assurance (QA) to compute “delivered dose” from linac logfile data. Methods: We added significant features to an IMRT/VMAT QA web application which is based on existing technologies (HTML5, Python, and Django). This tool interfaces with python, c-code libraries, and command line-based GPU applications to perform a MC-based IMRT/VMAT QA. The web app automates many complicated aspects of interfacing clinical DICOM and logfile data with cutting-edge GPU software to run a MC dose calculation. The resultant web app is powerful, easy to use, and is able to re-compute both plan dose (from DICOM data) and delivered dose (from logfile data). Both dynalog and trajectorylog file formats are supported. Users upload zipped DICOM RP, CT, and RD data and set the expected statistic uncertainty for the MC dose calculation. A 3D gamma index map, 3D dose distribution, gamma histogram, dosimetric statistics, and DVH curves are displayed to the user. Additional the user may upload the delivery logfile data from the linac to compute a 'delivered dose' calculation and corresponding gamma tests. A comprehensive PDF QA report summarizing the results can also be downloaded. Results: We successfully improved a web app for a GPU-based QA tool that consists of logfile parcing, fluence map generation, CT image processing, GPU based MC dose calculation, gamma index calculation, and DVH calculation. The result is an IMRT and VMAT QA tool that conducts an independent dose calculation for a given treatment plan and delivery log file. The system takes both DICOM data and logfile data to compute plan dose and delivered dose respectively. Conclusion: We sucessfully improved a GPU-based MC QA tool to allow for logfile dose calculation. The high efficiency and accessibility will greatly facilitate IMRT and VMAT QA.

07 ISOTOPE AND RADIATION SOURCES↗

Systems for the Intermodal Routing of Spent Nuclear Fuel

The safe and secure movement of spent nuclear fuel from shutdown and active reactor facilities to intermediate or long term storage sites may, in some instances, require the use of several modes of transportation to accomplish the move. To that end, a fully operable multi-modal routing system is being developed within Oak Ridge National Laboratory s (ORNL) WebTRAGIS (Transportation Routing Analysis Geographic Information System). This study aims to provide an overview of multi-modal routing, the existing state of the TRAGIS networks, the source data needs, and the requirements for developing structural relationships between various modes to create a suitable system for modeling the transport of spent nuclear fuel via a multimodal network. Modern transportation systems are comprised of interconnected, yet separate, modal networks. Efficient transportation networks rely upon the smooth transfer of cargoes at junction points that serve as connectors between modes. A key logistical impediment to the shipment of spent nuclear fuel is the absence of identified or designated transfer locations between transport modes. Understanding the potential network impacts on intermodal transportation of spent nuclear fuel is vital for planning transportation routes from origin to destination. By identifying key locations where modes intersect, routing decisions can be made to prioritize cost savings, optimize transport times and minimize potential risks to the population and environment. In order to facilitate such a process, ORNL began the development of a base intermodal network and associated routing code. The network was developed using previous intermodal networks and information from publicly available data sources to construct a database of potential intermodal transfer locations with likely capability to handle spent nuclear fuel casks. The coding development focused on modifying the existing WebTRAGIS routing code to accommodate intermodal transfers and the selection of prioritization constraints and modifiers to determine route selection. The limitations of the current model and future directions for development are discussed, including the current state of information on possible intermodal transfer locations for spent fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physics-Based Hazard Assessment for Critical Structures Near Large Earthquake Sources

We argue that for critical structures near large earthquake sources: (1) the ergodic assumption, recent history, and simplified descriptions of the hazard are not appropriate to rely on for earthquake ground motion prediction and can lead to a mis-estimation of the hazard and risk to structures; (2) a physics-based approach can address these issues; (3) a physics-based source model must be provided to generate realistic phasing effects from finite rupture and model near-source ground motion correctly; (4) wave propagations and site response should be site specific; (5) a much wider search of possible sources of ground motion can be achieved computationally with a physics-based approach; (6) unless one utilizes a physics-based approach, the hazard and risk to structures has unknown uncertainties; (7) uncertainties can be reduced with a physics-based approach, but not with an ergodic approach; (8) computational power and computer codes have advanced to the point that risk to structures can be calculated directly from source and site-specific ground motions. Spanning the variability of potential ground motion in a predictive situation is especially difficult for near-source areas, but that is the distance at which the hazard is the greatest. The basis of a “physical-based” approach is ground-motion syntheses derived from physics and an understanding of the earthquake process. This is an overview paper and results from previous studies are used to make the case for these conclusions. Our premise is that 50 years of strong motion records is insufficient to capture all possible ranges of site and propagation path conditions, rupture processes, and spatial geometric relationships between source and site. Predicting future earthquake scenarios is necessary; models that have little or no physical basis but have been tested and adjusted to fit available observations can only “predict” what happened in the past, which should be considered description as opposed to prediction. Finally, we have developed a methodology for synthesizing physics-based broadband ground motion that incorporates the effects of realistic earthquake rupture along specific faults and the actual geology between the source and site.

58 GEOSCIENCES↗

Compiler verification meets cross-language linking via data abstraction

Many real programs are written in multiple different programming languages, and supporting this pattern creates challenges for formal compiler verification. We describe our Coq verification of a compiler for a high-level language, such that the compiler correctness theorem allows us to derive partial-correctness Hoare-logic theorems for programs built by linking the assembly code output by our compiler and assembly code produced by other means. Our compiler supports such tricky features as storable cross-language function pointers, without giving up the usual benefits of being able to verify different compiler phases (including, in our case, two classic optimizations) independently. Here, the key technical innovation is a mixed operational and axiomatic semantics for the source language, with a built-in notion of abstract data types, such that compiled code interfaces with other languages only through axiomatically specified methods that mutate encapsulated private data, represented in whatever formats are most natural for those languages.

97 MATHEMATICS AND COMPUTING↗

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗

Application of Monte Carlo Methods to Perform Uncertainty and Sensitivity Analysis on Inverse Water-Rock Reactions with NETPATH

Methods were developed to quantify uncertainty and sensitivity for NETPATH inverse water-rock reaction models and to calculate dissolved inorganic carbon, carbon-14 groundwater travel times. The NETPATH models calculate upgradient groundwater mixing fractions that produce the downgradient target water chemistry along with amounts of mineral phases that are either precipitated or dissolved. Carbon-14 groundwater travel times are calculated based on the upgradient source-water fractions, carbonate mineral phase changes, and isotopic fractionation. Custom scripts and statistical code were developed for this study to facilitate modifying input parameters, running the NETPATH simulations, extracting relevant output, postprocessing the results, and producing graphs and summaries. The scripts read userspecified values for each constituent’s coefficient of variation, distribution, sensitivity parameter, maximum dissolution or precipitation amounts, and number of Monte Carlo simulations. Monte Carlo methods for analysis of parametric uncertainty assign a distribution to each uncertain variable, sample from those distributions, and evaluate the ensemble output. The uncertainty in input affected the variability of outputs, namely source-water mixing, phase dissolution and precipitation amounts, and carbon-14 travel time. Although NETPATH may provide models that satisfy the constraints, it is up to the geochemist to determine whether the results are geochemically reasonable. Two example water-rock reaction models from previous geochemical reports were considered in this study. Sensitivity analysis was also conducted to evaluate the change in output caused by a small change in input, one constituent at a time. Results were standardized to allow for sensitivity comparisons across all inputs, which results in a representative value for each scenario. The approach yielded insight into the uncertainty in water-rock reactions and travel times. For example, there was little variation in source-water fraction between the deterministic and Monte Carlo approaches, and therefore, little variation in travel times between approaches. Sensitivity analysis proved very useful for identifying the most important input constraints (dissolved-ion concentrations), which can reveal the variables that have the most influence on source-water fractions and carbon-14 travel times. Once these variables are determined, more focused effort can be applied to determining the proper distribution for each constraint. Second, Monte Carlo results for water-rock reaction modeling showed discrete and nonunique results. The NETPATH models provide the solutions that satisfy the constraints of upgradient and downgradient water chemistry. There can exist multiple, discrete solutions for any scenario and these discrete solutions cause grouping of results. As a result, the variability in output may not easily be represented by a single distribution or a mean and variance and care should be taken in the interpretation and reporting of results.

phase dissolution and precipitation amounts↗

Modeling study of a proposed field calibration source using K-40 and high-Z targets for sodium iodide detectors

Calibration sources based on the primordial isotope potassium-40 ( 40 K) have reduced controls on the source’s activity due to its terrestrial ubiquity and very low specific activity. Potassium–40’s beta emissions and 1,460.8 keV gamma ray can be used to induce K-shell fluorescence x rays in high-Z metals between 60 and 80 keV. A gamma ray calibration source that uses potassium chloride salt and a high-Z metal to create a two-point calibration for a sodium iodide field gamma spectroscopy instrument is thus proposed. The calibration source was designed in collaboration with the Sandia National Laboratory using the Monte Carlo N-Particle eXtended (MCNPX) transport code. Two methods of x-ray production were explored. First, a thin high-Z layer (HZL) was interposed between the detector and the potassium chloride-urethane source matrix. Second, bismuth metal powder was homogeneously mixed with a urethane binding agent to form a potassium chloride-bismuth matrix (KBM). The bismuth-based source was selected as the development model because it is inexpensive, nontoxic, and outperforms the high-Z layer method in simulation. As a result, based on the MCNPX studies, sealing a mixture of bismuth powder and potassium chloride into a thin plastic case could provide a light, inexpensive field calibration source.

40K↗

Exergy Analysis of Heating Options for Buildings in the Pacific Northwest

A framework for quantifying environmental impact in future energy codes is explored based on exergy analysis. An example exergy analysis is performed that quantifies the environmental impact of electric resistance heating with exergetic efficiency of only 13% compared to 30% and 40% for air source and ground source heat pumps. This provides motivation for removing electric resistance as a space heating option in future energy codes in the Northwest.

Dillon, Heather↗

SU-F-T-524: Investigation of the Dosimertric Benefits of Interchangeable Source Size of a Novel Rotating Gamma System

Purpose: Tremendous technological developments were made for conformal therapy techniques with linear accelerators, while less attention was paid to cobalt-60 units. The aim of the current study is to explore the dosimetric benefits of a novel rotating gamma ray system enhanced with interchangeable source sizes and multi-leaf collimator (MLC). Material and Methods: CybeRT is a novel rotating gamma ray machine with a ring gantry that ensures an iso-center accuracy of less than 0.3 mm. The new machine has a 70cm source axial distance allowing for improved penumbra compared to conventional machines. MCBEAM was used to simulate Cobalt-60 beams from the CybeRT head, while the MCPLAN code was used for modeling the MLC and for phantom/patient dose calculation. The CybeRT collimation will incorporate a system allowing for interchanging source sizes. In this work we have created phase space files for 1cm and 2cm source sizes. Evaluation of the system was done by comparing CybeRT beams with the 6MV beams in a water phantom and in patient geometry. Treatment plans were compared based on isodose distributions and dose volume histograms. Results: Profiles for the 1cm source were comparable to that from 6MV in the order of 6mm for 10×10 cm{sup 2} field size at the depth of maximum dose. This could ascribe to Cobalt-60 beams producing lowerenergy secondary electrons. Although, the 2cm source have a larger penumbra however it could be still used for large targets with proportionally increased dose rate. For large lung targets, the difference between cobalt and 6MV plans is clinically insignificant. Our preliminary results showed that interchanging source sizes will allow cobalt beams for volumetric arc therapy of both small lesions and large tumors. Conclusion: The CybeRT system will be a cost effective machine capable of performing advanced radiation therapy treatments of both small tumors and large target volumes.

60 APPLIED LIFE SCIENCES↗

Update to the Microcontroller Benchmark for Radiation Testing

LANL developed a benchmark of software code for radiation testing of microprocessors several years ago, and it was published under an open-source license on GitHub. Publishing the software is necessary for other researchers to adopt and implement this benchmark for radiation testing of other microprocessors to standardize test practices so that test data can be compared across different microprocessors. The original codes have been used several times by other organizations to test a wide range of microcontrollers and microprocessors. After several years of research, LANL is ready to update the benchmark. Changes include: 1. Addition of new codes that allow common software codes to be tested, 2. Addition of new codes that instrument more microprocessor circuitry, 3. Addition of input patterns that allow for a more compressive understanding of how the memory layout affects the sensitivity to radiation-induced faults and better use of automated test pattern generation standards, and 4. Modification of current codes for faster and more resilient detection, reporting and correction of radiation-induced faults. These codes have been tested by LANL researchers over the last few years, which has been published in the open literature. As the code base for the new benchmarks are stable, it is time to release the update to the GitHub repository, where the original codes were released.

Quinn, Heather↗