Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density model”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Composition-Based Density Model for High Level Waste Glasses

In this report, the Savannah River National Laboratory (SRNL) provides a first-principles model capable of predicting the density of high-level waste (HLW) glass based on the glass composition. The model relies on the additivity of the specific volume of bound glass oxides to obtain a quantitative evaluation of the approximate glass density. A database of glass densities and compositions was compiled from various sources including previous Defense Waste Processing Facility (DWPF) variability study glasses, HLW studies from other institutions, and literature references. There are a total of 1104 glass compositions in the database with densities ranging from 2.224 g/cm 3 to 3.997 g/cm 3 . The breadth of the compositional coverage space is reported as well. A predictive method was investigated in which the glass density was calculated as a function of the weighted partial molar volumes of the individual glass oxides and the molar mass of the glass. This method yielded a high correlation between the predicted density and the measured density. In addition, a linear regression of these predictions indicated that the predicted value was equal to the measured within the bounds of analytical uncertainty. The following conclusions were drawn from the results of these analyses: 1. Measured glass density is highly correlated with glass composition (R 2 ≈ 0.95). 2. The correlation between glass density and glass composition can be expressed as either: a. A weighted sum of the individual oxide glass-component densities, or b. A weighted sum of the individual oxide glass-component partial molar volumes. 3. A linear regression of the predicted density via partial molar volume method to the measured density yields a slope of 0.9975 and an intercept of -0.0007 – the uncertainty interval associated with these coefficients contain one and zero, respectively. 4. The fitted coefficients and associated uncertainty intervals of the linear regression of the predicted density via the partial molar volume method to the measured density suggest that the predicted theoretical density is, within analysis uncertainty, equal to the measured density of the glass. It is recommended that glass density measurements included with the variability study for each sludge batch be discontinued. In their place, SRNL will calculate a bounding glass density value for the frit/sludge system using the model. It is also recommended that DWPF implement the density model as part of the fissile loading determination portion of the SME acceptability process. Details for implementation are provided in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Composition-Based Density Model for High Level Waste Glasses

In this report, the Savannah River National Laboratory (SRNL) provides a first-principles model capable of predicting the density of high-level waste (HLW) glass based on the glass composition. The model relies on the additivity of the specific volume of bound glass oxides to obtain a quantitative evaluation of the approximate glass density.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

An Extensible High Energy Density Modeling Tool for Extreme Regimes

The overreaching goal of the project is to develop a simulation capability to model an entire High Energy Density (HED) target heated by an X-ray Free Electron Laser (XFEL) beam including the surrounding material with a focus on target material associated with subsequent shots. Our primary objective is to study target fratricide problems, where radiation/debris from the current shot degrade the target for subsequent shots. This work will have a large impact on research being conducted at HED facilities where the quality of data is directly related to the effective repetition rate of the driver, which is often limited by target fratricide problems. This project focuses on a stream of closely spaced droplet targets that are appropriate for new very high repetition rate sources such as the ~1 MHz LCLS-II upgraded XFEL. The capability of the PISALE (Pacific Island Structured-AMR with ALE) code is extended to model this class of problems. We summarize the capability and improvements made to the HED branch of the PISALE code. We discuss our simulations of target fratricide for small (5 micron) hydrogen droplets heated by an XFEL beam at SLAC. We also compare simulations of the initial dynamics of larger (50 micron) water droplets, heated by the same XFEL beam, with experimental data. We briefly discuss training, dissemination, and related activities associated with this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Exchange current density model for the contact-determined current-voltage behavior of solar cells

An analytic expression for the current–voltage [J(V)] behavior of a solar cell as limited by equilibrium exchange current densities of both carriers at both contacts is derived. The partial currents at both contacts to a generic semiconductor absorber are assumed to be linearly proportional to the excess carrier concentration at the interface with the contacts (e.g., as with Schottky-like contacts). The assumption that the quasi-Fermi levels in the absorber are approximately flat leads to an algebraic solution for the applied voltage as a function of current, which is inverted to obtain the analytic J(V) curve. The J(V) curve reveals distinct behavior associated with electrons and holes, separately, and allows for the determination of all critical performance parameters. In particular, it demonstrates how the characteristic features of the J(V) curve depend on the relative rate at which a particular carrier (electron or hole) is collected at one contact vs the other, rather than the relative rate of electron vs hole collection at a single contact. Furthermore, the model provides a unified explanation of how majority carrier extraction limitations cause nonideal J(V) behaviors such as S-shaped curves and dark/light crossover (i.e., failure of superposition). The efficacy and limitations of the model when applied to Schottky-type and doped semiconductor contacts are discussed. The work serves as a theoretical guide to scientists studying solar cells that are thought to be primarily limited by their contacts.

Solar cells↗
DOE OSTI2022

Defect Thermodynamics and Transport Properties of Proton Conducting Oxide BaZr 1–x Y x O 3–δ (x ≤ 0.1) Guided by Density Functional Theory Modeling

