Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust optimization”

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 1,099 records · Page 61

The Completion of a Geosynchronous Earth Orbit Survey with the Eugene Stansbery-Meter Class Autonomous Telescope

The Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT) is the primary optical sensor used by the NASA Orbital Debris Program Office (ODPO) to statistically characterize the geosynchronous Earth orbit (GEO) debris environment and support future Orbital Debris Engineering Model (ORDEM) releases. The ES-MCAT completed its first optical survey of the GEO region from 2020 to 2022. The primary goal of this survey was to autonomously collect and process GEO data with calculated photometric and astrometric uncertainties. A pointing plan was developed to provide uniform sampling within the region of interest (ROI) while accounting for predicted downtime due to insufficient observing conditions. Detections are autonomously correlated to the Space Surveillance Network (SSN) catalog to determine if objects are correlated targets (CTs) or uncorrelated targets (UCTs), the latter of which are of interest for modeling the GEO orbital debris environment. To assess the size detection sensitivity over time and monitor the general performance of the telescope’s optics and software, the optical throughput of the system and limiting magnitudes are evaluated on a routine basis. While the telescope’s ability to operate autonomously and remotely allowed for the GEO survey to continue throughout the COVID-19 pandemic, travel restrictions hampered routine cleaning of the optics during this time, and the primary mirror degraded enough to require recoating. The mirror was removed in 2022, concluding the first GEO survey. The primary mirror received a new coating designed to be more robust against the harsh environment surrounding the ES-MCAT’s location on Ascension Island, accounting for experience gained during operations over the first GEO survey. In early 2023, the recoated primary mirror was reinstalled, and the second GEO survey was initiated. The primary goal of the second GEO survey is to characterize the evolving GEO debris environment with updated optics, software, and pointing strategies while allowing for the inclusion of non-GEO regimes or those that are outside of the ROI. While the pointing method implemented in the first survey allowed for adequate coverage of the ROI over two years, it has been improved to include pointings that avoid the Moon’s position and the galactic plane to reduce software processing time and maximize the detection capabilities of fainter objects. This method also accounts for the changing weather patterns throughout the year and reduces coverage gaps in the ROI. Provided the success of the first two-year GEO survey using autonomous operations, the ODPO is actively collaborating with the United States Space Force (USSF) to make the ES-MCAT a contributing sensor to the SSN. This paper presents results from the first GEO survey including magnitude distributions and orbital parameters for CTs and UCTs. Details are provided for the automated processing pipeline and the optical system throughput for the previous and current primary mirror coatings. In addition, an updated strategy for the second GEO survey to optimize coverage over the ROI is discussed, as are preliminary results from the ongoing second survey.

Corbin Cruz↗

Assessment of Mars Atmospheric Temperature Retrievals from the Thermal Emission Spectrometer Radiances

Motivated by the needs of Mars data assimilation. particularly quantification of measurement errors and generation of averaging kernels. we have evaluated atmospheric temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) radiances. Multiple sets of retrievals have been considered in this study; (1) retrievals available from the Planetary Data System (PDS), (2) retrievals based on variants of the retrieval algorithm used to generate the PDS retrievals, and (3) retrievals produced using the Mars 1-Dimensional Retrieval (M1R) algorithm based on the Optimal Spectral Sampling (OSS ) forward model. The retrieved temperature profiles are compared to the MGS Radio Science (RS) temperature profiles. For the samples tested, the M1R temperature profiles can be made to agree within 2 K with the RS temperature profiles, but only after tuning the prior and error statistics. Use of a global prior that does not take into account the seasonal dependence leads errors of up 6 K. In polar samples. errors relative to the RS temperature profiles are even larger. In these samples, the PDS temperature profiles also exhibit a poor fit with RS temperatures. This fit is worse than reported in previous studies, indicating that the lack of fit is due to a bias correction to TES radiances implemented after 2004. To explain the differences between the PDS and Ml R temperatures, the algorithms are compared directly, with the OSS forward model inserted into the PDS algorithm. Factors such as the filtering parameter, the use of linear versus nonlinear constrained inversion, and the choice of the forward model, are found to contribute heavily to the differences in the temperature profiles retrieved in the polar regions, resulting in uncertainties of up to 6 K. Even outside the poles, changes in the a priori statistics result in different profile shapes which all fit the radiances within the specified error. The importance of the a priori statistics prevents reliable global retrievals based a single a priori and strongly implies that a robust science analysis must instead rely on retrievals employing localized a priori information, for example from an ensemble based data assimilation system such as the Local Ensemble Transform Kalman Filter (LETKF).

