Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional regression”

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

Assessment of Calibration Approaches for the Stress Relaxation Test

Advanced manufacturing technologies are enabling the next generation of superalloys for extreme temperatures, pressures, and longevity. There is a need for Accelerated Creep Test (ACT) methods to reduce the time needed to implement new creep resistant superalloys. Without the development of ACTs, the qualification of new superalloys using conventional creep tests can take 11+ years corresponding to the 100,000-hour service life of long-lived turbomachinery components. A refined stress relaxation test (SRT) is introduced to succeed real-time conventional creep tests (CCTs). The stress relaxation test generates extensive stress versus creep-strain-rate data within a single specimen. Since the specimen is held below the elastic limit and relatively little creep deformation has occurred, several SRTs at different temperatures can be performed using a single specimen. To obtain data, specimens are rapidly loaded in force-control to just below the elastic limit. The specimen is switched to displacement-control and held at a fixed displacement while stress relaxation occurs. This generates a stress versus time curve that must be converted into stress versus creep-strain-rate. The stress versus creep-strain-rate curves are utilized to calculate creep activation energy and produce a creep-strain-rate master curve. The objective of this study is to determine the best calibration approach for SRT: regression of stress, regression of creep strain, and the finite difference of creep strain. It is recommended that the differentiation be applied; however, finite difference applied to the raw data produces inconsistent calculations of the creep-strain-rate. In this study, a MATLAB algorithm is written to evaluate all three calibration options where, stress and creep deformation are regressed into natural logarithm functions. The resulting creep-strain-rate curves are compared to determine the most reliable and consistent method. Analyzing and evaluating the results the regression of stress and regression or creep strain calibration methods prove to provide accurate results. These methods will help replace CCTs by conducting less expensive ACTs, while increasing the overall efficiency of the development of new materials.

36 MATERIALS SCIENCE↗

On least squares collocation

It is shown that the least squares collocation approach to estimating geodetic parameters is identical to conventional minimum variance estimation. Hence, the least squares collocation estimator can be derived either by minimizing the usual least squares quadratic loss function or by computing a conditional expectation by means of the regression equation. When a deterministic functional relationship between the data and the parameters to be estimated is available, one can implement a least squares solution using the functional relation to obtain an equation of condition. It is proved the solution so obtained is identical to what is obtained through least squares collocation. The implications of this equivalance for the estimation of mean gravity anomalies are discussed.

Argentiero, P. D.↗

A Corrected Score Function Framework for Modelling Circadian Gene Expression

Many biological processes display oscillatory behaviour based on an approximately 24 h internal timing system specific to each individual. One process of particular interest is gene expression, for which several circadian transcriptomic studies have identified associations between gene expression during a 24 h period and an individual's health. A challenge with analysing data from these studies is that each individual's internal timing system is offset relative to the 24 h day-night cycle, where day–night cycle time is recorded for each collected sample. Laboratory procedures can accurately determine each individual's offset and determine the internal time of sample collection. However, these laboratory procedures are labour-intensive and expensive. Here, in this paper, we propose a corrected score function framework to obtain a regression model of gene expression given internal time when the offset of each individual is too burdensome to determine. A feature of this framework is that it does not require the probability distribution generating offsets to be symmetric with a mean of zero. Simulation studies validate the use of this corrected score function framework for cosinor regression, which is prevalent in circadian transcriptomic studies. Illustrations with data from three circadian transcriptomic studies further demonstrate that the proposed framework consistently mitigates bias relative to using a score function that does not account for this offset.

59 BASIC BIOLOGICAL SCIENCES↗

Design and analysis of machine learning exchange-correlation functionals via rotationally invariant convolutional descriptors

