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 109 records · Page 6

Simulation of Delamination Propagation in Composites Under High-Cycle Fatigue by Means of Cohesive-Zone Models

A damage model for the simulation of delamination propagation under high-cycle fatigue loading is proposed. The basis for the formulation is a cohesive law that links fracture and damage mechanics to establish the evolution of the damage variable in terms of the crack growth rate dA/dN. The damage state is obtained as a function of the loading conditions as well as the experimentally-determined coefficients of the Paris Law crack propagation rates for the material. It is shown that by using the constitutive fatigue damage model in a structural analysis, experimental results can be reproduced without the need of additional model-specific curve-fitting parameters.

Turon, Albert↗

Observations on Cost Modeling and Performance Measurement of Long Term Archives

This paper describes a prototype suite of Excel-based tools that could be used for estimating lifecycle costs for newly planned or modified long-term archival facilities. These tools may also prove valuable for monitoring the long-term performance of such facilities once operational. Cost estimation is by analogy, using statistical curve-fitting techniques across a database of comparable data activities. The database currently includes 29 operational data centers ranging from small (2 FTEs) to large (66 FTEs), and is readily expandable to include additional activities specifically involving data preservation and added value. Each comparable data center is described in terms of its staffing, throughput workload, archival and distribution requirements, levels of user service, overall complexity, degree of automation, and other data, comprising 94 distinct descriptors in all. The descriptors were developed by normalizing heterogeneous data from the various centers and mapping them into an Excel framework consistent with the OAIS reference model. A user-friendly tool is provided for generating input to and updating the comparables database. This tool can also be used to benchmark the performance (in terms of cost versus throughput) of an operational data center, and to update the staffing, cost and workload data on a periodic basis. The comparables database could thus provide a history of staffing and throughput over time, as a means of performance monitoring and providing feedback for continuous improvement. Ancillary tools are also provided for performing "what-if' cost exercises for planning purposes, and for graphical display of data and results. We provide a high-level description of the tools; present our experiences and observations on gathering the information and maintaining the database; and discuss how this tool set might be applied to long term archives.

Fontaine, Kathy↗

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗

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.↗

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↗

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↗

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↗

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.↗

Characterization of Collection Efficiency of the Common Research Model Midspan Wing Section in the IRT

This paper presents a preliminary study for the characterization of collection efficiency from icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center. A test method previously developed for measuring the attachment line maximum collection efficiency of a swept NACA 0012 airfoil model at zero angle of attack was applied to the leading-edge region of a 65%-scale version of the Common Research Model midspan wing section. A correlation for the stagnation line maximum collection efficiency as a function of the modified inertia parameter was obtained with LEWICE3D simulations utilizing a discrete number of drop diameters. It was then compared with the collection efficiency measurement data obtained in the IRT. For the experimental collection efficiency, two ice shape digitization procedures were utilized to extract 2-D chord-wise ice shape profiles, i.e., the Maximum Combined Cross Section or MCCS, and the Minimum Combined Cross Section or Min CCS, at selected span-wise locations from the 3-D scanned ice shapes. The preliminary result showed that the rime ice thickness method can be used to characterize the collection efficiency distribution, for conditions free of ice erosion, in the main ice shape region. However, a high-order statistical processing of the ice scan is needed for estimating the mean ice thickness in areas where feathers are prevalent. From the limited comparison of the experimental and LEWICE3D collection efficiency data in the CRM65 MS model leading edge area, it was shown that the attachment line maximum collection efficiency is reasonably estimated by the best curve-fit correlation in the range of modified inertia parameter tested. A tighter alignment control of the iced and cold clean model scans is needed to improve the comparison. Further evaluation of this correlation is recommended to assess its applicability for Common Research Model type swept wing icing scaling analysis.

Ice Scailing↗

Characterization of Collection Efficiency of the Common Research Model Midspan Wing Section in the IRT

This paper presents a preliminary study for the characterization of collection efficiency from icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center. A test method previously developed for measuring the attachment line maximum collection efficiency of a swept NACA 0012 airfoil model at zero angle of attack was applied to the leading-edge region of a 65%-scale version of the Common Research Model midspan wing section. A correlation for the stagnation line maximum collection efficiency as a function of the modified inertia parameter was obtained with LEWICE3D simulations utilizing a discrete number of drop diameters. It was then compared with the collection efficiency measurement data obtained in the IRT. For the experimental collection efficiency, two ice shape digitization procedures were utilized to extract 2-D chord-wise ice shape profiles, i.e., the Maximum Combined Cross Section or MCCS, and the Minimum Combined Cross Section or MinCCS, at selected span-wise locations from the 3-D scanned ice shapes. The preliminary result showed that the rime ice thickness method can be used to characterize the collection efficiency distribution, for conditions free of ice erosion, in the main ice shape region. However, a high-order statistical processing of the ice scan is needed for estimating the mean ice thickness in areas where feathers are prevalent. From the limited comparison of the experimental and LEWICE3D collection efficiency data in the CRM65 MS model leading edge area, it was shown that the attachment line maximum collection efficiency is reasonably estimated by the best curve-fit correlation in the range of modified inertia parameter tested. A tighter alignment control of the iced and cold clean model scans as well as the model angle of attack are needed to improve the comparison. Further evaluation of this correlation is recommended to assess its applicability for Common Research Model type swept wing icing scaling analysis.