Hoffman, Matthew J.↗

How Well Satellite Remote Sensing Can Inform Surface Level Ozone Production Sensitivity to Precursor Trace Gas Emissions?

Surface-level ozone (O3) pollution mitigation strategies rely on a clear understanding of the atmospheric chemical processes involving O3, NOx (nitric oxide (NO) + nitrogen dioxide (NO2)) and volatile organic compounds (VOCs). A robust spatiotemporal classification of O3 photochemical regimes with relative chemical sensitivity of local O3 formation to emission reductions of NOx (NOx-limited regime) versus VOCs (radical-limited regime) is required to design anthropogenic emission reduction policies in order to improve surface air quality. Numerous studies found that the ratio of satellite-retrieved Vertical Column Densities (VCDs) of formaldehyde (HCHO) to NO2 can diagnose O3 production sensitivity. However, uncertainty still remains about the accuracy of such estimates using satellite-based retrievals. This study conducts an investigation to identify how well satellite data can inform surface-level O3 sensitivity by validating satellite-based estimates against aircraft-based retrievals taken during the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) and Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaigns. For this study, we use HCHO and NO2 retrievals from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and high spatiotemporal resolution TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite. The Community Multiscale Air Quality (CMAQ) modelling system is used to provide trace gas vertical profiles to recalculate air mass factor for deriving consistent satellite-based VCDs. To demonstrate the ability of satellites to diagnose O3 production regimes, we validate satellite-derived O3 sensitivity regimes, along with the OMI and TROPOMI products applied, using HCHO and NO2 retrievals from the aircraft-based instruments: Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) spectrometer and Geostationary Coastal and Air Pollution Events Airborne Simulator (GCAS) instruments deployed during the LISTOS-2018 and OWLETS-2 campaigns. The results of this work will advance our fundamental knowledge of the spatiotemporal evolution of the complex O3-NOx-VOC chemistry, and will advance our understanding about the quality of remote-sensing, suborbital, and satellite data for accurately observing surface-level O3 production sensitivity.

Satellite↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

X-Band CubeSat Communication System Demonstration

Today's CubeSats mostly operate their communications at UHF- and S-band frequencies. UHF band is presently crowded, thus downlink communications are at lower data rates due to bandwidth limitations and are unreliable due to interference. This research presents an end-to-end robust, innovative, compact, efficient and low cost S-band uplink and X-band downlink CubeSat communication system demonstration between a balloon and a Near Earth Network (NEN) ground system. Since communication systems serve as umbilical cords for space missions, demonstration of this X-band communication system is critical for successfully supporting current and future CubeSat communication needs. This research has three main objectives. The first objective is to design, simulate, and test a CubeSat S- and X-band communication system. Satellite Tool Kit (STK) dynamic link budget calculations and HFSS Simulations and modeling results have been used to trade the merit of various designs for small satellite applications. S- and X-band antennas have been tested in the compact antenna test range at Goddard Space Flight Center (GSFC) to gather radiation pattern data. The second objective is simulate and test a CubeSat compatible X-band communication system at 12.5Mbps including S-band antennas, X-band antennas, Laboratory for Atmospheric and Space Physics (LASP) /GSFC transmitter and an S-band receiver from TRL-5 to TRL-8 by the end of this effort. Different X-band communication system components (antennas, diplexers, etc.) from GSFC, other NASA centers, universities, and private companies have been investigated and traded, and a complete component list for the communication system baseline has been developed by performing analytical and numerical analysis. This objective also includes running simulations and performing trades between different X-band antenna systems to optimize communication system performance. The final objective is to perform an end-to-end X-band CubeSat communication system demonstration between a balloon and/or a sounding rocket and a Near Earth Network (NEN) ground system. This paper presents CubeSat communication systems simulation results, analysis of X-band and S-band antennas and RF front-end components, transceiver design, analysis and optimization of space-to-ground communication performance, subsystem development, as well as the test results for an end-to-end X-band CubeSat communication system demonstration. The outcome of this work will be used to pave the way for next generation NEN-compatible X-band CubeSat communication systems to support higher data rates with more advanced modulation and forward error correction (FEC) coding schemes, and to support and attract new science missions at lower cost. It also includes an abbreviated concept of operations for CubeSat users to utilize the NEN, starting from first contact with NASA's communication network and continuing through on-orbit operations.

