Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local 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 55 records · Page 3

A Swing of the Pendulum: The Chemodynamics of the Local Stellar Halo Indicate Contributions from Several Radial Merger Events

We find that the chemical abundances and dynamics of APOGEE and GALAH stars in the local stellar halo are inconsistent with a scenario in which the inner halo is primarily composed of debris from a single massive, ancient merger event, as has been proposed to explain the Gaia-Enceladus/Gaia Sausage (GSE) structure. The data contain trends of chemical composition with energy that are opposite to expectations for a single massive, ancient merger event, and multiple chemical evolution paths with distinct dynamics are present. We use a Bayesian Gaussian mixture model regression algorithm to characterize the local stellar halo, and find that the data are fit best by a model with four components. We interpret these components as the Virgo Radial Merger (VRM), Cronus, Nereus, and Thamnos; however, Nereus and Thamnos likely represent more than one accretion event because the chemical abundance distributions of their member stars contain many peaks. Although the Cronus and Thamnos components have different dynamics, their chemical abundances suggest they may be related. We show that the distinct low- and high-α halo populations from Nissen & Schuster are explained by VRM and Cronus stars, as well as some in situ stars. Because the local stellar halo contains multiple substructures, different popular methods of selecting GSE stars will actually select different mixtures of these substructures, which may change the apparent chemodynamic properties of the selected stars. We also find that the Splash stars in the Solar region are shifted to higher v $\phi$ and slightly lower [Fe/H] than previously reported.

79 ASTRONOMY AND ASTROPHYSICS↗

Determining Latitudinal Extent of Energetic Electron Precipitation Using MEPED On-Board NOAA/POES

Energetic Electron Precipitation (EEP) from the plasma sheet and the radiation belts ionizes the polar lower thermosphere and mesosphere. EEP increases the production of NO x and HO x , which will catalytically destroy ozone, an important element of atmospheric dynamics. Therefore, measurement of the latitudinal extent of the precipitation boundaries is important in quantifying the atmospheric effects of the Sun-Earth interaction. This study uses measurements by the Medium Energy Proton Electron Detector (MEPED) of six NOAA/POES and EUMETSAT/METOP satellites from 2004 to 2014 to determine the latitudinal boundaries of EEP and their variability with geomagnetic activity and solar wind drivers. Variation of the boundaries for different electron energies and Magnetic Local Time (MLT) is studied. Regression analyses are applied to determine the best predictor variable based on solar wind parameters and geomagnetic indices. The highest correlation was found for the pressure-corrected Dst index when applying a linear regression model. A model of the equatorward EEP boundary is developed separately for three different energy channels, >43, >114, and >292 keV, and for 3 hour MLT sectors. For >43 keV EEP, 80% of the equatorward boundaries predicted by the model are within ±2.2° cgmlat. The model exhibits a solar cycle bias where it systematically exaggerates the equatorward movement of the EEP region during solar minimum. The highest accuracy of the model is found in periods dominated by corotating interaction regions/high speed solar wind streams. The result will be a key element for constructing a model of EEP variability to be applied in atmosphere climate models.

E. M. Babu↗

Estimation of Crop Gross Primary Production (GPP): Do Scaled (MODIS) Vegetation Indices Improve Performance? - 2

