Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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 1,171 records · Page 65

On-Line Fringe Tracking and Prediction at IOTA

The Infrared/Optical Telescope Array (IOTA) is a multi-aperture Michelson interferometer located on Mt. Hopkins near Tucson, Arizona. To enable viewing of fainter targets, an on-line fringe tracking system is presently under development at NASA Ames Research Center. The system has been developed off-line using actual data from IOTA, and is presently undergoing on-line implementation at IOTA. The system has two parts: (1) a fringe tracking system that identifies the center of a fringe packet by fitting a parametric model to the data; and (2) a fringe packet motion prediction system that uses characteristics of past fringe packets to predict fringe packet motion. Combined, this information will be used to optimize on-line the scanning trajectory, resulting in improved visibility of faint targets. Fringe packet identification is highly accurate and robust (99% of the 4000 fringe packets were identified correctly, the remaining 1% were either out of the scan range or too noisy to be seen) and is performed in 30-90 milliseconds on a Pentium II-based computer. Fringe packet prediction, currently performed using an adaptive linear predictor, delivers a 10% improvement over the baseline of predicting no motion.

Wilson, Edward↗

Multiwavelength Photometric Imaging of X-Ray and EUV Emission from Comet P/Tempel-Tuttle 1998

The unexpected discovery of x-ray emission from Comet Hyakutake in March 1996 (Lisse et ai. 1996) has produced a number of questions about the physical mechanism producing the radiation. The original detection and subsequent observations (Dennerl et ai. 1997, Mumma et al. 1997, Krasnopolsky et al. 1998, Owens et al. 1998. Lisse et al. 1999) have shown that the very soft (best fit thermal bremsstrahlung model kT approx. 0.2 keV) emission is due to an interaction between the solar wind and the comet's atmosphere. Using the results fiom the 15 comets detected to date in x-rays, we report on the latest results on cometary x-ray emission, including new results from Chandra and XMM. As-observed morphologies, spectra, and light curves will be discussed. Our emphasis will be on understanding the physical mechanism producing the emission, and using this to determine the nature of the cometary coma, the structure of the solar wind in the heliosphere, and the source of the local soft x-ray background. This work has been graciously supported by grants from the NASA Planetary Astronomy and Astrophysical Data Programs.

Lisse, Carey M.↗

Intrinsic Line Shape Measurements of the XRS Instrument on Astro-E2

The XRS instrument on the Astro-E2 observatory contains a substantially improved microcalorimeter array over the Astro-E mission. In addition to roughly a factor of 2 improvement in the detector resolution at 6 keV, the detector response is shown to be almost perfectly gaussian. We have made measurements of the detector response of the flight instrument, using a double crystal monochrometer at 4 and 8 keV, a 55-Fe internal conversion source, and x-ray induced fluorescence from a number of targets including Ti, Cu, and GaAs. The detector response has been measured to be entirely gaussian to at least 2 orders of magnitude down from the peak of the line or line complex. This is in sharp contrast to the results from the XRS on Astro-E where many channels exhibited excess counts on the high energy side of the spectral lines. Here we present details of the line shape measurement as well as the detector response as measured during the XRS ground calibration including details of the line fits and line models.

Porter, F. Scott↗

Three-dimensional organization of otolith-ocular reflexes in rhesus monkeys. III. Responses To translation