NEN↗

Enhancing GPM Passive and Combined Microwave Algorithms with Dynamic Surface Information for Drizzle Retrieval and Improved Precipitation Detection Over Land

Following the 2014 launch of the Global Precipitation Measurement Mission (GPM), an unprecedented combination of coincident active and passive microwave observations are available for state of the art precipitation retrieval. The GPM Combined Algorithm forms the backbone of this effort, optimizing geophysical variables for agreement with the full suite of multi-spectral information content. These combined retrievals are then utilized, along with a radiative transfer model, as a database applied for retrievals across a constellation of passive microwave radiometers of varying frequencies. By keeping such retrievals related through the transfer standard of the combined algorithm, level 3 products such as the Integrated Multi-satellitE Retrievals for GPM (IMERG) are able to provide consistent global products for users at the higher temporal resolution required for hydrological applications. In initial versions of the combined product, precipitation retrievals are carried out only in the presence of a signal from the active radar. As a result, light precipitation and drizzle below the threshold of DPR sensitivity are not included in any of the products down the chain from the constellation to IMERG. In this work, the effects of enhancing the retrievals with a surface emissivity and non-raining water vapor retrieval using the passive observations are explored. Over both ocean and land, the surface retrieval is used to identify areas with high probability of light precipitation and drizzle which is then quantified using techniques derived from the higher sensitivity CloudSat mission. Results indicate successful inclusion of drizzle in the retrievals that can then be included in the constellation databases, as well as improvement in passive microwave false positive precipitation signals over land in cases where surface scattering was misinterpreted as precipitation signal. The inclusion of the dynamic surface information also creates a more robust, radiometrically consistent retrieval scheme for process studies and hydrologic applications.

Ringerud, Sarah↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Bone Metabolism on ISS Missions

Spaceflight-induced bone loss is associated with increased bone resorption (1, 2), and either unchanged or decreased rates of bone formation. Resistive exercise had been proposed as a countermeasure, and data from bed rest supported this concept (3). An interim resistive exercise device (iRED) was flown for early ISS crews. Unfortunately, the iRED provided no greater bone protection than on missions where only aerobic and muscular endurance exercises were available (4, 5). In 2008, the Advanced Resistive Exercise Device (ARED), a more robust device with much greater resistance capability, (6, 7) was launched to the ISS. Astronauts who had access to ARED, coupled with adequate energy intake and vitamin D status, returned from ISS missions with bone mineral densities virtually unchanged from preflight (7). Bone biochemical markers showed that while the resistive exercise and adequate energy consumption did not mitigate the increased bone resorption, bone formation was increased (7, 8). The typical drop in circulating parathyroid hormone did not occur in ARED crewmembers. In 2014, an updated look at the densitometry data was published. This study confirmed the initial findings with a much larger set of data. In 42 astronauts (33 male, 9 female), the bone mineral density response to flight was the same for men and women (9), and those with access to the ARED did not have the typical decrease in bone mineral density that was observed in early ISS crewmembers with access to the iRED (Figure 1) (7). Biochemical markers of bone formation and resorption responded similarly in men and women. These data are encouraging, and represent the first in-flight evidence in the history of human space flight that diet and exercise can maintain bone mineral density on long-duration missions. However, the maintenance of bone mineral density through bone remodeling, that is, increases in both resorption and formation, may yield a bone with strength characteristics different from those that existed before space flight. Studies to assess bone strength after flight are underway at NASA, to better understand the results of bone remodeling. Studies are also underway to evaluate optimized exercise protocols and nutritional countermeasures. Regardless, there is clear evidence of progress being made to protect bone during spaceflight.

