Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Spectral Index and Quasi-Periodic Oscillation Frequency Correlation in Black Hole Sources: Observational Evidence of Two Phases and Phase Transition in Black Holes

Recent studies have shown that strong correlations are observed between the low frequencies (1-10 Hz) of quasi-periodic oscillations (QPOs) and the spectral power law index of several black hole (BH) candidate sources, in low (hard) states, steep power law (soft) states, and transitions between these states. The observations indicate that the X-ray spectra of such state (phases) show the presence of a power-law component and are sometimes related to simultaneous radio emission, indicating the probable presence of a jet. Strong QPOs (>20% rms) are present in the power density spectrum in the spectral range where the power-law component is dominant (i.e., 60%90%). This evidence contradicts the dominant, long-standing interpretation of QPOs as a signature of the thermal accretion disk. We present the data from the literature and our own data to illustrate the dominance of power-law index-QPO frequency correlations. We provide a model that identifies and explains the origin of the QPOs and how they are imprinted on the properties of the power-law flux component. We argue for the existence of a bounded compact coronal region that is a natural consequence of the adjustment of the Keplerian disk flow to the innermost sub-Keplerian boundary conditions near the central object and that ultimately leads to the formation of a transition layer (TL) between the adjustment radius and the innermost boundary. The model predicts two phases or states dictated by the photon upscattering produced in the TL: (1) a hard state, in which the TL is optically thin and very hot (kT approximately greater than 50 keV), producing photon upscattering via thermal Comptonization (the photon spectrum index Gamma approximates 1.7 for this state is dictated by gravitational energy release and Compton cooling in an optically thin shock near the adjustment radius), and (2) a soft state that is optically thick and relatively cold (kT approximately less than 5 keV the index for this state, Gamma approximates 2.8, is determined by soft-photon upscattering and photon trapping in a converging flow into the BH). In the TL model for the corona, the QPO frequency V(sub high) is related to the gravitational (close to Keplerian) frequency V(sub K) at the outer (adjustment) radius and v(sub low) is related to the TL's normal mode (magnetoacoustic) oscillation frequency v(sub MA) . The observed correlations between index and low and high QPO frequencies are readily explained in terms of this model. We also suggest a new method for evaluation of the BH mass using the index-frequency correlation.

Titarchuk, Lev↗
NASA NTRS2004

Spectral Index and Quasi-Periodic Oscillation Frequency Correlation in Black Hole (BH) Sources: Observational Evidence of Two Phases and Phase Transition in BHs

Recent studies have shown that strong correlations are observed between the low frequencies (1-10 Hz) of quasiperiodic oscillations (QPOs) and the spectral power law index of several Black Hole (BH) candidate sources, in low hard states, steep power-law (soft) states and in transition between these states. The observations indicate that the X-ray spectrum of such state (phases) show the presence of a power-law component and are sometimes related to simultaneous radio emission indicated the probable presence of a jet. Strong QPOs (less than 20% rms) are present in the power density spectrum in the spectral range where the power-law component is dominant ( i.e. 60-90% ). This evidence contradicts the dominant long standing interpretation of QPOs as a signature of the thermal accretion disk. We present the data from the literature and our own data to illustrate the dominance of power-law index-QPO frequency correlations. We provide a model, that identifies and explains the origin of the QPOs and how they are imprinted on the properties of power-law flux component. We argue the existence of a bounded compact coronal region which is a natural consequence of the adjustment of Keplerian disk flow to the innermost sub-Keplerian boundary conditions near the central object and that ultimately leads to the formation of a transition layer (TL) between the adjustment radius and the innermost boundary. The model predicts two phases or states dictated by the photon upscattering produced in the TL: (1) hard state, in which the TL is optically thin and very hot (kT approx. greater than 50 keV) producing photon upscattering via thermal Componization; the photon spectrum index Gamma appprox.1.5 for this state is dictated by gravitational energy release and Compton cooling in an optically thin shock near the adjustment radius; (2) a soft state which is optically thick and relatively cold (approx. less than 5 keV); the index for this state, Gamma approx. 2.8 is determined by soft-photon upscattering and photon trapping in converging flow into BH. In the TL model for corona the QPO frequency vnu(sub high) is related to the gravitational (close to Keplerian) frequency nu(sub K) at the outer (adjustment) radius and nu(sub low) is related to the TL s normal mode (magnetoacoustic) oscillation frequency nu(sub MA). The observed correlations between index and low and high QPO frequencies are readily explained in terms of this model. We also suggest a new method for evaluation of the BH mass using the index-frequency correlation.

