Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complexity bounding”

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.

446 records · Page 25

A Proposed Clock Synchronization Method for the Solar System Internet

Networked communications in space are necessary to achieve scalability in terms of the number of communicating nodes but also in terms of the overall system complexity. A key component to such a system is the ability to synchronize clocks, which is the focus of this paper. The so-called Solar System Internet (SSI) will be built upon Delay Tolerant Networking (DTN), which, in analogy to the Internet Protocol (IP), can be considered a suite of protocols necessary for networking in the space domain. Therefore, our goal is to extend this suite to include a DTN clock synchronization capability, analogous to the Network Time Protocol (NTP) used in the Internet. A motivating example of a network in space is NASA’s LunaNet, a vision for a multi-hop multi-path network extending to the moon wherein not all nodes will have direct connections to an authoritative reference clock. In this paper, we propose a general clock synchronization methodology and algorithm that could be used for LunaNet as well as more elaborate time-varying networks. In recent years, DTN has benefited from modeling efforts founded on the mathematical tool of sheaves. Here we continue this work to provide an approach to clock synchronization. Due to the time-varying nature of space networks, absolute consensus is not possible. However, the sheaf Laplacian provides a practical, distributed approach to approximating consensus by allowing data to diffuse through the network. In particular, the sheaf Laplacian is readily computable, lending our approach to implementation. Our approach is well suited to handle the difficulties of space networks. For instance, differences in clock accuracy mean certain nodes are more authoritative than others; we can account for these differences through hierarchies in the network, generalizing the strata in NTP. Furthermore, just as error estimation is an integral part of NTP, we are able to give concrete error bounds for our approach. Indeed, different applications (e.g., communications schedules, pointing, navigation, distributed science) will have different requirements, hence it is necessary to maintain clocks within a given tolerance. We outline some of the necessary steps to turn our approach into a practical network protocol that could be used in DTN, and we conclude the paper with suggestions for future research.

Michael Moy↗

Climate Change Impact Uncertainties for Maize in Panama: Farm Information, Climate Projections, and Yield Sensitivities

We present results from a pilot project to characterize and bound multi-disciplinary uncertainties around the assessment of maize (Zea mays) production impacts using the CERES-Maize crop model in a climate-sensitive region with a variety of farming systems (Panama). Segunda coa (autumn) maize yield in Panama currently suffers occasionally from high water stress at the end of the growing season, however under future climate conditions warmer temperatures accelerate crop maturation and elevated CO (sub 2) concentrations improve water retention. This combination reduces end-of-season water stresses and eventually leads to small mean yield gains according to median projections, although accelerated maturation reduces yields in seasons with low water stresses. Calibrations of cultivar traits, soil profile, and fertilizer amounts are most important for representing baseline yields, however sensitivity to all management factors is reduced in an assessment of future yield changes (most dramatically for fertilizers), suggesting that yield changes may be more generalizable than absolute yields. Uncertainty around General Circulation Model (GCM)s' projected changes in rainfall gain in importance throughout the century, with yield changes strongly correlated with growing season rainfall totals. Climate changes are expected to be obscured by the large inter-annual variations in Panamanian climate that will continue to be the dominant influence on seasonal maize yield into the coming decades. The relatively high (A2) and low (B1) emissions scenarios show little difference in their impact on future maize yields until the end of the century. Uncertainties related to the sensitivity of CERES-Maize to carbon dioxide concentrations have a substantial influence on projected changes, and remain a significant obstacle to climate change impacts assessment. Finally, an investigation into the potential of simple statistical yield emulators based upon key climate variables characterizes the important uncertainties behind the selection of climate change metrics and their performance against more complex process-based crop model simulations, revealing a danger in relying only on long-term mean quantities for crop impact assessment.

Maize↗

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

New Rover Conops with High-Performance Onboard Computing: Give Up Raw Data to Reduce Ops Cost and Do More Science