Smith, S. M.↗

Ground Vibration Testing Options for Space Launch Vehicles

New NASA launch vehicles will require development of robust systems in a fiscally-constrained environment. NASA, Department of Defense (DoD), and commercial space companies routinely conduct ground vibration tests as an essential part of math model validation and launch vehicle certification. Although ground vibration testing must be a part of the integrated test planning process, more affordable approaches must also be considered. A study evaluated several ground vibration test options for the NASA Constellation Program flight test vehicles, Orion-1 and Orion-2, which concluded that more affordable ground vibration test options are available. The motivation for ground vibration testing is supported by historical examples from NASA and DoD. The approach used in the present study employed surveys of ground vibration test subject-matter experts that provided data to qualitatively rank six test options. Twenty-five experts from NASA, DoD, and industry provided scoring and comments for this study. The current study determined that both element-level modal tests and integrated vehicle modal tests have technical merits. Both have been successful in validating structural dynamic math models of launch vehicles. However, element-level testing has less overall cost and schedule risk as compared to integrated vehicle testing. Future NASA launch vehicle development programs should anticipate that some structural dynamics testing will be necessary. Analysis alone will be inadequate to certify a crew-capable launch vehicle. At a minimum, component and element structural dynamic tests are recommended for new vehicle elements. Three viable structural dynamic test options were identified. Modal testing of the new vehicle elements and an integrated vehicle test on the mobile launcher provided the optimal trade between technical, cost, and schedule.

Patterson, Alan↗

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

Investigating the Simulink Auto-Coding Process

Model based program design is the most clear and direct way to develop algorithms and programs for interfacing with hardware. While coding "by hand" results in a more tailored product, the ever-growing size and complexity of modern-day applications can cause the project work load to quickly become unreasonable for one programmer. This has generally been addressed by splitting the product into separate modules to allow multiple developers to work in parallel on the same project, however this introduces new potentials for errors in the process. The fluidity, reliability and robustness of the code relies on the abilities of the programmers to communicate their methods to one another; furthermore, multiple programmers invites multiple potentially differing coding styles into the same product, which can cause a loss of readability or even module incompatibility. Fortunately, Mathworks has implemented an auto-coding feature that allows programmers to design their algorithms through the use of models and diagrams in the graphical programming environment Simulink, allowing the designer to visually determine what the hardware is to do. From here, the auto-coding feature handles converting the project into another programming language. This type of approach allows the designer to clearly see how the software will be directing the hardware without the need to try and interpret large amounts of code. In addition, it speeds up the programming process, minimizing the amount of man-hours spent on a single project, thus reducing the chance of human error as well as project turnover time. One such project that has benefited from the auto-coding procedure is Ramses, a portion of the GNC flight software on-board Orion that has been implemented primarily in Simulink. Currently, however, auto-coding Ramses into C++ requires 5 hours of code generation time. This causes issues if the tool ever needs to be debugged, as this code generation will need to occur with each edit to any part of the program; additionally, this is lost time that could be spent testing and analyzing the code. This is one of the more prominent issues with the auto-coding process, and while much information is available with regard to optimizing Simulink designs to produce efficient and reliable C++ code, not much research has been made public on how to reduce the code generation time. It is of interest to develop some insight as to what causes code generation times to be so significant, and determine if there are architecture guidelines or a desirable auto-coding configuration set to assist in streamlining this step of the design process for particular applications. To address the issue at hand, the Simulink coder was studied at a foundational level. For each different component type made available by the software, the features, auto-code generation time, and the format of the generated code were analyzed and documented. Tools were developed and documented to expedite these studies, particularly in the area of automating sequential builds to ensure accurate data was obtained. Next, the Ramses model was examined in an attempt to determine the composition and the types of technologies used in the model. This enabled the development of a model that uses similar technologies, but takes a fraction of the time to auto-code to reduce the turnaround time for experimentation. Lastly, the model was used to run a wide array of experiments and collect data to obtain knowledge about where to search for bottlenecks in the Ramses model. The resulting contributions of the overall effort consist of an experimental model for further investigation into the subject, as well as several automation tools to assist in analyzing the model, and a reference document offering insight to the auto-coding process, including documentation of the tools used in the model analysis, data illustrating some potential problem areas in the auto-coding process, and recommendations on areas or practices in the current Ramses model that should be further investigated. Several skills were required to be built up over the course of the internship project. First and foremost, my Simulink skills have improved drastically, as much of my experience had been modeling electronic circuits as opposed to software models. Furthermore, I am now comfortable working with the Simulink Auto-coder, a tool I had never used until this summer; this tool also tested my critical thinking and C++ knowledge as I had to interpret the C++ code it was generating and attempt to understand how the Simulink model affected the generated code. I had come into the internship with a solid understanding of Matlab code, but had done very little in using it to automate tasks, particularly Simulink tasks; along the same lines, I had rarely used shell script to automate and interface with programs, which I gained a fair amount of experience with this summer, including how to use regular expression. Lastly, soft-skills are an area everyone can continuously improve on; having never worked with NASA engineers, which to me seem to be a completely different breed than what I am used to (commercial electronic engineers), I learned to utilize the wealth of knowledge present at JSC. I wish I had come into the internship knowing exactly how helpful everyone in my branch would be, as I would have picked up on this sooner. I hope that having gained such a strong foundation in Simulink over this summer will open the opportunity to return to work on this project, or potentially other opportunities within the division. The idea of leaving a project I devoted ten weeks to is a hard one to cope with, so having the chance to pick up where I left off sounds appealing; alternatively, I am interested to see if there are any opening in the future that would allow me to work on a project that is more in-line with my research in estimation algorithms. Regardless, this summer has been a milestone in my professional career, and I hope this has started a long-term relationship between JSC and myself. I really enjoy the thought of building on my experience here over future summers while I work to complete my PhD at Missouri University of Science and Technology.

