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 1,063 records · Page 59

Analysis of the density field around a supersonic conical projectile using quantitative schlieren

Quantitative schlieren imaging is a flow measurement technique that is capable of measuring density fields throughout refractive flowfields. The technique was applied here to measure the density field surrounding supersonic conical projectiles in free flight. Shock waves attached to a supersonic conical projectile offer a simple geometry with the well-established Taylor–Maccoll analytical flow solution to which these experimental measurements were compared. The schlieren images recorded a projection of the index of refraction field surrounding the 10° half-angle cones which was converted to density first through an Abel inversion and then the Gladstone–Dale law. Three Abel inversion methods—two-point, three-point, and arbitrary ray axisymmetric projection (ARAP)—were applied to deconvolute the three-dimensional flow within the constrained axisymmetric flow field. The resulting reconstructed density profiles were compared to the Taylor–Maccoll solution, parameterized by cone geometry and Mach number. The experimental density fields demonstrated strong agreement with the theoretical profiles. Experimental consistency was confirmed across various projectile speeds, demonstrating quantitative schlieren’s capability to accurately reconstruct the density of the flow field, even within the resolution constraints imposed by high-speed imaging. In conclusion, an assessment of experimental uncertainties in the density reconstruction was performed.

Abel Inversion↗

High-field side scrape-off layer density profile measurements and implications for high-field side LHCD coupling in DIII-D

The high-field side (HFS) scrape-off layer (SOL) is an often under-diagnosed region of tokamak plasmas. Situated in a region with favorable curvature, the HFS SOL has minimal turbulence-induced radial transport compared to the low-field side. Using profile reflectometry, the HFS SOL density profile is measured with high temporal resolution for a wide range of DIII-D plasma discharges for the first time. The magnetic configuration, particularly the location of the secondary separatrix (SS) largely determines the HFS SOL density profile. Additionally, density perturbations induced by edge-localized modes (ELMs) are observed in the HFS SOL, and the impact of ELMs on the SOL density profile is characterized for a range of magnetic configurations. It is found that ELM-induced HFS SOL perturbations are generally localized to within the SS and ELM-induced changes to the HFS wall density are minimized in near double-null configurations. The characterization of the HFS SOL density profile on DIII-D is crucial for the effective coupling of the planned HFS lower hybrid current drive (LHCD) launcher. Experimental HFS SOL measurements are used as input to full-wave simulations of LHCD coupling. Furthermore, this work simulates the effect of the magnetic configuration on LHCD coupling. It also simulates LHCD coupling in ELMs and high q min discharges. During such discharges, the LHCD coupling is predicted to be resilient to ELM-induced SOL density perturbations.

high-field side↗

Modification of SOL profiles and fluctuations with line-average density and divertor flux expansion in TCV

A set of Ohmic density ramp experiments addressing the role of parallel connection length in modifying scrape off layer (SOL) properties has been performed on the TCV tokamak. The parallel connection length has been modified by varying the poloidal flux expansion f x . It will be shown that this modification does not influence neither the detachment density threshold, nor the development of a flat SOL density profile which instead depends strongly on the increase of the core line average density. The modification of the SOL upstream profile, with the appearance of what is generally called a density shoulder, has been related to the properties of filamentary blobs. Blob size increases with density, without any dependence on the parallel connection length both in the near and far SOL. The increase of the density decay length, corresponding to a profile flattening, has been related to the variation of the divertor normalized collisionality Λ div , showing that in TCV the increase of Λ div is not sufficient to guarantee the SOL upstream profile flattening.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental investigation of PDI bifurcation of lower hybrid waves during electron density ramp-up in EAST

Abstract The effect of parametric decay instability (PDI) on the current drive efficiency of 4.6 GHz lower hybrid (LH) waves in EAST is investigated experimentally, showing the PDI channel bifurcation of LH waves for the first time in EAST. First, experiments with three platforms of LH power were performed, achieving the LH power required for the PDI occurrence. Second, PDI bifurcation experiments were further carried out by ramping up the plasma electron density. The loop voltage increases with an increase in density, implying a decrease in the current driven by the LH wave. PDI bifurcation during electron density ramp-up was studied by analyzing the parallel refractive index ( n ∥ ) and the frequency spectrum broadening, which is measured with a radio frequency magnetic probe array recently installed close to the LH antenna. It is observed for the first time that they both first increase with density, then there is not much variation and a clear sideband in the frequency spectrum is also observed when the density is up to 4 × 10 19 m − 3 , suggesting a change in the PDI channel. Calculation of the mode growth rate driven by PDI shows that when the edge electron density is up to 1.9 × 10 18 m − 3 , the growth rate of the ion cyclotron quasi-mode (ICQM) will exceed that of the ion sound quasi-mode (ISQM), quantitively explaining that with an increase in density, the PDI channel partly transits from the ISQM to the ICQM channel. Studies provide a possible way to reduce the power deposition in the edge region and improve drive capability by means of mitigating PDI behavior.

