Search NASA⌕ Search

SEARCH · Search NASA

Results for “Global constraints”

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 289 records · Page 16

A Novel Multi-Spacecraft Interplanetary Global Trajectory Optimization Transcription

As the frontier of space exploration continues to advance, so does the design complexity of future interplanetary missions. One avenue of this increasing complexity includes a class of designs known as "Distributed Spacecraft Missions"; missions where multiple spacecraft coordinate to perform shared objectives. Current approaches for the global trajectory optimization of these Multi-Vehicle Missions (MVMs) are prone to shortcomings including laborious iterative design, considerable human-in-the-loop effort, treatment of the multi-vehicle problem as multiple separate trajectory optimization subproblems (resulting in suboptimal solutions where the whole is less than the sum of its parts), and poor handling of coordination objectives and constraints. There are only a handful of software platforms in existence capable of fully-automated, rapid, interplanetary mission and systems global optimization including the Parallel Global Multiobjective Optimizer (PaGMO), the Gravity Assisted Low-thrust Local Optimization Program (GALLOP), and the Evolutionary Mission Trajectory Generator (EMTG). However, none of these tools is capable of performing such tasks for MVM designs. The work outlined in this paper lays the groundwork for a technique to begin addressing these shortcomings. We present a fully-automated technique which frames interplanetary MVMs as Multi-Objective, Multi-Agent Hybrid Optimal Control Problems (MOMA HOCP). First, the basic functionality of this technique is validated on the single-vehicle problem of reproducing the Cassini interplanetary cruise.

Interplanetary↗

A Novel Multi-Spacecraft Interplanetary Global Trajectory Optimization Transcription

As the frontier of space exploration continues to advance, so does the design complexity of future interplanetary missions. One avenue of this increasing complexity includes a class of designs known as ``Distributed Spacecraft Missions"; missions where multiple spacecraft coordinate to perform shared objectives. Current approaches for the global trajectory optimization of these Multi-Vehicle Missions (MVMs) are prone to shortcomings including laborious iterative design, considerable human-in-the-loop effort, treatment of the multi-vehicle problem as multiple separate trajectory optimization subproblems (resulting in suboptimal solutions where the whole is less than the sum of its parts), and poor handling of coordination objectives and constraints. There are only a handful of software platforms in existence capable of fully-automated, rapid, interplanetary mission and systems global optimization including the Parallel Global Multiobjective Optimizer (PaGMO), the Gravity Assisted Low-thrust Local Optimization Program (GALLOP), and the Evolutionary Mission Trajectory Generator (EMTG). However, none of these tools is capable of performing such tasks for MVM designs. The work outlined in this paper lays the groundwork for a technique to begin addressing these shortcomings. We present a fully-automated technique which frames interplanetary MVMs as Multi-Objective, Multi-Agent Hybrid Optimal Control Problems (MOMA HOCP). First, the basic functionality of this technique is validated on the single-vehicle problem of reproducing the Cassini interplanetary cruise.

Napier, Sean W.↗

Channel Water Storage Anomaly: A New Remotely Sensed Quantity for Global River Analysis

River channels store large volumes of water globally, critically impacting ecological and biogeochemical processes. Despite the importance of river channel storage, there is not yet an observational constraint on this quantity. We introduce a 26-year record of entirely remotely sensed volumetric channel water storage (CWS) change on 26 major world rivers. We find mainstem volumetric CWS climatology amplitude (CA) represents an appreciable amount of basin-wide terrestrial water storage variability (median 2.78%, range 0.04%–12.54% across world rivers), despite mainstem rivers themselves represent an average of just 0.2% of basin area. We find that two global river routing schemes coupled with land surface models reasonably approximate CA (within ±50%) in only 11.5% (CaMa-Flood) and 30.7% (HyMap) of rivers considered. These findings demonstrate volumetric CWS is a useful quantity for assessing global hydrological model performance, and for advancing understanding of spatial patterns in global hydrology.

Stephen Coss↗

Optimizing GRACE/GRACE-FO Data and A priori Hydrological Knowledge for Improved Global Terrestial Water Storage Component Estimates