The three-dimensional (3-D) properties of the translational vestibulo-ocular reflexes (translational VORs) during lateral and fore-aft oscillations in complete darkness were studied in rhesus monkeys at frequencies between 0.16 and 25 Hz. In addition, constant velocity off-vertical axis rotations extended the frequency range to 0.02 Hz. During lateral motion, horizontal responses were in phase with linear velocity in the frequency range of 2-10 Hz. At both lower and higher frequencies, phase lags were introduced. Torsional response phase changed more than 180 degrees in the tested frequency range such that torsional eye movements, which could be regarded as compensatory to "an apparent roll tilt" at the lowest frequencies, became anticompensatory at all frequencies above approximately 1 Hz. These results suggest two functionally different frequency bandwidths for the translational VORs. In the low-frequency spectrum (<<0.5 Hz), horizontal responses compensatory to translation are small and high-pass-filtered whereas torsional response sensitivity is relatively frequency independent. At higher frequencies however, both horizontal and torsional response sensitivity and phase exhibit a similar frequency dependence, suggesting a common role during head translation. During up-down motion, vertical responses were in phase with translational velocity at 3-5 Hz but phase leads progressively increased for lower frequencies (>90 degrees at frequencies <0.2 Hz). No consistent dependence on static head orientation was observed for the vertical response components during up-down motion and the horizontal and torsional response components during lateral translation. The frequency response characteristics of the translational VORs were fitted by "periphery/brain stem" functions that related the linear acceleration input, transduced by primary otolith afferents, to the velocity signals providing the input to the velocity-to-position neural integrator and the oculomotor plant. The lowest-order, best-fit periphery/brain stem model that approximated the frequency dependence of the data consisted of a second order transfer function with two alternating poles (at 0.4 and 7.2 Hz) and zeros (at 0.035 and 3.4 Hz). In addition to clearly differentiator dynamics at low frequencies (less than approximately 0.5 Hz), there was no frequency bandwidth where the periphery/brain stem function could be approximated by an integrator, as previously suggested. In this scheme, the oculomotor plant dynamics are assumed to perform the necessary high-frequency integration as required by the reflex. The detailed frequency dependence of the data could only be precisely described by higher order functions with nonminimum phase characteristics that preclude simple filtering of afferent inputs and might be suggestive of distributed spatiotemporal processing of otolith signals in the translational VORs.

Non-NASA Center↗

Dependence of Ru2O3 Activity on Composition of Silicate Melts: Using Statistical Correlations to Infer Thermodynamic Behavior in the Melt

Understanding variations in activity with composition is an essential step in improving prediction of partition coefficients during magma evolution. Variations in activity with composition are complex and do not generally exhibit ideal behavior relative to a traditional melt-component set. Although deviations from component ideality can be modeled numerically by simply fitting to compositional variables (such as in a regular or subregular solution model), such models have not been particularly successful for describing variations in trace component activities. A better approach might be to try to identify components that do a better job of describing the behavior of the species in the melt. Electrochemical Measurement of Ru2O3 activities: Electrodes were inserted into silicate melt beads of various compositions (Table 1) suspended on Ptwire loops in a 1-atm gas mixing furnace. An electrical potential was imposed between the electrodes, the imposed potential increasing along a step ramp with a pulse imposed on each step (Fig. 1). Current flows between electrodes when electroactive species in the melt are oxidized or reduced at the electrodes. The resulting current was measured at the top and bottom of the voltage pulse, and the difference (the differential current) was plotted against potential. The peak of the resulting curve is related to the activity coefficient for the particular electroactive species (Ru2O3) in the melt [1, 2, 3]. A significant part of the nonideal contribution to activity is due not to intrinsic properties of the component in the melt, but to our ignorance about the state and mixing properties of the component in the melt.

Colson, R. O.↗

Water-Group Ion Distributions in the Mid-Cometosheath of Comet Halley

In the mid-cometosheath of comet Halley (1-2x10^5 km from the nucleus) the center- of-mass plasma frame is approximately the bulk flow velocity of the cometary ions, and the Alfven wave speed is an appreciable fraction of the flow speed. Here, the peaks of the water-group ion distributions observed by the Giotto Ion Mass Spectrometer are at velocities consistently below the expected pickup speed. It is shown that this effect is consistent with the scattering of the new pickup ions onto a bispherical shell distribution. The model does not fit the data inside similar 1.2x10^5 km however, possibly as a result of the growing importance of collisions or the presence of other processes such as scattering on obliquely-propagating magnetosonic waves.

Huddleston, D. E.↗

Comet Mineralogy as Inferred from Infrared Spectra of Comets

