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 595 records · Page 33

Calculation of Turbine Axial Thrust by Coupled CFD Simulations of the Main Flow Path and Secondary Cavity Flow in an SLI LOX Turbine

Each single reusable Space Launch Initiative (SLI) booster rocket is an engine operating at a record vacuum thrust level of over 730,000 Ibf using LOX and LH2. This thrust is more than 10% greater than that of the Delta IV rocket, resulting in relatively large LOX and LH2 turbopumps. Since the SLI rocket employs a staged combustion cycle the level of pressure is very high (thousands of psia). This high pressure creates many engineering challenges, including the balancing of axial-forces on the turbopumps. One of the main parameters in the calculation of the axial force is the cavity pressure upstream of the turbine disk. The flow in this cavity is very complex. The lack of understanding of this flow environment hinders the accurate prediction of axial thrust. In order to narrow down the uncertainty band around the actual turbine axial force, a coupled, unsteady computational methodology has been developed to simulate the interaction between the turbine main flow path and the cavity flow. The CORSAIR solver, an unsteady three- dimensional Navier-Stokes code for turbomachinery applications, was used to solve for both the main and the secondary flow fields. Turbine axial thrust values are presented in conjunction with the CFD simulation, together with several considerations regarding the turbine instrumentation for axial thrust estimations during test.

Dorney, D. J.↗

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↗

Noise Exposure Map of Urban Air Mobility

A noise exposure map is “a scaled geographic depiction of an airport, its noise exposure contours, noise-sensitive facilities, and land uses in the airport surrounding area” developed in accordance with the FAA’s 14 Code of Federal Regulation Part 150. This paper is the first to explore the applicability of airport noise exposure maps to Urban Air Mobility (UAM). The FAA’s airport noise compatibility planning program is first described. Then the applicability of the noise exposure map to Urban Air Mobility (UAM) is explored. Finally, new airspace infrastructure, including vertiport locations and UAM routes from NASA’s UAM engineering simulations, and local noise-sensitive facility locations and land use information were collected and processed to develop the noise exposure maps of UAM near the Dallas-Fort Worth area. The DNL noise contours resulting from a six-passenger electric quadrotor prototype vehicle are predicted using NASA’s AIRNOISEUAM software. The noise exposure maps of UAM are generated automatically using Python’s data analysis and visualization libraries. The results have applications for UAM’s noise compatibility planning, noise-reducing route planning, and vertiport location selection. (To hear the voice please utilize the video uploaded to the record)

Urban Air Mobility↗

A ROS-based Simulator for Testing the Enhanced Autonomous Navigation of the Mars 2020 Rover

In order to achieve the ambitious objectives of the Mars 2020 (M2020) mission, in particular the ability to autonomously traverse more challenging terrains more efficiently, new surface mobility software was developed for Enhanced Navigation (ENav). That decision was made early in the project, before most of the new surface flight software (FSW) existed, which created a need for a separate framework where the new navigation algorithms could be quickly prototyped and tested, before more realistic FSW-based testbeds became available. The JPL robotics team chose the Robot Operating System [1] (ROS) as the environment in which to test the new ENav algorithms. This made it possible to write the algorithms in the C language required by the FSW, so they could be directly ported over to the flight module later on, while leveraging all the C++ libraries and tools provided by ROS for simulation and testing. The ENav algorithms were developed as a separate C library, and stubs were used to replace any FSW-specific code, such as Event Reporting (EVRs) and data products (DPs). A ROS simulator was developed to generate a rich set of varied 3D terrains representative of the candidate Mars landing sites and simulate the physics of the rover motion, the point cloud perceived by the rover’s stereo vision system, and the new thinking-while-driving (TWD) navigation logic which directs the rover to drive autonomously to user-specified waypoints. To simulate the rover motion and perception, a ROS node was developed that uses a software library called HyperDrive Sim (HDSim), which is a wrapper for the Rover Sequencing and Visualization Program [2] (RSVP). That library provides roverterrain settling, realistic slip modelling, and camera rendering capability based on the rover’s NavCam machine vision models. To simulate the navigation logic, a ROS node was created that initializes and runs the ENav algorithms in a way that mimics the FSW execution, while also providing the capability to load and replay data products, including re-running the recorded inputs through the ENav algorithms for testing. An engineering Graphical User Interface (GUI) was also developed to visualize various elements, such as the rover pose during the drive, the simulated and perceived terrain, the selected local and global paths to the goal, the evaluated candidate paths and the reasons why they were rejected, the keep-in and keep-out zones (KIOZs), etc. Finally, an advanced Monte Carlo (MC) framework that can run many simulations in parallel on the Cloud and automatically generate reports that capture the key ENav performance metrics was developed to evaluate the system in a statisticallymeaningful way. This paper provides an overview of the ROSbased simulator used for testing the M2020 ENav algorithms.

