Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized least squared”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Calibration and evaluation of Skylab altimetry for geodetic determination of the geoid

The author has identified the following significant results. The Skylab altimeter experiment has proven the capability of the altimeter for measurement of sea surface topography. The geometric determination of the geoid/mean sea level from satellite altimetry is a new approach having significant applications in many disciplines including geodesy and oceanography. A generalized least squares collocation technique was developed for determination of the geoid from altimetry data. The technique solves for the altimetry geoid and determines one bias term for the combined effect of sea state, orbit, tides, geoid, and instrument error using sparse ground truth data. The influence of errors in orbit and a priori geoid values are discussed. Although the Skylab altimeter instrument accuracy is about plus or minus 1m, significant results were obtained in identification of large geoidal features such as over the Puerto Rico trench. Comparison of the results of several passes shows that good agreement exists between the general slopes of the altimeter geoid and the ground truth, and that the altimeter appears to be capable of providing more details than are now available with best known geoids.

Mourad, A. G.↗
DOE OSTI2022

An influence of manufacturing tolerances on Pin-Cell k-infinity of MOX fuel using data from the FUBILA experiment program

An influence of manufacturing tolerances (MTs) was evaluated on the pin-cell k-infinity of MOX fuel through the random sampling of CASMO5 calculations. The data from the FUBILA experiment program was used as the manufacturing parameters. For the uncertainties of element/isotope mass fractions, their covariance matrices were calculated by the generalized least square method. The total k-infinity uncertainty was 120-250 pcm (percent mille). From the breakdown of k-infinity uncertainty per materials, the MTs of fuel pellet had a dominant influence. The individual influences were also evaluated for element/isotope mass fractions, an inner/outer diameter, and a density. Those of element/isotope mass fractions were less than a few dozen pcm. Since the perturbations of fuel pellet diameter and density caused large variations to the total amount of heavy metals and effect of spatial self-shielding, they had a negative correlation with the k-infinity perturbation. On the AG3 (Al-Mg alloy) over-cladding, the inner/outer diameter perturbations had both negative and positive correlations. The negative one was due to decrease in amount of light water and the positive one was due to increase in amount of AG3. Since the effect of neutron slowing-down by light water had a dominant influence, the influence of the former negative correlation was larger than that of the latter positive correlation. It is, therefore, concluded that the consideration of the MTs that had a large influence on the effect of neutron slowing-down is important for precise quantification of the pin-cell k-infinity uncertainty of MOX fuel. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A semiparametric latent factor model for large scale temporal data with heteroscedasticity

Large scale temporal data have flourished in a vast array of applications, and their sophisticated structures, especially the heteroscedasticity among subjects with inter- and intra-temporal dependence, have fueled a great demand for new statistical models. In this paper, with covariate information, we consider a flexible model for large scale temporal data with subject-specific heteroscedasticity. Formally, the model employs latent semiparametric factors to simultaneously account for the subject-specific heteroscedasticity and the contemporaneous and/or serial correlations. The subject-specific heteroscedasticity is modeled as the product of the unobserved factor process and subject’s covariate effect, which is further characterized via additive models. For estimation, we propose a two-step procedure. First, the latent factor process and nonparametric loading are recovered through projection-based methods, and following, we estimate the regression components by approaches motivated from the generalized least squares. By scrupulously examining the non-asymptotic rates for recovering the factor process and its loading, we show the consistency and efficiency of estimated regression coefficients in the absence of prior knowledge of latent factor process and subject’s covariate effect. Here, the statistical guarantees remain valid even for finite time points that makes our method particularly appealing when the subjects significantly outnumber the observation time points. Using comprehensive simulations, we demonstrate the finite sample performance of our method, which corroborates the theoretical findings. Finally, we apply our method to a data set of air quality and energy consumption collected at 129 monitoring sites in the United States in 2015.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Evaluating 239 Pu(n,f) cross sections via machine learning using experimental data, covariances, and measurement features

In this paper, the neutron-induced 239 Pu fission cross section, 239 Pu(n,f), is evaluated from 1–20 MeV using experimental data and associated covariances while also considering information on the measurement, termed features here. For instance, methods to determine the background, sample backing material, or impurities in the sample, are explicitly taken into account in the evaluation process. To this end, outliers in the experimental data are identified with a modified version of the Hybrid Robust Support Vector Machine. In a second step, two machine learning methods (logistic regression with elastic net regularization and random forest regression with SHAP feature importance metric) are used to highlight measurement features that are common among many of the outlying data points. Based on this analysis, penalty uncertainties are added to the experimental covariances of outlying data points that have outlier measurement features and are put through the generalized-least-squares evaluation. The resulting evaluated mean values and covariances differ distinctly from those data evaluated without the penalty uncertainties. These results highlight that certain measurement features should be more closely examined.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Thermal adaptation of soil microbial growth traits in response to chronic warming