Here, we explore the potential of a data-driven approach to the design of exchange-correlation (xc) functionals. The approach, inspired by convolutional filters in computer vision and surrogate functions from optimization, utilizes convolutions of the electron density to form a feature space to represent local electronic environments and neural networks to map the features to the exchange-correlation energy density. These features are orbital free, and provide a systematic route to including information at various length scales. This work shows that convolutional descriptors are theoretically capable of an exact representation of the electron density, and proposes Maxwell-Cartesian spherical harmonic kernels as a class of rotationally invariant descriptors for the construction of machine learned functionals. The approach is demonstrated using data from the B3LYP functional on a number of small molecules containing C, H, O, and N along with a neural network regression model. The machine learned functionals are compared to standard physical approximations and the accuracy is assessed for the absolute energy of each molecular system as well as formation energies. The results indicate that it is possible to reproduce the exchange-correlation portion of B3LYP formation energies to within chemical accuracy using orbital-free descriptors with a spatial extent of 0.2 Å. The findings provide empirical insight into the spatial range of electron exchange, and suggest that the combination of convolutional descriptors and machine learning regression models is a promising framework for xc functional design, although challenges remain in obtaining training data and generating models consistent with pseudopotentials.

36 MATERIALS SCIENCE↗

Molecular dipole moment learning via rotationally equivariant derivative kernels in molecular-orbital-based machine learning

This study extends the accurate and transferable molecular-orbital-based machine learning (MOB-ML) approach to modeling the contribution of electron correlation to dipole moments at the cost of Hartree–Fock computations. A MOB pairwise decomposition of the correlation part of the dipole moment is applied, and these pair dipole moments could be further regressed as a universal function of MOs. The dipole MOB features consist of the energy MOB features and their responses to electric fields. An interpretable and rotationally equivariant derivative kernel for Gaussian process regression (GPR) is introduced to learn the dipole moment more efficiently. The proposed problem setup, feature design, and ML algorithm are shown to provide highly accurate models for both dipole moments and energies on water and 14 small molecules. To demonstrate the ability of MOB-ML to function as generalized density-matrix functionals for molecular dipole moments and energies of organic molecules, we further apply the proposed MOB-ML approach to train and test the molecules from the QM9 dataset. The application of local scalable GPR with Gaussian mixture model unsupervised clustering GPR scales up MOB-ML to a large-data regime while retaining the prediction accuracy. In addition, compared with the literature results, MOB-ML provides the best test mean absolute errors of 4.21 mD and 0.045 kcal/mol for dipole moment and energy models, respectively, when training on 110 000 QM9 molecules. The excellent transferability of the resulting QM9 models is also illustrated by the accurate predictions for four different series of peptides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A sparse Gaussian sigmoid basis function approximation of hyperspectral data for detection of solids

We define a new characterization of emissivity and reflectance curves for compositional exploitation of hyperspectral data. Here, our method decomposes each spectrum into a sparse set of Gaussian sigmoid components using penalized regression. Detection is based on the combination of Gaussian sigmoid components unique to a target material. Focusing on the presence of spectral upslopes and downslopes rather than spectral correlations makes detection more robust to both target variation and spectral variability from atmosphere and background encountered during the collection process. We present simulation studies that demonstrate the potential to reduce false positive rates without compromising sensitivity. Characterization of long–wave infrared (LWIR) experimental data validates our method using minerals of different particle sizes, measurement angles, and collection conditions.

97 MATHEMATICS AND COMPUTING↗

determining the mass of kepler-78b with nonparametric Gaussian process estimation

