Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Uncertainty Propagation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Atmospheric drag model calibrations for spacecraft lifetime prediction

Although solar activity prediction uncertainty normally dominates decay prediction error budget for near-Earth spacecraft, the effect of drag force modeling errors for given levels of solar activity needs to be considered. Two atmospheric density models, the modified Harris-Priester model and the Jacchia-Roberts model, to reproduce the decay histories of the Solar Mesosphere Explorer (SME) and Solar Maximum Mission (SMM) spacecraft in the 490- to 540-kilometer altitude range were analyzed. Historical solar activity data were used in the input to the density computations. For each spacecraft and atmospheric model, a drag scaling adjustment factor was determined for a high-solar-activity year, such that the observed annual decay in the mean semimajor axis was reproduced by an averaged variation-of-parameters (VOP) orbit propagation. The SME (SMM) calibration was performed using calendar year 1983 (1982). The resulting calibration factors differ by 20 to 40 percent from the predictions of the prelaunch ballistic coefficients. The orbit propagations for each spacecraft were extended to the middle of 1988 using the calibrated drag models. For the Jaccia-Roberts density model, the observed decay in the mean semimajor axis of SME (SMM) over the 4.5-year (5.5-year) predictive period was reproduced to within 1.5 (4.4) percent. The corresponding figure for the Harris-Priester model was 8.6 (20.6) percent. Detailed results and conclusions regarding the importance of accurate drag force modeling for lifetime predictions are presented.

Binebrink, A. L.↗
NASA NTRS2019

Forecasting the Ambient Solar Wind with Numerical Models. I. On the Implementation of an Operational Framework

The ambient solar wind conditions in interplanetary space and in the near-Earth environment are determined by activity on the Sun. Steady solar wind streams modulate the propagation behavior of interplanetary coronal mass ejections and are themselves an important driver of recurrent geomagnetic storm activity. The knowledge of the ambient solar wind flows and fields is thus an essential component of successful space weather forecasting. Here, we present an implementation of an operational framework for operating, validating, and optimizing models of the ambient solar wind flow on the example of Carrington Rotation 2077. We reconstruct the global topology of the coronal magnetic field using the potential field source surface model (PFSS) and the Schatten current sheet model (SCS) and discuss three empirical relationships for specifying the solar wind conditions near the Sun, namely the Wang–Sheeley (WS) model, the distance from the coronal hole boundary model (DCHB), and the Wang-Sheeley–Arge (WSA) model. By adding uncertainty in the latitude about the sub-Earth point, we select an ensemble of initial conditions and map the solutions to Earth by the Heliospheric Upwind eXtrapolation (HUX) model. We assess the forecasting performance from a continuous variable validation and find that the WSA model most accurately predicts the solar wind speed time series (RMSE≈83 km/s). We note that the process of ensemble forecasting slightly improves the forecasting performance of all solar wind models investigated. We conclude that the implemented framework is well suited for studying the relationship between coronal magnetic fields and the properties of the ambient solar wind flow in the near-Earth environment.

Reiss, Martin A.↗
NASA NTRS2009

Probabilistic Estimation of Critical Flaw Sizes in the Primary Structure Welds of the Ares I-X Launch Vehicle

The primary structure of the Ares I-X Upper Stage Simulator (USS) launch vehicle is constructed of welded mild steel plates. There is some concern over the possibility of structural failure due to welding flaws. It was considered critical to quantify the impact of uncertainties in residual stress, material porosity, applied loads, and material and crack growth properties on the reliability of the welds during its pre-flight and flight. A criterion--an existing maximum size crack at the weld toe must be smaller than the maximum allowable flaw size--was established to estimate the reliability of the welds. A spectrum of maximum allowable flaw sizes was developed for different possible combinations of all of the above listed variables by performing probabilistic crack growth analyses using the ANSYS finite element analysis code in conjunction with the NASGRO crack growth code. Two alternative methods were used to account for residual stresses: (1) The mean residual stress was assumed to be 41 ksi and a limit was set on the net section flow stress during crack propagation. The critical flaw size was determined by parametrically increasing the initial flaw size and detecting if this limit was exceeded during four complete flight cycles, and (2) The mean residual stress was assumed to be 49.6 ksi (the parent material s yield strength) and the net section flow stress limit was ignored. The critical flaw size was determined by parametrically increasing the initial flaw size and detecting if catastrophic crack growth occurred during four complete flight cycles. Both surface-crack models and through-crack models were utilized to characterize cracks in the weld toe.

