SEARCH · Engineering Papers

Results for “forward sensitivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Neutrino electromagnetic properties and the weak mixing angle at the LHC Forward Physics Facility

The LHC produces an intense beam of highly energetic neutrinos of all three flavors in the forward direction, and the Forward Physics Facility (FPF) has been proposed to house a suite of experiments taking advantage of this opportunity. In this study, we investigate the FPF’s potential to probe the neutrino electromagnetic properties, including neutrino millicharge, magnetic moment, and charge radius. We find that, due to the large flux of tau neutrinos at the LHC, the FPF detectors will be able to provide more sensitive constraints on the tau neutrino magnetic moment and millicharge than previous measurements at DONUT, by searching for excess in low recoil energy electron scattering events. We also find that, by precisely measuring the rate of neutral current deep inelastic scattering events, the FPF detectors have the potential to obtain the strongest experimental bounds on the neutrino charge radius for the electron neutrino, and one of the leading bounds for the muon neutrino flavor. The same signature could also be used to measure the weak mixing angle, and we estimate that sin 2 θ W could be measured to about 3% precision at a scale Q ∼ 10 GeV , shedding new light on the longstanding NuTeV anomaly. Published by the American Physical Society 2025

Abraham, Roshan Mammen (ORCID:0000000346783808)↗

DOE OSTI2024

Revenue Analysis for Energy Storage Systems in the United States

In this work we evaluate the potential revenue from energy storage using historical electricity prices, forward-looking projections of hourly electricity prices, and actual reported revenue. This analysis examines the impact of storage characteristics, specifically duration and round-trip efficiency, as well as locational elements of storage revenue within the current and projected U.S. power system. Figure ES-1 illustrates the revenue for a 1 MW storage system in seven market regions with durations range from 1 hour to 12 hours using both historical and forward-looking price data. The historical analysis covers more than 500 price nodes for each market region, while the forward-looking analysis includes balancing areas under different 10 scenarios of the electricity generation mix. The results indicate that the revenues consistently increase with duration, though the marginal value declines as duration grows. Moreover, the range of revenue depends on the system's operational location, and the electricity generation mix changes for future years. This range also widens with increased durations. In addition, the sensitivity analysis of round-trip efficiency reveals that as efficiency improves, system revenue increases, though the value of better round-trip efficiency declines as at higher efficiency levels.

25 ENERGY STORAGE↗

NASA NTRS2011

The RHESSI Microflare Height Distribution

We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6-10 keV count-rate when RHESSI's full sensitivity was available in order to find the smallest events (Christe et al., 2008). Between March 2002 and March 2007, a total of 25,006 events were found. Source locations were determined in the 4-10 keV, 10-15 keV, and 15-30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g (h) oc exp(-h/lambda) where lambda = 6.1 plus or minus 0.3 Mm is the scale height. A power law height distribution with a negative power law index, gamma = 3.1 plus or minus 0.3 is also consistent with the data. Interpreted as thermal loop top sources, these heights are compared to loops generated by a potential field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential field loop height distribution which may be a signature of the flare energization process.

Christe, P.↗

DOE OSTI2023

Noble gas insights into early impact delivery and volcanic outgassing to Earth's atmosphere: A limited role for the continental crust

Earth's atmosphere, crust and mantle have evolved together through continuous geochemical exchange throughout Earth's history. Constraints on the transport of volatile elements and compounds between these reservoirs are crucial for understanding how Earth could have stayed habitable for extended periods of time. In this work, we present a new forward model of He, Ne and Ar in the mantle, crust and atmosphere. We explore concentrations of noble gases at the end of accretion, bulk silicate Earth K/U ratios, crustal growth scenarios, and upper mantle processing rates throughout Earth's history. We search for parameter combinations that simultaneously satisfy observational constraints on present-day mantle 4 He/ 3 He and 40 Ar/ 36 Ar (sensitive to mantle outgassing and continental crust growth, which depletes the mantle of U, Th and K), atmospheric 20 Ne/ 22 Ne (which tracks the mix of outgassed vs. delivered Ne in the atmosphere), and atmospheric 40 Ar/ 36 Ar in the past and today (sensitive to volatile delivery, mantle outgassing and crustal growth). Leveraging this intertwined set of noble gas abundances and isotopic compositions yields new constraints on initial noble gas abundances in the mantle and on the proportions of atmospheric volatiles originating from delivery by impact degassing, mantle outgassing, and degassing of the continental crust. We find that atmospheric Ar isotopic evolution is primarily sensitive to the mantle processing rate history; the atmospheric Ar isotopic record should therefore not be used to reconstruct continental crust growth, but instead provides valuable insights into mantle processing rates. Our model predicts a measurably low 20 Ne/ 22 Ne ratio of ~9.7 in Archean atmospheric samples. Most of atmospheric primordial 36 Ar was directly delivered by chondritic bodies and not transferred to the atmosphere during an intense early episode of mantle outgassing. Nitrogen delivered by impact degassing could account for the present-day atmospheric nitrogen inventory.