Toupet, Olivier↗

The Solar Chromosphere/Corona Interface: FUV-EUV Observations and Modeling of Unresolved Coronal Funnels - I

A coronal funnel model, developed by Rabin (199l), was tested against a calibrated spectroheliogram recorded in 171 - 175 Angstrom bandpass. This image was recorded on board a sounding rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array, II (MSSTA II), MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center and the Lawrence Livermore National Laboratory, is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow bandpass, high resolution images (1 - 3 arc sec) at selected FUV, EUV and soft X-ray wavelengths (44 Angstroms - 1550 Angstroms). Using full-disk images centered at 1550 Angstroms (C IV) and 173 Angstroms (FE IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found, was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171 - 175 Angstrom bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the Chianti spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of FE IX plasma. The funnels are found to have constrictions, Gamma approx. 6 - 20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely, that a global network of funnels must cover approx. 70 - 95 % of the total solar surface area seen in the 171 - 175 Angstrom image, and = 45 % of the disk area seen in the 1550 Angstrom bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et af. (1974) and put forth in more detail by Gabriel (1976). Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by Raymond & Doyle (1981) for temperatures, 250,000 K = T = 650,000 K, based on full-disk observations made on board by SkyLab.

Martinez-Galarce, D. S.↗

An improved dataset for predicting mammal infecting viruses from genetic sequence information

There have been several attempts to develop machine learning (ML) models to identify human infecting viruses from their genomic sequences, with varying degrees of success. Direct comparison between models is problematic, because these models are typically trained and evaluated on different datasets with alternative data splitting schemes, features, and model performance metrics. In this paper we present a standardized dataset of mammal infecting and non-infecting viral pathogens, refined from the previous work of Mollentze et al. to include the latest literature evidence, roughly doubling the number of curated host-virus records available to the community, and new host target labels, primate and mammal. The new host labels were included for several reasons, including previous reports that classification performance is better at broader taxonomic ranks and the idea that there may be more data for primate infection that might serve as a suitable proxy for zoonotic potential and avoidance of false positives for human infection due to absence of evidence. On this dataset, we report the performance of eight machine learning models for predicting mammal-infecting viruses from their genomic sequences. We find that randomly assigning cases in our improved dataset to training/testing sets, when compared to the original assignments into training/testing in Mollentze et al., increases the overall average ROC AUC of prediction of human infection from 0.663 ± 0.070 to 0.784 ± 0.013, consistent with the reduction in phylogenetic distance between train and test sets (relative entropy change from 3.00 to 0.08). The broadest host category of mammal infection can be predicted most reliably at 0.850 ± 0.020. We share our improved dataset and code to enable standardized comparisons of machine learning methods to predict human host infections. Overall, we have presented preliminary evidence that classification of virus host infection is more tractable at higher taxonomic ranks, that unsurprisingly reducing the phylogenetic distance between training and test sets can improve predictive performance, that peptide kmer features appear to be harmful to out of sample model performance, and we are left with the question of whether models for virus host prediction can reasonably be expected to perform well in out of sample scenarios given the likelihood that viruses do not share a common ancestor. Consistent with this concern, when the data is resampled such that there is no overlap between viral families in training and test sets (relative entropy > 24), models perform no better than random chance at prediction of human infection regardless of whether kmers are included (ROC AUC 0.50 ± 0.08) or not (ROC AUC 0.50 ± 0.04).

59 BASIC BIOLOGICAL SCIENCES↗

Planetary Image Geometry Library

The Planetary Image Geometry (PIG) library is a multi-mission library used for projecting images (EDRs, or Experiment Data Records) and managing their geometry for in-situ missions. A collection of models describes cameras and their articulation, allowing application programs such as mosaickers, terrain generators, and pointing correction tools to be written in a multi-mission manner, without any knowledge of parameters specific to the supported missions. Camera model objects allow transformation of image coordinates to and from view vectors in XYZ space. Pointing models, specific to each mission, describe how to orient the camera models based on telemetry or other information. Surface models describe the surface in general terms. Coordinate system objects manage the various coordinate systems involved in most missions. File objects manage access to metadata (labels, including telemetry information) in the input EDRs and RDRs (Reduced Data Records). Label models manage metadata information in output files. Site objects keep track of different locations where the spacecraft might be at a given time. Radiometry models allow correction of radiometry for an image. Mission objects contain basic mission parameters. Pointing adjustment ("nav") files allow pointing to be corrected. The object-oriented structure (C++) makes it easy to subclass just the pieces of the library that are truly mission-specific. Typically, this involves just the pointing model and coordinate systems, and parts of the file model. Once the library was developed (initially for Mars Polar Lander, MPL), adding new missions ranged from two days to a few months, resulting in significant cost savings as compared to rewriting all the application programs for each mission. Currently supported missions include Mars Pathfinder (MPF), MPL, Mars Exploration Rover (MER), Phoenix, and Mars Science Lab (MSL). Applications based on this library create the majority of operational image RDRs for those missions. A Java wrapper around the library allows parts of it to be used from Java code (via a native JNI interface). Future conversions of all or part of the library to Java are contemplated.

Deen, Robert C.↗

Heading Toward Launch with the Integrated Multi-Satellite Retrievals for GPM (IMERG)

The Day-l algorithm for computing combined precipitation estimates in GPM is the Integrated Multi-satellitE Retrievals for GPM (IMERG). We plan for the period of record to encompass both the TRMM and GPM eras, and the coverage to extend to fully global as experience is gained in the difficult high-latitude environment. IMERG is being developed as a unified U.S. algorithm that takes advantage of strengths in the three groups that are contributing expertise: 1) the TRMM Multi-satellite Precipitation Analysis (TMPA), which addresses inter-satellite calibration of precipitation estimates and monthly scale combination of satellite and gauge analyses; 2) the CPC Morphing algorithm with Kalman Filtering (KF-CMORPH), which provides quality-weighted time interpolation of precipitation patterns following cloud motion; and 3) the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks using a Cloud Classification System (PERSIANN-CCS), which provides a neural-network-based scheme for generating microwave-calibrated precipitation estimates from geosynchronous infrared brightness temperatures. In this talk we summarize the major building blocks and important design issues driven by user needs and practical data issues. One concept being pioneered by the IMERG team is that the code system should produce estimates for the same time period but at different latencies to support the requirements of different groups of users. Another user requirement is that all these runs must be reprocessed as new IMERG versions are introduced. IMERG's status at meeting time will be summarized, and the processing scenario in the transition from TRMM to GPM will be laid out. Initially, IMERG will be run with TRMM-based calibration, and then a conversion to a GPM-based calibration will be employed after the GPM sensor products are validated. A complete reprocessing will be computed, which will complete the transition from TMPA.

