Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gamma Distribution Function”

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 37 records · Page 2

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the Prompt Fission Gamma Properties for Neutron Induced Fission of 235,238 U and 239 Pu

Herein we present details of the prompt fission gamma property evaluations included in the released ENDF/B-VIII.0 files. The average prompt gamma multiplicity, the average prompt gamma spectrum, the average total prompt gamma energy released and the prompt gamma multiplicity distributions, included in the ENDF/B-VIII.0 files, are presented for 235 U(n,f), 238 U(n,f), and 239 Pu(n,f) reactions, as a function of incident neutron energy up to 30 MeV. The evaluation is based on available experimental data and model calculations. Notable data include recent measurements by experimental groups at Los Alamos and Lawrence Livermore National Laboratories, at Geel and Budapest, at CEA labs, as well as seminal historical gamma production measurements by Drake at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Search For Dark Matter Signals in the GalacticHalo with HAWC

The High Altitude Water Cherenkov Observatory (HAWC) is a high energy (500 GeV to 100 TeV) gamma ray detector located in central Mexico. HAWC has both a wide (≈ 2 sr) field of view and near continuous duty cycle, making it ideal for unbiased sky surveys and analysis of highly extended sources. We make use of HAWC's survey abilities to search for dark matter signals originating from the Milky Way Galaxy. By taking an unbiased sample of the sky, we set characteristic upper limits on dark matter annihilation flux for various dark matter masses and annihilation channels as a function of declination. These characteristic limits are then combined with models of the dark matter spatial distribution in order to calculate the HAWC sensitivity to gamma-ray emission from both dark matter substructure and the main Galactic halo. We use this sensitivity to optimize a search for dark matter emission from the main Galactic halo.

79 ASTRONOMY AND ASTROPHYSICS↗

INSIGHTS INTO MARINE UNICELLULAR CYANOBACTERIAL AND NON-CYANOBACTERIAL DIAZOTROPHS THROUGH SINGLE-CELL ANALYSES

Nutrient availability affects primary productivity, and nitrogen (N) is the limiting nutrient over large regions of the ocean. Biological N fixation (BNF) is the conversion of atmospheric dinitrogen (N2) to bioavailable forms of N (ammonium and amino acids). A select group of Bacteria and Archaea, known as diazotrophs, are the only organisms capable of this conversion. BNF is a critical source of N in nutrient-depleted (oligotrophic) marine waters. As such, BNF studies are a fundamental component for better understanding primary productivity and carbon cycling, as well as parameterizing predictive future global ocean models. Unicellular diazotrophs are widespread in marine systems and, in some cases, important contributors of BNF. This dissertation investigates unicellular diazotrophs at the single-cell level, including: UCYN-A, a cyanobacterial diazotroph that lives in symbiosis with an algal host, non-cyanobacterial diazotrophs (NCDs), which are putative heterotrophs, as well as a specific NCD known as Gamma A, a commonly occurring organism throughout the tropics and subtropics that belongs to the class Gammaproteobacteria. The single-cell investigations used in these studies can provide high-resolution insights that are not possible from community-level analyses of BNF, such as the N2 fixation rate of individual cells and cell morphology. Singlecell nanoscale secondary ion mass spectrometry (nanoSIMS) analyses were used in chapter 2 and 3 to investigate the cell-specific N2 fixation rates of (2) UCYN-A in the viii cold waters of the Bering and Arctic Seas where N2 fixation was not previously believed to be possible, and (3) NCDs in the North Pacific Subtropical Gyre where we demonstrated the first direct evidence of NCD N2 fixation on open ocean particles. Chapter 4 focused on the uncultivated NCD, Gamma A, using geneFISH to visualize single cells in natural seawater samples. We successfully identified Gamma A by developing a fluorescent in situ hybridization (FISH) probe targeting the diazotrophic functional marker gene nifH, revealing Gamma A cells are most often attached to particles. Together these studies provide new insights into the distributions, activities, and physiologies of geographically widespread, unicellular diazotrophs, thus improving our understanding of their roles as sources of fixed N in the surface ocean.

Harding, K.↗

Vectorization techniques for probability distribution functions using VecCore

Probability distribution functions (PDFs) are very used in modeling random processes and physics simulations. Improving the performance of algorithms that generate many random numbers under complex PDFs is often a very challenging task when methods as direct functions are not available. In this work we present general strategies on how to vectorize some PDFs using VecCore library. We show the results for the Exponential, Gaussian, discrete Poisson and Gamma probability distributions.

Chaparro Amaro, Oscar R.↗

