Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust estimation”

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 883 records · Page 49

Using Gaia Data 2 to Constrain Local Dark Matter Density and Thin Dark Disk

We use stellar kinematics from the latest Gaia data release (DR2) to measure the local dark matter (DM) density ρDM in a heliocentric cylinder of radius R = 150 pc and half-height z = 200 pc. We also explore the prospect of using our analysis to estimate the DM density in local substructure by setting constraints on the surface density and scale height of a thin dark disk aligned with the baryonic disk and formed due to dissipative dark matter self-interactions. Performing the statistical analysis within a Bayesian framework for three types of tracers, we obtain ρDM = 0.016 ± 0.010 Mꙩ/p c 3 for A stars; early G stars give a similar result, while F stars yield a significantly higher value. For a thin dark disk, A stars set the strongest constraint: excluding surface densities (5–12) Mꙩ/pc 2 for scale heights below 100 pc with 95% confidence. The upper bound of this constraint implies . 1% of the Milky Way DM mass is present in a dissipative dark sector. Comparing our results with those derived using Tycho-Gaia Astrometric Solution (TGAS) data, we find that the uncertainty in our measurements of the local DM content is dominated by systematic errors that arise from assumptions of our dynamical analysis in the low z region. Furthermore, there will only be a marginal reduction in these uncertainties with more data in the Gaia era. We comment on the robustness of our method and discuss potential improvements for future work.

dark matter theory↗

USA Crop Yield Estimation with MODIS NDVI: Are Remotely Sensed Models Better Than Simple Trend Analyses?

Crop yield forecasting is performed monthly during the growing season by the United States Department of Agriculture’s National Agricultural Statistics Service. The underpinnings are long-established probability surveys reliant on farmers’ feedback in parallel with biophysical measurements. Over the last decade though, satellite imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) has been used to corroborate the survey information. This is facilitated through the Global Inventory Modeling and Mapping Studies/Global Agricultural Monitoring system, which provides open access to pertinent real-time normalized difference vegetation index (NDVI) data. Hence, two relatively straightforward MODIS-based modeling methods are employed operationally. The first model constitutes mid-season timing based on the maximum peak NDVI value, while the second is reflective of late-season timing by integrating accumulated NDVI over a threshold value. Corn model results nationally show the peak NDVI method provides a R^(2) of 0.88 and a coefficient of variation (CV) of 3.5%. The accumulated method, using an optimally derived 0.58 NDVI threshold, improves the performance to 0.93 and 2.7%, respectively. Both these models outperform simple trend analysis, which is 0.48 and 7.4%, correspondingly. For soybeans the R^(2) results of the peak NDVI model are 0.62, and 0.73 for the accumulated using a 0.56 threshold. CVs are 6.8% and 5.7%, respectively. Spring wheat’s R2performance with the accumulated NDVI model is 0.60 but just 0.40 with peak NDVI. The soybean and spring wheat models perform similarly to trend analysis. Winter wheat and upland cotton show poor model performance, regardless of method. Ultimately, corn yield forecasting derived from MODIS imagery is robust, and there are circumstances when forecasts for soybeans and spring wheat have merit too.

crop yield↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Training and Validation of Spectral Gap Filling Algorithm for Cpf-Ceres Intercalibration