Satellite remote sensing estimates of Gross Primary Production (GPP) have routinely been made using spectral Vegetation Indices (VIs) over the past two decades. The Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), the green band Wide Dynamic Range Vegetation Index (WDRVIgreen), and the green band Chlorophyll Index (CIgreen) have been employed to estimate GPP under the assumption that GPP is proportional to the product of VI and photosynthetically active radiation (PAR) (where VI is one of four VIs: NDVI, EVI, WDRVIgreen, or CIgreen). However, the empirical regressions between VI*PAR and GPP measured locally at flux towers do not pass through the origin (i.e., the zero X-Y value for regressions). Therefore they are somewhat difficult to interpret and apply. This study investigates (1) what are the scaling factors and offsets (i.e., regression slopes and intercepts) between the fraction of PAR absorbed by chlorophyll of a canopy (fAPARchl) and the VIs, and (2) whether the scaled VIs developed in (1) can eliminate the deficiency and improve the accuracy of GPP estimates. Three AmeriFlux maize and soybean fields were selected for this study, two of which are irrigated and one is rainfed. The four VIs and fAPARchl of the fields were computed with the MODerate resolution Imaging Spectroradiometer (MODIS) satellite images. The GPP estimation performance for the scaled VIs was compared to results obtained with the original VIs and evaluated with standard statistics: the coefficient of determination (R2), the root mean square error (RMSE), and the coefficient of variation (CV). Overall, the scaled EVI obtained the best performance. The performance of the scaled NDVI, EVI and WDRVIgreen was improved across sites, crop types and soil/background wetness conditions. The scaled CIgreen did not improve results, compared to the original CIgreen. The scaled green band indices (WDRVIgreen, CIgreen) did not exhibit superior performance to either the scaled EVI or NDVI in estimating crop daily GPP at these agricultural fields. The scaled VIs are more physiologically meaningful than original un-scaled VIs, but scaling factors and offsets may vary across crop types and surface conditions.

Vegetation Index↗

DEPRECATED AI-Batt-OS (Autonomous Identification of Battery Life Models - Open Source) [SWR 21-17]

DEPRECATED. This repository was archived by the owner on Jun 30, 2026. It is now read-only. Open source implementation of some of the methods utilized by AI-Batt, a battery lifetime modeling and analysis toolkit provided by the National Laboratory of the Rockies (NLR). This software demonstrates the use of bi-level optimization and symbolic regression techniques to semi-autonomously identify algebraic models predicting the capacity fade of lithium-ion batteries during calendar aging. Modeling the degradation of batteries is a complex task, due to the difficulty in separating the time-dependent and time-independent factors impacting cell level degradation, across multiple data series with different numbers of measurements and/or data quality. Bi-level optimization enables model parameters to be optimized to either the entire data set or to individual data series, allowing statistical disambiguation of global behaviors (data series independent) and local behaviors (data series dependent). Symbolic regression is used to automatically search for optimal low-dimesional models predicting the variation of locally optimized parameters versus time-independent experimental variables from millions of possible models, resulting in a more accurate and repeatable model identification process than is possible by a manual search. The provided tools also implement cross-validation and bootstrap resampling schemes, empowering statistical model comparison/selection and quantification of model uncertainties. An example script replicates the results from the manuscript "Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning", submitted to ECS. All code is written in MATLAB. Requires the Statistics and Machine Learning Toolbox. Contact Dr. Paul Gasper at Paul.Gasper@nlr.gov for any questions.

Gasper, Paul↗

Vehicle Localization in 3D World Coordinates Using Single Camera at Traffic Intersection

Optimizing traffic control systems at traffic intersections can reduce the network-wide fuel consumption, as well as emissions of conventional fuel-powered vehicles. While traffic signals have been controlled based on predetermined schedules, various adaptive signal control systems have recently been developed using advanced sensors such as cameras, radars, and LiDARs. Among these sensors, cameras can provide a cost-effective way to determine the number, location, type, and speed of the vehicles for better-informed decision-making at traffic intersections. In this research, a new approach for accurately determining vehicle locations near traffic intersections using a single camera is presented. For that purpose, a well-known object detection algorithm called YOLO is used to determine vehicle locations in video images captured by a traffic camera. YOLO draws a bounding box around each detected vehicle, and the vehicle location in the image coordinates is converted to the world coordinates using camera calibration data. During this process, a significant error between the center of a vehicle’s bounding box and the real center of the vehicle in the world coordinates is generated due to the angled view of the vehicles by a camera installed on a traffic light pole. As a means of mitigating this vehicle localization error, two different types of regression models are trained and applied to the centers of the bounding boxes of the camera-detected vehicles. The accuracy of the proposed approach is validated using both static camera images and live-streamed traffic video. Based on the improved vehicle localization, it is expected that more accurate traffic signal control can be made to improve the overall network-wide energy efficiency and traffic flow at traffic intersections.

