SEARCH · Engineering Papers

Results for “utility impact evaluation”

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

SourceResearch connectionFrom yearThrough year

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

504 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Diagnosing Aircraft Icing Potential from Satellite Cloud Retrievals

The threat for aircraft icing in clouds is a significant hazard that routinely impacts aviation operations. Accurate diagnoses and forecasts of aircraft icing conditions requires identifying the location and vertical distribution of clouds with super-cooled liquid water (SLW) droplets, as well as the characteristics of the droplet size distribution. Traditional forecasting methods rely on guidance from numerical models and conventional observations, neither of which currently resolve cloud properties adequately on the optimal scales needed for aviation. Satellite imagers provide measurements over large areas with high spatial resolution that can be interpreted to identify the locations and characteristics of clouds, including features associated with adverse weather and storms. This paper describes new techniques for interpreting cloud products derived from satellite data to infer the flight icing threat to aircraft. For unobscured low clouds, the icing threat is determined using empirical relationships developed from correlations between satellite imager retrievals of liquid water path and droplet size with icing conditions reported by pilots (PIREPS). For deep ice over water cloud systems, ice and liquid water content (IWC and LWC) profiles are derived by using the imager cloud properties to constrain climatological information on cloud vertical structure and water phase obtained apriori from radar and lidar observations, and from cloud model analyses. Retrievals of the SLW content embedded within overlapping clouds are mapped to the icing threat using guidance from an airfoil modeling study. Compared to PIREPS and ground-based icing remote sensing datasets, the satellite icing detection and intensity accuracies are approximately 90% and 70%, respectively, and found to be similar for both low level and deep ice over water cloud systems. The satellite-derived icing boundaries capture the reported altitudes over 90% of the time. Satellite analyses corresponding to the time and location of several recent aviation accidents and with icing PIREPS are also presented that reveal skill in identifying severe icing conditions. These results demonstrate the utility of satellite cloud retrievals for quantitatively diagnosing the potential for icing conditions on temporal and spatial scales that should be useful to the aviation community. Plans are being developed to deliver these new satellite products to the GOES-R Proving Ground in the near future so that they can be evaluated in operational applications.

Smith, William L., Jr.↗

NASA NTRS1976

An exploratory study to determine the integrated technological air transportation system ground requirements of liquid-hydrogen-fueled subsonic, long-haul civil air transports

An exploratory study was conducted to evaluate the feasibility and cost impact of implementing a liquid-hydrogen-fueled (LHz), air transport support system at Chicago-O’Hare International Airport. The study assumed an available supply of hydrogen to the airport boundary. It also assumed that the present widebody fleet at O’Hare was replaced by 400 passenger, LH2-fueled aircraft operating on the current frequency schedule and route network. The LH2 aircraft configurations used in this study were developed during a previous NASA-sponsored study. A baseline air terminal concept was developed which permitted airlines and the airport to operate JP- or LH2-fueled aircraft at common terminal gates. The concept included installation of a hydrogen liquefaction and storage facility on airport property, as well as the fuel distribution system. The capital investment and hydrogen-related operating costs to the airlines were estimated. The study concluded that the system would be technically and operationally feasible at O’Hare, and that economics would be a prime factor in future decisions regarding the use of LH2 as an air transport fuel. Research and technology recommendations were made to improve component efficiency/reliability. The study also concluded that additional investigation should be conducted on the implementation of a nationwide or worldwide LH2 air transportation system, to understand fully the operational and economic implications of the total system, before decisions could be made relative to committing to this option.

Hydrogen↗

NASA NTRS1999

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

NASA NTRS2018

CubeX: A Compact X-Ray Telescope Enables Both X-Ray Fluorescence Imaging Spectroscopy and Pulsar Timing Based Navigation

