Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Equation of State”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Liouvillian dynamics of the open Schwinger model: String breaking and kinetic dissipation in a thermal medium

Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Sublimation temperature sensor for temperature locale

The sublimation temperature sensor (or “sublime sensor”) provides a continuum of measurement locations in which certain maximum temperatures are achieved during a heat up/cool down cycle. A predetermined material is encapsulated within a vacuum-sealed, non-volatile long tube (i.e., both ends capped and L$\gg$ D). This assembly is then inserted and centered into a heated zone, such as a furnace, exhaust pipe, or reactor. As the temperature increases, the material will sublimate (i.e., a process of having both the solid and gaseous states of matter simultaneously present) and will begin to fill the void - moving outward in both directions toward the ends of the tube. Once beyond the elevated temperatures, the gas will de-sublimate (i.e., deposition) onto the inner wall of the tube. The desired result of the sensor is the ring of material that develops over a relatively short period of time. This material deposit can be equated with temperature at an exact location. There is no need to interpolate and/or extrapolate for the desired measurement. Accuracy has been recorded for temperature locations on the range of ±2 mm over a 1 m span. Likewise, the precision of the measurement is ±0.2% the overall sensor domain. Furthermore, individual tubes with unique materials and pressures can be bundled together to provide a complete temperature profile of the heated zone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Annual Summary Report (FY 2025) Performance Assessment for the Environmental Restoration Disposal Facility

Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) started accepting low-level radioactive, hazardous, and mixed wastes that were generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells or disposal areas and can accommodate future design expansions as needed. Currently, there are eight cells and two supercells in ERDF. Each supercell is the equivalent of two cells. During this reporting period, fiscal year (FY) 2025, extending from October 1, 2024, through September 30, 2025, approximately 50,653 metric tons (55,835 U.S. tons) of waste was disposed at ERDF. From ERDF inception through September 30, 2025, approximately 17.3 million metric tons (19.1 million U.S. tons) of waste has been disposed at ERDF, which equates to consumption of approximately 91.0% of the currently constructed disposal volume. According to the design of ERDF, the facility has the ability to be expanded as needed. In 2023, the Department of Energy approved the construction of supercell 11 to provide additional disposal capacity to meet the cleanup mission at the Hanford Site. The construction of supercell 11 began in FY 2025, and it is expected to take up to two years. Construction of the new cell will not affect ongoing waste disposal activities at ERDF.

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

A new-old approach to composite scalars with chiral fermion constituents

We develop a dynamical, Lorentz invariant theory of composite scalars in configuration space consisting of chiral fermions, interacting by the perturbative exchange of a massive “gluon” of coupling g 0 and mass M 0 2 (the coloron model). The formalism is inspired by, but goes beyond, old ideas of Yukawa and the Nambu-Jona-Lasinio (NJL) model. It yields a non-pointlike internal wave-function of the bound state, ϕ ( r ) , which satisfies a Schrödinger-Klein-Gordon (SKG) equation with eigenvalue μ 2 . For super-critical coupling, g 0 2 > g 0 c 2 , we have μ 2 < 0 leading to spontaneous symmetry breaking. The binding of chiral fermions is semiclassical, not loop-level as in NJL. The mass scale is determined by the interaction as in NJL. We mainly focus on the short-distance, large M 0 2 limit, yielding an NJL pointlike interaction, but the bound state internal wave-function, ϕ ( r → ) , remains spatially extended and dilutes ϕ ( 0 ) . This leads to power-law suppression of the induced Yukawa and quartic couplings and requires radically less fine-tuning of a hierarchy than does the NJL model. We include a discussion of loop corrections of the theory. A realistic top condensation model appears possible.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Effect of magnetic fields on Urca rates in neutron star mergers

