Search NASA⌕ Search

SEARCH · Search NASA

Results for “density”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cluster model of C 12 in the density functional theory framework

We employ the constrained density functional theory to investigate cluster phenomena for the 12 C nucleus. The proton and neutron densities are generated from the placement of three 4 He nuclei (α particles) geometrically. These densities are then used in a density constrained Hartree-Fock calculation that produces an antisymmetrized state with the same densities through energy minimization. In the calculations no a priori analytic form for the single-particle states is assumed and the full energy density functional is utilized. The geometrical scan of the energy landscape provides the ground state of 12 C as an equilateral triangular configuration of three αs with molecular bond like structures. The use of the nucleon localization function provides further insight to these configurations. One can conclude that these configurations are a hybrid between a pure mean-field and a pure α particle condensate. Furthermore, this development could facilitate density functional theory based fusion calculations with a more realistic 12 C ground state.

6 ≤ A ≤ 19↗
DOE OSTI2019

Selective laser melting of 304L stainless steel: Role of volumetric energy density on the microstructure, texture and mechanical properties

The role of volumetric energy density on the microstructural evolution, texture and mechanical properties of 304L stainless steel parts additively manufactured via selective laser melting process is investigated. 304L is chosen because it is a potential candidate to be used as a matrix in a metal matrix composite with nanoparticles dispersion for energy and high temperature applications. The highest relative density of 99 %±0.5 was achieved using a volumetric energy density of 1400 J/mm 3 . Both XRD analysis and Scheil simulation revealed the presence of a small trace of the delta ferrite phase, due to rapid solidification within the austenitic matrix of 304L. Here, a fine cellular substructure ranged between 0.4–1.8 μm, was detected across different energy density values. At the highest energy density value, a strong texture in the direction of [100] was identified. At lower energy density values, multicomponent texture was found due to high nucleation rate and the existing defects. Yield strength, ultimate tensile strength, and microhardness of samples with a relative density of 99 % were measured to be 540 ± 15 MPa, 660 ± 20 MPa and 254 ± 7 HV, respectively and higher than mechanical properties of conventionally manufactured 304L stainless steel. Heat treatment of the laser melted 304L at 1200 °C for 2 h, resulted in the nucleation of recrystallized equiaxed grains followed by a decrease in microhardness value from 233 ± 3 HV to 208 ± 8 HV due to disappearance of cellular substructure.

304L stainless steel↗
DOE OSTI2022

Turbine scale and siting considerations in wind plant layout optimization and implications for capacity density

Improvements in wind energy technology, reduced costs, and ambitious clean energy goals have led to projections of high wind contribution in coming years. Developing methodologies to design wind plants with a variety of siting constraints and turbine sizes helps enable high wind penetration, and gain a better understanding of how wind plants are sensitive to setback constraints and turbine design. In this paper, we present a two-step optimization method to simultaneously determine the optimal number of turbines and their locations in a wind plant domain divided into many small, discrete parcels. We present the optimized performance metrics of a wind plant optimized with different turbine sizes and ratings, and with different siting restrictions within the wind plant. Our results indicate that taller and larger turbines are more sensitive to increasing siting constraints. We also compare the optimal wind plant layouts and performance for wind plants optimized for minimum COE and maximum profit. Wind plants optimized for profit had 130%-190% of the capacity of plants optimized for COE, which demonstrates that the optimal results are greatly affected by the objective function, which should be carefully considered. Finally, in this paper we demonstrate the effect of increasing siting constraints on wind plant capacity density, and how the results change when different land areas are used to calculate capacity density. When using the entire wind plant boundary area to determine capacity density, increasing siting constraints decreases the capacity density. However, when we only use the available area (the area left after removing the siting constraints) to calculate the capacity density, increasing the siting constraints increases capacity density. This is a critical insight because of how capacity density is typically defined and used in research, and has important implications for assessment of technical potential and capacity expansion modeling, as well as future wind deployment potential.

17 WIND ENERGY↗
DOE OSTI2021

Effect of edge plasma density on hot spot in LHCD plasma in EAST

Hot spots are a serious challenge limiting long pulse operation with lower hybrid current drive (LHCD) in tokamaks. Here, in order to mitigate hot spots in the guard limiter and improve LHCD capability, the effect of edge plasma density on inducing hot spots and current drive has been studied in EAST. The temperature in the guard limiter of the LH antenna, inducing hot spot directly, increases with edge density and LH power. Studies show that the hot spot is mainly ascribed to the heat flux in front of LH antenna. Further simulation indicates that such spots correspond to the peak position of edge density due to local LH electric field. In addition, due to the stronger parametric instability (PI) behavior in the case of higher edge density, the current drive capability decreases with edge density. Strike-point splitting behaviour appears as density increase, in agreement with current profiles in the edge region and the reduction of total driven current, suggesting that more power is deposited in the edge region, which then contributes more to hot spot formation. These studies offer one possible idea to optimize the edge density so as to satisfy the coupling, mitigate the heat flux in the guard limiter, and improve current drive capability in fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toward a QUBO-Based Density Matrix Electronic Structure Method

