Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Modeling”

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 433 records · Page 24

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of feed composition on the production of off-gases during vitrification of simulated low-activity nuclear waste

During the vitrification of nuclear waste, hazardous and radioactive emissions are generated from the feed-to-glass conversion reactions, in addition to discharges from forced air bubbling and air inleakage. Although the major gaseous emissions are water vapor, nitrogen, and carbon dioxide, various monitored environmental pollutants are also released, such as nitrogen oxides or sulfur dioxide. In addition, reactions between organics and nitrates in the feed may also form products of incomplete combustion such as carbon monoxide and acetonitrile. Although off-gas emissions are commonly measured during both laboratory- and pilot-scale melter testing, no predictive tool is currently available to a priori estimate the composition of gaseous emissions during nuclear waste vitrification. This work forms a basis for the development of such predictive tool by measuring gas evolution from a broad range of simulated low-activity waste melter feeds using evolved gas analysis data and developing correlations between the feed and off-gas compositions. Using reaction stoichiometry and regression analysis, we demonstrate that next to the content of nitrogen and organic carbon in the feed, the gaseous emissions are affected by the feed reduction-oxidation conditions – the more the feed is reduced, the less nitrogen monoxide, and more carbon monoxide and acetonitrile evolves. In conclusion, the results presented in this work provide a first step towards reducing the amount of expensive physical melter testing and the regression analysis provides a simple tool for rapid optimization of feed composition with respect to off-gas composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for the 1/2 + intruder state in 35 P

The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number Z along N = 20 is of interest both in terms of better understanding the evolution of nuclear structure between spherical Ca 40 and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the N = 20 Island of Inversion, the n p n h (where n = 2 , 4 , 6 ) neutron excitations across a diminished N = 20 gap result in deformed and collective ground states, as observed in Mg 32 . In heavier isotones, n p n h excitations do not dominate in the ground states but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2 p 2 h x s 1 / 2 + state in P 35 , the only N = 20 isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the S 36 ( d , He 3 ) P 35 reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross section to populate such a configuration within a 2.5 to 3.6 MeV energy range provides a stringent constraint on the wave function compositions in both S 36 and P 35 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fine structure in the odd-odd proton emitter Tm 144

Axial symmetry breaking in 144 Tm is probed by examining its proton emission fine structure. The ground-state spin and parity in 144 Tm and daughter 143 Er are assigned unambiguously based on the corroboration of our calculations with the present data. Here we establish the first microscopic description of fine structure in odd-odd proton emitters, which is capable of resolving ambiguities present in the assignments of transitions in such nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled-cluster theory for strong entanglement in nuclei

Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Examining the justification for the introduction of a fermion localization function

Becke and Edgecombe suggested in 1990 a theoretical tool to describe electron localization in atoms and molecules, an idea which was borrowed by a large number of nuclear theorists since 2011 to describe nucleon localization in nuclear systems. I argue here that these arguments are highly questionable and cannot be used in interacting systems where effects beyond the naive mean field or the simple Hartree-Fock framework are important and the inclusion of correlations induced by particle interactions is necessary in order to introduce such a localization function. Furthermore, I also describe several aspects of the exchange and irreducible two-body density matrices, which depend on the character and strength of the two-particle interaction and which can be useful in justifying the derivation of an appropriate energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proton Upset Monte Carlo Simulation

The Proton Upset Monte Carlo Simulation (PROPSET) program calculates the frequency of on-orbit upsets in computer chips (for given orbits such as Low Earth Orbit, Lunar Orbit, and the like) from proton bombardment based on the results of heavy ion testing alone. The software simulates the bombardment of modern microelectronic components (computer chips) with high-energy (.200 MeV) protons. The nuclear interaction of the proton with the silicon of the chip is modeled and nuclear fragments from this interaction are tracked using Monte Carlo techniques to produce statistically accurate predictions.

O'Neill, Patrick M.↗

Exploring the sensitivity of charge-exchange ( p , n ) reactions to the neutron density distribution

The determination of the nuclear neutron properties suffers from uncontrolled uncertainties, which have attracted considerable attention recently, such as in the context of the PREX experiment. Our aim is to analyze the sensitivity of charge-exchange (p,n) reactions to the neutron density distribution ρ n and constrain the neutron characteristics in the nuclear structure models. Method: By combing the folding and the mean-field models, the nucleon-nucleus (N,A) potential can be obtained from the nuclear density distribution. Further, the (p,p) and (p,n) cross sections for 48 Ca and 208 Pb are calculated following the distorted-wave Born approximation method. As a result, compared with the (p,p) cross section, the effects of ρ n variation on the (p,n) cross section are significant, which is due to the impact of isovector properties. Based on the global folding model analyses of data, it is found that 48 Ca and 208 Pb have relatively large neutron skin thickness ΔR np . Results illustrate that the charge-exchange (p,n) reaction is a sensitive probe of ρ n . The results in this paper can offer useful guides for future experiments of neutron characteristics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Solar neutrino detection in liquid xenon detectors via charged-current scattering to excited states