Isospin-equilibrating weak processes, called “Urca” processes, are of fundamental importance in astrophysical environments like (proto-)neutron stars, neutron star mergers, and supernovae. In these environments, matter can reach high temperatures of tens of MeVs and be subject to large magnetic fields. We thus investigate Urca rates at different temperatures and field strengths by performing the full temperature and magnetic-fielddependent rate integrals for different equations of state. We find that the magnetic fields play an important role at temperatures of a few MeV, especially close to or below the direct Urca threshold, which is softened by the magnetic field. At higher temperatures, the effect of the magnetic fields can be overshadowed by the thermal effects. Finally, we observe that the magnetic field influences the neutron decay rates more strongly than the electron-capture rates, leading to a shift in the flavor equilibrium.

electroweak interactions in nuclear physics↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Isospin mixing and the cubic isobaric multiplet mass equation in the lowest T=2, A=32 quintet

The isobaric multiplet mass equation (IMME) is known to break down in the first T=2,A=32 isospin quintet. In this work we combine high-resolution experimental data with state-of-the-art shell-model calculations to investigate isospin mixing as a possible cause for this violation. The experimental data are used to validate isospin-mixing matrix elements calculated with newly developed shell-model Hamiltonians. Our analysis shows that isospin mixing with nonanalog T=1 states contributes to the IMME breakdown, making the requirement of an anomalous cubic term inevitable for the multiplet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Superconductivity out of a non-Fermi liquid: Free energy analysis

In this paper, we present an in-depth analysis of the condensation energy E c for a superconductor in a situation when superconductivity emerges out of a non-Fermi liquid due to pairing mediated by a massless boson. This is the case for electronic-mediated pairing near a quantum-critical point in metal, for pairing in SYK-type models, and for phonon-mediated pairing in the properly defined limit, when the dressed Debye frequency vanishes. We consider a subset of these quantum-critical models, in which the pairing in a channel with a proper spatial symmetry is described by an effective 0 + 1 dimensional model with the effective dynamical interaction V(Ω m ) = g¯ γ /|Ω m | γ , where γ is model-specific (the γ model). In previous papers, we argued that the pairing in the γ model is qualitatively different from that in a Fermi liquid, and the gap equation at T = 0 has an infinite number of topologically distinct solutions, Δ n (ω m ), where an integer n, running between 0 and infinity, is the number of zeros of Δ n (ω m ) on the positive Matsubara axis. This gives rise to the set of extrema of E c at E c,n , of which E c,0 is the global minimum. Here, the spectrum E c,n is discrete for a generic γ < 2 but becomes continuous at γ = 2–0. Here, we discuss in more detail the profile of the condensation energy near each E c,n and the transformation from a discrete to a continuous spectrum at γ → 2. We also discuss the free energy and the specific heat of the γ model in the normal state.

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

Ducted Assembly Steady-State Heat Transfer Software (DASSH): Theory Manual

The Ducted Assembly Steady-State Heat transfer software (DASSH) is being developed at Argonne National Laboratory to obtain a steady-state thermal fluids result for the coolant flow and temperature distribution in a hexagonally gridded reactor core with ducted assemblies. In DASSH, each assembly is assumed to contain a hexagonal bundle of wire-wrapped pins. DASSH divides the coolant in the pin bundle into subchannels on which it applies energy balances to determine coolant and duct temperatures. DASSH offers multiple models for inter assembly heat transfer and features the ability to calculate pin temperatures. This document describes the models and equations implemented in DASSH. It includes a description of the geometric system to which the models are applied and a derivation of the subchannel energy balance formulation is presented. It also reviews the available correlations and provides an overview of the process by which power distributions are generated. Finally, it covers the code execution procedure, discussing the main code objects and their interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Determination of asymptotic normalization coefficients for the channel 16 O→ α+ 12 C: excited state 16 O(0 + ; 6.05 MeV)

