Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Confidence interval”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Role of emission controls in reducing the 2050 climate change penalty for PM2.5 in China

Previous studies demonstrated that global warming can lead to deteriorated air quality even when anthropogenic emissions were kept constant, which has been called a climate change penalty on air quality. It is expected that anthropogenic emissions will decrease significantly in the future considering the aggressive emission control actions in China. However, the dependence of climate change penalty on the choice of emission scenario is still uncertain. To fill this gap, we conducted multiple independent model simulations to investigate the response of PM2.5 to future (2050) climate warming (RCP8.5) in China but with different emission scenarios, including the constant 2015 emissions, the 2050 CLE emissions (based on Current Legislation), and the 2050 MTFR emissions (based on Maximum Technically Feasible Reduction). For each set of emissions, we estimate climate change penalty as the difference in PM2.5 between a pair of simulations with either 2015 or 2050 meteorology. Under 2015 emissions, we find a PM2.5 climate change penalty of 1.43 μg m−3 in Eastern China, leading to an additional 35,000 PM2.5-related premature deaths [95% confidence interval (CI), 21,000-40,000] by 2050. However, the PM2.5 climate change penalty weakens to 0.24 μg m−3 with strict anthropogenic emission controls under the 2050 MTFR emissions, which decreases the associated PM2.5-related deaths to 17,000. The smaller MTFR climate change penalty contributes 14% of the total PM2.5 decrease when both emissions and meteorology are changed from 2015 to 2050, and 24% of total health benefits associated with this PM2.5 decrease in Eastern China. This finding suggests that controlling anthropogenic emissions can effectively reduce the climate change penalty on PM2.5 and its associated premature deaths, even though a climate change penalty still occurs even under MTFR. Strengthened controls on anthropogenic emissions are key to attaining air quality targets and protecting human health in the context of future global climate change.

Emission controls↗
NASA NTRSDate not supplied

Congruence of Clusters Defined By Whole Genome Sequencing and MALDI-TOF for Bacteria Isolated From Cleanrooms

Introduction: Oligotrophic conditions can render cleanrooms inhospitable to microbes. Despite these constraints, fungi and bacteria are frequently isolated from surfaces in astromaterials cleanrooms at the Johnson Space Center. Bacillus species are of particular concern because endospores belonging to this genus are resilient and can affect astromaterials. Current monitoring programs rely on 16S rRNA sequencing and the VITEK2 Compact system. These methods have limited power to resolve Bacillus species. Matrix-assisted laser desorption - time of flight mass spectrometry (MALDI-TOF MS), provides a rapid, low cost, method of identifying bacterial isolates and has a higher resolution than 16S rRNA sequencing, particularly for Bacillus species; however, few studies have compared this method to the industry gold standard, whole genome sequencing (WGS). Methods: Based on 16S rRNA classification, we selected 14 isolates for analysis with MALDI-TOF and WGS. Mass spectra were generated with MALDI-TOF MS and processed with custom scripts to identify clusters of closely related isolates and calculate a matrix of pairwise cosine similarity scores. Hybrid Illumina and Nanopore sequencing were used to generate draft genomes. Pairwise similarity scores were calculated from these genomes based on the average amino acid identity (AAI) predicted from single copy core genes. Congruence of clustering between these methods, was assessed by calculating adjusted Rand and Wallace coefficients. Results: Clusters of species generated from MALDI-TOF MS showed good agreement of phylotypes generated with WGS. Pairs of strains that were > 94% similar to each other in terms of predicted amino acid sequences consistently showed cosine similarities of mass spectra > 0.65 and, of the 9 clusters identified with WGS, 8 were identical with MALDI-TOF. This corresponds to an adjusted Rand index of 0.95 and a 95% confidence interval of 0.80 – 1.00 for adjusted Wallace coefficients. The only discordance was for a pair of isolates that were classified as Paenibacillus species. This pair showed relatively high similarity (0.84) in terms of MALDI-TOF MS but only 85% similarity in terms of AAI. Conclusion: This study shows that MALDI-TOF and WGS exhibit a similar ability to delineate Bacillus species isolated from cleanrooms and taxonomic units described by these two methods are consistent with one another. Since MALDI-TOF MS is low in cost and high in throughput, this approach appears to be an ideal option for routine microbial monitoring and identifying Bacillus species.

