Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constraint Checking”

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 109 records · Page 6

Checking Flight Rules with TraceContract: Application of a Scala DSL for Trace Analysis

Typically during the design and development of a NASA space mission, rules and constraints are identified to help reduce reasons for failure during operations. These flight rules are usually captured in a set of indexed tables, containing rule descriptions, rationales for the rules, and other information. Flight rules can be part of manual operations procedures carried out by humans. However, they can also be automated, and either implemented as on-board monitors, or as ground based monitors that are part of a ground data system. In the case of automated flight rules, one considerable expense to be addressed for any mission is the extensive process by which system engineers express flight rules in prose, software developers translate these requirements into code, and then both experts verify that the resulting application is correct. This paper explores the potential benefits of using an internal Scala DSL for general trace analysis, named TRACECONTRACT, to write executable specifications of flight rules. TRACECONTRACT can generally be applied to analysis of for example log files or for monitoring executing systems online.

temporal logic↗

Proceedings of the Second NASA Formal Methods Symposium

This publication contains the proceedings of the Second NASA Formal Methods Symposium sponsored by the National Aeronautics and Space Administration and held in Washington D.C. April 13-15, 2010. Topics covered include: Decision Engines for Software Analysis using Satisfiability Modulo Theories Solvers; Verification and Validation of Flight-Critical Systems; Formal Methods at Intel -- An Overview; Automatic Review of Abstract State Machines by Meta Property Verification; Hardware-independent Proofs of Numerical Programs; Slice-based Formal Specification Measures -- Mapping Coupling and Cohesion Measures to Formal Z; How Formal Methods Impels Discovery: A Short History of an Air Traffic Management Project; A Machine-Checked Proof of A State-Space Construction Algorithm; Automated Assume-Guarantee Reasoning for Omega-Regular Systems and Specifications; Modeling Regular Replacement for String Constraint Solving; Using Integer Clocks to Verify the Timing-Sync Sensor Network Protocol; Can Regulatory Bodies Expect Efficient Help from Formal Methods?; Synthesis of Greedy Algorithms Using Dominance Relations; A New Method for Incremental Testing of Finite State Machines; Verification of Faulty Message Passing Systems with Continuous State Space in PVS; Phase Two Feasibility Study for Software Safety Requirements Analysis Using Model Checking; A Prototype Embedding of Bluespec System Verilog in the PVS Theorem Prover; SimCheck: An Expressive Type System for Simulink; Coverage Metrics for Requirements-Based Testing: Evaluation of Effectiveness; Software Model Checking of ARINC-653 Flight Code with MCP; Evaluation of a Guideline by Formal Modelling of Cruise Control System in Event-B; Formal Verification of Large Software Systems; Symbolic Computation of Strongly Connected Components Using Saturation; Towards the Formal Verification of a Distributed Real-Time Automotive System; Slicing AADL Specifications for Model Checking; Model Checking with Edge-valued Decision Diagrams; and Data-flow based Model Analysis.

Munoz, Cesar↗

Kepler and Ground-Based Transits of the exo-Neptune HAT-P-11b

We analyze 26 archival Kepler transits of the exo-Neptune HAT-P-11b, supplemented by ground-based transits observed in the blue (B band) and near-IR (J band). Both the planet and host star are smaller than previously believed; our analysis yields Rp = 4.31 R xor 0.06 R xor and Rs = 0.683 R solar mass 0.009 R solar mass, both about 3 sigma smaller than the discovery values. Our ground-based transit data at wavelengths bracketing the Kepler bandpass serve to check the wavelength dependence of stellar limb darkening, and the J-band transit provides a precise and independent constraint on the transit duration. Both the limb darkening and transit duration from our ground-based data are consistent with the new Kepler values for the system parameters. Our smaller radius for the planet implies that its gaseous envelope can be less extensive than previously believed, being very similar to the H-He envelope of GJ 436b and Kepler-4b. HAT-P-11 is an active star, and signatures of star spot crossings are ubiquitous in the Kepler transit data. We develop and apply a methodology to correct the planetary radius for the presence of both crossed and uncrossed star spots. Star spot crossings are concentrated at phases 0.002 and +0.006. This is consistent with inferences from Rossiter-McLaughlin measurements that the planet transits nearly perpendicular to the stellar equator. We identify the dominant phases of star spot crossings with active latitudes on the star, and infer that the stellar rotational pole is inclined at about 12 deg 5 deg to the plane of the sky. We point out that precise transit measurements over long durations could in principle allow us to construct a stellar Butterfly diagram to probe the cyclic evolution of magnetic activity on this active K-dwarf star.

