Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Employment”

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 649 records · Page 36

Assessing the Calibration Differences in the Reflective Solar Bands of Terra MODIS and Landsat-7 Enhanced Thematic Mapper Plus

Long-term data records obtained from Earth observing sensors depend not only onthe calibration accuracy of individual sensors but also on the consistency across instruments andplatforms. Hence, sensor calibration intercomparison plays a vital role for a better understandingof various science products. The Moderate Resolution Imaging Spectroradiometer (MODIS)and enhanced thematic mapper plus (ETM+) on the Terra and Landsat 7 platforms have operatedsuccessfully since their launch, collecting measurements in the reflective solar and infrared partsof the spectrum. Terra MODIS has employed a reflectance-based calibration since beginning itsmission. In the case of ETM+, a radiance-based calibration was employed until recent years,when a reflectance-based calibration was introduced. Being in the AM constellation with lessthan 30 min difference in overpass times, near-simultaneous Earth scene measurements can beeffectively used to assess the calibration differences between the spectrally matching bands ofthese two instruments. The pseudoinvariant calibration sites (PICS) in the North African desertare widely used for on-orbit calibration and validation of satellite sensors. Four PICS from thisregion have been employed to assess the multitemporal reflectance differences. Correction forbidirectional reflectance, spectral response function mismatch, and impacts of atmosphericwater-vapor have been incorporated to provide an assessment of the long-term stability ofeach spectral band and reflectance differences amongst them. Results indicate that the spectralbands of both instruments show a long-term stability to within 2% from 2000 to 2017. Thetop-of-atmosphere reflectances between the two instruments postcorrection agree to within 4%.Also included in this paper is a detailed discussion of various parameters contributing to theuncertainties of this cross-calibration. The techniques presented in this paper can be furtherextended to perform similar intercomparison between Landsat 8 Operational Land Imager, AquaMODIS, and Suomi-NPP VIIRS.

Angal, Amit↗

TPSAS-NF1676L-19788-DND

The calibration of the historical AVHRR visible channels has always been hindered by the degrading NOAA satellite orbits. The complete AVHRR record spans over three decades making it useful for cloud, aerosol, and land use climate studies. Studies that monitor long-term changes in these properties require climate quality calibration. The AVHRR sensors do not have any onboard visible calibration, unlike the IR channels. Many AVHRR post-launch calibration methods have been published. Some of these are based on the pseudo invariant targets such as desert and polar ice and others employ direct calibration transfer from a well-calibrated sensor such as MODIS. Most historical studies only employed only one method and reconciling calibration differences from multiple studies was difficult. This study employs multiple calibration approaches, which are then merged according to their individual uncertainty. Since the visible spectral response functions are similar across AVHRR sensor and that the morning and afternoon orbits usually degrade in the same manner, an afternoon and morning orbiter during the MODIS era is chosen as the reference AVHRR sensor. The reference AVHRR sensor is inter-calibrated with the Aqua-MODIS Collection 6 standard using simultaneous nadir overpass (SNO) radiance pairs. Three pseudo invariant target approaches are utilized, deserts, polar ice and deep convective clouds. These sites are characterized using reference AVHRR reflectances as a function of solar zenith angle. The multiple calibration approach is then validated with non-reference AVHRR sensor during the MODIS era, using AVHRR and MODIS SNO measurements. For historical AVHRR satellite having consistent calibration across invariant targets then validates the approach. This calibration effort is in sponsored by the NOAA CDR program to provide climate quality visible calibration coefficients to the remote sensing community. Each of the calibration approaches are unique and when utilized in tandem can accurately calibrate the AVHRR record. AVHRR calibration methods cannot rely on repeatable annual orbits for calibration model development. Inconsistent results prompt code and invariant target examinations that cannot be determined using only single method approaches. Results from this study will be presented at the meeting.

David Doelling↗

Multiphysics Computational Analysis of a Perforated Plate Cooling Flow