Farnaz Mazhari↗
NASA NTRS2021

Development and Validation of an Empirical Ocean Color Algorithm with Uncertainties: A Case Study with the Particulate Backscattering Coefficient

We explored how algorithm (model) and in situ measurement (observation) uncertainties can effectively be incorporated into empirical ocean color model development and assessment. In this study we focused on methods for deriving the particulate backscattering coefficient at 555 nm, b(bp)(555)/(m). We developed a simple empirical algorithm for deriving b(bp)(555) as a function of a remote sensing reflectance line height (LH) metric. Model training was performed using a high-quality bio-optical dataset that contains coincident in situ measurements of the spectral remote sensing reflectances, R(rs)(λ)/(sr), and the spectral particulate backscattering coefficients, b(bp)(λ). The LH metric used is defined as the magnitude of Rrs(555) relative to a linear baseline drawn between R(rs)(490) and R(rs)(670). Using an independent validation dataset, we compared the skill of the LH-based model with two other models. We used contemporary validation metrics, including bias and mean absolute error (MAE), that were corrected for model and observation uncertainties. The results demonstrated that measurement uncertainties do indeed impact contemporary validation metrics such as mean bias and MAE. Zeta-scores and z-tests for overlapping confidence intervals were also explored as potential methods for assessing model skill.

ocean color↗
NASA NTRS2019

GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observedby LIGO and Virgo During the First and Second Observing Runs

We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above1M during the first and second observing runs of the advanced gravitational-wave detector network. During the first observing run (O1), from September 12, 2015 to January 19, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run (O2), which ran from November 30, 2016 to August 25, 2017, saw the first detection of gravitational waves from a binary neutron star in spiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818, and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between18.6þ3.2−0.7M and 84.4þ15.8−11.1M and range in distance between 320þ120−110 and 2840þ1400−1360Mpc. No neutron star–black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false-alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of 110−3840Gpc−3y−1 for binary neutron stars and 9.7−101Gpc−3y−1 for binary black holes assuming fixed population distributions and determine a neutron star–black hole merger rate 90% upper limit of 610Gpc−3y−1.

Astrophysics↗
NASA NTRSDate not supplied

Assessment of Crew Time for Maintenance and Repairs Activities for Lunar Surface Missions

NASA is currently evaluating different methods to predict how much time crewmembers will spend conducting repair and maintenance activities on future space missions. As mission scope and spacecraft architectures change, it will be necessary to understand how crew repair and maintenance timelines are impacted by mission operations and technology changes. Past work has been done using historical ISS data to accurately predict crew habitation and operation timelines, resulting in the development of NASA’s Exploration Crew Time Model (ECTM). However, understanding crew maintenance and repair requirements has posed a unique challenge due to the complexity of available datasets, the probabilistic nature of sub-system failures, and the impacts of reliability growth on failure rates. This paper presents a methodology to collect and condition empirical repair and maintenance time data from available data sets, to extrapolate from that data to estimate projected maintenance and repair times for a lunar Surface Habitat, and to assess how uncertainty in repair time could impact utilization time on the lunar surface. NASA International Space Station (ISS) maintenance and crew time data are logged into two central databases, the Maintenance Data Collection (MDC) and the Operations Planning Timeline Integration System (OPTimIS) respectively. Separately, each of these two datasets capture only portions of the complete set of data required to generate an accurate assessment of crew time spent on maintenance activities at a sub-system level. MDC provides a detailed catalog of failure events and an overview of the failure’s required maintenance and OPTimIS provides a description of crew activities and crew time durations dedicated to maintenance. To create a more useful crew time estimate for maintenance timelines, the authors developed a methodology to capture relevant data from each set and combine and utilize that data by linking crew time requirements to specific components. The authors compare the failure logs in the MDC to crew activity logs pulled from OPTimIS and then process the data to estimate required repair times for each failure event. Data is also classified by the outcome of each repair event, whether the failed component was replaced or whether it was repaired in place. The entire maintenance activity dataset is then categorized based on the class of failed component to allow for a statistically significant sample size for each class and to provide accurate crew time estimates for any components lacking relevant data. This resultant component repair time data can be used in the future to generate Mean Time To Repair (MTTR) estimates and confidence intervals for each class of component based on a probabilistic distribution of documented maintenance events. These improved MTTR values can then be applied to candidate element sub-system architectures, along with component Mean Time Between Failure (MTBF) data to generate distributions for potential required system crew repair time estimates for a given mission. Repair time distributions can then be used to develop more accurate crew schedules and to assess potential available utilization time.