A major portion of time during the tactical operation of Mars rovers is spent for selecting, prioritizing, and coordinating sciences and engineering activities such that they fit within resource constraints, including the downlink data volume, energy, and time. In particular, the downlink data volume constraint is getting particularly tighter in recent missions because modern instruments produce increasingly high data volume while the communication bandwidth is essentially bounded by the law of physics. Tactical operation would be substantially simplified, hence the operation cost could be reduced, if the data volume constraint is relaxed or even removed. In this abstract, we propose a new operation paradigm for achieving this goal. The key observation is that, both in science and engineering applications, the bit size of raw data is typically much greater than the volume of processed information that is needed for scientific or engineering analysis. For example, a full-resolution image from Mastcam-Z, the main science camera on Perseverance, is about 700 kB in volume and we downlinked 29,685 images up to Sol 243, totaling ~20 GB of data. But of course, scientists do not use every pixel of these images; what they really look for in the images are geological features, typically represented by specific geometric configurations or textures. An end product after processing hundreds of Mascam-Z images could be a single geological map summarizing the spatial distribution of the features. For another example, a 100-meter drive of Perseverance produces 7-12 MB of drive telemetry, which records every detail of the rover's motion at 8 Hz, including position, attitude, steering angles, encoder readings, motor currents and many other information. But what the ground engineers eventually pay attention to is the signs of anomaly, such as excessive motor currents or high slip; if a drive is nominal, the vast majority of this data is unused. What if, then, we process the raw data onboard and only downlink the processed data that is relevant to scientific or engineering analyses, such as a list of detected science features (with cropped images) or a list of potential signs of anomaly while driving? A major roadblock for such onboard, high-level information processing has been the onboard computational resource. RAD750, the main onboard computer of Perseverance, is obviously not sufficient for performing complex image or signal processing such as object detection, semantic segmentation, or anomaly detection. Interestingly, RAD750 is not the best processor that Perseverance has; Qualcomm's Snapdragon 801, a modern mobile processor, is on her Heli Base Station, a device for communicating with Mars Helicopter Ingenuity; also, Intel's Atom E3845 processors are on engineering cameras. In the reminder of this paper, we will introduce two particular uses cases of these high-performance co-processors (meaning auxiliary CPU, GPU, or other types of processors that are separate from the main processor that runs the main flight software) for lowering operation cost and accommodating more science activities for a given communication constraint.

Didier, A.↗

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗

Soluble Organic Matter (SOM) analysis of the Hayabusa2 samples: The first results

