Engineering PapersSearch

SEARCH · Engineering Papers

Results for “state estimation”

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 19 records

Exploration of an Adaptive Routine for Battery Modeling

The purpose of this document is to explore the use of adaptive routines in battery modeling. The adaptive routines consist of real-time state estimators combined with battery parameter model components that are adjusted in real-time as battery data becomes available. Several aspects are explored. It is shown that model parameter identification is possible for simple battery models using available input/output data measurements. The online system identification used is recursive least squares. Model identification may be combined with a state observer such as the extended Kalman filter or the unscented Kalman filter to form an adaptive model combined with state estimation. However, such a combination is found to be problematic due to uncertainty, observability and stability issues. This paper is organized as follows. Section 1 introduces adaptive routines and possible roles they play in battery modeling. In Section 2 real-time parameter identification is described with results based on battery data. Section 3 reviews various state estimators and results using a simple battery model. In Section 4 parameter identification and state estimation are combined to form an adaptive routine. Finally, in Section 5 conclusions are drawn and future work is suggested.

Adaptive

Mathematical Characterization of Battery Models

The purpose of this document is to demonstrate the use of the Extended Kalman Filter as a tool for battery state estimation and the estimation of battery state of charge. The mathematical details based on the equivalent circuit model are presented followed by an electrochemical engineering model. A simplified first-order model is used to demonstrate the procedure followed by second and third-order models. Next a simplified electrochemistry model is presented along with observer development. State observability is calculated for the simpler equivalent circuit models and the simplified electrochemistry model. An outline of the battery model parameter identification method is presented, and model performance based on experimental and flight data is demonstrated.

Battery

Remote Guidance of a Mothership-Daughtership Formation for Asteroid Exploration

This work presents an application of Model Predictive Control to maintain long-term stability for a three-satellite formation around the asteroid Apophis via remote guidance. The mothership and two daughterships perform tomographic characterization with a bistatic radar. The mothership estimates the states of all three spacecraft and determines the necessary control to maintain the desired antipodal daughtership states. An Extended Kalman Filter estimates the mothership states. A Batch Least Squares algorithm estimates the daughtership states and performs maneuver reconstruction for the mothership-administered thrust commands under thruster and gravitational uncertainty. The formulation maintains the antipodal angle with growing instability amidst these uncertainties.

Guidance

Online Prediction of Battery Discharge and Estimation of Parasitic Loads for an Electric Aircraft

Predicting whether or not vehicle batteries contain sufficient charge to support operations over the remainder of a given flight plan is critical for electric aircraft. This paper describes an approach for identifying upper and lower uncertainty bounds on predictions that aircraft batteries will continue to meet output power and voltage requirements over the remainder of a flight plan. Battery discharge prediction is considered here in terms of the following components; (i) online battery state of charge estimation; (ii) prediction of future battery power demand as a function of an aircraft flight plan; (iii) online estimation of additional parasitic battery loads; and finally, (iv) estimation of flight plan safety. Substantial uncertainty is considered to be an irremovable part of the battery discharge prediction problem. However, high-confidence estimates of flight plan safety or lack of safety are shown to be generated from even highly uncertain prognostic predictions.

Battery Discharge Prediction

Modular Battery Charge Controller