A multiphysics analysis of a film cooling flow was conducted using a computational fluid dynamics (CFD) solver in conjunction with a conjugate heat transfer (CHT) solver. The film cooling problem consisted of a square nozzle producing a heated subsonic jet over a perforated plate with three sections of 45 cooling holes each. Results from the multiphysics analysis were compared to a CFD-only approach and experimental results. Reynolds-averaged Navier-Stokes (RANS) computations were employed first, with and without the CHT solver. The multiphysics approach provided much closer agreement with experimentally measured plate temperatures than the CFD-only approach. Comparisons of computed velocities and temperatures were also made with experimental data. Hybrid RANS-large eddy simulation calculations employing the Delayed Detached Eddy Simulation (DDES) approach were also conducted within the multiphysics framework. The DDES approach was unable to resolve turbulent structures in the near wall regions. Hence, the RANS part of the DDES model was used in near wall regions, and produced surface temperatures along with velocity and temperature profiles close to the wall that were very similar to the RANS results. The DDES approach did resolve structures in the outer portion of the boundary layer especially near the end of the plate where the boundary layer was thicker. The key result was the improved surface temperatures produced by the multiphysics approach employing CHT over the fluid-only approach.

Computational Fluid Dynamics↗

Surface-Normal Active Flow Control for Lift Enhancement and Separation Mitigation for High-Lift Common Research Model

This study explores the feasibility of using small surface-normal jets near the flap trailing edge as active aerodynamic load control on the high-lift Common Research Model wing/-body configuration. Chimera Grid Tools are utilized to generate structured curvilinear overset grids, and the Reynolds-averaged Navier-Stokes solver OVERFLOW is employed to solve for the flow-field around the geometry. The so-called microjet is initially employed across both the inboard and the outboard flaps on the pressure-side near the trailing edge. It is shown that implementing the microjet on the inboard flap is more effective compared to implementing it on the outboard flap. This is because, prior to microjet implementation, the flow on the outboard flap exhibits extensive separation, while the flow on the inboard flap exhibits moderate separation. For microjet implementation across the inboard flap only, the relationship between momentum coefficient of the microjet and lift-enhancement is found to be∆CL'1.66√Cμfor the rangeCμ= 0.00−0.012. We show that implementing a microjet with a jet velocity ratio of one, which corresponds toCμ= 0.003, can shift the linear region of the lift curve by ∆CL= 0.08. The linear shift in the lift curve is significant for enhancing airplane performance such as increasing its payload. Microjet implementation effects on the drag coefficient are also investigated through a drag decomposition analysis. Further, we employ an induced drag analysis based on the spanwise load distribution and show that the microjet-related increase in pressure drag coefficient is dominated by the increase in the induced drag while microjet implementation reduces the form drag. These preliminary results show that favorable changes in aerodynamic performance can be achieved by using the surface-normal jets presented in this study.

ARMD↗

MER Opportunity dust-storm recovery operations and implications for future Mars surface missions

In June 2018, NASA’s Mars Exploration Rover Opportunity became engulfed in the most intense global dust storm observed in its 14-year mission and in Mars’ recorded history. Sapped of life-giving solar energy, Opportunity fell silent to ground operators on June 11, in what would be its final call home. Over the course of the next eight months, the MER team employed numerous recovery efforts and radiated over one thousand commands to wake the silent rover. Although unsuccessful, MER’s dust-storm recovery team changed the paradigm of dust-storm operations from Opportunity’s survival of a previous global dust storm in 2007. In this paper, the authors offer a glimpse into MER’s recovery efforts and lessons learned from Opportunity dust-storm operations for future solar-powered Mars surface missions. In the first section, the authors discuss the indicators of the approaching dust storm and the actions the team employed to reduce Opportunity’s power consumption and preserve available battery charge prior to loss of contact. The second section discusses the team’s recovery efforts until Opportunity’s declared End-of-Mission, the steps the team took to re-establish contact with Opportunity, the assumptions made during each step of the recovery process, and the commanding actions employed. Finally, the authors discuss the operational impacts of global dust storms on the safe and successful operation of solar-powered spacecraft on the Martian surface, and offer design recommendations for future solar-powered missions from the lessons learned during Opportunity’s 14-year mission and experience through two global dust storms.