Huffman, George J.↗

A method to generate heating and cooling schedules based on data from connected thermostats

Internet-connected thermostats are a promising new source of temperature and operational data in homes because they record inside temperatures, setpoints, and HVAC runtimes every five minutes. Over 20 million Internet-connected thermostats have already been installed in American homes. Data from about 20,000 connected thermostats were collected and organized by climate zone, number of occupants, floor area, and day type. A method was developed to create up to 40 representative temperature schedules which, together, can more accurately capture the diversity of heating and cooling behaviors. These results are suitable for input into schedules for building energy simulation models. This information enables more realistic simulations of American heating and cooling behavior, leading to more accurate estimates of energy consumption and savings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗

The Solar Chromosphere/Corona Interface. I. Far-Ultraviolet to Extreme-Ultraviolet Observations and Modeling of Unresolved Coronal Funnels

A coronal funnel model, developed by D. Rabin, was tested against a calibrated spectroheliogram recorded in the 170-1 75 A bandpass. This image was recorded on board a sounding-rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array II (MSSTA II). MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center, and the Lawrence Livermore National Laboratory' is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow-bandpass, high-resolution images (1 sec.- 3 sec.) at selected far-ultraviolet (FUV), extreme-ultraviolet (EUV), and soft X-ray wavelengths (44-1550 A). Using full disk images centered at 1550 A (C IV) and 173 A (Fe IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171-175 A bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the CHIANTI spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of Fe IX plasma. The funnels are found to have constrictions, Gamma approx. 6-20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely,that a global network of funnels must cover approx. 700 - 95% of the total solar surface area seen in the 171-175 A image, and a 5% of the disk area seen in the 1550 A bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et al. and put forth in more detail by Gabriel. Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by J. C. Raymond & J. G. Doyle for temperatures 250,000 K less than or = T less than or = 650,000 K, based on full-disk observations made on board Skylab.