Physics↗

Impact of ionization and transport on pedestal density structure in DIII-D and Alcator C-Mod

Abstract This paper investigates the role of ionization on the pedestal structure using both measurements and modeling for H-mode plasmas on DIII-D and Alcator C-Mod to enhance our ability to predict pedestal behavior in future pilot plants. The impact of the neutral penetration depth on the pedestal density is investigated using dimensionally matching hydrogen and deuterium DIII-D H-mode discharges at low and high electron density. The DIII-D Lyman- α diagnostic measurements show that hydrogen neutrals penetrate deeper inside the plasma on both the high field and low field side, while the pedestal electron density structure is similar for both isotopes. However, as the opaqueness increases we observe that the pedestal density gradient becomes stiff, similar to prior observations on DIII-D and C-Mod (Mordijck 2020 Nuclear Fusion 60 082006). In addition, these results also confirm prior measured and modeled poloidal asymmetries in neutral densities, indicating that to make transport predictions, 2D neutral modeling is necessary. The first direct validation of SOLPS-ITER for the measured brightness, emissivity and neutral densities for three different confinement regimes on C-Mod is introduced. The SOLPS-ITER model shows good agreement, within the constrains of the model for all regimes. In addition, a comparison of SOLPS-ITER modeling for DIII-D and C-Mod shows that as opaqueness increases, the role of divertor fueling and thus poloidal asymmetries in the neutral density profiles decreases. Based on these experimental and modeling results we estimate the size of a potential particle pinch using typical values for the diffusion coefficient for both DIII-D and C-Mod H-mode discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-linear dependence of ion heat flux on plasma density at the L–H transition of JET NBI-heated deuterium–tritium plasmas

Abstract Recent JET D–T campaigns opened the possibility of unique isotope studies to investigate the L–H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L–H transitions in the low and high-density branches of JET NBI-heated D–T plasmas. As discussed in the paper, L–H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux ( Q i ). We present the first power balance analysis of JET NBI-heated D–T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D–T Q i , with a lower Q i in D–T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L–H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Q i deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Q i and the power threshold minimum at JET. At low plasma density, NBI power dominates Q i , while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D–T NBI-heated plasmas at JET.

Vincenzi, P. (ORCID:0000000251564354)↗

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗

CryoTEN: efficiently enhancing cryo-EM density maps using transformers

Abstract Motivation Cryogenic electron microscopy (cryo-EM) is a core experimental technique used to determine the structure of macromolecules such as proteins. However, the effectiveness of cryo-EM is often hindered by the noise and missing density values in cryo-EM density maps caused by experimental conditions such as low contrast and conformational heterogeneity. Although various global and local map-sharpening techniques are widely employed to improve cryo-EM density maps, it is still challenging to efficiently improve their quality for building better protein structures from them. Results In this study, we introduce CryoTEN—a 3D UNETR++ style transformer to improve cryo-EM maps effectively. CryoTEN is trained using a diverse set of 1295 cryo-EM maps as inputs and their corresponding simulated maps generated from known protein structures as targets. An independent test set containing 150 maps is used to evaluate CryoTEN, and the results demonstrate that it can robustly enhance the quality of cryo-EM density maps. In addition, automatic de novo protein structure modeling shows that protein structures built from the density maps processed by CryoTEN have substantially better quality than those built from the original maps. Compared to the existing state-of-the-art deep learning methods for enhancing cryo-EM density maps, CryoTEN ranks second in improving the quality of density maps, while running >10 times faster and requiring much less GPU memory than them. Availability and implementation The source code and data are freely available at https://github.com/jianlin-cheng/cryoten.

Biochemistry & Molecular Biology↗

X-ray Studies of the Transformation from High-to Low-density Amorphous Water