The comprehensive information of global terrestrial water storage (TWS) components (soil moisture, groundwater, snow, surface water) is essential for effective assessment of water resource availability, climate variation, and disaster mitigation measures. Observational data provided by the Gravity Recovery And Climate Experiment (GRACE) and GRACE Follow-On satellite missions offer global TWS variation (ΔTWS) in terms of an integrated water column. However, GRACE spatial resolution is relatively coarse (i.e., 3°), and the vertically integrated value cannot be separated into ΔTWS components directly. This study demonstrates the feasibility to estimate ΔTWS components at any desired spatial-vertical resolution by effectively maintaining the native resolution of the employed hydrological knowledge. It utilizes a least-squares with constraints (LSC) approach to rigorously incorporate GRACE and GRACE-FO data and a priori hydrological knowledge, with the aim to improve global ΔTWS components’ accuracy and spatial resolution. The 3°×3° GRACE mascon derived ΔTWS data is disaggregated into the 0.5°×0.5° anomalous soil moisture storage (ΔSMS), groundwater storage (ΔGWS), snow water equivalent (ΔSWE), and surface water storage (ΔSWS) based on the covariance information obtained from the Community Atmosphere Biosphere Land Exchange (CABLE) and the PCRaster Global Water Balance (PCR-GLOBWB) models. Evaluation with different ground measurements and satellite products between 2002 and 2019 exhibits significantly improved accuracy in all individual ΔTWS components. This improvement is of particular note in ΔGWS and ΔSWS, where the LSC approach increases the globally averaged correlation values by approximately 0.13 and 0.05, respectively. Reliable prior knowledge leads to a more accurate ΔTWS component estimate, and the use of ensemble-mean knowledge yields the best result.

GRACE↗

Global circulation and distribution of hydrogen in thermosphere of Venus

The global density distribution of atomic H in the upper atmosphere of Venus is determined in terms of a global circulation model. It is shown that H produced on the dayside is efficiently advected to the nightside by the major gas CO2 where it is then convected to lower altitudes and recombines. A maximum night/day H density ratio of about 5 is derived, in which case exospheric return flow from night to day strongly limits the nightside density enhancement. A previously postulated constraint connecting the eddy diffusion coefficient with the escape flux of H in one dimensional models is no longer required when global circulation of H is considered.

Hartle, R. E.↗

Powered Descent Guidance with General Thrust-Pointing Constraints

The Powered Descent Guidance (PDG) algorithm and software for generating Mars pinpoint or precision landing guidance profiles has been enhanced to incorporate thrust-pointing constraints. Pointing constraints would typically be needed for onboard sensor and navigation systems that have specific field-of-view requirements to generate valid ground proximity and terrain-relative state measurements. The original PDG algorithm was designed to enforce both control and state constraints, including maximum and minimum thrust bounds, avoidance of the ground or descent within a glide slope cone, and maximum speed limits. The thrust-bound and thrust-pointing constraints within PDG are non-convex, which in general requires nonlinear optimization methods to generate solutions. The short duration of Mars powered descent requires guaranteed PDG convergence to a solution within a finite time; however, nonlinear optimization methods have no guarantees of convergence to the global optimal or convergence within finite computation time. A lossless convexification developed for the original PDG algorithm relaxed the non-convex thrust bound constraints. This relaxation was theoretically proven to provide valid and optimal solutions for the original, non-convex problem within a convex framework. As with the thrust bound constraint, a relaxation of the thrust-pointing constraint also provides a lossless convexification that ensures the enhanced relaxed PDG algorithm remains convex and retains validity for the original nonconvex problem. The enhanced PDG algorithm provides guidance profiles for pinpoint and precision landing that minimize fuel usage, minimize landing error to the target, and ensure satisfaction of all position and control constraints, including thrust bounds and now thrust-pointing constraints.

Carson, John M., III↗

The enigma of lunar magnetism

Current understandings of the nature and probable origin of lunar magnetism are surveyed. Results of examinations of returned lunar samples are discussed which reveal the main carrier of the observed natural remanent magnetization to be iron, occasionally alloyed with nickel and cobalt, but do not distinguish between thermoremanent and shock remanent origins, and surface magnetometer data is presented, which indicates small-scale magnetic fields with a wide range of field intensities implying localized, near-surface sources. A detailed examination is presented of orbital magnetometer and charged particle data concerning the geologic nature and origin of magnetic anomaly sources and the directional properties of the magnetization, which exhibit a random distribution except for a depletion in the north-south direction. A lunar magnetization survey with global coverage provided by a polar orbiting satellite is suggested as a means of placing stronger constraints on the origin of lunar crustal magnetization.

Hood, L. L.↗

Wrinkle ridges, reverse faulting, and the depth penetration of lithospheric stress in lunae planum, Mars