This paper describes the benefits of a miniaturized X-ray telescope payload in the context of a lunar mission. The first part describes the payload in detail, the second part summarizes a small satellite mission concept that utilizes its compact form factor and performance. The CubeX instrument can be used for both X-ray fluorescence (XRF) imaging spectroscopy and X-ray pulsar timing based navigation (XNAV). Using our recent technological advances in X-ray optics and sensors, CubeX combines high angular resolution (<1 arcminutes) Miniature Wolter-I X-ray optics (MiXO) with a common focal plane consisting of high spectral resolution (<150 eV at 1 keV) CMOS X-ray sensors and a high timing resolution (< 1 usec) SDD X-ray sensor. This novel combination of the instruments enables both XRF measurements and XNAV operations without moving parts. The high angular resolution of the MiXO opens a wide range of orbital configurations for observation. Given that performance, the instrument has unprecedented small volume (~1×1×6U), mass (<6 kg), and power (<9W) requirements and opens a wide range of applications for a variety of targets and missions including NEOs and Martian moons. In this paper we illustrate one potential application for a lunar mission concept: The elemental composition of the Moon holds keys to understanding the origin and evolution of both the Moon and the Earth. X-ray fluorescence (XRF), induced either by solar X-ray flux or energetic ions, carries decisive signatures of surface elemental composition. X-ray observations, therefore, give a unique, powerful diagnostic tool for remotely determining elemental abundances including major rock forming elements such as Mg, Al, Na, Si, Fe, and Ca. Through high-resolution XRF imaging spectroscopy, CubeX searches for small patches of elusive lower crust and mantle material excavated within and around impact craters. CubeX identifies regional compositional variations and allows straightforward comparison of elemental distributions with the surface topography from LRO and the gravity data from GRAIL. The elemental compositions of the lower crust and the mantle are sensitive to the conditions of the giant impact which led to the Moon's formation and the subsequent lunar magma ocean (LMO), and thus they are key missing pieces in understanding the formation and early evolution of the Moon. In between XRF observations, CubeX also leverages the technology of high resolution X-ray imaging and time series measurements to conduct XNAV operations and evaluate their performance. Deep space navigation is a critical issue for small planetary missions. XNAV can enable low-cost autonomous deep-space navigation, and has the potential to greatly assist, or even outperform, NASA's Deep Space Network (DSN) or ESA's European Space Tracking (ESTRACK). CubeX is designed to perform sequential observations of 3-4 millisecond pulsars (MSPs) to solve the spacecraft trajectory for absolute navigation, and explore the remote sensing capability of XNAV. In the presented mission concept, the Moon's relative proximity enables a straightforward evaluation of the XNAV performance through DSN.

x-ray spectroscopy↗

NASA NTRSDate not supplied

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗

NASA NTRS2010

Lessons Learned from Performance Testing of Humans in Spacesuits in Simulated Reduced Gravity

Introduction: The overarching objective of the Integrated Suit Test (IST) series is to evaluate suited human performance using reduced-gravity analogs and learn what aspects of an EVA suit system affect human performance. For this objective to be successfully achieved, the testing methodology should be valid and reproducible, and the partial-gravity simulations must be as accurate and realistic as possible. Objectives: To highlight some of the key lessons learned about partial-gravity analogs and testing methodology, and to suggest considerations for optimizing the effectiveness and quality of results of future tests. Methods: Performance testing of suited and unsuited subjects was undertaken in different reduced-gravity analogs including the Space Vehicle Mockup Facility s Partial Gravity Simulator (POGO), parabolic flight on the C-9 aircraft, underwater environments including NASA s Extreme Environment Mission Operations (NEEMO) and the Neutral Buoyancy Lab (NBL), and in field analogs including Desert Research and Technology Studies (RATS), the Haughton Mars Project (HMP), and the JSC Rock Pile. Subjects performed level walking, incline/decline walking, running, shoveling, picking up and transferring rocks, kneeling/standing, and task boards. Lessons Learned Analogs: No single analog will properly simulate all aspects of the true partial-gravity environment. The POGO is an ideal environment from the standpoint that there are no time limits or significant volumetric constraints, but it does have several limitations. It allows only 2 translational degrees of freedom (DOF) and applies true partial-gravity offload only through the subject s center of gravity (CG). Also, when a subject is doing non-stationary tasks, significant overhead inertia from the lift column seems to have a negative impact on performance. Parabolic flight allows full translational and rotational DOF and applies offload to all parts of the body, but the simulation lasts less than 30 seconds. When this is coupled with the volumetric constraints of the plane, both task selection and data collection options are significantly limited. The underwater environments also allow all 6 DOF and allow off-loading to be applied throughout the body, but the data collection capabilities are limited to little more than subjective ratings. In addition, water drag negatively affects performance of tasks requiring dynamic motion. Field analogs provide the ability to simulate lunar terrain and more realistic mission-like objectives, but all of them operate at 1-g, so suited human performance testing generally must utilize a reduced-mass or "mockup" suit, depending on study objectives. In general, the ground-based overhead-suspension partial-gravity analogs like POGO allow the most diverse data collection methods possible while still simulating partial gravity. However, as currently designed, the POGO has significant limitations. Design of the Active Response Gravity Offload System (ARGOS) has begun and is focusing on adding full x,y,z translational DOF, improved offload accuracy, increased lift capacity, and active control of the x and y axes to minimize offload system inertia. Additionally, a new gimbal is being designed to reduce mass and inertia and to be able to work with different suits, as the current gimbal only supports suited testing with the Mark III Technology Demonstrator Suit (MKIII).

