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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

HdsaLib

HdsaLib performs hyper-differential sensitivity analysis in optimization problems that are constrained by partial differential equation problems. The software aids researchers in understanding the impact of uncertainties on their optimization models. The tool has been applied for use in ice sheet modeling, fluid flow modeling, atmospheric tracer transport modeling, and subsurface flow modeling. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

van Bloemen Waanders, Bart [Sandia National Lab. (↗
DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Magnetic field induced hair structure in the charmonium gluon dissociation

We study the electromagnetic field effect on charmonium gluon dissociation in quark-gluon plasma. With the effective Hamiltonian derived from QCD multipole expansion under an external electromagnetic field, we first solve the two-body Schrödinger equation for a pair of charm quarks with mean field potentials for color and electromagnetic interactions and obtain the charmonium binding energies and wave functions, and then calculate the gluon dissociation cross section and decay width by taking the color electric and magnetic dipole interactions as perturbations above the mean field and employing Fermi’s golden rule. Considering the charmonium deformation in a magnetic field, the discrete Landau energy levels make the dissociation cross section grow hair, and the electric dipole channel is significantly changed, especially for the P -wave states χ c 0 and χ c ± . From our numerical calculation, the magnetic field strength e B = 5 m π 2 already changes the gluon dissociation strongly, which may indicate measurable effects in high-energy nuclear collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Neutron Next-Event Estimators Kinematics (Rev.2)

This paper reviews the kinematics of neutron elastic and inelastic scattering with moving and stationary targets for contributions to Neutron Next-Event Estimators (NEEs). NEEs are often used for simulating the response of detectors in locations that have few random walk particles. Contributions to NEEs from collisions differ from sampling the outgoing particle state in standard random walk collisions as the location of the estimator is fixed, thus fixing the scattering angle between the incoming and outgoing directions. To calculate contributions to NEEs, the outgoing energies and probabilities of scatter toward the estimator must be calculated. This paper presents the most general case that encompasses elastic and inelastic scattering with both moving and stationary targets. The moving target case is required for elastic scattering with thermal motion due to the free gas thermal treatment approximation. The stationary target equations are presented for comparison with equations presented in canonical Monte Carlo texts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Computational schemes for the Magnus expansion of the in-medium similarity renormalization group

The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schrödinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to 7MeV for ground-state energies and 0.5MeV for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.

39 ≤ A ≤ 58↗
DOE OSTI2023

Quarkyonic mean field theory

We discuss mean field theory of quarkyonic matter at zero temperature. We treat the nucleons with contact interactions in mean field approximation, discussing both vector and scalar mean field interactions. We treat the quarks without mean field vector interactions, but allow mass terms to be generated consistently from a scalar mean field consistent with the additive quark model for quark masses. Quarkyonic matter is composed of a shell of nucleons that under-occupy the total available phase space associated with the underlying quark degrees of freedom. Here, the fully occupied Fermi sphere beneath this shell of nucleons at high densities is thought of as quarks, but when this fully occupied distribution of states first appears, although the phase space is filled, the matter is at low density. For the transition between this low density and high density saturated matter, we advocate a dual description of the fully filled Fermi sea in terms of hadrons, and make a phenomenological hypothesis for the equation of state of this matter. We then proceed to an example where the mean field interactions are all vector and only associated with the nucleons, ignoring the effects of mass change associated with the scalar interactions. Except for the effects of Pauli blocking, the nucleons and quarks do not interact. To get a reasonable transition to quarkyonic matter the interaction of the quarks among themselves are assumed to be nonperturbative, and a simple phenomenological relation between quark Fermi energy and density is introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coupled Boltzmann transport equations of heavy quarks and quarkonia in quark-gluon plasma

We develop a framework of coupled transport equations for open heavy flavor and quarkonium states, in order to describe their transport inside the quark-gluon plasma. Our framework is capable of studying simultaneously both open and hidden heavy flavor observables in heavy-ion collision experiments and can account for both, uncorrelated and correlated recombination. Our recombination implementation depends on real-time open heavy quark and antiquark distributions. We carry out consistency tests to show how the interplay among open heavy flavor transport, quarkonium dissociation and recombination drives the system to equilibrium. We then apply our framework to study bottomonium production in heavy-ion collisions. We include Υ(1S), Υ(2S), Υ(3S), χb (1P) and χb (2P) in the framework and take feed-down contributions during the hadronic gas stage into account. Cold nuclear matter effects are included by using nuclear parton distribution functions for the initial primordial heavy flavor production. A calibrated 2 + 1 dimensional viscous hydrodynamics is used to describe the bulk QCD medium. We calculate both the nuclear modification factor R AA of all bottomonia states and the azimuthal angular anisotropy coefficient v 2 of the Υ(1S) state and find that our results agree reasonably with experimental measurements. Our calculations indicate that correlated cross-talk recombination is an important production mechanism of bottomonium in current heavy-ion experiments. The importance of correlated recombination can be tested experimentally by measuring the ratio of R AA ( χb (1P)) and R AA (Υ(2S)).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
DOE OSTI2020

Molecular bond stabilization in the strong-field dissociation of O 2 +

We theoretically examine the rotational and vibrational dynamics of O 2 + molecular ions exposed to intense, short laser pulses for conditions realized in contemporary pump-probe experiments. We solve the time-dependent Schrödinger equation within the Born-Oppenheimer approximation for an initial distribution of randomly aligned molecular ions. For fixed peak intensities, our numerical results show that total, angle-integrated O 2 + → O( 3 P) + O + ( 4 S 0 ) dissociation yields do not monotonically increase with increasing infrared-probe pulse duration. We find this pulse-duration-dependent stabilization to be consistent with the transient trapping of nuclear probability density in a light-induced (bond-hardening) potential-energy surface and robust against rotational excitation. We analyze this stabilization effect and its underlying bond-hardening mechanism (i) in the time domain, by following the evolution of partial nuclear probability densities associated with the dipole-coupled O 2 + (a 4 Π u ) and O 2 + (f 4 Π g ) cationic states, and (ii) in the frequency domain, by examining rovibrational quantum-beat spectra for the evolution of the partial nuclear probability densities associated with these states. In this work, our analysis reveals the characteristic timescale for the bond-hardening mechanism in O 2 + and explains the onset of bond stabilization for sufficiently long pulse durations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2020

Quarkyonic matter equation of state in beta-equilibrium

Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasiparticles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass (≳ 2 $M_⊙$) and the radius of a 1.4 $M_⊙$ star ($R_{1.4}$ ≲ 13.5 km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. Here, we propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Explicit Modeling of Pebble Temperature in the Porous-medium Framework for Pebble-bed Reactors Applications

In this study, an explicit model for pebble temperature calculations in the framework of a porous-medium approach has been developed and implemented in the SAM code. The solid-phase energy conservation equation, including the explicit modeling of pebble temperature, is a multi-scale model, which can predict the macroscopic (pebble-bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. Extensive and successful code verifications and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the previously implemented model, which assumed a local solid-phase thermal equilibrium that is not appropriate for pebble-bed reactor design and safety analysis applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Thermal Effects in Binary Neutron Star Mergers

Abstract We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nuclear fragmentation at the future electron-ion collider

We explore the potential of conducting low-energy nuclear physics studies, including nuclear structure and decay, at the future Electron-Ion Collider (EIC) at Brookhaven. By comparing the standard theory of electron-nucleus scattering with the equivalent photon method applied to Ultraperipheral Collisions (UPC) at the Large Hadron Collider (LHC) at CERN. In the limit of extremely high beam energies and small energy transfers, very transparent equations emerge. We apply these equations to analyze nuclear fragmentation in UPCs at the LHC and 𝑒𝐴 scattering at the EIC, demonstrating that the EIC could facilitate unique photonuclear physics studies. Furthermore, we have also shown that the fragmentation cross-sections at the EIC are about 1,000 times smaller than those at the LHC. At the LHC, the fragmentation of uranium nuclei displays characteristic double-hump mass distributions from fission events, while at the EIC, fragmentation is dominated by neutron emission and fewer few fission products, about 10,000 smaller number of events.

Bertulani, Carlos A. [East Texas A&M University, C↗
DOE OSTI2024

Analysis of static Wilson line correlators from lattice QCD at finite temperature with T -matrix approach

The thermodynamic T-matrix approach is used to study Wilson line correlators (WLCs) for a static quark-antiquark pair in the quark-gluon plasma (QGP). Selfconsistent results that incorporate constraints from the QGP equation of state can approximately reproduce WLCs computed in 2+1-flavor lattice-QCD (lQCD), provided the input potential exhibits less screening than in previous studies. Utilizing the updated potential to calculate pertinent heavylight T-matrices we evaluate thermal relaxation rates of heavy quarks in the QGP. We find a more pronounced temperature dependence for low-momentum quarks than in our previous results (with larger screening), which turns into a weaker temperature dependence of the (temperature-scaled) spatial diffusion coefficient, in fair agreement with the most recent lQCD data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