Deming, Drake↗

Two-Season Atacama Cosmology Telescope Polarimeter Lensing Power Spectrum

We report a measurement of the power spectrum of cosmic microwave background (CMB) lensing from two seasons of Atacama Cosmology Telescope polarimeter (ACTPol) CMB data. The CMB lensing power spectrum is extracted from both temperature and polarization data using quadratic estimators. We obtain results that are consistent with the expectation from the best-fit Planck CDM model over a range of multipoles L 80-2100, with an amplitude of lensing A(sub lens) = 1.06 +/- 0.15 stat +/- 0.06 sys relative to Planck. Our measurement of the CMB lensing power spectrum gives sigma 8 omega m(sup 0.25) = 0.643 +/- 0.054; including baryon acoustic oscillation scale data, we constrain the amplitude of density fluctuations to be sigma 8 = 0.831 +/- 0.053. We also update constraints on the neutrino mass sum. We verify our lensing measurement with a number of null tests and systematic checks, finding no evidence of significant systematic errors. This measurement relies on a small fraction of the ACTPol data already taken; more precise lensing results can therefore be expected from the full ACTPol data set.

Shewin, Blake D.↗

Exploring Network-Related Optimization Problems Using Quantum Heuristics

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics.One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Wang, Zhihui↗

Study network-related optimization problems using quantum alternating optimization ansatz

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics. One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Zhihui Wang↗

Models for Water Isotopes Constrained with Data from Crystal Face

During the year covered by this proposal we conducted work on several different topics, as reflected by our publications listed below. One major activity was to work with a group of about 10 scientists from around the country to prepare a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. Unfortunately budgetary constraints prevented these missions from moving forward. In conjunction with the group NASA Ames we built a new numerical model for deep convection and have applied that model to simulate the CRYSTAL isotope data. Our goal in particular has been to better understanding how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results are now in press in Geophys. Res. Lett. Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity ana isotope abundances with this model. The data suggest that isotopes do not provide a clear constraint on the mechanism by which water enters the stratosphere- whether by convection, or by slow ascent. Our work is relevant for the water isotope comparison experiments recently done by NASA. We are conducting numerical experiments related to this project to help understand that data.

Toon, Owen B.↗

Spaceborne Autonomous and Ground Based Relative Orbit Control for the TerraSAR-X/TanDEM-X Formation

TerraSAR-X (TSX) and TanDEM-X (TDX) are two advanced synthetic aperture radar (SAR) satellites flying in formation. SAR interferometry allows a high resolution imaging of the Earth by processing SAR images obtained from two slightly different orbits. TSX operates as a repeat-pass interferometer in the first phase of its lifetime and will be supplemented after two years by TDX in order to produce digital elevation models (DEM) with unprecedented accuracy. Such a flying formation makes indeed possible a simultaneous interferometric data acquisition characterized by highly flexible baselines with range of variations between a few hundreds meters and several kilometers [1]. TSX has been successfully launched on the 15th of June, 2007. TDX is expected to be launched on the 31st of May, 2009. A safe and robust maintenance of the formation is based on the concept of relative eccentricity/inclination (e/i) vector separation whose efficiency has already been demonstrated during the Gravity Recovery and Climate Experiment (GRACE) [2]. Here, the satellite relative motion is parameterized by mean of relative orbit elements and the key idea is to align the relative eccentricity and inclination vectors to minimize the hazard of a collision. Previous studies have already shown the pertinence of this concept and have described the way of controlling the formation using an impulsive deterministic control law [3]. Despite the completely different relative orbit control requirements, the same approach can be applied to the TSX/TDX formation. The task of TDX is to maintain the close formation configuration by actively controlling its relative motion with respect to TSX, the leader of the formation. TDX must replicate the absolute orbit keeping maneuvers executed by TSX and also compensate the natural deviation of the relative e/i vectors. In fact the relative orbital elements of the formation tend to drift because of the secular non-keplerian perturbations acting on both satellites. The goal of the ground segment is thus to regularly correct this configuration by performing small orbit correction maneuvers on TDX. The ground station contacts are limited due to the geographic position of the station and the costs for contact time. Only with a polar ground station a contact visibility is possible every orbit for LEO satellites. TSX and TDX use only the Weilheim ground station (in the southern part of Germany) during routine operations. This station allows two scheduled contact per day for the nominal orbit configuration, meaning that the satellite conditions can be checked with an interval of 12 hours. While this limitation is usually not critical for single satellite operations, the visibility constraints drive the achievable orbit control accuracy for a LEO formation if a ground based approach is chosen. Along-track position uncertainties and maneuver execution errors affect the relative motion and can be compensated only after a ground station contact.