Crew Time↗
NASA NTRSDate not supplied

Analysis of Pull Force Test Results for Crimped Connections

Crimped electrical contacts are reliable when strict process controls are followed during manufacturing and accompanied by continuous process verification through pull force testing. Cable and wire harness assemblies’ standards are developed and refined over time to provide the minimum pull force that a crimp contact must meet before it breaks from the wire. However, in practice the failures occur at a much higher tensile strength than the minimum required. The first section of this paper reviewed 780 pull force test results provided by NASA Centers that were collected and analyzed to determine how the data compare to NASA’s pre-existing requirements from the cable/harness standards NASA-STD-8739.4 and IPC/WHMA-A-620B-S. The measured tensile strength of most of the contact/conductor pairs (i.e., Contact/Wires or C/W) exceeded the minimum pull force values of NASA-STD-8739.4 and IPC/WHMA-A-620 standards by at least 100%. The C/W pair samples’ tensile strength values followed a normal distribution with an average tensile strength value that was at least 182% of the minimum requirement, and all the samples analyzed passed the pull force testing. In addition, the 95% confidence interval of the average tensile strength distributions for several C/W pairs was determined and plotted as error bars to show that the C/W pairs will meet and surpass the requirements. The frequency of pull force testing can be problematic for projects because of the cost and availability of spare contacts for the destructive test. It is possible to reduce the frequency of pull force testing if at the beginning of the production run, the conditions of the crimp tool and materials are verified, and the settings of the tool remain unchanged throughout the process. However, the project needs to understand and evaluate the impact to risk from reducing the frequency of testing prior to implementing process changes.

Alejandra Constante↗
NASA NTRSDate not supplied

Assessment of Crew Time for Maintenance and Repair Activities for Lunar Surface Missions

NASA is currently evaluating different methods to predict how much time crewmembers will spend conducting repair and maintenance activities on future space missions. As mission scope and spacecraft architectures change, understanding how crew repair and maintenance timelines are impacted by mission operations and technology changes is vital for future mission planning. Past work has been done using historical International Space Station (ISS) data to accurately predict crew habitation and operation timelines, resulting in the development of NASA’s Exploration Crew Time Model (ECTM). However, understanding crew maintenance and repair requirements has posed a unique challenge due to the complexity of available datasets, the probabilistic nature of sub-system failures, and the impacts of reliability growth on failure rates. This paper presents a methodology to collect and condition empirical repair and maintenance time data from available datasets, to extrapolate from that data to estimate projected maintenance and repair times for a lunar Surface Habitat (SH), and to assess how uncertainty in repair time could impact utilization time on the lunar surface. NASA ISS maintenance and crew time data are logged into two central databases: the Maintenance Data Collection (MDC) and the Operations Planning Timeline Integration System (OPTimIS). Separately, each of these two datasets capture only portions of the complete set of data required to generate an accurate assessment of crew time spent on maintenance activities at a sub-system level. To create a more useful crew time estimate for maintenance timelines, the authors developed a methodology to capture relevant data from each set and combine and utilize that data by linking crew time requirements to specific components. The authors compare the failure logs in the MDC to crew activity logs pulled from OPTimIS and then process the data to estimate required repair time for each failure and repair event. The entire maintenance activity dataset is then categorized based on the class of failed component to ensure a significant sample size for each class and accurate crew time estimates for any components lacking relevant data. This resultant component repair time data can be used in the future to generate Mean Time to Repair (MTTR) estimates and confidence intervals for each class of component based on a probabilistic distribution of documented maintenance events. These improved MTTR values can then be applied to candidate element sub-system architectures, along with component Mean Time Between Failure (MTBF) data to generate distributions for potential required system crew repair time estimates for a given mission. The authors applied these modeling methods to a case study of a crewed mission to the planned SH and produced expected corrective maintenance crew time distributions. The results produced an expected corrective maintenance crew time at over 24 hours per mission, and a maintenance crew time distribution that reflects the importance of planning for sufficient maintenance requirements each mission. Repair time distributions can then be used to develop more accurate crew schedules and to assess potential available utilization time.

