SEARCH · Engineering Papers
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.
The presentation is a status report of activities undertaken by the five laboratories (LANL, BNL, LBNL, PNNL, and LLNL) that are a part of the Energy-Dependent FPY Project Funded.
An extensive dataset of cumulative fission product yields has been generated under a joint collaboration between Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and the Triangle Universities Nuclear Laboratory (TUNL). The energy dependence of the cumulative yield for a select number of high-yield fission products has been measured using quasimonoenergetic neutrons with energies between 5.5 and 11.0 MeV. This is in addition to previously published data covering 0.5–4.5 and 14.8 MeV. The absolute number of fissions was determined during the irradiation period using dual-fission ionization chambers, and the fission products were measured postactivation by whole target γ-ray spectroscopy. Here this paper presents the absolute cumulative fission product yields as a function of incident neutron energy from the neutron-induced fission of 235 U, 238 U, and 239 Pu isotopes at four incident energies in the second-chance fission region and compares them with existing literature values. Corrections relevant to this collaboration's previously published data are also discussed.
The PERT card is used for first- and second-order perturbations of material densities, isotope densities, and isotope cross sections in the MCNP6.2 Monte Carlo code. Using the PERT card for firstand second-order sensitivities is straightforward. Reference 2 only looked at density perturbations with no energy dependence. This report looks at the sensitivity of a response to scattering cross sections as a function of energy and finds that there is a bug in the PERT capability. The SENSMG multigroup neutron sensitivity code was recently given the capability of computing sensitivities to user-specified reactions whose cross sections are available in a user-supplied NJOY output file. SENSMG uses the PARTISN multigroup discrete ordinates code for the transport. The NJOY-reaction capability is similar to that of SUSD3D, which also uses PARTISN for the multigroup transport. The PERT bug was identified by comparing MCNP6.2 outputs with SENSMG outputs. The next section of this report describes the test problem. Section III presents the MCNP6.2 PERT bug and proposes a solution. Section IV is a summary and conclusions. The SENSMG and MCNP input files are listed in the appendix.
Data from a survey of solar flare events undertaken by the UMD/MPI ULET telescope on the IMP-8 satellite during the period 1973-1977 yields examples of time and energy dependence of the abundances and spectra of the He, C, O and Fe heavy ions. Time variations are found in the O/He, O/C and Fe/O ratios which appear to be inconsistent with models based entirely on rigidity-dependent propagation in the interplanetary medium. It is speculated that such other factors as the abundance differences at the flare site on the sun, or features of the acceleration and release mechanisms, may play an important role in the variations discussed.
Here, we present a measurement of the first-order azimuthal anisotropy v 1 of deuterons from Au + Au collisions at $\sqrt{s_{NN}}$ = 7.7 , 11.5, 14.5, 19.6, 27, and 39 GeV recorded with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The energy dependence of the v 1 (y) slope, $dv_1/dy|_{y=0}$, for deuterons, where y is the rapidity, is extracted for semicentral collisions (10%–40% centrality) and compared with that of protons. While the v 1 (y) slopes of protons are generally negative for $\sqrt{s_{NN}}$ > 10 GeV , those for deuterons are consistent with zero, a strong enhancement of the v 1 (y) slope of deuterons is seen at the lowest collision energy (the largest baryon density) at $\sqrt{s_{NN}}$ = 7.7 GeV . In addition, we report the transverse momentum dependence of v 1 for protons and deuterons. The experimental results are compared with transport and coalescence models.
The cross section for coherent photonuclear production of J/ψ is presented as a function of the electromagnetic dissociation (EMD) of Pb. The measurement is performed with the ALICE detector in ultra-peripheral Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{{\textrm{s}}_{\textrm{NN}}}$ = 5.02 TeV. Cross sections are presented in five different J/ψ rapidity ranges within |y| < 4, with the J/ψ reconstructed via its dilepton decay channels. In some events the J/ψ is not accompanied by EMD, while other events do produce neutrons from EMD at beam rapidities either in one or the other beam direction, or in both. The cross sections in a given rapidity range and for different configurations of neutrons from EMD allow for the extraction of the energy dependence of this process in the range 17 < W γ Pb,n < 920 GeV, where W γ Pb,n is the centre-of-mass energy per nucleon of the γPb system. This range corresponds to a Bjorken-x interval spanning about three orders of magnitude: 1.1 × 10 —5 < x < 3.3 × 10 —2 . In addition to the ultra-peripheral and photonuclear cross sections, the nuclear suppression factor is obtained. These measurements point to a strong depletion of the gluon distribution in Pb nuclei over a broad, previously unexplored, energy range. These results, together with previous ALICE measurements, provide unprecedented information to probe quantum chromodynamics at high energies.
A Monte Carlo computational model has been developed which simulates atomic oxygen attack of protected polymers at defect sites in the protective coatings. The parameters defining how atomic oxygen interacts with polymers and protective coatings as well as the scattering processes which occur have been optimized to replicate experimental results observed from protected polyimide Kapton on the Long Duration Exposure Facility (LDEF) mission. Computational prediction of atomic oxygen undercutting at defect sites in protective coatings for various arrival energies was investigated. The atomic oxygen undercutting energy dependence predictions enable one to predict mass loss that would occur in low Earth orbit, based on lower energy ground laboratory atomic oxygen beam systems. Results of computational model prediction of undercut cavity size as a function of energy and defect size will be presented to provide insight into expected in-space mass loss of protected polymers with protective coating defects based on lower energy ground laboratory testing.
A comparison of light hypernuclei production, from Ultra-relativistic Quantum Molecular Dynamics (UrQMD)+coalescence and the thermal model, in heavy ion collisions over a wide range of beam energies and system sizes is presented. We find that both approaches provide generally similar results, with differences in specific details. Especially the ratios of hypertriton to Λ are affected by both the source radius Δr of the coalescence procedure as well as canonical effects. On the other hand, the double ratio S 3 is almost independent of canonical effects, which is in contrast to coalescence. Thus, both the beam energy dependence and centrality dependence of S 3 can be used to constrain the hypertriton source radius. To do so the currently available data are not yet sufficient. Elliptic flow is shown to be unaffected by the source size of the nuclei and an almost perfect mass scaling of the elliptic flow is observed. Furthermore, our predictions further suggest that the existence of the H-dibaryon(ΛΛ) seems ruled out by ALICE data.
In this study we focus on the radiation belt dynamics driven by the geomagnetic storms during September 2017. Besides the long-lasting three-belt structures of ultrarelativistic electrons (>2 MeV, existing for tens of days), which has been studied intensively during the Van Allen Probe era, it is found that magnetospheric electrons of hundreds of keVs can also have three-belt structures at similar L extent during storm time. Measurements of 500–800 keV electrons from MagEIS instrument onboard Van Allen Probes show double-peaked (L = 3.5 and 4.5, respectively) flux-versus-L-shell profile in the outer belt, which lasted for 2–3 days. During the time interval of such transient three-belt structure, the energy-versus-L spectrogram shows novel distributions differing from both “S-shaped” and “V-shaped” spectrograms reported previously. Such peculiar distribution also illustrates the energy-dependent occurrence of the three-belt profile. Furthermore, the gradual formation of “reversed energy spectrum” at L ~ 3.5 also indicates that hiss scattering inside the plasmapause contributed to the fast decay of sub-MeV remnant belt.
In the framework of the A Multi-Phase Transport (AMPT) model, the multi-particle azimuthal cumulant method is used to calculate the linear and mode-coupled contributions to the quadrangular flow harmonic (v4) and the mode-coupled response coefficient as functions of centrality in Au+Au collisions at sNN = 200, 39, 27 and 19.6 GeV. This study indicates that the linear and mode-coupled contributions to v4 are sensitive to beam energy change. Nevertheless, the correlations between different-order flow symmetry planes and the mode-coupled response coefficients show weak beam energy dependence. In addition, the presented results suggest that the experimental measurements that span a broad range of beam energies can be an additional constraint for the theoretical model calculations.
We report the first measurements of cumulants, up to 4th order, of deuteron number distributions and proton-deuteron correlations in Au+Au collisions recorded by the STAR experiment in phase-I of Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider. Deuteron cumulants, their ratios, and proton-deuteron mixed cumulants are presented for different collision centralities covering a range of center-of-mass energy per nucleon pair $\sqrt{s_{NN}} = 7.7$ to 200 GeV. It is found that the cumulant ratios at lower collision energies favor a canonical ensemble over a grand canonical ensemble in thermal models. An anti-correlation between proton and deuteron multiplicity is observed across all collision energies and centralities, consistent with the expectation from global baryon number conservation. The UrQMD model coupled with a phase-space coalescence mechanism qualitatively reproduces the collision-energy dependence of cumulant ratios and proton-deuteron correlations.
The linear and mode-coupled contributions to higher-order anisotropic flow are presented for Au+Au collisions at √ S NN = 27, 39, 54.4, and 200 GeV and compared to similar measurements for Pb+Pb collisions at the Large Hadron Collider (LHC). The coefficients and the flow harmonics' correlations, which characterize the linear and mode-coupled response to the lower-order anisotropies, indicate a beam energy dependence consistent with an influence from the specific shear viscosity (η/s). In contrast, the dimensionless coefficients, mode-coupled response coefficients, and normalized symmetric cumulants are approximately beam-energy independent, consistent with a significant role from initial-state effects. These measurements could provide unique supplemental constraints to (i) distinguish between different initial-state models and (ii) delineate the temperature (T) and baryon chemical potential (μ B ) dependence of the specific shear viscosity η/s(T, μ B ).
Introduction to the project: One of the Laboratory missions is to provide consistent, high-precision fission product; yield data critical for testing fission models and maintaining the safety and security of the nation’s nuclear weapons stockpile. However, high-quality, energy-differential fission product yield data is missing for certain actinides and neutron energies important to constrain the new U.S. Nuclear Data Program evaluation effort; It has been shown that the reactor antineutrino anomaly may be at least partially caused by roughly 20 fission products. The fission product yield data is missing or incomplete for many of these isotopes, thus it is necessary to accurately determine these values to better constrain the anomaly; Aitor Bracho is measuring very short-lived (seconds to minutes) fission product yields of 235 U and 239 Pu using monoenergetic neutron beams at E n = 60 and 560 keV; Aitor Bracho is using a direct approach utilizing a state-of-the-art rabbit transfer system, superior HPGe detector, and digital acquisition systems for fission decay measurements. The goals of this project: lop experimental capabilities and data analysis techniques to carry out the gamma-ray spectra analysis necessary for fission product yield calculation; Provide high-precision and energy-dependent fission product data supporting fission theory, neutrino physics, and applied physics.
We report the differential yields at mid-rapidity of the Breit-Wheeler process (𝛾𝛾 → 𝑒 + 𝑒 − ) in peripheral Au+Aucollisions at $\sqrt{s_{NN}}$ = 54.4 and 200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC), as a function of energy $\sqrt{s_{NN}}$, 𝑒 + 𝑒 − transverse momentum 𝑝 T , 𝑝$^{2}_{T}$, invariant mass 𝑀 𝑒𝑒 , and azimuthal angle. In the invariant mass range of 0.4 < 𝑀 𝑒𝑒 < 2.6GeV/𝑐 2 at low transverse momentum (𝑝 T < 0.15GeV/𝑐), the yields increase while the pair √⟨𝑝$^{2}_{T}$⟩ decreases with increasing $\sqrt{s_{NN}}$, a feature that is correctly predicted by the QED calculation. Here, the energy dependencies of the measured quantities are sensitive to the nuclear form factor, infrared divergence and photon polarization. The data are compiled and used to extract the charge radius of the Au nucleus.
An ionization spectrometer consisting of 12 layers of tungsten each 6.13 gm cm/2 thick and 7 iron modules each 70 gm cm/2 thick was flown in a balloon at an altitude of 6 gm cm/2 for 16 hours. The distributions of the first interaction points are used to determine the mean free path of the incident particles. At energies of 20 GeV the proton m.f.p. of 140 gm cm/2 is in agreement with that of other workers. The m.f.p.'s of heavy nuclei have been determined as a function of energy up to a total energy of 1000 GeV. These results are compared with the overlap model of nucleus-nucleus collisions. For C nuclei interacting with W the overlap parameter is -2.6 x 10/13 cm. The energy dependence of this parameter is discussed.
Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV γ-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to ∼100 GeV) and the Fermi Large Area Telescope to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current-sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second γ-ray light-curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the E field, as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.
The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.
Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.