58 GEOSCIENCES↗

DOE OSTI2023

Precision engineering of biological function with large-scale measurements and machine learning

As synthetic biology expands and accelerates into real-world applications, methods for quantitatively and precisely engineering biological function become increasingly relevant. This is particularly true for applications that require programmed sensing to dynamically regulate gene expression in response to stimuli. However, few methods have been described that can engineer biological sensing with any level of quantitative precision. Here, we present two complementary methods for precision engineering of genetic sensors: in silico selection and machine-learning-enabled forward engineering. Both methods use a large-scale genotype-phenotype dataset to identify DNA sequences that encode sensors with quantitatively specified dose response. First, we show that in silico selection can be used to engineer sensors with a wide range of dose-response curves. To demonstrate in silico selection for precise, multi-objective engineering, we simultaneously tune a genetic sensor’s sensitivity (EC 50 ) and saturating output to meet quantitative specifications. In addition, we engineer sensors with inverted dose-response and specified EC 50 . Second, we demonstrate a machine-learning-enabled approach to predictively engineer genetic sensors with mutation combinations that are not present in the large-scale dataset. We show that the interpretable machine learning results can be combined with a biophysical model to engineer sensors with improved inverted dose-response curves.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2018

Wide-field Infrared Survey Telescope (WFIRST): Composite Structure Verification for Operational Temperatures

The Wide-Field Infrared Survey Telescope (WFIRST) mission is the next large astrophysics observatory for NASA after the James Webb Space Telescope and is the top priority mission from the 2010 National Academy of Sciences' decadal survey. The WFIRST Optical Telescope Assembly (OTA) includes inherited composite support structures that were originally designed and tested for room temperature operation; however, the WFIRST mission will require operation at colder temperatures to achieve sufficient sensitivity for the infrared wavelengths. We will present the results and conclusions of testing completed at the coupon and engineering model level to verify that the inherited composite structures will maintain mechanical integrity and performance over the required temperature range. The testing included: (1) characterization testing of constituent material coupons, (2) thermal cycling and static load testing of a representative aft metering structure (AMS) and forward metering structure (FMS), and (3) thermal cycling and dynamic testing of a representative secondary mirror assembly (SMA).

metering↗

NASA NTRSDate not supplied

Scaling of Liquefaction Systems and the Effects of Gravity

Recent exploration initiatives both by NASA and others have produced interest in the liquefaction of fluids produced through in-situ processes on the surfaces of the Moon and Mars. Liquefaction of fluids in cryogenic temperatures is routinely done at large scale for terrestrial uses, but these processes do not necessarily scale down in an economical or technologically feasible manner for the anticipated rates of initial Lunar or Martian production plants. To understand appropriate processes and scaling parameters, various options should be considered. A multi-center team at NASA considered multiple different refrigeration cycles and refrigeration integration methodologies as well as how these might fit into early liquefaction plants. These studies resulted in the conclusion that integrating the tubing on the storage tank wall (for both structural and heat spreading reasons), preferably the tank which will actually use the fluid, is the nearest path forward. Given that assumption on the general framework of the liquefaction system, there are several different heat transfer mechanisms to consider from a scaling perspective. These include forced convection heat removal to the refrigeration system (or cryocooler), conduction through the tank wall heat exchanger, and condensation (with some natural convection) on the inner tank wall. Analysis of these mechanisms shows that while there is some sensitivity to gravitational level, that within the bounds of current interest (Lunar and Martian applications), this sensitivity does not dominate the liquefaction application. While there will be some effect, perhaps in a limiting manner as the tanks approach some high fill level, system level testing on the Earth should suffice for the performance prediction and demonstration of liquefaction operations as applicable to Lunar and Martian applications. However, as one approaches orbital conditions, this general approach to liquefaction will not be appropriate as the scaling of the heat transfer mechanisms is not appropriate. At this point, other approaches will need to be developed and demonstrated in the micro-gravity environment.