The Climate Absolute Radiance and Refractivity Observatory (CLARREO) Pathfinder (CPF) mission is set to launch an SI-traceable reflective solar (RS) spectrometer aboard the International Space Station to measure Earth-reflected solar radiation with a radiometric uncertainty of 0.3% (k=1). The CPF intercalibration team has devised a cutting-edge methodology to accurately transfer the benchmark CPF calibration reference to the shortwave (SW) channel (200-5000 nm) of the Clouds and the Earth’s Radiant Energy System (CERES) instrument. The spectral range of CPF measurements spans from 350-2300 nm, while the CERES SW channel measures the Earth-reflected broadband solar radiances between 200 nm to 5 μm. To conduct precise CPF-CERES intercalibration analysis, the CPF-like spectral radiances outside the CPF spectral range need to be estimated to match the CERES SW spectral range. In response, the team has developed a fast algorithm that leverages spectrally redundant information within the CPF-measured portion through principal component analysis (PCA) and utilizes pre-established spectral correlation relationships among wavelengths to extend the CPF spectrum below 350 nm and above 2300 nm. Our results show that the algorithm achieves excellent accuracy in generating the missing energy in the UV and IR portions. The RMS error in the UV region is less than 4.5x10-3 W/m2/sr/nm, while in the IR region, it is smaller than 8x10-5 W/m2/sr/nm. Our methodology was validated using measured EMIT radiance data, which covers the spectral range from 0.381 μm to 2.493 μm. We employed EMIT radiances within the wavelength range of 0.43 – 2.25 μm to generate radiances for both the shorter wavelength range (0.381 – 0.43 μm) and longer wavelength range (2.25 – 2.493 μm). The generated radiances agree very well with the measured EMIT radiances. The standard deviation in the integrated broadband radiances was about 0.1%, and the bias is less than 0.004% for over 1.5 million EMIT measured samples. These statistics show that the spectral gap filling algorithm is robust and effective in substantially reducing the spectral difference-induced uncertainty in the CPF-CERES intercalibration samples.

Qiguang Yang↗

Estimating the Single-Trial Characteristics of Event-Related Responses: Evaluation of the MCERP Algorithm

Single-trial event-related responses collected during the course of an experiment are typically averaged before analysis resulting in a rather crude picture of event-related brain dynamics. It has been quite clear for some time that these responses exhibit trial-to-trial variability: however, the computational techniques necessary to deal with such responses in noisy conditions have not been available. To this end we have developed the multiple-component, event-related potential model (mcERP), which assumes that the each event-related response consists of a sum of multiple evoked components each described by a stereotypical waveshape. These waveshapes are allowed to vary in amplitude and onset latency from trial to trial, which allows us to capture, to first-order, the trial-dependent variations in event-related brain dynamics. We have constructed many sets of synthetic data designed to simulate intracortical recordings from a 15 channel, linear-array multielectrode implanted acutely in V1 of an awake-behaving macaque undergoing visual stimulation with a red light flash. This synthetic data was used to characterize the performance of the mcERP algorithm. First we quantified the degree to which such trial-to-trial variability aids in the identification of multiple components, and we demonstrate that amplitude variability is a more important factor in component separation than latency variability. Second, we quantified the behavior of the algorithm under two distinct signal-to-noise ratio (SNR) conditions: Gaussian noise independently present in each channel, and highly correlated (1/f distributed), far-field noise presented identically in each channel of the array. The mcERP algorithm was found to be robust to noise accurately identifying all component waveshapes and their associated single-trial characteristics down to SNR levels of -20dB for Gaussian noise and -7dB for 1/f far-field noise. Comparisons of the performance of this algorithm with factor analysis (FA) and independent component analysis (ICA) will be described by Knuth et al. (SFN abstracts, 2002). In addition, the advantages of application of mcERP to real data will be described by Shah et al, (these abstracts, 2002: SFN abstracts, 2002).

Knuth, K. H.↗

Retrieval of Snow Water Equivalent by the Precipitation Imaging Package (PIP) over Northern Great Lakes

Performance of the Precipitation Imaging Package (PIP) for estimating the snow water equivalent (SWE) is evaluated through a comparative study with the collocated National Oceanic and Atmospheric Administration National Weather Service snow stake field measurements. The PIP together with a vertically pointing radar, a weighing bucket gauge, and a laser-optical disdrometer was deployed at the NWS Marquette, Michigan office building for a long-term field study supported by the National Aeronautics and Space Administration’s Global Precipitation Measurement mission Ground Validation program. The site was also equipped with a weather station. During the 2017-18 winter, PIP functioned nearly uninterrupted at frigid temperatures accumulating 2345.8 mm of geometric snow depth over a total of 499 hours. This long record consists of 30 events, and the PIP-retrieved and snow stake field measured SWE differed less than 15% in every event. For the windy events, the wind corrected stake field record agreed better with the PIP-retrieved SWE. Two of the major events with the longest duration and the highest accumulation are examined in detail. The particle mass with a given diameter was much lower during a shallow, colder, uniform lake-effect event than in the deep, less cold, and variable synoptic event. This study demonstrated that PIP is a robust instrument for operational use, and is reliable for deriving the bulk properties of falling snow.