Pai, Shantaram S.↗
NASA NTRS2004

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Method and Uncertainties - Part 1

A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating/drying profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and non-convective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud resolving model simulations, and from the Bayesian formulation itself. Synthetic rain rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in instantaneous rain rate estimates at 0.5 deg resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. These errors represent about 70-90% of the mean random deviation between collocated passive microwave and spaceborne radar rain rate estimates. The cumulative algorithm error in TMI estimates at monthly, 2.5 deg resolution is relatively small (less than 6% at 5 mm/day) compared to the random error due to infrequent satellite temporal sampling (8-35% at the same rain rate).

Olson, William S.↗
NASA NTRS2006

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Improved Method and Uncertainties - Part 1

A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and nonconvective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud-resolving model simulations, and from the Bayesian formulation itself. Synthetic rain-rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in TMI instantaneous rain-rate estimates at 0.5 -resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. Errors in collocated spaceborne radar rain-rate estimates are roughly 50%-80% of the TMI errors at this resolution. The estimated algorithm random error in TMI rain rates at monthly, 2.5deg resolution is relatively small (less than 6% at 5 mm day.1) in comparison with the random error resulting from infrequent satellite temporal sampling (8%-35% at the same rain rate). Percentage errors resulting from sampling decrease with increasing rain rate, and sampling errors in latent heating rates follow the same trend. Averaging over 3 months reduces sampling errors in rain rates to 6%-15% at 5 mm day.1, with proportionate reductions in latent heating sampling errors.

Olson, William S.↗
NASA NTRS1995

ROSAT observations of the gravitationally lensed system 0957+561

Observations of 0957+561 with the Einstein High Resolution Imager (HRI) in 1979 May and 1979 November and with the ROSAT HRI in 1991 May and 1992 October reveal large variations in the X-ray fluxes for images A and B which significantly differ from simultaneously observed changes of the corresponding optical continuum emission. Most striking was the apparent increase by over fivefold in the flux of image B from the late 1970's to the early 1900's. Further, in the 1990's the X-ray flux for image A increased by a factor of 1.7 +/- 0.3 and that of B by a factor of 1.9 +/- 0.2 between the two ROSAT observations (separated by 540 days), whereas optical measurements showed nearly no changes for A and B between these epochs. No significant changes in X-ray emission were observed on timescales of hours to days. Of we adopt 1.5 yr as the difference between the propagation time from the quasar to us via the two images ('time delay'), then the ratio of the X-ray flux of image B for the 1992 October epoch to that of A for the 1991 May epoch is 2.7 +/- 0.4. This ratio is significantly greater than the ratio of 1.05 +/- 0.03 observed in the optical R band, 0.75 in the broad line region (BLR), and 0.76 +/- 0.03 in the radio (VLBI lambda 13 cm core) for approximately the same epoch. The wavelength dependence in the ratio of the fluxes of the two images suggests that either microlensing may be significant for the X-ray band or the time delay is substanially different from 1.5 yr (and the intrinsic variation of the quasar emission were significant within an interval comparable to the uncertainty of the time delay), or both.

Chartas, G.↗
NASA NTRS2002

Algorithms for Autonomous GPS Orbit Determination and Formation Flying: Investigation of Initialization Approaches and Orbit Determination for HEO

This report summarizes the efforts to date in processing GPS measurements in High Earth Orbit (HEO) applications by the Colorado Center for Astrodynamics Research (CCAR). Two specific projects were conducted; initialization of the orbit propagation software, GEODE, using nominal orbital elements for the IMEX orbit, and processing of actual and simulated GPS data from the AMSAT satellite using a Doppler-only batch filter. CCAR has investigated a number of approaches for initialization of the GEODE orbit estimator with little a priori information. This document describes a batch solution approach that uses pseudorange or Doppler measurements collected over an orbital arc to compute an epoch state estimate. The algorithm is based on limited orbital element knowledge from which a coarse estimate of satellite position and velocity can be determined and used to initialize GEODE. This algorithm assumes knowledge of nominal orbital elements, (a, e, i, omega, omega) and uses a search on time of perigee passage (tau(sub p)) to estimate the host satellite position within the orbit and the approximate receiver clock bias. Results of the method are shown for a simulation including large orbital uncertainties and measurement errors. In addition, CCAR has attempted to process GPS data from the AMSAT satellite to obtain an initial estimation of the orbit. Limited GPS data have been received to date, with few satellites tracked and no computed point solutions. Unknown variables in the received data have made computations of a precise orbit using the recovered pseudorange difficult. This document describes the Doppler-only batch approach used to compute the AMSAT orbit. Both actual flight data from AMSAT, and simulated data generated using the Satellite Tool Kit and Goddard Space Flight Center's Flight Simulator, were processed. Results for each case and conclusion are presented.