Kepler-78b is a transiting planet that is 1.2 times the radius of Earth and orbits a young, active K dwarf every 8 hr. The mass of Kepler-78b has been independently reported by two teams based on radial velocity (RV) measurements using the HIRES and HARPS-N spectrographs. Due to the active nature of the host star, a stellar activity model is required to distinguish and isolate the planetary signal in RV data. Whereas previous studies tested parametric stellar activity models, we modeled this system using nonparametric Gaussian process (GP) regression. We produced a GP regression of relevant Kepler photometry. We then use the posterior parameter distribution for our photometric fit as a prior for our simultaneous GP + Keplerian orbit models of the RV data sets. We tested three simple kernel functions for our GP regressions. Based on a Bayesian likelihood analysis, we selected a quasi-periodic kernel model with GP hyperparameters coupled between the two RV data sets, giving a Doppler amplitude of 1.86 ± 0.25 m s{sup −1} and supporting our belief that the correlated noise we are modeling is astrophysical. The corresponding mass of 1.87{sub −0.26}{sup +0.27} M{sub ⊕} is consistent with that measured in previous studies, and more robust due to our nonparametric signal estimation. Based on our mass and the radius measurement from transit photometry, Kepler-78b has a bulk density of 6.0{sub −1.4}{sup +1.9} g cm{sup −3}. We estimate that Kepler-78b is 32% ± 26% iron using a two-component rock-iron model. This is consistent with an Earth-like composition, with uncertainty spanning Moon-like to Mercury-like compositions.

79 ASTRONOMY AND ASTROPHYSICS↗

SULI Technical Report: WIRE (Draft)

Local climate undoubtedly affects facility maintenance requirements; however, knowledge of the link between weather and facility condition over time is grounded in the 1990s and needs to be updated. WIRE aims to quantify the relationship between climate and infrastructure by calculating climate-based service life estimates to be used in BUILDER, a program the Department of Energy (DOE) uses for facility management. Phase One of this project is a case study on window-or-wall-mounted AC units. The independent variables influencing service life are load, maintenance quality, and effective age, with load calculated from three weather variables: temperature, humidity, and solar radiation. Given sufficient data, service life can be calculated as a function of these variables via regression. WIRE is currently collecting data. Preliminary estimated service life projections as a function of load derived from Lawrence Livermore National Laboratory (LLNL) data support the need for updated analysis on the relationship between climate and facility condition. With the groundwork laid, WIRE can generate deliverables immediately upon arrival of data. The methods used for this case study constitute a framework for deriving the relationship between weather variables and service life for other components.

42 ENGINEERING↗

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗

Proper orthogonal descriptors for efficient and accurate interatomic potentials

Here, we present the proper orthogonal descriptors for efficient and accuracy representation of the potential energy surface. The potential energy surface is represented as a many-body expansion of parametrized potentials in which the potentials are functions of atom positions and parameters. The proper orthogonal decomposition is employed to decompose the parametrized potentials into a set of proper orthogonal descriptors (PODs). Because of the rapid convergence of the proper orthogonal decomposition, relevant snapshots can be sampled exhaustively to represent the atomic neighborhood environment accurately with a small number of descriptors. The proper orthogonal descriptors are used to develop interatomic potentials by using a linear expansion of the descriptors and determining the expansion coefficients from a weighted least-squares regression against a density functional theory (DFT) training set. We present a comprehensive evaluation of the POD potentials on previously published DFT data sets comprising Li, Mo, Cu, Ni, Si, Ge, and Ta elements. The data sets represent a diverse pool of metals, transition metals, and semiconductors. The accuracy of the POD potentials are comparable to that of state-of-the-art machine learning potentials such as the spectral neighbor analysis potential (SNAP) and the atomic cluster expansion (ACE).

97 MATHEMATICS AND COMPUTING↗

Model selection and signal extraction using Gaussian Process regression

We present a novel computational approach for extracting localized signals from smooth background distributions. We focus on datasets that can be naturally presented as binned integer counts, demonstrating our procedure on the CERN open dataset with the Higgs boson signature, from the ATLAS collaboration at the Large Hadron Collider. Our approach is based on Gaussian Process (GP) regression — a powerful and flexible machine learning technique which has allowed us to model the background without specifying its functional form explicitly and separately measure the background and signal contributions in a robust and reproducible manner. Unlike functional fits, our GP-regression-based approach does not need to be constantly updated as more data becomes available. We discuss how to select the GP kernel type, considering trade-offs between kernel complexity and its ability to capture the features of the background distribution. We show that our GP framework can be used to detect the Higgs boson resonance in the data with more statistical significance than a polynomial fit specifically tailored to the dataset. Finally, we use Markov Chain Monte Carlo (MCMC) sampling to confirm the statistical significance of the extracted Higgs signature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deep learning: A guide for practitioners in the physical sciences