The Hayabusa2 spacecraft successfully collected the surface and possible sub-surface materials of the asteroid 162173 Ryugu. Ryugu is a C-type asteroid characterized by a low-albedo surface probably consisting of hydrous minerals and carbonaceous materials. [1] The direct optical and spectral analysis of the returned samples indicates that Ryugu material is dominated by hydrous carbonaceous chondrite-like matter (similar to CI chondrites) [2]. Since carbonaceous chondrites have generally yielded various types of organic matter, the collected Ryugu grains are expected to contain diverse types of organic compounds including bio-related molecules. The occurrence of organic compounds in the Ryugu surface will provide clues to the evolution of prebiotic molecules and their preservations associated with aqueous alteration of the primitive asteroid. The initial analysis of soluble organic matter (SOM) of the Hayabusa2-returned samples has been performed by an international team consisting of 32 members. Because the sample amount available for comprehensive SOM analyses is limited, and because the SOM is expected to be present as a complex mixture of various types of organic compounds with very small concentrations of each compound, high-sensitivity and high-resolution analytical techniques have been developed using carbonaceous meteorites [e.g. 3]. Two aggregate samples of the Ryugu grains (A106 and C107) were allocated for the solvent extractions and bulk carbon (C), hydrogen (H), nitrogen (N), sulfur (S) and oxygen (O) measurements. The A106 sample was collected during the 1st sampling in February 2019 and the C107 sample was collected during the 2nd sampling in July 2019 after the Small Carry Impactor (SCI) operation. They consist mainly of particles smaller than 1 mm in diameter, and each sample mass was 38-39 mg. They were firstly investigated spectroscopically in the near infrared wavelength range by the Stone Team prior to the solvent extractions. Other small grains (A0080 and C0057) were also allocated for this study to investigate the spatial distribution of organic compounds on the sample surface. The extraction and analytical measurements implemented by the SOM Team are summarized in Figure 1. Each powder sample was extracted sequentially with non-polar to polar solvents, i.e., hexane, dichloromethane (DCM), methanol (MeOH) and H2O, for non-targeted analysis to reveal the compound composition. Each solvent extract was analyzed by solution state nuclear magnetic resonance (NMR) spectroscopy [4], Fourier transform- ion cyclotron resonance/mass spectrometry (FT-ICR/MS) with ESI and APPI ionization [5] and by high-resolution mass spectroscopy using Orbitrap MS coupled with nano-liquid chromatography (nanoLC/Orbitrap MS) [6], and using two dimensional gas chromatography/mass spectrometry (GC×GC/MS). The extracted residues were passed to the Chemistry Team for further inorganic element analysis. The other powder sample was subjected to the hot water extraction for amino acid analyses including chiral isomer separation, which was performed by three-dimensional (3D) high-performance liquid chromatography (HPLC) with high-sensitivity fluorescence detection (FD) [7] and by HPLC/FD coupled with quadrupole-time of flight/mass spectrometry (QToF/MS) [8]. After the hot water extraction, the residue was split into two halves. One half was further extracted with hydrochloric acid (HCl) to analyze for amino acids in bound-form. The other half was sequentially extracted with DCM/MeOH (1/1) to analyze semi-polar compounds such as polycyclic aromatic hydrocarbons (PAHs) by GC/MS, followed by further extraction with formic acid to analyze polar heterocyclic compounds, and subsequent extraction with HCl to detect bound-form polar compounds. The extracted residues were passed to the IOM Team for the analysis of insoluble organic matter (IOM). Compound-specific stable isotope analyses will be performed using GC/combustion/isotope ratio mass spectrometry (GC/C/IRMS) if the compound concentration is high enough to enable such an isotopic measurement. All extraction procedures were performed on an ISO 6 (Class 100) clean bench inside an ISO 5 (Class 1000) clean room. Baked serpentine powder was also analyzed as a procedural blank. In situ organic compound analysis with the molecular imaging was performed using desorption electrospray ionization (DESI) equipped with Orbitrap MS [9, 10], followed by spatial imaging of organic compounds using ToF/secondary ion mass spectrometry (ToF/SIMS) [11]. The bulk chemical and isotopic compositions of CNS and HO were determined using nano-elemental analysis/isotope ratio mass spectrometry (nanoEA/IRMS) [12] and EA/pyrolysis/IRMS, respectively. We have identified a variety of indigenous organic compounds in the extracts of both A106 and C107 samples. The Ryugu grains host organic molecules under the high-vacuum and cosmic-ray irradiation environment of the asteroid surface. The analysis of extracted molecules is in progress, and the first results will be presented at the symposium.

Hiroshi Naraoka↗

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

The Utilization Profiles of the CCSDS Unified Space Link Protocol (USLP)