Icing Scaling↗

Improved Methods for Hyperelastic Material Characterizations Based on Uniaxial, Radial-Biaxial, and Strip-Biaxial Test Data

This report concerns the different analytical methods used for material characterizations of hyperelastic silicone rubber, a potential candidate elastomer for aerospace applications. To establish the constitutive equations for the silicon rubber, strain-energy gradients (material constants) must be determined from experiments. Earlier experiments by Blatz-Ko showed that the material constants of natural rubber determined from different stress fields were different. Therefore, in the current material characterization, uniaxial, equalbiaxial, and strip-biaxial tensile tests were conducted. Based on the resulting stress-stretch data, the material constants were determined by using a trendline on the stress-difference plots and energy plots, and also the MSC-Marc ® curve-fitting method. The results confirm that the material constants determined from uniaxial and biaxial stress fields are indeed different. The material constants determined from the equal-biaxial energy plots and from the MSC-Marc ® combined-curve-fitting method were found to be extremely close, and can be considered as universal material constants for establishing the constitutive equations for the silicone rubber under any stress field. The results of this report show that a uniaxially stretched specimen became slightly anisotropic due to crystallization of polymer molecular chains, and that equal-biaxial and strip-biaxial testing must be conducted to fully understand the material properties of hyperelastic materials.

William L Ko↗

Development and Validation of an Initial Electric Motor Rotor Vibration Model

An overview of the development of a low fidelity prediction tool for electric motor rotor resonant modes and frequencies. These predictions are useful for the prediction of UAM motor noise during the design phase. Finite Element Analysis of motor rotors has led to the development of a parameterized model to predict these modes shapes. This has in turn allowed to a curve-fit based model from parametric sweeps. These curves can predict resonant frequencies and shapes for 1-5kW class motors without running a simulation.

acoustics↗

Characterization of Large Drop Velocity in the NASA Icing Research Tunnel

This paper presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the test section. Some icing spray clouds with large drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow in the measurement area of the probe between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a pitot-static probe affixed on the mounting plate next to the stand. Those findings were later verified in the Icing Research Tunnel using a new pitot-static probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations were performed with the airflow moving towards the probe head with a Langmuir-D 7-bin distribution cloud containing large drops. The simulation results helped identify a critical drop-size threshold of 300 μm above which the velocities of larger drops are negligibly affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.

SLD↗

Comparison and Application of Battery Modeling Methods for Conceptual Electric Aircraft

Growing interest in design and optimization of electrified aircraft propulsion concepts prompts the need for accurate, flexible, and efficient methods to model battery systems. Presented in this paper are three battery modeling methods that have been used at NASA’s Glenn Research Center, each representing different mathematical or electrical approaches. Thévenin equivalent circuit, normalization, and curve-fitting methods are compared against battery cell test data for the X-57 Maxwell electric aircraft technology demonstrator. The methods are then applied in a simple multidisciplinary optimization context using NASA’s Six-Passenger Electric Quadrotor concept to determine their applicability and performance. The normalization method achieves the highest accuracy for steady and unsteady discharge rates with a voltage mean error percentage of 0.423% and 1.186%, respectively. Optimal quadrotor mission range between the models varies up to 0.5 nmi, identifying current battery modeling methods as a potentially significant contributor to mission analysis error. A set of relevant tools and techniques for conceptual battery modeling are identified in this paper, with conclusions made on the utility of each modeling approach for various design challenges.

power↗

Characterization of Large Drop Velocity in the NASA Icing Research Tunnel

This presentation presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the tunnel test section. Some icing spray clouds with large-sized drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a Pitot static probe affixed on the mounting plate. Those findings were later verified in the Icing Research Tunnel using a new Pitot probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations as airflow moving towards the probe head were run for a Langmuir-D 7-bin drop size distribution of a spray cloud with a nominally large value of medium volumetric diameter. The simulation results helped identify a critical drop-size threshold of 300 µm above which the velocities of larger drops are not affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.

SLD↗