Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive 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 253 records · Page 14

Attachment-Line Transition Due to Roughness on a 76 Degree Swept Cylinder at Mach 1.6

Experiments were conducted on a 76 deg. swept cylinder to establish the conditions for the attachment-line transition process with, and without, surface roughness in a low-disturbance ("quiet"), Mach number 1.6 flow. Local flow parameters were estimated from pressure measurements. These were in good agreement with predictions from an Euler code (CFL3D) and a boundary layer code (BL3D). Hot-wires and Schlieren photography were used to determine the state of the boundary layer. It was found that, for a near-adiabatic wall condition and a smooth surface, the attachment-line, boundary-layer remained laminar up to the highest attainable Reynolds number (R-bar of 790). Transition under the influence of trip wires was found to depend on wind-tunnel disturbance levels and the onset conditions have been established. Results suggest that current design practice, which is based upon data from conventional ("noisy") tunnels, may be highly conservative.

Coleman, Colin P.↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

Excitation of interstellar hydrogen chloride

We have computed new rate coefficients for the collisional excitation of HCl by He, in the close-coupled formalism and using an interaction potential determined recently by Willey, Choong, & DeLucia. Results have been obtained for temperatures between 10 K and 300 K. With the use of the infinite order sudden approximation, we have derived approximate expressions of general applicability which may be used to estimate how the rate constant for a transition (J to J prime) is apportioned among the various hyperfine states F prime of the final state J prime. Using these new rate coefficients, we have obtained predictions for the HCl rotational line strengths expected from a dense clump of interstellar gas, as a function of the HCl fractional abundance. Over a wide range of HCl abundances, we have found that the line luminosities are proportional to abundance(exp 2/3), a general result which can be explained using a simple analytical approximation. Our model for the excitation of HCl within a dense molecular cloud core indicates that the J = 1 goes to 0 line strengths measured by Blake, Keene, & Phillips toward the Orion Molecular Cloud (OMC-1) imply a fractional abundance n(HCl)/n(H2) approximately 2 x 10(exp -9), a value which amounts to only approximately 0.3% of the cosmic abundance of chlorine nuclei. Given a fractional abundance of 2 x 10(exp -9), the contribution of HCl emission to the total radiative cooling of a dense clump is small. For Orion, we predict a flux approximately 10(exp -19) W/sq cm for the HCl J = 3 goes to 2 line near 159.8 micrometers, suggesting that the strength of this line could be measured using the Infrared Space Observatory.

Neufild, David A.↗

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗

Autonomous Navigation of a Lunar Relay Using GNSS and Other Measurements

Many of the highest priority destinations at the Moon lack a continuous view of Earth, such as the lunar poles or lunar far side. Exploration of these sites will require spacecraft in cislunar space to relay communications and provide position, navigation, and timing (PNT) services. Accurate knowledge of relay position, velocity, and time is essential to these services. This paper describes a concept for a PNT Instrument being developed for the Lunar Communications Relay and Navigation Systems (LCRNS) Project. The instrument is intended as a payload that would enable autonomous, on-board, real-time navigation and timing using Global Navigation Satellite System (GNSS), optical navigation, and one-way measurements from Earth-based ground stations. Hardware-in-the-loop simulations using flight software are used to realistically characterize performance on hardware platforms with a path to flight. These results provide preliminary validation of the proposed PNT Instrument, demonstrate the benefits of augmenting GNSS with other measurements, and serve as an insightful reference for the design of future lunar missions, including those that will operate within the LunaNet framework of standards. This instrument concept relies on several technologies developed at NASA Goddard Space Flight Center (GSFC). For GNSS observables, the instrument relies on the high-altitude NavCube 3 mini (NC3m) GNSS receiver specifically designed for cislunar applications. The autoNGC system, which consists of flight software and a hardware platform, is responsible for fusing the observables using its extended Kalman filter, the Goddard Enhanced Onboard Navigation System (GEONS). Optical navigation observables are processed within autoNGC (“autonomous Navigation, Guidance, and Control”) using the Goddard Image Analysis & Navigation Tool (GIANT) which is also responsible for simulating high-fidelity images for test and analysis. In addition to describing the PNT Instrument and its components, the paper will present predicted performance based on simulation results. As a baseline, it will present GNSS-only hardware-in-the loop results using a NC3m test unit to process Spirent-simulated GPS signals in a potential lunar relay trajectory: a 12-hour elliptical frozen lunar orbit (ELFO). GEONS then processes the GPS pseudorange and time differenced carrier phase measurements to estimate and propagate the relay state (position, velocity, and time). These results extend previously published work that showed preliminary ELFO performance. Previous work has shown the importance of other measurement types, so additional simulations are performed which augment GNSS with ground station observables and several methods of optical navigation, including celestial navigation, limb-finding (e.g., observations of the lunar horizon), and terrain relative navigation (TRN). TRN involves correlating simulated predicted images of the lunar surface with actual imagery; misalignments of landmarks identified in each image are translated into relay state updates. TRN is valuable as a measurement of the relay’s state relative to the Moon, especially during GNSS outages or after maneuvers. One-way Pseudorange and Doppler measurements from Earth-based ground stations are also simulated. The full set of observables is processed using autoNGC. These simulations make use of autoNGC and NC3m test units, a lab atomic clock, and a pulse-per-second (PPS) generation and distribution system. This combination of subsystems, and the hardware platforms used in this analysis, represents a PNT Instrument that could be flown on a lunar relay. Results from the hardware-in-the-loop simulations presented in this paper provide a preliminary assessment of the achievable navigation performance of this instrument concept. PNT Instrument performance is compared to the GPS-only performance, and a discussion is provided on the apparent merits and challenges of each measurement type.

