Search NASA⌕ Search

SEARCH · Search NASA

Results for “Curve 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 505 records · Page 28

Distribution of Cost Growth in Robotic Space Science Missions

Cost growth characterization is a critical factor for effective cost risk analysis and project planning. This study analyzed low level budget changes in Jet Propulsion Laboratory-managed space science missions, which occurred during the development of the project. The data was then curve fit, according to cost distribution categories, to provide a reference set of distribution parameters with sufficient granularity to effectively model cost growth in robotic space science missions.

cost↗

Loss Factor Estimation Using the Impulse Response Decay Method on a Stiffened Structure

High-frequency vibroacoustic modeling is typically performed using energy-based techniques such as Statistical Energy Analysis (SEA). Energy models require an estimate of the internal damping loss factor. Unfortunately, the loss factor is difficult to estimate analytically, and experimental methods such as the power injection method can require extensive measurements over the structure of interest. This paper discusses the implications of estimating damping loss factors using the impulse response decay method (IRDM) from a limited set of response measurements. An automated procedure for implementing IRDM is described and then evaluated using data from a finite element model of a stiffened, curved panel. Estimated loss factors are compared with loss factors computed using a power injection method and a manual curve fit. The paper discusses the sensitivity of the IRDM loss factor estimates to damping of connected subsystems and the number and location of points in the measurement ensemble.

Cabell, Randolph↗

Translational Vestibulo-Ocular Reflexes During Off-Vertical Axis Rotation

The translational vestibulo-ocular reflex (tVOR) is an otolith-mediated response that stabilizes near vision during linear acceleration at higher frequencies where visually mediated reflexes are not adequate. The modulation of horizontal and vergence eye movements during Off-Vertical Axis Rotation (OVAR) are presumed to reflect the tVOR in response to the continuously varying linear acceleration in the interaural and nasooccipital axes, respectively. The purpose of this study was to examine the effect of frequency and fixation distance on the modulation of slow phase eye velocity (SPV) as further evidence that the tVOR is elicited during OVAR. Eighteen subjects were rotated about their longitudinal axis tilted by 30 deg off-vertical. Rotational velocities varied between 18 and 288 deg/sec corresponding to a frequency range of 0.05 to 0.8 Hz. Fixation distance was altered by asking subjects to imagine stationary targets that were briefly presented at 0.5, 1 and 2 m during some rotation cycles. The target flash was 40 msec in the nose-up position at eye level. Oculomotor responses were recorded in the dark using infrared binocular videography. Sinusoidal curve fits were used to derive amplitude, phase and bias velocity of the eye movements across multiple rotation cycles. Consistent with previous studies, the modulation of both horizontal and vergence SPV increased with stimulus frequency. The effect of fixation distance was negligible at lower frequencies. The modulation of horizontal and vergence SPV was; however, proportional to fixation distance during OVAR at 0.8 Hz. This increasing sensitivity and dependence on fixation distance of horizontal and vergence SPV during OVAR is consistent with tVOR characteristics measured during other types of linear motion. We conclude that the modulation of horizontal and vergence SPV will be diagnostically more useful at higher stimulus frequencies where the tVOR is more robust.

Wood, Scott J.↗

Smooth Phase Interpolated Keying

Smooth phase interpolated keying (SPIK) is an improved method of computing smooth phase-modulation waveforms for radio communication systems that convey digital information. SPIK is applicable to a variety of phase-shift-keying (PSK) modulation schemes, including quaternary PSK (QPSK), octonary PSK (8PSK), and 16PSK. In comparison with a related prior method, SPIK offers advantages of better performance and less complexity of implementation. In a PSK scheme, the underlying information waveform that one seeks to convey consists of discrete rectangular steps, but the spectral width of such a waveform is excessive for practical radio communication. Therefore, the problem is to smooth the step phase waveform in such a manner as to maintain power and bandwidth efficiency without incurring an unacceptably large error rate and without introducing undesired variations in the amplitude of the affected radio signal. Although the ideal constellation of PSK phasor points does not cause amplitude variations, filtering of the modulation waveform (in which, typically, a rectangular pulse is converted to a square-root raised cosine pulse) causes amplitude fluctuations. If a power-efficient nonlinear amplifier is used in the radio communication system, the fluctuating-amplitude signal can undergo significant spectral regrowth, thus compromising the bandwidth efficiency of the system. In the related prior method, one seeks to solve the problem in a procedure that comprises two major steps: phase-value generation and phase interpolation. SPIK follows the two-step approach of the related prior method, but the details of the steps are different. In the phase-value-generation step, the phase values of symbols in the PSK constellation are determined by a phase function that is said to be maximally smooth and that is chosen to minimize the spectral spread of the modulated signal. In this step, the constellation is divided into two groups by assigning, to information symbols, phase values that result in equal numbers of clockwise and counter-clockwise phase rotations for equally likely symbols. The purpose served by assigning phase values in this way is to prevent unnecessary generation of spectral lines and prevent net shifts of the carrier signal. In the phase-interpolation step, the smooth phase values are interpolated over a number, n, of consecutive symbols (including the present symbol) by means of an unconventional spline curve fit.