ABSTRACT Adaptation of soil microbes due to warming from climate change has been observed, but it remains unknown what microbial growth traits are adaptive to warming. We studied bacterial isolates from the Harvard Forest Long-Term Ecological Research site, where field soils have been experimentally heated to 5°C above ambient temperature with unheated controls for 30 years. We hypothesized that Alphaproteobacteria from warmed plots have (i) less temperature-sensitive growth rates; (ii) higher optimum growth temperatures; and (iii) higher maximum growth temperatures compared to isolates from control plots. We made high-throughput measurements of bacterial growth in liquid cultures over time and across temperatures from 22°C to 37°C in 2–3°C increments. We estimated growth rates by fitting Gompertz models to the growth data. Temperature sensitivity of growth rate, optimum growth temperature, and maximum growth temperature were estimated by the Ratkowsky 1983 model and a modified Macromolecular Rate Theory (MMRT) model. To determine evidence of adaptation, we ran phylogenetic generalized least squares tests on isolates from warmed and control soils. Our results showed evidence of adaptation of higher optimum growth temperature of bacterial isolates from heated soils. However, we observed no evidence of adaptation of temperature sensitivity of growth and maximum growth temperature. Our project begins to capture the shape of the temperature response curves, but illustrates that the relationship between growth and temperature is complex and cannot be limited to a single point in the biokinetic range. IMPORTANCE Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Bayesian Monte Carlo Evaluation Framework for Imperfect Nuclear Data

Bayesian evaluation of resolved resonance region (RRR) nuclear data has historically been carried out using the generalized least squares (GLS) formalism, as implemented in, e.g., SAMMY. We have recently developed a prototype of Bayesian Monte Carlo (BMC) evaluation framework, implemented using a Markov Chain Monte Carlo (MCMC) method with a Metropolis-Hastings (MH) acceptance criterion. This was done in order to remove the approximations underlying the conventional GLS evaluations, namely, the linear approximation, and the approximation that all probability density functions (PDFs) are of the normal kind. Recent works by others have used similar stochastic approaches to quantify cross section uncertainties from ENDF evaluated co-variances, and/or, from integral benchmark data, but those have not been conceived as an evaluation framework like the one presented here.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Temporal deposition of copper and zinc in the sediments of metal removal constructed wetlands

The objective of this study was to explore the effects of time, seasons, and total carbon (TC) on Copper (Cu) and Zinc (Zn) deposition in the surface sediments. This study was performed at the H-02 constructed wetland on the Savannah River Site (Aiken, SC, USA). Covering both warm (April-September) and cool (October-March) seasons, several sediment cores were collected twice a year from the H-02 constructed wetland cells from 2007 to 2013. Total concentrations of Cu and Zn were measured in the sediments. Concentrations of Cu and Zn (mean ± standard deviation) in the surface sediments over 7 years of operation increased from 6.0 ± 2.8 and 14.6 ± 4.5 mg kg -1 to 139.6 ± 87.7 and 279.3 ± 202.9 mg kg -1 dry weight, respectively. The linear regression model explained the behavior and the variability of Cu deposition in the sediments. On the other hand, using the generalized least squares extension with the linear regression model allowed for unequal variance and thus produced a model that explained the variance properly, and as a result, was more successful in explaining the pattern of Zn deposition. Total carbon significantly affected both Cu ( p = 0.047) and Zn ( p < 0.001). Time effect on Cu deposition was statistically significant ( p = 0.013), whereas Zn was significantly affected by the season ( p = 0.009).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A New Probability Table Methodology with Covariances on Probability Tables

Through this work, we test a new methodology of evaluating probability tables informed by a combination of physics models and experimental data. This methodology uses generalized least squares to produce a joint mean and covariance for average cross sections and probability tables. This methodology has been tested and validated using total Monte Carlo. Attempts were made to parameterize the probability table using a formulation provided in. Although the normal inverse Gaussian fit was successful for the total cross section probability distribution, the rational-function fit to the capture, inelastic, and fission cross sections was unsuccessful due to time constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

A generalized technique for using cones and dihedral angles in attitude determination, revision 1

Analytic development is presented for a general least squares attitude determination subroutine applicable to spinning satellites. The method is founded on a geometric approach which is completely divorced from considerations relating to particular types and configurations of onboard attitude sensors. Any mix of sensor measurements which can be first transformed (outside the program) to cone or dihedral angle data can be processed. A cone angle is an angle between the spin axis and a known direction line in space; a dihedral angle is an angle between two planes formed by the spin axis and each of two known direction lines. Many different kinds of sensor data can be transformed to these angles, which in turn constitute the actual program inputs, so that the subroutine can be applied without change to a variety of satellite missions. Either a constant or dynamic spin axis model can be handled. The program is also capable of solving for fixed biases in the input angles, in addition to the spin axis attitude solution.