47 OTHER INSTRUMENTATION↗

Effects of local compositional and structural disorder on vacancy formation in entropy-stabilized oxides from first-principles

Entropic stabilization has evolved into a strategy to create new oxide materials and realize novel functional properties engineered through the alloy composition. Achieving an atomistic understanding of these properties to enable their design, however, has been challenging due to the local compositional and structural disorder that underlies their fundamental structure-property relationships. Here, we combine high-throughput atomistic calculations and linear regression algorithms to investigate the role of local configurational and structural disorder on the thermodynamics of vacancy formation in (MgCoNiCuZn)O-based entropy-stabilized oxides (ESOs) and their influence on the electrical properties. We find that the cation-vacancy formation energies decrease with increasing local tensile strain caused by the deviation of the bond lengths in ESOs from the equilibrium bond length in the binary oxides. The oxygen-vacancy formation strongly depends on structural distortions associated with the local configuration of chemical species. Vacancies in ESOs exhibit deep thermodynamic transition levels that inhibit electrical conduction. By applying the charge-neutrality condition, we determine that the equilibrium concentrations of both oxygen and cation vacancies increase with increasing Cu mole fraction. Our results demonstrate that tuning the local chemistry and associated structural distortions by varying alloy composition acts an engineering principle that enables controlled defect formation in multi-component alloys.

36 MATERIALS SCIENCE↗

A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data

The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data.We developed and cross validated models to predict daily PM(sub 2.5) at a 1X 1 km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1 X 1 km grid predictions. We used mixed models regressing PM(sub 2.5) measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

Air pollution↗

Examining the Link Between Food Prices and Food Insecurity: a Multi-Level Analysis of Maize Price and Birthweight in Kenya

In developing countries where many poor people rely on rainfed, locally produced food for the majority of their caloric intake, shifts in climate and weather patterns can dramatically reduce agricultural productivity. The reduction in agricultural productivity reduces overall food availability and ultimately impacts food accessibility, putting millions of people at risk for malnutrition. In this project we focus on Kenya where roughly a third of households are food insecure. We examine the relationship of the price of maize and low birth weight to help quantify the impact of local food prices on one outcome of household food insecurity. Using spatially referenced data from recent Kenyan Demographic and Health Survey datasets, price data, livelihood information, and a remotely sensed-based measure of local growing season productivity, we develop a dataset linking pregnancies occurring from 2001 to 2008 to the spatially and temporally relevant maize price data. We construct several regression models to examine the impact of local maize prices and remotely sensed based estimates of crop production on infant birth weight - specifically low birth weight. The results of the models highlight the importance of including community crop production to evaluate maize price impacts on low birth weight outcomes. Also, because of the positive correlation between pre-pregnancy maize prices and birth weight, the results suggest that some households may benefit from high prices or that high prices may impact the number of conceptions. More generally, our work demonstrates that multilevel models that account for community-level variation are important for disentangling these complex relationships and can contribute to the discussion of how to design more effective food policies.

Low birth weight↗

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↗

Atomic resolution convergent beam electron diffraction analysis using convolutional neural networks

Two types of convolutional neural network (CNN) models, a discrete classification network and a continuous regression network, were trained to determine local sample thickness from convergent beam diffraction (CBED) patterns of SrTiO 3 collected in a scanning transmission electron microscope (STEM) at atomic column resolution. Acquisition of atomic resolution CBED patterns for this purpose requires careful balancing of CBED feature size in pixels, acquisition speed, and detector dynamic range. The training datasets were derived from multislice simulations, which must be convolved with incoherent source broadening. Sample thicknesses were also determined using quantitative high-angle annular dark-field (HAADF) STEM images acquired simultaneously. The regression CNN performed well on sample thinner than 35 nm, with 70% of the CNN results within 1 nm of HAADF thickness, and 1.0 nm overall root mean square error between the two measurements. The classification CNN was trained for a thicknesses up to 100 nm and yielded 66% of CNN results within one classification increment of 2 nm of HAADF thickness. Our approach depends on methods from computer vision including transfer learning and image augmentation.

