Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Laws”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Transfer and Flushing Evaluation for Potential Precipitated Solids in the 241-AP-106 to EMF Transfer Pipeline for DFLAW

Waste transfer pipeline flushing is necessary to maintain pipeline performance, but added flush fluid increases the nuclear waste inventory. During the Direct Feed Low-Activity Waste (DFLAW) mission, waste transfer lines from the Interim Low-Activity Waste Storage Tank (i.e., 241-AP-106) to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility must have the capability of being flushed to prevent accumulation of solids and to mitigate corrosion concerns. Operational experience at Hanford has shown that the solids precipitation occurs during processing operations from dilute Hanford liquids. Pacific Northwest National Laboratory was requested to evaluate the potential precipitated solid particle transport and flushing operational capabilities of the AP-106 through the Effluent Management Facility (EMF) low-point drain LAW feed pipeline. The evaluation results do not provide operational requirements, but rather provide a scoping basis for understanding the potential operational significance of solids precipitation in the pipeline. The salt and aluminum phase solids most likely to precipitate from the LAW feeds during cooling, evaporation, or mixing have particle densities that are estimated to range from 1.62 to 2.78 g/mL and a spherical particle size range of 8 to 2100 μm, depending on the solid phase considered. Application of the method to calculate the critical deposition velocity required by TFC-ENG-STD-26 at conservatively bounding estimates for the solids concentrations demonstrates that potentially precipitated solids may deposit on the bottom of the transfer pipe invert at the lower LAW feed pipeline flow rate, but the upper LAW feed pipeline flow rate will likely prevent deposition. The system pressure limit for the AP-106 through the EMF low-point drain LAW feed pipeline likely exceeds any pressure loss at a calculated critical deposition velocity even for the most adverse potential precipitated solid. However, should solids settle as a result of no-flow conditions, the flow capabilities of the AP-106 through the EMF low-point drain LAW feed pipeline are shown, based on the available literature, to be potentially inadequate for effective flushing operations to remove solids, even at solids concentrations below the maximum specified for LAW feed. Thus, stepwise solids accumulations over multiple transfers may be an issue, and eventual line plugging may occur. Further inadequacies may potentially be realized should the cohesive nature of certain precipitates be accounted for.

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

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations

We propose an unsupervised machine-learning checkpoint-restart (CR) algorithm for particle-in-cell (PIC) algorithms using Gaussian mixtures (GM). The algorithm compresses the particle population per spatial cell by constructing a velocity distribution function using GM. Particles are reconstructed at restart time by local resampling of the Gaussians. To guarantee fidelity of the CR process, we ensure the exact preservation of invariants such as charge, momentum, and energy for both compression and reconstruction stages, everywhere on the mesh. We also ensure the preservation of Gauss' law after particle reconstruction by exactly matching the density profile at restart time. As a result, the GM CR algorithm is shown to provide a clean, conservative restart capability while potentially affording orders of magnitude savings in input/output requirements. Here, we demonstrate the algorithm using a recently developed exactly energy- and charge-conserving PIC algorithm using both electrostatic and electromagnetic tests. The tests demonstrate not only a high-fidelity CR capability, but also its potential for enhancing the fidelity of the PIC solution for a given particle resolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Comparison of Linear Solvers for Resolving Flow in Three-Dimensional Discrete Fracture Networks

We compare various methods for resolving steady flow within three-dimensional discrete fracture networks, including direct methods, Krylov subspace methods with and without preconditioning, and multi-grid methods. We compared the performance of the methods based on compute times and scaling of the solution as a function of the number of grid nodes and log-variance of the hydraulic aperture. The methods are applied to three test cases: (a) variable density of networks with a truncated power-law distribution of fracture lengths, (b) a fixed network composed of monodisperse fracture sizes but varied permeability/aperture heterogeneity, (c) and a network based on field site in Nevada, US. We chose these cases to allow us to study the impact of the mesh size and flow properties, as well as to demonstrate our conclusions on a large-scale, realistic problem (more than 40 million mesh nodes). A direct solution using Cholesky factorization outperformed other methods for every example but was closely followed in performance by some algebraic multigrid (AMG) preconditioned Krylov subspace methods. Among the Krylov methods, conjugate gradients (CG) with an AMG preconditioner performs the best. Generally, Cholesky factorization is recommended, but CG with an AMG preconditioner may be suitable for very large problems beyond 40 million nodes where the entire linear system cannot reside in memory.

