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

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗
DOE OSTI2021

Quantum information approach to high energy interactions

High energy hadron interactions are commonly described by using a probabilistic parton model that ignores quantum entanglement present in the light-cone wave functions. Here, we argue that since a high energy interaction samples an instant snapshot of the hadron wave function, the phases of different Fock state wave functions cannot be measured—therefore the light-cone density matrix has to be traced over these unobservable phases. Performing this trace with the corresponding U(1) Haar integration measure leads to ‘Haar scrambling’ of the density matrix, and to the emergence of entanglement entropy. This entanglement entropy is determined by the Fock state probability distribution, and is thus directly related to the parton structure functions. As proposed earlier, at large rapidity η the hadron state becomes maximally entangled, and the entanglement entropy is S E ~η according to QCD evolution equations. When the phases of Fock state components are controlled, for example in spin asymmetry measurements, the Haar average cannot be performed, and the probabilistic parton description breaks down.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quasiclassical Gluon Fields and Low’s Soft Theorem at Small Momentum-Fraction x

In the high-energy limit, soft gluons can be approximately described by quasiclassical gluon fields. It is well known that the gluon field is a pure gauge field on the transverse plane at eikonal order. We derived the complete next-to-eikonal order solutions of the classical Yang-Mills equations for soft gluons in the dense nuclear regime. Utilizing these solutions, it is shown that Low’s soft theorem at small x can be obtained by considering off-diagonal matrix elements of quasiclassical chromoelectric field between single-gluon states in the dilute regime. We further propose on extending Low’s soft theorem at small x to incorporate the effects of gluon saturation in the dense regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The Neutron Skin of 48 Ca and 208 Pb: A Critical Analysis

We discuss the neutron skins of 48 Ca and 208 Pb. We review and critically examine modern predictions and empirical constraints, with special attention to the different interpretations of the findings from the PREX-II experiment and the recently reported value of the neutron skin in 48 Ca extracted from the CREX experiment. We argue that, in the spirit of the ab initio philosophy, the path to understanding the behavior of dense neutron-rich matter must not circumvent fundamental nuclear forces. Based only on that argument, a thick neutron skin in 208 Pb is highly unlikely.

chiral effective field theory↗
DOE OSTI2021

Structure of Stellar Remnants with Coupling to a Light Scalar

In this paper we study how a Yukawa coupling of the Standard Model fermions to a light scalar field effects the stellar structure of cold stellar remnants such as neutron stars. We elucidate the stellar structure phenomenology using a simple model of a massive scalar coupled to a single dominant fermion with no other interactions. For a broad scalar mass range ($10^{-10}\,\mathrm{eV}\ll m_\phi\ll10^3\,\mathrm{eV}$ for neutron stars) we show that the equation-of-state and stellar structure depends only the effective coupling $\mathfrak{g}=\frac{g_f\,m_f}{m_\phi}$, where $g_f$ is the Yukawa coupling, $m_f$ the fermion mass and $m_\phi$ is the scalar kinematic mass at nuclear densities. If $\mathfrak{g}>\mathcal{O}(1)$ the Yukawa coupled matter exhibits various anomalous behaviors including hydrodynamic instability, negative pressure, distinct phases (soft and hard) of matter with sharp phase boundaries between them and with vacuum. These anomalies can lead to stars consisting of only soft, only hard or hybrid of soft and hard matter. These stars can have varying signs of the slope of the mass-radius relation, anomalously large and small masses, gaps in allowed radii, multiple radii for the same mass, very thin crusts and radiate anomalously large amounts of energy when they form (in the form of neutrinos for neutron stars). To the extent that these anomalies have not and/or will not be observed limits the effective coupling to $\mathfrak{g}<\mathcal{O}(1)$. We argue this phenomenology is generic to realistic models of stars with Yukawa coupled matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Preliminary Implementation of Two-Dimensional Cartesian Solver in CTF-R