Crew Time↗
NASA NTRSDate not supplied

Variance Decomposition of MEDLI2 Reconstructed Heating Using Neural Networks

The Mars Entry, Descent, and Landing Instrumentation (MEDLI2) sensor suite collected data during entry of the Mars 2020 Perseverance rover into Mars’ atmosphere. This suite included a network of MEDLI2 Instrumented Sensor Plugs (MISPs). Each MISP was comprised of a cylinder made of Thermal Protection System (TPS) material with 1-3 embedded thermocouples (TCs), and it was flush mounted into the heatshield or backshell. Data from these in-depth TCs were used to reconstruct the aeroheating environment of the vehicle throughout entry. Surface heating was posed as an inverse problem, with the goal of estimating the surface heating by minimizing an objective function of the difference between MISP temperature measurements during flight and the temperature predictions derived from the Fully Implicit Ablation and Thermal response (FIAT) program. Given an aerothermal environment, FIAT calculates the material response and provides in-depth temperatures throughout the TPS material. To achieve the reverse, an internal tool called FIAT_Opt runs through multiple different environments until the output temperature at the TC depth closely matches the flight data. 95% confidence intervals on the reconstructed surface heating were obtained using Monte Carlo analysis, in which uncertainties in the thermocouple depth and the TPS material properties (e.g., density, thermal conductivity, heat capacity, emissivity) based on flight-lot material testing were included. A variance decomposition method using Sobol indices was employed to assess the sensitivity of the reconstructed peak heating to the TC placement and material property uncertainties. Variance decomposition was found to require tens of thousands of FIAT_Opt runs in order for the Sobol indices to converge. With a single FIAT_Opt run taking on the order of 40 minutes, the required number of computations would take months to complete, even if using multiple CPUs. To mitigate this problem, three machine learning models (ridge regression with cross-validation, random forest regression, and a deep neural network) were trained and tested using the 2000 Monte Carlo runs that were already completed. A subset of 1600 runs were used to train the model (i.e., training set), while the remaining 400 runs were used as the test set. The predictions from the deep neural network (DNN) on the test set showed nearly perfect agreement to the actual values computed with FIAT_Opt (R2 > 0.99). Using the DNN as a surrogate model, the variance decomposition using 50,000 runs was completed within minutes. The resulting Sobol indices showed that the reconstructed peak surface heating was most sensitive to the uncertainties in the thermal conductivity (ST = 0.37) and heat capacity (ST = 0.26). This method can be leveraged to provide requirements for material property measurements needed to improve the accuracy of surface heating prediction and ultimately lead to the reduction of design margins in the future. This presentation will include background on the MEDLI2 suite; the method used for inverse heating estimation; the way that material property uncertainties were accounted for using Monte Carlo analysis; a brief background on variance decomposition; the motivation for using machine learning in this context; how a neural network was trained on the data to enable variance decomposition in a fraction of the time; and the variance decomposition results for one of the MISPs.

Hannah Alpert↗
NASA NTRS2021

A Uniform Search for Nearby Planetary Companions to Hot Jupiters in TESS Data Reveals Hot Jupiters Are Still Lonely