Tectonic features on a planetary surface are commonly used as constraints on models to determine the state of stress at the time the features formed. Quantitative global stress models applied to understand the formation of the Tharsis province on Mars constrained by observed tectonics have calculated stresses at the surface of a thin elastic shell and have neglected the role of vertical structure in influencing the predicted pattern of surface deformation. Wrinkle ridges in the Lunae Planum region of Mars form a conentric pattern of regularly spaced features in the eastern and southeastern part of Tharsis; they are formed due to compressional stresses related to the response of the Martian lithosphere to the Tharsis bulge. As observed in the exposures of valley walls in areas such as the Kasei Valles, the surface plains unit is underlain by an unconsolidated impact-generated megaregolith that grades with depth into structurally competent lithospheric basement. The ridges have alternatively been hypothesized to reflect deformation restricted to the surface plains unit ('thin skinned deformation') and deformation that includes the surface unit, megaregolith and basement lithosphere ('thick skinned deformation'). We have adopted a finite element approach to quantify the nature of deformation associated with the development of wrinkle ridges in a vertically stratified elastic lithosphere. We used the program TECTON, which contains a slippery node capability that allowed us to explicitly take into account the presence of reverse faults believed to be associated with the ridges. In this study we focused on the strain field in the vicinity of a single ridge when slip occurs along the fault. We considered two initial model geometries. In the first, the reverse fault was assumed to be in the surface plains unit, and in the second the initial fault was located in lithospheric basement, immediately beneath the weak megaregolith. We are interested in the conditions underwhich strain in the surface layer and basement either penetrates or fails to penetrate through the megaregolith. We thus address the conditions required for an initial basement fault to propagate through the megaregolith to the surface, as well as the effect of the megareolith on the strain tensor in the vicinity of a fault that nucleates in the surface plains unit.

Zuber, M. T.↗

Terrestrial Environmental Variables Derived From EOS Platform Sensors

The three main objectives of the overall project were: 1. Adaptation of environmental constraint methods to take advantage of EOS sensors, specifically, MODIS, ASTER, and Landsat-7, in addition to the PM AVHRR observations 2. Refinement of environmental constraint methods based on fundamental scientific knowledge. 3. Assessment of spatial scaling patterns in environmental constraint measurements to evaluate the potential biases and errors that occur when estimating regional and global-scale NPP patterns with moderate to coarse satellite observations. These goals were modified because, on one hand, MODIS data did not become available until after the first year of the project and because of project staffing issues at the University of Maryland., The OSU portion of the project contained a modest amount of funding and responsibility compared to the University of Maryland and the University of Toledo.

Stadler, Stephen J.↗

Nonlinear Attitude Filtering Methods

The extended Kalman filter (EKF) is the workhorse of real-time spacecraft attitude estimation. Since the group SO3 of rotation matrices has dimension three, most attitude determination EKFs use lower- dimensional attitude parameterizations than the nine-parameter attitude matrix itself. The fact that all three- parameter representations of SO3 are singular or discontinuous for certain attitudes has led to extended discussions of constraints and attitude representations in EKFs. The most successful EKF uses a nonsingular parameterization for the global attitude, which necessarily has more than three parameters, while employing a three-component representation for the attitude errors. This filter has become known as the Multiplicative Extended Kalman Filter. These issues are now well understood, however, and the EKF has performed admirably in the vast majority of attitude determination applications. Nevertheless, poor performance or even divergence arising from the linearization implicit in the EKF has led to the development of nonlinear filters, most recently sigma point or unscented filters and particle filters.

Markley, F. Landis↗

Maintaining Atmospheric Mass and Water Balance Within Reanalysis

This report describes the modifications implemented into the Goddard Earth Observing System Version-5 (GEOS-5) Atmospheric Data Assimilation System (ADAS) to maintain global conservation of dry atmospheric mass as well as to preserve the model balance of globally integrated precipitation and surface evaporation during reanalysis. Section 1 begins with a review of these global quantities from four current reanalysis efforts. Section 2 introduces the modifications necessary to preserve these constraints within the atmospheric general circulation model (AGCM), the Gridpoint Statistical Interpolation (GSI) analysis procedure, and the Incremental Analysis Update (IAU) algorithm. Section 3 presents experiments quantifying the impact of the new procedure. Section 4 shows preliminary results from its use within the GMAO MERRA-2 Reanalysis project. Section 5 concludes with a summary.

IAU↗

Semi-analytic preliminary design of low-thrust missions

