Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariates”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying the drivers of ecosystem fluxes and water potential across the soil-plant-atmosphere continuum in an arid woodland

Dryland ecosystems occupy a vast swath of the terrestrial land surface and exert a sizeable impact on the cycling of carbon and water globally. These biomes are characterized by tightly coupled carbon and water cycles that respond rapidly to transient pulses in water availability. However, there exist many mechanistic uncertainties regarding the environmental drivers of, and linkages between, plant and ecosystem processes. Thus, drylands are often poorly represented in many vegetation and land surface models. An enhanced understanding of dryland ecosystem function is limited by the lack of long-term, co-located, and frequent measurements of plant and ecosystem processes. At a piñon-juniper woodland in southeastern Utah, USA, we collected a continuous dataset of meteorological conditions, soil water potential from surface to bedrock, tree water potential, and ecosystem carbon and water fluxes from eddy covariance. We found that predawn and midday tree water potential and daily ecosystem fluxes were highly sensitive to fluctuations in soil water availability, particularly in shallower layers, and that daytime variability in atmospheric drivers only loosely controlled these processes. The strong connections between shallow soil water potential, tree water potential, and ecosystem fluxes occurred because of the dominant role of precipitation pulses in driving vegetation activity, as even small pulses of moisture stimulated shallow soil water potential, tree water potential, and evapotranspiration for between 1 and 2 weeks. Carbon fluxes (net ecosystem exchange and gross primary productivity) were sensitive to precipitation pulses for longer, up to 3 weeks. Our results highlight that improved monitoring and sensing of shallow soil moisture can greatly enhance our understanding of dryland ecosystem function. Here, a better mechanistic understanding of the impacts of precipitation pulses is also needed to improve vegetation modeling of dryland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Remote sensing of seasonal variation of LAI and fAPAR in a deciduous broadleaf forest

Climate change is affecting the phenology of terrestrial ecosystems. In deciduous forests, phenology in leaf area index (LAI) is the primary driver of seasonal variation in the fraction of absorbed photosynthetically active radiation (fAPAR), which drives photosynthesis. Remote sensing has been widely used to estimate LAI and fAPAR. However, while many studies have examined both empirical and model-based relationships among LAI, fAPAR, and spectral vegetation indices (SVI) from remote sensing, few studies have systematically and empirically examined how relationships among these variables change over the growing season. In this study, we examine how and why seasonal-scale covariation differs among time series of remotely sensed SVIs and both LAI and fAPAR based on current understanding and theory. To do this we use newly available remote sensing data sets in combination with time series of in-situ measurements and a canopy radiative transfer model to analyze how seasonal variation in canopy and environmental conditions affect relationships among remotely sensed SVIs, LAI, and fAPAR at a temperate deciduous forest site in central Massachusetts. Our results show that accounting for seasonal variation in canopy shadowing, which is driven by variation in solar zenith angle, improved remote sensing-based estimates of LAI, fAPAR, and daily total APAR. Specifically, we show that the phenology of SVIs is strongly influenced by seasonal variation in near infrared (NIR) reflectance arising from systematic variation in the canopy shadow fraction that is independent of changes in LAI or fAPAR. Therefore, results of this work provide a refined basis for understanding how remote sensing can be used to monitor and model the phenology of LAI, fAPAR, APAR, and gross primary productivity in temperate deciduous forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of evapotranspiration and yield of maize: An inter-comparison among 41 maize models

Accurate simulation of crop water use (evapotranspiration, ET) can help crop growth models to assess the likely effects of climate change on future crop productivity, as well as being an aid for irrigation scheduling for today's growers. To determine how well maize (Zea mays L.) growth models can simulate ET, an initial inter-comparison study was conducted in 2019 under the umbrella of AgMIP (Agricultural Model Inter-Comparison and Improvement Project). Herein, we present results of a second inter-comparison study of 41 maize models that was conducted using more comprehensive datasets from two additional sites - Mead, Nebraska, USA and Bushland, Texas, USA. There were 20 treatment-years with varying irrigation levels over multiple seasons at both sites. ET was measured using eddy covariance at Mead and using large weighing lysimeters at Bushland. A wide range in ET rates was simulated among the models, yet several generally were able to simulate ET rates adequately. The ensemble median values were generally close to the observations, but a few of the models sometimes performed better than the median. Many of the models that did well at simulating ET for the Mead site did poorly for drier, windy days at the Bushland site, suggesting they need to improve how they handle humidity and wind. Additional variability came from the approaches used to simulate soil water evaporation. Fortunately, several models were identified that did well at simulating soil water evaporation, canopy transpiration, biomass accumulation, and grain yield. These models were older and have been widely used, which suggests that a larger number of users have tested these models over a wider range of conditions leading to their improvement. These revelations of the better approaches are leading to model improvements and more accurate simulations of ET.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗
DOE OSTI2025

