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 523 records · Page 29

Preliminary Modeling of CH3D From 4000 to 4550 cm (exp-1)

A new study of 12CH3D line positions and intensities was performed for the upper portion of the Enneadecad polyad between 4000 and 4550/cm. For this, FTIR spectra were recorded with D-enriched methane samples (at 80 K with a Bruker 125 IFS at 0.005/cm resolution and at 291 K with the McMath-Pierce FTS at 0.011/cm resolution, respectively). Line positions and intensities were retrieved by least square curve-fitting procedures and analyzed using the effective Hamiltonian and the effective dipole moment expressed in terms of irreducible tensor operators adapted to symmetric top molecules. Initially, only the cold spectrum was used to identify quantum assignments and predict 12CH3D relative intensities in this region. To assign higher quanta up to J equal 14, additional line positions and intensities were obtained from two room temperature spectra. In the final stage, measured intensities from both the cold and the room temperature data were normalized to corresponding values at 296 K and then averaged. Combining the two temperature datasets confirmed the assumed quantum assignments and also demonstrated the relative accuracies to be better than 70.0002/cm for line positions and at least +/-6% for intensities so that approx.1160 features were selected. Including additional assignments from the room temperature spectra alone permitted 1362 line intensities of 11 bands (involving 23 vibrational symmetry components) to be reproduced with an RMS of 9%. Over 4085 selected positions for 12 bands were modeled to 0.008/cm. Nevertheless a number of known assignments could not be modeled to within our experimental precisions. More work is needed to obtain a complete characterization of this complex polyad.

Nikitin, A. V.↗

A Model-Based Anomaly Detection Approach for Analyzing Streaming Aircraft Engine Measurement Data

This paper presents a model-based anomaly detection architecture designed for analyzing streaming transient aircraft engine measurement data. The technique calculates and monitors residuals between sensed engine outputs and model predicted outputs for anomaly detection purposes. Pivotal to the performance of this technique is the ability to construct a model that accurately reflects the nominal operating performance of the engine. The dynamic model applied in the architecture is a piecewise linear design comprising steady-state trim points and dynamic state space matrices. A simple curve-fitting technique for updating the model trim point information based on steadystate information extracted from available nominal engine measurement data is presented. Results from the application of the model-based approach for processing actual engine test data are shown. These include both nominal fault-free test case data and seeded fault test case data. The results indicate that the updates applied to improve the model trim point information also improve anomaly detection performance. Recommendations for follow-on enhancements to the technique are also presented and discussed.

Propulsion System Performance↗

Some Interesting Applications of Probabilistic Techiques in Structural Dynamic Analysis of Rocket Engines

Numerical and Analytical methods developed to determine damage accumulation in specific engine components when speed variation included. Dither Life Ratio shown to be well over factor of 2 for specific example. Steady‐State assumption shown to be accurate for most turbopump cases, allowing rapid calculation of DLR. If hot‐fire speed data unknown, Monte Carlo method developed that uses speed statistics for similar engines. Application of techniques allow analyst to reduce both uncertainty and excess conservatism. High values of DLR could allow previously unacceptable part to pass HCF criteria without redesign. Given benefit and ease of implementation, recommend that any finite life turbomachine component analysis adopt these techniques. Probability Values calculated, compared, and evaluated for several industry‐proposed methods for combining random and harmonic loads. Two new excel macros written to calculate combined load for any specific probability level. Closed form Curve fits generated for widely used 3(sigma) and 2(sigma) probability levels. For design of lightweight aerospace components, obtaining accurate, reproducible, statistically meaningful answer critical.

Brown, Andrew M.↗

Influence of Triply-Charged Ions and Ionization Cross-Sections in a Hybrid-PIC Model of a Hall Thruster Discharge

