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

Method for Making Measurements of the Post-Combustion Residence Time in a Gas Turbine Engine

A system and method of measuring a residence time in a gas-turbine engine is provided, whereby the method includes placing pressure sensors at a combustor entrance and at a turbine exit of the gas-turbine engine and measuring a combustor pressure at the combustor entrance and a turbine exit pressure at the turbine exit. The method further includes computing cross-spectrum functions between a combustor pressure sensor signal from the measured combustor pressure and a turbine exit pressure sensor signal from the measured turbine exit pressure, applying a linear curve fit to the cross-spectrum functions, and computing a post-combustion residence time from the linear curve fit.

Miles, Jeffrey H↗

Implementation of Stochastic Gradient Descent in an Automated Glow Peak Identification Software for Multiple Thermoluminescent Dosimeter Types

A glow-curve analysis code was previously developed in C++ to analyze thermoluminescent dosimeter glow curves using automated peak detection while applying a first-order kinetics model. A newer version of this code was implemented to improve the automated peak detection and curve fitting models. The Stochastic Gradient Descent Algorithm was introduced to replace the prior approach of taking first and second-order derivatives for peak detection. Additionally, early stopping mechanisms were invoked to improve the previously used Levenberg-Marquardt Algorithm employed for curve fitting. The two software versions were compared through glow curve analysis of different thermoluminescent dosimeter materials and calculation of the corresponding figures of merit. Altogether improvements were shown, namely an increase in the number of peaks detected and a reduction of the mean figure of merit by approximately 46%.

137Cs↗

Quantitative XPS of plutonium: Evaluation of the Pu4f peak shape, relative sensitivity factors and estimated detection limits

We report that high-resolution X-ray photoelectron spectra have been acquired from sputter cleaned and in situ oxidized samples of α-plutonium and a face-centred cubic δ-plutonium–gallium alloy. The differences in the Pu4f peak shape in these two materials have been investigated, and the poorly screened satellite peaks have been quantified. Curve fitting parameters for the Pu4p 3/2 , Pu4d 5/2 , Pu4f and Pu6p 3/2 photoelectron peaks are reported, and relative sensitivity factors have been determined. The Pu4f curve fit model has been applied to data acquired using different spectrometers and alloys. Examples of quantification of the plutonium spectra are provided and minimum detection limits are calculated for common impurities in plutonium metal.

36 MATERIALS SCIENCE↗

Nonlinear Acoustical Assessment of Precipitate Nucleation

The purpose of the present work is to show that measurements of the acoustic nonlinearity parameter in heat treatable alloys as a function of heat treatment time can provide quantitative information about the kinetics of precipitate nucleation and growth in such alloys. Generally, information on the kinetics of phase transformations is obtained from time-sequenced electron microscopical examination and differential scanning microcalorimetry. The present nonlinear acoustical assessment of precipitation kinetics is based on the development of a multiparameter analytical model of the effects on the nonlinearity parameter of precipitate nucleation and growth in the alloy system. A nonlinear curve fit of the model equation to the experimental data is then used to extract the kinetic parameters related to the nucleation and growth of the targeted precipitate. The analytical model and curve fit is applied to the assessment of S' precipitation in aluminum alloy 2024 during artificial aging from the T4 to the T6 temper.

Cantrell, John H.↗

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗

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↗

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

Calibration improvements to electronically scanned pressure systems and preliminary statistical assessment

Orifice-to-orifice inconsistencies in data acquired with an electronically-scanned pressure system at the beginning of a wind tunnel experiment forced modifications to the standard, instrument calibration procedures. These modifications included a large increase in the number of calibration points which would allow a critical examination of the calibration curve-fit process, and a subsequent post-test reduction of the pressure data. Evaluation of these data has resulted in an improved functional representation of the pressure-voltage signature for electronically-scanned pressures sensors, which can reduce the errors due to calibration curve fit to under 0.10 percent of reading compared to the manufacturer specified 0.10 percent of full scale. Application of the improved calibration function allows a more rational selection of the calibration set-point pressures. These pressures should be adjusted to achieve a voltage output which matches the physical shape of the pressure-voltage signature of the sensor. This process is conducted in lieu of the more traditional approach where a calibration pressure is specified and the resulting sensor voltage is recorded. The fifteen calibrations acquired over the two-week duration of the wind tunnel test were further used to perform a preliminary, statistical assessment of the variation in the calibration process. The results allowed the estimation of the bias uncertainty for a single instrument calibration; and, they form the precursor for more extensive and more controlled studies in the laboratory.