Ali Tokay↗

CHP-PRA Proof-of-Concept [Simulation] Sensitivity Assessment

An effort is underway to establish a Crew Health and Performance system (CHP) tradespace tool using a Probabilistic Risk Assessment (PRA) modeling and simulation system. The goal of the CHP-PRA effort is to provide a means of quantifying the integrated influence of CHP functions and capabilities on risk outcome metrics associated with health, performance, and long-term health. These metrics can then be used to establish potential risk-based trades on CHP system designed functionality and capabilities. Previously, our team demonstrated a proof-of-concept PRA approach that estimated the integrated influence of exercise countermeasures on 8 human system risks (Figure 1a) with outcomes associated with health and medical risk metrics. We reported that the change in the integrated relative health risk was small (Figure 1b) and that the small change in overall risk resulted from compounding and competing contribution levels of the individual risks. This interesting observation illustrates the emergent complexity of even straightforward representations of the human health and performance risk space and the ability of PRA models to capture this balance of global risk concerns. A key question that is not addressed in the initial analysis is “even though the global risk is relatively nominal, do any of the local risks become unacceptable?” In essence, we seek to determine what relative change in the human system risks are contributing to the relatively muted sensitivity of the proof-of-concept model combined risk assessments. Evaluations at the component risk level should elucidate if any individual risk reaches a high level that is subsequentially balanced by reductions in other areas. To further understand the relative changes in the component risks in the proof-of-concept model, and to elucidate how future refinements can be targeted, a means of establishing the contributions of the robustness of the proof-of-concept approach will be demonstrated.

astronaut health↗

Production and Characterization of Additively Manufactured Radiator Panels with Integral Branching Heat Pipes for High-Temperature Heat Rejection

Emerging concepts for fission surface power and nuclear electric propulsion necessitate lightweight, mechanically robust, and thermally efficient heat rejection radiators. State-of-the-art intermediate-temperature (~400 K) composite radiator assemblies have been developed based on titanium-water heat pipes bonded to metal, graphite, and carbon-fiber-based panels. NASA has identified a need for new radiator concepts that can operate at even higher temperatures (500 – 600 K), minimize thermal resistances and thermal stress failures at bond interfaces, and approach areal densities of 2 – 3 kg m -2 . To meet these needs, our team is developing additively manufactured (AM) radiator panels with integral branching wicking heat pipe networks. Water is selected as the working fluid for this temperature range. Based on simulations and thermal vacuum experiments, these branching embedded heat pipe networks can efficiently distribute heat over panels for finned surface efficiencies of η f >70% at T H = 500 K input heat. This paper first presents laser powder-bed fusion AM strategies to produce embedded porous structures for wicking heat pipes in Inconel 718 and titanium alloys (commercially pure and Ti-6Al-4V alloys). Post-build chemical and thermal treatments are described that yield hydrophilic wicking surfaces for operation with water. Transient rate-of-rise experiments with water and acetone are reported that yield estimates for AM wick porosity (ϵ), permeability (K), and effective pore radius (r pore ). Based on the wick characterization results, small prototype radiator panels (75 × 125 mm) with integrated heat pipe networks were manufactured. Heat rejection performance data are presented from cold thermal vacuum testing, with heat input temperatures up to ~510 K. Future efforts will focus on improving heat pipe performance, optimizing radiator mass, and evaluating larger panels to assess scalability.

thermal management↗

Production and Characterization of Additively Manufactured Radiator Panels With Integral Branching Heat Pipes for High-Temperature Heat Rejection