The sensitivity of xenon ionization rates to collision cross-sections is studied within the framework of a hybrid-PIC model of a Hall thruster discharge. A revised curve fit based on the Drawin form is proposed and is shown to better reproduce the measured crosssections at high electron energies, with differences in the integrated rate coefficients being on the order of 10% for electron temperatures between 20 eV and 30 eV. The revised fit is implemented into HPHall and the updated model is used to simulate NASA's HiVHAc EDU2 Hall thruster at discharge voltages of 300, 400, and 500 V. For all three operating points, the revised cross-sections result in an increase in the predicted thrust and anode efficiency, reducing the error relative to experimental performance measurements. Electron temperature and ionization reaction rates are shown to follow the trends expected based on the integrated rate coefficients. The effects of triply-charged xenon are also assessed. The predicted thruster performance is found to have little or no dependence on the presence of triply-charged ions. The fraction of ion current carried by triply-charged ions is found to be on the order of 1% and increases slightly with increasing discharge voltage. The reaction rates for the 0→III, I→III, and II→III ionization reactions are found to be of similar order of magnitude and are about one order of magnitude smaller than the rate of 0→II ionization in the discharge channel.

Hall Thrusters↗

A Model-Based Anomaly Detection Approach for Analyzing Streaming Aircraft Engine Measurement Data

This paper presents a model-based anomaly detection architecture designed for analyzing streaming transient aircraft engine measurement data. The technique calculates and monitors residuals between sensed engine outputs and model predicted outputs for anomaly detection purposes. Pivotal to the performance of this technique is the ability to construct a model that accurately reflects the nominal operating performance of the engine. The dynamic model applied in the architecture is a piecewise linear design comprising steady-state trim points and dynamic state space matrices. A simple curve-fitting technique for updating the model trim point information based on steadystate information extracted from available nominal engine measurement data is presented. Results from the application of the model-based approach for processing actual engine test data are shown. These include both nominal fault-free test case data and seeded fault test case data. The results indicate that the updates applied to improve the model trim point information also improve anomaly detection performance. Recommendations for follow-on enhancements to the technique are also presented and discussed.

Propulsion System Performance↗

Three Gravitationally Lensed Supernovae Behind Clash Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and A383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approx. 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approx. 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak-lensing maps of the clusters (in magnitude units, 2.5 log10 μ): 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

Three↗

A Methodology to Assess the Capability of Engine Designs to Meet Closed-loop Performance and Operability Requirements

Designing a closed-loop controller for an engine requires balancing trade-offs between performance and operability of the system. One such trade-off is the relationship between the 95% response time and minimum high-pressure compressor (HPC) surge margin (SM) attained during acceleration from idle to takeoff power. Assuming a controller has been designed to meet some specification on response time and minimum HPC SM for a mid-life (nominal) engine, there is no guarantee that these limits will not be violated as the engine ages, particularly as it reaches the end of its life. A characterization for the uncertainty in this closed-loop system due to aging is proposed that defines elliptical boundaries to estimate worst-case performance levels for a given control design point. The results of this characterization can be used to identify limiting design points that bound the possible con- troller designs yielding transient results that do not exceed specified limits in response time or minimum HPC SM. This characterization involves performing Monte Carlo simulation of the closed-loop system with controller constructed for a set of trial design points and developing curve fits to describe the size and orientation of each ellipse; a binary search procedure is then employed that uses these fits to identify the limiting design point. The method is demonstrated through application to a generic turbofan engine model in closed- loop with a simplified controller; it is found that the limit for which each controller was designed was exceeded by less than 4.76%. Extension of the characterization to another trade-off, that between the maximum high-pressure turbine (HPT) entrance temperature and minimum HPC SM, showed even better results: the maximum HPT temperature was estimated within 0.76%. Because of the accuracy in this estimation, this suggests another limit that may be taken into consideration during design and analysis. It also demonstrates the extension of the characterization to other attributes that contribute to the performance or operability of the engine. Metrics are proposed that, together, provide information on the shape of the trade-off between response time and minimum HPC SM, and how much each varies throughout the life cycle, at the limiting design points. These metrics also facilitate comparison of the expected transient behavior for multiple engine models.

engine control↗

A Methodology to Assess the Capability of Engine Designs to Meet Closed-Loop Performance and Operability Requirements

