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

Aspects of gravitational portals and freeze-in during reheating

We conduct a systematic investigation of freeze-in during reheating while taking care to include both direct and indirect production of dark matter (DM) via gravitational portals and inflaton decay. Direct production of DM can occur via gravitational scattering of the inflaton, while indirect production occurs through scattering in the Standard Model radiation bath. We consider two main contributions to the radiation bath during reheating. The first, which may dominate at the onset of the reheating process, is produced via gravitational scattering of the inflaton. The second (and more standard contribution) comes from inflaton decay. We consider a broad class of DM production rates parametrized as R χ ∝ T n + 6 / Λ n + 2 , and inflaton potentials with a power-law form V ( ϕ ) ∝ ϕ k about the minimum. We find the relic density produced by freeze-in for each contribution to the Standard Model bath for arbitrary k and n , and compare these with the DM density produced gravitationally by inflaton scattering. We find that freeze-in production from the gravitationally produced radiation bath can exceed that of the conventional decay bath and account for the observed relic density provided that m χ > T RH , with additional k - and n -dependent constraints. For each freeze-in interaction considered, we also find m χ - and T RH -dependent limits on the beyond the Standard Model scale, Λ , for which gravitational production will exceed ordinary freeze-in production. Published by the American Physical Society 2025

Henrich, Stephen E. (ORCID:0000000323669519)↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

TeV Neutrinos and Hard X-Rays from Relativistic Reconnection in the Corona of NGC 1068

Abstract The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is dn p / dE p ∝ E p − 1 below the break and dn p / dE p ∝ E p − s above the break, with IceCube neutrino observations suggesting s ≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.

Astronomy & Astrophysics↗
DOE OSTI2024

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanisms of water oxidation on heterogeneous catalyst surfaces

Water oxidation, an essential step in photosynthesis, has attracted intense research attention. Understanding the reaction pathways at the electrocatalyst/water interface is of great importance for the development of water oxidation catalysts. How the water is oxidized on the electrocatalyst surface by the positive charges is still an open question. This review summarizes current advances in studies on surface chemistry within the context of water oxidation, including the intermediates, reaction mechanisms, and their influences on the reaction kinetics. The Tafel analyses of some electrocatalysts and the rate-laws relative to charge consumption rates are also presented. Moreover, how the multiple charge transfer relies on the intermediate coverage and the accumulated charge numbers is outlined. Finally, the intermediates and rate-determining steps on some water oxidation catalysts are discussed based on density functional theories.

14 SOLAR ENERGY↗
DOE OSTI2020

Toward fusion plasma scenario planning for NSTX-U using machine-learning-accelerated models

One of the most promising devices for realizing power production through nuclear fusion is the tokamak. To maximize performance, it is preferable that tokamak reactors achieve advanced operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic (MHD) stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high \emph{fusion gain} --- the ratio of produced fusion power to the external power provided through the plasma boundary. Precise and robust control of the evolution of the plasma boundary shape as well as the spatial distribution of the plasma current, density, temperature, and rotation will be essential to achieving and maintaining such scenarios. The complexity of the evolution of tokamak plasmas, arising due to nonlinearities and coupling between various parameters, motivates the use of model-based control algorithms that can account for the system dynamics. In this work, a learning-based accelerated model trained on data from the National Spherical Torus Experiment Upgrade (NSTX-U) is employed to develop planning and control strategies for regulating the density and temperature profile evolution around desired trajectories. The proposed model combines empirical scaling laws developed across multiple devices with neural networks trained on empirical data from NSTX-U and a database of first-principles-based computationally intensive simulations. The reduced execution time of the accelerated model will enable practical application of optimization algorithms and reinforcement learning approaches for scenario planning and control development. An initial demonstration of applying optimization approaches to the learning-based model is presented, including a strategy for mitigating the effect of leaving the finite validity range of the accelerated model. The approach shows promise for actuator planning between experiments and in real-time.

machine learning↗
DOE OSTI2020

A priori analysis of a power-law sub-filter scale mixing model for transported PDF modeling of high Karlovitz turbulent premixed flames