Norcross, Jason R.↗

NASA NTRS2016

Hybrid Propulsion In-Situ Resource Utilization Test Facility Results for Performance Characterization

Hybrid propulsion presents a promising alternative to conventional systems for in-space propulsion applications using In-Situ Propellant Production (ISPP). A hybrid propulsion test facility has been built at the NASA Jet Propulsion Laboratory (JPL) with the capability of evaluating ISPP oxidizers (gaseous mixtures of O2 and CO2) and modifying fuel type and chamber geometry. This work is the result of a continuing e ort at JPL to characterize the burn characteristics and performance of different propellant combinations and motor configurations in order to re ne and support preliminary designs. Twelve tests using gaseous oxygen and two paraffin based fuels (Black Paraffin and SP1X) from the 2015 campaign are presented. Different reduction techniques are compared in an e ort to determine the uncertainty associated with deriving a regression rate law and evaluating c* efficiency with the collected data.Hybrid motors are mechanically and chemically simple propulsion systems that have the potential forlong-term storability with low gross mass. The oxidizer (such as gaseous O2) can be produced in-situ,signicantly reducing the landed mass. The solid, inert fuel is stored in the combustion chamber in its most dense state, further reducing the total structural mass and volume. Considerable eorts have been focusedon ISPP for Environmental Control and Life Support Systems (ECLSS) and human exploration of Mars.At optimum oxidizer to fuel ratio (O/F), about 70% of the total propellant mass could be produced in-situ.Paran based fuels exhibit a fast regression rate behavior due to the formation of a thin liquid layer on thefuel surface.1 A high burn rate is necessary in order to reach a medium to high thrust range. In order toachieve an O/F ratio close to optimum performance (between 2.1 and 2.7 for the O2/Paran combination),the fuel mass ow needs to be suciently large compared to the oxidizer ow. The slow regression rateof classical hybrid rocket fuels leads to long fuel grain designs or complicated multi-port structures. Singleport, high regression rate fuels have the potential to reduce residual propellant mass and have a grainlength to diameter ratio (L/D) more favorable for packaging. However, hybrid motors remain at a fairly lowTechnology Readiness Level (TRL), and solutions to accurately measure and increase performance are stilllacking, especially for small-scale motors that are constrained to tight packaging.A total of 12 hotres were performed and analyzed in this campaign. All of these tests were conductedat pressures above the critical pressure of paran wax, approximately 6.7 bar (97 psi).2 Three dierentmotor congurations were explored while conserving a single cylindrical port fuel grain. Additionally, sincethrottling is a crucial advantage of hybrid rocket propulsion systems, low oxidizer mass ows (60 - 100 g/s)and uxes (5 - 30 g/cm2/s) were investigated to determine the performance impact in this regime.

Mechentel, F. S.↗

NASA NTRSDate not supplied

James Webb Space Telescope Navigation Optimization Challenges