Everhart, Joel L.↗

Measurement of Residual Flexibility for Substructures Having Prominent Flexible Interfaces

Verification of a dynamic model of a constrained structure requires a modal survey test of the physical structure and subsequent modification of the model to obtain the best agreement possible with test data. Constrained-boundary or fixed-base testing has historically been the most common approach for verifying constrained mathematical models, since the boundary conditions of the test article are designed to match the actual constraints in service. However, there are difficulties involved with fixed-base testing, in some cases making the approach impractical. It is not possible to conduct a truly fixed-base test due to coupling between the test article and the fixture. In addition, it is often difficult to accurately simulate the actual boundary constraints, and the cost of designing and constructing the fixture may be prohibitive. For use when fixed-base testing proves impractical or undesirable, alternate free-boundary test methods have been investigated, including the residual flexibility technique. The residual flexibility approach has been treated analytically in considerable detail and has had limited frequency response measurements for the method. This concern is well-justified for a number of reasons. First, residual flexibilities are very small numbers, typically on the order of 1.0E-6 in/lb for translational diagonal terms, and orders of magnitude smaller for off-diagonal values. This poses difficulty in obtaining accurate and noise-free measurements, especially for points removed from the excitation source. A second difficulty encountered in residual measurements lies in obtaining a clean residual function in the process of subtracting synthesized modal data from a measured response function. Inaccuracies occur since modes are not subtracted exactly, but only to the accuracy of the curve fits for each mode; these errors are compounded with increasing distance from the excitation point. In this paper, the residual flexibility method is applied to a simple structure in both test and analysis. Measured and predicted residual functions are compared, and regions of poor data in the measured curves are described. It is found that for accurate residual measurements, frequency response functions having prominent stiffness lines in the acceleration/force format are needed. The lack of such stiffness lines increases measurement errors. Interface drive point frequency respose functions for shuttle orbiter payloads exhibit dominant stiffness lines, making the residual test approach a good candidate for payload modal tests when constrained tests are inappropriate. Difficulties in extracting a residual flexibility value from noisy test data are discussed. It is shown that use of a weighted second order least-squares curve fit of the measured residual function allows identification of residual flexibility that compares very well with predictions for the simple structure. This approach also provides an estimate of second order residual mass effects.

Tinker, Michael L.↗

Multiwavelength Pyrometry To Correct For Reflections

Computerized curve fitting yields more and better information on thermal radiation. Multiwavelength pyrometry involves measurement of spectrum of thermal radiation emitted by and reflected from specimen. Auxiliary source enables determination of spectral reflectance of specimen. Spectral reflectance and temperature of specimen obtained from spectral measurements by nonlinear least-squares curve-fitting routine.

Ng, Daniel↗

Column atmospheric water vapor retrievals from airborne imaging spectrometer data

High-spatial-resolution column atmospheric water vapor amounts were derived from spectral data collected by the Airborne Visible Infrared Imaging Spectrometer (AVIRIS). The quantitative derivation is made by curve fitting observed spectra with calculated spectra in the 1.14- and 0.94-micron water-vapor band absorption regions with a nonlinear least-squares technique. The precision of the retrieved column water vapor is approximately 5 percent. The derived column water vapor amounts are independent of the absolute surface reflectance. Curve fitting of spectra near 1 micron from areas covered with vegetation indicates that both the amount of atmospheric water vapor and the moisture content of vegetation can be retrieved simultaneously. It should be possible to measure column water vapor over land areas from satellite altitude with the proposed high-resolution imaging spectrometer or even the moderate-resolution imaging spectrometer.

Gao, Bo-Cai↗

Status of Cycle 23 Forecasts

Forecasts for the amplitude of cycle 23 that were reported prior to the start of, the cycle covered a full range of, values from very small to very large. A forecast reached by the consensus of a panel of forecasters convened at the time of minimum in 1996 [Joselyn et al., 1997] suggested that this cycle would be much larger than average with the smoothed International Sunspot Number reaching a maximum of 160 +/- 30 in the middle of the year 2000. A recent survey of solar cycle prediction techniques [Hathaway et al., 1999] found that the two most reliable techniques for forecasting the cycle prior to its start give similar predictions for this cycle's maximum -- 154 +/- 26 and 153 +/- 33. Curve-fitting and regression techniques can be used with some confidence now that cycle 23 is well underway. These techniques indicate a more modest sunspot cycle with a maximum of 112 + 10 -- only slightly larger than average. The current (May 2000) prediction using the combined predictions from both precursors and curve-fitting gives a cycle amplitude of about 135 +/- 20. This is within the errors given by the consensus and the precursor technique predictions but very close to their lower bounds. Inspection of the characteristics of cycle 23 thus far indicates that in some respects this cycle is an outlier but in other respects it is quite normal.