Here we report about the structural evolution during the conversion from high-density amorphous ices at ambient pressure to the low-density state. Using high-energy X-ray diffraction we have monitored the transformation by following in reciprocal space the structure factor SOO(Q) and derive in real space the pair distribution function gOO(r). Heating equilibrated high-density amorphous ice (eHDA) at a fast rate (4 K/min), the transition to the low-density form occurs very rapidly, while domains of both high- and low-density coexist. On the other hand, the transition in the case of unannealed HDA (uHDA) and very-high-density amorphous ice (VHDA) is more complex and of continuous nature. The direct comparison of eHDA and uHDA indicates that the molecular structure of uHDA contains a larger amount of tetrahedral motives. The different crystallization behavior of the derived low-density amorphous (LDA) states is interpreted as emanating from an increased tetrahedral coordination present in uHDA.

amorphous ice↗

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗

Phosphorus allotropes: Stability of black versus red phosphorus re-examined by means of the van der Waals inclusive density functional method

We evaluate the energetic stabilities of white, red, and black allotropes of phosphorus using density functional theory (DFT) and hybrid functional methods, van der Waals (vdW) corrections (DFT+vdW and hybrid+vdW), vdW density functionals, and random phase approximation (RPA). In this work, we find that stability of black phosphorus over red-V (i.e., the violet form) is not ubiquitous among these methods, and the calculated enthalpies for the reaction phosphorus (red-V)→ phosphorus (black) are scattered between –20 and 40 meV/atom. With local density and generalized gradient approximations, and hybrid functionals, mean absolute errors (MAEs) in densities of P allotropes relative to experiments are found to be around 10%–25%, whereas with vdW-inclusive methods, MAEs in densities drop below ~5%. While the inconsistency among the density functional methods could not shed light on the stability puzzle of black versus red phosphorus, comparison of their accuracy in predicting densities and the supplementary RPA results on relative stabilities indicate that opposite to the common belief, black and red phosphorus are almost degenerate, or the red-V (violet) form of phosphorus might even be the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vortices in low-density neutron matter and cold Fermi gases

Cold gas experiments can be tuned to achieve strongly-interacting regimes such as that of low-density neutron matter found in neutron-stars' crusts. We report T=0 diffusion Monte Carlo results (i) for the ground state of both spin-1/2 fermions with short-range interactions and low-density neutron matter in a cylindrical container, and (ii) properties of these systems with a vortex line excitation. We calculate the equation of state for cold atoms and low-density neutron matter in the bulk systems, and we contrast it to our results in the cylindrical container. We compute the vortex line excitation energy for different interaction strengths, and we find agreement between cold gases and neutron matter for very low densities. We also calculate density profiles, which allow us to determine the density depletion at the vortex core, which depends strongly on the short-ranged interaction in cold atomic gases, but it is of ≈25% for neutron matter in the density regimes studied in this work. Our results can be used to constrain neutron matter properties by using measurements from cold Fermi gases experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Auxiliary function approach for determining symmetry energy at suprasaturation densities

Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear energy density optimization: Shell structure

Nuclear density functional theory is the only microscopical theory that can be applied throughout the entire nuclear landscape. Its key ingredient is the energy density functional. Here, we propose a new parametrization UNEDF 2 of the Skyrme energy density functional. The functional optimization is carried out using the POUNDERS optimization algorithm within the framework of the Skyrme Hartree-Fock-Bogoliubov theory. Compared to the previous parametrization UNEDF 1, restrictions on the tensor term of the energy density have been lifted, yielding a very general form of the energy density functional up to second order in derivatives of the one-body density matrix. In order to impose constraints on all the parameters of the functional, selected data on single-particle splittings in spherical doubly-magic nuclei have been included into the experimental dataset. The agreement with both bulk and spectroscopic nuclear properties achieved by the resulting UNEDF 2 parametrization is comparable with UNEDF 1. While there is a small improvement on single-particle spectra and binding energies of closed shell nuclei, the reproduction of fission barriers and fission isomer excitation energies has degraded. As compared to previous UNEDF parametrizations, the parameter confidence interval for UNEDF 2 is narrower. In particular, our results overlap well with those obtained in previous systematic studies of the spin-orbit and tensor terms. UNEDF 2 can be viewed as an all-around Skyrme EDF that performs reasonably well for both global nuclear properties and shell structure. However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally. These results suggest that the standard Skyrme energy density has reached its limits, and significant changes to the form of the functional are needed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Higher-order nonequilibrium term: Effective power density quantifying evolution towards or away from local thermodynamic equilibrium