We present the results of a uniform search for additional planets around all stars with confirmed hot Jupiters observed by the Transiting Exoplanet Survey Satellite (TESS) in its Cycle 1 survey of the southern ecliptic hemisphere. Our search comprises 184 total planetary systems with confirmed hot Jupiters with Rp > 8 R⊕ and orbital period <10 days. The Transit Least Squares algorithm was utilized to search for periodic signals that may have been missed by other planet search pipelines. While we recovered 169 of these confirmed hot Jupiters, our search yielded no new statistically validated planetary candidates in the parameter space searched (P < 14 days). A lack of planet candidates nearby hot Jupiters in the TESS data supports results from previous transit searches of each individual system, now down to the photometric precision of TESS. This is consistent with expectations from a high-eccentricity migration formation scenario, but additional formation indicators are needed for definitive confirmation. We injected transit signals into the light curves of the hot Jupiter sample to probe the pipeline's sensitivity to the target parameter space, finding a dependence proportional to ${R}_{p}^{2.32}{P}^{-0.88}$ for planets within 0.3 ≤ Rp ≤ 4 R⊕ and 1 ≤ P ≤ 14 days. A statistical analysis accounting for this sensitivity provides a median and 90% confidence interval of ${7.3}_{-7.3}^{+15.2} \% $ for the rate of hot Jupiters with nearby companions in this target parameter space. This study demonstrates how TESS uniquely enables comprehensive searches for nearby planetary companions to nearly all the known hot Jupiters.

Benjamin J Hord↗
NASA NTRS2022

Global Hotspots of Salt Marsh Change and Carbon Emissions

Salt marshes provide ecosystem services such as carbon sequestration, coastal protection, sea-level-rise (SLR) adaptation and recreation. SLR, storm events, drainage and mangrove encroachment are known drivers of salt marsh loss. However, the global magnitude and location of changes in salt marsh extent remains uncertain. Here we conduct a global and systematic change analysis of Landsat satellite imagery from the years 2000–2019 to quantify the loss, gain and recovery of salt marsh ecosystems and then estimate the impact of these changes on blue carbon stocks. We show a net salt marsh loss globally, equivalent to an area double the size of Singapore (719 km^2), with a loss rate of 0.28% year^(−1) from 2000 to 2019. Net global losses resulted in 16.3 (0.4–33.2, 90% confidence interval) Tg CO_2 e year^(−1) emissions from 2000 to 2019 and a 0.045 (−0.14–0.115) Tg CO_2 e year^(−1) reduction of carbon burial. Russia and the USA accounted for 64% of salt marsh losses, driven by hurricanes and coastal erosion. Our findings highlight the vulnerability of salt marsh systems to climatic changes such as SLR and intensification of storms and cyclones.

Carbon Monitoring↗
NASA NTRSDate not supplied

Component-level Performance and Mass Sensitivity Analysis of NEP MW-class Power System

Nuclear electric propulsion (NEP) is a promising option towards enabling missions to Mars and is an area of interest for NASA’s Space Nuclear Propulsion project. This project is currently investigating technology development opportunities for an NEP vehicle. Physics-based modeling can be used in the early stages of technology development to gain understanding of the effects of technology and performance assumptions on the system performance and mass. This information can then inform technology maturation planning for near term development. By using a Brayton power conversion model and vehicle mass model for megawatt class NEP applications, a sensitivity analysis is performed to assess the impact of individual components’ performance on the power conversion system performance and system mass. A Monte Carlo simulation is also used to determine the variability in system mass based on uncertainty within the modeling parameters. The sensitivity analysis shows a high sensitivity to power conversion inlet temperature, compressor inlet temperature, and recuperator performance. A Monte Carlo analysis suggests a range of -10% to +15% for a 90% confidence interval on system mass based on the uncertainties in the model inputs.

Nuclear electric propulsion↗
NASA NTRS2019

Simultaneous X-Ray and Infrared Observations of Sagittarius A*'s Variability