Sub-channel codes are one of the the modeling and simulation tools used for thermal-hydraulic analysis of nuclear reactors. A few examples of such sub-channel codes are the COolant Boiling in Rod Arrays (COBRA) family of codes. The approximations that are used to simplify the fluid conservation equations into sub-channel form, mainly that of axially-dominated flow, lead to noticeable limitations on sub-channels solvers for problems with significant flow in lateral directions. In this report, a two-dimensional Cartesian solver is developed and implemented within CTF-R, which is the residual solver in the North Carolina State University version of COBRA-TF (CTF). The new solver will enable CTF to simulate flow that is not axially-dominated. The appropriate Cartesian forms of the conservation equations are derived and implemented in the solver. Once the conservation equations are established, the process of constructing the matrix system was altered to solve a two-dimensional staggered grid system. A simple case was used to test that the two-dimensional Cartesian solver is accurate. The test problem does not include any source terms or flow in the lateral direction. The results show that the solver was able to run the simple case and converge to a steady-state solution. Future work will focus on testing existing capabilities by using test cases that include transients and equation cross-terms. Future work will also include adding additional capabilities such as enabling the solver to include cases with source terms and three dimensional cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

MAGNET Scaling and Methodology

The purpose of this study was to analyze the heat transfer of the Microreactor Agile Non-nuclear Experimental Testbed (MAGNET) within the Dynamic Energy Transport and Integration Laboratory (DETAIL) and develop scaling equations and models to couple with other systems. Hierarchical Two-Tiered Scaling (H2TS) methodologies were applied to DETAIL’s MAGNET facility to scale and project data sets while conserving the observed behavior based on first principles. The MAGNET system was successfully scaled using H2TS and multiple system parameters were determined or calculated from experimental data including steady state and transient data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Long-Range Transverse-Momentum Correlations and Radial Flow in Pb-Pb Collisions at the LHC

This Letter presents measurements of long-range transverse-momentum correlations using a new observable, 𝑣 0⁡ (𝑝 T ), serving as a probe of event-by-event radial-flow fluctuations, the underlying radial expansion, and the medium’s properties in heavy-ion collisions. Results are reported for inclusive charged particles, pions, kaons, and protons across various centrality intervals in Pb-Pb collisions at $\sqrt{𝑠_{\textrm{NN}}}$ = 5.02 TeV, recorded by the ALICE detector. A pseudorapidity-gap technique, similar to that used in anisotropic-flow studies, is employed to suppress short-range correlations. At low 𝑝 T , a characteristic mass ordering consistent with hydrodynamic collective flow is observed. At higher 𝑝 T (>3 GeV/𝑐), protons exhibit larger 𝑣 0 ⁡(𝑝 T ) than pions and kaons, in agreement with expectations from quark-recombination models. Comparisons to viscous hydrodynamic calculations with varying bulk viscosity and equation of state demonstrate the sensitivity of the 𝑣 0 ⁡(𝑝 T ) observable to these key medium properties. The findings establish 𝑣 0 ⁡(𝑝 T ) as a valuable addition to the set of observables used in Bayesian analyses for extracting the transport properties and constraining the equation of state of strongly interacting matter, while also helping to systematically explore its sensitivity and impact within such global studies.

Acharya, S. [INFN] (ORCID:0000000292135329)↗
DOE OSTI2024

A Reduced-Order Model of a Nuclear Power Plant with Thermal Power Dispatch

This paper presents reduced-order modeling of thermal power dispatch (TPD) from a pressurized water reactor (PWR) for providing heat to nearby heat consuming industrial processes that seek to take advantage of nuclear heat to reduce carbon emissions. The reactor model includes the neutronics of the reactor core, thermal–hydraulics of the primary coolant cycle, and a three-lump model of the steam generator (SG). The secondary coolant cycle is represented with quasi-steady state mass and energy balance equations. The secondary cycle consists of a steam extraction system, high-pressure and low-pressure turbines, moisture separator and reheater, high-pressure and low-pressure feedwater heaters, deaerator, feedwater and condensate pumps, and a condenser. The steam produced by the SG is distributed between the turbines and the extraction steam line (XSL) that delivers steam to nearby industrial processes, such as production of clean hydrogen. The reduced-order simulator is verified by comparing predictions with results from separate validated steady-state and transient full-scope PWR simulators for TPD levels between 0% and 70% of the rated reactor power. All simulators indicate that the flow rate of steam in the main steam line and turbine systems decrease with increasing TPD, which causes a reduction in PWR electric power generation. The results are analyzed to assess the impact of TPD on system efficiency and feedwater flow control. Due to the simplicity of the proposed reduced-order model, it can be scaled to represent a PWR of any size with a few parametric changes. In the future, the proposed reduced-order model will be integrated into a power system model in a digital real-time simulator (DRTS) and physical hardware-in-the-loop simulations.

