Search NASA⌕ Search

SEARCH · Search NASA

Results for “reach code”

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 541 records · Page 30

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗

Spreader-Bar Radiation Detection System Enhancements: A Modeling and Simulation Study

An investigation of radiation detectors on crane-mounted platforms was performed to understand the limitations and optimal capability of these types of systems. Crane-mounted systems are primarily focused on crane-mounted radiation detector systems, where the detectors are mounted on a spreader-bar platform. Spreader bars are attached with cables to a crane at seaports to load and unload containers from cargo ships. The bars can spread from 20 feet to 53 feet using hydraulics, allowing a single bar to move a variety of container lengths. Spreader-bars are also used on straddle carriers, which move containers around the laydown yard in the port. Spreader bar radiation detection (SBRD) systems are commercially available from at least two companies. These systems employ a set of spectroscopic sensors placed on the underside of a gantry crane and held near the top of a cargo container for detection. These systems have been previously tested and some issues and limitations observed. Radiation detectors mounted on spreader-bars are single sided systems (observing cargo only from the top) and it is therefore challenging to obtain the same sensitivity as double-sided systems (for example, portal monitors) especially when considering the possible amounts of intervening cargo that could be present. The approximately exponential attenuation of gamma rays in intervening dense cargo presents an inherent limitation to detection. However, SBRD systems hold detectors closer to the container, have a longer dwell time, and may have lower background (over the ship); all of these aspects help mitigate the single-sided issue. In addition, there are operational considerations that make spreader-bar radiation systems attractive. These include scanning for radiation during other operations (loading, unloading or container movement) which minimizes additional time for scanning, incorporating the radiation detection into existing equipment to minimize initial cost, and in the case of ship unloading, scanning the containers before reaching the laydown yard. This study is focused on researching optimized configurations of SBRD systems, including number of detectors, locations, and best use of detector data, to determine a reasonable limit on the capability of these types of systems. The results will help support decisions on where and how these systems might be effectively deployed, and could be used for future specification development for SBRD acquisition. The next phase of this research and investigation into SBRD systems and quantification of the possible improved capability with an enhanced SBRD is the focus of this report. Supporting this investigation, detailed models of an enhanced SBRD system were created using a 3D radiation transport code, and models were run to simulate the data from sources and background for a realistic cargo-loading port environment. The enhanced SBRD had increased detector material (numbers of detectors) compared to the commercial systems, with a distribution providing more uniform coverage. Measurements using small industrial sources were made with an enhanced SBRD system developed by Sandia National Laboratory to support model validation to establish credibility for the simulations of these threat and background scenarios. A significant simulation effort was undertaken to simulate sources of interest within a 40-foot cargo container. In order to estimate the performance of a generalized SBRD system in scenarios including background and cargo, a highly enriched uranium (HEU) test source (with 3x10 -8 232 U by weight) within a wood cargo model was simulated for a 40-foot cargo container. This source was placed at a variety of locations, and the flux from the source was multiplied as necessary to create estimated minimal detectable quantities (MDQ) figures for various quantities of HEU. Simulated data sets were produced corresponding to a complete container offloading scenario for a large number of statistical realizations. Significant work was performed on both the algorithm development and the container cargo distributions entering the US, and the result of this modeling, validation, and simulation effort was a set of estimates of the performance of a variety of radiation detection algorithms in this cargo-screening environment. In the case of algorithms, a variety of contemporary radiation detection and identification approaches were employed, and the optimal approaches to selecting time and location subsets of detector data were explored. In all cases, algorithm performance was compared, and the best approaches were suggested along with the best expected performance from the enhanced SBRD. In general, template-matching schemes employing a limited set of regression parameters tend to provide best results. The improvements of an enhanced SBRD with increased detector material and improved algorithms result in an effective MDQ of 10 kg of HEU in ~50 percent of the imported cargo in 40-foot containers, assuming a random source location. This is approximately two times the capability of the estimated performance of the commercial systems. For an MDQ of the 10 kg HEU source located in the worst location (bottom of the 40-foot container), the volume of cargo scanned drops to 10% for the enhanced SBRD, again a factor of two or more larger than the commercial systems. Enhancing the SBRD systems by adding more detectors, providing a more uniform coverage, and incorporating more advanced algorithms will result in a significant increase in capability. Although this investigation provided some indication of the performance capabilities of an enhanced SBRD system, additional work is needed to estimate the capability for other container configurations, such as the 20 and 45 foot containers. Since the 20-foot container weight limit is similar to the 40-foot, more dense cargo is typically shipped in 20-foot containers. This will result in a smaller percentage of cargo scanned for the 20-foot containers (for a particular source quantity) compared to the 40-foot containers. In addition to other configuration, other source material needs to be investigated to understand the capabilities for the entire threat space, where only HEU is considered in this work.