Titarchuk, Lev↗
NASA NTRS1992

Near-resonant excitation and propagation of eccentric density waves by external forcing

An overview is presented of the astronomical evidence that relatively massive, distended, gaseous disks form as a natural by-product of the process of star formation, and also the numerical evidence that SLING-amplified eccentric modes in the outer parts of such disks can drive one-armed spiral density waves in the inner parts by near-resonant excitation and propagation. An ordinary differential equation (ODE) of the second order that approximately governs the nonlocalized forcing of waves in a disk satisfying Lindblad resonance almost everywhere is derived. When transformed and appended with an extra model term, this ODE implies, for free waves, the usual asymptotic results of the WKBJ dispersion relationship and the propagation Goldreich-Tremaine (1978) formula for the resonant torque exerted on a localized Lindblad resonance. An analytical solution is given for the rate of energy and angular momentum transfer by nonlocalized near-resonant forcing in the case when the disk has power-law dependences on the radius of the surface density and temperature.

Ostriker, Eve C.↗
NASA NTRS1976

Measurements of uranium mass confined in high density plasmas

An X-ray absorption method for measuring the amount of uranium confined in high density, rf-heated uranium plasmas is described. A comparison of measured absorption of 8 keV X-rays with absorption calculated using Beer Law indicated that the method could be used to measure uranium densities from 3 times 10 to the 16th power atoms/cu cm to 5 times 10 to the 18th power atoms/cu cm. Tests were conducted to measure the density of uranium in an rf-heated argon plasma with UF6 infection and with the power to maintain the discharge supplied by a 1.2 MW rf induction heater facility. The uranium density was measured as the flow rate through the test chamber was varied. A maximum uranium density of 3.85 times 10 to the 17th power atoms/cu cm was measured.

Stoeffler, R. C.↗
NASA NTRS2006

Electron Velocity Distribution Function in Magnetic Clouds in the Solar Wind

We present a study of the kinetic properties of the electron velocity distribution functions within magnetic clouds, since they are the dominant thermal component. The study is based on high time resolution data from the GSFC WIND/SWE electron spectrometer and the Berkeley 3DP electron plasma instruments. Recent studies on magnetic clouds have shown observational evidence of anti-correlation between the total electron density and electron temperature, which suggest a polytrope law P(sub e) = alpha(Nu(sub e) (sup gamma)) for electrons with the constant gamma approximates 0.5 < 1. This anti-correlation and small polytropic gamma-values is interpreted in the context of the presence of highly non-Maxwellian electron distributions (i.e. non-thermal) within magnetic clouds. These works suggested that the non-thermal electrons can contribute as much as 50% of the total electron pressure within magnetic clouds. We have revisited some of the magnetic cloud events previously studied and attempted to quantify the nature of the non-thermal electrons by modeling the electron velocity distribution function using a kappa distribution function to characterize the kinetic non-thermal effects. If non-thermal tail effects are the source for the anti-correlation between the moment electron temperature and density and if the kappa distribution is a reasonable representative model of non-thermal effects, then the electron velocity distribution within magnetic clouds should show indication for small K-values when gamma < 1.

Nieves-Chinchil, Teresa↗
DOE OSTI2021

Robust d-Wave Superconductivity in the Square-Lattice t–J Model