Designing a closed-loop controller for an engine requires balancing trade-offs between performance and operability of the system. One such trade-off is the relationship between the 95 percent response time and minimum high-pressure compressor (HPC) surge margin (SM) attained during acceleration from idle to takeoff power. Assuming a controller has been designed to meet some specification on response time and minimum HPC SM for a mid-life (nominal) engine, there is no guarantee that these limits will not be violated as the engine ages, particularly as it reaches the end of its life. A characterization for the uncertainty in this closed-loop system due to aging is proposed that defines elliptical boundaries to estimate worst-case performance levels for a given control design point. The results of this characterization can be used to identify limiting design points that bound the possible controller designs yielding transient results that do not exceed specified limits in response time or minimum HPC SM. This characterization involves performing Monte Carlo simulation of the closed-loop system with controller constructed for a set of trial design points and developing curve fits to describe the size and orientation of each ellipse; a binary search procedure is then employed that uses these fits to identify the limiting design point. The method is demonstrated through application to a generic turbofan engine model in closed-loop with a simplified controller; it is found that the limit for which each controller was designed was exceeded by less than 4.76 percent. Extension of the characterization to another trade-off, that between the maximum high-pressure turbine (HPT) entrance temperature and minimum HPC SM, showed even better results: the maximum HPT temperature was estimated within 0.76 percent. Because of the accuracy in this estimation, this suggests another limit that may be taken into consideration during design and analysis. It also demonstrates the extension of the characterization to other attributes that contribute to the performance or operability of the engine. Metrics are proposed that, together, provide information on the shape of the trade-off between response time and minimum HPC SM, and how much each varies throughout the life cycle, at the limiting design points. These metrics also facilitate comparison of the expected transient behavior for multiple engine models.

systems analysis↗

High-Performance Reaction Wheel Optimization for Fine-Pointing Space Platforms: Minimizing Induced Vibration Effects on Jitter Performance plus Lessons Learned from Hubble Space Telescope for Current and Future Spacecraft Applications

The Hubble Space Telescope (HST) applies large-diameter optics (2.5-m primary mirror) for diffraction-limited resolution spanning an extended wavelength range (approx. 100-2500 nm). Its Pointing Control System (PCS) Reaction Wheel Assemblies (RWAs), in the Support Systems Module (SSM), acquired an unprecedented set of high-sensitivity Induced Vibration (IV) data for 5 flight-certified RWAs: dwelling at set rotation rates. Focused on 4 key ratios, force and moment harmonic values (in 3 local principal directions) are extracted in the RWA operating range (0-3000 RPM). The IV test data, obtained under ambient lab conditions, are investigated in detail, evaluated, compiled, and curve-fitted; variational trends, core causes, and unforeseen anomalies are addressed. In aggregate, these values constitute a statistically-valid basis to quantify ground test-to-test variations and facilitate extrapolations to on-orbit conditions. Accumulated knowledge of bearing-rotor vibrational sources, corresponding harmonic contributions, and salient elements of IV key variability factors are discussed. An evolved methodology is presented for absolute assessments and relative comparisons of macro-level IV signal magnitude due to micro-level construction-assembly geometric details/imperfections stemming from both electrical drive and primary bearing design parameters. Based upon studies of same-size/similar-design momentum wheels' IV changes, upper estimates due to transitions from ground tests to orbital conditions are derived. Recommended HST RWA choices are discussed relative to system optimization/tradeoffs of Line-Of-Sight (LOS) vector-pointing focal-plane error driven by higher IV transmissibilities through low-damped structural dynamics that stimulate optical elements. Unique analytical disturbance results for orbital HST accelerations are described applicable to microgravity efforts. Conclusions, lessons learned, historical context/insights, and perspectives on future applications are given; these previously unpublished data and findings represents a valuable resource for fine-pointing spacecraft or space-based platforms using RWAs, Control Moment Gyros (CMGs), Momentum Wheels, or other ball-bearing-based rotational units.

Hasha, Martin D.↗

Data Analysis & Statistical Methods for Command File Errors

This paper explains current work on modeling for managing the risk of command file errors. It is focused on analyzing actual data from a JPL spaceflight mission to build models for evaluating and predicting error rates as a function of several key variables. We constructed a rich dataset by considering the number of errors, the number of files radiated, including the number commands and blocks in each file, as well as subjective estimates of workload and operational novelty. We have assessed these data using different curve fitting and distribution fitting techniques, such as multiple regression analysis, and maximum likelihood estimation to see how much of the variability in the error rates can be explained with these. We have also used goodness of fit testing strategies and principal component analysis to further assess our data. Finally, we constructed a model of expected error rates based on the what these statistics bore out as critical drivers to the error rate. This model allows project management to evaluate the error rate against a theoretically expected rate as well as anticipate future error rates.