Ben Ashman↗

Pointing system modeling for performance prediction

The paper considers the modeling and analysis of the JPL Precision Pointing and Tracking System, a gyro stabilized platform pointing system using optical target tracking. Sensor and actuator models are described which realistically represent actual hardware, and a practical controller and state estimator is developed through modern control techniques. The effect of the actuator nonlinearities on the system performance evaluation is discussed along with the modeling of the noise sources in the gyro. Finally, the limitations of implementing the control algorithms on the spacecraft computer, and the resulting tradeoffs are presented.

Boussalis, D.↗

Binary Oscillatory Crossflow Electrophoresis

Electrophoresis has long been recognized as an effective analytic technique for the separation of proteins and other charged species, however attempts at scaling up to accommodate commercial volumes have met with limited success. In this report we describe a novel electrophoretic separation technique - Binary Oscillatory Crossflow Electrophoresis (BOCE). Numerical simulations indicate that the technique has the potential for preparative scale throughputs with high resolution, while simultaneously avoiding many problems common to conventional electrophoresis. The technique utilizes the interaction of an oscillatory electric field and a transverse oscillatory shear flow to create an active binary filter for the separation of charged protein species. An oscillatory electric field is applied across the narrow gap of a rectangular channel inducing a periodic motion of charged protein species. The amplitude of this motion depends on the dimensionless electrophoretic mobility, alpha = E(sub o)mu/(omega)d, where E(sub o) is the amplitude of the electric field oscillations, mu is the dimensional mobility, omega is the angular frequency of oscillation and d is the channel gap width. An oscillatory shear flow is induced along the length of the channel resulting in the separation of species with different mobilities. We present a model that predicts the oscillatory behavior of charged species and allows estimation of both the magnitude of the induced convective velocity and the effective diffusivity as a function of a in infinitely long channels. Numerical results indicate that in addition to the mobility dependence, the steady state behavior of solute species may be strongly affected by oscillating fluid into and out of the active electric field region at the ends of the cell. The effect is most pronounced using time dependent shear flows of the same frequency (cos((omega)t)) flow mode) as the electric field oscillations. Under such conditions, experiments indicate that solute is drawn into the cell from reservoirs at both ends of the cell leading to a large mass build up. As a consequence, any initially induced mass flux will vanish after short times. This effect was not captured by the infinite channel model and hence numerical and experimental results deviated significantly. The revised model including finite cell lengths and reservoir volumes allowed quantitative predictions of the time history of the concentration profile throughout the system. This latter model accurately describes the fluxes observed for both oscillatory flow modes in experiments using single protein species. Based on the results obtained from research funded under NASA grant NAG-8-1080.S, we conclude that binary separations are not possible using purely oscillatory flow modes because of end effects associated with the cos((omega)t) mode. Our research shows, however, that a combination of cos(2(omega)t) and steady flow should lead to efficient separation free of end effects. This possibility is currently under investigation.

Molloy, Richard F.↗

Use of Remotely-Sensed Rainfall Data in Global Modeling and Data Assimilation