Martinez-Galarce, Dennis S.↗

CalWave - TEAMER UMaine 2021 Tank Testing

This is the data submission to the teamer report for RFTS1 during which tank testing at UMaine W2 tank occurred. This submission includes the relevant information and recorded data of CalWave's UMaine wave tank testing conducted under TEAMER Request for Support Round 1. During the testing CalWave assessed an autonomous control approach for changing PTO gain parameters based on "on board" sea state detection. The data submitted includes an overview of normalized PTO gains based on sea state values and for a set of 8 relevant sea states the WEC motion. All data was collected on CalWave's SCADA system and UMaine's motion tracking system. The data submitted suffices the creation of figures as shown in the PDF report submitted alongside the data and the submission includes excel sheets (parameter gains), MATLAB structures (absorber motion), and the mentioned PDF report.

16 TIDAL AND WAVE POWER↗

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗

Full Length Assembly Testing in PELICAN (Final Report)

In support of the development of the U.S. Department of Energy (DOE) Versatile Test Reactor (VTR), a thermal hydraulics test facility was constructed to generate experimental measurement of the pressure drop across a single full-scale assembly containing prototypic axial reflectors, fuel, and plena components. Constructed and operated at Argonne National Laboratory, the Pressure drop Experimental Loop for Investigations of Core Assemblies in Nuclear reactors (PELICAN) facility was designed to achieve hydraulic conditions identical to those anticipated for a full-scale fuel assembly located in the VTR core in the region with the highest flow rate. Using water as surrogate for liquid sodium, the flow loop was operated at elevated temperatures and pressures to match the thermophysical properties of liquid sodium and ensure matching Reynolds and Euler numbers. The measurement objectives for data generated from this test facility was driven primarily by the validation needs for code calculations and simulations of the reference VTR core. These objectives focused on the need to validate pressure drop results across the various segments of the fuel assembly as they relate directly to the pumping power and safety of the reactor. Presented in this report are experimental results and analytical comparisons based on testing of a full-length assembly in PELICAN. Housed within a hexagonal test section extending 3.4 m in length, the tested assembly features a prototypic lower reflector, grid plates, wire-wrapped rod bundle, upper reflector, and exit region. The rod bundle extends over 1.5 m in length and contains 217 individual wire-wrapped rods with dimensions that best reflect the reference VTR design. The as-tested bundle assembly was fabricated using 316 stainless steel 0.25-inch (6.35-mm) diameter rods wrapped with 0.04-inch (1.016-mm) diameter wire at a helical pitch of 10.51 inch (26.6 cm). Details of the method for in-house wire-wrapping, assembly, and installation are provided later in this report. Experimental measurements of pressure drop at 19 positions along the test assembly were recorded for a range of flow conditions, with special attention paid to key locations within the assembly, including component inlet and outlet, transition, and wire-wrapped rod bundle regions. Testing conditions were based on 110°C water with flow rates ranging from 50 to 450 GPM (3 to 27 kg/s) at the inlet of the test assembly generating Reynolds numbers and velocities up to ~8.0×10 4 and ~7.8 m/s, respectively, within the rod bundle region. Non-dimensional values for the friction factor were then calculated based on these experimental measurements and compared against those predicted by various analytical correlations available from open literature. Predictions by the upgraded Cheng and Todreas, Rehme, and Novendstern correlations fell within 4% to those values measured experimentally.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fast 2D Bicephalous Convolutional Autoencoder for Compressing 3D Time Projection Chamber Data