Accurate modeling of mixing in large-eddy simulation (LES) / transported probability density function (PDF) modeling of turbulent combustion remains an outstanding issue. The issue is particularly salient in turbulent premixed combustion under extreme conditions such as high-Karlovitz number Ka. Here, the present study addresses this issue by conducting an a priori analysis of a power-law scaling based mixing timescale model for the transported PDF model. A recently produced DNS dataset of a high-Ka turbulent jet flame is used for the analysis. A power-law scaling is observed for a scaling factor used to model the sub-filter scale mixing timescale in this high-Ka turbulent premixed DNS flame when the LES filter size is much greater than the characteristic thermal thickness of a laminar premixed flame. The sensitivity of the observed power-law scaling to the different viewpoints (local or global) and to the different scalars for the data analysis is examined and the dependence of the model parameters on the dimensionless numbers Ka and Re (the Reynolds number) is investigated. Different model formulations for the mixing timescale are then constructed and assessed in the DNS flame. The proposed model is found to be able to reproduce the mixing timescale informed by the high-Ka DNS flame significantly better than a previous model.

42 ENGINEERING↗
DOE OSTI2026

Electric accumulation of millicharged particles

A terrestrial population of millicharged particles that interact significantly with normal matter can arise if they make up a dark matter subcomponent or if they are light enough to be produced in cosmic-ray air showers. Such particles thermalize to terrestrial temperatures through repeated scatters with normal matter in Earth’s environment. We show that a simple electrified shell (for example, a Van de Graaff generator) functions as an efficient accumulator of such room-temperature millicharged particles, parametrically enhancing their local density by as much as 12 orders of magnitude. This can be used to boost the sensitivity of any detector housed in the shell’s interior, such as ion traps and tests of Coulomb’s law. In a companion paper, we apply this specifically to Cavendish tests of Coulomb’s law, and we show that a well-established setup can probe a large region of unexplored parameter space, with sensitivity to the irreducible density of millicharged particles generated from cosmic rays that outperforms future accelerator searches for sub-GeV masses.

Berlin, Asher [Fermilab] (ORCID:0000000211561482)↗
DOE OSTI2022

Pseudocapacitive Charge Storage in Electrochromic Transition-Metal Oxide Thin Films

Electrochromic pseudocapacitive transition-metal oxide materials, such as tungsten oxide, which combine fast response, high energy density, and optical effects, can play a significant role as energy storage materials. Here we investigate the electrochemical kinetics of thin films of tungsten oxide, which turn transparent to sky-blue color in the lithiated state due to the reduction of W 6+ to W 5+ . In this study, we investigated the charge density, charge transfer, ion diffusion, and interfacial behavior upon Li + insertion/de-insertion in WO 3 . The electrochromic thin film's pseudocapacitive and electrical double layer mechanism was differentiated based on the power-law. Faradaic diffusion-controlled process dominates over the surface capacitive behavior at scan rates below 40 mV s –1 . These films exhibit an areal charge density of around 100 mC cm –2 and a capacitance of 80 mF cm –2 , superior to most comparable electrochromic materials and supercapacitors. This work combines electrochromics and energy storage properties and provides a fundamental understanding of pseudocapacitive and electrochromic mechanisms in WO 3 .

25 ENERGY STORAGE↗
DOE OSTI2022

IER-501 CED-3b: Experiment Execution Summary for the Pulsed-Neutron Die-Away Experimental Testbed for Thermal Scattering Law Benchmarks (PNDA)

This report summarizes the experiments performed for IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments were conducted in two separate two-day campaigns: one campaign for high-density polyethylene and another for Lucite. All experiments were performed at LLNL. They will become high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high-priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). The report presents the PNDA design and its equipment, and it documents the experimental die-away curves. It presents the characterization measurements that have been performed to-date. Importantly, mass spectrometry measurements to determine sample impurities were not yet completed. These will be included in the final benchmark. The report gives the fitted decay constants for the die-away curves of each target sample. It also provides the data for the die-away curves in the appendices. The HDPE experiments examined twelve targets of varying size. The Lucite measurements included ten targets of varying size.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