Ardaens, J. S.↗

Cloud Condensation in Titan's Lower Stratosphere

A 1-D condensation model is developed for the purpose of reproducing ice clouds in Titan's lower stratosphere observed by the Composite Infrared Spectrometer (CIRS) onboard Cassini. Hydrogen cyanide (HCN), cyanoacetylene (HC3N), and ethane (C2H6) vapors are treated as chemically inert gas species that flow from an upper boundary at 500 km to a condensation sink near Titan's tropopause (-45 km). Gas vertical profiles are determined from eddy mixing and a downward flux at the upper boundary. The condensation sink is based upon diffusive growth of the cloud particles and is proportional to the degree of supersaturation in the cloud formation regIOn. Observations of the vapor phase abundances above the condensation levels and the locations and properties of the ice clouds provide constraints on the free parameters in the model. Vapor phase abundances are determined from CIRS mid-IR observations, whereas cloud particle sizes, altitudes, and latitudinal distributions are derived from analyses of CIRS far-IR observations of Titan. Specific cloud constraints include: I) mean particle radii of2-3 J.lm inferred from the V6 506 cm- band of HC3N, 2) latitudinal abundance distributions of condensed nitriles, inferred from a composite emission feature that peaks at 160/cm , and 3) a possible hydrocarbon cloud layer at high latitudes, located near an altitude of 60 km, which peaks between 60 and 80 cm l . Nitrile abundances appear to diminish substantially at high northern latitudes over the time period 2005 to 2010 (northern mid winter to early spring). Use of multiple gas species provides a consistency check on the eddy mixing coefficient profile. The flux at the upper boundary is the net column chemical production from the upper atmosphere and provides a constraint on chemical pathways leading to the production of these compounds. Comparison of the differing lifetimes, vapor phase transport, vapor phase loss rate, and particle sedimentation, sheds light on temporal stability of the clouds.

Romani, Paul N.↗

Multiply scaled constrained nonlinear equation solvers

To improve the numerical stability of nonlinear equation solvers, a partitioned multiply scaled constraint scheme is developed. This scheme enables hierarchical levels of control for nonlinear equation solvers. To complement the procedure, partitioned convergence checks are established along with self-adaptive partitioning schemes. Overall, such procedures greatly enhance the numerical stability of the original solvers. To demonstrate and motivate the development of the scheme, the problem of nonlinear heat conduction is considered. In this context the main emphasis is given to successive substitution-type schemes. To verify the improved numerical characteristics associated with partitioned multiply scaled solvers, results are presented for several benchmark examples.

Padovan, Joe↗

A Framework for a Supervisory Expert System for Robotic Manipulators with Joint-Position Limits and Joint-Rate Limits

