Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Global brain activity and its coupling with cerebrospinal fluid flow is related to tau pathology

Abstract INTRODUCTION Factors responsible for the deposition of pathological tau in the brain are incompletely understood. This study links macroscale tau deposition in the human brain to cerebrospinal fluid (CSF) flow dynamics using resting‐state functional magnetic resonance imaging (rsfMRI). METHODS Low‐frequency (< 0.1 Hz) resting‐state global brain activity is coupled with CSF flow and potentially reflects CSF dynamics‐related clearance. We examined the correlation between rsfMRI measures of CSF inflow and global activity (gBOLD–CSF coupling) as a predictor, interacting with amyloid beta (Aβ), of tau and cortical thickness (dependent variables) across Alzheimer's Disease Neuroimaging Initiative (ADNI) participants from cognitively unimpaired through mild cognitive impairment (MCI) and Alzheimer's disease (AD). RESULTS Tau deposition in Aβ+ participants, accompanied by cortical thinning and cognitive decline, is associated with decreased gBOLD–CSF coupling. Tau mediates the relationship between coupling and thickness. DISCUSSION Findings suggest that resting‐state global brain activity and CSF movements comodulate Alzheimer's tau deposition, presumably related to CSF clearance. Highlights A non‐invasive functional magnetic resonance imaging (fMRI) assessment of a CSF clearance‐related process is carried out. Global brain activity is coupled with CSF inflow in human fMRI during resting state. Global fMRI–CSF coupling is correlated with tau in Alzheimer's disease (AD). This coupling measure is also associated with cortical thickness, mediated by tau.

Neurosciences & Neurology↗
DOE OSTI2022

Pervasive alterations to snow-dominated ecosystem functions under climate change

Climate change projections consistently demonstrate that warming temperatures and dwindling seasonal snowpack will elicit cascading effects on ecosystem function and water resource availability. Despite this consensus, little is known about potential changes in the variability of ecohydrological conditions, which is also required to inform climate change adaptation and mitigation strategies. Considering potential changes in ecohydrological variability is critical to evaluating the emergence of trends, assessing the likelihood of extreme events such as floods and droughts, and identifying when tipping points may be reached that fundamentally alter ecohydrological function. Using a single-model Large Ensemble with sophisticated terrestrial ecosystem representation, we characterize projected changes in the mean state and variability of ecohydrological processes in historically snow-dominated regions of the Northern Hemisphere. Widespread snowpack reductions, earlier snowmelt timing, longer growing seasons, drier soils, and increased fire risk are projected for this century under a high-emissions scenario. In addition to these changes in the mean state, increased variability in winter snowmelt will increase growing-season water deficits and increase the stochasticity of runoff. Thus, with warming, declining snowpack loses its dependable buffering capacity so that runoff quantity and timing more closely reflect the episodic characteristics of precipitation. This results in a declining predictability of annual runoff from maximum snow water equivalent, which has critical implications for ecosystem stress and water resource management. So our results suggest that there is a strong likelihood of pervasive alterations to ecohydrological function that may be expected with climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Qumode transfer between continuous- and discrete-variable devices

Transferring quantum information between different types of quantum hardware is crucial for integrated quantum technology. In particular, converting information between continuous variable (CV) and discrete variable (DV) devices enables many applications in quantum networking, quantum sensing, quantum machine learning, and quantum computing. This paper addresses the transfer of CV-encoded information between CV and DV devices. We present a resource-efficient method for encoding CV states and implementing CV gates on DV devices, as well as two measurement-based protocols for transferring CV states between CV and DV devices. The success probability of the transfer protocols depends on the measurement outcome and can be increased to near-deterministic values by adding ancillary qubits to the DV devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

High-Resolution LiDAR Observations for Coupled Effects of Dry-Air Entrainment and Haze-Cloud Interactions on Cloud Vertical Structure