Nelson, Robert W.↗

Improving Computational Efficiency of Prognostics Algorithms in Resource-Constrained Settings

In engineering and aerospace applications, it is vital to operational success to have insight into the expected performance and health of physical systems. The field of prognostics and health management provides quantitative methods for monitoring, predicting, and managing system health. Prognostics algorithms can be employed to assess the current state of a system, propagate the system throughout time, and predict potential anomalies or failures that may occur. While they can provide accurate prediction results, effective prognostics algorithms can be challenging to use in resource-constrained settings due to computational limitations and high computational latency, leading to obsolete predictions. Thus, computationally efficient and accurate algorithms are necessary for future remaining useful life predictions. In this work, we implement new algorithmic approaches for prediction, quantitatively compare them via a battery degradation use-case, and provide recommendations of potential improvements to a prognostics framework. One approach to prediction is through sampling, whereby the current state of a physical system is sampled many times and each sample is propagated forward until failure is reached, resulting in a distribution of failure values. To improve the efficiency of this process, we implemented five new algorithmic approaches to prediction, including three distinct sampling methods (standard Monte Carlo, Quasi-Monte Carlo, and Latin Hypercube Sampling), a variable time step algorithm, and a variable sample size algorithm. To compare the algorithms, we employ a variety of metrics designed specifically to analyze both computational efficiency and model accuracy. Our metrics include accuracy to compare the average predicted value to ground truth, mean absolute deviation to illustrate dispersion, specific percentile error to describe accuracy within a user-defined risk tolerance, and code run-time. To quantitatively analyze our results, we employ a use-case of degradation of a Lithium-ion battery. We use an electrochemistry-based model to describe the current health state of the battery, and implement our prediction algorithms to propagate forward in time until end-of-discharge (EOD) is reached. Notably, through this work it was found that none of our sampling approaches had a significant impact on computational efficiency or model accuracy in predicting EOD of the battery. We find that while the sampling methods are unique, the distributions they generate are similar, ultimately producing final predictions that are nearly identical. In exploring the effect of the time step within the prediction algorithm, we found that prediction accuracy was highly dependent on the time step used, and that implementing a variable time step within a particular prediction may provide an increase in computational efficiency while also maintaining prediction accuracy. Finally, implementing a variable sample size also affected prediction, and our results show that tuning both the magnitude and timing of the sample size adjustment can result in improved computation speed and maintained prediction accuracy. Taken together, our findings highlight the challenge of performing prognostics in resource-constrained settings, and illustrate the potential of developing new prediction algorithms to improve computational efficiency.

prognostics↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

Extended Stroke and Miniaturized Reverse-Operation DTE Thermal Switches

This paper describes two advanced thermal switches under development that are offshoots of the recently flight-qualified 2500:1 reverse-operation differential thermal expansion (DTE) thermal switch (ROD-TSW). The first is the extended stroke ROD-TSW (ES-ROD-TSW) while the second is the miniaturized ROD-TSW (mini-ROD-TSW). Each employs multiple diametrically nested DTE stages in contrast to the single-stage ROD-TSW. The ES-ROD-TSW employs two stages and a negative CTE material known as Allvar and its main purpose is to extend its stroke for use in non-vacuum planetary environments. The mini-ROD-TSW employs four stages to miniaturize its size for use in compact lunar instruments and rovers. The ES-ROD-TSW was TVAC tested at JPL and its ON/OFF ratio is 13000:1 based on a measured ON conductance of 6.89 W/K and a calculated OFF conductance of 0.00053 W/K. The mini-ROD-TSW is targeted to have an ON conductance of 1.5 W/K and an OFF conductance of 0.0015 W/K, both verified by analysis. Test verification of the mini-ROD-TSW remains to be carried out. The paper will provide design/performance details of each device.

Farmer, Jeffery T.↗

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

Global 3D Data Visualization and Analysis Platform With Advanced Machine Learning Capabilities in Support of Lunar Exploration