Using generalized logarithmic spirals to approximate low-thrust trajectories, a new strategy for the design of low-thrust gravity-assist transfers has been developed. Each transfer leg is defined by a semi-analytic model, and its solution is equivalent to a hybrid Lambert’s problem. The method is suitable for approximating both flyby and rendezvous transfer legs. A branch and prune algorithm is used to generate a collection of initial guesses for further optimization. The analytic nature of the low-thrust model simplifies the pruning step, since dynamical and operational constraints (like maximum thrust or total v) can be imposed easily. The solutions obtained with the global search algorithm can be post-processed, filtered, and ranked according to various criteria. This is where the versatility of the method resides, because changing the selection criteria does not require a new search. Selected candidates are then optimized further, in order to generate actual low-thrust orbits. Two mission design examples are presented: an asteroid deflection mission using a kinetic impactor, and a rendezvous mission to Jupiter. These examples are used to analyze the convergence of the optimization stage, in particular how far from the optimal solution the initial guesses are.

Park, Ryan S.↗

Continental-pelagic carbonate partitioning and the global carbonate-silicate cycle

A carbonate-silicate geochemical cycle model is developed and used to explore dynamic and climatic consequences of constraints on shallow-water carbonate burial and possible carbon loss to the mantle associated with sea-floor subduction. The model partitions carbonate deposition between shallow-water and deep-water environments and includes carbon fluxes between the mantle and lithosphere. When total lithospheric carbonate mass is constant, there are two stable steady states, one in which the carbonate burial flux is mostly continental and another in which it is mostly pelagic. The continental steady state is characterized by a low metamorphic CO2 flux to the atmosphere and predominantly shallow-water carbonate burial. The pelagic steady state is characterized by a high metamorphic CO2 flux and predominantly deep-water carbonate burial. For reasonable parameter values, when total lithospheric carbonate mass is allowed to vary, the model oscillates between predominantly continental and predominantly pelagic modes. Model results suggest that carbonate deposition patterns established during the Cenozoic may be pushing the Earth system from the continental to the pelagic mode on a time scale of 10(8) yr, with a possible consequent order-of-magnitude increase in the metamorphic CO2 flux to the atmosphere.

NASA Program Exobiology↗

Constraints on Early Mars Evolution and Dichotomy Origin from Relaxation Modeling of Dichotomy Boundary in the Ismenius Region

The Martian dichotomy is a global feature separating the northern and southern hemispheres. The 3.5 - 4 Gyr old feature is manifested by a topographic difference of 2-6 km and crustal thickness difference of approx. 15 - 30 km between the two hemispheres. In the Ismenius region, sections of the boundary are characterized by a single scarp with a slope of approx. 20 deg. - 23 deg. and are believed to be among the most well preserved parts of the dichotomy boundary. The origin of the dichotomy is unknown. Endogenic hypotheses do not predict the steep slopes (scarps) of the dichotomy boundary. Exogenic models for forming the northern lowlands by impact cratering, associate the scarps along the dichotomy boundary with craters' rims, but are not globally consistent with the topography and gravity. In order to better understand the origin of the Martian dichotomy, it is necessary to know if the steep scarps along the boundary represent the original shape of the dichotomy. Smrekar et al. presented evidence showing that the boundary scarp in Ismenius is a fault along which the highland crust was down faulted. We test whether the relaxation process could produce faulting along the dichotomy boundary and examine the crustal and mantle conditions that would allow for faulting to occur within 1 Gyr and preserve the long wavelength topography over another 3 Gyr. We approach the problem by a combination of numerical and semi-analytical modeling. We test different viscosity profiles and crustal thicknesses by comparing our modeled magnitude, location and timing of plastic strain and displacements to detailed geologic observations in the Ismenius region.

Guest, A.↗

Advancing Solar Irradiance Measurement for Climate-Related Studies: Accurate Constraint on Direct Aerosol Radiative Effect (DARE)

Earth's climate is driven primarily by solar radiation. As summarized in various IPCC reports, the global average of radiative forcing for different agents and mechanisms, such as aerosols or CO2 doubling, is in the range of a few W/sq m. However, when solar irradiance is measured by broadband radiometers, such as the fleet of Eppley Precision Solar Pyranometers (PSP) and equivalent instrumentation employed worldwide, the measurement uncertainty is larger than 2% (e.g., WMO specification of pyranometer, 2008). Thus, out of the approx. 184 W/sq m (approx.263 W/sq m if cloud-free) surface solar insolation (Trenberth et al. 2009), the measurement uncertainty is greater than +/-3.6 W/sq m, overwhelming the climate change signals. To discern these signals, less than a 1 % measurement uncertainty is required and is currently achievable only by means of a newly developed methodology employing a modified PSP-like pyranometer and an updated calibration equation to account for its thermal effects (li and Tsay, 2010). In this talk, we will show that some auxiliary measurements, such as those from a collocated pyrgeometer or air temperature sensors, can help correct historical datasets. Additionally, we will also demonstrate that a pyrheliometer is not free of the thermal effect; therefore, comparing to a high cost yet still not thermal-effect-free "direct + diffuse" approach in measuring surface solar irradiance, our new method is more economical, and more likely to be suitable for correcting a wide variety of historical datasets. Modeling simulations will be presented that a corrected solar irradiance measurement has a significant impact on aerosol forcing, and thus plays an important role in climate studies.