36 MATERIALS SCIENCE↗

Virtual sensing of wind turbine hub loads and drivetrain fatigue damage, Virtuelle Sensoren für die Messung von Hauptwellenlasten und Ermüdungsschäden im Antriebstrang von Windenergieanlagen

Abstract: This paper presents a Digital Twin for virtual sensing of wind turbine aerodynamic hub loads, as well as monitoring the accumulated fatigue damage and remaining useful life in drivetrain bearings based on measurements of the Supervisory Control and Data Acquisition (SCADA) and the drivetrain condition monitoring system (CMS). The aerodynamic load estimation is realized with data-driven regression models, while the estimation of local bearing loads and damage is conducted with physics-based, analytical models. Field measurements of the DOE 1.5 research turbine are used for model training and validation. The results show low errors of 6.4% and 1.1% in the predicted damage at the main and the generator side high-speed bearing respectively. Zusammenfassung: In diesem Aufsatz wird ein digitaler Zwilling für Windenergieanlagen vorgestellt, welcher die virtuelle Erfassung der Hauptwellenlasten und die Zustandsüberwachung von Ermüdungschäden und der verbleibende Nutzungsdauer der Antriebsstranglager ermöglicht. Der digital Zwilling nutzt Messdaten des Supervisory Control and Data Acquisition (SCADA) Systems und des Zustandsüberwachungssystems des Antriebsstranges (CMS). Die Berechnung der Hauptwellenlasten ist mit datenbasierten Regressionsmodellen umgesetzt, während die Berechnung der Lagerkräfte und der Ermüdungsschaden mit physikbasierten Modelle durchgeführt wird. Für die Modellentwicklung und -validierung werden Feldmessdaten der DOE 1.5MW Turbine eingesetzt. Die Abweichungen in den Ermüdungsschäden am Hauptwellenlager und am Generatorwellenlager betragen lediglich 6,4% beziehungsweise 1,1%.

17 WIND ENERGY↗

Soil Water Percolation Flux, April 2017 to March 2019, BR-Ma2, Manaus

Soil water percolation flux measurements derived from six passive wick flux meters across three topographic positions (valley, slope and plateau). Data has been processed to compute the mean for each topographic position where two flux meters were installed at each topographic position. Missing data represented a few percent of the total data set and were filled using two approaches. In the event one of the flux meters malfunctioned, we used data from the other sensor from the same topographic position without averaging. When both sensors were not working, such as during the period from September to October 2017 for the slope and valley and July 2018 for the plateau, data was filled using a nonlinear regression model that was based on local precipitation data and coefficients that were empirically derived for each topographic position. The metadata tab included in the .csv has additional information on locations where the sensors were installed, and other installation/maintenance details. Contact ksolander@lanl.gov if you need to use this dataset for additional information.

54 ENVIRONMENTAL SCIENCES↗

Temperature Increase Reduces Global Yields of Major Crops in Four Independent Estimates

Wheat, rice, maize, and soybean provide two-thirds of human caloric intake. Assessing the impact of global temperature increase on production of these crops is therefore critical to maintaining global food supply, but different studies have yielded different results. Here, we investigated the impacts of temperature on yields of the four crops by compiling extensive published results from four analytical methods: global grid-based and local point-based models, statistical regressions, and field-warming experiments. Results from the different methods consistently showed negative temperature impacts on crop yield at the global scale, generally underpinned by similar impacts at country and site scales. Without CO2 fertilization, effective adaptation, and genetic improvement, each degree-Celsius increase in global mean temperature would, on average, reduce global yields of wheat by 6.0%, rice by 3.2%, maize by 7.4%, and soybean by 3.1%. Results are highly heterogeneous across crops and geographical areas, with some positive impact estimates. Multi-method analyses improved the confidence in assessments of future climate impacts on global major crops and suggest crop- and region-specific adaptation strategies to ensure food security for an increasing world population.

Zhao, Chuang↗

The metallicity’s fundamental dependence on both local and global galactic quantities