Emerging concepts for fission surface power and nuclear electric propulsion necessitate lightweight, mechanically robust, and thermally efficient heat rejection radiators. State-of-the-art intermediate-temperature (~400 K) composite radiator assemblies have been developed based on titanium-water heat pipes bonded to metal, graphite, and carbon-fiber-based panels. NASA has identified a need for new radiator concepts that can operate at even higher temperatures (500 – 600 K), minimize thermal resistances and thermal stress failures at bond interfaces, and approach areal densities of 2 – 3 kg m -2 . To meet these needs, our team is developing additively manufactured (AM) radiator panels with integral branching wicking heat pipe networks. Water is selected as the working fluid for this temperature range. Based on simulations and thermal vacuum experiments, these branching embedded heat pipe networks can efficiently distribute heat over panels for finned surface efficiencies of η f >70% at T H = 500 K input heat. This paper first presents laser powder-bed fusion AM strategies to produce embedded porous structures for wicking heat pipes in Inconel 718 and titanium alloys (commercially pure and Ti-6Al-4V alloys). Post-build chemical and thermal treatments are described that yield hydrophilic wicking surfaces for operation with water. Transient rate-of-rise experiments with water and acetone are reported that yield estimates for AM wick porosity (ϵ), permeability (K), and effective pore radius ( r pore ). Based on the wick characterization results, small prototype radiator panels (75 × 125 mm) with integrated heat pipe networks were manufactured. Heat rejection performance data are presented from cold thermal vacuum testing, with heat input temperatures up to ~510 K. Future efforts will focus on improving heat pipe performance, optimizing radiator mass, and evaluating larger panels to assess scalability.

nuclear electric propulsion↗

Real-World Testing of LiDAR-Inertial Based Navigation and Mapping for Precision Landing

The fusion of LiDAR and inertial measurements duringspacecraft descent and landing can be used to estimate alander’s navigation state and map the terrain below. Together,these data products can be used to enable safe and preciselanding on celestial bodies for which a priori orbital reconnaissanceis insufficient for hazard detection and avoidance.Unlike camera images used in visual terrain relative navigation,LiDAR scans are insensitive to changes in illumination; as aresult, the technique can be used to land in poorly lit areas,or at times of day when the lighting conditions are incongruentwith existing orbital imagery. In this paper, we extend previouswork in which we introduced a factor graph based smoothingapproach for LiDAR-inertial navigation and mapping. Whereasthe algorithms were previously tested on simulated data, thispaper presents testing on real-world data. Data from theAutonomous Landing Hazard Avoidance Technology (ALHAT)airplane flight tests in the Yucca Flats and Death Valley in2009 (FT3), the Morpheus vertical take off and landing flighttests at Kennedy Space Center in 2014 (FT6), and the landingof Perseverance and Ingenuity on Mars in 2021 (M2020) wereused to evaluate algorithm performance. In this paper, weextend our LiDAR-inertial technique to work with a variety ofranging technologies: single point laser altimetry (FT3), denseflash LiDAR (FT6), and six-beam radar (M2020). A thoroughperformance analysis for all three datasets is presented. Datasetpreparation, improvements in algorithm robustness, and outlierrejection, which were necessitated by the transition to realworlddata, are discussed.

Trawny, Nikolas↗

The Repair Maintenance and Fabrication Facility in the Common Habitat Architecture

The Common Habitat Architecture seeks to increase the habitability of long-duration human spaceflight systems. A key aspect of this is vehicle survivability. Missions beyond low Earth orbit need onboard capabilities for Repair, Maintenance, And Fabrication (RMAF) to overcome potential contingency scenarios. Strategies employed in historic human spaceflight such as redundancy management, reliability, sparing, orbital replacement units, and aborts may be insufficient by themselves. Based on subject matter input, a list of 53 critical failures defining a set of incidents that can render a key spacecraft subsystem inoperable were generated. A subsequent analysis found that a robust in-space RMAF system capable of performing 14 key functions can potentially repair a subsystem plagued by any of these failures. An ancillary benefit is this capability may provide psychological benefits to the crew, by enabling greater self-sufficiency in earth-independent problem solving. A basic RMAF facility has been defined for the Common Habitat, consisting of five workstations. A work bench and computer workstation provide a multipurpose horizontal work surface, computing interface, and tools storage. A CNC machining center provides a subtractive manufacturing capability for metals and plastics. A multi-material 3D printing facility provides additive manufacturing capabilities for plastic, metals, and printed electronics. A welding facility is used for joining metal components where a higher strength is needed than can be achieved with fasteners or adhesives. A glovebox facility is used to perform work that is too hazardous for any of the other workstations. This may include hardware brought in from outside the spacecraft that could potentially contaminate the cabin environment. Forward work includes considering the accommodation of additional manufacturing processes not modeled in the current system, assessing the ability of systems to operate in partial gravity and microgravity environments, incorporation of the system into the Common Habitat Computer Aided Design (CAD) model, bottoms-up mass estimating, and a crew time analysis.