The purpose of this paper is to identify the utilization profiles for interfacing the Data Protocol Sublayer using the Unified Space Link Protocols (USLP) (reference 1) with the space link coding procedures as specified in the CCSDS Coding & Synchronization Blue Books (references 2 through 5), used in both telecommand and telemetry applications. This paper describes how the USLP Protocol utilizes the coding and synchronization sublayer to support: a. Direct to Earth (DTE) telemetry links for engineering and science data b. Direct to Earth (DTE) telemetry links for very high rate science data c. Direct from Earth (DFE) command, sequencing and flight software loads d. Space to Space Links (Proximity) utilized by orbiters for data exchange to/from surface bound assets. The CCSDS has divided the functions of the Data Link Layer into two sublayers: the Data Link Protocol Sublayer (DLP-SL) and the Coding and Synchronization Sublayer (CS-SL). The Data Link Protocol Sublayer (DLP-SL) interfaces to the users, accepting the data that is to be transported, on the sending side of the link, and delivering that data on the receiving end. The Transfer Frame is the data unit that is transferred across the Data Link Protocol Sublayer and the Coding and Synchronization Sublayer boundary. The Coding and Synchronization Sublayer (CS-SL) provides the encoding, randomization, and frame synchronization functions that prepares the USLP Transfer Frame for transport across the space link. The CS-SL is divided into 2 processes: 1) The Frame Interface Processes (FIP) performs the interface functions required to prepare the data for delivery to the Coding/Decoding Process (CDP). This process includes prepending a Frame Start Marker to the provided frame, when management has designated that the frame is not to be aligned to the codeblock or when there is no block code used. 2) The Coding/Decoding Process (CDP) performs the forward error correction processes that are used to optimize the performance of the link and minimize the error rate. The CDP creates the symbol stream that is delivered to the Physical Layer. The transfer of the USLP transfer frames across different types of space links is the focus of this paper. The Protocol Data Unit (PDU) that is passed in both directions between the Data Link Protocol Sublayer (DLP-SL) and Coding and Synchronization Sublayer (CS-SL) is the transfer frame. The USLP frame structure provides flexibility that can be constrained by the functions utilized within the CS-SL that prepare the transfer frame for transit. For example, the USLP transfer frame contains a length field that enables the frame to be of variable length but CS-SL under certain conditions may constrain the frame to be fixed in length. This paper describes 5 operational modes available for use by the Data Link Layer to provide data exchange across the USLP space link. These modes are different because different operational requirements apply to vastly different types of space links and thus the communications implementation requirements differ. The environmental issues include the power or energy available, the distance between the end points of the link, the complexity of the equipment available at those end points, the atmospheric conditions and radiometric frequency selection. The CS-SL utilizes different forward error correcting codes supported by specific operational modes to configure the data for transit. This paper describes all of the operational modes in a series of data models which decompose the functionality between the Data Link Protocol Sublayer and the Coding and Synchronization sublayer. The operational modes described are: 1. Uncoded Mode: has been used for short links that contain significant available power to provide an acceptable frame error rate. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 2. Convolutional Only Mode: is currently the prime forward error correction coding used for the proximity links. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 3. Variable Length Frame Aligned to Variable Length Codeblock (TC): is used for Direct from Earth links were power levels are high and the simple, least complex code i.e., the BCH code is used. This mode has been in use since the early 1970s. The BCH code is a short code and the decoder is easy to implement. 4. Fixed Length Frame Aligned to Fixed Length Codeblock (AOS/TM): was introduced when the concatenated Convolutional and Reed-Solomon Code was formulated to provide significant reduction in link data error rate and the ability to determine if there was an error in the decoded codeblock. The frame is aligned to the codeblock so that there is a one to one relationship of frame errors to codeblock errors without additional error detection coding being added. This mode requires the protocol frames to be the exact size of the message portion of the codeblock. 5. Frames Unaligned to Fixed Length Codeblocks (Currently used for very high rates and space to space links): This mode is currently used for missions that have a very high data rate that can be controlled adaptively as the environment changes and as the next generation operating mode for the proximity link. This mode from a coded data stream point of view is exactly like that described in 4. above, except that the frame need not be aligned to the codeblock. There is no requirement on frame length when using this mode. Thus when using USLP it can be used to support links that require short or long frames. There is also no mandatory requirement that frames cannot be separated by idle data reducing the tight data rate connection requirements between the data link protocol sublayer and the coding & synchronization sublayer. In conclusion, how these operational modes can be put to use in mission operational scenarios is described for Direct from Earth links (DFE), Direct to Earth links (DTE), and Proximity links.

Greenberg, E.↗

A General Closed-Form Solution for the Lunar Reconnaissance Orbiter (LRO) Antenna Pointing System