Cryogenics↗

DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

DOE OSTI2020

Isopropanol dehydration reaction rate kinetics measurement using H 2 O time histories

H 2 O formation during thermal decomposition of isopropanol was measured behind reflected shock waves at temperatures ranging from 1127 to 1621 K at an average pressure of 1.42 atm using a laser absorption technique. Of the five modern chemical kinetics models used to compare the H 2 O time histories, the model from Li et al (Combust. Flame 2019;207:171-185) showed the best overall agreement. Sensitivity and rate of production analyses using the Li et al model (as well as those from AramcoMech 3.0, CRECK, and Togbe et al (Energy Fuel. 2011;25:676-683)) showed unimolecular dehydration of isopropanol, iC 3 H 7 OH ⇌ H 2 O + C 3 H 6 (R1), is nearly the sole reaction controlling H2O production at early times, allowing for an a priori measurement of the forward rate constant k 1 . The Arrhenius expression k1 (s –1 ) = 2.60 × 10 13 exp(-31 120 K/T) was determined to represent best the data from this study. According to the models and previous experimental investigations, the pressure investigated is well within the high-pressure limit (HPL) for this reaction. Additional, higher-pressure experiments also confirmed the HPL assumption. An uncertainty analysis was performed by varying secondary reactions within their uncertainties and examining their effect on the overall prediction, establishing an uncertainty within ±20% for all but the highest temperature cases, which have a maximum uncertainty of ±40%. Experiments conducted with a radical trapper, toluene, showed little influence from radical chemistry, suggesting this estimated uncertainty is fairly conservative. Experimental data from Heyne et al (Z Phys Chem. 2015;229:881-907) were found to be in good agreement with the rate measurements from this study and, therefore, a second Arrhenius expression, k 1 (s –1 ) = 2.11 × 10 13 exp(-30 820 K/T), was found to represent both datasets well. This second expression has a larger temperature range of 976-1621 K. Here the present study provides the first high-temperature data collected for this reaction, adding to the limited data available for isopropanol in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1974

Preliminary performance estimates of a highly maneuverable remotely piloted vehicle

A computerized synthesis program has been used to assess the effects of various vehicle and mission parameters on the performance of a highly maneuverable remotely piloted vehicle (RPV) for the air-to-air combat role. The configuration used in the study is a trapezoidal-wing and body concept, with forward-mounted stabilizing and control surfaces. The study mission consists of an outbound cruise, an acceleration phase, a series of subsonic and supersonic turns, and a return cruise. Performance is evaluated in terms of both the required vehicle weight to accomplish this mission and combat effectiveness as measured by turning and acceleration capability. The report describes the synthesis program, the mission, the vehicle, and the results of sensitivity and trade studies.

Nelms, W. P., Jr.↗

NASA NTRS2000

Test Cases for a Rectangular Supercritical Wing Undergoing Pitching Oscillations