This report addresses the problem of path planning and control of robotic manipulators which have joint-position limits and joint-rate limits. The manipulators move autonomously and carry out variable tasks in a dynamic, unstructured and cluttered environment. The issue considered is whether the robotic manipulator can achieve all its tasks, and if it cannot, the objective is to identify the closest achievable goal. This problem is formalized and systematically solved for generic manipulators by using inverse kinematics and forward kinematics. Inverse kinematics are employed to define the subspace, workspace and constrained workspace, which are then used to identify when a task is not achievable. The closest achievable goal is obtained by determining weights for an optimal control redistribution scheme. These weights are quantified by using forward kinematics. Conditions leading to joint rate limits are identified, in particular it is established that all generic manipulators have singularities at the boundary of their workspace, while some have loci of singularities inside their workspace. Once the manipulator singularity is identified the command redistribution scheme is used to compute the closest achievable Cartesian velocities. Two examples are used to illustrate the use of the algorithm: A three link planar manipulator and the Unimation Puma 560. Implementation of the derived algorithm is effected by using a supervisory expert system to check whether the desired goal lies in the constrained workspace and if not, to evoke the redistribution scheme which determines the constraint relaxation between end effector position and orientation, and then computes optimal gains.

Mutambara, Arthur G. O.↗

Tests of General Relativity With the Binary Black Hole Signals from the LIGO-Virgo Catalog GWTC-1

The detection of gravitational waves by Advanced LIGO and Advanced Virgo provides an opportunity to test general relativity in a regime that is inaccessible to traditional astronomical observations and laboratory tests. We present four tests of the consistency of the data with binary black hole gravitational waveforms predicted by general relativity. One test subtracts the best-fit waveform from the data and checks the consistency of the residual with detector noise. The second test checks the consistency of the low- and high-frequency parts of the observed signals. The third test checks that phenomenological deviations introduced in the waveform model (including in the post-Newtonian coefficients) are consistent with 0. The fourth test constrains modifications to the propagation of gravitational waves due to a modified dispersion relation, including that from a massive graviton. We present results both for individual events and also results obtained by combining together particularly strong events from the first and second observing runs of Advanced LIGO and Advanced Virgo, as collected in the catalog GWTC-1. We do not find any inconsistency of the data with the predictions of general relativity and improve our previously presented combined constraints by factors of 1.1 to 2.5. In particular, we bound the mass of the graviton to be 𝑚𝑔≤4.7×10 −23 eV/𝑐 2 (90% credible level), an improvement of a factor of 1.6 over our previously presented results. Additionally, we check that the four gravitational-wave events published for the first time in GWTC-1 do not lead to stronger constraints on alternative polarizations than those published previously.

B P Abbott↗

Timing analysis by model checking

The safety of modern avionics relies on high integrity software that can be verified to meet hard real-time requirements. The limits of verification technology therefore determine acceptable engineering practice. To simplify verification problems, safety-critical systems are commonly implemented under the severe constraints of a cyclic executive, which make design an expensive trial-and-error process highly intolerant of change. Important advances in analysis techniques, such as rate monotonic analysis (RMA), have provided a theoretical and practical basis for easing these onerous restrictions. But RMA and its kindred have two limitations: they apply only to verifying the requirement of schedulability (that tasks meet their deadlines) and they cannot be applied to many common programming paradigms. We address both these limitations by applying model checking, a technique with successful industrial applications in hardware design. Model checking algorithms analyze finite state machines, either by explicit state enumeration or by symbolic manipulation. Since quantitative timing properties involve a potentially unbounded state variable (a clock), our first problem is to construct a finite approximation that is conservative for the properties being analyzed-if the approximation satisfies the properties of interest, so does the infinite model. To reduce the potential for state space explosion we must further optimize this finite model. Experiments with some simple optimizations have yielded a hundred-fold efficiency improvement over published techniques.

Naydich, Dimitri↗

Optimal Network-Topology Design

Candidate network designs tested for acceptability and cost. Optimal Network Topology Design computer program developed as part of study on topology design and analysis of performance of Space Station Information System (SSIS) network. Uses efficient algorithm to generate candidate network designs consisting of subsets of set of all network components, in increasing order of total costs and checks each design to see whether it forms acceptable network. Technique gives true cost-optimal network and particularly useful when network has many constraints and not too many components. Program written in PASCAL.

Li, Victor O. K.↗

Self-adaptive predictor-corrector algorithm for static nonlinear structural analysis

