Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Record”

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 415 records · Page 23

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↗

Multispectral and thermal surface imagery and surface elevation mosaics (camspec-air)

This dataset contains high resolution image products (orthomosaics) acquired from midsized uncrewed aerial systems, which have been processed for value added quality. The instrument itself, a multispectral imager, the Altum by Micasense, captures 6 spectral bands (red, green blue, NIR, red edge, and LWIR/thermal1) as radiance, which is converted to reflectance. The code used to develop these images first uses tools from the Micasense python library2 to apply dark level corrections, row gradient corrections, and radiometric corrections. Next it uses the processing API from Agisoft Metashape software to align and mosaic the processed imagery, following the processes developed by the USGS' structure from motion workflow documentation3. Captures at different altitudes (recorded in MSL) produce an orthomosaic, a tif image containing information related to the 6 spectral bands, and a digital elevation model (DEM), a tif image containing information related to the elevation of the surveyed terraine. Metadata included in every image can be used to extract lat, lon, and reflectance values. 1https://www.arm.gov/publications/tech_reports/handbooks/doe-sc-arm-tr-281.pdf 2https://micasense.github.io/imageprocessing/MicaSense%20Image%20Processing%20Setup.html 3https://pubs.usgs.gov/of/2021/1039/ofr20211039.pdf

54 ENVIRONMENTAL SCIENCES↗

Multispectral and thermal surface imagery and surface elevation mosaics - SGP July 2022

This data set contains high-resolution image products (orthomosaics) acquired from midsized uncrewed aerial systems that have been processed for value-added quality. The instrument itself, a multispectral imager, the Altum by Micasense, captures six spectral bands (red, green blue, NIR, red edge, and LWIR/thermal1) as radiance, which is converted to reflectance via custom code. The code used to develop these images first uses tools from the Micasense Python library2 to apply dark level corrections, row gradient corrections, and radiometric corrections. Next, it uses the processing API from Agisoft Metashape software to align and mosaic the processed imagery, following the processes developed by the USGS' structure from motion workflow documentation.3 Captures at different altitudes (recorded in MSL) produce an orthomosaic, a tif image containing information related to the six spectral bands, and a digital elevation model (DEM), a tif image containing information related to the elevation of the surveyed terraine. Metadata included in every image can be used to extract lat, lon, and reflectance values. 1 https://www.arm.gov/publications/tech_reports/handbooks/doe-sc-arm-tr-281.pdf 2 https://micasense.github.io/imageprocessing/MicaSense%20Image%20Processing%20Setup.html 3 https://pubs.usgs.gov/of/2021/1039/ofr20211039.pdf

54 ENVIRONMENTAL SCIENCES↗

Analysis of Buzz in a Supersonic Inlet

A dual-stream, low-boom supersonic inlet designed for use on a small, Mach 1.6 aircraft was tested experimentally in the 8- by 6-Foot Supersonic Wind Tunnel (SWT) at the NASA Glenn Research Center (GRC). The tests showed that the inlet had good recovery and stable operation over large mass flow range. The inlet went into buzz at mass flows well below that needed for engine operation, and the experiments generated a wealth of data during buzz. High frequency response pressure measurements and high-speed schlieren videos were recorded for many buzz events. The objective of the present work was to use computational fluid dynamics (CFD) to predict some of the experimental data taken during buzz, compare those predictions to the experimental data, and to use both datasets to explain the physics of the buzz cycle. The calculations were done with the Wind-US CFD code using a second-order time-accurate differencing scheme and the SST turbulence model. Computed Mach number contours were compared with schlieren images, and ensemble-averaged unsteady pressures were compared to data. The results showed that the buzz cycle consisted partly of spike buzz, an unsteady oscillation of the main shock at the spike tip while the inlet pressure dropped, and partly of choked flow while the inlet repressurized. Most of the results could be explained by theory proposed by Dailey in 1954, but did not support commonly used acoustic resonance explanations.

Chima, Rodrick V.↗

Data and code from: Multivariate bayesian regression model for predicting disposed ash composition at U.S. coal fired power stations