Correlation Analysis↗

Characterization of DBD Plasma Actuators Performance without External Flow : Thrust-Voltage Quadratic Relationship in Logarithmic Space for Sinusoidal Excitation - Part I

We present results of thrust measurements of Dielectric Barrier Discharge (DBD) plasma actuators. We have used a test setup, measurement, and data processing methodology that we developed in prior work. The tests were conducted with High Density Polyethylene (HDPE) actuators of three thicknesses. The applied voltage driving the actuators was a pure sinusoidal waveform. The test setup was suspended actuators with a partial liquid interface. The tests were conducted at low ambient humidity. The thrust was measured with an analytical balance and the results were corrected for anti-thrust to isolate the plasma generated thrust. Applying this approach resulted in smooth and repeatable data. It also enabled curve fitting that yielded quadratic relations between the plasma thrust and voltage in log-log space at constant frequencies. The results contrast power law relationships developed in literature that appear to be a rough approximation over a limited voltage range.

Dielectric Barrier Discharge↗

Atomic Scattering Factor of the ASTRO-H (Hitomi) SXT Reflector Around the Gold's L Edges

The atomic scattering factor in the energy range of 11.2 - 15.4 keV for the ASTRO-H Soft X-ray Telescope (SXT) is reported. The large effective area of the SXT makes use of photon spectra above 10 keV viable, unlike most other X-ray satellites with total-reflection mirror optics. Presence of gold's L-edges in the energy band is a major issue, as it complicates the function of the effective area. In order to model the area, the reflectivity measurements in the 11.2 - 15.4 keV band with the energy pitch of 0.4 - 0.7 eV were made in the synchrotron beam-line Spring-8 BL01B1. We obtained atomic scattering factors f1 and f2 by the curve fitting to the reflectivities of our witness sample. The edges associated with the L-I, II, and III transitions are identified, of which the depths are found to be roughly 60 shallower than those expected from the Henkes atomic scattering factor.

ASTRO-H Soft X-ray↗

Isolating Added Mass Load Components of CPAS Main Clusters

The current simulation for the Capsule Parachute Assembly System (CPAS) lacks fidelity in representing added mass for the 116 ft Do ringsail Main parachute. The availability of 3-D models of inflating Main canopies allowed for better estimation the enclosed air volume as a function of time. This was combined with trajectory state information to estimate the components making up measured axial loads. A proof-of-concept for an alternate simulation algorithm was developed based on enclosed volume as the primary independent variable rather than drag area growth. Databases of volume growth and parachute drag area vs. volume were developed for several flight tests. Other state information was read directly from test data, rather than numerically propagated. The resulting simulated peak loads were close in timing and magnitude to the measured loads data. However, results are very sensitive to data curve fitting and may not be suitable for Monte Carlo simulations. It was assumed that apparent mass was either negligible or a small fraction of enclosed mass, with little difference in results.

Ray, Eric S.↗

MODIS Solar Diffuser On-Orbit Degradation Characterization Using Improved SDSM Screen Modeling

The Solar Diffuser (SD) is used for the MODIS reflective solar bands (RSB) calibration. An on-board Solar Diffuser Stability Monitor (SDSM) tracks the degradation of its on-orbit bi-directional reflectance factor (BRF). To best match the SDSM detector signals from its Sun view and SD view, a fixed attenuation screen is placed in its Sun view path, where the responses show ripples up to 10%, much larger than design expectation. Algorithms have been developed since the mission beginning to mitigate the impacts of these ripples. In recent years, a look-up-table (LUT) based approach has been implemented to account for these ripples. The LUT modeling of the elevation and azimuth angles is constructed from the detector 9 (D9) of SDSM observations in the MODIS early mission. The response of other detectors is normalized to D9 to reduce the ripples observed in the sun-view data. The accuracy of all detectors degradation estimation depends on how well the D9 approximated. After multiple years of operation (Terra: 16 years; Aqua: 14 years), degradation behavior of all detectors can be monitored by their own. This paper revisits the LUT modeling and proposes a dynamic scheme to build a LUT independently for each detector. Further refinement in the Sun view screen characterization will be highlighted to ensure the degradation estimation accuracy. Results of both Terra and Aqua SD on-orbit degradation are derived from the improved modeling and curve fitting strategy.