Gualdoni, Matthew J.↗

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

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

thermal management↗

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

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

nuclear electric propulsion↗

Revisiting the Effectiveness Of HCHO/NO2 Ratios for Inferring Ozone Sensitivity to Its Precursors Using High Resolution Airborne Remote Sensing Observations in A High Ozone Episode During the KORUS-AQ Campaign

The nonlinear chemical processes involved in ozone production (P(O3)) have necessitated using proxy indicators to convey information about the primary dependence of P(O3) on volatile organic compounds (VOCs) or nitrogen oxides (NOx). In particular, the ratio of remotely sensed columns of formaldehyde (HCHO) to nitrogen dioxide (NO2) has been widely used for studying O3 sensitivity. Previous studies found that the errors in retrievals and the incoherent relationship between the column and the near-surface concentrations are a barrier in applying the ratio in a robust way. In addition to these obstacles, we provide calculational-observational evidence, using an ensemble of 0-D photochemical box models constrained by DC-8 aircraft measurements on an ozone event during the Korea-United States Air Quality (KORUS-AQ) campaign over Seoul, to demonstrate the chemical feedback of NO2 on the formation of HCHO is a controlling factor for the transition line between NOx-sensitive and NOx-saturated regimes. A fixed value (~2.7) of the ratio of the chemical loss of NOx (LNOx) to the chemical loss of HO2+RO2 (LROx) perceptibly differentiates the regimes. Following this value, data points with a ratio of HCHO/NO2 less than 1 can be safely classified as NOx-saturated regime, whereas points with ratios between 1 and 4 fall into one or the other regime. We attribute this mainly to the HCHO-NO2 chemical relationship causing the transition line to occur at larger (smaller) HCHO/NO2 ratios in VOC-rich (VOC-poor) environments. We then redefine the transition line to LNOx/LROx~2.7 that accounts for the HCHO-NO2 chemical relationship leading to HCHO = 3.7 × (NO2 – 1.14 × 1016 molec.cm-2). Although the revised formula is locally calibrated (i.e., requires for readjustment for other regions), its mathematical format removes the need for having a wide range of thresholds used in HCHO/NO2 ratios that is a result of the chemical feedback. Therefore, to be able to properly take the chemical feedback into consideration, the use of HCHO = a × (NO2 – b) formula should be preferred to the ratio in future works. We then use the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument to study O3 sensitivity in Seoul. The unprecedented spatial (250 × 250 m2) and temporal (~every 2 h) resolutions of HCHO and NO2 observations form the sensor enhance our understanding of P(O3) in Seoul; rather than providing a crude label for the entire city, more in-depth variabilities in chemical regimes are observed that should be able to inform mitigation strategies correspondingly.