Emission from Sagittarius A* is highly variable at both X-ray and infrared (IR) wavelengths. Observations over the last ∼20 yr have revealed X-ray flares that rise above a quiescent thermal background about once per day, while faint X-ray flares from Sgr A* are undetectable below the constant thermal emission. In contrast, the IR emission of Sgr A* is observed to be continuously variable. Recently, simultaneous observations have indicated a rise in IR flux density around the same time as every distinct X-ray flare, while the opposite is not always true (peaks in the IR emission may not be coincident with an X-ray flare). Characterizing the behavior of these simultaneous X-ray/IR events and measuring any time lag between them can constrain models of Sgr A* ʼs accretion flow and the flare emission mechanism. Using 100+ hours of data from a coordinated campaign between the Spitzer Space Telescope and the Chandra X-ray Observatory, we present results of the longest simultaneous IR and X-ray observations of Sgr A* taken to date. The cross-correlation between the IR and X-ray light curves in this unprecedented data set, which includes four modest X-ray/IR flares, indicates that flaring in the X-ray may lead the IR by approximately 10–20 min with 68% confidence. However, the 99.7% confidence interval on the time-lag also includes zero, i.e., the flaring remains statistically consistent with simultaneity. Long-duration and simultaneous multi-wavelength observations of additional bright flares will improve our ability to constrain the flare timing characteristics and emission mechanisms, and must be a priority for Galactic Center observing campaigns.

H. Boyce↗
NASA NTRS2022

Global Premature Mortality By Dust and Pollution PM 2.5 Estimated From Aerosol Reanalysis of the Modern-Era Retrospective Analysis for Research and Applications, Version 2

This study quantifies global premature deaths attributable to long-term exposure of ambient PM 2.5 , or PM 2.5 -attributable mortality, by dust and pollution sources. We used NASA’s Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) aerosol reanalysis product for PM 2.5 and the cause-specific relative risk (RR) from the integrated exposure-response (IER) model to estimate global PM2.5-attributable mortality for five causes of deaths, namely ischaemic heart disease (IHD), cerebrovascular disease (CEV) or stroke, lung cancer (LC), chronic obstructive pulmonary disease (COPD), and acute lower respiratory infection (ALRI). The estimated yearly global PM 2.5 -attributable mortality in 2019 amounts to 2.89 (1.38–4.48) millions, which is composed of 1.19 (0.73–1.84) million from IHD, 1.01 (0.35–1.55) million from CEV, 0.29 (0.11–0.48) million from COPD, 0.23 (0.14–0.33) million from ALRI, and 0.17 (0.04–0.28) million from LC (the numbers in parentheses represent the estimated mortality range due corresponding to RR spread at the 95% confidence interval). The mortality counts vary with geopolitical regions substantially, with the highest number of deaths occurring in Asia. China and India account for 40% and 23% of the global PM 2.5 -attributable deaths, respectively. In terms of sources of PM 2.5 , about 22% of the global all-cause PM 2.5 -attributable deaths are caused by desert dust. The largest dust attribution is 37% for ALRI. The relative contributions of dust and pollution sources vary with the causes of deaths and geographical regions. Enforcing air pollution regulations to transfer areas from PM 2.5 nonattainment to PM2.5 attainment can have great health benefits. Being attainable with the United States air quality standard (AQS) of 15 μg/m 3 globally would have avoided nearly 40% or 1.2 million premature deaths. The most recent update of PM 2.5 guideline from 10 to 5 μg/m 3 by the World Health Organization (WHO) would potentially save additional one million lives. Our study highlights the importance of distinguishing aerodynamic size from geometric size in accurately assessing the global health burden of PM 2.5 and particularly for dust. A use of geometric size in diagnosing dust PM 2.5 from the model simulation, a common approach in current health burden assessment, could overestimate the PM 2.5 level in the dust belt by 40–170%, leading to an overestimate of global all-cause mortality by 1 million or 32%.

PM2.5↗
NASA NTRS2022

Bayesian Rules of Thumb: Robust Uncertainty Quantification in Early Project Cost Estimation