Borah, Deva K.↗

Platform for Post-Processing Waveform-Based NDE

Signal- and image-processing methods are commonly needed to extract information from the waves, improve resolution of, and highlight defects in an image. Since some similarity exists for all waveform-based nondestructive evaluation (NDE) methods, it would seem that a common software platform containing multiple signal- and image-processing techniques to process the waveforms and images makes sense where multiple techniques, scientists, engineers, and organizations are involved. NDE Wave & Image Processor Version 2.0 software provides a single, integrated signal- and image-processing and analysis environment for total NDE data processing and analysis. It brings some of the most useful algorithms developed for NDE over the past 20 years into a commercial-grade product. The software can import signal/spectroscopic data, image data, and image series data. This software offers the user hundreds of basic and advanced signal- and image-processing capabilities including esoteric 1D and 2D wavelet-based de-noising, de-trending, and filtering. Batch processing is included for signal- and image-processing capability so that an optimized sequence of processing operations can be applied to entire folders of signals, spectra, and images. Additionally, an extensive interactive model-based curve-fitting facility has been included to allow fitting of spectroscopy data such as from Raman spectroscopy. An extensive joint-time frequency module is included for analysis of non-stationary or transient data such as that from acoustic emission, vibration, or earthquake data.

Roth, Don J.↗

Template Matching Approach to Signal Prediction

A new approach to signal prediction and prognostic assessment of spacecraft health resolves an inherent difficulty in fusing sensor data with simulated data. This technique builds upon previous work that demonstrated the importance of physics-based transient models to accurate prediction of signal dynamics and system performance. While models can greatly improve predictive accuracy, they are difficult to apply in general because of variations in model type, accuracy, or intended purpose. However, virtually any flight project will have at least some modeling capability at its disposal, whether a full-blown simulation, partial physics models, dynamic look-up tables, a brassboard analogue system, or simple hand-driven calculation by a team of experts. Many models can be used to develop a predict, or an estimate of the next day s or next cycle s behavior, which is typically used for planning purposes. The fidelity of a predict varies from one project to another, depending on the complexity of the simulation (i.e. linearized or full differential equations) and the level of detail in anticipated system operation, but typically any predict cannot be adapted to changing conditions or adjusted spacecraft command execution. Applying a predict blindly, without adapting the predict to current conditions, produces mixed results at best, primarily due to mismatches between assumed execution of spacecraft activities and actual times of execution. This results in the predict becoming useless during periods of complicated behavior, exactly when the predict would be most valuable. Each spacecraft operation tends to show up as a transient in the data, and if the transients are misaligned, using the predict can actually harm forecasting performance. To address this problem, the approach here expresses the predict in terms of a baseline function superposed with one or more transient functions. These transients serve as signal templates, which can be relocated in time and space against the signal background. One then has the ability to reconstruct a signal regardless of the precise timing of the transients. During operation, one applies the actual start times of spacecraft activities as they occur, and produces a reconstructed, accurate predict in real-time. This general approach is valid under two important conditions. First, the transients themselves must be time-invariant. Second, the transients must be reasonably consistent with respect to different operating points. Both of these assumptions are generally valid, but in the case of a complicated system with numerous types of overlapping transients, this approach may not be effective. Fortunately, there tend to be few transients in spacecraft telemetry of sensor quantities or low-level health and status information because these signals rarely reflect multiple different types of operation. Furthermore, if the predict is at least reasonably close to actual operation, the shift either in time, or in the operating point when the transient occurs - is likely to be small. The proposed approach considers three ways to recognize a transient. The first and most reliable is to use a different signal that identifies operating mode - often the transient will be correlated to a change in operating mode, and this change can usually be detected positively from discrete signals in spacecraft telemetry. If there is no useful mode signal, the second option is to identify the transient template by hand, and detect the onset of the transient using a curve-fitting approach. Finally, if one elects not to choose by hand, there is an option for automatic selection. This technique has been applied to sensor data from several JPL missions and industrial applications, demonstrating an improvement in accurate prediction of future behavior and early detection of system problems. This technique is applicable to practically any time-varying, quantitative sensor measurement.

Mackey, Ryan↗

Method and apparatus for shape and end position determination using an optical fiber