ABSTRACT We study the scaling relations between gas-phase metallicity, stellar mass surface density (Σ*), star formation rate surface density (ΣSFR), and molecular gas surface density ($\Sigma _{{\rm H}_2}$) in local star-forming galaxies on scales of a kpc. We employ optical integral field spectroscopy from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, and ALMA data for a subset of MaNGA galaxies. We use partial correlation coefficients and Random Forest regression to determine the relative importance of local and global galactic properties in setting the gas-phase metallicity. We find that the local metallicity depends primarily on Σ* (the resolved mass–metallicity relation, rMZR), and has a secondary anticorrelation with ΣSFR (i.e. a spatially resolved version of the ‘Fundamental Metallicity Relation’, rFMR). We find that $\Sigma _{{\rm H}_2}$ is less important than ΣSFR in determining the local metallicity. This result indicates that gas accretion, resulting in local metallicity dilution and local boosting of star formation, is unlikely to be the primary origin of the rFMR. The local metallicity depends also on the global properties of galaxies. We find a strong dependence on the total stellar mass (M*) and a weaker (inverse) dependence on the total SFR. The global metallicity scaling relations, therefore, do not simply stem out of their resolved counterparts; global properties and processes, such as the global gravitational potential well, galaxy-scale winds and global redistribution/mixing of metals, likely contribute to the local metallicity, in addition to local production and retention.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-Driven Learning of Nonautonomous Systems

In this work, we present a numerical framework for recovering unknown nonautonomous dynamical systems with time-dependent inputs. To circumvent the difficulty presented by the nonautonomous nature of the system, our method transforms the solution state into piecewise integration of the system over a discrete set of time instances. The time-dependent inputs are then locally parameterized by using a proper model, for example, polynomial regression, in the pieces determined by the time instances. This transforms the original system into a piecewise parametric system that is locally time invariant. We then design a deep neural network structure to learn the local models. Once the network model is constructed, it can be iteratively used over time to conduct global system prediction. We provide theoretical analysis of our algorithm and present a number of numerical examples to demonstrate the effectiveness of the method.

97 MATHEMATICS AND COMPUTING↗

Data-Driven Learning of Non-Autonomous Systems

We present a numerical framework for recovering unknown non-autonomous dynamical systems with time-dependent inputs. To circumvent the difficulty presented by the non-autonomous nature of the system, our method transforms the solution state into piecewise integration of the system over a discrete set of time instances. The time-dependent inputs are then locally parameterized by using a proper model, for example, polynomial regression, in the pieces determined by the time instances. This transforms the original system into a piecewise parametric system that is locally time invariant. We then design a deep neural network structure to learn the local models. Once the network model is constructed, it can be iteratively used over time to conduct global system prediction. We provide theoretical analysis of our algorithm and present a number of numerical examples to demonstrate the effectiveness of the method.

97 MATHEMATICS AND COMPUTING↗

SRBench++: Principled Benchmarking of Symbolic Regression With Domain-Expert Interpretation

Symbolic regression searches for analytic expressions that accurately describe studied phenomena. The main promise of this approach is that it may return an interpretable model that can be insightful to users, while maintaining high accuracy. The current standard for benchmarking these algorithms is SRBench, which evaluates methods on hundreds of datasets that are a mix of real-world and simulated processes spanning multiple domains. At present, the ability of SRBench to evaluate interpretability is limited to measuring the size of expressions on real-world data, and the exactness of model forms on synthetic data. In practice, model size is only one of many factors used by subject experts to determine how interpretable a model truly is. Furthermore, SRBench does not characterize algorithm performance on specific, challenging sub-tasks of regression such as feature selection and evasion of local minima. In this work, we propose and evaluate an approach to benchmarking SR algorithms that addresses these limitations of SRBench by 1) incorporating expert evaluations of interpretability on a domain-specific task, and 2) evaluating algorithms over distinct properties of data science tasks. We evaluate 12 modern symbolic regression algorithms on these benchmarks and present an in-depth analysis of the results, discuss current challenges of symbolic regression algorithms and highlight possible improvements for the benchmark itself.

97 MATHEMATICS AND COMPUTING↗