Machine learning is finding increasingly broad applications in the physical sciences. This most often involves building a model relationship between a dependent, measurable output, and an associated set of controllable, but complicated, independent inputs. We present a tutorial on current techniques in machine learning—a jumping-off point for interested researchers to advance their work. We focus on deep neural networks with an emphasis on demystifying deep learning. We begin with background ideas in machine learning and some example applications from current research in plasma physics. We discuss supervised learning techniques for modeling complicated functions, beginning with familiar regression schemes, and then advancing to more sophisticated deep learning methods. We also address unsupervised learning and techniques for reducing the dimensionality of input spaces. Along the way, we describe methods for practitioners to help ensure that their models generalize from their training data to as-yet-unseen test data. We describe classes of tasks—predicting scalars, handling images, and fitting time-series—and prepare the reader to choose an appropriate technique. We finally point out some limitations to modern machine learning and speculate on some ways that practitioners from the physical sciences may be particularly suited to help.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cycle length dependence of stellar magnetic activity and solar cycle 23

Solar cycle (SC) 23 was extraordinarily long with remarkably low magnetic activity. We have investigated whether this is a common behavior of solar-type stars. From the Ca ii H and K line intensities of 111 stars observed at Mount Wilson Observatory from 1966 to 1991, we have retrieved data of all 23 G-type stars and recalculated their cycle lengths using the damped least-squares method for the chromospheric activity index S as a function of time. A regression analysis was performed to find relations between the derived cycle length, P{sub avg}, and the index for excess chromospheric emission, R{sub HK}{sup ′}. As a noteworthy result, we found a segregation between young and old solar-type stars in the cycle length-activity correlation. We incorporated the relation for the solar-type stars into the previously known rule for stellar chromospheric activity and brightness to estimate the variation of solar brightness from SC 22 to SC 23 as (0.12 ± 0.06)%, much higher than the actual variation of total solar irradiance (TSI) ≤0.02%. We have then examined solar spectral irradiance (SSI) to find a good phase correlation with a sunspot number in the wavelength range of 170–260 nm, which is close to the spectral range effective in heating the Earth’s atmosphere. Therefore, it appears that SSI rather than TSI is a good indicator of the chromospheric activity, and its cycle length dependent variation would be more relevant to the possible role of the Sun in the cyclic variation of the Earth’s atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Degradation Mechanisms and Lifetime Prediction for Lithium-Ion Batteries -- A Control Perspective: Preprint

Predictive models of Li-ion battery lifetime must consider a multiplicity of electrochemical, thermal, and mechanical degradation modes experienced by batteries in application environments. To complicate matters, Li-ion batteries can experience different degradation trajectories that depend on storage and cycling history of the application environment. Rates of degradation are controlled by factors such as temperature history, electrochemical operating window, and charge/discharge rate. We present a generalized battery life prognostic model framework for battery systems design and control. The model framework consists of trial functions that are statistically regressed to Li-ion cell life datasets wherein the cells have been aged under different levels of stress. Degradation mechanisms and rate laws dependent on temperature, storage, and cycling condition are regressed to the data, with multiple model hypotheses evaluated and the best model down-selected based on statistics. The resulting life prognostic model, implemented in state variable form, is extensible to arbitrary real-world scenarios. The model is applicable in real-time control algorithms to maximize battery life and performance. We discuss efforts to reduce lifetime prediction error and accommodate its inevitable impact in controller design.

25 ENERGY STORAGE↗

Linking leaf dark respiration to leaf traits and reflectance spectroscopy across diverse forest types

