Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 343 records · Page 19

The saturated activities of Na-22, Mn-54, and Co-56 and the depth of sampling of soils

A multivariate linear regression study shows that the saturated activities of Na-22, Co-56, and Mn-54 are functions of the chemical composition and the reciprocal of the cube root of the sample weight. A procedure for taking the temporal variation in these short-lived radionuclides into account is described. The regression results are used to estimate the average depths of sampling of soil material. The data indicate that all 28 samples were in the same conditions of exposure when recovered as they had been for the last approximately 10 years. The saturation activity of Mn-54 is a very weak function of time.

Keith, J. E.↗

Modeling the Echelle Spectra Continuum with Alpha Shapes and Local Regression Fitting

Continuum normalization of echelle spectra is an important data analysis step that is difficult to automate. Polynomial fitting requires a reasonably high-order model to follow the steep slope of the blaze function. However, in the presence of deep spectral lines, a high-order polynomial fit can result in ripples in the normalized continuum that increase errors in spectral analysis. Here, we present two algorithms for flattening the spectrum continuum. The Alpha-shape Fitting to Spectrum algorithm is completely data driven, using an alpha shape to obtain an initial estimate of the blaze function. The Alpha-shape and Lab Source Fitting to Spectrum algorithm incorporates a continuum constraint from a laboratory source reference spectrum for the blaze function estimation. These algorithms are tested on a simulated spectrum, where we demonstrate improved normalization compared to polynomial regression for continuum fitting. We show an additional application, using the algorithms for mitigation of spatially correlated quantum efficiency variations and fringing in the charge-coupled device detector of the EXtreme PREcision Spectrometer.

79 ASTRONOMY AND ASTROPHYSICS↗

Remote quantification of Cm(III) and HNO 3 by fluorescence spectroscopy and chemometrics

A unique approach to remotely quantify Cm(III) (0–100 µg mL −1 ) in HNO 3 (1–12 M) using steady-state laser fluorescence spectroscopy and multivariate regression models was developed. Photoluminescence is amenable to remote measurements using fiber-optic cables and is sensitive to numerous lanthanide and actinide species. In-line measurements can provide feedback to support complex processing in harsh environments (e.g., hot cells) to help guide and optimize radiochemical separations. In this work, Cm(III) spectra were acquired remotely in a glove box as a function of HNO 3 concentration to better understand spectral characteristics and evaluate the utility of multivariate regression models in this system. Furthermore, the Cm(III) fluorescence peak shape, width, position, and intensity changed significantly as a function of HNO 3 concentration, likely because of the displacement of emission quenching inner-sphere water molecules and complexation with nitrate ions. Despite significant covariance and nonlinearity in the data, a D-optimal design strategy successfully minimized training set sample size and was used to build effective partial least squares regression models for Cm(III) and HNO 3 concentrations without a priori knowledge of solution conditions. Chemometrics for modeling complex fluorescence spectra are promising and may find widespread applicability for online analysis in numerous chemical systems found in the nuclear field.

Actinide↗

Regression Model Term Selection for the Analysis of Strain-Gage Balance Calibration Data

The paper discusses the selection of regression model terms for the analysis of wind tunnel strain-gage balance calibration data. Different function class combinations are presented that may be used to analyze calibration data using either a non-iterative or an iterative method. The role of the intercept term in a regression model of calibration data is reviewed. In addition, useful algorithms and metrics originating from linear algebra and statistics are recommended that will help an analyst (i) to identify and avoid both linear and near-linear dependencies between regression model terms and (ii) to make sure that the selected regression model of the calibration data uses only statistically significant terms. Three different tests are suggested that may be used to objectively assess the predictive capability of the final regression model of the calibration data. These tests use both the original data points and regression model independent confirmation points. Finally, data from a simplified manual calibration of the Ames MK40 balance is used to illustrate the application of some of the metrics and tests to a realistic calibration data set.

Ulbrich, Norbert Manfred↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Characterization of the oncogenic function of centromere protein F in hepatocellular carcinoma