Tsay, Si-Chee↗

Assessment of Anthropogenic and Climatic Impacts on the Global Carbon Cycle Using a 3-D Model Constrained by Isotopic Carbon Measurements and Remote Sensing of Vegetation

Our original proposal called for improved modeling of the terrestrial biospheric carbon cycle, specifically using biome-specific process models to account for both the energy and water budgets of plant growth, to facilitate investigations into recent changes in global atmospheric CO2 abundance and regional distribution. The carbon fluxes predicted by these models were to be incorporated into a global model of CO2 transport to establish large-scale regional fluxes of CO2 to and from the terrestrial biosphere subject to constraints imposed by direct measurements of atmospheric CO2 and its 13C/12C isotopic ratio. Our work was coordinated with a NASA project (NASA NAGW-3151) at the University of Montana under the direction of Steven Running, and was partially funded by the Electric Power Research Institute. The primary objective of this project was to develop and test the Biome-BGC model, a global biological process model with a daily time step which simulates the water, energy and carbon budgets of plant growth. The primary product, the unique global gridded daily land temperature, and the precipitation data set which was used to drive the process model is described. The Biome-BGC model was tested by comparison with a simpler biological model driven by satellite-derived (NDVI) Normalized Difference Vegetation Index and (PAR) Photosynthetically Active Radiation data and by comparison with atmospheric CO2 observations. The simple NDVI model is also described. To facilitate the comparison with atmospheric CO2 observations, a three-dimensional atmospheric transport model was used to produce predictions of atmospheric CO2 variations given CO2 fluxes owing to (NPP) Net Primary Productivity and heterotrophic respiration that were produced by the Biome-BGC model and by the NDVI model. The transport model that we used in this project, and errors associated with transport simulations, were characterized by a comparison of 12 transport models.

Keeling, Charles D.↗

The SMAP Level 4 Carbon PRODUCT for Monitoring Terrestrial Ecosystem-Atmosphere CO2 Exchange

The NASA Soil Moisture Active Passive (SMAP) mission Level 4 Carbon (L4_C) product provides model estimates of Net Ecosystem CO2 exchange (NEE) incorporating SMAP soil moisture information as a primary driver. The L4_C product provides NEE, computed as total respiration less gross photosynthesis, at a daily time step and approximate 14-day latency posted to a 9-km global grid summarized by plant functional type. The L4_C product includes component carbon fluxes, surface soil organic carbon stocks, underlying environmental constraints, and detailed uncertainty metrics. The L4_C model is driven by the SMAP Level 4 Soil Moisture (L4_SM) data assimilation product, with additional inputs from the Goddard Earth Observing System, Version 5 (GEOS-5) weather analysis and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. The L4_C data record extends from March 2015 to present with ongoing production. Initial comparisons against global CO2 eddy flux tower measurements, satellite Solar Induced Canopy Florescence (SIF) and other independent observation benchmarks show favorable L4_C performance and accuracy, capturing the dynamic biosphere response to recent weather anomalies and demonstrating the value of SMAP observations for monitoring of global terrestrial water and carbon cycle linkages.

SMAP↗

The Role of CO2 in Aqueous Alteration of Ultra-Mafic Rocks and the Formation of MF-,FE-Rich Aqueous Solutons on Early Mars

An adequate understanding of water on Mars that moves beyond the simplistic "warmwet" vs. "cold-dry" dichotomy must include strong constraints on the variables: water/rock ratio, time, temperature, and chemical composition. By constraining these variables first on local, then regional and global scales we will be capable of precisely targeting landed missions to definitively understand the history of water on Mars and the possible existence of life. Data from remote sensing of Mars, landed missions, and martian meteorites indicate that secondary minerals formed from aqueous fluids on Mars are predominately Fe- and Mg-rich. The unique Mg-, Fe-rich carbonates in the ALH 84001 meteorite provide an excellent opportunity to provide strong constraints on an Fe-, Mg-rich aqueous system on early Mars. This work seeks to use the unusual chemical compositions of the ALH 84001 carbonates as a constraint for the composition of their formation fluid. These constraints can be used to better understand aqueous processes at a critical time in martian history.

Niles, Paul B.↗