crane-mounted↗

Triggering Collapse of the Presolar Dense Cloud Core and Injecting Short-lived Radioisotopes with a Shock Wave. V. Nonisothermal Collapse Regime

Recent meteoritical analyses support an initial abundance of the short-lived radioisotope (SLRI) {sup 60}Fe that may be high enough to require nucleosynthesis in a core-collapse supernova, followed by rapid incorporation into primitive meteoritical components, rather than a scenario where such isotopes were inherited from a well-mixed region of a giant molecular cloud polluted by a variety of supernovae remnants and massive star winds. This paper continues to explore the former scenario, by calculating three-dimensional, adaptive mesh refinement, hydrodynamical code (FLASH 2.5) models of the self-gravitational, dynamical collapse of a molecular cloud core that has been struck by a thin shock front with a speed of 40 km s{sup −1}, leading to the injection of shock front matter into the collapsing cloud through the formation of Rayleigh–Taylor fingers at the shock–cloud intersection. These models extend the previous work into the nonisothermal collapse regime using a polytropic approximation to represent compressional heating in the optically thick protostar. The models show that the injection efficiencies of shock front materials are enhanced compared to previous models, which were not carried into the nonisothermal regime, and so did not reach such high densities. The new models, combined with the recent estimates of initial {sup 60}Fe abundances, imply that the supernova triggering and injection scenario remains a plausible explanation for the origin of the SLRIs involved in the formation of our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

The HUNTER Dynamic Human Reliability Analysis Tool: Coupling an External Plant Code

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a framework to support dynamic human reliability analysis (HRA) with the aim to develop standalone software to perform the dynamic HRA calculations. Within the HRA, human actions in nuclear power plants (NPPs) are predicated by plant states, and human actions influence the plant. In other words, plant operations are necessarily recursive, and it becomes challenging to model complex human-plant interactions. Consequently, we have linked two software simulations that complement those shortcomings. RELAP5-3D—the Reactor Excursion and Leak Analysis Program (RELAP; Aumiller, Tomlinson, and Bauer 2001) is the foundational thermal-hydraulic software used to model nuclear systems. Using RELAP5-3D, we have simulated the plant operations proceeding according to procedures developed to address emergent situations in NPPs. Plant operations include various actions such as the operator checking plant parameters, as well as actions that are continuously performed over time until a specific parameter reaches certain criteria. That means that HUNTER and RELAP5-3D exchange information with each other and should be carried out simultaneously over time. To simulate plant operations, which represent the actual operator checks of plant parameters and corresponding manual control actions, changes in plant status are identified through simulation and performed according to the criteria and order of the procedure. Thus, the goal of coupling HUNTER with RELAP5-3D is to facilitate synchronous coupling, where human and plant models provide iterative feedback loops that drive the course of actions. The advantage of coupling with RELAP5-3D to serve as the external environment module in HUNTER is the ability to customize the plant model and streamline for particular model applications. In this paper, we will address the key features of the coupling and the coupling structure built to perform the feedback loops.

99 GENERAL AND MISCELLANEOUS↗

A WISER way for simulating the performance of gratings