08 HYDROGEN↗
DOE OSTI2024

Perspectives on Next Steps in Arms Control and Non-Proliferation

Much has happened in arms control over the last year as this chapter has been in development. The prospects for arms control as a viable national security tool have, if anything, grown darker. Russia has falsely equated tearing down the remaining arms control edifice with creating potential leverage vis-à-vis its actions in Ukraine. China has failed to update its talking points on being a responsible arms control actor with its massive nuclear force expansion. The United States has stated clearly that it is open for business on arms control and risk reduction, but it unfortunately has no customers in Russia and China. It remains unclear what specifically the United States is selling in terms of what it wants in a potential agreement. Competitive dynamics and an atmosphere of great power competition are driving all parties further away from the table.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Spatial and momentum imaging of the pion and kaon

We determine the leading Fock state light front wave functions (LFWFs) of the pion and kaon via light front projections of the covariant Bethe-Salpeter wave function. Using these LFWFs we study the multidimensional images of the valence quarks in the pion and kaon that are provided by their generalized parton distribution functions (GPDs) and transverse momentum dependent parton distribution functions (TMDs). Moments of the GPDs are taken to obtain the electromagnetic and gravitational form factors of the pion and kaon, and comparisons to available experimental and lattice data are made. Highlights from this study include predictions that the mean-squared impact parameter for the quarks in the pion and kaon are: ( b T 2 ) u π = 0.11 fm 2 , ( b T 2 ) s K = 0.08 fm 2 , and ( b T 2 ) u K = 0.13 fm 2 , and therefore the s quark in the kaon is much closer to the center of transverse momentum than the u quark. From the electromagnetic and gravitational form factors we find that the light-cone energy radii are about 60% smaller than the light-cone charge radii for each quark sector in the pion and kaon. A quantitative measure of the importance of the leading Fock state is obtained via comparison with a full Dyson-Schwinger equations calculation (containing an infinite tower of Fock states) for the pion form factor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Creation Synthetic Data to Train a Digital Twin to Predict Reactor Operations

Understanding techniques to strengthen the nuclear safeguards regime is crucial in preventing nuclear proliferation due to recent advancements in the nuclear energy industry such as Generation IV reactors and microreactors. Prior to the construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant's reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. A digital twin is defined as a virtual model that works in unison to represent a physical asset, with a transference of data between the virtual and physical assets [1]. This work serves as validation for training a digital twin on synthetic data fabricated via means of Serpent reactor physics and point kinetics equations simulations. In this case this work is based on parameters of Idaho State University's AGN-201 reactor. The synthetic data can then be utilized to train machine learning models in the future to further investigate the utility of these methods. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and decide whether these methods should be used in the future to inform inspectors of a reactor's proliferation potential.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

The Completion of Surrogate Testing for Low-Temperature TEX and a Look Towards the Future

To address the mounting need for below room temperature nuclear data validation, the Low-Temperature Thermal Epithermal eXperiments (LT-TEX) have been designed. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Neutron stars with exceptionally light QCD axions

We present a comprehensive study of axion condensed neutron stars that arise in models of an exceptionally light axion that couples to quantum chromodynamics (QCD). These axions solve the strong-charge-parity (𝐶⁢𝑃) problem, but have a mass-squared lighter than that due to QCD by a factor of 𝜖 < 1. Inside dense matter, the axion potential is altered, and much of the matter in neutron stars resides in the axion condensed phase where the strong-𝐶⁢𝑃 parameter 𝜃 = 𝜋 and 𝐶⁢𝑃 remains a good symmetry. In these regions, masses and interactions of nucleons are modified, in turn changing the equation of state (EOS), structure, and phenomenology of the neutron stars. We take the first steps toward the study of the EOS of neutron star matter at 𝜃 = 𝜋 within chiral effective field theory and use relativistic mean field theory to deduce the resulting changes to nuclear matter and the neutron star low-density EOS. We derive constraints on the exceptionally light axion parameter space based on observations of the thermal relaxation of accreting neutron stars, isolated neutron star cooling, and pulsar glitches, excluding the region up to 5 ×10 −7 ≲ 𝜖 ≲ 0.2 for 𝑚 𝑎 ≳ 2 × 10 −9 eV. We comment on potential changes to the neutron star mass-radius relationship, and discuss the possibility of novel, nuclear-density compact objects with 𝜃 = 𝜋 that are stabilized not by gravity but by the axion potential.

