Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Python codes”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Search for New Isotope Production Pathways

The isotope group at Jefferson Lab is carrying out R&D for producing medically interesting radioisotopes, especially those with theranostic (therapeutic and diagnostic) attributes. Here the search for viable production mechanisms has been expanded to multi-step reactions where a daughter is produced from the target and decays into a medically interesting granddaughter radioisotope. It is difficult to find efficient production routes when investigating both the initial excitation reaction as well as the decay routes leading to medically interesting isotopes. The overall goal of this project is to create a structured code in Python to find these decay routes by automatically exploring the large number of isotopes and their possible decay modes. The program structure is described, and preliminary results are presented.

Dabill, L.↗
DOE OSTI2022

Development and Commercialization of an IDAES-Based Power Plant Performance Monitoring and Optimization System

DOE and NETL have created an advanced, open-source computational platform through the Institute for the Design of Advanced Energy Systems (IDAES). The IDAES platform is a very extensive modeling environment that can be used for a broad range of power plant and process applications. MapEx Software is developing and commercializing a software application that makes it easier to set-up and run IDAES-based analyses. The MapEx-developed software application will replace the need for custom Python language coding with a user-friendly graphical user interface where the user can construct a flowsheet diagram of the IDAES model by inserting icons representing the plant equipment onto the screen. This application will make the implementation of modeling and optimization of existing fossil-fired power plants more straight-forward and less time-consuming. The effort focuses on performance monitoring and optimization of plant operations for the existing coal-fired power plant fleet but is built on a structure that allows expansion into the broad range of applications where IDAES methods may be applied.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented here is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameterfile driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Advanced Precipitation and Boundary Layer Data Products Derived from ARM Radar Wind Profilers

This research project was successful in delivering on four main objectives. First, software was developed to accurately calculate 915-MHz radar wind profiler (RWP) spectrum moments from the recorded Doppler velocity power spectra. Second, software was developed to calibrate the RWP reflectivity factor using collocated surface disdrometer observations. Third, the Python processing code was documented and given to the ARM Infrastructure to produce ARM ‘b level’ calibrated RWP products. Fourth, calibrated RWP products were uploaded to the ARM Archive as PI Products for 10 years of SGP RWP observations and for GoAmazon and TRACER field campaign RWP observations. In addition to working with RWP observations, this research project also worked with KAZR observations to distinguish insects from boundary layer clouds to help improve the ARSCL cloud mask product. The PI worked with senior and early career ARM funded scientists at BNL exploring how to include calibrated RWP moments into future versions of the ARSCL product.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Advancing 1D thermo-hydraulic tools for large cryogenic facilities

The Cryogenic Division at Fermilab develops large-scale cryogenic systems for particle accelerators and superconducting test facilities. To support design and diagnostics, a Python-based code was created to calculate pressure drops in components such as valves and pipes. This paper presents recent enhancements to the code, including new heat transfer functions that improve the accuracy of thermal and hydraulic predictions. The first application models pressure and temperature changes in PIP-II relief pipes, aiding pipe sizing and protecting relief valves. The second example analyzes heat load evolution in a pipe carrying sub-atmospheric helium, helping interpret temperature sensor data and understand gas return behavior to cold compressors. These improvements significantly expand the tool’s capabilities, offering a practical resource for designing and troubleshooting cryogenic systems under dynamic thermal and flow conditions.

Beckwith, Rosalyn [Fermilab]↗
DOE OSTI2025

Computational Prediction of Infrasound Arrival Times and Directions from Stationary and Moving Impulsive Sources