Highlights: •Overexpression of CENPF is frequently detected in HCC. •Upregulation of CENPF serves as an independent prognosis factor in HCC patients. •CENPF functions as an oncogene in HCC by promoting cell G2/M transition. -- Abstract: Centromere protein F (CENPF) is an essential nuclear protein associated with the centromere-kinetochore complex and plays a critical role in chromosome segregation during mitosis. Up-regulation of CENPF expression has previously been detected in several solid tumors. In this study, we aim to study the expression and functional role of CENPF in hepatocellular carcinoma (HCC). We found CENPF was frequently overexpressed in HCC as compared with non-tumor tissue. Up-regulated CENPF expression in HCC was positively correlated with serum AFP, venous invasion, advanced differentiation stage and a shorter overall survival. Cox regression analysis found that overexpression of CENPF was an independent prognosis factor in HCC. Functional studies found that silencing CENPF could decrease the ability of the cells to proliferate, form colonies and induce tumor formation in nude mice. Silencing CENPF also resulted in the cell cycle arrest at G2/M checkpoint by down-regulating cell cycle proteins cdc2 and cyclin B1. Our data suggest that CENPF is frequently overexpressed in HCC and plays a critical role in driving HCC tumorigenesis.

60 APPLIED LIFE SCIENCES↗

Inverse regression for ridge recovery: a data-driven approach for parameter reduction in computer experiments

Parameter reduction can enable otherwise infeasible design and uncertainty studies with modern computational science models that contain several input parameters. In statistical regression, techniques for sufficient dimension reduction (SDR) use data to reduce the predictor dimension of a regression problem. A computational scientist hoping to use SDR for parameter reduction encounters a problem: a computer prediction is best represented by a deterministic function of the inputs, so data comprised of computer simulation queries fail to satisfy the SDR assumptions. To address this problem, here we interpret SDR methods sliced inverse regression (SIR) and sliced average variance estimation (SAVE) as estimating the directions of a ridge function, which is a composition of a low-dimensional linear transformation with a nonlinear function. Within this interpretation, SIR and SAVE estimate matrices of integrals whose column spaces are contained in the ridge directions’ span; we analyze and numerically verify convergence of these column spaces as the number of computer model queries increases. Moreover, we show example functions that are not ridge functions but whose inverse conditional moment matrices are low-rank. Consequently, the computational scientist should beware when using SIR and SAVE for parameter reduction, since SIR and SAVE may mistakenly suggest that truly important directions are unimportant.

97 MATHEMATICS AND COMPUTING↗

Multinomial logistic regression modelling of obesity and overweight among primary school students in a rural area of Negeri Sembilan

Multinomial logistic regression is widely used to model the outcomes of a polytomous response variable, a categorical dependent variable with more than two categories. The model assumes that the conditional mean of the dependent categorical variables is the logistic function of an affine combination of predictor variables. Its procedure gives a number of logistic regression models that make specific comparisons of the response categories. When there are q categories of the response variable, the model consists of q-1 logit equations which are fitted simultaneously. The model is validated by variable selection procedures, tests of regression coefficients, a significant test of the overall model, goodness-of-fit measures, and validation of predicted probabilities using odds ratio. This study used the multinomial logistic regression model to investigate obesity and overweight among primary school students in a rural area on the basis of their demographic profiles, lifestyles and on the diet and food intake. The results indicated that obesity and overweight of students are related to gender, religion, sleep duration, time spent on electronic games, breakfast intake in a week, with whom meals are taken, protein intake, and also, the interaction between breakfast intake in a week with sleep duration, and the interaction between gender and protein intake.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How analytic choices can affect the extraction of electromagnetic form factors from elastic electron scattering cross section data