Steady and unsteady measured pressures for a Rectangular Supercritical Wing (RSW) undergoing pitching oscillations have been presented. From the several hundred compiled data points, 27 static and 36 pitching oscillation cases have been proposed for computational Test Cases to illustrate the trends with Mach number, reduced frequency, and angle of attack. The wing was designed to be a simple configuration for Computational Fluid Dynamics (CFD) comparisons. The wing had an unswept rectangular planform plus a tip of revolution, a panel aspect ratio of 2.0, a twelve per cent thick supercritical airfoil section, and no twist. The model was tested over a wide range of Mach numbers, from 0.27 to 0.90, corresponding to low subsonic flows up to strong transonic flows. The higher Mach numbers are well beyond the design Mach number such as might be required for flutter verification beyond cruise conditions. The pitching oscillations covered a broad range of reduced frequencies. Some early calculations for this wing are given for lifting pressure as calculated from a linear lifting surface program and from a transonic small perturbation program. The unsteady results were given primarily for a mild transonic condition at M = 0.70. For these cases the agreement with the data was only fair, possibly resulting from the omission of viscous effects. Supercritical airfoil sections are known to be sensitive to viscous effects (for example, one case cited). Calculations using a higher level code with the full potential equations have been presented for one of the same cases, and with the Euler equations. The agreement around the leading edge was improved, but overall the agreement was not completely satisfactory. Typically for low-aspect-ratio rectangular wings, transonic shock waves on the wing tend to sweep forward from root to tip such that there are strong three-dimensional effects. It might also be noted that for most of the test, the model was tested with free transition, but a few points were taken with an added transition strip for comparison. Some unpublished results of a rigid wing of the same airfoil and planform that was tested on the pitch and plunge apparatus mount system (PAPA) showed effects of the lower surface transition Strip on flutter at the lower subsonic Mach numbers. Significant effects of a transition strip were also obtained on a wing with a thicker supercritical section on the PAPA mount system. Both of these flutter tests on the PAPA resulted in very low reduced frequencies that may be a factor in this influence of the transition strip. However, these results indicate that correlation studies for RSW may require some attention to the estimation of transition location to accurately treat viscous effects. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. An overview of the model and tests is given and the standard formulary for these data is listed. Sample data points are presented in both tabular and graphical form. A complete tabulation and plotting of all the Test Cases is given. Only the static pressures and the real and imaginary parts of the first harmonic of the unsteady pressures are available. All the data for the test are available in electronic file form. The Test Cases are also available as separate electronic files.

Bennett, Robert M.↗

NASA NTRS2019

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach Using SMOS L-Band Microwave Data

The uncertainty of soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several (quasi-)operational and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different literature-based RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters, and that the use of multi-angular and multi-polarization TB does not necessarily lead to superior retrievals compared to retrievals based on TB data at a single incidence angle and polarization. The second part of this paper evaluates ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM. As expected for bounded variables, the spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term mean bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread of 0.037 cu m/cu m approximates the actual time series unbiased root-mean-square-difference of 0.042 cu m/cu m between ensemble mean retrievals and in situ reference data across the reference sites. However, the temporal variability in the anomaly ensemble spread reveals higher-order biases in the retrieval error, which should be accounted for when characterizing retrieval error.

Jan Quets↗

DOE OSTI2021

Practical considerations for amino acid isotope analysis

Over the last few decades, isotopic analysis of amino acids at the compound- and position-specific levels has been rapidly advancing across diverse fields. As these techniques progress, evaluation of isotopic fractionation associated with sample workup is essential. This critical review of analytical methods through the lens of isotope geochemistry provides a benchmark for researchers across disciplines seeking to make compound- and position-specific amino acid isotope measurements. We focus on preparation, acid hydrolysis, clean-up, derivatization, separation, and C, H, N, and S isotope measurement. Despite substantial customizability across these steps, the following general recommendations should maximize recovery while minimizing isotopic fractionation. Samples should be freeze-dried and stored anoxically at ≤ –20 °C prior to conventional acid hydrolysis (6N HCl, 110 °C, 20–24 h, anoxic), which suffices for many residues. Both gas and liquid chromatographic (GC and LC, respectively) techniques are well-established and separate about 15 amino acids; LC bypasses the need for derivatization, while GC provides higher sensitivity. Furthermore, when derivatization is needed, n-acetyl and alkoxycarbonyl esters provide the most reproducible C isotope ratios. For compound-specific analyses, online GC–IRMS and LC–IRMS systems offer the easiest workflow, but EA–IRMS enables potential multi-element isotope analysis. Emerging techniques like high-resolution mass spectrometry are also promising for multi-element analysis and recover position-specific isotopic information. Looking forward to the next decade of innovation, isotope geochemists and ecologists can improve amino acid isotope analysis by focusing on streamlining multi-element analysis and standardizing calibration practices across laboratories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Evaluating spectral cloud effective radius retrievals from the Enhanced MODIS Airborne Simulator (eMAS) during ORACLES