This study demonstrates the high-resolution profiling of cloud microphysics in a laboratory chamber using Time-Correlated Single Photon Counting (TCSPC) LiDAR. We present a novel retrieval method to derive vertical extinction (𝜎) profiles, constrained by in situ measurements, to diagnose responses to dry-air entrainment. In clean clouds, the LiDAR signals and retrieved 𝜎 remain relatively uniform, with entrainment effects confined to the upper layer. In contrast, polluted clouds exhibit strong vertical variability and a transition to a water-vapor-limited state. Entrainment significantly enhances the haze number concentration (𝑁ℎ), particularly near the bottom, creating highly height-dependent extinction profiles for polluted clouds. Our results highlight the capability of high-resolution LiDAR in capturing fine-scale vertical inhomogeneities. This approach provides a robust framework for quantifying how aerosol loading modulates entrainment sensitivity, offering new insights into the transition between buffered and water-vapor-limited regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of a point kinetics subroutine for Molten Salt reactors in RELAP5-3D

RELAP5-3D has historically focused on accurate modeling of design-basis and beyond-design-basis accidents in solid fuel light-water reactors. The assumption used in solid fuel reactors that delayed neutron precursors remain at the location in which they were produced is inherently incorrect for a molten salt reactor with flowing fuel. Therefore, new methods for modeling MSR kinetics were required. Subroutine used in RELAP5-3D were updated to analyze point kinetics in a molten salt reactor during both steady state and transient analyses. The phenomena related to MSR kinetics that were added to RELAP5-3D include arrays to store previous values for independent and dependent variables, arrays to add terms to existing kinetics equations, and a new reactivity bias term. A custom MATLAB assessment code was used to establish a suite of verification tests. This capability is undergoing final documentation and testing for public release with a future version of RELAP5-3D.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

The value of fusion energy to a decarbonized United States electric grid

Fusion could be a part of future decarbonized electricity systems, but it will need to compete with other technologies. In particular, pulsed tokamak plants have a unique operational mode, and evaluating which characteristics make them economically competitive can help select between design pathways. Using a capacity expansion and operations model, here we determined cost thresholds for pulsed tokamaks to reach a range of penetration levels in a future decarbonized US Eastern Interconnection. The required capital cost to reach a fusion capacity of 100 GW varied from $2,700 to $7,500 kW -1 , and the equilibrium penetration increases rapidly with decreasing cost. The value per unit power capacity depends on the variable operational cost and on the cost of its competition, particularly fission, much more than on the pulse cycle parameters. These findings can therefore provide initial cost targets for fusion more generally in the United States.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Measurement report: Role of organic coating and chemical composition on ice nucleation potential of atmospheric particles in European Arctic

Understanding the ice nucleation (IN) potential of Arctic aerosols is critical for predicting their influence on cloud formation and water cycles in this vulnerable region. This study investigates the role of particle composition, organic coatings, and aerosol sources in modulating ice nucleating particle (INPs) abundance across five aerosol samples collected at the Gruvebadet Observatory Station in Ny-Ålesund, Svalbard. The IN potential of Arctic aerosol particles was studied by investigating chemical, morphological, and INP abundance measurements. Single-particle analyses revealed distinct differences in mixing state, organic volume fraction (OVF), and organic coating morphology across samples. OVF distributions were linked to particle origin, with marine-influenced Na-rich particles often exhibiting thin organic coatings, while long-range transported particles showed thicker organic coatings. Biogenic contributions, though variable, were linked to heat-sensitive INPs, suggesting a role for labile biological macromolecules under certain meteorological conditions. Spearman rank correlation analysis between particle composition and immersion-mode INP concentrations at two freezing temperatures indicated that organic-rich and Na-rich particles were positively associated with enhanced INP abundance. However, discrepancies in INP abundance were observed for particles with thicker organic coatings, where the morphological configuration of the organic material may play a role. The results highlight that Arctic INP variability is governed not only by chemical composition but also by the morphological configuration of organic material, which can either enhance or inhibit ice nucleation depending on its abundance, distribution, thickness, and mixing state. These findings underscore the combined influence of source regions, atmospheric processing, and organic–inorganic interactions in shaping Arctic aerosol freezing behavior.