Axelrad, Penina↗
NASA NTRS1994

Measuring rainwater content by radar using propagation differential phase shift

While radars measure several quantities closely coupled to the rainfall rate, for frequencies less than 15 GHz, estimates of the rainwater content W are traditionally computed from the radar reflectivity factor Z or the rate of attenuation A--quantities only weakly related to W. Consequently, instantaneous point estimates of W using Z and A are often erroneous. A more natural, alternative parameter for estimating W at these frequencies is the specific polarization propagation differential phase shift phi(sub DP), which is a measure of the change in the difference between phases of vertically (V) and horizontally (H) polarized waves with increasing distance from a radar. It is now well known that W is nearly linearly related to phi(sub DP) divided by (1 - reversed R), where reversed R is the mass-weighted mean axis ratio of the raindrops. Unfortunately, such relations are not widely used in part because measurements of phi(sub DP) are scarce but also because one must determine reversed R. In this work it is shown that this parameter can be estimated using the differential reflectivity (Z(sub H)/Z(sub V) at 3 GHz. An alternative technique is suggested for higher frequencies when the differential reflectivity becomes degraded by attenuation. While theory indicates that it should be possible using phi(sub DP) to estimate W quite accurately, measurement errors increase the uncertainty to +/- 18%-35% depending on reversed R. While far from ideal, it appears that these estimates are likely to be considerably more accurate than those deduced using currently available methods.

Jameson, A. R.↗
NASA NTRS1997

The Wardle Instability in Interstellar Shocks: Gas Temperture and Line Emission - 2

We have modeled the gas temperature structure in unstable C-type shocks and obtained predictions for the resultant CO and H2 rotational line emissions, using numerical simulations of the Wardle instability. Our model for the thermal balance of the gas includes ion-neutral frictional heating; compressional heating; radiative cooling due to rotational and ro-vibrational transitions of the molecules CO, H2O, and H2; and gas-grain collisional cooling. We obtained results for the gas temperature distribution in-and H2 and CO line emission from-shocks of neutral Alfvenic Mach number 10 and velocity 20 or 40 km/ s in which the Wardle instability has saturated. Both two- and three-dimensional simulations were carried out for shocks in which the preshock magnetic field is perpendicular to the shock propagation direction, and a two-dimensional simulation was carried out for the case in which the magnetic field is obliquely oriented with respect to the shock propagation direction. Although the Wardle instability profoundly affects the density structure behind C-type shocks, most of the shock-excited molecular line emission is generated upstream of the region where the strongest effects of the instability are felt. Thus the Wardle instability has a relatively small effect on the overall gas temperature distribution in-and the emission-line spectrum from-C-type shocks, at least for the cases that we have considered. In none of the cases that we have considered thus far did any of the predicted emission-line luminosities change by more than a factor of 2.5, and in most cases the effects of instability were significantly smaller than that. Slightly larger changes in the line luminosities seem likely for three-dimensional simulations of oblique shocks, although such simulations have yet to be carried out and lie beyond the scope of this study. Given the typical uncertainties that are always present when model predictions are compared with real astronomical data, we conclude that Wardle instability does not imprint any clear observational signature on the shock-excited CO and H2 line strengths. This result justifies the use of one-dimensional steady shock models in the interpretation of observations of shock-excited line emission in regions of star formation. Our three-dimensional simulations of perpendicular shocks revealed the presence of warm filamentary structures that are aligned along the magnetic field, a result that is of possible relevance to models of water maser emission from C-type shocks.

Neufeld, David A.↗
NASA NTRSDate not supplied

Early Information Parameter-Set Analysis for Satellite Close Approaches using Machine Learning

Understanding orbital mechanics is essential in space flight and navigation applications, and leveraging modern force models for flight path projection remains an important aspect in space mission design and operation. However, force models do not capture all the dynamics or perturbations in the space environment and thus are subject to errors in predicting the state vectors. The further out the predicted miss distance between spacecraft is from the time of closest approach (TCA), the larger the propagated errors in the predicted miss distance at TCA is. The dependency on these force models for spacecraft flight state prediction calls for a more reliable method that can quantify, or even reduce, these propagated errors. With recent advances in the field artificial intelligence, specifically in machine and deep learning algorithms, a model that implements these approaches can improve on the modern force model approach. The goal for this work is to provide an early-information decision-making threshold, in order to prioritize risk assessment implementation, given the ongoing increase of space objects. In analyzing the relationship of several parameters from conjunction data messages(CDMs) and solar information, early information becomes viable in miss distance prediction with unsupervised learning techniques, which learn the parameters that are linked together with miss distance and probability of collision (Pc) variables. Another approach implemented for identifying relationships within CDMs is supervised learning, in which a shallow neural network binary classifier learns to distinguish events with Pc values¡108. These parameters detected in the unsupervised process are then applied to a regression neural network, which predicts the miss distance at TCA for a specific event within a given uncertainty bound. For the regression neural network, a Long Short Term Memory (LSTM) neural network is implemented, which yields memory about each time step in an event. Using an LSTM network, the model learns to predict miss distance within 0.2km of the value measured at TCA. Although there is a limited amount of "close miss" data to train a network, the network learns to associate parameters, like large energy dissipation rates with the secondary object, with an elevated Pc

Brianna I. Robertson↗
NASA NTRS1988

Model-scale sound propagation experiment

The results of a scale model propagation experiment to investigate grazing propagation above a finite impedance boundary are reported. In the experiment, a 20 x 25 ft ground plane was installed in an anechoic chamber. Propagation tests were performed over the plywood surface of the ground plane and with the ground plane covered with felt, styrofoam, and fiberboard. Tests were performed with discrete tones in the frequency range of 10 to 15 kHz. The acoustic source and microphones varied in height above the test surface from flush to 6 in. Microphones were located in a linear array up to 18 ft from the source. A preliminary experiment using the same ground plane, but only testing the plywood and felt surfaces was performed. The results of this first experiment were encouraging, but data variability and repeatability were poor, particularly, for the felt surface, making comparisons with theoretical predictions difficult. In the main experiment the sound source, microphones, microphone positioning, data acquisition, quality of the anechoic chamber, and environmental control of the anechoic chamber were improved. High-quality, repeatable acoustic data were measured in the main experiment for all four test surfaces. Comparisons with predictions are good, but limited by uncertainties of the impedance values of the test surfaces.

Willshire, William L., Jr.↗
NASA NTRS2015

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation that is based on a coupling between inertial navigation and flush air data sensing pressure measurements. In this approach, the full navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to directly estimate atmospheric winds and density using a nonlinear weighted least-squares algorithm. The approach uses a high fidelity model of atmosphere stored in table-look-up form, along with simplified models of that are propagated along the trajectory within the algorithm to provide prior estimates and covariances to aid the air data state solution. Thus, the method is essentially a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere and winds are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the discrete-time observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content to the system. The algorithm is then applied to the design of the pressure measurement system for the Mars 2020 mission. The pressure port layout is optimized to maximize the observability of atmospheric states along the trajectory. Linear covariance analysis is performed to assess estimator performance for a given pressure measurement uncertainty. The results indicate that the new tightly-coupled estimator can produce enhanced estimates of atmospheric states when compared with existing algorithms.

Karlgaard, Christopher D.↗
NASA NTRS1986

Detonation wave compression in gas turbines

A study was made of the concept of augmenting the performance of low pressure ratio gas turbines by detonation wave compression of part of the flow. The concept exploits the constant volume heat release of detonation waves to increase the efficiency of the Brayton cycle. In the models studied, a fraction of the compressor output was channeled into detonation ducts where it was processed by transient transverse detonation waves. Gas dynamic studies determined the maximum cycling frequency of detonation ducts, proved that upstream propagation of pressure pulses represented no problems and determined the variations of detonation duct output with time. Mixing and wave compression were used to recombine the combustor and detonation duct flows and a concept for a spiral collector to further smooth the pressure and temperature pulses was presented as an optional component. The best performance was obtained with a single firing of the ducts so that the flow could be re-established before the next detonation was initiated. At the optimum conditions of maximum frequency of the detonation ducts, the gas turbine efficiency was found to be 45 percent while that of a corresponding pressure ratio 5 conventional gas turbine was only 26%. Comparable improvements in specific fuel consumption data were found for gas turbines operating as jet engines, turbofans, and shaft output machines. Direct use of the detonation duct output for jet propulsion proved unsatisfactory. Careful analysis of the models of the fluid flow phenomena led to the conclusion that even more elaborate calculations would not diminish the uncertainties in the analysis of the system. Feasibility of the concept to work as an engine now requires validation in an engineering laboratory experiment.

Wortman, A.↗
NASA NTRS2012

A Geometric Analysis to Protect Manned Assets from Newly Launched Objects - COLA Gap Analysis

A safety risk was identified for the International Space Station (ISS) by The Aerospace Corporation following the launch of GPS IIR-20 (March 24, 2009), when the spent upper stage of the launch vehicle unexpectedly crossed inside the ISS notification box shortly after launch. This event highlighted a 56-hour vulnerability period following the end of the launch Collision Avoidance (COLA) process where the ISS would be unable to react to a conjunction with a newly launched object. Current launch COLA processes screen each launched object across the launch window to determine if an object's nominal trajectory is predicted to pass within 200 km of the ISS (or any other manned/mannable object), resulting in a launch time closure. These launch COLA screens are performed from launch through separation plus I 00 minutes. Once the objects are in orbit, they are cataloged and evaluated as part of routine on-orbit conjunction assessment processes. However, as the GPS IIR-20 scenario illustrated, there is a vulnerability period in the time line between the end of launch COLA coverage and the beginning of standard on-orbit COLA assessment activities. The gap between existing launch and on-orbit COLA processes is driven by the time it takes to track and catalog a launched object, identify a conjunction, and plan and execute a collision avoidance maneuver. For the ISS, the total time required to accomplish an of these steps is 56 hours. To protect human lives, NASA/JSC has requested that an US launches take additional steps to protect the ISS during this "COLA gap" period. The uncertainty in the state of a spent upper stage can be quite large after all bums are complete and all remaining propellants are expelled to safe the stage. Simply extending the launch COLA process an additional 56 hours is not a viable option as the 3-sigma position uncertainty will far exceed the 200 km miss-distance criterion. Additionally, performing a probability of collision (Pc) analysis over this period is also not practical due to the limiting effects of these large orbit state uncertainties. An estimated upper bound for Pc for a typical spent upper stage if nominally aligned for a direct broadside collision with the ISS is only on the order of 10-6. For a smaller manned object such as a Soyuz capsule, the risk level decreases to an order of 10'8 . In comparison, the Air Force Range policy (AFI 91-217) for launch COLAs would only eliminate launch opportunities when conjunctions with objects exceed a Pc of 10'5 This paper demonstrates a conservative geometry-based methodology that may be used to determine if launch opportunities pose a threat to the ISS during the COLA gap period. The NASA Launch Services Program at Kennedy Space Center has developed this COLA gap analysis method and employed it fQr three NASA missions to identify potential ISS conjunctions and corresponding launch window closures during the 56-hour at-risk period. In the analysis, for each launch opportunity, the nominal trajectory of the spent upper stage and the orbit state of the ISS are propagated over the 56 hour period. Each time the upper stage crosses the orbit plane of the ISS, the relative radial and argument of latitude separations are calculated. A window cutout is identified if these separation differences fall within a mission-specific violation box, which is determined from the evaluation of a Monte Carlo dispersions analysis that quantifies the potential variation in the upper stage radial and argument of latitude differences. This paper details the results of these analyses and their impacts to each mission.

Hametz, Mark E.↗
NASA NTRSDate not supplied

Navigation Performance of the BioSentinel Deep Space CubeSat Mission

The BioSentinel mission was recently launched aboard the SLS launch vehicle (LV) as part of the Artemis- 1 campaign. The BioSentinel navigation team successfully tracked and guided the spacecraft through a lunar gravity assist to its destination Earth-trailing heliocentric orbit. This 6U CubeSat carries live yeast cells to analyze the effects of radiation at large distances from Earth, becoming the first biological payload in Deep Space. Prelaunch activities included mission design updates, orbit determination rehearsals and the development of a tracking schedule in coordination with the Artemis-1 payload office and the Deep Space Network (DSN). An important influence on the trajectories of Artemis I secondaries was the uncertainty associated with deployment from the Interim Cryogenic Propulsion System (ICPS), the upper stage of the SLS LV. The ICPS was rotating at a rate of 1 rpm; there was also an uncertainty in the spin axis attitude, which translated into an unknown clock angle of deployment. The variability in this angle and magnitude of deployment implied the existence of a non-negligible risk of a lunar impact, which was evaluated for various potential launch dates. We present the results of Monte Carlo analyses and compute the pertinent maneuvers to avoid it. In addition, we present a comparison with the actual deployment once the mission launched by reconstructing our trajectory with tracking data. On November 16th 2022 BioSentinel successfully deployed from ICPS and the navigation team started to receive 2-way Doppler and Sequential Ranging data from the DSN. We processed early data to try to obtain a first ephemeris using Initial Orbit Determination (IOD) methods such as the least squares. Soon after deployment, the spacecraft was tumbling and entered safe mode, creating a period where the tracking data were sparse. The mission team recovered the spacecraft and after four tracking passes, we solved for a first ephemeris that was sent to the DSN for better tracking of the spacecraft. After propagating this first ephemeris solution, we determined that we avoided impact with a margin of a few hundred km from the lunar surface. More tracking data over the next few days (from DSN as well as ESA antennas) allowed for a more refined orbit solution predicting a periselene altitude of 406 km and a lunar eclipse lasting 36.5 minutes. Therefore, BioSentinel operators aborted any correction maneuvers. This periselene altitude also gave us the necessary energy to achieve a heliocentric orbit. The next challenge was due to the necessary adjustments in our orbit determination method due to the large energy boost resulting from the lunar flyby. After a series of tracking passes we were able to get a nominal solution that resulted into a stable trajectory. This paper discusses in detail the navigation performance using the X-band IRIS transponder, as well as the challenges and lessons learned prior to and during this deep space, CubeSat mission.

Andres Dono Perez↗
NASA NTRS2000

Burning in Outer Space: Microgravity

A better understanding of combustion can lead to significant technological advances, such as less polluting, more fuel-efficient vehicles. Unfortunately, gravity can interfere with the study of combustion. Gravity drags down gases that are cooler- and, therefore, denser-than heated gases. This movement mixes the fuel and the oxidizer substance that promotes burning. Because of this mixing, an observer cannot necessarily distinguish what is happening as a result of the natural combustion process and what is caused, by the pull of gravity. To remove this uncertainty, scientists can conduct experiments that simulate the negation of gravity through freefall. This condition is known as a microgravity environment. A micro-gravity experiment may take place in a chamber that is dropped down a hole or from a high-speed drop tower. The experiment also be conducted in an airplane or a rocket during freefall in a parabolic flight path. This method provides less than a minute of microgravity at most. An experiment that requires the prolonged absence of gravity may necessitate the use of an orbiting spacecraft as a venue. However, access to an orbital laboratory is difficult to acquire. High-end computing centers such as the NCCS can provide a practical alternative to operating in microgravity. Scientists can model phenomena such as combustion without gravitys observational interference. The study of microgravity combustion produces important benefits beyond increased observational accuracy. Certain valuable materials that are produced through combustion can be formed with a more uniform crystal structure-and, therefore, improved structural quality-when the pull of gravity is removed. Furthermore, understanding how fires propagate in the absence of gravity can improve fire safety aboard spacecraft.

Matkowsky, Bernard↗
NASA NTRS2002

Pulse Detonation Rocket MHD Power Experiment

A pulse detonation research engine (MSFC (Marshall Space Flight Center) Model PDRE (Pulse Detonation Rocket Engine) G-2) has been developed for the purpose of examining integrated propulsion and magnetohydrodynamic power generation applications. The engine is based on a rectangular cross-section tube coupled to a converging-diverging nozzle, which is in turn attached to a segmented Faraday channel. As part of the shakedown testing activity, the pressure wave was interrogated along the length of the engine while running on hydrogen/oxygen propellants. Rapid transition to detonation wave propagation was insured through the use of a short Schelkin spiral near the head of the engine. The measured detonation wave velocities were in excess of 2500 m/s in agreement with the theoretical C-J velocity. The engine was first tested in a straight tube configuration without a nozzle, and the time resolved thrust was measured simultaneously with the head-end pressure. Similar measurements were made with the converging-diverging nozzle attached. The time correlation of the thrust and head-end pressure data was found to be excellent. The major purpose of the converging-diverging nozzle was to configure the engine for driving an MHD generator for the direct production of electrical power. Additional tests were therefore necessary in which seed (cesium-hydroxide dissolved in methanol) was directly injected into the engine as a spray. The exhaust plume was then interrogated with a microwave interferometer in an attempt to characterize the plasma conditions, and emission spectroscopy measurements were also acquired. Data reduction efforts indicate that the plasma exhaust is very highly ionized, although there is some uncertainty at this time as to the relative abundance of negative OH ions. The emission spectroscopy data provided some indication of the species in the exhaust as well as a measurement of temperature. A 24-electrode-pair segmented Faraday channel and 0.6 Tesla permanent magnet assembly were then installed on Marshall Space Flight Center's (MSFC's) rectangular channel pulse detonation research engine. Magnetohydrodynamic (MHD) electrical power extraction experiments were carried out for a range of load impedances in which cesium hydroxide seed (dissolved in methanol) was sprayed into the gaseous oxygen/hydrogen propellants. Positive power extraction was obtained, but preliminary analysis of the data indicated that the plasma electrical conductivity is lower than anticipated and the near-electrode voltage drop is not negligible. It is believed that the electrical conductivity is reduced due to a large population of negative OH ions. This occurs because OH has a strong affinity for capturing free electrons. The effect of near-electrode voltage drop is associated with the high surface-to-volume ratio of the channel (1-inch by 1-inch cross-section) where surface effects play a dominant role. As usual for MHD devices, higher performance will require larger scale devices. Overall, the gathered data is extremely valuable from the standpoint of understanding plasma behavior and for developing empirical scaling laws.