For most comets, infrared (IR) spectroscopy (remote sensing) is the method through which we diagnose the mineralogy and size distribution of dust in their comae. The shape and contrast of the IR spectral features depend on the particle size: optically active minerals (absorbing of visible and near-IR solar photons) and submicron solid grains or highly porous (> 90% vacuum) grains primarily contribute to the shapes of the observed resonances. Comet mineralogies typically are determined by fitting thermal emission models of ensembles of discrete mineral grains to observed IR spectral energy distributions. The absorptivities (Q-abs) and scattering efficiencies (Q-scat) of the discrete mineral grains are computed using Mie scattering, Maxwell-Garnet mixing, Discrete Dipole Approximation, and Multi-Layered Sphere codes. These techniques when applied to crystalline minerals, specifically olivine (Mg_x, Fe_1-x)2 Si04, x>0.9, require the use of ellipsoidal shaped particles with elongated axial ratios or hollow spheres to produce the shapes of the resonances observed both from comet comae and laboratory samples. The wavelength positions of the distinct resonances from submicron-radii crystalline silicates, as well as their thermal equilibrium temperatures, constrain the crystalline olivine to have a relatively high Mg-content (x>0.9, or Fo>90). Only resonances computed for submicron Mg-rich crystalline olivine and crystalline orthopyroxene match the observed IR spectral features. However, this has led to the interpretation that micron-radii and larger crystals are absent from comet comae. Furthermore, the mass fraction of silicate crystals is dependent upon whether just the submicron portion of the size distribution is being compared or the submicron crystals compare to the aggregates of porous amorphous silicates that are computationally tractable as porous spheres. We will discuss the Deep Impact results as examples of these challenges to interpreting mid-IR spectra of short-period comets.

Wooden, Diane H.↗

A New Look at Some Solar Wind Turbulence Puzzles

Some aspects of solar wind turbulence have defied explanation. While it seems likely that the evolution of Alfvenicity and power spectra are largely explained by the shearing of an initial population of solar-generated Alfvenic fluctuations, the evolution of the anisotropies of the turbulence does not fit into the model so far. A two-component model, consisting of slab waves and quasi-two-dimensional fluctuations, offers some ideas, but does not account for the turning of both wave-vector-space power anisotropies and minimum variance directions in the fluctuating vectors as the Parker spiral turns. We will show observations that indicate that the minimum variance evolution is likely not due to traditional turbulence mechanisms, and offer arguments that the idea of two-component turbulence is at best a local approximation that is of little help in explaining the evolution of the fluctuations. Finally, time-permitting, we will discuss some observations that suggest that the low Alfvenicity of many regions of the solar wind in the inner heliosphere is not due to turbulent evolution, but rather to the existence of convected structures, including mini-clouds and other twisted flux tubes, that were formed with low Alfvenicity. There is still a role for turbulence in the above picture, but it is highly modified from the traditional views.

Roberts, Aaron↗

A Revised Model of Jupiter's Inner Electron Belts: Updating the Divine Radiation Model

In 1983, Divine presented a comprehensive model of the Jovian charged particle environment that has long served as a reference for missions to Jupiter. However, in situ observations by Galileo and synchrotron observations from Earth indicate the need to update the model in the inner radiation zone. Specifically, a review of the model for 1 MeV < E < 100 MeV trapped electrons suggests that, based on the new synchrotron observations, the pitch angle distributions within L < 4 need to be updated by introducing two additional components: one near the Jovian magnetic equator and one at high magnetic latitudes. We report modifications to the model that reproduce these observations. The new model improves the fit to synchrotron emission observations and remains consistent with the original fit to the in situ Pioneer and Voyager data. Further modifications incorporating observations from the Galileo and Cassini spacecraft will be reported in the future.

missions↗

Workflow Agents vs. Expert Systems: Problem Solving Methods in Work Systems Design