The National Aeronautics and Space Administration s (NASA) Lunar Reconnaissance Orbiter (LRO) launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle into a direct insertion trajectory to the Moon LRO, designed, built, and operated by the NASA Goddard Space Flight Center in Greenbelt, MD, is gathering crucial data on the lunar environment that will help astronauts prepare for long-duration lunar expeditions. During the mission s nominal life of one year its six instruments and one technology demonstrator will find safe landing site, locate potential resources, characterize the radiation environment and test new technology. To date, LRO has been operating well within the bounds of its requirements and has been collecting excellent science data images taken from the LRO Camera Narrow Angle Camera (LROC NAC) of the Apollo landing sites have appeared on cable news networks. A significant amount of information on LRO s science instruments is provided at the LRO mission webpage. LRO s Attitude Control System (ACS), in addition to controlling the orientation of the spacecraft is also responsible for pointing the High Gain Antenna (HGA). A dual-axis (or double-gimbaled) antenna, deployed on a meter-long boom, is required to point at a selected Earth ground station. Due to signal loss over the distance from the Moon to Earth, pointing precision for the antenna system is very tight. Since the HGA has to be deployed in spaceflight, its exact geometry relative to the spacecraft body is uncertain. In addition, thermal distortions and mechanical errors/tolerances must be characterized and removed to realize the greatest gain from the antenna system. These reasons necessitate the need for an in-flight calibration. Once in orbit around the moon, a series of attitude maneuvers was conducted to provide data needed to determine optimal parameters to load onboard, which would account for the environmental and mechanical errors at any antenna orientation. The nominal geometry for the HGA involves an outer gimbal axis that is exactly perpendicular to the inner gimbal axis, and a target direction that is exactly perpendicular to the outer gimbal axis. For this nominal geometry, closed-form solutions of the desired gimbal angles are simple to get for a desired target direction specified in the spacecraft body fame. If the gimbal axes and the antenna boresight are slightly misaligned, the nominal closed-form solution is not sufficiently accurate for computing the gimbal angles needed to point at a target. In this situation, either a general closed-form solution has to be developed for a mechanism with general geometries, or a correction scheme has to be applied to the nominal closed-form solutions. The latter has been adopted for Solar Dynamics Observatory (SDO) as can be seen in Reference 1, and the former has been used for LRO. The advantage of the general closed-form solution is the use of a small number of parameters for the correction of nominal solutions, especially in the regions near singularities. Singularities here refer to cases when the nominal closed-form solutions have two or more solutions. Algorithm complexity, however, is the disadvantage of the general closed-form solution.

Shah, Neerav↗

Exercise Equipment Usability Assessment for a Deep Space Concept Vehicle

With international aspirations to send astronauts to deep space, the world is now faced with the complex problem of keeping astronauts healthy in unexplored hostile environments for durations of time never before attempted by humans. The great physical demands imparted by space exploration compound the problem of astronaut health, as the astronauts must not only be healthy, but physically fit upon destination arrival in order to perform the scientific tasks required of them. Additionally, future deep space exploration necessitates the development of environments conducive to long-duration habitation that would supplement propulsive vehicles. Space Launch System (SLS) core stage barrel sections present large volumes of robust structure that can be recycled and used for long duration habitation. This assessment will focus on one such conceptual craft, referred to as the SLS Derived Habitat (SLS-DH). Marshall Space Flight Center's (MSFC) Advanced Concepts Office (ACO) has formulated a high-level layout of this SLS-DH with parameters such as floor number and orientation, floor designations, grid dimensions, wall placement, etc. Yet to be determined, however, is the layout of the exercise area. Currently the SLS-DH features three floors laid out longitudinally, leaving 2m of height between the floor and ceilings. This short distance between levels introduces challenges for proper placement of exercise equipment such as treadmills and stationary bicycles, as the dynamic envelope for the 95th percentile male astronauts is greater than 2m. This study aims to assess the optimal equipment layout and sizing for the exercise area of this habitat. Figure 1 illustrates the layout of the DSH concept demonstrator located at MSFC. The exercise area is located on the lower level, seen here as the front half of the level occupied by a crew member. This small volume does not allow for numerous or bulky exercise machines, so the conceptual equipment has been limited to a treadmill and stationary bicycle. With the most current treadmill aboard the International Space Station (ISS), the Combined Operational Load-Bearing External Resistance Treadmill (COLBERT), being located in an International Standard Payload Rack (ISPR), the bottom of the conceptual treadmill features a height of 38in. Making the treadmill flush with the floor would be impossible in this rack configuration, as the distance from the outer wall of the spacecraft to the bottom floor would be too shallow. From preliminary sizing, the 38in required for the bottom of the treadmill combined with a 78in operational envelope for a 95th percentile may not be accommodated in the exercise area in a vertical orientation. Figure 2 demonstrates the volume required (in maroon) for an ISPR-bound treadmill in the concept demonstrator. Early indications as seen in this figure indicate that the crew members would contact the ceiling in such an arrangement. An assessment will be conducted to evaluate various orientations of exercise equipment in the concept demonstrator. Orientations to be tested include putting the bottom of the treadmill on the wall, having the treadmill at an angle in the floor both horizontally and vertically, and having a shorter (non-rack bound) treadmill in a vertical orientation on the floor. This assessment will yield findings regarding sizing of the area and how well participants feel they could exercise in such an environment. Due to the restrictions of assessing a microgravity vehicle in a normal-gravity environment, simulations in MSFC's Virtual Environments Lab (VEL) may be necessary. Final deliverables will include recommendations regarding the location and size of possible exercise equipment aboard the SLS-Derived DSH.