Density functional theory-based thermodynamic modeling was performed to determine the effect of humidity and H 2 /O 2 gas pressure on the defect chemistry and transport properties of proton conducting oxide BaZr 1-x Y x O 3-δ (x ≤ 0.1) in the temperature range of 800–1200K, relevant for solid oxide fuel/electrolysis cell applications. The first-principles charge defect analysis was carried out to obtain the defect energetic ensembles as a function of Fermi level as well as the dependence of the hydration and oxidation reaction energies on complex defect-dopant configurations. It is shown that oxygen vacancies introduced to compensate for the extra Ba vacancies may cause a slight increase in the concentration of hydroxyl and proton species upon hydration of Ba deficient BaZr 1-x Y x O 3-δ , while the overall proton diffusivity is predicted to decrease with increased proton hopping barriers upon trapping near the Ba vacancies. The developed defect model is further demonstrated to be able to describe the bulk defect chemistry and transport properties of BaZr 0.9 Y 0.1 O 3-δ under the solid oxide cell operating conditions. Finally, the theoretical framework developed in this work further allows the inclusion of configurational, energetic and electronic characteristics of the defects that can be used as a complementary tool to experimental measurements for testing various mechanistic pathways or to provide essential mechanistic data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Defect Thermodynamics and Transport Properties of Perovskite and Fluorite Materials for Solid-Oxide and Proton Conducting Oxide Cells Evaluated Based on Density Functional Theory Modeling

Density functional theory based defect thermodynamic modeling was performed to determine the effect of humidity and H2/O2 gas pressure on various defect chemistry and transport properties of perovskite and fluorite oxides for solid-oxide and proton-conducting-oxide cell applications, with inclusion of the electronic-conducting oxides (as electrodes) and insulating oxides (as electrolytes). Automatic defect generation workflow and first-principles charged defect analysis were implemented on NETL Joule supercomputer for modeling defect equilibria and transport properties of insulating oxides as electrolytes in SOCs and proton-conducting ceramic cells. A GNU Octave defect model subroutines were developed to facilitate defect modeling of electronic conducting oxides in a wide range of operating conditions guided by modeling and experiments. The model includes the hydride defect formation reaction under reducing conditions and allows to incorporate nonstoichiometry effects on the defect thermodynamic parameters. The developed model serves as a platform to facilitate fundamental understanding of the defect thermodynamics in SOC oxide materials and can be used as a novel tool in computational materials screening for SOC and other energy applications.

Lee, Yueh-Lin↗
DOE OSTI2023

Modeling of the cold electron plasma density for radiation belt physics

This review focusses strictly on existing plasma density models, including ionospheric source models, empirical density models, physics-based and machine-learning density models. This review is framed in the context of radiation belt physics and space weather codes. The review is limited to the most commonly used models or to models recently developed and promising. A great variety of conditions is considered such as the magnetic local time variation, geomagnetic conditions, ionospheric source regions, radial and latitudinal dependence, and collisional vs. collisionless conditions. These models can serve to complement satellite observations of the electron plasma density when data are lacking, are for most of them commonly used in radiation belt physics simulations and can improve our understanding of the plasmasphere dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Uniform modelling of the stellar density of thirteen tidal streams within the Galactic halo

ABSTRACT We present the results of fitting a flexible stellar stream density model to a collection of thirteen streams around the Milky Way, using photometric data from DES, DECaLS, and Pan-STARRS. We construct density maps for each stream and characterize their tracks on the sky, width, and distance modulus curves along the length of each stream. We use these measurements to compute lengths and total luminosities of streams and identify substructures. Several streams show prominent substructures, such as stream broadening, gaps, large deviations of stream tracks, and sharp changes in stream densities. Examining the group of streams as a population, as expected we find that streams with globular cluster progenitors are typically narrower than those with dwarf galaxy progenitors, with streams around 100 pc wide showing overlap between the two populations. We also note the average luminosity of globular cluster streams is significantly lower than the typical luminosity of intact globular clusters. The likely explanation is that observed globular cluster streams preferentially come from lower luminosity and lower density clusters. The stream measurements done in a uniform manner presented here will be helpful for more detailed stream studies such as identifying candidate stream members for spectroscopic follow up and stellar stream dynamical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Projected density matrix embedding theory with applications to the two-dimensional Hubbard model

Density matrix embedding theory (DMET) is a quantum embedding theory for strongly correlated systems. From a computational perspective, one bottleneck in DMET is the optimization of the correlation potential to achieve self-consistency, especially for heterogeneous systems of large size. We introduce a new approach, called projected DMET (p-DMET), which achieves self-consistency without needing to optimize the correlation potential. We demonstrate the performance of p-DMET on the two-dimensional Hubbard model.

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

Diffusion-Model-Assisted Supervised Learning of Generative Models for Density Estimation