Free-moving Quantitative Gamma-ray Imaging

The ability to map and estimate the activity of radiological source distributions in unknown three-dimensional environments has applications in the prevention and response to radiological accidents or threats as well as the enforcement and verification of international nuclear non-proliferation agreements. Such a capability requires well-characterized detector response functions, accurate time-dependent detector position and orientation data, a digitized representation of the surrounding 3D environment, and appropriate image reconstruction and uncertainty quantification methods. We have previously demonstrated 3D mapping of gamma-ray emitters with free-moving detector systems on a relative intensity scale using a technique called Scene Data Fusion (SDF). Here we characterize the detector response of a multi-element gamma-ray imaging system using experimentally benchmarked Monte Carlo simulations and perform 3D mapping on an absolute intensity scale. We present experimental reconstruction results from hand-carried and airborne measurements with point-like and distributed sources in known configurations, demonstrating quantitative SDF in complex 3D environments.

47 OTHER INSTRUMENTATION↗

Examination of the Radiolytic Stability of the GENIORS Solvent

The INL radiolysis test loop system has been employed to evaluate the effects of gamma radiolysis the stability and performance of the GENIORS solvent extraction process. This report describes details and results of the experiments performed and documents completion of milestone M3FT-21IN030101024 “Testing of Radiolytic Stability of m-TDDGA Based Solvent at INL”. The purpose of the proposed dynamic irradiation tests to be performed is to evaluate the radiolytic stability of the GENIORS process solvent for An/Ln recovery under process-like conditions. Recent research efforts within the European Union GENIORS project have focused upon demonstrating additional process improvements. These improvements include simplification of the solvent composition, use of reagents composed of only C, H, O, and N atoms (the CHON principle), identification of a diglycolamide derivative with more favorable Pu loading behavior, and use of an aqueous soluble stripping reagent that is not sulfonated. As a result, a new solvent utilizing dimethyl-N,N,N’,N’-tetradecyl-diglycolamide (cis-mTDDGA) as the An/Ln extractant and the aqueous soluble 2,6-bis[1-(propan-1-ol)-1,2,3-triazol-4yl]-pyridine (PTD) as the An/Ln stripping reagent was proposed. The impact of gamma radiolysis upon the efficacy of the GENIORS process solvent was evaluated by determining the americium, europium, and cerium distribution ratios as a function of absorbed dose. The measured distribution ratios demonstrated that the GENIORS solvent did undergo radiolytic degradation. However, at absorbed doses corresponding to realistic process conditions, the extraction and stripping performance of the solvent remains acceptable. Compositional analysis of the irradiated solvent revealed a decrease in the concentration of cis-mTDDGA, which is consistent with the observed changes in metal distribution ratios as a function of absorbed dose.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗

Relativistic Shear Boundary Layer and the Gamma-Ray Emission of GW170817

We present an analysis of the radiation characteristics of kinetic shear boundary layers created by relativistic plasma jets. Using a model of electromagnetic field data based on particle-in-cell simulations of an electron–ion plasma, we solve the motion of individual test electrons and compute their instantaneous radiated power and peak frequency. By analyzing a large number of test electrons in this manner, we find two distinct electron populations present around the shear boundary layer. The most highly radiative electrons execute looping motion due to crossed electric and magnetic fields as they are accelerated along the bulk flow of the jet, and eventually cross the shear boundary interface at steep angles. Electrons that never cross the shear boundary interface radiate much less energy as a group. Summing over all of the highly radiative electrons, we compute the distribution of the total radiated energy as a function of the angle relative to the bulk flow. Furthermore, this result has important potential implications for the observed radiation output of short gamma-ray bursts viewed at large angles from the jet axis, such as the neutron star merger event GW170817/GRB 170817A.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Approaches for Calculating and Aggregating Detection Probabilities for Nuclear Material Diversion