Hathaway, David H.↗

Machine Learning Calabi–Yau Metrics

We apply machine learning to the problem of finding numerical Calabi–Yau metrics. Building on Donaldson's algorithm for calculating balanced metrics on Kähler manifolds, we combine conventional curve fitting and machine-learning techniques to numerically approximate Ricci-flat metrics. We show that machine learning is able to predict the Calabi–Yau metric and quantities associated with it, such as its determinant, having seen only a small sample of training data. Using this in conjunction with a straightforward curve fitting routine, we demonstrate that it is possible to find highly accurate numerical metrics much more quickly than by using Donaldson's algorithm alone, with our new machine-learning algorithm decreasing the time required by between one and two orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wing Shape Sensing from Measured Strain

A new two-step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least-squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to interpolate and extrapolate the deflection and slope of the entire structure through the use of the System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular plate wing. The theory is then applied to test data from a cantilevered swept-plate wing model. Computed results are compared with finite element results, results using another strainbased method, and photogrammetry data. For the computational model under an aeroelastic load, maximum deflection errors in the fore and aft, lateral, and vertical directions are -3.2%, 0.28%, and 0.09%, respectively; and maximum slope errors in roll and pitch directions are 0.28% and -3.2%, respectively. For the experimental model, deflection results at the tip are shown to be accurate to within 3.8% of the photogrammetry data and are accurate to within 2.2% in most cases. In general, excellent matching between target and computed values are accomplished in this study. Future refinement of this theory will allow it to monitor the deflection and health of an entire aircraft in real time, allowing for aerodynamic load computation, active flexible motion control, and active induced drag reduction.

shape sensing↗

Union through UNITY: Cosmology with 2000 SNe Using a Unified Bayesian Framework

Type Ia supernovae (SNe Ia) were instrumental in establishing the acceleration of the Universe’s expansion. By virtue of their combination of distance reach, precision, and prevalence, they continue to provide key cosmological constraints, complementing other cosmological probes. Individual SN surveys cover only over about a factor of 2 in redshift, so compilations of multiple SN data sets are strongly beneficial. We assemble an up-to-date “Union” compilation of 2087 cosmologically useful SNe Ia from 24 data sets (“Union3”). We take care to put all SNe on the same distance scale and update the light-curve fitting with SALT3 to use the full rest-frame optical. Over the next few years, the number of cosmologically useful SNe Ia will increase by more than a factor of 10, and keeping systematic uncertainties subdominant will be more challenging than ever. We discuss the importance of treating outliers, selection effects, light-curve shape/color populations/standardization relations, unexplained dispersion, and heterogeneous observations simultaneously. We present an updated Bayesian framework, called UNITY1.5 (Unified Nonlinear Inference for Type-Ia cosmologY), that incorporates significant improvements in our ability to model selection effects, standardization, and systematic uncertainties compared to earlier analyses. As an analysis byproduct, we also recover the posterior of the SN-only peculiar-velocity field, although we do not interpret it in this work. We compute updated cosmological constraints with Union3 and UNITY1.5, finding weak 1.7σ–2.6σ tension with flat cold dark matter and possible evidence for thawing dark energy (w0 > − 1, wa < 0). We release our SN distances, light-curve fits, and UNITY1.5 framework to the community.

Rubin, David↗

Frequency domain state-space system identification

An algorithm for identifying state-space models from frequency response data of linear systems is presented. A matrix-fraction description of the transfer function is employed to curve-fit the frequency response data, using the least-squares method. The parameters of the matrix-fraction representation are then used to construct the Markov parameters of the system. Finally, state-space models are obtained through the Eigensystem Realization Algorithm using Markov parameters. The main advantage of this approach is that the curve-fitting and the Markov parameter construction are linear problems which avoid the difficulties of nonlinear optimization of other approaches. Another advantage is that it avoids windowing distortions associated with other frequency domain methods.

Chen, Chung-Wen↗

Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

Johnston, Christopher O.↗

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