Here, we present a supervised learning framework of training generative models for density estimation. Generative models, including generative adversarial networks (GANs), normalizing flows, and variational auto-encoders (VAEs), are usually considered as unsupervised learning models, because labeled data are usually unavailable for training. Despite the success of the generative models, there are several issues with the unsupervised training, e.g., requirement of reversible architectures, vanishing gradients, and training instability. To enable supervised learning in generative models, we utilize the score-based diffusion model to generate labeled data. Unlike existing diffusion models that train neural networks to learn the score function, we develop a training-free score estimation method. This approach uses mini-batch-based Monte Carlo estimators to directly approximate the score function at any spatial-temporal location in solving an ordinary differential equation (ODE), corresponding to the reverse-time stochastic differential equation (SDE). This approach can offer both high accuracy and substantial time savings in neural network training. Once the labeled data are generated, we can train a simple, fully connected neural network to learn the generative model in the supervised manner. Compared with existing normalizing flow models, our method does not require the use of reversible neural networks and avoids the computation of the Jacobian matrix. Compared with existing diffusion models, our method does not need to solve the reverse-time SDE to generate new samples. As a result, the sampling efficiency is significantly improved. We demonstrate the performance of our method by applying it to a set of 2D datasets as well as real data from the University of California Irvine (UCI) repository.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multiplicative noise and the diffusion of conserved densities

Stochastic fluid dynamics governs the long time tails of hydrodynamic correlation functions, and the critical slowing down of relaxation phenomena in the vicinity of a critical point in the phase diagram. In this work we study the role of multiplicative noise in stochastic fluid dynamics. Multiplicative noise arises from the dependence of transport coefficients, such as the diffusion constants for charge and momentum, on fluctuating hydrodynamic variables. We study long time tails and relaxation in the diffusion of a conserved density (model B), and a conserved density coupled to the transverse momentum density (model H). Careful attention is paid to fluctuation-dissipation relations. We observe that multiplicative noise contributes at the same order as non-linear interactions in model B, but is a higher order correction to the relaxation of a scalar density and the tail of the stress tensor correlation function in model H.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Density-Matrix Model for Photon-Driven Transport in Quantum Cascade Lasers

In this report we develop a time-dependent density-matrix model to study photon-assisted (PA) electron transport in quantum cascade lasers. The Markovian equation of motion for the density matrix in the presence of an optical field is solved for an arbitrary field amplitude. Level-broadening terms emerge from microscopic Hamiltonians and supplant the need for empirical parameters that are often employed in related approaches. We show that, in quantum cascade lasers with diagonal design, photon resonances have a pronounced impact on electron dynamics around and above the lasing threshold, an effect that stems from the large spatial separation between the upper and lower lasing states. With the inclusion of PA tunneling, the calculated current density and output power are in good agreement with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DeePKS-kit: A package for developing machine learning-based chemically accurate energy and density functional models

We introduce DeePKS-kit, an open-source software package for developing machine learning based energy and density functional models. DeePKS-kit is interfaced with PyTorch, an open-source machine learning library, and PySCF, an ab initio computational chemistry program that provides simple and customized tools for developing quantum chemistry codes. It supports the DeePHF and DeePKS methods. In addition to explaining the details in the methodology and the software, we also provide an example of developing a chemically accurate model for water clusters.

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

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Dislocation Density-Based Model for Flow Behavior of a Near-α Titanium Alloy Considering Effects of Initial Lamellar Thickness

A dislocation density-based constitutive model that considered the effects of initial lamellar thickness was proposed to predict the flow behavior of a near-α high-temperature Ti-5.4Al-3.7Sn-3.3Zr-0.5Mo-0.4Si alloy. The evolution equations of dislocation density were improved by introducing the influence of the dynamic globularization of lamellar microstructure with different initial thicknesses. The evolution equations and Hall-Petch behavior were coupled into the constitutive model, which was composed of athermal and thermal stresses. According to the established constitutive equations, the variation and contributions of corresponding stress components were analyzed. The calculated dislocation densities were compared with the experimental values obtained by electron backscattered diffraction data to verify the model. Moreover, the flow stresses of the studied alloy with different initial lamellar microstructures under various deformation conditions were predicted. Good agreements between the predicted and experimental results were obtained. Additionally, steady-state flow behavior was further analyzed and the relationship between the calculated steady-state flow stress and the equiaxed α-phase size was quantitatively obtained. The results indicated the validity of the proposed constitutive model for predicting the flow behavior of the alloy with different initial lamellar microstructures over a range of deformation temperatures and strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An SYK-inspired model with density–density interactions: Spectral & wave function statistics, Green’s function and phase diagram

The Sachdev–Ye–Kitaev (SYK) model is a rare example of a strongly-interacting system that is analytically tractable. Tractability arises because the model is largely structureless by design and therefore artificial: while the interaction is restricted to two-body terms, interaction matrix elements are “randomized” and therefore the corresponding interaction operator does not commute with the local density. Unlike conventional density–density-type interactions, the SYK-interaction is, in this sense, not integrable. We here investigate a variant of the (complex) SYK model, which restores this integrability. It features a randomized single-body term and a density–density-type interaction. We present numerical investigations suggesting that the model exhibits two integrable phases separated by several intermediate phases including a chaotic one. The chaotic phase carries several characteristic SYK-signatures including in the spectral statistics and the frequency scaling of the Green’s function and therefore should be adiabatically connected to the non-Fermi liquid phase of the original SYK model. Thus, our model Hamiltonian provides a bridge from the SYK-model towards microscopic realism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