Soft X-ray monochromators for synchrotron radiation sources have been continuously developed over the years, improving energy resolution and stability. Great effort has been made in improving the surface quality of the optics involved, reaching values permitting diffraction-limited images. Still, one problem has not been solved, nor fully understood, yet: groove placing errors. Nowadays, these are one of the major factors responsible for resolving the power reduction of diffraction-gratings-based X-ray monochromators. Despite decades of use of gratings, there is not yet a universally established method for predicting and simulating the effect of groove placing errors on monochromator performance. This is especially important in the new scenario of high-coherent X-ray sources, i.e. diffraction-limited storage rings and free-electron lasers. To address this problem, in this article an approach based on WISER (Wavefront propagatIon Simulation codE libRary) is presented. WISER is a physical optics simulation package, also available in the user-friendly Orange Synchrotron Radiation Suite – OASYS. Even though it was originally conceived to assess the focusing performance of X-ray mirrors in the presence of height defects, it perfectly simulates the performance of a periodic (or quasi-periodic) structure like a diffraction grating. In this article, the way to use WISER and its application to a specific case, e.g. the design of a monochromator for the upgrade of a beamline at the Advanced Light Source, are shown. A simple rule for estimating how well the grooves are placed on a grating, based on calculation of the Strehl ratio, is also presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Bonded Bracket Assmebly for Frameless Solar Panels

In February 2011 the US Department of Energy announced their new Sunshot Initiative. The Sunshot goal is to reduce the total cost of solar energy systems by about 75 percent before the end of the decade. The DOE estimated that a total installed cost of $$1 per watt for photovoltaic systems would be equivalent to 5-6¢/kilowatt hour (kWh) for energy available from the grid. The DOE also estimated that to meet the $1 per watt goal, PV module costs would need to be reduced to $ .50 per watt, balance of systems costs would need to be reduced to $.40 per watt, and power electronic costs would need to reach $.10 per watt. To address the BOS balance of systems cost component of the $1 per watt goal, the DOE announced a funding opportunity called (BOS-X) Extreme Balance of System Hardware Cost Reductions. The DOE identified eight areas within the total BOS costs: 1) installation labor, 2) installation materials, 3) installation overhead and profit, 4) tracker, 5) permitting and commissioning, 6) site preparation, 7) land acquisition, 8) sales tax. The BOS-X funding announcement requested applications in four specific topics;Topic 1: Transformational Building Integrated Photovoltaic (BIPV) Modules; Topic 2: Roof and Ground Mount Innovations; Topic 3: Transformational Photovoltaic System Designs; and Topic 4: Development of New Wind Load Codes for PV Systems.The application submitted by ARaymond Tinnerman reflected the requirements listed in Topic #2, Roof and Ground Mount Innovations. The goal of topic #2 was to develop technologies that would result in the extreme reduction of material and labor costs associated with applications that require physical connections and attachments to roof and ground mount structures. The topics researched in this project included component cost reduction, labor reduction, weight reduction, wiring innovations, and alternative material utilization. The project objectives included; 1) The development of an innovative quick snap bracket assembly that would be bonded to frameless PV modules for commercial rooftop installations; 2) The development of a composite pultruded rail to replace traditional racking materials; 3) In partnership with a roofing company, pilot the certification of a commercial roof to be solar panel compliant, eliminating the need for structural analysis and government oversight resulting in significantly decreased permitting costs; and 4) Reduce the sum of all cost impacts in topic #2 from a baseline total of $2.05/watt to $.34/watt.

14 SOLAR ENERGY↗

Study on Water-ingress Accident during the Initial Core and the Running-in Phase of HTR-PM - 18429

Water ingress accident is one of the typical accidents with potential hazard to the high temperature gas-cooled reactor (HTGR), which will cause reactivity introduction, as well as the chemical corrosion of the graphite fuel elements and reflector structure material. Besides, increase of the primary pressure may result in the opening of the safety valves, consequently leading the release of radioactive isotopes. Therefore high attention has been paid to water ingress phenomenon during the design of the HTGR to verify its inherent safety characteristics. The online continuous discharge and multi-pass of the fuel elements is another important characteristics for the pebble-bed HTGR which can effectively increase the reactor operation availability. From the initial core, a long period of time, namely the running-in phase, is needed to reach the final long-term stable equilibrium core. During this period of time, the core states are totally different, including the number of the fuel elements, the fuel enrichment, the accumulation of the fission product, the temperature feedback coefficient of the reactivity, and so on. In this paper, based on the design of 250 MW Pebble-bed Module High Temperature gas-cooled Reactor (HTR-PM) and a candidate proposal for its initial core and running-in phase, considering three typical states during the running-in phase, the design basis accident of a double-ended guillotine break of one heating tube has been analyzed with the TINTE code, and the results are compared with that of the equilibrium core. Some important parameters during the transient process, including the power, the temperature of the fuel element, the pressure of the primary coolant, as well as the graphite corrosion rate, are studied. To a certain extent, the results can verity the design of the initial core and the running-in phase, and also further prove the safety feature of the HTR. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Impact of the uncertainty in α-captures on {sup 22}Ne on the weak s-process in massive stars