Precipitation observations derived from microwave sensors available from the Tropical Rainfall Measuring Mission (TRMM) and the proposed Global Precipitation Mission (GPM) can provide crucial information needed for improving global modeling, data assimilation, and numerical weather prediction. New methodologies are being developed at NASA to make effective use of this new data type in these applications. Currently, global analyses contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. We show that assimilating 6-h averaged TRMM rainfall retrievals improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the upper tropospheric moisture in the analysis produced by the Goddard Earth Observing System (GEOS) Data Assimilation System. The improved analysis also leads to improved short-range forecasts in the tropics. The above results were obtained using a variational assimilation procedure that uses rainfall observations to derive moisture and temperature tendency corrections every 6 hours to compensate for errors arising from imperfect initial conditions and deficiencies in the model physics. We will describe a developmental path towards using space-borne rainfall data to empirically estimate and correct for state-dependent systematic errors in parameterized model physics. The study provides a demonstration of the potential of using remote-sensed rainfall data from microwave instruments to improve the 4-dimensional global datasets for climate analysis and numerical weather prediction.

Hou, Arthur↗

Flight test evaluation of predicted light aircraft drag, performance, and stability

A technique was developed which permits simultaneous extraction of complete lift, drag, and thrust power curves from time histories of a single aircraft maneuver such as a pullup (from V sub max to V sub stall) and pushover (to sub V max for level flight.) The technique is an extension to non-linear equations of motion of the parameter identification methods of lliff and Taylor and includes provisions for internal data compatibility improvement as well. The technique was show to be capable of correcting random errors in the most sensitive data channel and yielding highly accurate results. This technique was applied to flight data taken on the ATLIT aircraft. The drag and power values obtained from the initial least squares estimate are about 15% less than the 'true' values. If one takes into account the rather dirty wing and fuselage existing at the time of the tests, however, the predictions are reasonably accurate. The steady state lift measurements agree well with the extracted values only for small values of alpha. The predicted value of the lift at alpha = 0 is about 33% below that found in steady state tests while the predicted lift slope is 13% below the steady state value.

Smetana, F. O.↗

Transport properties of monatomic carbon. II - Contributions from excited electronic states

The aim of the present analysis is to verify the conventional assumption that the transport properties of atoms and molecules in excited electronic states are nearly the same as in the ground state. Reasonably exact calculations of the transport properties of excited-state carbon atoms are obtained on the basis of available experimental and theoretical data on the interactions between carbon atoms in excited electronic states. A comparison with results obtained with several models that have been used to estimate the transport properties of excited species indicates that the assumption under consideration is correct, but that the models do not predict very accurately the transport properties of excited species.

Biolsi, L.↗

High temperature furnace modeling and performance verifications

Analytical, numerical, and experimental studies were performed on two classes of high temperature materials processing sources for their potential use as directional solidification furnaces. The research concentrated on a commercially available high temperature furnace using a zirconia ceramic tube as the heating element and an Arc Furnace based on a tube welder. The first objective was to assemble the zirconia furnace and construct parts needed to successfully perform experiments. The 2nd objective was to evaluate the zirconia furnace performance as a directional solidification furnace element. The 3rd objective was to establish a data base on materials used in the furnace construction, with particular emphasis on emissivities, transmissivities, and absorptivities as functions of wavelength and temperature. A 1-D and 2-D spectral radiation heat transfer model was developed for comparison with standard modeling techniques, and were used to predict wall and crucible temperatures. The 4th objective addressed the development of a SINDA model for the Arc Furnace and was used to design sample holders and to estimate cooling media temperatures for the steady state operation of the furnace. And, the 5th objective addressed the initial performance evaluation of the Arc Furnace and associated equipment for directional solidification. Results of these objectives are presented.

Smith, James E., Jr.↗

Scenario for Hollow Cathode End-Of-Life

Recent successful hollow cathode life tests have demonstrated that lifetimes can meet the requirements of several space applications. However, there are no methods for assessing cathode lifetime short of demonstrating the requirement. Previous attempts to estimate or predict cathode lifetime were based on relatively simple chemical depletion models derived from the dispenser cathode community. To address this lack of predicative capability, a scenario for hollow cathode lifetime under steady-state operating conditions is proposed. This scenario has been derived primarily from the operating behavior and post-test condition of a hollow cathode that was operated for 28,000 hours. In this scenario, the insert chemistry evolves through three relatively distinct phases over the course of the cathode lifetime. These phases are believed to correspond to demonstrable changes in cathode operation. The implications for cathode lifetime limits resulting from this scenario are examined, including methods to assess cathode lifetime without operating to End-of- Life and methods to extend the cathode lifetime.