High-energy large-scale particle colliders produce data at high speed in the order of 1 terabytes per second in nuclear physics and petabytes per second in high energy physics. Developing real-time data compression algorithms to reduce such data at high throughput to fit permanent storage has drawn increasing attention. Specifically, at the newly constructed sPHENIX experiment at the Relativistic Heavy Ion Collider (RHIC), a time projection chamber is used as the main tracking detector, which records particle trajectories in a volume of three-dimensional (3D) cylinder. The resulting data are usually very sparse with occupancy around 10.8%. Such sparsity presents a challenge to conventional learning-free lossy compression algorithms, such as SZ, ZFP, and MGARD. The 3D convolutional neural network (CNN)-based approach, Bicephalous Convolutional Autoencoder (BCAE), outperforms traditional methods both in compression rate and reconstruction accuracy. BCAE can also utilize the computation power of graphical processing units suitable for deployment in a modern heterogeneous highperformance computing environment. This work introduces two BCAE variants: BCAE++ and BCAE-2D. BCAE++ achieves a 15% better compression ratio and a 77% better reconstruction accuracy measured in mean absolute error compared with BCAE. BCAE-2D treats the radial direction as the channel dimension of an image, resulting in a 3× speedup in compression throughput. In addition, we demonstrate an unbalanced autoencoder with a larger decoder can improve reconstruction accuracy without significantly sacrificing throughput. Lastly, we observe both the BCAE++ and BCAE-2D can benefit more from using half-precision mode in throughput (76 - 79% increase) without loss in reconstruction accuracy. The source code and links to data and pretrained models can be found at https://github.com/BNL-DAQ-LDRD/NeuralCompression_v2

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Suomi National Polar-Orbiting Partnership (SNPP): Continuing NASA Research and Applications