Massive stars at solar metallicity contribute to the production of heavy elements with atomic masses between A = 60 and A = 90 via the so-called weak s-process (which takes place during core He and shell C burning phases). Furthermore, recent studies have shown that rotation boosts the s-process production in massive stars at low metallicities, with a production that may reach the barium neutron-magic peak. These results are very sensitive to neutron source and neutron poison reaction rates. For the weak s-process, the main neutron source is the reaction {sup 22}Ne(α,n){sup 25}Mg, which is in competition with {sup 22}Ne(α,γ){sup 26}Mg. The uncertainty of both rates strongly affects the nucleosynthesis predictions from stellar model calculations. In this study, we investigate the impact of the uncertainty in α-captures on {sup 22}Ne on the s-process nucleosynthesis in massive stars both at solar and at very low metallicity. For this purpose, we post-process, with the Nugrid mppnp code, non-rotating and rotating evolutionary models 25M{sub ⊙} stars at two different metallicities: Z = Z{sub ⊙} and Z = 10{sup −5}Z{sub ⊙}, respectively. Our results show that uncertainty of {sup 22}Ne(α,n){sup 25}Mg and {sup 22}Ne(α,γ){sup 26}Mg rates have a significant impact on the final elemental production especially for metal poor rotating models. Beside uncertainties in the neutron source reactions, for fast rotating massive stars at low metallicity we revisit the impact of the neutron poisoning effect by the reaction chain {sup 16}O(n,γ){sup 17}O(α,γ){sup 21}Ne, in competition with the {sup 17}O(α,n){sup 20}Ne, recycling the neutrons captured by {sup 16}O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collective mode Brownian dynamics: A method for fast relaxation of statistical ensembles

Sampling equilibrium configurations of correlated systems of particles with long relaxation times (e.g., polymeric solutions) using conventional molecular dynamics and Monte Carlo methods can be challenging. This is especially true for systems with complicated, extended bond network topologies and other interactions that make the use and design of specialized relaxation protocols infeasible. We introduce a method based on Brownian dynamics simulations that can reduce the computational time it takes to reach equilibrium and draw decorrelated samples. Importantly, the method is completely agnostic to the particle configuration and the specifics of interparticle forces. In particular, we develop a mobility matrix that excites non-local, collective motion of N particles and can be computed efficiently in $\mathcal{O}(N)$ time. Particle motion in this scheme is computed by integrating the overdamped Langevin equation with an Euler–Maruyama scheme, in which Brownian displacements are drawn efficiently using a low-rank representation of the mobility matrix in position and wave space. We demonstrate the efficacy of the method with various examples from the realm of soft condensed matter and release a massively parallel implementation of the code as a plugin for the open-source package HOOMD-blue [J. A. Anderson et al., J. Comput. Phys. 227, 5342 (2008) and J. Glaser et al., Comput. Phys. Commun. 192, 97 (2015)] which runs on graphics processing units

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamical Aspects of the Quark Gluon Plasma

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of tapering structures on the characteristics of a coaxial-waveguide gyrotron backward-wave oscillator

This study analyzes the characteristics of a gyrotron backward-wave oscillator (gyro-BWO) with a longitudinally tapered coaxial-waveguide by using a single-mode, self-consistent nonlinear code. Simulation results indicate that although tapering the inner wall or the outer wall can significantly raise the start-oscillation current, the former is more suitable for mode selection than the latter because an increase of the start-oscillation current by a tapered inner wall heavily depends on the chosen C value (i.e., the average ratio of the outer radius to the inner radius over the axial waveguide length). Selective suppression of the competing mode by tapering the inner wall is numerically demonstrated. Moreover, efficiency of the coaxial gyro-BWO is increased by tapering the outer wall. Properly down-tapering the outer wall ensures that the coaxial gyro-BWO can reach a maximum efficiency over twice that with a uniform one.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