A multiphase selfadaptive predictor corrector type algorithm was developed. This algorithm enables the solution of highly nonlinear structural responses including kinematic, kinetic and material effects as well as pro/post buckling behavior. The strategy involves three main phases: (1) the use of a warpable hyperelliptic constraint surface which serves to upperbound dependent iterate excursions during successive incremental Newton Ramphson (INR) type iterations; (20 uses an energy constraint to scale the generation of successive iterates so as to maintain the appropriate form of local convergence behavior; (3) the use of quality of convergence checks which enable various self adaptive modifications of the algorithmic structure when necessary. The restructuring is achieved by tightening various conditioning parameters as well as switch to different algorithmic levels to improve the convergence process. The capabilities of the procedure to handle various types of static nonlinear structural behavior are illustrated.

Padovan, J.↗

Influence of Nucleation Mechanisms on the Radiative Properties of Deep Convective Clouds and Subvisible Cirrus in CRYSTAL/FACE

During the past few years we have conducted work on several different topics, as reflected by our publications. As one of the Co-Project scientists for The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL FACE) we worked to help design the mission and then conduct it in the field. Another major activity during the past two years has been to pull together various groups to formulate plans for follow on missions to CRYSTAL FACE. We organized a workshop at the University of Colorado during the summer of 2003 to assess the best locations for future missions. Working with a group of about 10 scientists from around the country we prepared a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. In conjunction with the group at NASA Ames we have helped build a new numerical model for deep convection and have applied that model to simulate the CRYSTAL data. Our goal in particular has been to better understand how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results have just been submitted to Geophys. Res. Lett (Smith et al., 2005). Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity and isotope abundances with this model. We have also been heavily involved in trying to improve our understanding of nitric acid condensation on ice. Gao et al (2004) have shown that water supersaturations above ice occur when the atmosphere is supersaturated with respect to nitric acid trihydrate. As one of the co-authors of that work, we suggested the mechanism that may explain why this is occurring. Essentially, ice does not like to grow near unit supersaturation, but does so because the water molecules can find sites on the ice surface to attach themselves to before they fly off the ice surface. This phenomena was well known in the 1960s when it was a source of debate about whether condensation and evaporation coefficients for ice would be the same. Evaporation does not require any molecular orientation, while condensation does, so it was possible that the coefficients would differ. They don't differ because the water molecules rapidly move across the surface and find places to attach. Nitric acid may be occupying these preferred sites and therefore the water molecules can't find a desirable place to attach. We anticipate that this research will be the subject of laboratory work during the coming few years. Another possibility that has been suggested is that cubic ice is forming in clouds. We have measured the vapor pressure of cubic ice, and plan to publish that result in the next few months. We have also been working on additional aspects of the condensation of nitric acid on ice. With Y. Kondo we studied the condensation of NOy on ice using the SOLVE data. Gamblin et al. have continued this work. The CRYSTAL NOy and HNO, groups have shown th their data can be fit using standard Langmuir isotherms as suggested in some, but not all, laboratory studies. We have found in the SOLVE data set that this is not the case. Moreover some laboratory studies show there are important kinetic effects that may be occurring in the atmosphere limiting the transfer of nitric acid to the ice. The SOLVE data seem consistent with these studies. We are currently re-analyzing the CRYSTAL data to look for these kinetic effects. There are a number of implications of these studies. One of the more interesting is that the nitric acid coating on ice can be used as a cloud clock to determine how long the cloud parcel has been in existence. We have also been involved with several laboratory studies. We have worked to improve the database on ice optical constants, which are critical for remote sensing. We have also studied the ways in which ice nucleates on clays. We suspect now that the standard theories used for depositional ice nucleation are completely incorrect. Further work will be needed to develop a new theory.

Toon, Owen B.↗

CO2 Insulation for Thermal Control of the Mars Science Laboratory