Systems engineers often make use of cost Rules ofThumb in order to estimate cost during early phases of projectformulation. These Rules of Thumb typically take the form ofa sequence of percentages over which a total cost is allocatedacross NASA WBS elements. Rules of Thumb can then be usedto extrapolate cost from one or more known WBS elements tothe remaining unknown WBS elements, assisting early projectformulation architecture studies (such as those in JPL’s Team Xand A Team).A number of issues can arise when generating and using costRules of Thumb. For example, many records of project costsconsist of incomplete data. Typical methods of dealing withincomplete cost allocation data include (a) ignoring missionswith incomplete data, or (b) taking averages of the non-zero percentagesacross missions, but both of these methods can result inbiased estimates if the existence of incomplete data correlateswith total mission cost or any particular WBS element. Anothercommon example is cost reported in one or more incorrect WBSelements. This is especially prevalent in smaller missions whereit is more common for engineers to perform tasks that fall underthe purview of multiple WBS elements.Furthermore, a Rule of Thumb estimate is typically reported asa point estimate; there is no reported uncertainty around thepercentages used to generate an allocation. Even in the rarecase in which confidence intervals around mean percentages areprovided, there may be positive or negative correlations betweenWBS elements which can skew estimates.Here we attempt to address these problems by formulatingprobabilistic Rules of Thumb in which a distribution of allocationschemes, rather than a single allocation scheme, is generated.We use a bootstrap imputation method to simultaneouslyaccount for uncertainty in the missing data while using allavailable information contained in the dataset. The imputeddatasets are then input into a multivariate Bayesian modelwhich accounts for correlations between WBS elements andproperly accounts for uncertainty in the final Rule of Thumbpercentages and predictions. We describe the mathematicalmodel and provides snippets of R code utilizing the brms(Bayesian Regression Models using Stan) package. To illustratethis model, we generate a Bayesian Level 2 WBS Cost Rule ofThumb for MIDEX (Medium-Class Explorers) missions withdata extracted from NASA’s CADRe. We then compare thismethod’s performance with the classical Rule of Thumb method.

Hooke, Melissa A↗
NASA NTRSDate not supplied

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗
NASA NTRS2023

A Sensitive Search for Supernova Emission Associated with the Extremely Energetic and Nearby GRB 221009A

We report observations of the optical counterpart of the long gamma-ray burst GRB 221009A. Due to the extreme rarity of being both nearby (z = 0.151) and highly energetic (E ϒ,iso ≥ 10 54 erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of γ-ray properties. Adopting a phenomenological power-law model for the afterglow and host galaxy estimates from high-resolution Hubble Space Telescope imaging, we use Bayesian model comparison techniques to determine the likelihood of an associated supernova (SN) contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence (K Bayes =10 1.2 for the presence of an additional component arising from an associated SN, SN 2022xiw, and find that it must be substantially fainter (<67% as bright at the 99% confidence interval) than SN 1998bw. Given the large and uncertain line-of-sight extinction, we attempt to constrain the SN parameters (M Ni , M ej and E KE ) under several different assumptions with respect to the host galaxy's extinction. We find properties that are broadly consistent with previous GRB-associated SNe: M Ni = 0.05–0.25 M ⊙ , M ej = 3.5–11.1 M ⊙ , and E KE = (1.6–5.2) × 10 52 erg. We note that these properties are weakly constrained due to the faintness of the SN with respect to the afterglow and host emission, but we do find a robust upper limit on M Ni of M Ni < 0.36 M ⊙ . Given the tremendous range in isotropic gamma-ray energy release exhibited by GRBs (seven orders of magnitude), the SN emission appears to be decoupled from the central engine in these systems.

Core-collapse supernovae↗
NASA NTRS2023

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗
NASA NTRSDate not supplied

Assessing Long-Term Health (LTH) Outcomes in Astronauts

NASA has linked potential Long-Term Health (LTH) risks of astronauts to their spaceflight experience, including the effects of space radiation, microgravity, and stressors such as isolation, sleep deprivation, and stress. Therefore, assessing the LTH outcomes for astronauts exposed to these hazards becomes a critical investigation in ensuring any deleterious effects are minimal. Prior research conducted into LTH outcomes of astronauts has not revealed increases in non-accidental mortality rates when compared to the U.S. general population, analog healthy cohorts, and professional athletes. One study did report an increased incidence of non-fatal cardiovascular disease events (Hazard Ratio=2.41, 95% Confidence Interval 1.26 to 4.63). Additional evidence indicates that astronauts may have increases in incidences of melanoma, prostate, and hematologic cancers, though the former is consistent with rates observed among aircraft pilots and the latter two are likely caused by detection bias due to increased screening on former astronauts.

long term health↗