Werking, R. D.↗
NASA NTRS1975

A random search algorithm for laboratory computers

The small laboratory computer is ideal for experimental control and data acquisition. Postexperimental data processing is often performed on large computers because of the availability of sophisticated programs, but costs and data compatibility are negative factors. Parameter optimization can be accomplished on the small computer, offering ease of programming, data compatibility, and low cost. A previously proposed random-search algorithm ('random creep') was found to be very slow in convergence. A method is proposed (the 'random leap' algorithm) which starts in a global search mode and automatically adjusts step size to speed convergence. A FORTRAN executive program for the random-leap algorithm is presented which calls a user-supplied function subroutine. An example of a function subroutine is given which calculates maximum-likelihood estimates of receiver operating-characteristic parameters from binary response data. Other applications in parameter estimation, generalized least squares, and matrix inversion are discussed.

Curry, R. E.↗
NASA NTRS1977

Improved ground truth geoid for the GEOS-3 calibration area

The purpose of this investigation is to develop methods and procedures are reported for computing a detailed geoid to be used as geodetic ground truth for the calibration and verification of GEOS-3 altimeter data. The technique developed is based on rectifying the best available detailed geoid so that the rectified geoid will have correct scale, orientation, shape and position with respect to the geocenter. The approach involved the development of a mathematical model based on a second degree polynomial, in rectangular Cartesian coordinates, describing the geoid undulations at the control stations. A generalized least squares solution was obtained for the polynomial which describes the variation of the undulation differences between the control stations geoid and the gravimetric geoid. Three rectified geoid were determined. These geoids correspond to three sets of tracking station data: (1) WFC/C-band data; (2) GSFC/C-band data; and (3) OSU-275 data. The absolute accuracy of these rectified geoids is linearly correlated with the uncertainties of the tracking station coordinates and, to a certain extent, with those of the detailed geoid being rectified.

Mourad, A. G.↗
NASA NTRS1979

Error in total ozone measurements arising from aerosol attenuation

A generalized least squares method for deducing both total ozone and aerosol extinction spectrum parameters from Dobson spectrophotometer measurements was developed. An error analysis applied to this system indicates that there is little advantage to additional measurements once a sufficient number of line pairs have been employed to solve for the selected detail in the attenuation model. It is shown that when there is a predominance of small particles (less than about 0.35 microns in diameter) the total ozone from the standard AD system is too high by about one percent. When larger particles are present the derived total ozone may be an overestimate or an underestimate but serious errors occur only for narrow polydispersions.

Thomas, R. W. L.↗
NASA NTRS1986

A recursive solution for a fading memory filter derived from Kalman filter theory

A simple recursive solution for a class of fading memory tracking filters is presented. A fading memory filter provides estimates of filter states based on past measurements, similar to a traditional Kalman filter. Unlike a Kalman filter, an exponentially decaying weight is applied to older measurements, discounting their effect on present state estimates. It is shown that Kalman filters and fading memory filters are closely related solutions to a general least squares estimator problem. Closed form filter transfer functions are derived for a time invariant, steady state, fading memory filter. These can be applied in loop filter implementation of the Deep Space Network (DSN) Advanced Receiver carrier phase locked loop (PLL).

Statman, J. I.↗
NASA NTRS2002

Experimental Investigations of Non-Stationary Properties In Radiometer Receivers Using Measurements of Multiple Calibration References

Radiometers must be periodically calibrated because the receiver response fluctuates. Many techniques exist to correct for the time varying response of a radiometer receiver. An analytical technique has been developed that uses generalized least squares regression (LSR) to predict the performance of a wide variety of calibration algorithms. The total measurement uncertainty including the uncertainty of the calibration can be computed using LSR. The uncertainties of the calibration samples used in the regression are based upon treating the receiver fluctuations as non-stationary processes. Signals originating from the different sources of emission are treated as simultaneously existing random processes. Thus, the radiometer output is a series of samples obtained from these random processes. The samples are treated as random variables but because the underlying processes are non-stationary the statistics of the samples are treated as non-stationary. The statistics of the calibration samples depend upon the time for which the samples are to be applied. The statistics of the random variables are equated to the mean statistics of the non-stationary processes over the interval defined by the time of calibration sample and when it is applied. This analysis opens the opportunity for experimental investigation into the underlying properties of receiver non stationarity through the use of multiple calibration references. In this presentation we will discuss the application of LSR to the analysis of various calibration algorithms, requirements for experimental verification of the theory, and preliminary results from analyzing experiment measurements.

Racette, Paul↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