Introduction: The science goals for NASA’s Artemis program include: a) Understanding the character and origin of lunar polar volatiles, b) Conducting experimental science in the lunar environment and c) Investigating and mitigating exploration risks [1]. The permanently shadowed regions (PSRs) on the Lunar south pole are expected to host large quantities of water-ice and volatiles that are important for sustainable Lunar exploration [2]. There are several missions such as onboard Korea Pathfinder Lunar Orbiter (KPLO: Korean name Danuri) with onboard ShadowCam camera [3], Astrobotic Peregrine Mission One [4], and other efforts underway to obtain high resolution topographic, minerals, volatiles and other information on the moon. We envision a need in immediate future for platforms to integrate these data sets, provide rendering and visualization capabilities in the context of a 3D Lunar globe for easier information access and analysis. NASA's Celestial Mapping System (CMS) [5] is developed to address the need for 3D tools for planetary science investigations, mission planning, in-situ operations, in a 3D-first design constructed around a unified view of a planetary globe. At present CMS provides many critical functionalities that include: 1) equipment planning and optimized placement on Lunar surface 2) line-of-sight (LOS) analysis 3) powerful measurement tools based on 3D terrain with realistic 3D models to represent rovers, astronauts and equipment 4) visualization of derived mapping products (e.g. resource maps), and 5) a data engine for hosting new observations that are not available in other contemporary lunar data tools [5, 7]. Planetary Data Ingestion: CMS can consume and analyze data from locally hosted and external third party sources. It is compatible with Open Geospatial Consortium (OGC) data and file standards and currently integrates datasets from the Astrogeology Science Center of USGS. This includes global and local data acquired from NASA (LRO, Clementine, Lunar Orbiter) and JAXA (SELENE/Kaguya), with capability of integrating more datasets. In addition, users can specify other WMS-hosted data endpoints, which CMS can then query and stream data from automatically. To set-up an automated process for ingestion and accurate rendering, visualization and analysis of external 3rd party planetary datasets within CMS, we initiated the process of ingesting unique dataset of super-enhanced images of the permanently shadowed regions (PSRs) at the lunar poles which were produced by the Hyper-effective nOise Removal U-net Software (HORUS) tool [8]. This tool was developed to enhance the extremely low-light images of the interior of PSRs and provide the ability to see within these regions at and discern surface features (i.e. boulders and craters) down to 3 meters in size. We focused on the Nobile region on the Lunar south pole, selected site for VIPER mission and stitched several images to create a high-resolution map within one of the PSR of Nobile crater. Figure 1 shows the dark PSR zone form the original NAC layer of LRO as the base layer (left image) and the illuminated areas within that crater (center) which was created by ingesting and merging several of HORUS generated images. At present we employ a semi-automated process to ensure spatial accuracy and merger of several overlapping zones. However, we are in the process of completely automating this process by employing AI based techniques that would rank, sort, and stack the images based on their information density. The georectification of the images would employ selected features. Analysis on Ingested Planetary Datasets: Once an external planetary data-set is successfully ingested, georectified and merged seamlessly as a data-layer; CMS’ numerous analysis tools can be used on this data. A Line of Sight (LOS) tool has been developed for CMS which analyzes terrain profiles and obstructions to determine visibility for remote observers [5,6]. Figure 1 (right image) shows the viewshed analysis on the same PSR in the Nobile region. The yellow pin shows the observer location outside the PSR. The yellow area shows the visible part of PSR. The obstructed area with no visibility for the observer is shown in red. The Measurements tool allows the user to take area and distance measurements of features on the terrain using various shapes. Measurement type can be specified in a number of ways: Line, Path, Polygon, Circle, Ellipse, Square, Rectangle or Freehand. Once the shape is specified, elevation information can then be extracted along each of these shapes. Figure 2 (left) shows the measurements performed on a crater n illuminated PSR in Nobile region. The equipment placement tool allows the user to place a 3D equipment model at a desired location and analyze its coverage area. The equipment placement tool is coupled with LOS to determine the coverage. Figure 2 (right) shows an equipment placed on the Lunar terrain and it’s coverage area. The red rays are blocked sight lines and the green rays are non-obstructed sight lines with the cyan lines showing the point of intersection with the terrain. More details are provided in the video demonstrations in Reference 5. Overcoming Polar Distortions: 3D geospatial applications exhibit significant distortions in polar imagery due to several reasons: 1) distortions in the source imagery, 2) incompatible tessellation algorithms at the poles, and 3) map projections. We are leveraging new tessellation algorithms and reprojecting data using projections that are better suited for Lunar poles. The goal is to seamlessly switch to polar projections while maintaining 3D view and navigation.