Explainable machine learning to quantify the value of proximal remote sensing in latent energy flux estimation

Proximal remote sensing has the potential to provide critical information on vegetation biophysical factors that can predict land-atmosphere exchange of water and energy. Latent energy (LE) flux is traditionally estimated using process-based models which rely on vegetation parameters that change during the growing season. Data-driven models have the potential to address these issues by offering flexible predictor selection and more efficient utilization of the information in predictor sets. These models require careful choice of predictors to avoid redundancy and allow robust cross-validation. In this study we present a systematic and comprehensive evaluation of machine learning (ML) models to assess the capability of meteorological and proximal sensing data for predicting LE at a half-hourly temporal resolution across multiple growing seasons for an agricultural system. The results presented here demonstrate that a model using four environmental predictors in combination with two proximal sensing variables can capture 88 % of the variability in LE. ML models using only three predictors (one meteorological and two proximal remote sensing) captured 81 % of LE variability, offering the best trade-off between performance and complexity. An ML model utilizing only two predictors, one proximal remote sensing variable and downwelling radiation, captured 77 % of LE variability. These results demonstrate the power of proximal remote sensing and meteorological observations to estimate land-atmosphere water vapor exchange, providing a solution where more direct methods such as eddy covariance are not available and for evaluations of agronomic management and genotypic variations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Post-Combat-Injury Opioid Prescription and Alcohol Use Disorder in the Military

Previous studies have identified combat exposure and combat traumatic experience as problematic drinking risk factors. Increasing evidence suggests that opioid use increases the risk of alcohol use disorder. This study investigated the association between opioid prescription use after injury and (1) alcohol use disorder and (2) severity of alcohol use disorder among deployed military servicemembers. Deidentified health records data of 9,029 deployed servicemembers from a retrospective cohort study were analyzed. Data were randomly selected from the Department of Defense Trauma Registry and included servicemembers with combat injuries during deployment in Iraq or Afghanistan (2002–2016). Pharmacy records and International Classification of Diseases, Ninth and Tenth Revision diagnosis codes were used. Three groups were identified (no opioid prescription use, nonpersistent opioid prescription use, and persistent opioid prescription use) and were compared on the basis of alcohol use disorder risk using Cox proportional hazard models. Data analyses were performed in 2021. Of the 9,029 servicemembers with combat injury, 2,262 developed alcohol use disorder (1,322 developed severe alcohol use disorder). Compared with no opioid prescription use, increased alcohol use disorder risk was associated with persistent opioid prescription use, with a hazard ratio of 1.13 (95% CI=1.02, 1.26). After covariate adjustment, increased risk remained statistically significant (hazards ratio=1.24; 95% CI=1.10, 1.39). There was no significant difference in alcohol use disorder risk between no opioid prescription use and nonpersistent opioid prescription use. The risk of severe alcohol use disorder did not vary by opioid use among servicemembers with alcohol use disorder diagnosis. Furthermore, the findings of the study suggest that the incidence of alcohol use disorder was higher among injured servicemembers with persistent opioid prescription use than among those without opioid use. If replicated in prospective studies, the findings highlight the need for clinicians to consider the current and history of alcohol use of patients in initiating treatment involving opioids.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

An uncertainty quantification method relevant to material test reactors

Within material test reactor calculations, energy dependent flux and reaction rate uncertainties are typically not quantified when performing as-run analyses to determine the neutron field experienced by the experiment. When high fidelity Monte-Carlo codes are used in such analyses, straight forward methods to calculate output uncertainties are not available, instead expert opinion is used to postulate computational uncertainties. New methods to propagate uncertainties through these high fidelity simulations are available when sufficient computational power is available. A tool is developed here for sampling any part of an MCNP input from random distributions to determine output uncertainties based on those inputs. Another tool is developed to sample nuclear data cross-section in ACE format using multi-group nuclear data covariances. The Total Monte-Carlo Method and Gesellschaft für Anlagen-und Reaktorsicherheit method (GRS) are implemented and compared to one another as well as MCNP sensitivity and uncertainty calculations. The methods were applied to the Godiva critical sphere k-eigenvalue, the UAM pin-cell benchmark energy dependent flux and reaction rates, and the Advanced Test Reactor energy dependent flux within an experimental location. Furthermore, the two methods agree well, with GRS allowing for an order of magnitude speedup for reaction rate uncertainty calculations and several orders of magnitude for eigenvalue uncertainty calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