KORUS-AQ; formaldehyde; NOx; VOC; OMI; GeoTASO; oz↗

A Prototype Early Planetary Organic Processor Assembly (OPA) Based on Dual-Stage Anaerobic Membrane Bioreactor (AnMBR) for Fecal and Food Waste Treatment and Resource Recovery

Long-duration, deep-space exploration and habitation missions demand robust and reliable technologies to ensure crew health, safety, and mission success. Local food production will be essential for crew nutrition and morale. However, at $10,000/lb, the payload costs and mass/volume limitations to transport and provide the necessary resources, including fertilizer, for an anticipated 30-month mission become challenging over time. For mission success and sustainability, the Environmental Control and Life Support System (ECLSS) of the near future will need to recover resources from all “waste” sources and be near-closed loop. Organic wastes (e.g., fecal and food) offer a renewable source of C, N, P, water and other trace elements necessary to sustain crop production. However, these high solidwastes are often difficult to treatdue to factors includingheterogeneity, complexity, high organic strength, and the presence of pathogens. To date, there is no flight-ready technology capable of treating mixed organic wastes, creating a technology gap for future space missions. To address this need, a prototype Organic Processor Assembly (OPA) was developed through collaboration between the University of South Florida (USF) and NASA’s Kennedy Space Center (KSC). The OPA is based on the anaerobic membrane bioreactor (AnMBR), a hybrid technology coupling high-rate anaerobic digestion with membrane filtration. The system is designed for an early planetary base (EPB) scenario to aid in closing the resource recovery loop and decreaseresupply dependence. This presentation discusses initial research pertaining to: 1) design challenges in maximizing hydraulic/organic throughputandReliability, Availability, Maintainability, and Safety(RAMS)while minimizingmass and volume; 2) create capabilities for treating simulated high solids waste under steady and non-steady state conditions; and 3) measuresolids performance parameter(s). Future research and development pertaining to further optimization on system operation, performance, and expanded treatment capabilities arepresented.

Bioreactor↗

Assessing and Advancing the Potential of Quantum Computing: A NASA Case Study

Quantum computing is one of the most enticing computational paradigms with the potential to revolutionize diverse areas of future-generation computational systems. While quantum computing hardware has advanced rapidly, from tiny laboratory experiments to quantum chips that can outperform even the largest supercomputers on specialized computational tasks, these noisy- intermediate scale quantum (NISQ) processors are still too small and non-robust to be directly useful for any real-world applications. In this paper, we describe NASA’s work in assessing and advancing the potential of quantum computing. We discuss advances in algorithms, both near- and longer-term, and the results of our explorations on current hardware as well as with simulations, including illustrating the benefits of algorithm-hardware codesign in the NISQ era. This work also includes physics-inspired classical algorithms that can be used at application scale today. We discuss innovative tools supporting the assessment and advancement of quantum computing, and describe improved methods for simulating quantum systems of various types on high performance computing systems that incorporate realistic error models. We provide an overview of recent methods for benchmarking, evaluating, and characterizing quantum hardware for error mitigation, computational purposes.

quantum computing↗

Clinical Decision Support Project