The James Webb Space Telescope (JWST) is a NASA flagship mission launched on December 25, 2021. The early orbit phase was highlighted by three midcourse correction burns that were performed to maneuver the vehicle into a libration orbit at the second Sun-Earth-Moon (SEM) libration point (L2). During the coast to L2, the vehicle’s primary observatory and sunshield were deployed. Once the tennis-court-sized sunshield was unfurled and tensioned into place, the large area exposed to solar radiation pressure (SRP) dictated that the SRP force model would need to account for the vehicle’s geometry, reflective properties, and orientation relative to the Sun. Following the insertion of JWST into its L2 libration orbit on January 24, 2022, the commissioning phase for the vehicle’s observatory commenced along with the first cycle of routine station-keeping maneuvers. The NASA Goddard Space Flight Center’s Flight Dynamics Facility (FDF) provides navigation services for the JWST mission. The FDF serves as the prime or backup navigation operations center for more than 30 active missions spanning a wide array of flight regimes performing orbit determination, maneuver planning, trajectory optimization, and tracking data evaluation. Definitive orbit determination for JWST is performed by the FDF using an Extended Kalman Filter (EKF) to estimate the vehicle’s trajectory and reflective properties using Deep Space Network (DSN) TRK-2-34 tracking data and spacecraft attitude telemetry. For the purposes of orbit prediction and station-keeping maneuver targeting, the sensitivity of the SRP force to the vehicle’s orientation requires that predictive attitude information must be incorporated to generate accurate orbit predictions and maneuver plans. Predictive attitude plans are provided to the FDF by the JWST Spacecraft Operations Center (SOC). Short-term (28 day) orbit predictions are propagated using a Short Range Attitude Plan (SRAP) to model the vehicle’s future attitude states for up to a week. SRAP attitude data is highly reliable, as it reflects finalized planned attitude states that the vehicle will be commanded to attain during the ensuing week. A Long Range Attitude Plan (LRAP) can be implemented to model predicted attitude states beyond one week, but these predictions are tentative and subject to revision due to changes in science operations plans. Instead, a conservative approach of applying a Sun-Pointing Neutral (SPN) attitude configuration is utilized for long-term (2 year) orbit predictions where reliable planned attitude data is not available. SPN attitude mode aligns the net SRP force along the JWST-to-Sun vector where it becomes independent of the vehicle’s angle of rotation about this vector, defined as the Sun yaw. Given the dynamic instability of the libration orbit, station-keeping thrust must be applied in either the sunward or anti-sunward direction to balance the resulting orbit. The attitude constraints of the vehicle also impose limits on the available pointing directions for station-keeping thrusters. The thrusters cannot be aligned with the optimal pointing direction for sunward maneuvers, rendering sunward maneuvers to be significantly less fuel efficient than anti-sunward maneuvers. Consequently, station-keeping maneuvers must be targeted to balance the orbit while ensuring that the next maneuver will also be performed in the anti-sunward direction to optimize propellant usage. When targeting a station-keeping maneuver, the attitude and SRP modeling configuration which is applied to the predicted post-maneuver trajectory will dictate the direction of the subsequent station-keeping maneuver, assuming nominal propulsion system performance. If the predicted post-maneuver attitude states result in an SRP model which under-predicts the cumulative SRP impact on the post-maneuver orbit, a targeted anti-sunward maneuver will be larger than necessary and the next maneuver will need to be executed sunward in order to compensate. In contrast, an anti-sunward maneuver which is targeted using an SRP model which over-predicts the cumulative SRP impact will achieve station-keeping while ensuring that the next maneuver will likewise be performed anti-sunward. For this reason, station-keeping maneuvers are targeted while applying the SPN attitude mode, as this mode entails the largest possible SRP area cross-section and therefore a larger modeled cumulative SRP impact. This paper documents the NASA Goddard Space Flight Center’s FDF support for JWST on-orbit operations and the analysis projects undertaken utilizing the experiences and data accumulated throughout the first full year of routine science operations. The results of these analysis efforts have been used to implement improvements to orbit prediction accuracy and maneuver efficiency which have the potential to prolong the lifespan of JWST to continue to conduct ground-breaking infrared astronomy for decades to come.

Flight Dynamics↗

NASA NTRSDate not supplied

Rational Design of Nanoplasmonic Array Geometries for Biosensing