This report addresses the need to predict infrasound signal arrival times and back azimuths at monitoring stations, enabling more focused and efficient searches within recorded waveform data. The primary challenge is estimating expected signal arrival windows for stationary and moving acoustic sources, such as chemical explosions, volcanic eruptions, meteoroids, and spacecraft re-entry events. To address this challenge, a reproducible methodology is described that uses simplified propagation speeds for boundary layer, tropospheric, stratospheric, and thermospheric atmospheric waveguides. While the Python source code itself is not freely available, this document provides detailed, step-by-step instructions, and equations enabling users to replicate and adapt the method independently. The method reliably predicts signal arrival intervals and back azimuths, thereby supporting rapid detection and accurate interpretation of infrasound events. Results demonstrate that this method effectively identifies plausible signal arrival intervals and directions, facilitating faster event detection and more reliable interpretation. This methodology directly supports atmospheric monitoring, planetary defense, and forensic analysis of explosive atmospheric events.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

The Impact of Operated Window Shading on Visual Comfort, Non-Visual Health, and Energy Demand from Electric Lighting

This paper introduces a simulation-based workflow to assess annual circadian lighting performance, glare risk and energy demand for an array of seated view positions under operated shading and electric lighting conditions. This workflow uses a combination of 9-band and 81-band simulations to represent daylight and electric lighting respectively. Annual daylight-driven calculations rely on a python-based code provided in the LARK simulation engine which was used to run 180-degree high dynamic range (HDR) renderings with 9-bands of spectral resolution. The results reveal the impact of shading systems on equivalent melanopic lux (EML) and daylight glare probability (DGP) over time and the energy demand (kWh) from electric lighting systems to supplement eye-level recommendations from the WELL v2 Building Standard.

Rockcastle, Siobhan↗
DOE OSTI2025

Data for "A study of resistive peeling-ballooning modes across spherical tokamaks"

The data in this repository correspond to the figures of the paper "A study of resistive peeling-ballooning modes across spherical tokamaks" published in the journal Plasma Physics and Controlled Fusion. The figures are created based on simulations with the M3D-C1 and NIMROD codes using Python post-processing tools.

edge localized modes↗
DOE OSTI2020

Constraining the Halo Mass of Damped Lyα Absorption Systems (DLAs) at z = 2–3.5 Using the Quasar-CMB Lensing Cross-correlation

We study the cross-correlation of damped Lyα systems (DLAs) and their background quasars, using the most updated DLA catalog and the Planck 2018 CMB lensing convergence field. Our measurement suggests that the DLA bias b DLA is smaller than 3.1, corresponding to log (M/M ⊙ h -1 ) ⩽ 12.3 at a confidence of 90%. These constraints are broadly consistent with Alonso et al. and previous measurements by cross-correlation between DLAs and the Lyα forest (e.g., Font-Ribera et al.; Prez-Rfols et al.). Further, our results demonstrate the potential of obtaining a more precise measurement of the halo mass of the high-redshift sources using next generation CMB experiments with a higher angular resolution. The python-based codes and data products of our analysis are available at https://github.com/LittleLin1999/CMB-lensingxDLA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Calibrating radar wind profiler reflectivity factor using surface disdrometer observations

Abstract. This study uses surface disdrometer reflectivity factor estimates to calibrate the vertical and off-vertical pointing radar beams produced by an ultra high frequency (UHF) band radar wind profiler (RWP) deployed at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) Central Facility in northern Oklahoma from April 2011 through July 2019. The methodology consists of five steps. First, the recorded Doppler velocity power spectra are adjusted to account for Nyquist velocity aliasing and coherent integration filtering effects. Second, the spectrum moments are calculated. The third step increases the signal-to-noise ratio (SNR) due to inflated noise power estimates during convective rain events that cause SNR to be biased low. The fourth step determines the RWP calibration constant for one radar beam (called the “reference” beam) by comparing uncalibrated RWP reflectivity factors at 500 m above the ground to 1 min resolution surface disdrometer reflectivity factors. The last step uses the calibrated reference beam reflectivity factor to calibrate the other radar beams during precipitation. There are two key findings. The RWP sensitivity decreased by approximately 3 to 4 dB yr−1 as the hardware aged. This drift was slow enough that the reference calibration constant can be estimated over 3-month intervals using episodic rain events. The calibrated moments are available on the DOE ARM data archive, and the Python processing code is available on public repositories.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Towards Plume Impingement Modeling in Space Environments