This dataset contains the code and data files needed for implementation of a Multivariate Bayesian Regression model, described in Jin et al. (2025), for the historical prediction of the chemical composition of disposed coal ash at U.S. coal fired power plants as a function of annualized coal purchase data. The integrated coal supply data file (CoalSupplyDataset.csv) represents a compilation of monthly fuel purchase records for the period 1973-2022 at major U.S. power stations. These records were obtained from the U.S. Energy Information Administration. The CSV file also contains, for each coal purchase record, the coal region of the mine as defined by the U.S. Geological Survey. Data entry errors and data gaps in the EIA records were corrected as described in Jin et al. This CSV file represents the integrated coal supply data after corrections were made. The model structure and fitting parameters are encoded in pickle file format (Bayesian.pkl). The model was developed with the coal supply data and coal ash composition data, apportioned according to the Stratified Shuffle Split for training and testing subsets. The model was built using Python and the PyMC library. Reference Publication: Jin, Z.; Huang, J.; Hower, J.C.; Hsu-Kim, H.(2025). Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites. Environmental Science & Technology.

Coal ash composition↗

GOES-R Geostationary Lightning Mapper Performance Specifications and Algorithms

The Geostationary Lightning Mapper (GLM) is a single channel, near-IR imager/optical transient event detector, used to detect, locate and measure total lightning activity over the full-disk. The next generation NOAA Geostationary Operational Environmental Satellite (GOES-R) series will carry a GLM that will provide continuous day and night observations of lightning. The mission objectives for the GLM are to: (1) Provide continuous, full-disk lightning measurements for storm warning and nowcasting, (2) Provide early warning of tornadic activity, and (2) Accumulate a long-term database to track decadal changes of lightning. The GLM owes its heritage to the NASA Lightning Imaging Sensor (1997- present) and the Optical Transient Detector (1995-2000), which were developed for the Earth Observing System and have produced a combined 13 year data record of global lightning activity. GOES-R Risk Reduction Team and Algorithm Working Group Lightning Applications Team have begun to develop the Level 2 algorithms and applications. The science data will consist of lightning "events", "groups", and "flashes". The algorithm is being designed to be an efficient user of the computational resources. This may include parallelization of the code and the concept of sub-dividing the GLM FOV into regions to be processed in parallel. Proxy total lightning data from the NASA Lightning Imaging Sensor on the Tropical Rainfall Measuring Mission (TRMM) satellite and regional test beds (e.g., Lightning Mapping Arrays in North Alabama, Oklahoma, Central Florida, and the Washington DC Metropolitan area) are being used to develop the prelaunch algorithms and applications, and also improve our knowledge of thunderstorm initiation and evolution.

Mach, Douglas M.↗

Experiment and simulation of high-speed gas jet penetration into a semicircular fluidized bed

This work marks the third in a series of experiments that were in a semi-circular, gas-fluidized bed with side jets. In this work, the particles are (nominally) 1 mm ceramic beads. The bed is operated just at and slightly above and below the minimum fluidization velocity and additional fluidization is provided by two high-speed gas located on the sides of the bed near the flat, front face of the unit. Two primary measurements are taken: high-speed video recording of the front of the bed and bed pressure drop from a tap in the back of the bed. Particle Image Velocimetry (PIV) is used to determine particle motion, characterized as a mean Froude number, from the high-speed video. A CFD-DEM model of the bed is presented using the recently released MFIX-Exa code. Four model subvariants are considered using two methods of representing the jets and two drag models, both of which are calibrated to exactly match the experimentally measured minimum fluidization velocity. Although it is more difficult to determine the jet penetration depths in a straightforward manner as in the previous works using Froude number contours, the CFD-DEM results compare quite well to the PIV measurements, particularly for submodel flow Syam. Unfortunately, the good agreement of the solids-phase is overshadowed by significant disagreement in the gas-phase data. Specifically, the predicted time averaged standard deviation of the pressure drop is found to be over an order of magnitude larger than measured. Due to the low value of the measurements, just 1% of the mean bed pressure drop, it seems possible that the data is in error. On the other hand, the model may not be accurately capturing pressure attenuation through an under-fluidized region in the back of the bed. Without the possibility additional experiments to test the validity of the data, this work is simply being reported “as is” without being able to indicate which, either the simulation or the experiment, is more correct.<br>

Fullmer, William D.↗

NEMA-Phase Compliant Traffic Signal Controller Module in SUMO