Habitat↗

Thermal Infrared Remote Sensing for Analysis of Landscape Ecological Processes: Methods and Applications

Thermal Infrared (TIR) remote sensing data can provide important measurements of surface energy fluxes and temperatures, which are integral to understanding landscape processes and responses. One example of this is the successful application of TIR remote sensing data to estimate evapotranspiration and soil moisture, where results from a number of studies suggest that satellite-based measurements from TIR remote sensing data can lead to more accurate regional-scale estimates of daily evapotranspiration. With further refinement in analytical techniques and models, the use of TIR data from airborne and satellite sensors could be very useful for parameterizing surface moisture conditions and developing better simulations of landscape energy exchange over a variety of conditions and space and time scales. Thus, TIR remote sensing data can significantly contribute to the observation, measurement, and analysis of energy balance characteristics (i.e., the fluxes and redistribution of thermal energy within and across the land surface) as an implicit and important aspect of landscape dynamics and landscape functioning. The application of TIR remote sensing data in landscape ecological studies has been limited, however, for several fundamental reasons that relate primarily to the perceived difficulty in use and availability of these data by the landscape ecology community, and from the fragmentation of references on TIR remote sensing throughout the scientific literature. It is our purpose here to provide evidence from work that has employed TIR remote sensing for analysis of landscape characteristics to illustrate how these data can provide important data for the improved measurement of landscape energy response and energy flux relationships. We examine the direct or indirect use of TIR remote sensing data to analyze landscape biophysical characteristics, thereby offering some insight on how these data can be used more robustly to further the understanding and modeling of landscape ecological processes.

Quattrochi, Dale A.↗

Robust Exploration and Commercial Missions to the Moon Using NTR LANTR Propulsion and Lunar-Derived Propellants

The nuclear thermal rocket (NTR) has frequently been identified as a key space asset required for the human exploration of Mars. This proven technology can also provide the affordable access through cislunar space necessary for commercial development and sustained human presence on the Moon. In his post-Apollo Integrated Space Program Plan (1970-1990), Wernher von Braun, proposed a reusable nuclear thermal propulsion stage (NTPS) to deliver cargo and crew to the Moon to establish a lunar base before undertaking human missions to Mars. The NTR option was selected by von Braun because it was a demonstrated technology capable of generating both high thrust and high specific impulse (Isp 900 s) twice that of todays best chemical rockets. In NASAs Mars Design Reference Architecture (DRA) 5.0 study, the crewed Mars transfer vehicle used three 25 klbf Pewee engines the smallest and highest performing engine tested in the Rover program along with graphite composite fuel. Smaller, lunar transfer vehicles consisting of a NTPS using three approximately 16.5 klbf Small Nuclear Rocket Engines (SNREs), an in-line propellant tank, plus the payload can enable a variety of reusable lunar missions. These include cargo delivery and crewed lunar landing missions. Even weeklong tourism missions carrying passengers into lunar orbit for a day of sightseeing and picture taking are possible. The NTR can play an important role in the next phase of lunar exploration and development by providing an affordable in-space lunar transportation system (LTS) that can allow initial outposts to evolve into settlements supported by a variety of commercial activities such as in-situ propellant production used to supply strategically located propellant depots and transportation nodes. The utilization of iron-rich volcanic glass or lunar polar ice (LPI) deposits (each estimated at billions of metric tons) for propellant production can significantly reduce the launch mass requirements from Earth and can enable reusable, surface-based lunar landing vehicles (LLVs) using liquid oxygen/hydrogen (LOX/LH2) chemical rocket engines. Afterwards, LOX/LH2 propellant depots can be established in lunar equatorial and polar orbits to supply the LTS. At this point a modified version of the conventional NTR called the LOX-augmented NTR, or LANTR would be introduced into the LTS allowing bipropellant operation and leveraging the mission benefits of refueling with lunar-derived propellants for Earth return. The bipropellant LANTR engine utilizes the large divergent section of its nozzle as an afterburner into which oxygen is injected and supersonically combusted with nuclear preheated hydrogen emerging from the engines choked sonic throat essentially scramjet propulsion in reverse. By varying the oxygen-to-hydrogen mixture ratio, LANTR engines can operate over a range of thrust and Isp values while the reactor core power level remains relatively constant. Eventually, a LANTR-based LTS can enable a rapid commuter shuttle with one-way trip times to and from the Moon ranging from 36 to 24 hours. Even if only 1 of the extracted propellant from identified volcanic glass and polar ice deposits were available for use in lunar orbit, such a supply could support daily commuter flights to the Moon for many thousands of years! An evolutionary mission architecture is outlined and a variety of lunar missions and transfer vehicle designs are examined, along with the increasing demands on propellant production as mission complexity increases. A comparison of vehicle features and engine operating characteristics, for both NTR and LANTR engines, is also provided along with a brief discussion on the propellant production issues associated with using volcanic glass and LPI as source material.