As NASA plans for exploration missions into deep space, significant challenges are realized due to the distance from Earth. Beside the effects of microgravity and radiation exposure, the astronauts face the additional constraints of isolation, lack of resupply, increasingly difficult evacuation, and delayed and disrupted communication with ground-based medical care providers. These constraints require a paradigm shift from current medical care where crews rely on the real-time communications with ground-based medical care providers toward Earth-independent medical operations for astronaut medical care. Medical expertise and decision-making are ground-based for current International Space Station and planned Lunar missions. However, a deep space exploration crew will need to autonomously perform the detection, diagnosis, treatment, and prevention of medical conditions. One approach to provide Earth-independent medical operations is to augment the requisite knowledge, skills, and abilities (KSAs) of a time-constrained crew—operating under stressful conditions, combatting fatigue, and facing a potential medical crisis—with a robust clinical decision support system (CDSS). The CDSS is envisioned as an integrated, software-based tool deployed on a laptop computer or handheld device. The CDSS will assist the crew and ground support when interacting with knowledge/data bases (e.g. records, pharmacy, schedule), instrumentation (e.g. imaging, physiological monitoring devices), and habitat (e.g. wellness system, task performance system) and vehicle systems (e.g. environmental system, communication system). In addition, the human interface will employ a context-based approach that accounts for the crew’s situation. Thus, extraneous and clinically/operationally non-relevant information are reduced to avoid an increase in cognitive load. The framework of an ideal spaceflight CDSS is to include core and advanced analytical features that maintain a flexible platform for integrating new technology in the future. The Exploration Medical Capability (ExMC) Element of the Human Research Program (HRP) is expanding the boundaries of space medical systems to advance the care of astronauts on future exploration missions beyond low Earth orbit by actively identifying and testing next-generation medical care and crew health maintenance technologies. The Clinical Decision Support (CDS) project addressed ap Medical-701 within the Inflight Medical Conditions risk: “We need to increase inflight medical capabilities and identify new capabilities that (a) maximize benefit and/or (b) reduce “costs” on human system/mission/vehicle resources.” Though mass, volume, and power will face increasing constraints, the projected computational capabilities of spacecraft systems will increase exponentially as information technology advances in this decade and beyond. Hence, data, software, and computational resources will play an essential and synergistic role in maintaining crew health, wellness, and performance in deep space missions. The focus of the CDS project was to develop recommended requirements for an in-vehicle CDSS that acts as a ‘virtual assistant’ for delivering optimal health, performance, and medical care during exploration missions. In fiscal year 2022 (FY22), the CDS project was chartered to baseline and/or revise all CDS project related documentation and update the CDS project model to include the revised CDSS Concept of Operations, revised systems-based modeling language (SysML) activity diagrams, and baseline requirements. The focus of this presentation will be an overview of the CDS products and CDS model content.

Decision Support↗

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi ]Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post ]landing segment of the Flight Software (FSW). The vehicle crew is the prime input for touchdown detection, followed by an autonomous FSW algorithm, and finally a strictly time based backup timer. RCS thrusters must be inhibited before submersion in water to protect against possible damage due to firing these jets under water. In addition, neglecting to declare touchdown will not allow the vehicle to transition to post ]landing activities such as activating the Crew Module Up ]righting System (CMUS), resulting in possible loss of communication and difficult recovery. A previous AIAA paper gAssessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module h concluded that a strictly Inertial Measurement Unit (IMU) based detection method using an acceleration spike algorithm had the highest safety margins and shortest detection times of other methods considered. That study utilized finite element simulations of vehicle splashdown, generated by LS ]DYNA, which were expanded to a larger set of results using a Kriging surface fit. The study also used the Decelerator Systems Simulation (DSS) to generate flight dynamics during vehicle descent under parachutes. Proto ]type IMU and FSW MATLAB models provided the basis for initial algorithm development and testing. This paper documents an in ]depth trade study, using the same dynamics data and MATLAB simulations as the earlier work, to further develop the acceleration detection method. By studying the combined effects of data rate, filtering on the rotational acceleration correction, data persistence limits and values of acceleration thresholds, an optimal configuration was determined. The lever arm calculation, which removes the centripetal acceleration caused by vehicle rotation, requires that the vehicle angular acceleration be derived from vehicle body rates, necessitating the addition of a 2nd order filter to smooth the data. It was determined that using 200 Hz data directly from the vehicle IMU outperforms the 40 Hz FSW data rate. Data persistence counter values and acceleration thresholds were balanced in order to meet desired safety and performance. The algorithm proved to exhibit ample safety margin against early detection while under parachutes, and adequate performance upon vehicle splashdown. Fall times from algorithm initiation were also studied, and a backup timer length was chosen to provide a large safety margin, yet still trigger detection before CMUS inflation. This timer serves as a backup to the primary acceleration detection method. Additionally, these parameters were tested for safety on actual flight test data, demonstrating expected safety margins.

Daum, Jared↗