Maps↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

Global 3D Data Visualization and Analysis Platform with Advanced Machine Learning Capabilities in Support of Lunar Exploration

Introduction: The science goals for NASA’s Artemis program include: a) Understanding the character and origin of lunar polar volatiles, b) Conducting experimental science in the lunar environment and c) Investigating and mitigating exploration risks. The permanently shadowed regions (PSRs) on the Lunar south pole are expected to host large quantities of water-ice and volatiles that are important for sustainable Lunar exploration. There are several missions such as onboard Korea Pathfinder Lunar Orbiter (KPLO: Korean name Danuri) with onboard ShadowCam camera, Astrobotic Peregrine Mission One [4], and other efforts underway to obtain high resolution topographic, minerals, volatiles and other information on the moon. We envision a need in immediate future for platforms to integrate these data sets, provide rendering and visualization capabilities in the context of a 3D Lunar globe for easier information access and analysis. NASA's Celestial Mapping System (CMS) is developed to address the need for 3D tools for planetary science investigations, mission planning, in-situ operations, in a 3D-first design constructed around a unified view of a planetary globe. At present CMS provides many critical functionalities that include: 1) equipment planning and optimized placement on Lunar surface 2) line-of-sight (LOS) analysis 3) powerful measurement tools based on 3D terrain with realistic 3D models to represent rovers, astronauts and equipment 4) visualization of derived mapping products (e.g. resource maps), and 5) a data engine for hosting new observations that are not available in other contemporary lunar data tools. Planetary Data Ingestion: CMS can consume and analyze data from locally hosted and external third party sources. It is compatible with Open Geospatial Consortium (OGC) data and file standards and currently integrates datasets from the Astrogeology Science Center of USGS. This includes global and local data acquired from NASA (LRO, Clementine, Lunar Orbiter) and JAXA (SELENE/Kaguya), with capability of integrating more datasets. In addition, users can specify other WMS-hosted data endpoints, which CMS can then query and stream data from automatically. To set-up an automated process for ingestion and accurate rendering, visualization and analysis of external 3rd party planetary datasets within CMS, we initiated the process of ingesting unique dataset of super-enhanced images of the permanently shadowed regions (PSRs) at the lunar poles which were produced by the Hyper-effective nOise Removal U-net Software (HORUS) tool. This tool was developed to enhance the extremely low-light images of the interior of PSRs and provide the ability to see within these regions at and discern surface features (i.e. boulders and craters) down to 3 meters in size. We focused on the Nobile region on the Lunar south pole, selected site for VIPER mission and stitched several images to create a high-resolution map within one of the PSR of Nobile crater. Figure 1 shows the dark PSR zone form the original NAC layer of LRO as the base layer (left image) and the illuminated areas within that crater (center) which was created by ingesting and merging several of HORUS generated images. At present we employ a semi-automated process to ensure spatial accuracy and merger of several overlapping zones. However, we are in the process of completely automating this process by employing AI based techniques that would rank, sort, and stack the images based on their information density. The georectification of the images would employ selected features. Analysis on Ingested Planetary Datasets: Once an external planetary data-set is successfully ingested, georectified and merged seamlessly as a data-layer; CMS’ numerous analysis tools can be used on this data. A Line of Sight (LOS) tool has been developed for CMS which analyzes terrain profiles and obstructions to determine visibility for remote observers. Figure 1 (right image) shows the viewshed analysis on the same PSR in the Nobile region. The yellow pin shows the observer location outside the PSR. The yellow area shows the visible part of PSR. The obstructed area with no visibility for the observer is shown in red. The Measurements tool allows the user to take area and distance measurements of features on the terrain using various shapes. Measurement type can be specified in a number of ways: Line, Path, Polygon, Circle, Ellipse, Square, Rectangle or Freehand. Once the shape is specified, elevation information can then be extracted along each of these shapes. Figure 2 (left) shows the measurements performed on a crater n illuminated PSR in Nobile region. The equipment placement tool allows the user to place a 3D equipment model at a desired location and analyze its coverage area. The equipment placement tool is coupled with LOS to determine the coverage. Figure 2 (right) shows an equipment placed on the Lunar terrain and it’s coverage area. The red rays are blocked sight lines and the green rays are non-obstructed sight lines with the cyan lines showing the point of intersection with the terrain. More details are provided in the video demonstrations in Reference 5. Overcoming Polar Distortions: 3D geospatial applications exhibit significant distortions in polar imagery due to several reasons: 1) distortions in the source imagery, 2) incompatible tessellation algorithms at the poles, and 3) map projections. We are leveraging new tessellation algorithms and reprojecting data using projections that are better suited for Lunar poles. The goal is to seamlessly switch to polar projections while maintaining 3D view and navigation.