Scientists often try to incorporate prior knowledge into their regression algorithms, such as a particular analytic behavior or a known value at a kinematic endpoint. Unfortunately, there is often no unique way to make use of this prior knowledge, and thus, different analytic choices can lead to very different regression results from the same set of data. In this paper, to illustrate this point in the context of the proton electromagnetic form factors, we use the Mainz elastic data with its 1422 cross section points and 31 normalization parameters. Starting with a complex unbound non-linear regression, we will show how the addition of a single theory-motivated constraint removes an oscillation from the magnetic form factor and shifts the extracted proton charge radius. We then repeat both regressions using the same algorithm, but with a rebinned version of the Mainz dataset. These examples illustrate how analytic choices, such as the function that is being used or even the binning of the data, can dramatically affect the results of a complex regression. These results also demonstrate why it is critical when using regression algorithms to have either a physical model in mind or a firm mathematical basis

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian and LASSO Regressions for Comparative Permeability Modeling of Sandstone Reservoirs

In this research, Bayesian model averaging (BMA) and least absolute shrinkage and selection operator regression (LASSO) algorithms were adopted for the core permeability modeling as a function of well log and core measurements. More specifically, the core permeability (dependent factor) was modeled given the well log and core data (independent variables) and then predicted in non-cored intervals of a well in a sandstone formation. The BMA is a stochastic linear modeling and a Bayesian parameter selection. Among 50 linear models generated as a function of the independent variables and depicted in Occam’s window, the best model of the highest posterior probability is determined. The best model has the optimal subset parameters that most influences the response factor. The core permeability modeling was again conducted using the LASSO algorithm, which adopts a completely different way of modeling and subset selection by using the penalized least-squared equation. Both BMA and LASSO resulted in a very accurate prediction of the core permeability by achieving perfect measured and calculated permeability matching. Matches for the two algorithms were quantified by computing the adjusted R-squared and root-mean-squared prediction error. In addition, results of the BMA and LASSO were compared to the conventional multiple linear regression (LM). The LASSO algorithm led to an accurate matching similar to LM, but slightly better than BMA. Therefore, both BMA and LASSO represent integrated procedures for accurate permeability estimation similar to conventional regression analysis.

58 GEOSCIENCES↗

Predictive Models for Regional Hepatic Function Based on 99mTc-IDA SPECT and Local Radiation Dose for Physiologic Adaptive Radiation Therapy

Purpose: High-dose radiation therapy (RT) for intrahepatic cancer is limited by the development of liver injury. This study investigated whether regional hepatic function assessed before and during the course of RT using 99mTc-labeled iminodiacetic acid (IDA) single photon emission computed tomography (SPECT) could predict regional liver function reserve after RT. Methods and Materials: Fourteen patients treated with RT for intrahepatic cancers underwent dynamic 99mTc-IDA SPECT scans before RT, during, and 1 month after completion of RT. Indocyanine green (ICG) tests, a measure of overall liver function, were performed within 1 day of each scan. Three-dimensional volumetric hepatic extraction fraction (HEF) images of the liver were estimated by deconvolution analysis. After coregistration of the CT/SPECT and the treatment planning CT, HEF dose–response functions during and after RT were generated. The volumetric mean of the HEFs in the whole liver was correlated with ICG clearance time. Three models, dose, priori, and adaptive models, were developed using multivariate linear regression to assess whether the regional HEFs measured before and during RT helped predict regional hepatic function after RT. Results: The mean of the volumetric liver HEFs was significantly correlated with ICG clearance half-life time (r=−0.80, P<.0001), for all time points. Linear correlations between local doses and regional HEFs 1 month after RT were significant in 12 patients. In the priori model, regional HEF after RT was predicted by the planned dose and regional HEF assessed before RT (R=0.71, P<.0001). In the adaptive model, regional HEF after RT was predicted by regional HEF reassessed during RT and the remaining planned local dose (R=0.83, P<.0001). Conclusions: 99mTc-IDA SPECT obtained during RT could be used to assess regional hepatic function and helped predict post-RT regional liver function reserve. This could support individualized adaptive radiation treatment strategies to maximize tumor control and minimize the risk of liver damage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Physiological Factors Contributing to Postflight Changes in Functional Performance