Rhodes, Brooke M.↗

Cloud Computing Methods for Near Rectilinear Halo Orbit Trajectory Design

Complicated mission design problems require innovative computational solutions. As spacecraft depart from a proposed Gateway in a Near Rectilinear Halo Orbit (NRHO), recontact analysis is required to avoid risk of collision and ensure safe operations. Escape dynamics from NRHOs are governed by multiple gravitational bodies, yielding a trajectory design space that is exhaustively large. This paper summarizes the recontact analysis for departure from the NRHO and describes how the Deep Space Trajectory Explorer (DSTE) trajectory design software incorporates high performance cloud computing to compute and visualize the orbit design space. Recent focus on exploration missions to cislunar space has kindled accelerated interest in multibody orbit solutions. Trajectory analysis in the presence of multiple gravity fields is complex, and innovative computational tools are needed to simplify complicated design spaces, to generate large quantities of data quickly, and to visualize the output for user accessibility. The Gateway mission is a prime example. The Gateway1 is proposed as a human outpost in deep space. The current baseline orbit for the Gateway is a Near Rectilinear Halo Orbit (NRHO) near the Moon.2 The NRHO exists in a regime that experiences the gravitational effects of the Earth and the Moon simultaneously, complicating orbit analysis. The mission design process benefits greatly from updated computational tools for multibody missions like the Gateway. As an example, consider the problem of assessing the risk of collision in an NRHO. As a staging location to missions to the lunar surface and beyond the Earth-Moon system, the Gateway will experience spacecraft and other objects regularly arriving and departing. Departing objects potentially include spent logistics modules, visiting crew vehicles, debris objects, wastewater particles, and cubesats. Each departure is governed by the dynamics of the Gateway orbit and the surrounding dynamical environment. Over time, any unmaintained object in such an orbit eventually departs due to the small instabilities associated with the NRHOs. A separation maneuver speeds the departure from the NRHO, but the effects of the maneuver on the spacecraft behavior depend on the location, magnitude, and direction of the burn. Escape dynamics from the NRHO with regard to these maneuver options open up an enormous potential trajectory design space where subtle changes in input can produce dramatically large changes in the results. Any departing object must avoid recontacting the Gateway as it leaves the lunar vicinity, and a recontact analysis thus involves a significant number of computations and extensive output data. To explore the dynamics of this extensive design space, the Deep Space Trajectory Explorer3 (DSTE) trajectory design software incorporates new High Performance Computing (HPC) services and novel interactive visualizations. This paper details the HPC and cloud infrastructure techniques that are implemented in the DSTE, applying the new capabilities to analysis of recontact risk with the Gateway in NRHO. NEAR RECTILINEAR HALO ORBITS The Gateway is planned to fly in a lunar NRHO as its baseline orbit. The NRHO families of orbits are subsets of the larger halo families, which originate from planar orbits near the L1 and L2 libration points; the Earth-Moon L2 halo family appears in Figure 1. Each halo orbit is perfectly periodic in the Circular Restricted 3-Body Problem (CR3BP) and becomes a quasi-periodic orbit in a higher fidelity ephemeris force model. The NRHOs are defined as those members of the halo family with bounded stability properties;2 they pass near the Moon at perilune and are nearly polar. Families exist with apolunes located both above the lunar north pole and above the lunar south pole; the Gateway is planned to reside in a southern L2 NRHO in a 9:2 resonance with the lunar synodic period. The 9:2 NRHO is characterized by a period of about 6.5 days, a perilune radius of about 3,500 km, and an apolune radius of about 71,000 km; it is strongly affected by the gravity of both the Earth and the Moon simultaneously. This NRHO offers extended communications with assets on the south pole of the Moon,4 as well as low-cost orbit maintenance and attitude control,5 favorable eclipse avoidance properties,6 and inexpensive transfers from Earth and to other destinations.5,7 The NRHO portion of the southern L2 halo family is highlighted in black in Figure 1, and the 9:2 NRHO appears in blue.