58 GEOSCIENCES↗
DOE OSTI2026

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗
DOE OSTI2022

On the origin of multihole oxygen evolution in haematite photoanodes

The oxygen evolution reaction (OER) plays a crucial role in (photo)electrochemical devices that use renewable energy to produce synthetic fuels. Recent measurements on semiconducting oxides have found a power law dependence of the OER rate on surface hole density, suggesting a multihole mechanism. In this study, using transient photocurrent measurements, density functional theory simulations and microkinetic modelling, we have uncovered the origin of this behaviour in haematite. We show here that the OER rate has a third-order dependence on the surface hole density. We propose a mechanism wherein the reaction proceeds by accumulating oxidizing equivalents through a sequence of one-electron oxidations of surface hydroxy groups. The key O–O bond formation step occurs by the dissociative chemisorption of a hydroxide ion involving three oxyl sites. At variance with the case of metallic oxides, the activation energy of this step is weakly dependent on the surface hole coverage, leading to the observed power law.

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

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Higgs-portal dark matter in Brane-world cosmology

Abstract The Higgs-portal scalar dark matter (DM) model is a simple extension of the Standard Model (SM) to incorporate a DM particle to the SM, where a $$Z_2$$ Z 2 -odd real scalar field is introduced as a DM candidate. We consider this DM model in the context of 5-dimensional brane-world cosmology, where our 3-dimensional space is realized as a hyper-surface embedded in 4-dimensional space. In the setup, all the SM and DM fields reside on the hyper-surface while graviton lives in the bulk. We consider two well-known brane-world cosmologies, namely, the Randall–Sundrum (RS) and the Gauss–Bonnet (GB) brane-world cosmologies, in which the standard Big Bang cosmology is reproduced at low temperatures below the so-called “transition temperature” while at high temperatures the expansion law of the universe is significantly modified. Such a non-standard expansion law directly impacts the prediction for the relic density of the Higgs-portal DM. We investigate the brane-world cosmological effects and identify the allowed model parameter region by combining the constraints from the observed DM relic density, and the direct and indirect DM detection experiments. It is well-known that only DM masses in the vicinity of half the Higgs boson mass are allowed in the Higgs-portal scalar DM model. We find that the allowed parameter region becomes more severely constrained and even disappears in the RS cosmology, while the GB cosmological effect significantly enlarges the allowed region. Upon discovering Higgs-portal DM, we can determine transition temperature in the GB brane-world cosmology.

Liu, Taoli↗
DOE OSTI2024

DualSPHysics-INL

Funded by the DOE's Bioenergy Technology Office through the Feedstock-Conversion Interface Consortium, INL researchers developed this code to model biomass freestock flow in various handling equipment, such as hoppers and augers. Built on top of an existing open-source code DualSPHysics (https://dual.sphysics.org/) , the enrichment includes modification of the mass conservation equation that switching tracking density to void ratio, adding a hypoplastic constitutive law to better capture the flow physics of this type of material and implementing a novel boundary condition that can handle the dynamic contact between material and equipment. All implementation were realized via Nvidia CUDA, so GPU accelaration can be leveraged to signifiantly speed up the computational process.