Spacecraft design↗

Dense Feature Tracking of Atmospheric Winds with Deep Optical Flow

Atmospheric winds are a key physical phenomenon impacting natural hazards, energy transport, ocean currents, large-scale circulation, and ecosystem fluxes. Observing winds is a complex process and presents a large gap in NASA’s Earth Observation System. Atmospheric motion vectors (AMVs) aim to fill this gap by making numerical estimates of cloud movement between sequences of multi-spectral satellite images, tracking clouds and water vapor. Recent imaging hardware and software advancements have enabled the use of numerical optical flow techniques to produce accurate and dense vector fields outperforming traditional methods. This work presents WindFlow as the first machine learning based system for feature tracking atmospheric motion using optical flow. Due to the lack of large-scale satellite-based observations, we leverage high-resolution numerical simulations from NASA's GEOS-5 Nature Run to perform supervised learning and transfer to satellite images. We demonstrate that our approach using deep learning based optical flow scales to ultra-high-resolution images of size 2881x5760 with less than 1 m/s bias and 2.5 m/s average error. Four network and learning architectures are compared and it is found that recurrent all-pairs field transforms (RAFT) produces the lowest errors on all metrics for wind speed and direction. Results on held out numerical outputs shows RAFT's good performance in each of the spatial, temporal, and physical dimensions. A comparison between WindFlow and an operational AMV product against rawinsonde observations show that RAFT transfers across simulations and thermal infrared satellite observations. This work shows that machine learning based optical flow is an efficient approach to generating robust feature tracking for AMVs consistently over large regions.

Atmospheric winds↗

Development of a Supervisory Tool for Fault Detection and Diagnosis of DC Electric Power Systems with the Application of Deep Space Vehicles

This dissertation formulates the problem of fault detection and diagnosis of DC electric power systems for the application of autonomous spacecraft. The ability to accurately identify and isolate failures in the electrical power system is critical to ensure the reliability of a spacecraft. This problem becomes more pronounced during deep space missions that lack the ability to monitor from ground control. The current state of electrical power system fault supervision is insufficient to guarantee highly reliable and robust operation. To solve this issue, a combination of model-based and rules-based techniques are used in a hierarchical framework to improve the diagnostic performance of the spacecraft electrical power system. Noise, disturbances, and modeling errors are considered in the design of the method. Practical considerations related to the hardware and software are discussed for the flight application. A wide array of failure types are simulated in a series of experiments to assess the functionality of the design. The experiments showed that the methods used improved the diagnostic capability of the autonomous system while taking into account the limitations attributed to flight software requirements. The significance of this study is to provide a framework capable of advanced diagnostics of an electrical power system with little to no interaction from a human operator.

Fault Detection and Diagnosis↗