A new 181 Ta neutron resolved resonance region evaluation

A new 181 Ta neutron resolved resonance region evaluation has been performed from the thermal energy range up to approximately 2.5 keV. The R-matrix SAMMY code was used with the Reich–Moore approximation to evaluate resonance parameters from several experimental data sets. A Monte Carlo approach was used for resonance spin assignments and generating 59 small fictitious resonance levels which were shown to improve the cumulative level, Porter-Thomas, and Wigner distributions as compared to theoretical predictions. Covariance information was also generated for the entire resolved resonance region. Finally, the positive impact of the new evaluation was validated through benchmark calculations which were sensitive to the 181 Ta cross section and showed improvement in the reactivity bias for several benchmark cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Relativistic levels of mesic atoms

We revisit the derivation of the covariant two-body scalar-fermion equation with a Coulomb interaction, presented in a previous paper. We show that it can be given the formal aspect of a Dirac equation, but for the fact that the eigenvalue is also contained in one of the coefficients and thus it is not linearly included. The discussion of the boundary value problem is therefore different, although some properties of the Dirac equation can be recovered in an approximation bringing back to the concept of reduced mass. We discuss a mixed analytic-numerical method of solution which allows to obtain very accurate results and we calculate the lowest levels, states and QED corrections for pionic and kaonic atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

General tensor structure for electron scattering in terms of invariant responses

The general structure of semi-inclusive polarized electron scattering from polarized spin-1/2 targets is developed for use at all energy scales, from modest-energy nuclear physics applications to use in very high energy particle physics. The leptonic and hadronic tensors that enter in the formalism are constructed in a general covariant way in terms of kinematic factors that are frame dependent but model independent and invariant response functions which contain all of the model-dependent dynamics. In the process of developing the general problem the relationships to the conventional responses expressed in terms of the helicity components of the exchanged virtual photon are presented. For semi-inclusive electron scattering with polarized electrons and polarized spin-1/2 targets one finds that 18 invariant response functions are required, each depending on four Lorentz scalar invariants. Additionally it is shown how the semi-inclusive cross sections are related via integrations over the momentum of the selected coincidence particle and sums over open channels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Uranium integral fission product yields for a spectrally-shaped 14.1 MeV neutron source at the National Ignition Facility

This paper describes the experimental results for an energy tuning assembly created to modify the National Ignition Facility deuterium–tritium fusion neutron source into a notional thermonuclear and prompt fission neutron spectrum, which has applications in integral measurements, nuclear data benchmarks, and radiation effects on microelectronics. The Monte Carlo neutron transport utilized MCNP5 to estimate the ETA-modified fluence using the ENDF-B/VIII.0 and IRDFF-II continuous energy nuclear data libraries, and SCALE Sampler was used to estimate the systematic nuclear data covariance using ENDF-B/VII.1 and IRDFF-II in a 252-group structure. The experiment fielded eight activation foils and a highly enriched uranium sample. This provided fifteen reaction channels that were used in a forward-fit comparison to the modeled results and to unfold the neutron spectrum using STAYSL. Gamma-ray spectrometry was performed on the activation and highly enriched uranium foils, and the reduced Χ 2 between the modeled and experimental values was 1.21. The results from the STAYSL unfold, reduced Χ 2 = 1.62, indicated that the modeled neutron spectrum was achieved and the systematic nuclear data uncertainty associated with the neutron transport and activation product cross sections was representative of the experiment. Integral cumulative fission product yield data were collected for 37 mass chains with a combination of gamma-ray spectrometry and radiochemical analysis. Fission product analysis was generally in agreement with two models using a semi-empirical fit and the General Observables of Fission code, with the exception of mass chains 88, 109, 111, 112, 113, 129, 139, 142, 144, 151, and 156.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Gaussian process regression constrained by boundary value problems

We develop a framework for Gaussian processes regression constrained by boundary value problems. The framework may be applied to infer the solution of a well-posed boundary value problem with a known second-order differential operator and boundary conditions, but for which only scattered observations of the source term are available. Scattered observations of the solution may also be used in the regression. The framework combines co-kriging with the linear transformation of a Gaussian process together with the use of kernels given by spectral expansions in eigenfunctions of the boundary value problem. Furthermore, it benefits from a reduced-rank property of covariance matrices. We demonstrate that the resulting framework yields more accurate and stable solution inference as compared to physics-informed Gaussian process regression without boundary condition constraints.