After 30 years human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to act as a home base for Lunar exploration. A large space station must be assembled while in orbit, using a “piecemeal” approach. In this context, it means that the different modules will arrive at different times and attach to what is already in service. The ISS provides an excellent example of this approach, and the proposed Lunar Gateway will undergo a similar assembly process. This assembly is achieved via “docking” maneuvers between modules, which are made possible by sequential firings of the onboard reaction control system (RCS) thrusters. They work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking approach. The issue is that the gas from these thrusters' forms flow structures described as “plumes” and can impinge onto the outer surfaces of the space station, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion. All mechanisms that can damage the space station and must be avoided. Both permanent and visiting modules will have these RCS thruster exhaust impingement problems. Accurately and efficiently modeling these plume is an involved multi-physics calculation but also an important tool when designing the control algorithms of approaching modules. This poster presents progress towards this simulation on two fronts. First is the verification of OpenFOAM for rarefied plume impingement calculations by direct comparison to published DAC cases. Second is the estimation of plume impingement strikes over the time scale of an entire docking event. This is done by using a simple plume source flow model and a prescribed motion visualizer—coded in Python. Together these tools help push NASA's capabilities for simulating these plume impingement effects.

Rarefied Flows↗
NASA NTRSDate not supplied

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.

Aeroacoustics↗
DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Low-Temperature Geothermal Play Fairway Analysis for the Denver Basin

This project is part of a national initiative to showcase the benefits of incorporating low-temperature geothermal resource assessment into the deployment of geothermal heating, combined heat and power (CHP), and geothermal direct-use (GDU) technologies. The initiative was established to accelerate the country's decarbonization efforts by identifying potential for low-temperature geothermal resource utilization (<150 degrees C, e.g., CHP and GDU) in selected sedimentary basins with numerous population centers. The play fairway analysis (PFA) methodologies in this study were adapted from previous PFA investigations of sedimentary basin geothermal play types (SBGPTs) that evaluated the potential for low-temperature resources (<150 degrees C). Workflows, relevant datasets, a new Python library, and common and composite geological criteria maps are utilized to develop low-temperature geothermal resource favorability maps for the Denver Basin, a sedimentary basin spanning Colorado, Nebraska, and Wyoming. The replication of these methodologies in other SBGPTs can evaluate potential for low-temperature resources. To facilitate future assessment of low-temperature geothermal resources in SBGPTs, this project provides PFA workflows, data, tools, and favorability maps that will ultimately support the utilization of low-temperature geothermal resources in sedimentary basins.

combined heat and power↗
DOE OSTI2025

Advancing 1D thermo-hydraulic tools for large cryogenic facilities

The Cryogenic Division at Fermilab develops large-scale cryogenic systems for particle accelerators and superconducting magnet test facilities. To support design and diagnostics, a Python-based code was created to calculate pressure drops in components such as valves and pipes. This paper presents recent enhancements to the code, including new heat transfer functions that improve the accuracy of thermal and hydraulic predictions. The first application models pressure and temperature changes in Proton Improvement Plant PIP-II relief pipes, aiding pipe sizing and protecting relief valves. The second example analyses heat load evolution in a pipe carrying sub-atmospheric helium, helping interpret temperature sensor data and understand gas return behavior to cold compressors. These improvements significantly expand the tool's capabilities, offering a practical resource for designing and troubleshooting cryogenic systems under dynamic thermal and flow conditions.

Beckwith, R. [Fermilab]↗
DOE OSTI2021

Tiered Energy Use in Buildings (TEB)

This tool creates a set of building energy simulations and organizes them according to a tier-based model representing circuits that will be on or off depending on a microgrid controller (not included in this software). This Python based code uses the open-source RCBuildingSimulator that is under an amended MIT license. 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. SAND2021-2528 O

Villa, Daniel↗