Satellite remote sensing retrievals of cloud effective radius (CER) are widely used for studies of aerosol–cloud interactions. Such retrievals, however, rely on forward radiative transfer (RT) calculations using simplified assumptions that can lead to retrieval errors when the real atmosphere deviates from the forward model. Here, coincident airborne remote sensing and in situ observations obtained during NASA's ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) field campaign are used to evaluate retrievals of CER for marine boundary layer stratocumulus clouds and to explore impacts of forward RT model assumptions and other confounding factors. Specifically, spectral CER retrievals from the Enhanced MODIS Airborne Simulator (eMAS) and the Research Scanning Polarimeter (RSP) are compared with polarimetric retrievals from RSP and with CER derived from droplet size distributions (DSDs) observed by the Phase Doppler Interferometer (PDI) and a combination of the Cloud and Aerosol Spectrometer (CAS) and the Two-Dimensional Stereo Probe (2D-S). The sensitivities of the eMAS and RSP spectral retrievals to assumptions about the DSD effective variance (CEV) and liquid water complex index of refraction are explored. CER and CEV inferred from eMAS spectral reflectance observations of the backscatter glory provide additional context for the spectral CER retrievals. The spectral and polarimetric CER retrieval agreement is case dependent, and updating the retrieval RT assumptions, including using RSP polarimetric CEV retrievals as a constraint, yields mixed results that are tied to differing sensitivities to vertical heterogeneity. Moreover, the in situ cloud probes, often used as the benchmark for remote sensing CER retrieval assessments, themselves do not agree, with PDI DSDs yielding CER values 1.3–1.6 µm larger than CAS and with CEV roughly 50 %–60 % smaller than CAS. Implications for the interpretation of spectral and polarimetric CER retrievals and their agreement are discussed.