The authors built and tested a stochastic simulation to estimate achieved detection probabilities (DPs) on a stratum basis, over a tailorable range of diverted amounts from 0 to 2 significant quantities (SQ), using typical International Atomic Energy Agency (IAEA) inspection data: i.e., SQ in stratum, number of items, number of gross/partial/bias defect measurements conducted, and realistic relative standard deviation (RSD) values for typical IAEA verification measurements. For bulk strata, the model calculates achieved DP at 0.01 SQ diversion increments; for item strata, the model calculates DP using the smallest realistic diversion increment (e.g., a plate, pin, or coupon). After successfully benchmarking against IAEA deterministic models, the simulation was used to test the sensitivity of DP to certain standard assumptions and selected input parameters. First, the equal defect assumption was tested; the results suggest significant complexity in the effectiveness of partial defect measurements. Next, the authors explored the sensitivity of DP to the assumed RSD of attribute tests. Then, the authors compared non-normal models for instrument performance (e.g., logistic, step, or arbitrary functions) to the typical results from a normal distribution (characterized by RSD). This last comparison was supplemented with experimentally derived performance data for an HM-5 gamma spectrometer. The HM-5 was used to make enrichment measurements on both LEU and HEU MTR fuel elements as plates were removed, and the results fit with logistic and step curves and applied in the simulation. These stochastic DP results were compared to DP estimates from a deterministic model assuming a normal curve and typical RSD, yielding insights that could improve effectiveness in the field. These early results illustrate the potential of stochastic models to better understand achieved DP and to improve safeguards effectiveness.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Characterization of sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements with 3D hybrid particle-in-cell code

Hard x-rays produced by intense laser-produced fast electrons interacting with solids are a vital source for producing radiographs of high-density objects and implosion cores for inertial confinement fusion. Accurate calculation of hard x-ray sources requires a three-dimensional (3D) simulation geometry that fully models the electron transport dynamics, including electron recirculation and the generation of absolute photon yields. To date, 3D simulations of laser-produced bremsstrahlung photons over tens of picoseconds and code benchmarking have not been performed definitively. In this study, we characterize sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements for refluxing and non-refluxing targets using the 3D hybrid particle-in-cell (PIC), Large Scale Plasma code. Bremsstrahlung radiation and escaped electron data were obtained by focusing a 50-TW Leopard laser (15 J, 0.35 ps, 2 × 10 19 W/cm 2 ) on a 100-μm-thick Cu foil and a Cu with a large plastic backing (Cu–CH target). Data for both the Cu and Cu–CH targets were reproduced for simulations with a given set of electron parameters. Comparison of the simulations revealed that the hard x-ray emission from the Cu target was significantly longer in duration than that from the Cu–CH target. Here, the benchmarked hybrid PIC code could prove to be a powerful tool in the design and optimization of time- and angular-dependent bremsstrahlung sources for flash x-ray and gamma-ray radiography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extended Interacting Boson Model with Configuration Mixing Description For 98 Mo

An extended version of the interacting boson model (IBM), encompassing two-particle–two-hole (2p-2h) excitations originating from below a valence shell and one-particle–one-hole (1p-1h) excitations extending to the upper next-nearest major shell, is developed to investigate cross-shell excitation effects in the transitional nucleus 98 Mo. This extended IBM framework is utilized to reproduce experimentally observed properties such as low-lying positive-parity level energies, 𝐵⁡(𝐸⁢2) values, and the electric quadrupole moment of the 2$^{+}_{1}$ state, 𝜌 2 ⁡(𝐸⁢0, 0$^{+}_{2}$ → 0$^{+}_{g}$), as well as the strength functions of the isoscalar giant monopole resonance (ISGMR) and the isoscalar giant quadrupole resonance (ISGQR) in 98 Mo. It is shown that the extended consistent-𝑄 IBM, incorporating configuration mixing of both 2p-2h excitations from below the valence shell and 1p-1h excitations to the upper next-nearest major shell, describes the low-lying structural properties and the ISGMR and ISGQR strength distributions of 98 Mo quite well, with the ISGMR and the ISGQR strength distributions comparable to the results of the self-consistent quasiparticle random-phase approximation with Skyrme interactions.

90 ≤ A ≤ 149↗

ResNet and CycleGAN for pulse shape discrimination of He-4 detector pulses: Recovering pulses conventional algorithms fail to label unanimously