Phillips, Sean M.↗

Effect of salinity on the preservation of biomarkers in hypersaline microbial mat kerogens

Kerogen is the residue that remains after the minerals and organic compounds have been extracted from a rock using acid and organic solvents, respectively [1]. It is not only the largest organic carbon sink in the geological record but the best-preserved fraction due to its recalcitrance and ability to retain syngenetic source information during deposition [2]. Due to the complex, macromolecular, cross-linked polymeric structure of kerogen, it is difficult to analyse intact, often requiring breakdown into smaller subunits through thermal (pyrolysis) or chemical treatment (chemolysis) [3]. Saline and hypersaline environments host diverse microbial communities, which can trap and bind organic matrices such as cellular components (e.g., extracellular polymeric substances, cell wall, peptidoglycan, sheath material, etc.) with sediments, forming layered, laminated structures called microbial mats. Here we present depth profiles of two hypersaline cores through a microbial mat with differing total salinities. We looked at the composition of organic compounds released from kerogen during high temperature, hydrogen-assisted, pyrolysis to understand the effect of salinity on the composition and timing of their incorporation [4]. Field Site: The Exportadora del Sal S. A. (ESSA) saltern of Guerrero Negro (Baja California Sur, Mexico) lies within a typical sabhka environment adjacent to Laguna Ojo de Liebre. Note the spatial extent of the 13 major ponds of the ESSA system. The GN ponds thrive under well-monitored and stable conditions with steady state accretion and degradation of mat layers (around 0.5–1.4 cm per year). Organic mats typically attain up to 10 cm in thickness, depending on the location, which represents about 60 years of growth since these ponds were established. Sedimentation rate is rapid in these non-lithifying and non-mineralising mats, 1–1.5 cm per year. This environment is ideal for investigating in situ lipid preservation and diagenesis without the influence of constantly changing physical parameters. Methods: P4n5 and P5AB cores were collected as ~9 cm and 6.5 cm fresh cores, respectively, in June 2001 and September 2010, respectively from the ESSA saltern. The P4n5 core (9-9.1% salinity) was taken from Pond 4 near 5 and sub-divided into 10 layers. The P5AB core (11.2% salinity) was taken from Pond 5A near 5B and sub-divided into 8 layers. Lyophilised microbial mat powders were solvent-extracted using a modified Bligh-Dyer method and analysed by gas chromatography-mass spectrometry (GC-MS). Additionally, a subset of layers—layer 4 and layer 7—which represented active biomass and sediment-processed, respectively, were subjected to mild-acid methanolysis. The respective residues from conventional solvent extraction and acid methanolysis were loaded with a molybdenum sulfide catalyst and placed into a stainless-steel reactor. Samples were run on the hydropyrolysis set up (heated to 500 °C with 13-15 MPa of H2 pressure), extracted, and analysed by GC-MS. Results: Molecular profiles from the two ponds were significantly different and reflect the contributions of different photosynthetic and respiratory microbial communities to preserved organic matter. For example, the higher salinity P5AB core showed evidence of greater archaeal inputs (similarly observed in lab culture experiments [5]) as well as potential kerogen-bound carotenoids/carotenoid rearranged products or fragments. A hydrophobic emulsion formed during the acid methanolysis processing of layer 4 from both P4n5 and P5AB and layer 7 of P5AB. Hydropyrolysis of this hydrophobic residues released a greater abundance of polycyclic lipids relative to the pre-extracted control. Implications and Future Work: Catalytic hydropyrolysis of kerogen can rapidly generate abundant saturated pyrolysate products from the bound biomarker pool without altering the structures or stereochemistries of the products [6]. Chemical processing of solvent-extractable residues indicated that cellular matrices such as extracellular polymeric substances (EPS) may play a role in the sequestration of polycyclic lipid biomarkers such as steranes and hopanes. Their higher relative abundance compared to the pre-extracted control in the higher salinity layers and cores preliminarily aligns with this hypothesis, although further experiments are required to confirm this. It has been well-documented that EPS plays an important role in mineralisation, specifically carbonate formation, in microbial ecosystems. EPS can enhance calcium carbonate precipitation by providing diffusion-limited sites that create alkalinity gradients in response to microbial processes [7]. It has been experimentally demonstrated that salinity influences the total amount of EPS (both loosely- and tightly-bound) which increases with increasing salinity [8]. The higher salinity at Guerrero Negro switches the microbial population from filamentous Microcoleus to being dominated by Phormidium, Oscillatoria, and unicellular cyanobacteria. We demonstrated that lipid binding into kerogen via strong covalent linkages occurs at the very earliest stages of sedimentary diagenesis. We are currently investigating the influence of salinity on organic preservation through experimental and modelling approaches. Understanding how key biomarkers transform into preserved organic matter (i.e., from precursor biolipids to bound geolipids) in brine ecosystems will aid the search for organic biosignatures on other planetary bodies, especially Icy Moons and modern Mars.