A method of determining the shape of an unbound optical fiber includes collecting strain data along a length of the fiber, calculating curvature and bending direction data of the fiber using the strain data, curve-fitting the curvature and bending direction data to derive curvature and bending direction functions, calculating a torsion function using the bending direction function, and determining the 3D shape from the curvature, bending direction, and torsion functions. An apparatus for determining the 3D shape of the fiber includes a fiber optic cable unbound with respect to a protective sleeve, strain sensors positioned along the cable, and a controller in communication with the sensors. The controller has an algorithm for determining a 3D shape and end position of the fiber by calculating a set of curvature and bending direction data, deriving curvature, bending, and torsion functions, and solving Frenet-Serret equations using these functions.

Moore, Jason P.↗

Dynamic Analysis of Sounding Rocket Pneumatic System Revision

The recent fusion of decades of advancements in mathematical models, numerical algorithms and curve fitting techniques marked the beginning of a new era in the science of simulation. It is becoming indispensable to the study of rockets and aerospace analysis. In pneumatic system, which is the main focus of this paper, particular emphasis will be placed on the efforts of compressible flow in Attitude Control System of sounding rocket.

Armen, Jerald↗

Flow Property Measurement Using Laser-Induced Fluorescence in the NASA Ames Interaction Heating Facility

The spectroscopic diagnostic technique of two photon absorption laser-induced fluorescence (TALIF) of atomic species has been applied to single-point measurements of velocity and static temperature in the NASA Ames Interaction Heating Facility (IHF) arc jet. Excitation spectra of atomic oxygen and nitrogen were recorded while scanning a tunable dye laser over the absorption feature. Thirty excitation spectra were acquired during 8 arc jet runs at two facility operating conditions; the number of scans per run varied between 2 and 6. Curve fits to the spectra were analyzed to recover their Doppler shifts and widths, from which the flow velocities and static temperatures, respectively, were determined. An increase in the number of independent flow property pairs from each as-measured scan was obtained by extracting multiple lower-resolution scans. The larger population sample size enabled the mean property values and their uncertainties for each run to be characterized with greater confidence. The average plus or minus 2 sigma uncertainties in the mean velocities and temperatures for all 8 runs were plus or minus 1.4% and plus or minus 11%, respectively.

Grinstead, Jay Henderson↗

Comparison of Ablation Predictions for Carbonaceous Materials Using CEA and JANAF-Based Species Thermodynamics

In most previous work at NASA Ames Research Center, ablation predictions for carbonaceous materials were obtained using a species thermodynamics database developed by Aerotherm Corporation. This database is derived mostly from the JANAF thermochemical tables. However, the CEA thermodynamics database, also used by NASA, is considered more up to date. In this work, the FIAT code was modified to use CEA-based curve fits for species thermodynamics, then analyses using both the JANAF and CEA thermodynamics were performed for carbon and carbon phenolic materials over a range of test conditions. The ablation predictions are comparable at lower heat fluxes where the dominant mechanism is carbon oxidation. However, the predictions begin to diverge in the sublimation regime, with the CEA model predicting lower recession. The disagreement is more significant for carbon phenolic than for carbon, and this difference is attributed to hydrocarbon species that may contribute to the ablation rate.

Milos, Frank S.↗

A Model for Hydrogen Thermal Conductivity and Viscosity Including the Critical Point

In order to conduct a thermal analysis of heat transfer to liquid hydrogen near the critical point, an accurate understanding of the thermal transport properties is required. A review of the available literature on hydrogen transport properties identified a lack of useful equations to predict the thermal conductivity and viscosity of liquid hydrogen. The tables published by the National Bureau of Standards were used to perform a series of curve fits to generate the needed correlation equations. These equations give the thermal conductivity and viscosity of hydrogen below 100 K. They agree with the published NBS tables, with less than a 1.5 percent error for temperatures below 100 K and pressures from the triple point to 1000 KPa. These equations also capture the divergence in the thermal conductivity at the critical point

Wagner, Howard A.↗

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↗

Receive Mode Analysis and Design of Microstrip Reflectarrays