The controller modules in SUMO use a stage-based control structure. A phase is defined as a stage of all allowed movements at a time instance. However, traffic signal controllers used in North America widely use National Electrical Manufacturers Association (NEMA) phase definition. A NEMA phase is defined by a certain flow movement at an intersection. At one time, more than one NEMA phase could happen together as long as they do not conflict with each other. We can visualize the NEMA phases and timings in Ring-and-Barrier structured NEMA diagrams. For one controller, only one phase from a ring can be activated at a time. Phases from different rings could be activated together as long as they are not from the different sides of a barrier. When a controller is operated in fixed-time control mode, we can model the NEMA phase timing as a corresponding stage-based control timing without any issues. When introducing actuation into the signal control, a Ring-and-Barrier structured traffic signal controller can be more flexible than stage-based controller by allowing different possible phase combinations. We made two efforts in modeling Ring-and-Barrier structured controllers in SUMO. One is to translate a NEMA phases timing into SUMO-readable phases and timings as an additional file for SUMO. This translation worked well for fixed-time control. To model actuated control and coordinated actuated control, we augmented the SUMO source code by adding a Ring-and-Barrier structured controller module. This module could implement traffic signal timing from controllers using NEMA phases. We also augmented TraCI to be able to set new NEMA phase timings during simulations. We examined the Ring-and-Barrier structured traffic signal controller module by both visually observing the simulation animations and the simulation records. The developed control module can model the generalized Ring-and-Barrier structured traffic signal timing that is used in North America. SEE: https://github.com/eclipse/sumo/blob/main/src/microsim/traffic_lights/NEMAController.cpp

Wang, Qichao↗

TOGA/COARE AMMR 1992 data processing

The complete set of Tropical Ocean and Global Atmosphere (TOGA)/Coupled Ocean Atmosphere Response Experiment (COARE) flight data for the 91.65 GHz Airborne Meteorological Radiometer (AMMR92) contains data from nineteen flights: two test flights, four transit flights, and thirteen experimental flights. The data flight occurred between December 16, 1992 and February 28, 1993. Data collection from the AMMR92 during the first ten flights of TOGA/COARE was performed using the executable code TSK30041. These are IBM PC/XT programs used by the NASA Goddard Space Flight Center (GSFC). During one flight, inconsistencies were found during the operation of the AMMR92 using the GSFC data acquisition system. Consequently, the Georgia Tech (GT) data acquisition system was used during all successive TOGA/COARE flights. These inconsistencies were found during the data processing to affect the recorded data as well. Errors are caused by an insufficient pre- and post-calibration setting period for the splash-plate mechanism. The splash-plate operates asynchronusly with the data acquisition system (there is no position feedback to the GSFC or GT data system). This condition caused both the calibration and the post-calibration scene measurement to be corrupted on a randomly occurring basis when the GSFC system was used. This problem did not occur with the GT data acquisition system due to sufficient allowance for splash-plate settling. After TOGA/COARE it was determined that calibration of the instrument was a function of the scene brightness temperature. Therefore, the orientation error in the main antenna beam of the AMMR92 is hypothesized to be caused by misalignment of the internal 'splash-plate' responsible for directing the antenna beam toward the scene or toward the calibration loads. Misalignment of the splash-plate is responsible for 'scene feedthrough' during calibration. Laboratory investigation at Georgia Tech found that each polarization is affected differently by the splash-plate alignment error. This is likely to cause significant and unique errors in the absolute calibration of each channel.

Kunkee, D. B.↗

Dynamic Deformation Measurements of an Aeroelastic Semispan Model

The techniques used to acquire, reduce, and analyze dynamic deformation measurements of an aeroelastic semispan wind tunnel model are presented. Single-camera, single-view video photogrammetry (also referred to as videogrammetric model deformation, or VMD) was used to determine dynamic aeroelastic deformation of the semispan 'Models for Aeroelastic Validation Research Involving Computation' (MAVRIC) model in the Transonic Dynamics Tunnel at the NASA Langley Research Center. Dynamic deformation was determined from optical retroreflective tape targets at five semispan locations located on the wing from the root to the tip. Digitized video images from a charge coupled device (CCD) camera were recorded and processed to automatically determine target image plane locations that were then corrected for sensor, lens, and frame grabber spatial errors. Videogrammetric dynamic data were acquired at a 60-Hz rate for time records of up to 6 seconds during portions of this flutter/Limit Cycle Oscillation (LCO) test at Mach numbers from 0.3 to 0.96. Spectral analysis of the deformation data is used to identify dominant frequencies in the wing motion. The dynamic data will be used to separate aerodynamic and structural effects and to provide time history deflection data for Computational Aeroelasticity code evaluation and validation.