Background: Molecular diagnostics provide early and accurate diagnosis, which is essential for the prevention and treatment of infectious as well as chronic diseases. These tests are designed to detect disease-specific bioanalytes such as nucleic acid (DNA or RNA) or protein (antigens, antibodies) biomarkers. In the context of infectious disease diagnosis, nucleic acid-based detection methods are known to provide more specific and sensitive results. Here, the presence of a unique sequence belonging to the pathogenic genomic material is targeted to identify species, organism, genera and/or antimicrobial resistant gene markers. The majority of the common nucleic acid based diagnostic techniques require amplification (polymerase chain reaction, isothermal amplification etc.) of the pathogenic genetic material prior to detection impacting diagnostic speed, complexity, and cost thereby limiting ease of use. Thus, the development of simplified nucleic acid-based diagnostics that can be even used in resource-poor settings may hugely benefit patients across the globe. Nanopath is a molecular diagnostics company utilizing a solid-state nanosensor to enable sequence-specific detection of target nucleic acids without the need of amplification. These nanostructures enable ultra-sensitive biomarker detection using geometric, feature-dependent properties highly dependent on the local dielectric environment, allowing them to be sensitive to low concentration binding events. This paper describes an application of this approach to provide highly relevant clinical information within a single doctor’s office visit. Intro: The Nanopath team is in collaboration with NASA (National Aeronautics and Space Administration) and NIST (National Institute of Standards and Technology) to push the bounds of the fundamental physics associated with their biosensing platform. The ability of metals to support electromagnetic surface waves gives rise to surface plasmons when optically illuminated. This property, and its strong sensitivity to changes in the local refractive index, allows for the use of metal nanoparticles as ultra-sensitive transducers. In prior work by members of this team, ensembles of randomly oriented nanoparticles (i.e., colloidal nanorods dispersed on chip) were employed for sequence-specific nucleic acid sensing (1-3). While these particle sensors have the advantage of rapid fabrication, they suffer from low sensitivity and quality factor due to the random particle dispersity. In contrast, in this study we employ ordered array nanoparticle ensembles which can be used to improve sensor sensitivity and figure-of-merit. Study Methods Overview: In this talk, we detail the results of sensing experiments and computational simulations to outline a rational design of the structure of these plasmonic nanoparticle arrays for biomolecular sensing. Through simulation and experiment, we iteratively tailor nanostructure dimension to provide high quality signal and large resonance shifts upon modeled nucleic acid binding. In particular, full-wave electromagnetic simulations were conducted using Lumerical photonic simulation software in which periodic boundary conditions were applied in the x- and y- dimensions for each of the nanoplasmonic sensor geometries. To simulate the resonance response to changes in the bulk solution in contact with the sensor surface, the refractive index of the surrounding media was changed appropriately. Nucleic acid hybridization events were modeled using either using spherical structures approximating the relevant radius of genomic material as estimated by polymer models, or as conformal layers with the known refractive indices for nucleic acids. On the basis of initial simulations, nanosensors were fabricated using traditional electron-beam lithography protocols at NIST. To evaluate consensus between simulations and experiments, bulk sensing experiments were carried out in which the resonance peaks were obtained by submerging the sensors in refractive index standards. Key nanosensor characteristics including resonance peak locations, resonance peak shifts as a function of refractive index, and figure of merit (FOM) of extinction curves were examined between the experimental and simulation results prior to proceeding with simulations on additional geometries and more complex solution conditions, and further device fabrication. This iterative process is repeated toward a rational design of nanoplasmonic array geometries for biosensing optimizing response for targeted disease detection. In summary, this study puts forth a methodology for rational design and characterization of regularly spaced nanoparticle arrays for optics-based biosensing. The results of this study will allow for more informed design of nanostructure geometries towards sequence-specific nucleic acid detection. These improved designs have the potential to improve clinical sensitivity and limit-of-detection across disease indication.

sensor↗

NASA NTRS2017

Computational Modeling to Limit the Impact Displays and Indicator Lights Have on Habitable Volume Operational Lighting Constraints

