Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density perturbation”

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

Experimental and theoretical evidence of universality in superfluid vortex reconnections

The minimum separation between reconnecting vortices in fluids and superfluids obeys a universal scaling law with respect to time. The prereconnection and the postreconnection prefactors of this scaling law are different, a property related to irreversibility and to energy transfer and dissipation mechanisms. In the present work, we determine the temperature dependence of these prefactors in superfluid helium from experiments and a numeric model which fully accounts for the coupled dynamics of the superfluid vortex lines and the thermal normal fluid component. At all temperatures, we observe a pre- and postreconnection asymmetry similar to that observed in other superfluids and in classical viscous fluids, indicating that vortex reconnections display a universal behavior independent of the small-scale regularizing dynamics. We also numerically show that each vortex reconnection event represents a sudden injection of energy in the normal fluid. Finally we argue that in a turbulent flow, these punctuated energy injections can sustain the normal fluid in a perturbed state, provided that the density of superfluid vortices is large enough.

Reconnections↗
DOE OSTI2022

DIII-D research advancing the physics basis for optimizing the tokamak approach to fusion energy

DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ~8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co-I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anomalous isotope effect on the optical bandgap in a monolayer transition metal dichalcogenide semiconductor

Isotope effects have received increasing attention in materials science and engineering because altering isotopes directly affects phonons, which can affect both thermal properties and optoelectronic properties of conventional semiconductors. However, how isotopic mass affects the optoelectronic properties in 2D semiconductors remains unclear because of measurement uncertainties resulting from sample heterogeneities. Here, we report an anomalous optical bandgap energy red shift of 13 (±7) milli–electron volts as mass of Mo isotopes is increased in laterally structured 100 MoS 2 - 92 MoS 2 monolayers grown by a two-step chemical vapor deposition that mitigates the effects of heterogeneities. This trend, which is opposite to that observed in conventional semiconductors, is explained by many-body perturbation and time-dependent density functional theories that reveal unusually large exciton binding energy renormalizations exceeding the ground-state renormalization energy due to strong coupling between confined excitons and phonons. The isotope effect on the optical bandgap reported here provides perspective on the important role of exciton-phonon coupling in the physical properties of two-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Auger spectroscopy beyond the ultra-short core-hole relaxation time approximation

Abstract We present a time-dependent computational approach to study Auger electron spectroscopy (AES) beyond the ultra-short core-hole relaxation time approximation and, as a test case, we apply it to the paradigmatic example of a one-dimensional Mott insulator represented by a half-filled Hubbard chain. The Auger spectrum is usually calculated by assuming that, after the creation of a core-hole, the system thermalizes almost instantaneously. This leads to a relatively simple analytical expression that uses the ground-state with a core-hole as a reference state and ignores all the transient dynamics related to the screening of the core-hole. In this picture, the response of the system can be associated to the pair spectral function. On the other hand, in our numerical calculations, the core hole is created by a light pulse, allowing one to study the transient dynamics of the system in terms of the pulse duration and in the non-perturbative regime. Time-dependent density matrix renormalization group calculations reveal that the relaxation process involves the creation of a polarization cloud of doublon excitations that have an effect similar to photo-doping. As a consequence, there is a leak of spectral weight to higher energies into what otherwise would be the Mott gap. For longer pulses, these excited states, mostly comprised of doublons, can dominate the spectrum. By changing the duration of the light-pulse, the entire screening process can be resolved in time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Core States of Neutron Stars from Anatomizing Their Scaled Structure Equations

Given an Equation of State (EOS) for neutron star (NS) matter, there is a unique mass–radius sequence characterized by a maximum mass ${M}_{\mathrm{NS}}^{\max}$ at radius $R$ max . We first show analytically that the ${M}_{\mathrm{NS}}^{\max}$ and $R$ max scale linearly with two different combinations of the NS central pressure $P$ $c$ and energy density $ε$ $c$ , by dissecting perturbatively the dimensionless Tolman–Oppenheimer–Volkoff (TOV) equations governing NS internal variables. The scaling relations are then verified via 87 widely used and rather diverse phenomenological as well as 17 microscopic NS EOSs with/without considering hadron–quark phase transitions and hyperons, by solving numerically the original TOV equations. The EOS of the densest NS matter allowed before it collapses into a black hole is then obtained. Using the universal ${M}_{\mathrm{NS}}^{\max}$ and $R$ max scalings and Neutron Star Interior Composition Explorer and XMM-Newton mass–radius observational data for PSR J0740+6620, a very narrow constraining band on the NS central EOS is extracted directly from the data for the first time, without using any specific input EOS model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Diii-D Research Advancing The Physics Basis For Optimizing The Tokamak Approach To Fusion Energy

