Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Numerical methods for reduction of topside ionograms

Several alternative methods for solving the group height equation are presented. Three of these are now in operation at Ames Research Center and use data contained in a single ionogram trace. From the data an electron density profile N(h) is computed. If the ionogram also exhibits other traces, reverse ionogram traces are computed, using the N(h) profile, for comparison with the redundant data. When agreement is poor, the initial data trace is reinterpreted, another N(h) profile computed, and the reverse traces generated once again. This process is repeated until a desired degree of consistency is achieved. To reduce the necessity for human intervention and eliminate decision making required in conjunction with the preceding methods, a method is proposed that accepts as input, all data from a single ionogram. In general, no electron density function will satisfy these data exactly, but a best N(h) profile can be computed. Finally, a method is described that eliminates the need to assume that the ionosphere is spherically stratified. Horizontal gradients in electron density are detected and accounted for by processing several ionograms from the same satellite pass simultaneously. This idea is derived as an extension of one of the basic methods.

Mcculley, L.↗
DOE OSTI2023

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗
DOE OSTI2024

Lower hybrid drift instability in nonthermal plasmas

Lower hybrid drift instability (LHDI) is driven by the cross-field current and operates in the vicinity of the lower-hybrid frequency, between the ion- and electron-gyro frequencies, and with wavelengths between the electron and ion thermal gyro radii. The free energy source that drives this instability resides in the density gradient associated with an inhomogeneous plasma. The existing literature on LHDI assumes that the charged particle distribution function is given by a Maxwellian form, but the space plasma is pervasively observed to feature nonthermal characteristics. This paper extends the theory of LHDI to nonthermal plasmas. The generalized theory of LHDI is, thus, applicable to various space plasma environments characterized by nonthermal plasma velocity distribution functions.

Arya, Neetasha (ORCID:0009000236075605)↗
DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗
DOE OSTI2023

Accurate Kohn-Sham auxiliary system from the ground-state density of solids

The Kohn-Sham (KS) system is an auxiliary system whose effective potential is unknown in most cases. It is in principle determined by the ground-state density and it has been found numerically for some low-dimensional systems by inverting the KS equations starting from a given accurate density. For solids, only approximate results are available. In this work, we determine accurate exchange-correlation (xc) potentials for Si and NaCl using the ground-state densities obtained from auxiliary field quantum Monte Carlo calculations. We show that these xc potentials can be rationalized as an ensemble of a few local functions of the density, whose form depends on the specific environment and can be well characterized by the gradient of the density and the local kinetic energy density. Further, the KS band structure can be obtained with high accuracy. The true KS band gap turns out to be larger than the prediction of the local density approximation, but significantly smaller than the measurable photoemission gap, which confirms previous estimates. Finally, our findings show that the conjecture that very different xc potentials can lead to very similar densities and other KS observables is true also in solids, which questions the meaning of details of the potentials and, at the same time, confirms the stability of the KS system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Adsorption and Growth of Organic Materials on Metal Surfaces

This work investigated extensively the role of van der Waals (vdWs) interactions, when incorporated self-consistently into density functional theory (DFT). The investigations follow two avenues: 1) screening of several vdWs schemes at the generalized gradient approximation (GGA), using adsorption of small and large organic molecules on a variety of metal surfaces; and 2) simulating the adsorption of 2D materials such as silicene and phosphorene on metal surfaces. This second avenue is done in connection were subject of experimental investigations. These investigations led us to the conclusion that, if the inclusion of the vdWs interactions improved the nature of the adsorption and brought some excellent agreements with the experimental observations of the adsorption geometries, there is still a lot to do when it comes to the changes to the electronic structure upon adsorption.

2D materials.↗
NASA NTRS1993

Characterization of the thermal conductivity for Advanced Toughened Uni-piece Fibrous Insulations

Advanced Toughened Uni-piece Fibrous Insulations (TUFI) is discussed in terms of their thermal response to an arc-jet air stream. A modification of the existing Ames thermal conductivity program to predict the thermal response of these functionally gradient materials is described in the paper. The modified program was used to evaluate the effect of density, surface porosity, and density gradient through the TUFI materials on the thermal response of these insulations. Predictions using a finite-difference code and calculated thermal conductivity values from the modified program were compared with in-depth temperature measurements taken from TUFI insulations during short exposures to arc-jet hypersonic air streams.