NASA has demonstrated an interest in improving astronaut health and performance through the installment of a new lighting countermeasure on the International Space Station. The Solid State Lighting Assembly (SSLA) system is designed to positively influence astronaut health by providing a daily change to light spectrum to improve circadian entrainment. Unfortunately, existing NASA standards and requirements define ambient light level requirements for crew sleep and other tasks, yet the number of light-emitting diode (LED) indicators and displays within a habitable volume is currently uncontrolled. Because each of these light sources has its own unique spectral properties, the additive lighting environment ends up becoming something different from what was planned or researched. Restricting the use of displays and indicators is not a solution because these systems provide beneficial feedback to the crew. The research team for this grant used computer-based computational modeling and real-world lighting mockups to document the impact that light sources other than the ambient lighting system contribute to the ambient spectral lighting environment. In particular, the team was focused on understanding the impacts of long-term tasks located in front of avionics or computer displays. The team also wanted to understand options for mitigating the changes to the ambient light spectrum in the interest of maintaining the performance of a lighting countermeasure. The project utilized a variety of physical and computer-based simulations to determine direct relationships between system implementation and light spectrum. Using real-world data, computer models were built in the commercially available optics analysis software Zemax Optics Studio(c). The team also built a mockup test facility that had the same volume and configuration as one of the Zemax models. The team collected over 1200 spectral irradiance measurements, each representing a different configuration of the mockup. Analysis of the data showed a measurable impact on ambient light spectrum. This data showed that obvious design techniques exist that can be used to bind the ambient light spectrum closer to the planned spectral operating environment for the observer's eye point. The following observations should be considered when designing an operational environment that is dominated by computer displays. When more light is directed into the field of view of the observer, the greater the impact it will make on various human factors issues that depend on spectral shape and intensity. Because viewing angle has a large part to play in the amount of light flux on the crewmember's retina, beam shape, combined with light source location is an important factor for determining percent probable incident flux on the observer from any combination of light sources. Computer graphics design and display lumen output are major factors influencing the amount of spectrally intense light projected into the environment and in the viewer's direction. Use of adjustable white point display software was useful only if the predominant background color was white and if it matched the ambient light system's color. Display graphics that used a predominantly black background had the least influence on unplanned spectral energy projected into the environment. Percent reflectance makes a difference in total energy reflected back into an environment, and within certain architectural geometries, reflectance can be used to control the amount of a light spectrum that is allowed to perpetuate in the environment. Data showed that room volume and distance from significant light sources influence the total spectrum in a room. Smaller environments had a homogenizing effect on total light spectrum, whereas light from multiple sources in larger environments was less mixed. The findings indicated above should be considered when making recommendations for practice or standards for architectural systems. The ambient lighting system, surface reflectance, and display and indicator implementation all factor into the users' spectral environment. A variety of low-cost solutions exist to mitigate the impact of light from non-architectural lighting systems, and much potential for system automation and integration of display systems with the ambient environment. This team believes that proper planning can be used to avoid integration problems and also believes that human-in-the-loop evaluations, real-world test and measurement, and computer modeling can be used to determine how changes to a process, display graphics, and architecture will help maintain the planned spectral operating lighting environment.

Clark, T. A.↗

NASA NTRS2017

Biosignatures in the Context of Low Energy Flux

Many of the features that are thought of as biosignatures - including the mediation of chemical and physical processes with speed, specificity, and selectivity - result directly or indirectly from life's unique capability to mediate and direct energy flux. As such, it is important to consider the impact that differences in energy flux may have on the quantity and quality of evidence for life. Earth differs from every other body in our solar system in the magnitude of biologically-usable energy flux into a liquid water environment. On a global basis, the capture of light energy into photosynthesis and the flux of chemical energy represented in the products of that photosynthesis (organic material + O2) are about six and four orders of magnitude larger, respectively, than the flux of energy represented in geochemical sources. Our conception of what an inhabited world "looks like" and our intuition about how to search for life are based in this high-energy context. Energy fluxes on worlds beyond Earth may be better approximated by the million-fold smaller flux provided to Earth's biosphere by geochemical sources. As a result, the nature, abundance, and quality of evidence for life that could be expected on an inhabited extraterrestrial world within our solar system may differ profoundly from that found on Earth. Understanding this potential difference in quantitative terms provides important context for the formulation of life detection strategies. The influence of energy flux on biosignatures can be evaluated through reference to the two basic purposes into which life partitions energy flux: (1) Life expends energy to sustain existing biomass in a metabolic steady state (metabolically functional but non-growing). The formal representation of this relationship in the traditional microbiology literature equates biomass directly with energy flux. The direct implication is that worlds having lower energy flux will have correspondingly lower potential to support biomass. Life detection strategies that directly target extant organisms should therefore be prepared to encounter average biomass densities that may be many orders of magnitude smaller than those found in most of Earth's surface environments (2) Life expends energy to synthesize new biomass. An end-member case in which new biomass is created at the energy-limited rate and the corresponding cells are immediately destroyed (so that the energy partitioned to cell maintenance is minimized) establishes an upper bound on the rate at which biological material can enter a bulk global pool. For a specified bulk concentration [i] of any particular biological compound, i, or for biologically produced matter overall, this synthesis rate, R (sub i), defines a characteristic time scale tau (sub i) equals [i] divided by R (sub i). tau (sub i) can be thought of as (a) the minimum time required for biosynthesis to yield a specific bulk concentration (e.g., a detection threshold) of i, and (b) the average residence time of i within a bulk pool when [i] is held in steady state through a balance between biosynthesis and attrition by physical, chemical, or biological consumption. tau (sub i) becomes an important quantity in considering the potential utility of enantiomeric excess (as a product of homochiral biosynthesis) as a biosignature. Spontaneous racemization of amino acids acts to "erase" the signature of homochiral synthesis over time scales that may range from hundreds to hundreds of thousands of years, depending on temperature. For environments in which low energy flux translates to low rates of biosynthesis, including the synthesis of homochiral amino acids, amino acid residence times in pools having detectable concentrations may compare to or significantly exceed the time scale for racemization. This and similar consequences of long residence times should be considered in the formulation of life detection strategies based on detection of biologically-produced species. Fluxes of biologically-useful energy on potentially habitable worlds within our solar system are, at present, not well constrained. Improving such constraint has the potential to inform priorities in the formulation and targeting of search-for-life strategies, based on the implications of energy flux for the abundance and quality of biosignatures overall, and in specific categories.