During the 1980s, a community of artificial intelligence researchers became interested in formalizing problem solving methods as part of an effort called "second generation expert systems" (2nd GES). How do the motivations and results of this research relate to building tools for the workplace today? We provide an historical review of how the theory of expertise has developed, a progress report on a tool for designing and implementing model-based automation (Brahms), and a concrete example how we apply 2nd GES concepts today in an agent-based system for space flight operations (OCAMS). Brahms incorporates an ontology for modeling work practices, what people are doing in the course of a day, characterized as "activities." OCAMS was developed using a simulation-to-implementation methodology, in which a prototype tool was embedded in a simulation of future work practices. OCAMS uses model-based methods to interactively plan its actions and keep track of the work to be done. The problem solving methods of practice are interactive, employing reasoning for and through action in the real world. Analogously, it is as if a medical expert system were charged not just with interpreting culture results, but actually interacting with a patient. Our perspective shifts from building a "problem solving" (expert) system to building an actor in the world. The reusable components in work system designs include entire "problem solvers" (e.g., a planning subsystem), interoperability frameworks, and workflow agents that use and revise models dynamically in a network of people and tools. Consequently, the research focus shifts so "problem solving methods" include ways of knowing that models do not fit the world, and ways of interacting with other agents and people to gain or verify information and (ultimately) adapt rules and procedures to resolve problematic situations.

Clancey, William J.↗

Ground Operations of the ISS GNC Babb-Mueller Atmospheric Density Model

The ISS GNC system was updated recently with a new software release that provides onboard state determination capability. Prior to this release, only the Russian segment maintained and propagated the onboard state, which was periodically updated through Russian ground tracking. The new software gives the US segment the capability for maintaining the onboard state, and includes new GPS and state vector propagation capabilities. Part of this software package is an atmospheric density model based on the Babb-Mueller algorithm. Babb-Mueller efficiently mimics a full analytical density model, such as the Jacchia model. While lacchia is very robust and is used in the Mission Control Center, it is too computationally intensive for use onboard. Thus, Babb-Mueller was chosen as an alternative. The onboard model depends on a set of calibration coefficients that produce a curve fit to the lacchia model. The ISS GNC system only maintains one set of coefficients onboard, so a new set must be uplinked by controllers when the atmospheric conditions change. The onboard density model provides a real-time density value, which is used to calculate the drag experienced by the ISS. This drag value is then incorporated into the onboard propagation of the state vector. The propagation of the state vector, and therefore operation of the BabbMueller algorithm, will be most critical when GPS updates and secondary state vector sources fail. When GPS is active, the onboard state vector will be updated every ten seconds, so the propagation error is irrelevant. When GPS is inactive, the state vector must be updated at least every 24 hours, based on current protocol. Therefore, the Babb-Mueller coefficients must be accurate enough to fulfill the state vector accuracy requirements for at least one day. A ground operations concept was needed in order to manage both the on board Babb-Mueller density model and the onboard state quality. The Babb-Mueller coefficients can be determined operationally in two ways. The first method is to calibrate the coefficients in real-time, where a set of custom coefficients is generated for the real-time atmospheric conditions. The second approach is to generate pre-canned sets of coefficients that encompass the expected atmospheric conditions over the lifetime of the vehicle. These predetermined sets are known as occurrences. Even though a particular occurrence will not match the true atmospheric conditions, the error will be constrained by limiting the breadth of each occurrence. Both methods were investigated and the advantages and disadvantages of each were considered. The choice between these implementations was a trade-off between the additional accuracy of the real-time calibration and the simpler development for the approach using occurrences. The operations concept for the frequency of updates was also explored, and depends on the deviation in solar flux that still achieves the necessary accuracy of the coefficients. This was determined based on historical solar flux trends. This analysis resulted in an accurate and reliable implementation of the Babb-Mueller coefficients and how flight controllers use them during realtime operations.

Brogan, Jonathan↗

Adaptive Sampling of Time Series During Remote Exploration