Pulse shape discrimination (PSD) capable detectors, such as He-4, that respond to neutron and gamma-ray 7 interactions have a threshold deposited energy value below which n/γ discrimination vanishes when using 8 conventional PSD algorithms. Recent attempts in applying supervised learning based artificial neural 9 networks for PSD use the pulses in the separated regions to train the networks so they can be used to classify 10 another set of separated pulses. In doing so, pulses previously indistinguishable are not recovered for 11 classification, which would have increased the number of neutron and gamma-ray pulses that could be used 12 for further analysis. Assuming the reason why conventional PSD algorithms have unseparated regions is 13 because the parameter space of the algorithms fail to capture the intrinsic (but subtle) distinguishing 14 behavior of some of the neutron and gamma-ray pulses, a cycle-consistent generative adversarial network 15 (CycleGAN) was trained to amplify those differences and extract well separated neutron and gamma-ray 16 clusters. Results show that, once the network is trained with pulses from separated and unseparated regions, 17 it was able to transform the pulses in the unseparated region to improve the PSD. Subsequent n/γ 18 classification was performed using deep residual network (ResNet) that takes pulses with 512 data points 19 as an input. Two different ResNets were explored – simple ResNet and modified ResNet which takes 20 segmented pulse inputs in the first layer and the corresponding time axis values in the last hidden layer. 21 The later approach enables the network to extract time correlated pulse features to enhance its ability to 22 capture the pulse behaviors relevant for PSD. Although it achieves slightly lower accuracy, 99.41% versus 23 99.89%, based on simply counting the number of correct n/γ labels assigned, compared to the simple 24 ResNet, the modified ResNets architecture was able to decreases the cross-entropy loss function by half, 25 which implies that the correct n/γ labels assigned are less likely to be accidental. PSD parameter 26 distributions based on n/γ classification by ResNet before and after transforming unseparated pulses using 27 CycleGAN show that by enhancing the separation between neutrons and gamma-rays, the transformation 28 helps improve the performance of classifier networks that are trained using labeled dataset. The 29 enhancement of neutron and gamma-ray separation by the CycleGAN increased the PSD figure of merit 30 (FOM) by up to 70% in some regions. Here, the results show that, if a given detector achieves clear separation 31 between neutron and gamma-ray pulses in any energy region, such neural network approaches can help 32 lower the energy threshold for the separation and increasing the number of neutron and gamma-ray pulses 33 that can be used for further analysis.

4He↗

Extracting the gamma-ray source-count distribution below the Fermi-LAT detection limit with deep learning

We reconstruct the extra-galactic gamma-ray source-count distribution, or dN/dS, of resolved and unresolved sources by adopting machine learning techniques. Specifically, we train a convolutional neural network on synthetic 2-dimensional sky-maps, which are built by varying parameters of underlying source-counts models and incorporate the Fermi-LAT instrumental response functions. The trained neural network is then applied to the Fermi-LAT data, from which we estimate the source count distribution down to flux levels a factor of 50 below the Fermi-LAT threshold. We perform our analysis using 14 years of data collected in the (1,10) GeV energy range. The results we obtain show a source count distribution which, in the resolved regime, is in excellent agreement with the one derived from cataloged sources, and then extends as dN/dS ~ S -2 in the unresolved regime, down to fluxes of 5 · 10 -12 cm -2 s -1 . The neural network architecture and the devised methodology have the flexibility to enable future analyses to study the energy dependence of the source-count distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of the differential cross sections of the production of Z + jets and $\gamma$ + jets and of Z boson emission collinear with a jet in pp collisions at $ \sqrt{s} $ = 13 TeV

Measurements of the differential cross sections of Z + jets and γ + jets production, and their ratio, are presented as a function of the boson transverse momentum. Measurements are also presented of the angular distribution between the Z boson and the closest jet. The analysis is based on pp collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 35.9 fb$^{−1}$ recorded by the CMS experiment at the LHC. The results, corrected for detector effects, are compared with various theoretical predictions. In general, the predictions at higher orders in perturbation theory show better agreement with the measurements. This work provides the first measurement of the ratio of the differential cross sections of Z + jets and γ + jets production at 13 TeV, as well as the first direct measurement of Z bosons emitted collinearly with a jet.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Flat-spectrum Radio Quasars and BL Lacs Dominate the Anisotropy of the Unresolved Gamma-Ray Background

We analyze the angular power spectrum (APS) of the unresolved gamma-ray background (UGRB) emission and combine it with the measured properties of the resolved gamma-ray sources of the Fermi-LAT 4FGL catalog. Our goals are to dissect the composition of the gamma-ray sky and to establish the relevance of different classes of source populations of active galactic nuclei in determining the observed size of the UGRB anisotropy, especially at low energies. We find that, under physical assumptions for the spectral energy distribution, i.e., by using the 4FGL catalog data as a prior, two populations are required to fit the APS data, namely flat-spectrum radio quasars at low energies and BL Lacs at higher energies. The inferred luminosity functions agree well with the extrapolation of the flat-spectrum radio quasar and BL Lac ones obtained from the 4FLG catalog. We use these luminosity functions to calculate the UGRB intensity from blazars, finding a contribution of 20% at 1 GeV and 30% above 10 GeV. Finally, bounds on an additional gamma-ray emission due to annihilating dark matter are also derived.

79 ASTRONOMY AND ASTROPHYSICS↗