hypersaline↗

Rocket/Radar Sporadic-E Experiment Conducted during the El Coqui 2 Campaign

In order to investigate the complex electrodynamics and neutral-plasma coupling inherent to sporadic-E layers in the earth's mid-latitude ionosphere, a series of rocket/radar experiments were planned as part of the NASA El Coqui H Campaign from Tortuguero Launch Range, Puerto Rico, in March-April, 1998. The rocket experiments consisted of two pairs of "mother-daughter" payloads with limited apogees so that the payloads "hovered" in the sporadic-E region (95-125 km). Each payload pair included vector DC and AC electric field detectors, a highly accurate flux-gate DC magnetometer, an ion mass spectrometer, an ionization gauge, and spaced-electric field receivers to measure the wavelength and phase velocity of the unstable plasma waves. Separate rockets were included to simultaneously carry aloft TMA trails to measure the neutral wind and its velocity shear, believed responsible for the sporadic-E layer formation. In addition to the rocket experiments, incoherent scatter radar measurements of plasma density and drift velocity were gathered almost every night during the 3 week campaign. Continuous VHF backscatter radar operations were carried out from a site near Salinas, Puerto Rico, where 3-m backscatter echoes were observed associated with sporadic-E and other types of low altitude ionospheric layers. Other radars that operated during the campaign included an HF backscatter system near Ponce, Puerto Rico, and a second VHF backscatter radar set up near Aguadila Puerto Rico. On 24 March 1998, one of the instrumented rockets was launched, attaining an apogee of 129 km. The payloads successfully pierced an intense sporadic-E layer observed by both the Arecibo radar and the in-situ density and ion mass spectrometer probes. In-situ DC electric fields revealed very low (about 1-2 mV/m) ambient fields with small amplitude structures of the same order. No high frequency (short scale) waves were observed, consistent with the VHF backscatter observations at the time of the launch. An overview of the observations will be presented.

Pfaff, R.F.↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