Density matrix electronic structure theory is used in many quantum chemistry methods to “alleviate” the computational cost that arises from directly using wave functions. Although density matrix based methods are computationally more efficient than wave function based methods, significant computational effort is involved. Because the Schrödinger equation needs to be solved as an eigenvalue problem, the time-to-solution scales cubically with the system size in mean-field type approaches such as Hartree–Fock and density functional theory and is solved as many times in order to reach charge or field self-consistency. We hereby propose and study a method to compute the density matrix by using a quadratic unconstrained binary optimization (QUBO) solver. This method could be useful to solve the problem with quantum computers and, more specifically, quantum annealers. Our proposed approach is based on a direct construction of the density matrix using a QUBO eigensolver. We explore the main parameters of the algorithm focusing on precision and efficiency. We show that, while direct construction of the density matrix using a QUBO formulation is possible, the efficiency and precision have room for improvement. Moreover, calculations performed with quantum annealing on D-Wave’s new Advantage quantum computer are compared with results obtained with classical simulated annealing, further highlighting some problems of the proposed method. Finally, we also suggest alternative methods that could lead to a more efficient QUBO-based density matrix construction.

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

Local Spin Density Approximation Strongly Improved by a Better-Informed Local Scaling of Its Self-Interaction Correction

The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.

Approximation↗
DOE OSTI2020

Two-Dimensional Molecular Charge Density Waves in Single-Layer-Thick Islands of a Dirac Fermion System

Charge density waves have been intensely studied in inorganic materials such as transition metal dichalcogenides however its counterpart in organic materials has yet to be explored in detail. Here we report the finding of a robust two dimensional charge density waves in molecular layers formed by α-(BEDT-TTF) 2 -I 3 on a Ag(111) surface. Low temperature scanning tunnelling microscopy images of a multi-layer thick α-(BEDT-TTF) 2 -I 3 on Ag(111) substrate reveal coexistence of 5a0 x 5a0 and √31 a0 x √31 a0 R9° charge density wave patterns commensurate with the underlying molecular lattice at 80 K. Both charge density wave patterns remain in nano-size molecular islands with just a single constituent molecular-layer thickness at 80 K and 5 K. Local tunneling spectroscopy measurements reveal the variation of the gap from 244 meV to 288 meV between the maximum and minimum charge density wave locations. Density functional theory calculations further confirm a vertical positioning of BEDT-TTF molecules in the molecular layer. While the observed charge density wave patterns are stable for the defect sites, they can be reversely switched for one molecular lattice site by means of inelastic tunnelling electron energy transfer with the electron energies exceeding 400 meV using a scanning tunneling microscope manipulation scheme.

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

Room temperature depinning of the charge-density waves in quasi-two-dimensional 1T-TaS 2 devices

Here we report on the depinning of nearly commensurate charge-density waves in 1T-TaS 2 thin films at room temperature. A combination of the differential current–voltage measurements with the low-frequency noise spectroscopy provides unambiguous means for detecting the depinning threshold field in quasi-2D materials. The depinning process in 1T-TaS 2 is not accompanied by an observable abrupt increase in electric current—in striking contrast to depinning in the conventional charge-density-wave materials with quasi-1D crystal structure. We explained it by the fact that the current density from the charge-density waves in the 1T-TaS 2 devices is orders of magnitude smaller than the current density of the free carriers available in the discommensuration network surrounding the commensurate charge-density wave islands. The depinning fields in 1T-TaS 2 thin-film devices are several orders of magnitude larger than those in quasi-1D van der Waals materials. Obtained results are important for the proposed applications of the charge-density wave devices in electronics.

2D materials↗
DOE OSTI2023

Density limits as disruption forecasters for spherical tokamaks

Fusion power output from spherical tokamaks would benefit from increased confined plasma density, but there exists a limit on the density before confinement is lost and the plasma current is disrupted. This density limit has long been characterized by a simple, global Greenwald limit proportional to the plasma current and inversely proportional to the cross sectional area of the plasma. It is shown that in the database of discharges from the National Spherical Tokamak Experiment (NSTX) and Mega Ampere Spherical Tokamak (MAST) spherical tokamaks, the likelihood of disruption does increase above the Greenwald limit, and especially in the plasma current rampdown phase. The physics of the density limit has been recently theoretically explored through local criteria. Several of these are tested using the disruption event characterization and forecasting (DECAF TM ) code for their potential effectiveness as disruption warning signals. For a limited set of NSTX discharges, a local island power balance criteria was found to be less reliable, presently, than the Greenwald limit. For this work, an empirical critical edge line density and a boundary turbulent transport limit were both tested for MAST-U, which has an electron density profile measurement with high spatial resolution in the outer part of the plasma. Both were found to have similar dependencies on key plasma parameters. In a limited set of MAST-U discharges that appear to disrupt due to rising density at values under the Greenwald limit, crossing of the boundary turbulent transport limit occurred close to the time of disruption. Finally, these limits were evaluated for their potential use in real-time, and it was found that with the necessary real-time inputs and with refinement through further testing, these limits could be implemented in a real-time disruption forecasting system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effect of resistivity on poloidal asymmetries in electron density in the edge and scrape off layer in 3D full-F gyro-fluid simulations with FELTOR