Lata, Nurun Nahar [Pacific Northwest National Labo↗
DOE OSTI2025

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗
DOE OSTI2022

Continuous conditional generative adversarial networks for data-driven solutions of poroelasticity with heterogeneous material properties

Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Competition between magnetic order and charge localization in Na 2 IrO 3 thin crystal devices thin crystal devices

Spin orbit assisted Mott insulators such as sodium iridate (Na 2 IrO 3 ) have been an important subject of study in recent years. In these materials, the interplay of electronic correlations, spin-orbit coupling, crystal field effects, and a honeycomb arrangement of ions bring exciting ground states, predicted in the frame of the Kitaev model. The insulating character of Na 2 IrO 3 has hampered its integration to an electronic device, desirable for applications, such as the manipulation of quasiparticles interesting for topological quantum computing. Here we show through electronic transport measurements supported by angle-resolved photoemission spectroscopy (ARPES) experiments, that electronic transport in Na 2 IrO 3 is ruled by variable range hopping and it is strongly dependent on the magnetic ordering transition known for bulk Na 2 IrO 3 , as well as on external electric fields. Furthermore, electronic transport measurements allow us to deduce a value for the localization length and the density of states in our Na 2 IrO 3 thin crystal devices, and offer an alternative approach to study insulating 2D-materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A comparison of model validation approaches for echo state networks using climate model replicates

As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. Here, in this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A directional microcrack damage theory for brittle solids based on continuous hyperplasticity

We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.

Materials Science↗
DOE OSTI2025

Strength and ductility of additively manufactured 316L stainless steel: Impact of neutron irradiation and data variability

Here, this article presents the mechanical properties of additively manufactured (AM) 316L stainless steel processed via the laser powder bed fusion (LPBF) method, focusing on the effects of neutron irradiation on mechanical properties and the variability in strength and ductility data. The rapid melting-solidification process and multiple heating-cooling cycles inherent in LPBF typically result in a fine, metastable microstructure with significant local variability. AM 316L builds of varying thicknesses were fabricated, and SS-J3 miniature tensile specimens were machined from six different locations. These specimens were irradiated with fast neutrons to doses of 2 and 10 dpa at target temperatures of 300 °C and 600 °C. Post-irradiation tensile tests were conducted at room temperature, 300 °C, and 600 °C. Compared to conventional 316L stainless steel, AM 316L exhibited higher initial strength but lower ductility. Irradiation at 300 °C caused significant hardening and prompt necking at yield, with limited uniform ductility, although embrittlement was not observed up to 10 dpa. While neutron irradiation, particularly at 600 °C, increased the variability in strength and ductility data, no clear dependence of mechanical properties on build thickness or sampling location was found—contrary to the conventional perception that AM materials may exhibit high property variability. Furthermore, we observed that the variability in property data for LPBF-processed 316L was relatively low compared to that of wrought 316L stainless steel. This reduced variability in AM 316L steel may be attributed to its highly metastable, stress-containing microstructure, which is discussed in the context of general tensile property variations.

Additively manufactured 316L stainless steel↗
DOE OSTI2025

Impact of amorphous pockets on displacement damage evolution in silicon

Silicon has long been known to exhibit amorphization in response to heavy particle bombardment. For doses below the total amorphization threshold, partial amorphization is observed in the form of scattered amorphous pockets. While extensive research has gone into modeling the formation and evolution of amorphous pockets in response to irradiation, no studies yet investigate their impact on the evolution of other damage such as interstitial supersaturation and clustering. In this study, we survey the impact of amorphous pockets on defect evolution in silicon when treated as static sinks. MD is first used to show that amorphous pockets provide energetically favorable sites for point defects relative to the crystalline bulk, supporting the hypothesis that they act as sinks. A 0-D cluster dynamics model is then constructed, taking an interstitial clustering model from the literature and including amorphous pockets as a sink species. We conduct our survey for temperatures between 30 and 400 °C and sink strengths between 1 to 6 x 10 10 cm −2 . Both implantation- and radiation-induced damage states are investigated using interstitial and vacancy concentrations as initial condition variables. We find that, due to the differing migration rates of the interstitial and the vacancy, amorphous pockets have a non-monotonic impact on the final damage state depending on the effective sink strength of the amorphous pockets, resulting in increased damage formation in regimes of intermediate amorphization. In conclusion, this result emphasizes the important role of amorphous pockets in governing the evolution of damage in partially amorphized crystalline materials.

36 MATERIALS SCIENCE↗