Astronauts experience alterations in multiple physiological systems due to exposure to the microgravity conditions of space flight. These physiological changes include sensorimotor disturbances, cardiovascular deconditioning and loss of muscle mass and strength. These changes might affect the ability of crewmembers to perform critical mission tasks immediately after landing on lunar and Martian surfaces. To date, changes in functional performance have not been systematically studied or correlated with physiological changes. To understand how changes in physiological function impact functional performance an interdisciplinary pre/postflight testing regimen (Functional Task Test, FTT) has been developed that systematically evaluates both astronaut postflight functional performance and related physiological changes. The overall objectives of the FTT are to: Develop a set of functional tasks that represent critical mission tasks for Constellation. Determine the ability to perform these tasks after flight. Identify the key physiological factors that contribute to functional decrements. Use this information to develop targeted countermeasures. The functional test battery was designed to address high priority tasks identified by the Constellation program as critical for mission success. The set of functional tests making up the FTT include the: 1) Seat Egress and Walk Test, 2) Ladder Climb Test, 3) Recovery from Fall/Stand Test, 4) Rock Translation Test, 5) Jump Down Test, 6) Torque Generation Test, and 7) Construction Activity Board Test. Corresponding physiological measures include assessments of postural and gait control, dynamic visual acuity, fine motor control, plasma volume, orthostatic intolerance, upper and lower body muscle strength, power, fatigue, control and neuromuscular drive. Crewmembers will perform both functional and physiological tests before and after short (Shuttle) and long-duration (ISS) space flight. Data will be collected on R+0 (Shuttle only), R+1, R+6 and R+30. Using a multivariate regression model we will identify which physiological systems contribute the most to impaired performance on each functional test. This will allow us to identify the physiological systems that play the largest role in decrement in functional performance. Using this information we can then design and implement countermeasures that specifically target the physiological systems most responsible for the altered functional performance associated with space flight.

Bloomberg, J. J.↗

Calibration and Data Analysis of the MC-130 Air Balance

Design, calibration, calibration analysis, and intended use of the MC-130 air balance are discussed. The MC-130 balance is an 8.0 inch diameter force balance that has two separate internal air flow systems and one external bellows system. The manual calibration of the balance consisted of a total of 1854 data points with both unpressurized and pressurized air flowing through the balance. A subset of 1160 data points was chosen for the calibration data analysis. The regression analysis of the subset was performed using two fundamentally different analysis approaches. First, the data analysis was performed using a recently developed extension of the Iterative Method. This approach fits gage outputs as a function of both applied balance loads and bellows pressures while still allowing the application of the iteration scheme that is used with the Iterative Method. Then, for comparison, the axial force was also analyzed using the Non-Iterative Method. This alternate approach directly fits loads as a function of measured gage outputs and bellows pressures and does not require a load iteration. The regression models used by both the extended Iterative and Non-Iterative Method were constructed such that they met a set of widely accepted statistical quality requirements. These requirements lead to reliable regression models and prevent overfitting of data because they ensure that no hidden near-linear dependencies between regression model terms exist and that only statistically significant terms are included. Finally, a comparison of the axial force residuals was performed. Overall, axial force estimates obtained from both methods show excellent agreement as the differences of the standard deviation of the axial force residuals are on the order of 0.001 % of the axial force capacity.

Booth, Dennis↗

Low-level arsenic exposure via drinking water consumption and female fecundity - A preliminary investigation