Chen, H.↗

Reflectivity Around the Gold L-Edges of X-Ray Reflector of the Soft X-Ray Telescope Onboard ASTRO-H

We report the atomic scattering factor in the 11.215.4 keV for the ASTRO-H Soft X-ray Telescope (SXT) obtained in the ground based measurements. The large effective area of the SXT covers above 10 keV. In fact, the flight data show the spectra of the celestical objects in the hard X-ray band. In order to model the area, the reflectivity measurements in the 11.2-15.4 keV band with the energy pitch of 0.4-0.7 eV were made in the synchrotron beamline Spring-8 BL01B1. We obtained atomic scattering factors f1 and f2 by the curve fitting to the reflectivities of our witness sample. The edges associated with the golds L-I, II, and III transitions are identified, of which the depths are found to be roughly 60 percent shallower than those expected from the Henke's atomic scattering factor.

X-rays↗

Evaluation of Sentinel-1A Data For Above Ground Biomass Estimation in Different Forests in India

Use of remote sensing data for mapping and monitoring of forest biomass across large spatial scales can aid in addressing uncertainties in carbon cycle. Earlier, several researchers reported on the use of Synthetic Aperture Radar (SAR) data for characterizing forest structural parameters and the above ground biomass estimation. However, these studies cannot be generalized and the algorithms cannot be applied to all types of forests without additional information on the forest physiognomy, stand structure and biomass characteristics. The radar backscatter signal also saturates as forest parameters such as biomass and the tree height increase. It is also not clear how different polarizations (VV versus VH) impact the backscatter retrievals in different forested regions. Thus, it is important to evaluate the potential of SAR data in different landscapes for characterizing forest structural parameters. In this study, the SAR data from Sentinel-1A has been used to characterize forest structural parameters including the above ground biomass from tropical forests of India. Ground based data on tree density, basal area and above ground biomass data from thirty-eight different forested sites has been collected to relate to SAR data. After the pre-processing of Sentinel 1-A data for radiometric calibration, geo-correction, terrain correction and speckle filtering, the variability in the backscatter signal in relation tree density, basal area and above biomass density has been investigated. Results from the curve fitting approach suggested exponential model between the Sentinel-1A backscatter versus tree density and above ground biomass whereas the relationship was almost linear with the basal area in the VV polarization mode. Of the different parameters, tree density could explain most of the variations in backscatter. Both VV and VH backscatter signals could explain only thirty and thirty three percent of variation in above biomass in different forest sites of India. Results also suggested saturation of the Sentinel-1A backscatter signal around hundred tonnes per hectare for VV polarization and one hundred and forty five tonnes per hectare for VH polarization. The presentation will highlight the above results in addition to potentials and limitations of Sentinel-1A data for retrieving forest structural parameters. Also, background information on different forest types of India, biomass variations and forest type mapping efforts in the region will be presented.

Data↗

Diagnostic Techniques to Elucidate the Aerodynamic Performance of Acoustic Liners

In support of Topic A.2.8 of NASA NRA NNH10ZEA001N, the University of Florida (UF) has investigated the use of flow field optical diagnostic and micromachined sensor-based techniques for assessing the wall shear stress on an acoustic liner. Stereoscopic particle image velocimetry (sPIV) was used to study the velocity field over a liner in the Grazing Flow Impedance Duct (GFID). The results indicate that the use of a control volume based method to determine the wall shear stress is prone to significant error. The skin friction over the liner as measured using velocity curve fitting techniques was shown to be locally reduced behind an orifice, relative to the hard wall case in a streamwise plane centered on the orifice. The capacitive wall shear stress sensor exhibited a linear response for a range of shear stresses over a hard wall. PIV over the liner is consistent with lifting of the near wall turbulent structure as it passes over an orifice, followed by a region of low wall shear stress.

June, Jason↗

Unsteady Aerodynamic Force Sensing from Strain Data

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm.

aerodynamic load sensing↗