A common approach to assess the nature of energy conversion in a classical fluid or plasma is to compare power densities of the various possible energy conversion mechanisms. A leading research area is quantifying energy conversion for systems that are not in local thermodynamic equilibrium (LTE), as is common in a number of fluid and plasma systems. Here we introduce the “higher-order nonequilibrium term” (HORNET) effective power density, which quantifies the rate of change of departure of a phase space density from LTE. It has dimensions of power density, which allows for quantitative comparisons with standard power densities. We employ particle-in-cell simulations to calculate HORNET during two processes, magnetic reconnection and decaying kinetic turbulence in collisionless magnetized plasmas, that inherently produce non-LTE effects. We investigate the spatial variation of HORNET and the time evolution of its spatial average. By comparing HORNET with power densities describing changes to the internal energy (pressure dilatation, Pi − D , and divergence of the vector heat flux density), we find that HORNET can be a significant fraction of these other measures (8% and 35% for electrons and ions, respectively, for reconnection; up to 67% for both electrons and ions for turbulence), meaning evolution of the system towards or away from LTE can be dynamically important. Applications to numerous plasma phenomena are discussed. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles equation-of-state table of silicon and its effects on high-energy-density plasma simulations

Using density-functional theory–based molecular-dynamics simulations, we have investigated the equation of state for silicon in a wide range of plasma density and temperature conditions of ρ=0.001–500g/cm 3 and T=2000–10 8 K. With these calculations, we have established a first-principles equation-of-state (FPEOS) table of silicon for high-energy-density (HED) plasma simulations. When compared with the widely used SESAME-EOS model (Table 3810), we find that the FPEOS-predicted Hugoniot is ~20% softer; for off-Hugoniot plasma conditions, the pressure and internal energy in FPEOS are lower than those of SESAME EOS for temperatures above T ≈ 1–10 eV (depending on density), while the former becomes higher in the low- T regime. The pressure difference between FPEOS and SESAME 3810 can reach to ~50%, especially in the warm-dense-matter regime. Implementing the FPEOS table of silicon into our hydrocodes, we have studied its effects on Si-target implosions. When compared with the one-dimensional radiation-hydrodynamics simulation using the SESAME 3810 EOS model, the FPEOS simulation showed that (1) the shock speed in silicon is ~10% slower; (2) the peak density of an in-flight Si shell during implosion is ~20% higher than the SESAME 3810 simulation; (3) the maximum density reached in the FPEOS simulation is ~40% higher at the peak compression; and (4) the final areal density and neutron yield are, respectively, ~30% and ~70% higher predicted by FPEOS versus the traditional simulation using SESAME 3810. All of these features can be attributed to the larger compressibility of silicon predicted by FPEOS. Furthermore, these results indicate that an accurate EOS table, like the FPEOS presented here, could be essential for the precise design of targets for HED experiments.

36 MATERIALS SCIENCE↗

Accuracy of ab initio electron correlation and electron densities in vanadium dioxide

Here, diffusion quantum Monte Carlo results are used as a reference to analyze properties related to phase stability and magnetism in vanadium dioxide computed with various formulations of density functional theory. We introduce metrics related to energetics, electron densities and spin densities that give us insight on both local and global variations in the antiferromagnetic M1 and R phases. Importantly, these metrics can address contributions arising from the challenging description of the 3d orbital physics in this material. We observe that the best description of energetics between the structural phases does not correspond to the best accuracy in the charge density, which is consistent with observations made recently by Medvedev et~al. in the context of isolated atoms. However, we do find evidence that an accurate spin density connects to correct energetic ordering of different magnetic states in VO$_2$, although local, semilocal, and meta-GGA functionals tend to erroneously favor demagnetization of the vanadium sites. The recently developed SCAN functional stands out as remaining nearly balanced in terms of magnetization across the M1-R transition and correctly predicting the ground state crystal structure. In addition to ranking current density functionals, our reference energies and densities serve as important benchmarks for future functional development.

36 MATERIALS SCIENCE↗

Nanoscale mapping of the electron density at Al grain boundaries and correlation with grain-boundary energy

Grain-boundary (GB) structures and energies are often calculated and have revealed correlations between the GB energy and change in electron density at the GB. In this work, the plasmon peak in valence electron energy-loss spectroscopy (VEELS) was used to determine the variation in electron density across four well-characterized GBs in Al, spanning a range of known GB energies. The results show that the plasmon energy is lower at the GB than in the adjacent grains due to a decrease in electron density, and the GB energy increases proportional to the density decrease. The decrease in electron density also extends further into adjacent grains with increasing GB energy, extending beyond the geometric changes, or physical width, revealed by electron microscopy. Plasmon damping also increases with increasing GB energy, indicative of increasing disruption of the electron density with increasing GB energy. These results demonstrate that VEELS can be a valuable tool for detecting small electron density changes at GBs, and this change clearly influences, and is correlated to, the GB energy.

36 MATERIALS SCIENCE↗