The Suomi National Polar-orbiting Partnership (SNPP) satellite was successfully launched into a polar orbit on October 28, 2011 carrying 5 remote sensing instruments designed to provide data to improve weather forecasts and to increase understanding of long-term climate change. SNPP provides operational continuity of satellite-based observations for NOAA's Polar-orbiting Operational Environmental Satellites (POES) and continues the long-term record of climate quality observations established by NASA's Earth Observing System (EOS) satellites. In the 2003 to 2011 pre-launch timeframe, NASA's SNPP Science Team assessed the adequacy of the operational Raw Data Records (RDRs), Sensor Data Records (SDRs), and Environmental Data Records (EDRs) from the SNPP instruments for use in NASA Earth Science research, examined the operational algorithms used to produce those data records, and proposed a path forward for the production of climate quality products from SNPP. In order to perform these tasks, a distributed data system, the NASA Science Data Segment (SDS), ingested RDRs, SDRs, and EDRs from the NOAA Archive and Distribution and Interface Data Processing Segments, ADS and IDPS, respectively. The SDS also obtained operational algorithms for evaluation purposes from the NOAA Government Resource for Algorithm Verification, Independent Testing and Evaluation (GRAVITE). Within the NASA SDS, five Product Evaluation and Test Elements (PEATEs) received, ingested, and stored data and performed NASA's data processing, evaluation, and analysis activities. The distributed nature of this data distribution system was established by physically housing each PEATE within one of five Climate Analysis Research Systems (CARS) located at either at a NASA or a university institution. The CARS were organized around 5 key EDRs directly in support of the following NASA Earth Science focus areas: atmospheric sounding, ocean, land, ozone, and atmospheric composition products. The PEATES provided the system level interface with members of the NASA SNPP Science Team and other science investigators within each CARS. A sixth Earth Radiation Budget CARS was established at NASA Langley Research Center (NASA LaRC) to support instrument performance, data evaluation, and analysis for the SNPP Clouds and the Earth's Radiant Budget Energy System (CERES) instrument. Following the 2011 launch of SNPP, spacecraft commissioning, and instrument activation, the NASA SNPP Science Team evaluated the operational RDRs, SDRs, and EDRs produced by the NOAA ADS and IDPS. A key part in that evaluation was the NASA Science Team's independent processing of operational RDRs and SDRs to EDRs using the latest NASA science algorithms. The NASA science evaluation was completed in the December 2012 to April 2014 timeframe with the release of a series of NASA Science Team Discipline Reports. In summary, these reports indicated that the RDRs produced by the SNPP instruments were of sufficiently high quality to be used to create data products suitable for NASA Earth System science and applications. However, the quality of the SDRs and EDRs were found to vary greatly when considering suitability for NASA science. The need for improvements in operational algorithms, adoption of different algorithmic approaches, greater monitoring of on-orbit instrument calibration, greater attention to data product validation, and data reprocessing were prominent findings in the reports. In response to these findings, NASA, in late 2013, directed the NASA SNPP Science Team to use SNPP instrument data to develop data products of sufficiently high quality to enable the continuation of EOS time series data records and to develop innovative, practical applications of SNPP data. This direction necessitated a transition of the SDS data system from its pre-launch assessment mode to one of full data processing and production. To do this, the PEATES, which served as NASA's data product testing environment during the prelaunch and early on-orbit periods, were transitioned to Science Investigator-led Processing Systems (SIPS). The distributed data architecture was maintained in this new system by locating the SIPS at the same institutions at which the CARS and PEATES were located. The SIPS acquire raw SNPP instrument Level 0 (i.e. RDR) data over the full SNPP mission from the NOAA ADS and IDPS through the NASA SDS Data Distribution and Depository Element (SD3E). The SIPS process those data into NASA Level 1, Level 2, and global, gridded Level 3 standard products using peer-reviewed algorithms provided by members of the NASA Science Team. The SIPS work with the NASA SNPP Science Team in obtaining enhanced, refined, or alternate real-time algorithms to support the capabilities of the Direct Readout Laboratory (DRL). All data products, algorithm source codes, coefficients, and auxiliary data used in product generation are archived in an assigned NASA Distributed Active Archive Center (DAAC).

remote sensing↗

GBaTSv2: a revised synthesis of the likely basal thermal state of the Greenland Ice Sheet

The basal thermal state (frozen or thawed) of the Greenland Ice Sheet is under-constrained due to few direct measurements, yet knowledge of this state is becoming increasingly important to interpret modern changes in ice flow. The first synthesis of this state relied on inferences from widespread airborne and satellite observations and numerical models, for which most of the underlying datasets have since been updated. Further, new and independent constraints on the basal thermal state have been developed from analysis of basal and englacial reflections observed by airborne radar sounding. Here we synthesize constraints on the Greenland Ice Sheet's basal thermal state from boreholes, thermomechanical ice-flow models that participated in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6; Coupled Model Intercomparison Project Phase 6), IceBridge BedMachine Greenland v4 bed topography, Making Earth Science Data Records for Use in Research Environments (MEaSUREs) Multi-Year Greenland Ice Sheet Velocity Mosaic v1 and multiple inferences of a thawed bed from airborne radar sounding. Most constraints can only identify where the bed is likely thawed rather than where it is frozen. This revised synthesis of the Greenland likely Basal Thermal State version 2 (GBaTSv2) indicates that 33 % of the ice sheet's bed is likely thawed, 40 % is likely frozen and the remainder (28 %) is too uncertain to specify. The spatial pattern of GBaTSv2 is broadly similar to the previous synthesis, including a scalloped frozen core and thawed outlet-glacier systems. Although the likely basal thermal state of nearly half (46 %) of the ice sheet changed designation, the assigned state changed from likely frozen to likely thawed (or vice versa) for less than 6 % of the ice sheet. This revised synthesis suggests that more of northern Greenland is likely thawed at its bed and conversely that more of southern Greenland is likely frozen, both of which influence interpretation of the ice sheet's present subglacial hydrology and models of its future evolution. The GBaTSv2 dataset, including both code that performed the analysis and the resulting datasets, is freely available at https://doi.org/10.5281/zenodo.6759384 (MacGregor, 2022).

Joseph A. MacGregor↗