CHARACTERIZING EXOPLANETARY ATMOSPHERES THROUGH INFRARED POLARIMETRY

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties, and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g., due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1%, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line of sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes.

79 ASTRONOMY AND ASTROPHYSICS↗

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗

Progress in the Measurement of Prompt Neutron Output in Neutron-Induced Fission of 239 Pu: The Chi-Nu Project

Knowledge of the prompt neutron output from neutron-induced fission of 239 Pu is essential to nuclear defense and energy applications. Because the experimental data base is inadequate to constrain models of the spectrum of neutrons from fission, we have embarked on a multi-year program of new measurements of the prompt fission neutron spectrum as a function of the energy of the incident neutron that induces the fission. This is a joint LANL-LLNL project with participation by four university groups in advanced detector development. In the past year much progress has been achieved. The building housing the experiment has been completed and the flight paths including one for this experiment are now constructed. A parallel-plate avalanche counter (PPAC) containing 235 U is installed and working in the beam line. Neutron detectors for both low (<1 MeV) and high (> 0.5 MeV) neutrons have been assembled. The data acquisition system based on waveform digitizers has undergone extensive development and now can handle 10 PPAC channels and 10 neutron detectors simultaneously. In addition, a detector test station centered on a 252 Cf spontaneous fission chamber has been commissioned so that we can characterize our neutron detectors in a well understood experimental environment. Analysis of experiments at this facility has helped us understand the performance of several types of neutron detectors. In addition, the analysis work uncovered that the “standard” 252 Cf fission neutron spectrum in MCNPX disagrees significantly from that presently accepted by the international standards community. One year ago we submitted a report (LA-UR-11-05576) on the accuracies that could be reached under certain assumptions. Several sources of possible systematic uncertainties were identified, and procedures to control these uncertainties were suggested. This report is included as an appendix to this report. One element missing from these estimates was the contribution of backgrounds. It could not be quantified because, at that time, the building and flight path for the experiments had not been constructed. We now have made some measurements of the background and can now add this component to the accuracy projections. This report summarizes the progress in the past year. It also discusses measurements reported in the literature and some factors that should be considered when using these data to optimize parameters in nuclear model codes used in constructing evaluated nuclear data bases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

IAEA CRP on HTGR Uncertainties in Modeling: Assessment of Phase I Lattice to Core Model Uncertainties

Best-estimate plus uncertainty analysis of reactors is replacing the traditional conservative (stacked uncertainty) method for safety and licensing analysis. To facilitate uncertainty analysis applications, a comprehensive approach and methodology must be developed and applied. High temperature gas cooled reactors (HTGRs) have several features that require techniques not used in light-water reactor analysis (e.g., coated-particle design and large graphite quantities at high temperatures). The International Atomic Energy Agency has therefore launched the Coordinated Research Project on HTGR Uncertainty Analysis in Modeling to study uncertainty propagation in the HTGR analysis chain. The benchmark problem defined for the prismatic design is represented by the General Atomics Modular HTGR 350. The main focus of this report is the compilation and discussion of the results obtained for various permutations of Exercise I 2c and the use of the cross section data in Exercise II 1a of the prismatic benchmark, which is defined as the last and first steps of the lattice and core simulation phases, respectively. The report summarizes the Idaho National Laboratory (INL) best estimate results obtained for Exercise I 2a (fresh single-fuel block), Exercise I 2b (depleted single-fuel block), and Exercise I 2c (super cell) in addition to the first results of an investigation into the cross section generation effects for the super-cell problem. The two dimensional deterministic code known as the New ESC based Weighting Transport (NEWT) included in the Standardized Computer Analyses for Licensing Evaluation (SCALE) 6.1.2 package was used for the cross section evaluation, and the results obtained were compared to the three dimensional stochastic SCALE module KENO VI. The NEWT cross section libraries were generated for several permutations of the current benchmark super-cell geometry and were then provided as input to the Phase II core calculation of the stand alone neutronics Exercise II 1a. The steady state core calculations were simulated with the INL coupled-code system known as the Parallel and Highly Innovative Simulation for INL Code System (PHISICS) and the system thermal-hydraulics code known as the Reactor Excursion and Leak Analysis Program (RELAP) 5 3D using the nuclear data libraries previously generated with NEWT. It was observed that significant differences in terms of multiplication factor and neutron flux exist between the various permutations of the Phase I super-cell lattice calculations. The use of these cross section libraries only leads to minor changes in the Phase II core simulation results for fresh fuel but shows significantly larger discrepancies for spent fuel cores. Furthermore, large incongruities were found between the SCALE NEWT and KENO VI results for the super cells, and while some trends could be identified, a final conclusion on this issue could not yet be reached. This report will be revised in mid 2016 with more detailed analyses of the super-cell problems and their effects on the core models, using the latest version of SCALE (6.2). The super-cell models seem to show substantial improvements in terms of neutron flux as compared to single-block models, particularly at thermal energies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Transfer Request from Rutgers University to University of Illinois Urbana Champaign

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons will be exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges must be developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients will be calculated across beam energies. New techniques will be developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project will provide essential guidance to the upcoming Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report on Assessment of Advanced Reactor Fuel Performance