DIII-D physics research addresses critical challenges for operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-Ip steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs. ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. With this physics understanding and suite of actuators, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q=10 in ITER at 9 MA was coupled to a detached divertor, and a near Super H-mode optimized-pedestal scenario with co-Ip beam injection, was coupled to a radiative divertor. The Hybrid core scenario was achieved directly, without the need for internal flux pumping, using off-axis current drive actuators, and a controller to avoid mode locking and disruptions in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

Fenstermacher, M.↗
DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Interference with gravitational instability: Hot and fuzzy dark matter

Wave or fuzzy dark matter produced with high momenta behaves in many ways like hot particle dark matter while also possessing seemingly different phenomenology due to wave interference. We develop wave perturbation theory to show that white noise density fluctuations generated by the interference of high-momenta waves are gravitationally unstable in the usual way during matter domination above the free-streaming scale and stabilize below the free-streaming scale, much like the analogous effects for massive neutrinos in hot dark matter. We verify and illustrate these effects in the density power spectra of Newtonian Schrödinger-Poisson simulations. In the cosmological context, this would cause a gradual suppression of the initial white noise isocurvature perturbations below the free-streaming scale at matter radiation equality, unlike cold dark matter isocurvature fluctuations, and virial stability of dark matter halos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Carrier density oscillation in the photoexcited semiconductor

Abstract The perturbation of a semiconductor from the thermodynamic equilibrium often leads to the display of non-linear dynamics and formation of spatiotemporal patterns due to the spontaneous generation of competing processes. Here, we use scanning ultrafast electron microscopy to show that the transport of hot carriers in the strongly excited semiconductor slows down by turning into an oscillatory process; this is evidenced by the expansions and contractions in the second moment of the distribution. We attribute such a response to the electric field generated by the spatial separation of photo-excited electrons and holes under intrinsic and photo-induced fields; we then introduce a transport model that mimics the experimental observation. Our finding provides a direct imaging evidence for the electrostatic oscillation of hot carriers in the highly excited semiconductor and offers new insights into their dynamics in space and time.

Physics↗
DOE OSTI2022

Spin relaxation rate for heavy quarks in weakly coupled QCD plasma

We compute the relaxation rate of the spin density of heavy quarks in a perturbative QCD plasma to leading-log order in the coupling constant g. The spin relaxation rate Γ s in spin hydrodynamics is shown to be Γ s ~ g 4 log(1/g)T(T/M) 2 in the heavy-quark limit T/M $\ll$ 1, which is smaller than the relaxation rate of other non- hydrodynamic modes by additional powers of T/M. We demonstrate three different methods to evaluate the spin relaxation rate: 1) the Green-Kubo formula in the spin hydrodynamic regime, 2) the spin density correlation function in the strict hydrodynamic limit, and 3) quantum kinetic theory of the spin distribution function in momentum space. We highlight the interesting differences between these methods, while they are ultimately connected to each other by the underlying Ward-Takahashi identity for the non-conserved spin density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Large scale structure reconstruction with short-wavelength modes

Large scale density modes are difficult to measure because they are sensitive to systematic observational errors in galaxy surveys, but we can study them indirectly by observing their impact on small scale perturbations. Cosmological perturbation theory predicts that second-order density inhomogeneities are a convolution of a short- and a long-wavelength mode. This arises physically because small scale structures grow at different rates depending on the large scale environment in which they reside. This induces an off-diagonal term in the two-point statistics in Fourier space that we use as the basis for a quadratic estimator for the large scale field. Here, we demonstrate that this quadratic estimator works well on an N-body simulation of size (2.5 h -1 Gpc) 3 . In particular, the quadratic estimator successfully reconstructs the long-wavelength modes using only small-scale information. This opens up novel opportunities to study structure on the largest observable scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Development of real-time density feedback control on MAST-U in L-mode

In this paper we report on the development and demonstration of density feedback control for MAST-U. Sinusoidal perturbations are used to measure the frequency response from a deuterium gas valve (actuator) to line-integrated core electron density measured by the interferometer (sensor). In the frequency range relevant for control design, only two system-identification experiments were needed to regress a first-order dynamic model. This control-oriented model informs the offline design of a proportional integral controller with the established loop-shaping controller design method. After offline verification of the controller implementation, control is demonstrated by experimentally tracking a staircase reference for the line-integrated electron density. This paper demonstrates the efficiency of controller design using system-identification and loop-shaping, providing reliable density control for MAST-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