Jin, Wencheng [Idaho National Laboratory (INL), Id↗
DOE OSTI2024

Chemical-specific Parameters Dataset

The chemical-specific parameters dataset is searchable for physicochemical information for multiple chemicals simultaneously. After selecting chemicals of interest and the desired parameters, the RAIS will generate a table containing the values, chosen according to an established hierarchy. Results can be downloaded in Excel format. Over 40 parameters are available, including melting point, boiling point, density, density, vapor pressure, water solubility, and Henry’s Law constants. Thirteen primary sources are used to populate the dataset of chemical-specific parameters. These values should be used in cancer risk and noncancer hazard assessments for the calculation of preliminary remediation goals (PRGs), hazard characterization, and transport modeling. Users can select up to 1000 chemicals per query. The dataset supports environmental risk assessments, regulatory decision-making, and environmental planning with tools for benchmarking against risk-based standards. This structured approach ensures a robust evaluation of environmental risks tailored to regulatory needs.

Dolislager, Fred [Oak Ridge National Laboratory (O↗
DOE OSTI2025

Ferromagnets, a new anomaly, instantons, and (noninvertible) continuous translations

We discuss a large class of classical field theories with continuous translation symmetry. In the quantum theory, a new anomaly explicitly breaks this translation symmetry to a discrete symmetry. Furthermore, this discrete translation symmetry is extended by a d – 2-form global symmetry. All these theories can be described as U(1) gauge theories where Gauss law states that the system has nonzero charge density. Special cases of such systems can be phrased as theories with a compact phase space. Examples are ferromagnets and lattices in the lowest Landau level. In some cases, the broken continuous translation symmetry can be resurrected as a noninvertible symmetry. We clarify the relation between the discrete translation symmetry of the continuum theory and the discrete translation symmetry of an underlying lattice model. Our treatment unifies, clarifies, and extends earlier works on the same subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Many-body quantum dynamics slows down at low density

We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

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↗
DOE OSTI2020

Extracting a mixing parameter from 2D radiographic imaging of variable-density turbulent flow

We extract a suitably averaged fluctuating density from the two-dimensional radiographic image of a flow. The X-ray attenuation is given by the Beer–Lambert law which exponentially damps the incident beam intensity by a factor proportional to the density, opacity and thickness of the target. By making reasonable assumptions for the mean density, opacity and effective thickness of the target flow, we estimate the density fluctuation contribution to the attenuation. The extracted density fluctuations averaged across the thickness of the flow in the direction of the beam may be used to form the density–specific-volume correlation $b$. In a statistical description of variable-density turbulence, $b$ quantifies the degree of mixedness. The ability to extract a measure of mixedness from experimental data would be a powerful tool that could be used in the validation of mix models. The scheme proposed is tested for DNS data computed for variable density buoyancy-driven mixing. We quantify the deficits in the extracted value of $b$ due to target thickness, Atwood number and modeled signal noise. This analysis justifies using the proposed scheme to infer the mix parameter from thin targets at moderate to low Atwood numbers. Furthermore, to illustrate how the scheme might be used in a practical problem, we demonstrate its application to a radiographic image of counter-shear flow obtained from experiments at the National Ignition Facility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Diagonal entropy and topological phase transitions in extended Kitaev chains

We investigate the diagonal entropy for ground states of the extended Kitaev chains with extensive pairing and hopping terms. The systems contain rich topological phases equivalently represented by topological invariant winding numbers and Majorana zero modes. Both the finite size scaling law and block scaling law of the diagonal entropy are studied, which indicates that the diagonal entropy demonstrates volume effect. The parameter of volume term is regarded as the diagonal entropy density, which can identify the critical points of symmetry-protected topological phase transitions efficiently in the studied models, even for those with higher winding numbers. The formulation of block scaling law and the capability of diagonal entropy density in detecting topological phase transitions are independent of the chosen bases. In order to manifest the advantage of diagonal entropy, we also calculate the global entanglement, which cannot show clear signatures of the topological phase transitions. This work provides a new quantum-informatic approach to characterize the feature of the topologically ordered states and may motivate a deep understanding of the quantum coherence and diagonal entropy in various condensed matter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Liquid-like dynamics in a solid-state lithium electrolyte

Superionic materials represent a regime intermediate between the crystalline and liquid states of matter. Despite the considerable interest in potential applications for solid-state batteries or thermoelectric devices, it remains unclear whether the fast ionic diffusion observed in superionic materials reflects liquid-like dynamics or whether the hops of mobile ions are inherently coupled to more conventional lattice phonons. Here we reveal a crossover from crystalline vibrations to relaxational dynamics of ionic diffusion in the superionic compound Li 6 PS 5 Cl, a candidate solid-state electrolyte. By combining inelastic and quasi-elastic neutron-scattering measurements with first-principles-based machine-learned molecular dynamics simulations, we found that the vibrational density of states in the superionic state strongly deviates from the quadratic behaviour expected from the Debye law of lattice dynamics. The superionic dynamics emerges from overdamped phonon quasiparticles to give rise to a linear density of states characteristic of instantaneous normal modes in the liquid state. Further, we showed that the coupling of lattice phonons with a dynamic breathing of the Li+ diffusion bottleneck enables an order-of-magnitude increase in diffusivity. Thus, our results shed insights into superionics for future energy storage and conversion technologies.

25 ENERGY STORAGE↗