Stewart, David A.↗
DOE OSTI2021

Cryogenic operation of planar ultrananocrystalline diamond field emission source in SRF injector

Here, we report a demonstration of electron beam generation in a superconducting radio frequency (SRF) Teraelectronvolt Energy Superconducting Linear Accelerator 1.3 GHz gun equipped with a field emission cathode and operated at 2 K. The cathode is a submicrometer film of nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] deposited atop a Nb RRR300 cathode plug that is 4 mm in diameter. Measurements of the output current showed that it increased exponentially as a function of the gun gradient. Our results demonstrate a feasible path toward simplified fully cryogenic SRF injector technology. At a maximal gradient of 0.9 MV/m, a current density of 22 μA/cm2 was obtained [calculated as a 2.75 μA current over the total cathode surface area of π(2 mm)2]. One important finding is that the electron emitter made of (N)UNCD, a material long been known as a highly efficient field emission material, demonstrated a record low turn-on gradient of 0.6 MV/m. A hypothesis explaining this behavior is proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Interplanetary stream magnetism - Kinematic effects

The particle density and the magnetic-field intensity and direction are calculated for volume elements of the solar wind as a function of the initial magnetic-field direction and the initial speed gradient. It is assumed that the velocity is constant and radial. These assumptions are approximately valid between about 0.1 and 1.0 AU for many streams. Time profiles of the particle density, field intensity, and velocity are calculated for corotating streams, neglecting effects of pressure gradients. The compression and rarefaction of the magnetic field depend sensitively on the initial field direction. By averaging over a typical stream, it is found that the average radial field intensity is inversely proportional to the square of the heliocentric distance, whereas the average intensity in the direction of the planets' motion does not vary in a simple way, consistent with deep space observations. Changes of field direction may be very large, depending on the initial angle; but when the initial angle at 0.1 AU is such that the base of the field line corotates with the sun, the spiral angle is the preferred direction at 1 AU. The theory is also applicable to nonstationary flows.

Burlaga, L. F.↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1986

Analyses of convective flows and spatial gradients in energetic ion observations

Based on the measured anisotropy of three-dimensional energetic ion distributions obtained from the Medium Energy Particle Instrument experiment of the ISEE 1 spacecraft, bulk flow, epsilon x B convection, and pressure gradients in magnetospheric plasmas are measured with two techniques. In the first technique, using a Compton-Getting transform into the appropriate rest frame in an iterative fashion and converging to the unique velocity that removes the first-order anisotropy from the flux angular distribution, the actual measured angular distribution, and spectrum at each angle, are employed. The second technique uses a first moment of the distribution function in which the shift in velocity space of the distribution function is obtained at two phase space density values. The convection contribution to the velocity shift is assumed to be perpendicular to the pressure gradient contribution if an MHD generator condition is assumed, rendering the two shifts separable, and a gradient scale and a lower limit to the convection electric field are derived in one calculation.

Mitchell, D. G.↗
NASA NTRS1984

Density-functional calculations of the surface tension of liquid Al and Na

Calculations of the surface tensions of liquid Al and Na are described using the full ionic density functional formalism of Wood and Stroud (1983). Surface tensions are in good agreement with experiment in both cases, with results substantially better for Al than those found previously in the gradient approximation. Preliminary minimization with respect to surface profile leads to an oscillatory profile superimposed on a nearly steplike ionic density disribution; the oscillations have a wavellength of about a hardsphere diameter.

Stroud, D.↗
NASA NTRS1987

A thermodynamic model of the Hg(0.8)Cd(0.2)Te-iodine transport system. I - Te-saturated source material. II - Source material composition within the homogeneity range

A thermochemical analysis of the Hg(0.8)Cd(0.2)Te-iodine vapor transport system is presented, and theoretical calculations of diffusion-controlled mass transport rates are made. The predicted mass fluxes are compared with experimental data obtained from transport experiments under vertical, stabilizing conditions reported earlier and with results of additional transport experiments conducted during the present study. Experimental mass transport rate studies of the transport system for fixed amount of excess Hg as a function of transport agent pressure are presented. The mass fluxes are determined for the vertical, stabilizing orientation of the density gradient relative to the gravitational vector. In order to compare experimental mass transport rates with computed values, the thermochemical analysis is extended to take the formation of Hg vacancies in the above compound into account along with their effect on the partial pressure of the system.

Chandra, D.↗