Graves, Sharon S.↗

Calculation of Dose Deposition in 3D Voxels by Heavy Ions and Simulation of gamma-H2AX Experiments

The biological response to high-LET radiation is different from low-LET radiation due to several factors, notably difference in energy deposition and formation of radiolytic species. Of particular importance in radiobiology is the formation of double-strand breaks (DSB), which can be detected by -H2AX foci experiments. These experiments has revealed important differences in the spatial distribution of DSB induced by low- and high-LET radiations [1,2]. To simulate -H2AX experiments, models based on amorphous track with radial dose are often combined with random walk chromosome models [3,4]. In this work, a new approach using the Monte-Carlo track structure code RITRACKS [5] and chromosome models have been used to simulate DSB formation. At first, RITRACKS have been used to simulate the irradiation of a cubic volume of 5 m by 1) 450 1H+ ions of 300 MeV (LET 0.3 keV/ m) and 2) by 1 56Fe26+ ion of 1 GeV/amu (LET 150 keV/ m). All energy deposition events are recorded to calculate dose in voxels of 20 m. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. Many differences are found in the spatial distribution of dose voxels for the 56Fe26+ ion. The track structure can be distinguished, and voxels with very high dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and indicate clustered energy deposition, which may be responsible for complex DSB. In the second step, assuming that DSB will be found only in voxels where energy is deposited by the radiation, the intersection points between voxels with dose > 0 and simulated chromosomes were obtained. The spatial distribution of the intersection points is similar to -H2AX foci experiments. These preliminary results suggest that combining stochastic track structure and chromosome models could be a good approach to understand radiation-induced DSB and chromosome aberrations.

Plante, I.↗

Retrieving Atmospheric Temperature and Moisture Profiles from NPP CRIS/ATMS Sensors Using Crimss EDR Algorithm

The Joint Polar Satellite System (JPSS) is a U.S. National Oceanic and Atmospheric Administration (NOAA) mission in collaboration with the U.S. National Aeronautical Space Administration (NASA) and international partners. The NPP Cross-track Infrared Microwave Sounding Suite (CrIMSS) consists of the infrared (IR) Crosstrack Infrared Sounder (CrIS) and the microwave (MW) Advanced Technology Microwave Sounder (ATMS). The CrIS instrument is hyperspectral interferometer, which measures high spectral and spatial resolution upwelling infrared radiances. The ATMS is a 22-channel radiometer similar to Advanced Microwave Sounding Units (AMSU) A and B. It measures top of atmosphere MW upwelling radiation and provides capability of sounding below clouds. The CrIMSS Environmental Data Record (EDR) algorithm provides three EDRs, namely the atmospheric vertical temperature, moisture and pressure profiles (AVTP, AVMP and AVPP, respectively), with the lower tropospheric AVTP and the AVMP being JPSS Key Performance Parameters (KPPs). The operational CrIMSS EDR an algorithm was originally designed to run on large IBM computers with dedicated data management subsystem (DMS). We have ported the operational code to simple Linux systems by replacing DMS with appropriate interfaces. We also changed the interface of the operational code so that we can read data from both the CrIMSS science code and the operational code and be able to compare lookup tables, parameter files, and output results. The detail of the CrIMSS EDR algorithm is described in reference [1]. We will present results of testing the CrIMSS EDR operational algorithm using proxy data generated from the Infrared Atmospheric Sounding Interferometer (IASI) satellite data and from the NPP CrIS/ATMS data.

Liu, X.↗

Climate model postprocessing

The development of new postprocessing software of the climate modeling group is summarized. Code, test, and perform simulations with global general circulation models are described. The models improve understanding and ability to predict the vagaries of weather and climate. To learn from and utilize the model results, it is necessary to create elaborate postprocessing software to allow analysis of the large volume of data produced. The models produce sigma history tapes. The sigma history records are interpolated to pressure history records, which are written on a pressure history tape. The model results are analyzed on pressure surfaces, with snap shots or time averages.

Abeles, J.↗

Prospects for gamma-ray line observations of individual supernovae