In this contribution, we analyze the poloidal asymmetries in electron density in the edge and scrape-off layer (SOL) in a COMPASS-sized, diverted tokamak with the 3D full-F, isothermal, electromagnetic, and gyro-fluid model FELTOR. The study is performed for different simulations that span over 2 orders of magnitude in resistivity. The poloidal asymmetries of density are evaluated in the edge and the SOL relative to the outer midplane (OMP), where the highest densities are usually found. In the closed magnetic surfaces, the relative poloidal asymmetry with respect to the OMP is not larger than 20% independently of the plasma resistivity. For the open field lines in the SOL, the relative density asymmetry can range from 55% for the highest resistivity to around 40% for the lowest. The lowest densities inside the separatrix are found between the inner midplane and the top of the magnetic configuration, away from the X-point. In the SOL it is usually close to the X-point in the high field side. The observations in the closed field lines are consistent with ballooning transport but in the SOL oppose experimental evidence for high-density plasmas. In conclusion, this indicates the necessity for more complex physics to reproduce the experimental observations in the SOL for higher density plasmas, such as neutral-plasma dynamics and realistic divertor conditions.

3D full-f gyrofluid↗
DOE OSTI2021

Enhanced particle transport events approaching the density limit of the J-TEXT tokamak

Enhanced particle transport events are discovered and analyzed as the density limit of the J-TEXT tokamak is approached. Edge shear layer collapse is observed and the ratio of Reynolds power to turbulence production decreases. Simultaneously, the divergence of turbulence internal energy flux (i.e. turbulence spreading) increases, indicating that shear layer collapse triggers an outward spreading event. Studies of correlations show that the enhanced particle transport events are quasi-coherent, and manifested primarily in density fluctuations which exhibit positive skewness. Furthermore, electron adiabaticity emerges as the critical parameter which signals transport event onset. For α < 0.35 as density approaches the Greenwald density, both turbulence spreading and density fluctuations rise rapidly. Taken together, these results elucidate the connections between edge shear layer, density fluctuations, particle transport events, turbulence spreading and plasma edge cooling as the density limit is approached.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation and analysis of lower-hybrid-current-drive density limit in EAST

Lower hybrid current drive (LHCD) experiments with line-averaged density up to ~5.1 × 10 19 m –3 were performed in EAST L-mode plasmas. When the line-averaged density rises above a critical value, the hard x-ray (HXR) emission falls to the noise level, indicating that the LHCD density limit is encountered. The experimental results show that the LHCD density limit can be increased with higher wave source frequency (ƒ 0 ) and higher magnetic field (B t ). Although a higher LHCD density limit is obtained by a higher magnetic field for both 2.45 GHz and 4.6 GHz waves, the results show a stronger dependence on the magnetic field for the 4.6 GHz case. Analysis suggests that, for normal operation with a relatively low magnetic field (1.6 T $\leqslant$ B t $\leqslant$ 2.5 T) on EAST, the dominant mechanisms responsible for the LHCD density limit are different between the 2.45 GHz and 4.6 GHz waves. The wave accessibility plays a more significant role during 4.6 GHz LHCD experiments, while parasitic losses due to parametric decay instability (PDIs) dominate the accessibility issue in the 2.45 GHz case. Collisional loss in the scrape-off layer (SOL) may explain the 4.6 GHz result when combined with the accessibility limit at high density and low temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The role of isotope mass on neutral fueling and density pedestal structure in the DIII-D tokamak

Abstract Experimental measurements on DIII-D of hydrogen neutral penetration lengths ( λ n 0 ) on the high field side (HFS) are longer by a factor of 2 than for deuterium consistent with the thermal velocity ratio for neutrals at the same temperature ( v t h H / v t h D = 2 ) . This ratio is constant for both low and high pedestal electron density. At low pedestal density ( n e ∼ 4 × 10 19 m − 3 ) , the neutral penetration length is greater than the density pedestal width for both isotopes, and the additional 40% increase of neutral penetration in hydrogen widens the pedestal by the same amount. As the density pedestal height increases ( n e ∼ 6 × 10 19 m − 3 ) , the neutral penetration lengths drop below the density pedestal widths for both isotopes, and the increased penetration of hydrogen has no increased effect on the pedestal width compared to deuterium. Extrapolating to future reactor-relevant high electron density pedestals, the isotope-mass change in neutral fueling on the HFS from the deepest neutral penetration of hydrogen, to the shortest neutral penetration of tritium will be negligible (0.2–0.4 cm) in comparison to estimates of the density pedestal width (6–8.5 cm).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