The National Aeronautics and Space Administration (NASA) is sending a large (>850 kg) rover as part of the Mars Science Laboratory (MSL) mission to Mars in 2011. The rover's primary power source is a Multi-Mission Radioisotope Thermoelectric Generator (MMRTG) that generates roughly 2000 W of heat, which is converted to approximately 110 W of electrical power for use by the rover electronics, science instruments, and mechanism-actuators. The large rover size and extreme thermal environments (cold and hot) for which the rover is designed for led to a sophisticated thermal control system to keep it within allowable temperature limits. The pre-existing Martian atmosphere of low thermal conductivity CO2 gas (8 Torr) is used to thermally protect the rover and its components from the extremely cold Martian environment (temperatures as low as -130 deg C). Conventional vacuum based insulation like Multi Layer Insulation (MLI) is not effective in a gaseous atmosphere, so engineered gaps between the warm rover internal components and the cold rover external structure were employed to implement this thermal isolation. Large gaps would lead to more thermal isolation, but would also require more of the precious volume available within the rover. Therefore, a balance of the degree of thermal isolation achieved vs. the volume of rover utilized is required to reach an acceptable design. The temperature differences between the controlled components and the rover structure vary from location to location so each gap has to be evaluated on a case-by-case basis to arrive at an optimal thickness. For every configuration and temperature difference, there is a critical thickness below which the heat transfer mechanism is dominated by simple gaseous thermal conduction. For larger gaps, the mechanism is dominated by natural convection. In general, convection leads to a poorer level of thermal isolation as compared to conduction. All these considerations play important roles in the optimization process. A three-step process was utilized to design this insulation. The first step is to come up with a simple, textbook based, closed-form equation assessment of gap thickness vs. resultant thermal isolation achieved. The second step is a more sophisticated numerical assessment using Computational Fluid Dynamics (CFD) software to investigate the effect of complicated geometries and temperature contours along them to arrive at the effective thermal isolation in a CO2 atmosphere. The third step is to test samples of representative geometries in a CO2 filled chamber to measure the thermal isolation achieved. The results of these assessments along with the consistency checks across these methods leads to the formulation of design-guidelines for gap implementation within the rover geometry. Finally, based on the geometric and functional constraints within the real rover system, a detailed design that accommodates all these factors is arrived at. This paper will describe in detail this entire process, the results of these assessments and the final design that was implemented.

CO2↗

Space Shuttle Day-of-Launch Trajectory Design Operations

A top priority of any launch vehicle is to insert as much mass into the desired orbit as possible. This requirement must be traded against vehicle capability in terms of dynamic control, thermal constraints, and structural margins. The vehicle is certified to specific structural limits which will yield certain performance characteristics of mass to orbit. Some limits cannot be certified generically and must be checked with each mission design. The most sensitive limits require an assessment on the day-of-launch. To further minimize vehicle loads while maximizing vehicle performance, a day-of-launch trajectory can be designed. This design is optimized according to that day s wind and atmospheric conditions, which increase the probability of launch. The day-of-launch trajectory design and verification process is critical to the vehicle s safety. The Day-Of-Launch I-Load Update (DOLILU) is the process by which the National Aeronautics and Space Administration's (NASA) Space Shuttle Program tailors the vehicle steering commands to fit that day s environmental conditions and then rigorously verifies the integrated vehicle trajectory s loads, controls, and performance. This process has been successfully used for almost twenty years and shares many of the same elements with other launch vehicles that execute a day-of-launch trajectory design or day-of-launch trajectory verification. Weather balloon data is gathered at the launch site and transmitted to the Johnson Space Center s Mission Control. The vehicle s first stage trajectory is then adjusted to the measured wind and atmosphere data. The resultant trajectory must satisfy loads and controls constraints. Additionally, these assessments statistically protect for non-observed dispersions. One such dispersion is the change in the wind from the last measured balloon to launch time. This process is started in the hours before launch and is repeated several times as the launch count proceeds. Should the trajectory design not meet all constraint criteria, Shuttle would be No-Go for launch. This Shuttle methodology is very similar to other unmanned launch vehicles. By extension, this method would likely be employed for any future NASA launch vehicle. This paper will review the Shuttle s day-of-launch trajectory optimization and verification operations as an example of a more generic application of day-of-launch design and validation. With Shuttle s retirement, it is fitting to document the current state of this critical process and capture lessons learned to benefit current and future launch vehicle endeavors.

Harrington, Brian E.↗