This work deals with the challenge of online adaptive data collection in a time series. A remote sensor or explorer agent adapts its rate of data collection in order to track anomalous events while obeying constraints on time and power. This problem is challenging because the agent has limited visibility (all its datapoints lie in the past) and limited control (it can only decide when to collect its next datapoint). This problem is treated from an information-theoretic perspective, fitting a probabilistic model to collected data and optimizing the future sampling strategy to maximize information gain. The performance characteristics of stationary and nonstationary Gaussian process models are compared. Self-throttling sensors could benefit environmental sensor networks and monitoring as well as robotic exploration. Explorer agents can improve performance by adjusting their data collection rate, preserving scarce power or bandwidth resources during uninteresting times while fully covering anomalous events of interest. For example, a remote earthquake sensor could conserve power by limiting its measurements during normal conditions and increasing its cadence during rare earthquake events. A similar capability could improve sensor platforms traversing a fixed trajectory, such as an exploration rover transect or a deep space flyby. These agents can adapt observation times to improve sample coverage during moments of rapid change. An adaptive sampling approach couples sensor autonomy, instrument interpretation, and sampling. The challenge is addressed as an active learning problem, which already has extensive theoretical treatment in the statistics and machine learning literature. A statistical Gaussian process (GP) model is employed to guide sample decisions that maximize information gain. Nonsta tion - ary (e.g., time-varying) covariance relationships permit the system to represent and track local anomalies, in contrast with current GP approaches. Most common GP models are stationary, e.g., the covariance relationships are time-invariant. In such cases, information gain is independent of previously collected data, and the optimal solution can always be computed in advance. Information-optimal sampling of a stationary GP time series thus reduces to even spacing, and such models are not appropriate for tracking localized anomalies. Additionally, GP model inference can be computationally expensive.

Thompson, David R.↗

Measurements of DSD Second Moment Based on Laser Extinction

Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through rain (or dust in the rocket case) yields an estimate of the 2nd moment of the particle cloud, and rainfall drop size distribution (DSD) in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rainfall make in direct measurements of the DSD. Most common of these instruments are the rainfall rate gauge measuring the 1 1/3 th moment, (when using a D(exp 2/3) dependency on terminal velocity). Instruments that scatter microwaves off of hydrometeors, such as the WSR-880, vertical wind profilers, and microwave disdrometers, measure the 6th moment of the DSD. By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain, yield a measurement of the DSD 2nd moment, using the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required, depending on the intensity of the rainfall rate. For moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. A photo-detector could replace the camera, for automated installations. In order to spatially correlate the 2nd moment measurements to a collocated disdrometer or tipping bucket, the laser's beam path can be reflected multiple times using mirrors to restrict the spatial extent of the measurement. In cases where a disdrometer is not available, complete DSD estimates can be produced by parametric fitting of DSD model to the 2nd moment data in conjunction with tipping bucket data. In cases where a disdrometer is collocated, the laser extinction technique may yield a significant improvement to insitu disdrometer validation and calibration strategies

Lane, John E.↗

Improved CPAS Photogrammetric Capabilities for Engineering Development Unit (EDU) Testing

This paper focuses on two key improvements to the photogrammetric analysis capabilities of the Capsule Parachute Assembly System (CPAS) for the Orion vehicle. The Engineering Development Unit (EDU) system deploys Drogue and Pilot parachutes via mortar, where an important metric is the muzzle velocity. This can be estimated using a high speed camera pointed along the mortar trajectory. The distance to the camera is computed from the apparent size of features of known dimension. This method was validated with a ground test and compares favorably with simulations. The second major photogrammetric product is measuring the geometry of the Main parachute cluster during steady-state descent using onboard cameras. This is challenging as the current test vehicles are suspended by a single-point attachment unlike earlier stable platforms suspended under a confluence fitting. The mathematical modeling of fly-out angles and projected areas has undergone significant revision. As the test program continues, several lessons were learned about optimizing the camera usage, installation, and settings to obtain the highest quality imagery possible.

Ray, Eric S.↗

Lunar Regolith Excavation Competition

The Lunar Regolith Excavation Competition is a new competition that needs graphics, logos, rules, as well as an arena. Although this is the first year of the competition, the competition is modeled after an existing competition, the Centennial Lunar Excavator Challenge. This competition however is aimed at college students. This makes the challenge identifying key aspects of the original competition and modeling them to fit into an easier task, and creating exciting advertisement that helps encourage participation. By using a youth focus group, young insight, as well as guiding advice from experts in the field, hopefully an arena can be designed and built, rules can be molded and created to fit, and alluring graphics can be printed to bring about a successful first year of the Lunar Regolith Excavation Competition.