Traditionally microstrip or printed reflectarrays are designed using the transmit mode technique. In this method, the size of each printed element is chosen so as to provide the required value of the reflection phase such that a collimated beam results along a given direction. The reflection phase of each printed element is approximated using an infinite array model. The infinite array model is an excellent engineering approximation for a large microstrip array since the size or orientation of elements exhibits a slow spatial variation. In this model, the reflection phase from a given printed element is approximated by that of an infinite array of elements of the same size and orientation when illuminated by a local plane wave. Thus the reflection phase is a function of the size (or orientation) of the element, the elevation and azimuth angles of incidence of a local plane wave, and polarization. Typically, one computes the reflection phase of the infinite array as a function of several parameters such as size/orientation, elevation and azimuth angles of incidence, and in some cases for vertical and horizontal polarization. The design requires the selection of the size/orientation of the printed element to realize the required phase by interpolating or curve fitting all the computed data. This is a substantially complicated problem, especially in applications requiring a computationally intensive commercial code to determine the reflection phase. In dual polarization applications requiring rectangular patches, one needs to determine the reflection phase as a function of five parameters (dimensions of the rectangular patch, elevation and azimuth angles of incidence, and polarization). This is an extremely complex problem. The new method employs the reciprocity principle and reaction concept, two well-known concepts in electromagnetics to derive the receive mode analysis and design techniques. In the "receive mode design" technique, the reflection phase is computed for a plane wave incident on the reflectarray from the direction of the beam peak. In antenna applications with a single collimated beam, this method is extremely simple since all printed elements see the same angles of incidence. Thus the number of parameters is reduced by two when compared to the transmit mode design. The reflection phase computation as a function of five parameters in the rectangular patch array discussed previously is reduced to a computational problem with three parameters in the receive mode. Furthermore, if the beam peak is in the broadside direction, the receive mode design is polarization independent and the reflection phase computation is a function of two parameters only. For a square patch array, it is a function of the size, one parameter only, thus making it extremely simple.

Rengarajan, Sembiam↗

Non-Contact Thrust Stand Calibration Method for Repetitively-Pulsed Electric Thrusters

A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.

Wong, Andrea R.↗

Routines for Computing Pressure Drops in Venturis

A set of computer-program routines has been developed for calculating pressure drops and recoveries of flows through standard venturis, nozzle venturis, and orifices. Relative to prior methods used for such calculations, the method implemented by these routines offers greater accuracy because it involves fewer simplifying assumptions and is more generally applicable to wide ranges of flow conditions. These routines are based on conservation of momentum and energy equations for real nonideal fluids, the properties of which are calculated by curve-fitting subroutines based on empirical properties data. These routines are capable of representing cavitating, choked, non-cavitating, and unchoked flow conditions for liquids, gases, and supercritical fluids. For a computation of flow through a given venturi, nozzle venturi, or orifice, the routines determine which flow condition occurs: First, they calculate a throat pressure under the assumption that the flow is unchoked or non-cavitating, then they calculate the throat pressure under the assumption that the flow is choked or cavitating. The assumption that yields the higher throat pressure is selected as the correct one.

de Quay, Laurence↗

Insulator Surface Charge as a Function of Pressure: Theory and Simulation

A two-phase equilibrium model was developed to explain the discontinuous surface charge decay versus atmospheric pressure of insulators that had been charged triboelectrically. The two-phase model is an electrostatic form of the Langmuir Isotherm for ions adsorbed on a surface in equilibrium with ions in the gas phase. In this paper, the model was extended to account for vibrational states of the adsorbed surface ions via the vibrational partition function. An analysis is performed that rules out Paschen discharge as the cause of the discharge observed. Also, a numerical simulation is performed using NWChem to calculate the adsorption energies of ions on insulator surfaces for comparison to curve fit adsorption energies developed from the model and experimental data.

Hogue, Michael D.↗

Pulsar Polar Cap and Slot Gap Models: Confronting Fermi Data

Rotation-powered pulsars are excellent laboratories for studying particle acceleration as well as fundamental physics of strong gravity, strong magnetic fields and relativity. I will review acceleration and gamma-ray emission from the pulsar polar cap and slot gap. Predictions of these models can be tested with the data set on pulsars collected by the Large Area Telescope on the Fermi Gamma-Ray Telescope over the last four years, using both detailed light curve fitting and population synthesis.

Harding, Alice K.↗

A Standard Law for the Equatorward Drift of the Sunspot Zones

The latitudinal location of the sunspot zones in each hemisphere is determined by calculating the centroid position of sunspot areas for each solar rotation from May 1874 to June 2012. When these centroid positions are plotted and analyzed as functions of time from each sunspot cycle maximum there appears to be systematic differences in the positions and equatorward drift rates as a function of sunspot cycle amplitude. If, instead, these centroid positions are plotted and analyzed as functions of time from each sunspot cycle minimum then most of the differences in the positions and equatorward drift rates disappear. The differences that remain disappear entirely if curve fitting is used to determine the starting times (which vary by as much as 8 months from the times of minima). The sunspot zone latitudes and equatorward drift measured relative to this starting time follow a standard path for all cycles with no dependence upon cycle strength or hemispheric dominance. Although Cycle 23 was peculiar in its length and the strength of the polar fields it produced, it too shows no significant variation from this standard. This standard law, and the lack of variation with sunspot cycle characteristics, is consistent with Dynamo Wave mechanisms but not consistent with current Flux Transport Dynamo models for the equatorward drift of the sunspot zones.

Hathaway, David H.↗