20-Years of Atmospheric Temperature, Water Vapor, Cloud, and Surface Temperature Anomalies and Trends Derived From Operational Hyperspectral Ir Sounders

Hyperspectral IR sounders such as AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C provide high-quality atmospheric temperature, water, vapor, and greenhouse gas vertical profiles. Additionally, they provide atmospheric cloud properties, surface emissivity, and surface skin temperatures. We have developed two algorithms which can consistently derive these products from multiple IR sounders. The first one is a Single Field-of-view Sounder Atmospheric Product (SIFSAP) algorithm and the second one is a Climate Fingerprinting Sounder Product (ClimFiSP) algorithm. Compared to current operational AIRS and CrIS Level-2 (L2) algorithms, which perform one retrieval for each 3 by 3 field of views (FOVs) using a cloud-clearing approach, the SiFSAP algorithm, on the other hand, performs one retrieval for each FOV using an all-sky optimal estimation approach. The SiFSAP algorithm retrieves all the above-mentioned atmosphere and surface properties simultaneously including cloud properties with 3-time higher spatial resolution and 9-times more products. The core of the SiFSAP algorithm is an accurate and fast Principal Component-based Radiative Transfer Model (PCRTM), which can calculate hyperspectral radiance spectra under both clear and cloudy conditions. The PCRTM was developed in the past decade using consistent reference line-by-line radiative transfer model and spectroscopy for hyperspectral sounders such as AIRS, CrIS, IASI, NAST-I, and S-HIS. The SiFSAP retrieval algorithm also uses the same climatology a priori and associated covariances, which makes it ideal for generating high quality products for both weather and climate applications. Climate products are typically derived by performing spatial and temporal averaging of L2 products. It is a time-consuming process to generate L2 data products since AIRS, CrIS, and IASI have millions of observations each day with thousands of spectral channels for each observation. Additionally, differences in L2 retrieval algorithms for different satellite sensors can lead to errors in the climate products. Our ClimFiSP algorithm, which performs retrievals from spatiotemporally averaged L1 hyperspectral radiances directly, will be orders of magnitude faster than traditional method. he ClimFiSP algorithm uses consistent radiative kernels and a robust spectral fingerprinting method. It provides accurate data climate data fusion products from multiple satellite sensors. We have applied this method to both AIRS and CrIS (on SNPP and on NOAA 20) data and generated two decades climate data records for atmospheric temperature, water vapor, cloud, trace gases, and surface skin temperature. Both SiFSAP and ClimFiSP will be available at NASA GES DISC data center for public access.

Xu Liu↗

Constraining the Climate and Ocean PH of the Early Earth with a Geological Carbon Cycle Model

The early Earth's environment is controversial. Climatic estimates range from hot to glacial, and inferred marine pH spans strongly alkaline to acidic. Better understanding of early climate and ocean chemistry would improve our knowledge of the origin of life and its coevolution with the environment. Here, we use a geological carbon cycle model with ocean chemistry to calculate self-consistent histories of climate and ocean pH. Our carbon cycle model includes an empirically justified temperature and pH dependence of seafloor weathering, allowing the relative importance of continental and seafloor weathering to be evaluated. We find that the Archean climate was likely temperate (0-50deg C) due to the combined negative feedbacks of continental and seafloor weathering. Ocean pH evolves monotonically from 6.6(sup + 0.6) (sub -0.4) (2 sigma) at 4.0 Ga to 7.0(sup +0.7) (sub-0.5) (2 sigma) at the Archean-Proterozoic boundary, and to 7.9(sup +0.1) (sub -0.2) (2 sigma) at the Proterozoic-Phanerozoic boundary. This evolution is driven by the secular decline of pCO2, which in turn is a consequence of increasing solar luminosity, but is moderated by carbonate alkalinity delivered from continental and seafloor weathering. Archean seafloor weathering may have been a comparable carbon sink to continental weathering, but is less dominant than previously assumed, and would not have induced global glaciation. We show how these conclusions are robust to a wide range of scenarios for continental growth, internal heat flow evolution and outgassing history, greenhouse gas abundances, and changes in the biotic enhancement of weathering.

Krissansen-Totton, Joshua↗