Liles, Cassandra↗

Using CFD Techniques to Predict Slosh Force Frequency and Damping Rate

Resonant effects and energy dissipation due to sloshing fuel inside propellant tanks are problems that arise in the initial design of any spacecraft or launch vehicle. A faster and more reliable method for calculating these effects during the design stages is needed. Using Computational Fluid Dynamics (CFD) techniques, a model of these fuel tanks can be created and used to predict important parameters such as resonant slosh frequency and damping rate. This initial study addresses the case of free surface slosh. Future studies will focus on creating models for tanks fitted with propellant management devices (PMD) such as diaphragms and baffles.

Marsell, Brandon↗

Automating Hyperspectral Data for Rapid Response in Volcanic Emergencies

In a volcanic emergency, time is of the essence. It is vital to quantify eruption parameters (thermal emission, effusion rate, location of activity) and distribute this information as quickly as possible to decision-makers in order to enable effective evaluation of eruption-related risk and hazard. The goal of this work was to automate and streamline processing of spacecraft hyperspectral data, automate product generation, and automate distribution of products. Visible and Short-Wave Infrared Images of volcanic eruption in Iceland in May 2010." class="caption" align="right">The software rapidly processes hyperspectral data, correcting for incident sunlight where necessary, and atmospheric transmission; detects thermally anomalous pixels; fits data with model black-body thermal emission spectra to determine radiant flux; calculates atmospheric convection thermal removal; and then calculates total heat loss. From these results, an estimation of effusion rate is made. Maps are generated of thermal emission and location (see figure). Products are posted online, and relevant parties notified. Effusion rate data are added to historical record and plotted to identify spikes in activity for persistently active eruptions. The entire process from start to end is autonomous. Future spacecraft, especially those in deep space, can react to detection of transient processes without the need to communicate with Earth, thus increasing science return. Terrestrially, this removes the need for human intervention.

Davies, Ashley G.↗

Characterization of New Hard X-ray Cataclysmic Variables

Aims. We aim at characterizing a sample of nine new hard X-ray selected Cataclysmic Variable (CVs), to unambiguously identify them as magnetic systems of the Intermediate Polar (IP) type. Methods. We performed detailed timing and spectral analysis by using X-ray, and simultaneous UV and optical data collected by XMM−Newton, complemented with hard X-ray data provided by INTEGRAL and Swift. The pulse arrival time were used to estimate the orbital periods. The broad band X-ray spectra were fitted using composite models consisting of different absorbing columns and emission components. Results. Strong X-ray pulses at the White Dwarf (WD) spin period are detected and found to decrease with energy. Most sources are spin-dominated systems in the X-rays, though four are beat dominated at optical wavelengths. We estimated the orbital period in all system (except for IGR J16500-3307), providing the first estimate for IGRJ08390-4833, IGRJ18308-1232, and IGR J18173-2509. All X-ray spectra are multi-temperature. V2069 Cyg and RX J0636+3535 poses a soft X-ray optically thick component at kT approx. 80 eV. An intense K (sub alpha) Fe line at 6.4 keV is detected in all sources. An absorption edge at 0.76 keV from OVII is detected in IGR J08390-4833. The WD masses and lower limits to the accretion rates are also estimated. Conclusions. We found all sources to be IPs. IGR J08390-4833, V2069 Cyg, and IGR J16500-3307 are pure disc accretors, while IGR J18308-1232, IGR J1509-6649, IGR J17195-4100, and RX J0636+3535 display a disc-overflow accretion mode. All sources show a temperature gradient in the post-shock regions and a highly absorbed emission from material located in the pre-shock flow which is also responsible for the X-ray pulsations. Reflection at the WD surface is likely the origin of the fluorescent iron line. There is an increasing evidence for the presence of a warm absorber in IPs, a feature that needs future exploration. The addition of two systems to the subgroup of soft X-ray IPs confirms a relatively large (∼ 30%) incidence.

hard X-ray↗