IER-501 CED-3b: Experiment Execution Summary for the Pulsed-Neutron Die-Away Experimental Testbed for Thermal Scattering Law Benchmarks (PNDA)

This report summarizes the experiments done for IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments were conducted in two separate two-day campaigns: one campaign for high-density polyethylene and another for Lucite. All experiments were performed at Lawrence Livermore Laboratory. They will become high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). The report presents the PNDA design and its equipment, and it documents the experimental die-away curves. It presents the characterization measurements that have been performed to-date. Importantly, mass spectrometry measurements to determine sample impurities were not yet completed. These will be included in the final benchmark. The report gives the fitted decay constants for the die-away curves of each target sample. It also provides the data for the die-away curves in the appendices. The HDPE experiments examined twelve targets of varying size. The Lucite measurements included ten targets of varying size. Experimenters in the campaign were Daniel Siefman, William Zywiec, and Ruby Araj.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Spontaneous Symmetry–Breaking of Nonequilibrium Steady–States Caused by Nonlinear Electrical Transport

Negative differential resistance (NDR) in certain materials has been attributed to spontaneous emergence of symmetry-breaking electrical current density localization from a previously homogeneous distribution, which is postulated to occur due to the nonequilibrium thermodynamic force of minimization of entropy production. However, this phenomenon has not been quantitatively predicted based on intrinsic material properties and an applied electrical stimulus. Herein an instability criterion is derived for localization of current density and temperature from a thermal fluctuation in a parallel conductor model of a thin film that is subject to Newton's law of cooling. The conditions for steady–state electro-thermal localization is predicted, verifying a decrease in entropy production upon localization. Electro-thermal localization accompanied by a decrease of entropy production is confirmed in a multiphysics simulation of current flow in a thin film. The instability criterion predicts conditions for spontaneous current density localization, relating symmetry breaking fundamentally to dynamical instability via Local Activity theory.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Liquid-to-gas transfer of sodium in a liquid cathode glow discharge

Abstract Plasma-liquid interactions have been extensively studied with a focus on the transport of reactive species from the plasma to the liquid phase and their induced liquid phase chemistry and resulting applications. While solute transfer from the liquid to the gas phase in plasmas has been widely used in analytical chemistry, the underlying processes remain relatively unexplored. We report spatially and temporally resolved absolute density measurements of sodium in a plasma with a NaCl solution cathode using two-photon absorption laser induced fluorescence (TaLIF). The observed non-linear increase in sodium density with solution conductivity is shown to correlate with droplet generation as visualized by Mie scattering. The findings are explained by droplet generation by electrospray induced by Taylor cone formation as underpinning mechanism for the introduction of sodium in the plasma. An analytical sheath model combined with a scaling law shows an increase in electric field force with solution conductivity that is consistent with the observed non-linear increase in sodium density in the plasma with solution conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Cholesterol modulates membrane elasticity via unified biophysical laws

Cholesterol and lipid unsaturation underlie a balance of opposing forces that features prominently in adaptive cell responses to diet and environmental cues. These competing factors have resulted in contradictory observations of membrane elasticity across different measurement scales, requiring chemical specificity to explain incompatible structural and elastic effects. Here, we demonstrate that – unlike macroscopic observations – lipid membranes exhibit a unified elastic behavior in the mesoscopic regime between molecular and macroscopic dimensions. Using nuclear spin techniques and computational analysis, we find that mesoscopic bending moduli follow a universal dependence on the lipid packing density regardless of cholesterol content, lipid unsaturation, or temperature. Our observations reveal that compositional complexity can be explained by simple biophysical laws that directly map membrane elasticity to molecular packing associated with biological function, curvature transformations, and protein interactions. The obtained scaling laws closely align with theoretical predictions based on conformational chain entropy and elastic stress fields. These findings provide unique insights into the membrane design rules optimized by nature and unlock predictive capabilities for guiding the functional performance of lipid-based materials in synthetic biology and real-world applications.

Kumarage, Teshani [Virginia Polytechnic Inst. and ↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

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

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. Here, we address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