Asymptotic normalization coefficients (ANC) determine the overall normalization of cross sections of peripheral radiative capture reactions. In the present paper, we treat the ANC C for the virtual decay 16 O(0 + ; 6.05 MeV)→ α+ 12 C(g.s.), the known values of which are characterized by a large spread (0.29 - 1.65) × 10 3 fm -1/2 . The ANC C is found by analytic continuation in energy of the α 12 C s-wave scattering amplitude, known from the phase-shift analysis of experimental data, to the pole corresponding to the 16O bound state and lying in the unphysical region of negative energies. To determine C, two different methods of analytic continuation are used. In the first method, the scattering data are approximated by the sum of polynomials in energy in the physical region and then extrapolated to the pole. The best way of extrapolation is chosen on the basis of the exactly solvable model. Within the second approach, the ANC C is found by solving the Schrödinger equation for the two-body α 12 C potential, the parameters of which are selected from the requirement of the best description of the phase-shift analysis data at a fixed experimental binding energy of 16 O(0 + ; 6.05 MeV) in the α+ 12 C channel. The values of the ANC C obtained within these two methods lie in the interval (886–1139) fm -1/2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

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

Is the trace anomaly at its minimum value at neutron star centers?

While the equation of state (EOS) 𝑃⁡(𝜖) of neutron star (NS) matter has been extensively studied, the EOS-parameter 𝜙 =𝑃/𝜖 or equivalently the dimensionless trace anomaly Δ =1/3 −𝜙, which quantifies the balance between pressure 𝑃 and energy density 𝜖, remains far less explored, especially in NS cores. Its bounds and density profile carry crucial information about the nature of superdense matter. Physically, the EOS-parameter 𝜙 represents the mean stiffness of matter accumulated from the stellar surface up to a given density. Based on the intrinsic structure of the Tolman-Oppenheimer-Volkoff equations, we show that 𝜙 decreases monotonically outward from the NS center, independent of any specific input NS EOS model. Furthermore, observational evidence of a peak in the SSS density profile near the center effectively rules out a valley and a subsequent peak in the radial profile of 𝜙 at similar densities, reinforcing its monotonic decrease. Finally, these model-independent relations impose strong constraints on the near-center behavior of the EOS-parameter 𝜙, particularly demonstrating that the mean stiffness (or equivalently Δ) reaches a local maximum (minimum) at the center.

nuclear astrophysics↗
DOE OSTI2019

Toward the implementation of self-powered, wireless, real-time reactor power sensing

The successful development of real-time, self-powered, wireless reactor power sensing technology would enable key operating parameters of every fuel assembly in a commercial reactor core to be continuously monitored, and potentially be utilized for applications such as closed-loop reactor control. This paper discusses the expansion of successful work on the application of Self Powered Neutron Detectors (SPND) for reactor power monitoring and the wireless transmission of these detector signals using advanced vacuum micro-electronic (VME) transmitter technology developed by Westinghouse in collaboration with the Penn State University Radiation Science and Engineering Center. Results demonstrate that this technology may be successfully applied via long-term testing in the Penn State TRIGA reactor. SPND are inherently slower than, for example, ion chambers due to their dependence on radioactive decay to generate a signal. Additional research developed delay compensation algorithms to accelerate the response of these detectors for use in near-real-time closed-loop reactor control. Initial results indicate that practical signal acceleration can be obtained digitally or with an equivalent analog circuit. Future research will apply an inverse model of the sensor’s balance differential equations to more-effectively accelerate the SPND signal. Use of inverse models will also prove to be adaptive in the sense that time-dependent burn-up of these detectors, due to neutron activation of the detector material, can be modeled and applied in the acceleration algorithm.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Khuri-Treiman equations for 3$π$ decays of particles with spin

Khuri-Treiman equations have proven to be a useful theoretical tool in the analysis of three-body decays, especially into the 3$π$ final state. In this work we present in full detail the necessary generalization of the formalism to study the decays of particles with arbitrary spin, parity, and charge conjugation. To this extent, we find it most convenient to work with helicity amplitudes instead of the so-called invariant amplitudes, especially when dealing with the unitarity relations. The isobar expansions in the three possible ( s-, t-, and u- ) final channels are related with the appropriate crossing matrices. We pay special attention to the kinematical singularities and constraints of the helicity amplitudes, showing that these can be derived by means of the crossing matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