Unravelling competing orders emergent in doped Mott insulators and their interplay with unconventional superconductivity is one of the major challenges in condensed matter physics. Here, to explore the possible superconducting state in a doped Mott insulator, we study the square-lattice t-J model with both the nearest-neighbor and next-nearest-neighbor electron hoppings and spin interactions. By using the state-of-the-art density matrix renormalization group calculation with imposing charge U(1) and spin SU(2) symmetries on the six-leg cylinders, we establish a quantum phase diagram including three phases: a stripe charge density wave phase, a superconducting phase without static charge order, and a superconducting phase coexistent with a weak charge stripe order. Crucially, we demonstrate that the superconducting phase has a power-law pairing correlation that decays much slower than the charge density and spin correlations, which is a quasi-1D descendant of the uniform d-wave superconductor in two dimensions. These findings reveal that enhanced charge and spin fluctuations with optimal doping is able to produce robust d-wave superconductivity in doped Mott insulators, providing a foundation for connecting theories of superconductivity to models of strongly correlated systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Plasma irregularities in the comet's tail

Scintillation theory is invoked to explain fluctuations in radio intensity observed during occultation of the extragalactic radio source PKS 2025-15 by the plasma tail of comet 1973 XII on Jan. 5, 1975. Plasma irregularities and turbulence in the tail of the comet (Kohoutek 1973f) are fitted to a Gaussian spectrum and to a Kolmogorov power-law spectrum in analyzing the scintillation data. The rms fluctuation of electron density in the cometary tail is reported at 80 electrons per cu mm, the inner scale of the fluctuation at 800 km, and the largest scale of fluctuation at possibly 400,000 km. A hump in the comet power-law spectrum is noted. Use of the power spectrum of electron density fluctuations to predict the power spectrum of magnetic field fluctuations for irregularities associated with hydromagnetic turbulence is recommended.

Lee, L. C.↗
NASA NTRS2021

The Field Substellar Mass Function Based on the Full-sky 20 pc Census of 525 L, T, and Y Dwarfs

We present final Spitzer trigonometric parallaxes for 361 L, T, and Y dwarfs. We combine these with prior studies to build a list of 525 known L, T, and Y dwarfs within 20 pc of the Sun, 38 of which are presented here for the first time. Using published photometry and spectroscopy as well as our own follow-up, we present an array of color–magnitude and color–color diagrams to further characterize census members, and we provide polynomial fits to the bulk trends. Using these characterizations, we assign each object a T(eff) value and judge sample completeness over bins of T(eff) and spectral type. Except for types ≥T8 and T(eff) < 600 K, our census is statistically complete to the 20 pc limit. We compare our measured space densities to simulated density distributions and find that the best fit is a power law (dN/dM ∝ M^(-α) with α = 0.6 ± 0.1. We find that the evolutionary models of Saumon & Marley correctly predict the observed magnitude of the space density spike seen at 1200 K < T(eff) < 1350 K, believed to be caused by an increase in the cooling timescale across the L/T transition. Defining the low-mass terminus using this sample requires a more statistically robust and complete sample of dwarfs ≥Y0.5 and with T(eff) < 400 K. We conclude that such frigid objects must exist in substantial numbers, despite the fact that few have so far been identified, and we discuss possible reasons why they have largely eluded detection.

J. Davy Kirkpatrick↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
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 OSTI2019

Features of near and far scrape-off layer heat fluxes on the Wendelstein 7-X inboard limiters

The power fall-off length is used as a characteristic dimension to describe the power exhaust in the scrape-off layer (SOL) in magnetic confinement fusion devices. Measurements from the limiter campaign of the stellarator Wendelstein 7-X with inboard limiters are presented with the first-time characterization of the limiter heat loads. Two fall-off regimes are found with a narrow near SOL with power fall-off length in millimeters and a wider far SOL with fall-off lengths of several centimeters. An attempt is made to describe both regimes with a scaling law for plasma with different heating powers and densities. Finally, the results confirm the major geometry effects of the connection length on the heat transport predicted by 3D modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