High level arsenic exposure is associated with reproductive toxicity in experimental and observational studies; however, few data exist to assess risks at low levels. Even less data are available to evaluate the impact of low level arsenic exposure on human fecundity. Our aim in this pilot study was a preliminary evaluation of associations between low level drinking water arsenic contamination and female fecundity. This retrospective study was conducted among women previously recruited to a hospital-based case-control study of spontaneous pregnancy loss in Timiá¹£ County, Romania. Women (n=94) with planned pregnancies of 5–20 weeks gestation completed a comprehensive physician-administered study questionnaire and reported the number of menstrual cycles attempting to conceive as the time to pregnancy (TTP). Drinking water samples were collected from residential drinking water sources and we determined arsenic levels using hydride generation-atomic absorption spectrometry (HG-AAS). Multivariable Cox-proportional hazards regression with Efron approximation was employed to evaluate TTP as a function of drinking water arsenic concentrations among planned pregnancies, adjusted for covariates. There was no main effect for drinking water arsenic exposure, yet the conditional probability for pregnancy was modestly lower among arsenic exposed women with longer TTPs, relative to women with shorter TTPs, and relative to unexposed women. For example, 1 µg/L average drinking water arsenic conferred 5%, 8%, and 10% lower likelihoods for pregnancy in the 6th, 9th, and 12th cycles, respectively (P=0.01). While preliminary, our results suggest that low level arsenic contamination in residential drinking water sources may further impair fecundity among women with longer waiting times; however, this hypothesis requires confirmation by a future, more definitive study. - Highlights: • We assessed low level drinking water arsenic as a predictor of fecundability. • Arsenic did not affect time to pregnancy among women conceiving quickly. • Arsenic increased time to pregnancy among women taking longer to conceive. • Low level drinking water arsenic may adversely impact women with lower fecundity.

54 ENVIRONMENTAL SCIENCES↗

Insights into Prismatic Loop Formation in Irradiated Fe–Cr Alloys from Hypothesis-Driven Active Learning and Causal Analysis

Neutron and electron irradiation experimental studies conducted on body-centered cubic Fe and Fe–Cr alloys have established two prismatic dislocation loop populations, which have Burgers vectors of either a/2$\langle$111$\rangle$ or a$\langle$100$\rangle$. Here, the loop formation depends on factors such as dose (D), dose rate (D rt ), temperature (T), chromium content (Cr%), and other alloying elements. Hence, it is important to understand how irradiation-induced dislocation loops evolve conditional upon the loop characteristics, such as loop density (DD), average loop size d̅, and irradiation parameters (D, D rt , T, and irradiation type), which is still an active area of research. To understand these complex structure–property relationships, machine learning (ML) is employed in a three-step approach. This includes imputing missing data with a k-nearest neighbor, generating functionalized features, and assessing feature importance with random forest classification and regression. Physics-based features are incorporated in a hypothesis-driven active learning scheme to overcome data unavailability challenges. Insights obtained from ML models (i) to categorize dislocation loop types, show the highest correlation with d̅; (ii) Log(DD), obtained through mathematical formulations involving D, Cr%, d̅, and T (e.g., Log(DD) ~ D + exp(-Cr%) + 1/d̅ and log(DD) ~ D + exp(-Cr%) + 1/T). Hypothesis-driven active learning is able to predict Log(DD) in which the experimental date is not known. Causal models verify cause–effect relationships for dislocation loop classification and irradiation factors in FeCr alloys.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Earthquake catalog-based machine learning identification of laboratory fault states and the effects of magnitude of completeness

Here, machine learning regression can predict macroscopic fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. Here we show that a similar approach is successful using event catalogs derived from the continuous data. Our methods are applicable to catalogs of arbitrary scale and magnitude of completeness. We investigate how machine learning regression from an event catalog of laboratory earthquakes performs as a function of the catalog magnitude of completeness. We find that strong model performance requires a sufficiently low magnitude of completeness, and below this magnitude of completeness, model performance saturates.

58 GEOSCIENCES↗

Kinetic modelling for zinc (II) ions biosorption onto Luffa cylindrica

The biosorption of Zinc (II) ions onto a biomaterial - Luffa cylindrica has been studied. This biomaterial was characterized by elemental analysis, surface area, pore size distribution, scanning electron microscopy, and the biomaterial before and after sorption, was characterized by Fourier Transform Infra Red (FTIR) spectrometer. The kinetic nonlinear models fitted were Pseudo-first order, Pseudo-second order and Intra-particle diffusion. A comparison of non-linear regression method in selecting the kinetic model was made. Four error functions, namely coefficient of determination (R{sup 2}), hybrid fractional error function (HYBRID), average relative error (ARE), and sum of the errors squared (ERRSQ), were used to predict the parameters of the kinetic models. The strength of this study is that a biomaterial with wide distribution particularly in the tropical world and which occurs as waste material could be put into effective utilization as a biosorbent to address a crucial environmental problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