Three activities have been performed in FY 2019 to support the Advanced Fuel Cycle Campaign, which include assessment of annular fuel performance in an ultra-high burnup sodium-cooled reactor (SFR), modeling of annular fuel behavior during a transient scenario, and feasibility of a sustainable closed fuel cycle utilizing melt-refining process. Annular fuel performance parameters, including temperature profile, swelling strain, and cumulative damage fraction (CDF), in a 3000 MWt Sustainable Sodium-cooled Fast Reactor (SSFR) were assessed using the BISON fuel performance code. The SSFR is once-through ultra-long-life core based on a 34-batch fuel management scheme with 1.5-year cycle length. The fuel residence time in the SSFR is 51 years and the average discharged burnup is about 30%. Annular fuel has potential to irradiate for ultra-high burnup without major fuel-cladding mechanical interaction. This feature is achieved by using a low smeared density (~55%). In the SSFR design, however, 75% smeared density and a thick fuel pin were preferred to maintain criticality for 51 years. Due to the high smeared density and thick fuel pin configuration, the annular fuel was predicted to suffer premature cladding failure before reaching the discharge burnup of the SSFR. The current fuel creep model predicts a prominent creep strain due to gravity over the half-century irradiation, which may significantly affect the evolution of CDF. Thus, in order to use the annular fuel in an ultra-high burnup core such as the SSFR, a low smeared density and thin annular fuel is strongly recommended.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SU-E-T-148: Efficient Verification Method for Modulated Electron Radiotherapy Treatment Plans

Purpose: For shallow tumors, modulated electron radiotherapy (MERT) promises a reduction of dose to distal organs at risk. At our institution a framework was developed in order to create treatment plans for MERT employing both forward and inverse optimization. In this work, an efficient quality assurance (QA) process is established. Methods: Treatment plans for three different tumor sites were created using an inverse optimization. These plans consist of 6–12 segments and energies between 6 and 18 MeV. An already established QA process for photon IMRT plans is now extended to additionally handle MERT plans. First, the dose distributions are calculated in a homogenous water phantom. For this task a dedicated Monte Carlo (MC) framework for MERT is used. Second, the segments are applied on a stand-alone amorphous silicon electronic portal imaging device (EPID) using a source-to-surface distance of 70 cm. This device was calibrated for electron beams in a previous work. An in-house developed analysis software, is then utilized for comparisons and evaluation of the measured and calculated dose distributions. Results: For all three plans the calculated dose distributions agree well with the measured ones. Using a 2D gamma comparison (2% of dose max/2 mm and 10% dose threshold) passing rates >98% are achieved.The dose calculation for each plan on the water phantom, using voxels of 0.2×0.2×0.2 cm{sup 3}, takes at maximum 30 min on a single core Pentium 2.66 GHz system with 6 GB RAM, to reach a statistical uncertainty of 2% (1 std. dev.). Conclusion: An already established QA procedure for IMRT photon plans was applied for MERT. The dedicated MC framework and the use of EPID measurements allow an efficient QA procedure in a clinical environment. This work was supported by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