Maps↗

Alternate Methods of Model Reduction to Avoid Dynamic Modal Truncation Error

Loads analysis is traditionally performed using dynamically reduced models, which provides the benefit to reduce run time. If the reduced model frequency cutoff is not chosen appropriately, the model will lack the dynamic content required to fully represent the response of the non-reduced model. Standard guidance for reduced model frequency content, provided in NASA-STD-5002, is to solve fixed base modes up to a minimum of 1.5x the model frequency content of interest and to employ static modal truncation methods such as residual vectors, the mode acceleration method, and the residual flexibility method to account for the truncated flexibility of the missing modes. Fixed base modes require the selection of a set of degrees of freedom to be constrained which, if not properly selected, may affect the accuracy of the reduced model by excluding some of the dynamic characteristic of the full model. In this case, the standard NASA guidance would be insufficient, but it may not be readily apparent that a portion of the reduced model response is missing. This error was encountered during an independent verification and validation (IV&V) effort, where it was observed that the resulting dynamic response was lower than the inline analysis. In this specific case, despite following the standard NASA model reduction guidelines in the selection of the frequency cutoff, the inline model still did not fully capture the necessary dynamic content. As part of the IV&V, an alternate reduction methodology was employed using an unconstrained mode acceleration method. The original model initially performed a constrained reduction to twice the frequency content of interest before doing a free-free run, employing the mode acceleration method to account for the truncated modes. In contrast in the IV&V, the unconstrained model reduced directly to the needed frequency content of the free-free run, avoiding any interactions between constraints and dynamic content. To verify the model, the original reduction methodology was used to generate a series of Hurty-Craig-Bampton reductions, each with a higher frequency cutoff than the previous. The results were shown to converge once the frequency cutoff increased past eight times the frequency content of interest. At the same frequency cutoff, the results of the Hurty-Craig-Bampton model converged with the results of the unconstrained mode acceleration model. This comparative study provided confidence that the results of the unconstrained modal acceleration reduced model were correct and that the Hurty-Craig-Bampton needed to increase its frequency cutoff to fully capture the dynamic response.

Erin Simmons↗

Alternate Methods of Model Reduction to Avoid Dynamic Modal Truncation Error