42 ENGINEERING↗
DOE OSTI2025

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗
DOE OSTI2025

Estimating soybean yields from high-temporal-resolution multi-source data using deep learning

Accurate and timely crop yield prediction is crucial for ensuring food security and maintaining stable agricultural markets. In recent years, there has been a surge in interest in leveraging high-temporal-resolution, multi-source data for effective crop growth monitoring and yield estimation. A notable challenge arises from the difficulty in capturing the intricate interactions between variables across different time steps within these high-temporal-resolution time series datasets. This complexity hinders the reliable extraction of yield information from voluminous and often noisy datasets, especially during periods of extreme weather events. Here, in this study, we propose an Attention and Graph Isomorphism Network-enhanced Bi-directional Long Short-Term Memory network (AGB-LSTM) for estimating county-level soybean yield in the United States. This model integrates a diverse set of remote sensing data, including Near-Infrared Reflectance of Vegetation (NIRv), Sun-Induced chlorophyll Fluorescence (SIF), and Gross Primary Productivity (GPP), along with environmental covariates. The AGB-LSTM effectively leverages information related to crop yield from high-temporal-resolution time series data (5-days), achieving an accuracy of R²= 0.67 and rRMSE = 14.46%. This approach significantly outperforms traditional machine learning methods such as Random Forest (RF) (R²= 0.52, rRMSE = 17.36%) and Bi-LSTM (R²= 0.58, rRMSE = 16.17%). Sensitivity experiments with different time steps and ranges demonstrated that our model could accurately and stably predict yields 1 to 2 months before harvest. Moreover, data with a finer temporal resolution consistently improved prediction performance, resulting in an approximately 20% increase in and an approximately 20% decrease in rRMSE compared to using monthly composites. We also evaluated the robustness of the model under extreme climate events and observed strong performance (R²= 0.50, rRMSE = 21.32%). Finally, yield mapping for major soybean-producing regions in North America in 2023 revealed spatial patterns that closely matched USDA yield reports. Our findings suggest that the AGB-LSTM model is a promising and effective method for estimating yield and has notable potential for global crop yield forecasting.

Deep learning↗
DOE OSTI2020

Optimal sensor network design for multi-scale, time-varying differential algebraic equation systems: Application to an entrained-flow gasifier refractory brick

An algorithm for optimal sensor network design for multi-scale, time-varying differential algebraic equation systems with non-separable dynamics is presented. As the process is time-varying, an integral normalized posterior error covariance of a multi-scale filter is minimized to obtain the optimal sensor locations. For reducing the computational cost, an adaptive sampling rate approach is considered for the slowly-varying variables. The algorithm is applied to a smart refractory brick with embedded sensors as part of an entrained-flow gasifier. Thermistors and interdigital capacitors are considered as candidate measurement technologies for estimating temperature and slag penetration profile along the gasifier wall. Finally, when the optimal set of sensors obtained from the algorithm is used for estimating temperature and slag penetration profiles in a multi-scale Kalman filter framework, satisfactory estimates are obtained despite high measurement noise and model mismatch.

42 ENGINEERING↗
DOE OSTI2021

A comparison between residential relocation timing of Sydney and Chicago residents: A Bayesian survival analysis

We report that understanding households' behaviour in residential relocation timing is of great importance in the field of transport engineering and economics. This research aims to develop a residential relocation model by considering the potential dynamic impacts of other households' decisions and variables, including economic and demographic attributes, housing features, intra-household decision-making structures, travel mode choice, and other life-course attributes. A multivariate parametric survival model with both fixed and time-varying covariates is developed. To the best of the authors' knowledge, this study is the first paper in the literature of residential relocation timing to propose the use of a Bayesian model in contrast to the widely used classic frequentist approach and have conducted a discussion on its advantages. An emerging residential relocation dataset collected for two cities in Australia and the USA (Sydney and Chicago cities) has been used, which covers residence, vehicle ownership, occupation, education, economic and demographic attributes of respondents. A comprehensive comparison between the results of two cities and a comparison between two Bayesian and frequentist approaches are made. This study confirms the impact of life-course variables, intra-household decision-making behaviours, and sociodemographic attributes on home mobility. According to the model outputs, the accelerating or decelerating impact of explanatory variables on the relocation timing has been almost the same in the two cities. The Bayesian model was confirmed to have some advantages over the frequentist model, including being straightforward to interpret, availability of making inferences on the results, and ease of handling complex models, and optimisation convergence complexities.

99 GENERAL AND MISCELLANEOUS↗