We investigate the prospects for real-time detection of solar neutrinos via the charged-current neutrino-nucleus scattering process in liquid xenon time projection chambers. We use a nuclear shell model, benchmarked with experimental data, to calculate the cross sections for populating specific excited states of the cesium nuclei produced by neutrino capture on Xe 131 and Xe 136 . The shell model is further used to compute the decay schemes of the low-lying 1 + excited states of Cs 136 , for which there is sparse experimental data. We explore the possibility of tagging the characteristic deexcitation γ rays/conversion electrons using two techniques: spatial separation of their energy deposits using event topology and their time separation using delayed coincidence. The efficiencies in each case are evaluated within a range of realistic detector parameters. We find that the topological signatures are likely to be dominated by radon backgrounds, but that a delayed-coincidence signature from long-lived states predicted in Cs 136 may enable background-free detection of CNO neutrino interactions in next-generation experiments with smaller uncertainty than current measurements. We also estimate the sensitivity as a function of exposure for detecting the solar-temperature-induced line shift in Be 7 neutrino emission, which may provide a new test of solar models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. Refined nuclear interaction models in these energy regimes will also be valuable for other applications, such as measurements of reactor, solar, and atmospheric neutrinos. This manuscript discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Projection algorithm for state preparation on quantum computers

We present an efficient method to prepare states of a many-body system on quantum hardware, first isolating individual quantum numbers and then using time evolution to isolate the energy. Our method in its simplest form requires only one additional auxiliary qubit. The total time evolved for an accurate solution is proportional to the ratio of the spectrum range of the trial state to the gap to the lowest excited state, a substantial improvement over other projection algorithms, and the accuracy increases exponentially with the time evolved. Isolating the quantum numbers is efficient because of the known eigenvalues and increases the gap thus shortening the propagation time required. The success rate of the algorithm, or the probability of producing the desired state, is a simple function of measurement times and phases and is dominated by the square overlap of the original state to the desired state. We present examples from the nuclear shell model and the Heisenberg model. In conclusion, we compare this algorithm to previous algorithms for short evolution times and discuss potential further improvements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Halo effective field theory analysis of one-neutron knockout reactions of Be 11 and C 15

Background: One-nucleon knockout reactions provide insightful information on the single-particle structure of nuclei. When applied to one-neutron halo nuclei, they are purely peripheral, suggesting that they could be properly modeled by describing the projectile within a halo effective field theory (halo-EFT). Purpose: We reanalyze the one-neutron knockout measurements of 11 Be and 15 C —both one-neutron halo nuclei—on beryllium at about 60 MeV/nucleon. We consider halo-EFT descriptions of these nuclei which already provide excellent agreement with breakup and transfer data. Method: Here we include a halo-EFT description of the projectile within an eikonal-based model of the reaction and compare its outcome to existing data. Results: Excellent agreement with experiment is found for both nuclei. The asymptotic normalization coefficients inferred from this comparison confirm predictions from ab initio nuclear-structure calculations and values deduced from transfer data. Conclusions: Halo-EFT can be reliably used to analyze one-neutron knockout reactions measured for halo nuclei and test predictions from state-of-the-art nuclear structure models on these experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influence of solar heating and precipitation scavenging on the simulated lifetime of post-nuclear war smoke

The behavior of smoke injected into the atmosphere by massive fires that might follow a nuclear war was simulated. Studies with a three-dimensional global atmospheric circulation model showed that heating of the smoke by sunlight would be important and might produce several effects that would decrease the efficiency with which precipitation removes smoke from the atmosphere. The heating gives rise to vertical motions that carry smoke well above the original injection height. Heating of the smoke also causes the tropopause, which is initially above the smoke, to reform below the heated smoke layer. Smoke above the tropopause is physically isolated from precipitation below. Consequently, the atmospheric residence time of the remaining smoke is greatly increased over the prescribed residence times used in previous models of nuclear winter.

Malone, R. C.↗

The harmonic oscillator and nuclear physics

The three-dimensional harmonic oscillator plays a central role in nuclear physics. It provides the underlying structure of the independent-particle shell model and gives rise to the dynamical group structures on which models of nuclear collective motion are based. It is shown that the three-dimensional harmonic oscillator features a rich variety of coherent states, including vibrations of the monopole, dipole, and quadrupole types, and rotations of the rigid flow, vortex flow, and irrotational flow types. Nuclear collective states exhibit all of these flows. It is also shown that the coherent state representations, which have their origins in applications to the dynamical groups of the simple harmonic oscillator, can be extended to vector coherent state representations with a much wider range of applicability. As a result, coherent state theory and vector coherent state theory become powerful tools in the application of algebraic methods in physics.

Rowe, D. J.↗

From light neutrino decay phenomenology to muon neutrino cross-section measurement at MINERvA experiment

This thesis tackles two essential topics of neutrino physics: neutrino decay and neutrino cross-section measurement. First, the invisible and visible neutrino decay is analyzed through a phenomenological approach, considering future long-baseline neutrino experiments such as DUNE and a hypothetical neutrino beam toward the ANDES laboratory. The study takes into account the ν μ and ν e disappearance and appearance, respectively, for both FHC and RHC flux modes. The results showed a negligible matter effect for DUNE but significantly more notable at ANDES. At 90% C.L., the sensitivity to the decay parameter α3 can be as small as 2 × 10 -6 eV 2 for a chosen coupling. The impact of neutrino decay in the determination of θ 23 and δ CP were also shown. Second, the double-differential cross-section measurement for ν μ -carbon interactions with three-momentum transfer |$\vec{q}$| < 1.2 GeV obtained with medium energy exposures in the NuMI beam at MINERvA experiment are reported. The measurement is presented as a function of |$\vec{q}$| and E avail and reviews different interaction models and nuclear effects along quasi-elastic to resonance processes to define a new model for a better agreement. The double differential cross sections are compared to the MnvTunes, GENIE, and NuWro predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