Loads analysis is traditionally performed using dynamically reduced models, which provides the benefit to reduce run time. If the reduced model frequency cutoff is not chosen appropriately, the model will lack the dynamic content required to fully represent the response of the non-reduced model. Standard guidance for reduced model frequency content, provided in NASA-STD-5002, is to solve fixed base modes up to a minimum of 1.5x the model frequency content of interest and to employ static modal truncation methods such as residual vectors, the mode acceleration method, and the residual flexibility method to account for the truncated flexibility of the missing modes. Fixed base modes require the selection of a set of degrees of freedom to be constrained which, if not properly selected, may affect the accuracy of the reduced model by excluding some of the dynamic characteristic of the full model. In this case, the standard NASA guidance would be insufficient, but it may not be readily apparent that a portion of the reduced model response is missing. This error was encountered during an independent verification and validation (IV&V) effort, where it was observed that the resulting dynamic response was lower than the inline analysis. In this specific case, despite following the standard NASA model reduction guidelines in the selection of the frequency cutoff, the inline model still did not fully capture the necessary dynamic content. As part of the IV&V, an alternate reduction methodology was employed using an unconstrained mode acceleration method. The original model initially performed a constrained reduction to twice the frequency content of interest before doing a free-free run, employing the mode acceleration method to account for the truncated modes. In contrast in the IV&V, the unconstrained model reduced directly to the needed frequency content of the free-free run, avoiding any interactions between constraints and dynamic content. To verify the model, the original reduction methodology was used to generate a series of Hurty-Craig-Bampton reductions, each with a higher frequency cutoff than the previous. The results were shown to converge once the frequency cutoff increased past eight times the frequency content of interest. At the same frequency cutoff, the results of the Hurty-Craig-Bampton model converged with the results of the unconstrained mode acceleration model. This comparative study provided confidence that the results of the unconstrained modal acceleration reduced model were correct and that the Hurty-Craig-Bampton needed to increase its frequency cutoff to fully capture the dynamic response.

Erin Simmons↗

Crowd-Sourced Technology Challenge for Improving Visual Color Detection of Hydrazine and Monomethylhydrazine Vapors in Spacecraft Environments

NASA currently uses a visual colorimetric detection method for potential hydrazine, monomethylhydrazine (MMH), or unsymmetrical dimethylhydrazine (UDMH) contamination in the International Space Station. Astronauts exposed to propellants or their residues during extravehicular activities may transfer contaminants into the airlock. The colorimetric detection method employs the Contamination Detection Kit (CDK), which uses a potassium tetrachloroaurate redox reaction with the propellant hydrazine vapors and a color comparison card to determine airborne concentrations. Seeking ideas for improvement, the NASA Tournament Lab (NTL) crowdsourced a way to tackle the challenge of detecting hydrazine and MMH vapors using colorimetric detection methods. This Rid the Rocket competition drew over 200 participants and 20 submissions from around the world proposing innovative ways to develop a new chemical colorimetric detection method for hydrazine and MMH vapors on spacecraft. Using a phased approach to evaluate contestants, NASA eventually narrowed the field to five finalists from the United States, Romania, Taiwan, and India. Concept papers and hardware submissions were judged on feasibility, creativity, and ability to detect hydrazine and MMH vapors before being sent to the NASA White Sands Test Facility for laboratory evaluation. Finalists employed variations of sampling methods and color-detection chemistry using a variety of sampling pumps and indicator pads or solutions—including those employing potassium or hydrogen tetrachloroaurate, para -dimethylaminobenzaldehyde (PDAB), and modifiers including sodium metasilicate and cetyltrimethylammonium bromide—to enhance gold nanoparticle formation and surface plasmon resonance (SPR) resulting in visual blue to purple color development. This paper presents a summary of the crowdsourced submissions and results of laboratory testing.

Crowd-Sourced↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

U.S. Aerospace and Aviation Industry: A State-By-State Analysis

President George W. Bush and the Congress created the Commission on the Future of the United States Aerospace Industry to evaluate the current and future health of the industry as well as the challenges that lie ahead for the U.S. workforce and the economy. To accomplish our mission, we commissioned a study on the economic impact of the aerospace industry nationally and on a state-by-state basis, using the best available U.S. government data. This study sought to define the core of the aerospace industry. The resulting data represents that very core those jobs specifically tied to commercial and civilian aerospace. U.S. Aerospace and Aviation: A State-by-State Analysis examines the civilian and commercial aerospace and aviation industry by employment, wages, payroll, and establishments. The report found that the U.S. civilian and commercial aerospace and aviation industry employed over 2 million workers in 2001.

STATES (UNITED STATES)↗