The gamma-ray line emission from individual type I and type II supernovae are studied using numerical simulations and photon propagation codes to predict the flux levels and line shapes. For both types, the gamma ray lines with the highest flux from an individual event are the 0.847 and 1.238 MeV lines from the Ni-56 to Co-56 to Fe-56 decay chain. For type I supernovae, the 0.847 MeV line peaks at about 70 days after event onset. The historical record indicates an approximate discovery rate of once in 10 years for balloon-borne instruments, once in two to three years for the Gamma-Ray Observatory, and once in one to two years for a proposed space mission. The 0.847 MeV line flux from type II supernovae peaks at about 600 days after event onset at a low level which restricts observations to the events in the Galaxy and its nearest neighbors. The expected line shape is narrower than for type I supernovae.

Gehrels, Neil↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

Lossless compression of AVIRIS data: Comparison of methods and instrument constraints

A family of lossless compression methods, allowing exact image reconstruction, are evaluated for compressing Airborne Visible/Infrared Imaging Spectrometers (AVIRIS) image data. The methods are used on Differential Pulse Code Modulation (DPCM). The compressed data have an entropy of order 6 bits/pixel. A theoretical model indicates that significantly better lossless compression is unlikely to be achieved because of limits caused by the noise in the AVIRIS channels. AVIRIS data differ from data produced by other visible/near-infrared sensors, such as LANDSAT-TM or SPOT, in several ways. Firstly, the data are recorded at a greater resolution (12 bits, though packed into 16-bit words). Secondly, the spectral channels are relatively narrow and provide continuous coverage of the spectrum so that the data in adjacent channels are generally highly correlated. Thirdly, the noise characteristics of the AVIRIS are defined by the channels' Noise Equivalent Radiances (NER's), and these NER's show that, at some wavelengths, the least significant 5 or 6 bits of data are essentially noise.

Roger, R. E.↗

Magnetically Actuated Propellant Orientation Experiment, Controlling Fluid Motion With Magnetic Fields in a Low-Gravity Environment

This report details the results of a series of fluid motion experiments to investigate the use of magnets to orient fluids in a low-gravity environment. The fluid of interest for this project was liquid oxygen (LO2) since it exhibits a paramagnetic behavior (is attracted to magnetic fields). However, due to safety and handling concerns, a water-based ferromagnetic mixture (produced by Ferrofluidics Corporation) was selected to simplify procedures. Three ferromagnetic fluid mixture strengths and a nonmagnetic water baseline were tested using three different initial fluid positions with respect to the magnet. Experiment accelerometer data were used with a modified computational fluid dynamics code termed CFX-4 (by AEA Technologies) to predict fluid motion. These predictions compared favorably with experiment video data, verifying the code's ability to predict fluid motion with and without magnetic influences. Additional predictions were generated for LO2 with the same test conditions and geometries used in the testing. Test hardware consisted of a cylindrical Plexiglas tank (6-in. bore with 10-in. length), a 6,000-G rare Earth magnet (10-in. ring), three-axis accelerometer package, and a video recorder system. All tests were conducted aboard the NASA Reduced-Gravity Workshop, a KC-135A aircraft.

Martin, J. J.↗

Measurement of local convective heat transfer coefficients from a smooth and roughened NACA-0012 airfoil: Flight test data

Wind tunnels typically have higher free stream turbulence levels than are found in flight. Turbulence intensity was measured to be 0.5 percent in the NASA Lewis Icing Research Tunnel (IRT) with the cloud making sprays off and around 2 percent with cloud making equipment on. Turbulence intensity for flight conditions was found to be too low to make meaningful measurements for smooth air. This difference between free stream and wing tunnel conditions has raised questions as to the validity of results obtained in the IRT. One objective of these tests was to determine the effect of free stream turbulence on convective heat transfer for the NASA Lewis LEWICE ice growth prediction code. These tests provide in-flight heat transfer data for a NASA-0012 airfoil with a 533 cm chord. Future tests will measure heat transfer data from the same airfoil in the Lewis Icing Research Tunnel. Roughness was obtained by the attachment of small, 2 mm diameter hemispheres of uniform size to the airfoil in three different patterns. Heat transfer measurements were recorded in flight on the NASA Lewis Twin Otter Icing Research Aircraft. Measurements were taken for the smooth and roughened surfaces at various aircraft speeds and angles of attack up to four degrees. Results are presented as Frossling number versus position on the airfoil for various roughnesses and angles of attack.

Newton, James E.↗