ENERGY FLUX↗

NASA NTRSDate not supplied

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

NASA NTRS2007

An Experimental System for a Global Flood Prediction: From Satellite Precipitation Data to a Flood Inundation Map

Floods impact more people globally than any other type of natural disaster. It has been established by experience that the most effective means to reduce the property damage and life loss caused by floods is the development of flood early warning systems. However, advances for such a system have been constrained by the difficulty in estimating rainfall continuously over space (catchment-. national-, continental-. or even global-scale areas) and time (hourly to daily). Particularly, insufficient in situ data, long delay in data transmission and absence of real-time data sharing agreements in many trans-boundary basins hamper the development of a real-time system at the regional to global scale. In many countries around the world, particularly in the tropics where rainfall and flooding co-exist in abundance, satellite-based precipitation estimation may be the best source of rainfall data for those data scarce (ungauged) areas and trans-boundary basins. Satellite remote sensing data acquired and processed in real time can now provide the space-time information on rainfall fluxes needed to monitor severe flood events around the world. This can be achieved by integrating the satellite-derived forcing data with hydrological models, which can be parameterized by a tailored geospatial database. An example that is a key to this progress is NASA's contribution to the Tropical Rainfall Measuring Mission (TRMM), launched in November 1997. Hence, in an effort to evolve toward a more hydrologically-relevant flood alert system, this talk articulates a module-structured framework for quasi-global flood potential naming, that is 'up to date' with the state of the art on satellite rainfall estimation and the improved geospatial datasets. The system is modular in design with the flexibility that permits changes in the model structure and in the choice of components. Four major components included in the system are: 1) multi-satellite precipitation estimation; 2) characterization of land surface including digital elevation from NASA SRTM, topography-derived hydrologic parameters such as flood direction. flow accumulation, basin, and river network etc.; 3) spatially distributed hydrological models to infiltrate rainfall and route overland runoff; and 4) an implementation interface to relay thc input data to the models and display the flood inundation results to the users and decision-makers. Early results appear reasonable in terms of location and frequency of events. Case studies of this experimental system are evaluated with surface runoff data and other river monitoring systems. such as Dartmouth Flood Observatory's "Surface Water Watch" array of river reaches that are measured daily via other satellite remote sensing data. A major outcome of this progress will be the availability of a global overview of flood alerts that should consequently improve the performance of Decision Support System. We expect these developments in utilization of satellite remote sensing technology to offer a practical solution to the challenge of building a cost-effective early warning system for data scarce and under-developed areas.

Adler, Robert↗

NASA NTRS1998

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗

NASA NTRS2003

Catalyzed Combustion In Micro-Propulsion Devices: Project Status