Leaf dark respiration (R dark ), an important yet rarely quantified component of carbon cycling in forest ecosystems, is often simulated from leaf traits such as the maximum carboxylation capacity (V cmax ), leaf mass per area (LMA), nitrogen (N) and phosphorus (P) concentrations, in terrestrial biosphere models. However, the validity of these relationships across forest types remains to be thoroughly assessed. Here, in this study, we analyzed R dark variability and its associations with V cmax and other leaf traits across three temperate, subtropical and tropical forests in China, evaluating the effectiveness of leaf spectroscopy as a superior monitoring alternative. We found that leaf magnesium and calcium concentrations were more significant in explaining cross-site R dark than commonly used traits like LMA, N and P concentrations, but univariate trait–R dark relationships were always weak (r 2 ≤ 0.15) and forest-specific. Although multivariate relationships of leaf traits improved the model performance, leaf spectroscopy outperformed trait–R dark relationships, accurately predicted cross-site R dark (r 2 = 0.65) and pinpointed the factors contributing to R dark variability. Our findings reveal a few novel traits with greater cross-site scalability regarding R dark , challenging the use of empirical trait–R dark relationships in process models and emphasize the potential of leaf spectroscopy as a promising alternative for estimating R dark , which could ultimately improve process modeling of terrestrial plant respiration.

59 BASIC BIOLOGICAL SCIENCES↗

A Metagenomic Meta-analysis Reveals Functional Signatures of Health and Disease in the Human Gut Microbiome

While recent research indicates that human health is affected by the gut microbiome, the functional mechanisms that underlie host-microbiome interactions remain poorly resolved. Metagenomic clinical studies can address this problem by revealing specific microbial functions that stratify healthy and diseased individuals. To improve our understanding of the relationship between the gut microbiome and health, we conducted the first integrative functional analysis of nearly 2,000 publicly available fecal metagenomic samples obtained from eight clinical studies. We identified characteristics of the gut microbiome that associate generally with disease, including functional alpha-diversity, beta-diversity, and beta-dispersion. Using regression modeling, we identified specific microbial functions that robustly stratify diseased individuals from healthy controls. Many of these functions overlapped multiple diseases, suggesting a general role in host health, while others were specific to a single disease and may indicate disease-specific etiologies. Our results clarify potential microbiome-mediated mechanisms of disease and reveal features of the microbiome that may be useful for the development of microbiome-based diagnostics.

59 BASIC BIOLOGICAL SCIENCES↗

Application of visible and near-infrared spectroscopy to classification of Miscanthus species

Here, the feasibility of visible and near infrared (NIR) spectroscopy as tool to classify Miscanthus samples was explored in this study. Three types of Miscanthus plants, namely, M. sinensis , M. sacchariflorus and M. fIoridulus , were analyzed using a NIR spectrophotometer. Several classification models based on the NIR spectra data were developed using line discriminated analysis (LDA), partial least squares (PLS), least squares support vector machine regression (LSSVR), radial basis function (RBF) and neural network (NN). The principal component analysis (PCA) presented rough classification with overlapping samples, while the models of Line_LSSVR, RBF_LSSVR and RBF_NN presented almost same calibration and validation results. Due to the higher speed of Line_LSSVR than RBF_LSSVR and RBF_NN, we selected the line_LSSVR model as a representative. In our study, the model based on line_LSSVR showed higher accuracy than LDA and PLS models. The total correct classification rates of 87.79 and 96.51% were observed based on LDA and PLS model in the testing set, respectively, while the line_LSSVR showed 99.42% of total correct classification rate. Meanwhile, the lin_LSSVR model in the testing set showed correct classification rate of 100, 100 and 96.77% for M. sinensis , M. sacchariflorus and M. fIoridulus , respectively. The lin_LSSVR model assigned 99.42% of samples to the right groups, except one M. fIoridulus sample. The results demonstrated that NIR spectra combined with a preliminary morphological classification could be an effective and reliable procedure for the classification of Miscanthus species.

59 BASIC BIOLOGICAL SCIENCES↗