Meyer, Kerry [NASA Goddard Space Flight Center (GS↗

NASA NTRS2011

Standardization of Spore Inactivation Method for PMA-PhyloChip Analysis

In compliance with the Committee on Space Research (COSPAR) planetary protection policy, National Aeronautics and Space Administration (NASA) monitors the total microbial burden of spacecraft as a means for minimizing the inadvertent transfer of viable contaminant microorganisms to extraterrestrial environments (forward contamination). NASA standard assay-based counts are used both as a proxy for relative surface cleanliness and to estimate overall microbial burden as well as to assess whether forward planetary protection risk criteria are met for a given mission, which vary by the planetary body to be explored and whether or not life detection missions are present. Despite efforts to reduce presence of microorganisms from spacecraft prior to launch, microbes have been isolated from spacecraft and associated surfaces within the extreme conditions of clean room facilities using state of the art molecular technologies. Development of a more sensitive method that will better enumerate all viable microorganisms from spacecraft and associated surfaces could support future life detection missions. Current culture-based (NASA standard spore assay) and nucleic-acid-based polymerase chain reaction (PCR) methods have significant shortcomings in this type of analysis. The overall goal of this project is to evaluate and validate a new molecular method based on the use of a deoxyribonucleic acid (DNA) intercalating agent propidium monoazide (PMA). This is used in combination with DNA microarray (PhyloChip) which has been shown to identify very low levels of organisms on spacecraft associated surfaces. PMA can only penetrate the membrane of dead cells. Once penetrated, it intercalates the DNA and, upon photolysis using visible light it produces stable DNA monoadducts. This allows DNA to be unavailable for further PCR analysis. The specific aim of this study is to standardize the spore inactivation method for PMA-PhyloChip analysis. We have used the bacterial spores Bacillus subtilis 168 (standard laboratory isolate) as a test organism.

spore inactivation↗

DOE OSTI2024

Review of data-driven models for quantifying load shed by non-residential buildings in the United States

Shifting and shedding power demand in buildings can be cost-effective techniques for grids to function reliably and for end users to earn compensation. Grid operators reimburse customers in proportion to the quantity of load shed. Simple data-driven methods are used to quantify this shed, which is the difference between a measured load during the event and modeled "baseline" that would have occurred in absence of the event. These methods have evolved over the years and in many cases have been integrated with building physics, to make them a hybrid between physics based and empirical models. However, there is no comprehensive analysis that provides guidance to building operators, grid operators and researchers in selecting appropriate models based on their specific needs and available data. Here, this work aims to fill this gap by critically assessing the performance of baseline models put forward from the year 2000 through 2023. The literature reviewed includes reports generated by grid operators, reports from national laboratories and academic journal articles. The work outlines modeling features like the inputs, training period, estimation method, adjustments to fine tune the predictions and metrics to evaluate the performance. A comprehensive list of 50 models has been provided. For each model, the study explores the applicability of the model to weather sensitive buildings, variability in the building profile, timing of the event, and whether the building reduces energy consumption before an event. The work identifies the situations in which a particular model works and draws lessons based on evidence of performance. Finally, recommendations to aid in model selection are given.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Investigation of Correlation Methods for Use in Criticality Safety

Although their adoption by practitioners has been limited, the introduction of similarity indices in criticality safety was a major step forward in reducing the reliance on expert judgement in discerning applicable experiments for the validation of new appliations in criticality safety analyses. Similarity indices have been successfully employed in bias trending and data assimilation techniques, but it is often unclear which acceptance criteria should be used. In their 2004 paper, Broadhead et al. specify the most widely used similarity parameter, ck, as an acceptance cutoff at 0.9. (Broadhead et al., ”Sensitivity and Uncertainty-Based Criticality Safety Validation Techniques,” Nucl. Sci. Eng. 146, 340–366, 2004). Experiments with a ck < 0.9 are often not considered applicable for code validation. This heuristic is based on quantitative studies and engineering judgement, but in some cases, experiments with ck < 0.9 can be used to accurately estimate computational bias. This suggests that further analysis is needed to determine what components of ck are driving applicability and accuracy in bias estimation. For cases in which applicable experiments may not be available (as is the case with UF6 transport canisters), understanding what distinguishes experiments in providing adequate bias estimates aside from just the similarity index is particularly necessary. To further the goal to better interpret ck values, several visualization tools were developed to assist in the investigation of which components of ck are driving applicability.

ck↗

NASA NTRS2009

Hardware Implementation of Serially Concatenated PPM Decoder

A prototype decoder for a serially concatenated pulse position modulation (SCPPM) code has been implemented in a field-programmable gate array (FPGA). At the time of this reporting, this is the first known hardware SCPPM decoder. The SCPPM coding scheme, conceived for free-space optical communications with both deep-space and terrestrial applications in mind, is an improvement of several dB over the conventional Reed-Solomon PPM scheme. The design of the FPGA SCPPM decoder is based on a turbo decoding algorithm that requires relatively low computational complexity while delivering error-rate performance within approximately 1 dB of channel capacity. The SCPPM encoder consists of an outer convolutional encoder, an interleaver, an accumulator, and an inner modulation encoder (more precisely, a mapping of bits to PPM symbols). Each code is describable by a trellis (a finite directed graph). The SCPPM decoder consists of an inner soft-in-soft-out (SISO) module, a de-interleaver, an outer SISO module, and an interleaver connected in a loop (see figure). Each SISO module applies the Bahl-Cocke-Jelinek-Raviv (BCJR) algorithm to compute a-posteriori bit log-likelihood ratios (LLRs) from apriori LLRs by traversing the code trellis in forward and backward directions. The SISO modules iteratively refine the LLRs by passing the estimates between one another much like the working of a turbine engine. Extrinsic information (the difference between the a-posteriori and a-priori LLRs) is exchanged rather than the a-posteriori LLRs to minimize undesired feedback. All computations are performed in the logarithmic domain, wherein multiplications are translated into additions, thereby reducing complexity and sensitivity to fixed-point implementation roundoff errors. To lower the required memory for storing channel likelihood data and the amounts of data transfer between the decoder and the receiver, one can discard the majority of channel likelihoods, using only the remainder in operation of the decoder. This is accomplished in the receiver by transmitting only a subset consisting of the likelihoods that correspond to time slots containing the largest numbers of observed photons during each PPM symbol period. The assumed number of observed photons in the remaining time slots is set to the mean of a noise slot. In low background noise, the selection of a small subset in this manner results in only negligible loss. Other features of the decoder design to reduce complexity and increase speed include (1) quantization of metrics in an efficient procedure chosen to incur no more than a small performance loss and (2) the use of the max-star function that allows sum of exponentials to be computed by simple operations that involve only an addition, a subtraction, and a table lookup. Another prominent feature of the design is a provision for access to interleaver and de-interleaver memory in a single clock cycle, eliminating the multiple clock-cycle latency characteristic of prior interleaver and de-interleaver designs.

Moision, Bruce↗