In recent years, there has been a tendency toward shrinking the size of spacecraft. New classes of spacecraft called micro-spacecraft have been defined by their mass, power, and size ranges. Spacecraft in the range of 20 to 100 kg represent the class most likely to be utilized by most small sat users in the near future. There are also efforts to develop 10 to 20 kg class spacecraft for use in satellite constellations. More ambitious efforts will be to develop spacecraft less than 10 kg, in which MEMS fabrication technology is required. These new micro-spacecraft will require new micro-propulsion technology. Although micro-propulsion includes electric propulsion approaches, the focus of this proposed program is micro-chemical propulsion which requires the development of microcombustors. As combustors are scaled down, the surface to volume ratio increases. The heat release rate in the combustor scales with volume, while heat loss rate scales with surface area. Consequently, heat loss eventually dominates over heat release when the combustor size becomes smaller, thereby leading to flame quenching. The limitations imposed on chamber length and diameter has an immediate impact on the degree of miniaturization of a micro-combustor. Before micro-combustors can be realized, such a difficulty must be overcome. One viable combustion alternative is to take advantage of surface catalysis. Micro-chemical propulsion for small spacecraft can be used for primary thrust, orbit insertion, trajectory-control, and attitude control. Grouping micro-propulsion devices in arrays will allow their use for larger thrust applications. By using an array composed of hundreds or thousands of micro-thruster units, a particular configuration can be arranged to be best suited for a specific application. Moreover, different thruster sizes would provide for a range of thrust levels (from N s to mN s) within the same array. Several thrusters could be fired simultaneously for thrust levels higher than the basic units, or in a rapid sequence in order to provide gradual but steady low-g acceleration. These arrays of micro-propulsion systems would offer unprecedented flexibility and redundancy for satellite propulsion and reaction control for launch vehicles. A high-pressure bi-propellant micro-rocket engine is already being developed using MEMS technology. High pressure turbopumps and valves are to be incorporated onto the rocket chip . High pressure combustion of methane and O2 in a micro-combustor has been demonstrated without catalysis, but ignition was established with a spark. This combustor has rectangular dimensions of 1.5 mm by 8 mm (hydraulic diameter 3.9 mm) and a length of 4.5 mm and was operated at 1250 kPa with plans to operate it at 12.7 MPa. These high operating pressures enable the combustion process in these devices, but these pressures are not practical for pressure fed satellite propulsion systems. Note that the use of these propellants requires an ignition system and that the use of a spark would impose a size limitation to this micro-propulsion device because the spark unit cannot be shrunk proportionately with the thruster. Results presented in this paper consist of an experimental evaluation of the minimum catalyst temperature for initiating/supporting combustion in sub-millimeter diameter tubes. The tubes are resistively heated and reactive premixed gases are passed through the tubes. Tube temperature and inlet pressure are monitored for an indication of exothermic reactions and composition changes in the gases.

Sung, C. J.↗

NASA NTRS2016

The Impact of New Trends in Satellite Launches on Orbital Debris Environment

The main goal of this study is to examine the impact of new trends in satellite launch activities on the orbital debris environment and collision risk. Starting from the launch of the first artificial satellite in 1957, space borne technology has become an indispensable part of our lives. More than 6,000 satellites have been launched into Earth orbit. Though the annual number of satellites launched stayed flat for many decades, the trend has recently changed. The satellite market has been undergoing a major evolution with new space companies replacing the traditional approach of deploying a few large, complex and costly satellites with an approach to use a multitude of smaller, less complex and cheaper satellites. This new approach creates a sharp increase in the number of satellites and so the historic trends are no longer representative. As a foundation for this study, a scenario for satellite deployments based on the publicly announced future satellite missions has been developed. These constellation-deploying companies include, but are not limited to, Blacksky, CICERO, EROS, Landmapper, Leosat, Northstar, O3b, OmniEarth, OneWeb, Orbcomm, OuterNet, PlanetIQ, Planet Labs, Radarsat, RapidEye Next Generation, Sentinel, Skybox, SpaceX, and Spire. Information such as the annual number of launches, the number of orbital planes to be used by the constellation, as well as apogee, perigee, inclination, spacecraft mass and area were included or approximated. Besides the production of satellites, a widespread ongoing effort to enhance orbital injection capabilities will allow delivery of more spacecraft more accurately into Earth orbits. A long list of companies such as Microcosm, Rocket Lab, Firefly Space Systems, Sierra Nevada Corporation and Arca Space Corporation are developing new launch vehicles dedicated for small satellites. There are other projects which intend to develop interstages with propulsive capabilities which will allow the deployment of satellites into their desired orbits beyond the restrictions of the launch vehicle used. These near future orbital injection technologies are also covered in the developed scenario. Using the above-mentioned background information, this study aims to examine how the orbital debris environment will be affected from the new dynamics of the emerging space markets. We developed a simulation tool that is capable of propagating the objects in a given deployment scenario with variable-sized time-steps as small as one second. Over the course of the run, the software also detects collisions; additional debris objects are then created according to the NASA breakup model and are fed back into the simulation framework. Examining the simulation results, the total number of particles to accumulate in different orbits can be monitored and the number of conjunctions can be tracked to assess the collision risks. The simulation makes it possible to follow the short- and long-term effects of a particular satellite or constellation on the space environment. Likewise, the effects of changes in the debris environment on a particular satellite or constellation can be evaluated. It is authors hope that the results of this paper and further utilization of the developed simulation tool will assist in the investigation of more accurate deorbiting metrics to replace the generic 25-year disposal guidelines, as well as to guide future launches toward more sustainable and safe orbits.

New Space↗

NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

NASA NTRS2013

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