A new approach to masterless, distributed, digital-charge control for batteries requiring charge control has been developed and implemented. This approach is required in battery chemistries that need cell-level charge control for safety and is characterized by the use of one controller per cell, resulting in redundant sensors for critical components, such as voltage, temperature, and current. The charge controllers in a given battery interact in a masterless fashion for the purpose of cell balancing, charge control, and state-of-charge estimation. This makes the battery system invariably fault-tolerant. The solution to the single-fault failure, due to the use of a single charge controller (CC), was solved by implementing one CC per cell and linking them via an isolated communication bus [e.g., controller area network (CAN)] in a masterless fashion so that the failure of one or more CCs will not impact the remaining functional CCs. Each micro-controller-based CC digitizes the cell voltage (V(sub cell)), two cell temperatures, and the voltage across the switch (V); the latter variable is used in conjunction with V(sub cell) to estimate the bypass current for a given bypass resistor. Furthermore, CC1 digitizes the battery current (I1) and battery voltage (V(sub batt) and CC5 digitizes a second battery current (I2). As a result, redundant readings are taken for temperature, battery current, and battery voltage through the summation of the individual cell voltages given that each CC knows the voltage of the other cells. For the purpose of cell balancing, each CC periodically and independently transmits its cell voltage and stores the received cell voltage of the other cells in an array. The position in the array depends on the identifier (ID) of the transmitting CC. After eight cell voltage receptions, the array is checked to see if one or more cells did not transmit. If one or more transmissions are missing, the missing cell(s) is (are) eliminated from cell-balancing calculations. The cell-balancing algorithm is based on the error between the cell s voltage and the other cells and is categorized into four zones of operation. The algorithm is executed every second and, if cell balancing is activated, the error variable is set to a negative low value. The largest error between the cell and the other cells is found and the zone of operation determined. If the error is zero or negative, then the cell is at the lowest voltage and no balancing action is needed. If the error is less than a predetermined negative value, a Cell Bad Flag is set. If the error is positive, then cell balancing is needed, but a hysteretic zone is added to prevent the bypass circuit from triggering repeatedly near zero error. This approach keeps the cells within a predetermined voltage range.

Button, Robert

Performance of Two Battery Prognostic Applications used by Two Octocopters for Safe Low Altitude Autonomous Flight Operations

This paper addresses the problem of building trust in online predictions of the remaining available flying time for two different electric Unmanned Aerial Vehicles (eUAVs) powered by lithium-ion-polymer batteries. Flight tests for various automation research missions for the two vehicles were monitored using two on-board battery health management applications to make predictions of the remaining flying time (RFT) for each eUAV and to predict the state of the battery. Playback of the voltage, current and temperature profiles of the battery discharge were used to assess the accuracy of the estimation of the voltage and the charge states of the models as well as the estimate of the RFT. The reference ground truth values were the observed landing time and the measured battery pack resting pack voltage 20 minutes after the flight. The predicted RFT, state of charge (SoC), and state of energy (SoE) were compared with the observed results. Noise values of one standard deviation from the mean values of the internal charge states of the battery model during a reference run were used to vary the states during simulation. One application used an equivalent circuit model of the electrical dynamics of the battery pack, and the other application used a reduced-order electrochemistry model. The variation of the model state components was compared to the variation in the estimate of the RFT and the variation in the SoE to estimate a confidence factor. Variation in the estimates caused by factors affecting the off-line laboratory parameter identification experiments is considered. Variation in the estimates due to environmental factors are discussed.

Assurance

Incremental Learning for Passive Microwave Precipitation Retrievals using Advanced Technology Microwave Sounder

Spaceborne passive microwave (PMW) radiometry is central to global precipitation monitoring, yet retrieval uncertainties remain substantial, particularly for cross-track sounders whose variable footprints and channel configurations are optimized for atmospheric temperature and moisture profiling rather than precipitation. Consequently, existing operational products often exhibit angular-dependent biases, limited effective swath utilization, unrealistic rainfall probability distributions, and systematic misclassification of precipitation phase. These limitations are further compounded by the scarcity of globally accurate and representative precipitation observations, as training data from the Dual-frequency Precipitation Radar (DPR) and the Cloud Profiling Radar (CPR) are spatially sparse, lack uniform global coverage, and exhibit heterogeneous error characteristics across precipitation regimes. To address these challenges, this study presents a supervised retrieval algorithm that incrementally trains an ensemble of extreme gradient-boosted decision trees by augmenting base learners with pre-training on reanalysis data and post-training on coincident DPR and CPR observations matched with the Advanced Technology Microwave Sounder (ATMS). By transferring prior information from reanalysis to posterior constraints from radar observations and adopting a sequential detection–estimation strategy for precipitation phase and rate retrieval, the proposed approach yields retrievals across the full ATMS swath that are largely free from persistent deficiencies in current Global Precipitation Measurement (GPM) passive microwave operational products. In particular, the method resolves bimodal artifacts in rainfall retrievals and mitigates systematic high-latitude snowfall biases, including overestimation across the Arctic and underestimation across the Antarctic. Validation against independent Multi-Radar Multi-Sensor (MRMS) data over the Contiguous United States (CONUS) further demonstrates improved performance in precipitation phase detection and rate estimation relative to both reanalysis and current GPM PMW products.

Mahyar Garshasbi

Control of Thermo-Acoustics Instabilities: The Multi-Scale Extended Kalman Approach

"Multi-Scale Extended Kalman" (MSEK) is a novel model-based control approach recently found to be effective for suppressing combustion instabilities in gas turbines. A control law formulated in this approach for fuel modulation demonstrated steady suppression of a high-frequency combustion instability (less than 500Hz) in a liquid-fuel combustion test rig under engine-realistic conditions. To make-up for severe transport-delays on control effect, the MSEK controller combines a wavelet -like Multi-Scale analysis and an Extended Kalman Observer to predict the thermo-acoustic states of combustion pressure perturbations. The commanded fuel modulation is composed of a damper action based on the predicted states, and a tones suppression action based on the Multi-Scale estimation of thermal excitations and other transient disturbances. The controller performs automatic adjustments of the gain and phase of these actions to minimize the Time-Scale Averaged Variances of the pressures inside the combustion zone and upstream of the injector. The successful demonstration of Active Combustion Control with this MSEK controller completed an important NASA milestone for the current research in advanced combustion technologies.

Dzu K Le

X-Ray Polarization of Z-Type Neutron Star Low-Mass X-Ray Binaries II. Spectropolarimetric Analysis

IXPE has provided for the first time detailed energy- and time-resolved X-ray polarimetry of Z-type neutron star low-mass X-ray binaries (NS-LMXBs) as they move along their color-color diagrams (CCDs). These sources can reach the highest polarization observed for NS-LXMBs in the 2–8 keV range when they move along the horizontal branch. In a previous paper, we characterized the spectral state of a sample of Z-sources using the CCD and estimated the polarization with model-independent analysis. Here, we present detailed spectropolarimetric analysis for each source on each branch using data from IXPE, NICER, and NuSTAR. The continuum X-ray emission of all the sources is well described with a combination of thermal accretion disk emission plus a harder Comptonized component. In addition, reflection features, in particular the relativistically broadened Fe Kα line, are observed for our sources, except GX 5–1. For most of the sources and branches, the main contribution to the X-ray emission and polarization is due to Comptonization: moving from the horizontal branch (HB) to the normal branch (NB), the polarization degree (PD) in the 2–8 keV band varies from about 6% to 3–4%, while the PD is loosely constrained in the flaring branch (FB), due to the shorter exposures. These PD values are significantly higher than theoretical expectations for typical spreading or boundary layer configurations. The polarization of the disk is generally lower (below 3%) but still higher than predictions for an electron scattering-dominated, plane-parallel atmosphere above the disk observed at the corresponding inclination. Moreover, the polarization angle (PA) of the disk seems to be significantly misaligned and not perpendicular to that of Comptonization. We find no correlation between the polarization signal and the inclination, nor with the contribution of reflected photons throughout the Z-track.

accretion

The Growth of Zeolites A, X and Mordenite in Space

Zeolites are a class of crystalline aluminosilicate materials that form the backbone of the chemical process industry worldwide. They are used primarily as adsorbents and catalysts and support to a significant extent the positive balance of trade realized by the chemical industry in the United States (around $19 billion in 1991). The magnitude of their efforts can be appreciated when one realizes that since their introduction as 'cracking catalysts' in the early 1960's, they have saved the equivalent of 60 percent of the total oil production from Alaska's North Slope. Thus the performance of zeolite catalysts can have a profound effect on the U.S. economy. It is estimated that a 1 percent increase in yield of the gasoline fraction per barrel of oil would represent a savings of 22 million barrels of crude oil per year, representing a reduction of $400 million in the United States' balance of payments. Thus any activity that results in improvement in zeolite catalyst performance is of significant scientific and industrial interest. In addition, due to their 'stability,' uniformity, and, within limits, their 'engineerable' structures, zeolites are being tested as potential adsorbents to purify gases and liquids at the parts-per-billion levels needed in today's electronic, biomedical, and biotechnology industries and for the environment. Other exotic applications, such as host materials for quantum-confined semiconductor atomic arrays, are also being investigated. Because of the importance of this class of material, extensive efforts have been made to characterize their structures and to understand their nucleation and growth mechanisms, so as to be able to custom-make zeolites for a desired application. To date, both the nucleation mechanics and chemistry (such as what are the 'key' nutrients) are, as yet, still unknown for many, if not all, systems. The problem is compounded because there is usually a 'gel' phase present that is assumed to control the degree of supersaturation, and this gel undergoes a continuous 'polymerization' type reaction during nucleation and growth. Generally, for structure characterization and diffusion studies, which are useful in evaluating zeolites for improving yield in petroleum refining as well as for many of the proposed new applications (e.g., catalytic membranes, molecular electronics, chemical sensors) large zeolites (greater than 100 to 1000 times normal size) with minimum lattice defects are desired. Presently, the lack of understanding of zeolite nucleation and growth precludes the custom design of zeolites for these or other uses. It was hypothesized that the microgravity levels achieved in an orbiting spacecraft could help to isolate the possible effects of natural convection (which affects defect formation) and minimize sedimentation, which occurs since zeolites are twice as dense as the solution from which they are formed. This was expected to promote larger crystals by allowing growing crystals a longer residence time in a high-concentration nutrient field. Thus it was hypothesized that the microgravity environment of Earth orbit would allow the growth of large, more defect-free zeolite crystals in high yield.

A Sacco, Jr

Enabling Interoperability in Earth System Digital Twins (ESDT): Integrating Observations, Models, and AI for Actionable Insights Through NASA'S Intelligent Systems Technology Program

NASA’s Intelligent Systems Technology Program (IST) is driving a paradigm shift in Earth science through the development of Earth System Digital Twins (ESDT). These integrated information systems create a dynamic "digital replica" of the Earth by harmonizing continuous, multi-source observations with high-fidelity models and state-of-the-art artificial intelligence (AI) that enable “What now?”, “What next?”, and “What if?” scenario building. These scenarios are reflected in NASA IST’s series of ESDTs, from the Coastal Zone Digital Twin that integrates complex data on the current state of the Chesapeake Bay to the Terrestrial Environmental Rapid-Replication and Assimilation Hydrometeorological (TerraHydro) AI-based ESDT that forecasts water movement across Earth’s surface, to the Agriculture Land Information System (AgLIS) which can be used to assess optimal planting dates and crop yield estimates. By bridging the gap between vast data archives and actionable insights, these projects enable a system-of-systems approach to understanding complex, interacting Earth processes. This poster will highlight recent innovations and future directions from NASA’s ESDT initiatives: Continuous Data Assimilation & Multi-Source Fusion. A core requirement of the ESDT work is the transition from static models to dynamic "living" replicas. This involves creating frameworks for the continual assimilation of near-real-time data from uncoordinated, heterogeneous sources, including satellite observations and airborne assets, and ground-based Internet of Things (IoT) sensors. These systems link design, operational status, and environmental data, ensuring the digital twin accurately reflects the current state of the physical Earth system. High-Fidelity Hybrid Modeling & Computational Acceleration to enable interactive "what-if" explorations, programs are moving beyond traditional, slow physical solvers by developing fast surrogate machine learning models and Deep Generative Models (DGMs). These hybrid approaches use neural networks to emulate complex physics, such as cloud feedback or ocean dynamics, at a fraction of the original computing cost, often leveraging advanced hardware like Graphics Processing Units (GPUs) to achieve the necessary scale. Federated Ecosystems & Interoperable Frameworks rather than building isolated tools, NASA IST is moving toward federated ESDTs and reusable analytic collaborative frameworks. This theme focuses on interoperability standards and common ontologies that allow specialized digital twins to interact and share data. This system-of-systems architecture supports multi-discipline investigations, such as analyzing how upstream watershed changes impact downstream urban flooding or how wildfire emissions affect regional air quality. By leveraging these advancements, ESDTs empower researchers and decision-makers to conduct real-time analysis and run complex hypothetical scenarios, ultimately improving our understanding of Earth’s evolving systems and informing critical real-world applications.

Earth System

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan

Lunar Electrical Conductivity, Permeability and Temperature From Apollo Magnetometer Experiments

Magnetometers have been deployed at four Apollo sites on the Moon to measure remanent and induced lunar magnetic fields. Measurements from this network of instruments have been used to calculate the electrical conductivity, temperature, magnetic permeability, and iron abundance of the lunar interior. The measured lunar remanent fields range from 3 gammas (γ) minimum at the Apollo 15 site to 327 γ maximum at the Apollo 16 site. Simultaneous magnetic field and solar plasma pressure measurements show that the remanent fields at the Apollo 12 and 16 sites interact with, and are compressed by, the solar wind. Remanent fields at Apollo 12 and Apollo 16 are increased 16 γ and 32 γ, respectively, by a solar plasma bulk pressure increase of 1.5 x 10 -7 dynes/cm 2 . Global lunar fields due to eddy currents, induced in the lunar interior by magnetic transients, have been analyzed to calculate an electrical conductivity profile for the Moon. From nightside magnetometer data in the solar wind it has been found that deeper than 170 km into the Moon the conductivity rises from 3 x 10 -4 mhos/m to 10 -2 mhos/m at 1000-km depth. Recent analysis of data obtained in the geomagnetic tail, in regions free of complicating plasma effects, yields results consistent with nightside values. Conductivity profiles have been used to calculate the lunar temperature for an assumed lunar material of olivine. In the outer layer (~ 170 km thick) the temperature rises to 110° C, after which it gradually increases with depth to 1500° C at a depth of ~ 1000 km. Simultaneous measurements by magnetometers on the lunar surface and in orbit around the moon are used to construct a whole-moon hysteresis curve, from which the global lunar magnetic permeability is determined to be μ = 1.012 ± 0.006. The corresponding global induced dipole moment is 2 x 10 18 gauss-cm 3 for typical inducing fields of 10 -4 gauss in the lunar environment. Lunar free iron abundance corresponding to the global permeability is determined to be 2.5 ± 2.0 wt. %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior. For a free iron/orthopyroxene lunar composition the total iron content is calculated to be 12.8 ± 1.0 wt. %; for a free iron/olivine composition, total iron content is 5.5 ± 1.2 wt. %. Other lunar models with an iron core and with a shallow iron-rich layer are also discussed in light of the measured global lunar permeability. Velocities and thicknesses of the Earth's magnetopause and bow shock have been estimated from simultaneous magnetometer measurements. Average speeds are determined to be about 50 km/s for the magnetopause and 70 km/s for the bow shock, although there are large variations in the measurement for any particular boundary crossing. Corresponding measured boundary thicknesses average about 2300 km for the magnetopause and 1400 km for the bow shock.

P Dyal

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff

Equivalent System Mass Comparison of ECLSS CO2 Removal Technologies

This project employs equivalent system mass (ESM) analyses to evaluate and compare state-of-the-art spacecraft CO2 removal systems. The ESM methodology converts sizing characteristics such as volume, power, and cooling requirements into a unified metric of mass, allowing for simple quantitative comparison of equipment impacts on a system level. Comparison of six high TRL technologies – the Carbon Dioxide Removal System (CDRS), Thermal Amine Scrubber (TAS), Carbon Dioxide Removal Assembly (CDRA), Four-Bed CO2 Scrubber (4BCO2), CO2 and Humidity Control Swing Bed (CHC), and the Carbon Dioxide Removal by Ionic Liquid Sorbent (CDRILS) – was achieved using an existing ESM spreadsheet tool that formerly ranked such technologies for consideration on Gateway. Fidelity of the tool was increased through revision with the most up-to-date sizing and performance data for each technology. ESM values were then estimated using NASA-supplied Lunar, Martian, and low-Earth orbit habitat infrastructure costs. This paper details the ESM approach used, lists major performance variables for each CO2 removal option, and summarizes conclusions of the technology comparison. The results are intended to act as a guide to assist programs in their selection of CO2 technologies best suited for a specific vehicle or mission.

Madeleine C Oliver

A Comparison of Doppler Lidar Wind Sensors for Earth-Orbit Global Measurement Applications

Now, there are four Doppler lidar configurations which are being promoted for the measurement of tropospheric winds: (1) the coherent CO2 Lidar, operating in the 9 micrometer region using a pulsed, atmospheric pressure CO2 gas discharge laser transmitter, and heterodyne detection; (2) the coherent Neodymium doped YAG or Glass Lidar, operating at 1.06 micrometers, using flashlamp or diode laser optical pumping of the solid state laser medium, and heterodyne detection; (3) the Neodymium doped YAG/Glass Lidar, operating at the doubled frequency (at 530 nm wavelength), again using flashlamp or diode laser pumping of the laser transmitter, and using a high resolution tandem Fabry-Perot filter and direct detection; and (4) the Raman shifted Xenon Chloride Lidar, operating at 350 nm wavelength, using a pulsed, atmospheric pressure XeCl gas discharge laser transmitter at 308 nm, Raman shifted in a high pressure hydrogen cell to 350 nm in order to avoid strong stratospheric ozone absorption, also using a high resolution tandem Fabry-Perot filter and direct detection. Comparisons of these four systems can include many factors and tradeoffs. The major portion of this comparison is devoted to efficiency. Efficiency comparisons are made by estimating the number of transmitted photons required for a single pulse wind velocity estimate of + or - 1 m/s accuracy in the middle troposphere, from an altitude of 800 km, which is assured to be reasonable for a polar orbiting platform.

Robert T Menzies

Equivalent System Mass Comparison of ECLSS CO 2 Removal Technologies

This project employs equivalent system mass (ESM) analyses to evaluate and compare state-of-the-art spacecraft CO 2 removal systems. The ESM methodology converts sizing characteristics such as volume, power, and cooling requirements into a unified metric of mass, allowing for simple quantitative comparison of equipment impacts on a system level. Comparison of six high TRL technologies – the Carbon Dioxide Removal System (CDRS), Thermal Amine Scrubber (TAS), Carbon Dioxide Removal Assembly (CDRA), Four-Bed CO 2 Scrubber (4BCO2), CO 2 and Humidity Control Swing Bed (CHC), and the Carbon Dioxide Removal by Ionic Liquid Sorbent (CDRILS) – was achieved using an existing ESM spreadsheet tool that formerly ranked such technologies for consideration on Gateway. Fidelity of the tool was increased through revision with the most up-to-date sizing and performance data for each technology. ESM values were then estimated using NASA-supplied Lunar, Martian, and low-Earth orbit habitat infrastructure costs. This paper details the ESM approach used, lists major performance variables for each CO 2 removal option, and summarizes conclusions of the technology comparison. The results are intended to act as a guide to assist programs in their selection of CO 2 technologies best suited for a specific vehicle or mission.

Madeleine Oliver

The thermal engine of Venus: Implications from the impact crater distribution

Questions concerning the state of Venus' volcanological activity are currently under investigation. These questions were raised by the imaging radar and altimetric data obtained by the Magellan spacecraft, which began mapping Venus in Sept. 1990. Clues about the answers to these questions are provided by the nature and distribution of impact craters on the surface of Venus. Statistical tests show that the distribution of craters on Venus cannot be distinguished from a completely spatially random population. Further, the majority of craters look relatively pristine; only about 5 percent of the observed craters appear to be embayed or flooded by lavas. The simplest interpretation of these observations is a model in which the impact craters lie on a surface that has been undisturbed by volcanological and tectonic processes since the time the surface was formed. Counting up the number of craters and estimating the rate at which meteoroids strike the surface then gives an age for the venusian surface of 500 million years; in this model the volcanic activity since that time has been essentially nil. This scenario is described as a 'production' model, or perhaps more appropriately, a 'catastrophic' model.

Phillips, Roger J.