axions↗
DOE OSTI2025

Perturbative treatment of nonlocal chiral interactions in auxiliary-field diffusion Monte Carlo calculations

Nuclear many-body systems, ranging from nuclei to neutron stars, are some of the most interesting physical phenomena in our universe, and quantum Monte Carlo (QMC) approaches are among the most accurate many-body methods currently available to study them. In recent decades, interactions derived from chiral effective field theory (EFT) have been widely adopted in the study of nuclear many-body systems. One drawback of the QMC approach is the requirement that the nuclear interactions need to be local, whereas chiral EFT interactions usually contain nonlocalities. In this work, we leverage the capability of computing second-order perturbative corrections to the ground-state energy in order to develop a self-consistent approach to including nonlocal operators in QMC calculations. In conclusion, we investigate both the deuteron and the neutron-matter equation of state in order to show the robustness of our technique and pave the way for future QMC calculations at higher orders in the EFT, where nonlocal operators cannot be avoided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

On the nonadiabatic collisional quenching of OH(A) by H 2 : a four coupled quasi-diabatic state description

A four-state diabatic potential energy matrix (DPEM), H d , for the description of the nonadiabatic quenching of OH(A 2 Σ + ) by collisions with H 2 is reported. The DPEM is constructed as a fit to adiabatic energies, energy gradients, and derivative couplings obtained exclusively from multireference configuration interaction wave functions. A four-adiabatic-electronic-state representation is used in order to describe all energetically accessible regions of the nuclear coordinate space. Partial permutation-inversion symmetry is incorporated into the representation. Here, the fit is based on electronic structure data at 42 882 points, described by over 1.6 million least squares equations with a root mean square (mean unsigned) error of 178(83) cm –1 . Comparison of ab initio and H d determined minima, saddle points, and energy minimized points on C 2v , C s , C ∞v , and C 1 (noncoplanar) portions of two conical intersection seams are used to establish the accuracy of the H d .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

New Yield Estimates for Nuclear Detonations Over Water

One of the primary methods used by the United States to estimate the yields of atmospheric nuclear tests performed in the 1940’s, 50’s, and 60’s was to film each detonation with a high-speed camera, measure the blast radius as a function of time, and then use Taylor’s spherical blast wave equation to calculate the source energy (i.e., the yield). For surface detonations, the yield in Taylor’s equation was assumed to be twice as large as the actual yield to account for the energy reflected from the surface, assuming the ground was a perfect reflector. Recent simulations of surface detonations over water (referred to as barge shots) have shown that Taylor’s equation must be modified further to account for material entrainment and heat losses. We have computed this lost energy using a high-order radiation-hydrodynamics code to determine the reduction in blast radius vs. time. By analyzing simulation data in the same manner as film data, we have determined that a yield correction factor of 1.27± 0.04 should be applied to the barge shots.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

A Semi-Empirical Density Law for Ternary, Homogeneous PuCl 3 /HCl/H 2 O Solutions

Nuclear material operations pose unique hazards that are not encountered in other chemical, energy, or manufacturing industries. One of these hazards is the potential for a nuclear criticality accident when handling fissile isotopes such as 235 U and 239 Pu. These hazards are particularly high when fissile material is dissolved in solution as the neutron behaviors of the system can change rapidly with the physical and chemical changes accessible in solution. Current estimates of solution density used for criticality safety are outdated and hinder fissile material handling. Developing new estimates for these safety calculations requires experimental characterization and the derivation of empirical density models. We have derived a density law describing PuCl 3 /HCl/H 2 O solutions from experimental data characterizing solution density. Density data was treated using a Pitzer-derived eight-parameter equation, defining density as a function of analyte concentrations, temperature, and interactions between these variables. The model is predictive across the concentration and temperature ranges from which it was derived. The potential effects of varying oxidation states of plutonium, which are easily accessible in aqueous media, on the bulk solution density of the ternary system were also investigated. The resulting Pitzer-derived density law was applied to a nuclear criticality safety model, and the impact of the experimental characterization of solution density relative to previous estimates was demonstrated to be significant and suggest that the current approach to estimating density in nuclear criticality safety calculations may lead to overly conservative controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