Litchford, Ron J.↗
NASA NTRS1992

Multiple scattering lidar returns from stratus clouds

Multiple scattering lidar returns from stratus clouds were measured using a multi-field-of-view (MFOV) lidar operating at 1.054 microns. The detector consists of four concentric silicon photodiodes which define half-angle fields of view (FOV) of 3.8, 12.5, 25 and 38.5 mrad. The central FOV receives the conventional lidar signal, while the outer FOV's receive only multiply scattered contributions. The ratios of the signals in the outer FOV's to the signal in the central FOV is an indication of the lateral spreading of the scattered component of the laser pulse as it propagates through the cloud. MFOV lidar returns from stratus clouds measured between October 1991 and March 1992 can be divided into two distinct types, those with large multiple scattering ratios and those with small ratios. An example of each type of return is shown. Both measurements were made at an elevation angle of 45 degrees. Clouds with small multiple scattering signals probably have a high concentration of much larger particles on the order of hundreds of micrometers in size. This is a typical size range for suspended ice crystals or precipitations. Stratus clouds often have a high concentration of ice crystals even when there is no precipitation. Large ice crystals would give smaller signals in the outer FOV's because much of the scattered intensity is contained in a narrow diffraction peak with an angular width on the order of milliradians. The result is that for a given extinction, many more orders of scattering are required for the laser pulse to spread out. So far we have not been able to do simulations of MFOV lidar returns from ice crystal clouds because of uncertainties about the phase function of the crystals, particularly the magnitude of the backscatter peak at 180 degrees. On two occasions, MFOV lidar returns measured just prior to snowfall, showed a striking vertical profile. An example is shown. From 900 to 1300m, the multiply scattered signals are negligible compared to the return in the central FOV. Abruptly, at 1300m strong signals begin in the outer FOVs. These results could be explained by the presence of a cloud layer composed of water droplets at a range of 1300m (or 900m above ground) with snow precipitations forming at the bottom of the cloud. By measuring the polarization ratio of lidar signals, Pal and Carswell have observed snow precipitating from a cloud layer composed of water droplets. Based on these results, we believe that analysis of multiple scattering lidar returns can yield important information on the phase of cloud particles.

Hutt, D. L.↗