Sarver-Verhey, Timothy R.↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

Soil Moisture Data Assimilation

Accurate knowledge of soil moisture at the continental scale is important for improving predictions of weather, agricultural productivity and natural hazards, but observations of soil moisture at such scales are limited to indirect measurements, either obtained through satellite remote sensing or from meteorological networks. Land surface models simulate soil moisture processes, using observation-based meteorological forcing data, and auxiliary information about soil, terrain and vegetation characteristics. Enhanced estimates of soil moisture and other land surface variables, along with their uncertainty, can be obtained by assimilating observations of soil moisture into land surface models. These assimilation results are of direct relevance for the initialization of hydro-meteorological ensemble forecasting systems. The success of the assimilation depends on the choice of the assimilation technique, the nature of the model and the assimilated observations, and, most importantly, the characterization of model and observation error. Systematic differences between satellite-based microwave observations or satellite-retrieved soil moisture and their simulated counterparts require special attention. Other challenges include inferring root-zone soil moisture information from observations that pertain to a shallow surface soil layer, propagating information to unobserved areas and downscaling of coarse information to finer-scale soil moisture estimates. This chapter summarizes state-of-the-art solutions to these issues with conceptual data assimilation examples, using techniques ranging from simplified optimal interpolation to spatial ensemble Kalman filtering. In addition, operational soil moisture assimilation systems are discussed that support numerical weather prediction at ECMWF and provide value-added soil moisture products for the NASA Soil Moisture Active Passive mission.

radar backscatter↗

Nonequilibrium Stagnation-Line Radiative Heating for Fire II

This paper presents a detailed analysis of the shock-layer radiative heating to the Fire II vehicle using a new air radiation model and a viscous shock-layer flowfield model. This new air radiation model contains the most up-to-date properties for modeling the atomic-line, atomic photoionization, molecular band, and non-Boltzmann processes. The applied viscous shock-layer flowfield analysis contains the same thermophysical properties and nonequilibrium models as the LAURA Navier-Stokes code. Radiation-flowfield coupling, or radiation cooling, is accounted for in detail in this study. It is shown to reduce the radiative heating by about 30% for the peak radiative heating points, while reducing the convective heating only slightly. A detailed review of past Fire II radiative heating studies is presented. It is observed that the scatter in the radiation predicted by these past studies is mostly a result of the different flowfield chemistry models and the treatment of the electronic state populations. The present predictions provide, on average throughout the trajectory, a better comparison with Fire II flight data than any previous study. The magnitude of the vacuum ultraviolet (VUV) contribution to the radiative flux is estimated from the calorimeter measurements. This is achieved using the radiometer measurements and the predicted convective heating. The VUV radiation predicted by the present model agrees well with the VUV contribution inferred from the Fire II calorimeter measurement, although only when radiation-flowfield coupling is accounted for. This agreement provides evidence that the present model accurately models the VUV radiation, which is shown to contribute significantly to the Fire II radiative heating.

Johnston, Christopher O.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

A Coupled-Cluster Study of XON (X=H, F, Cl), and the XON (left and right arrow) XNO Transition States

The XON molecules (X=H, F, and Cl) have been studied using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with a double polarized triple-zeta basis set. The equilibrium geometries, dipole moments, harmonic vibrational frequencies and infrared intensities have been predicted. The X-O bond distance is shown to be abnormally long for X=F and Cl, and this is attributed to the degree of ionic bonding and the stability of NO(+). Based on Mulliken population analyses, it is shown that there is a significant degree of X(-), ON(+) ionic bonding character for FON and ClON, whereas for HON the ionic character is reduced and best described as H(+), NO(-). The stability of the XON molecule relative to the XNO isomer is shown to increase in the order HON<FON<ClON, although even ClON is 27.0 plus or minus 1.0 kcal/mol (0 K) less stable than ClNO. The XON (left and right arrow) XNO transition states are also investigated. FON and ClON possess the lowest barrier heights to isomerization (7.2 plus or minus 1.0 kcal/mol at 0 K), but these are sufficiently large to suggest that isomerization should not be rapid at low temperatures.

Lee, Timothy↗

Determination of Eros Physical Parameters for Near Earth Asteroid Rendezvous Orbit Phase Navigation

Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Miller, J. K.↗