Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate dynamic programming”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Improved scaling of temperature-accelerated dynamics using localization

While temperature-accelerated dynamics (TAD) is a powerful method for carrying out non-equilibrium simulations of systems over extended time scales, the computational cost of serial TAD increases approximately as N{sup 3} where N is the number of atoms. In addition, although a parallel TAD method based on domain decomposition [Y. Shim et al., Phys. Rev. B 76, 205439 (2007)] has been shown to provide significantly improved scaling, the dynamics in such an approach is only approximate while the size of activated events is limited by the spatial decomposition size. Accordingly, it is of interest to develop methods to improve the scaling of serial TAD. As a first step in understanding the factors which determine the scaling behavior, we first present results for the overall scaling of serial TAD and its components, which were obtained from simulations of Ag/Ag(100) growth and Ag/Ag(100) annealing, and compare with theoretical predictions. We then discuss two methods based on localization which may be used to address two of the primary “bottlenecks” to the scaling of serial TAD with system size. By implementing both of these methods, we find that for intermediate system-sizes, the scaling is improved by almost a factor of N{sup 1/2}. Some additional possible methods to improve the scaling of TAD are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Energy & Environmental Research Center (EERC)–U.S. Department of Energy (DOE) Joint Program on Research and Development for Fossil Energy-Related Resources

The EERC teamed with the DOE National Energy Technology Laboratory to conduct complementary research and development under the Joint Program on Research and Development for Fossil Energy-Related Resources (Cooperative Agreement No. DE-FC26-08NT43291). This 10-year program merged two previous 10-year agreements—DE-FC26-98FT40320 and DE-FC26-98FT40321—into a single follow-on agreement. The Fossil Energy cooperative agreement between DOE and the EERC was initiated in 1983 when the EERC was defederalized. This dynamic 35-year partnership has evolved to further research in areas that address DOE Office of Fossil Energy program goals while serving the missions and exploiting the strengths of both organizations through this agreement and a new 5-year Fossil Energy-Related Resources cooperative agreement established in 2015. As the 10-year program has ended, the purpose of this report is to summarize the activity and significant outcomes that were accomplished by the EERC. This program is an excellent model for research, development, demonstration, and commercialization partnerships between government, industry, and the applied science and engineering communities to bring cutting-edge science closer to commercial application. The overarching goal was to support DOE Fossil Energy goals by advancing the scientific knowledge and technical development essential to ensuring future sustainable supplies of affordable energy and clean water, protecting and restoring the environment, and reducing dependence on foreign energy sources, thereby increasing U.S. energy security. The strategic objective was to advance continued use of domestic fossil fuels as a mainstay of U.S. energy production by making fossil energy systems nonpolluting and more efficient, capturing and sequestering greenhouse gases, and integrating the use of fossil and renewable energy sources into the energy mix. The research achievements from this partnership have exceeded the goals and objectives by contributing substantially to the development of science-based energy and environmental policy, educational foundations for science and technology, and commercialization and international marketing of energy technologies. The program comprised 75 individual projects, which included basic/applied fossil energy research projects involving no cost share and development, demonstration, and commercialization projects (many derived from applied research projects) that included nonfederal cost share of approximately 40% (the contractually required minimum level of cumulative nonfederal cost share was 29%). In all, 32,798,698 of federal funding was leveraged with industry cash and in-kind cost share of 21,697,129 from 80 nonfederal partners for a total of 54,495,827 of funded research that 1) advanced the development of unconventional tight oil resources in the Bakken Formation; 2) demonstrated the application of CO 2 Storage and EOR in tight oil formations; 3) developed advanced carbon capture technology; 4) addressed key application and infrastructure concerns for commercial deployment of CO 2 ; 5) demonstrated innovative, integrated water management strategies for energy development projects; 6) developed novel pollution control technologies for mercury and hazardous air pollutants; and 7) performed world-class outreach and dissemination of research results to industry partners and other stakeholders. These outcomes speak to the valuable return on the substantial investment in this program.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2010

QM02 Strength Measurement

In late April, Paul Emma reported that his orbit fitting program could find a reasonably good fit only if the strength of QM02 was changed from design value of -5.83 kG to -6.25 kG - a strength change of 7.3%. In late May, we made a focal length measurement of QM02 by turning off all focusing optics between YC07 and BPMS1 (in the spectrometer line) except for QM02 and adjusted the strength of QM02 so that vertical kicks by YC07 did not produce any displacements at BPMS1 (see Figure 1). The result was quoted in the LCLS elog was that QM02 appeared to 6% too weak, and approximately agreed with Paul's observation. The analysis used for the entry in the log book was based on the thin lens approximation and used the following numbers: Distance YC07 to QM02 - 5.128 m; Distance QM02 to BPMS1 - 1.778 m; and Energy - 135 MeV. These distances were computed from the X,Z coordinates given the on the large plot of the Injector on the wall of the control room. On review of the MAD output file coordinates, it seems that the distance used for QM02 to BPMS1 is not 1.778 m. The correct value is Distance, center of QM02 to BPMS1 - 1.845 m. There may be a typo on the wall chart values for the coordinates of BPMS1, or perhaps there was a misinterpretation of edge versus center of QM02. In any case, the effect of this change is that the thin lens estimate changes from 6% too weak to 9% too weak. At John Galayda's suggestion, we looked into the thin lens versus thick lens approximation. A Mathematica program was written to solve for the K value of the QM02, in the thick lens approximation, that provides point to point focusing from YC07 to BPMS1, and to compare this number with the value obtained using the thin lens approximation. The length of QM02 used in the thick lens calculation is the effective length determined by magnetic measurements of 0.108 m. The result of the Mathematica calculation is that the thin lens approximation predicts less magnet strength is required to produce the same focusing by about 1.3%. When both the distance correction and the thick lens approximation are taken into account, the result is: K{sub thick} - K{sub mm}/K{sub thick} = -7.6% where K{sub mm} is the value obtained from magnetic measurements and K{sub thick} is the value of K obtained from the focal length measurement in the thick lens approximation. That is, QM02 acts weaker than it was measured by magnetic measurements by 7.6%. This is remarkably close to Paul's original estimate. The unexpected weakness could in principle be due to several things: shorted turns, a current calibration error, magnetic measurement error; it could even be due to the presence of an gradient from QM01, which is of opposite sign and right next to QM02, despite it being set to zero current (although it was not DAC-zeroed). Plans have been implemented to remove and replace QM02 at the earliest ROD.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Thermal Neutron Scattering Law Evaluation for Zirconium Carbide and Critical Mass Calculations

Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Towards a unified picture of the water self-ions at the air-water interface: a density functional theory perspective

The propensities of the water self ions, H3O+ and OH- , for the air-water interface has implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O+ and/or OH- prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs. interfacial behavior of H3O+ and OH- that employs forces derived from density functional theory with a generalized gradient approximation exchangecorrelation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O+ as a function of the position of the ion in a 215-molecule water slab. The PMF is flat, suggesting that H3O+ has equal propensity for the air-water interface and the bulk. We compare the PMF for H3O+ to our previously computed PMF for OH- adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs. the bulk are connected with interfacial propensity. We find that the solvation shell of H3O+ is only slightly dependent on its position in the water slab, while OH- partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions. DJT was supported by National Science Foundation grant CHE-0909227. CJM was supported by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated for the Department of Energy by Battelle. The potential of mean force required resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC05-00OR22725. The remaining simulations and analysis used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. at at Lawrence Berkeley National Laboratory. MDB is grateful for the support of the Linus Pauling Distinguished Postdoctoral Fellowship Program at PNNL.

Baer, Marcel D.↗
DOE OSTI2015

Charge Recombination, Transport Dynamics, and Interfacial Effects in Organic Solar Cells

The need for renewable sources of energy is well known. Conversion of sunlight to electricity using solar cells is one of the most important opportunities for creating renewable energy sources. The research carried out under DE-FG02-08ER46535 focused on the science and technology of “Plastic” solar cells comprised of organic (i.e. carbon based) semiconductors. The Bulk Heterojunction concept involves a phase separated blend of two organic semiconductors each with dimensions in the nano-meter length scale --- one a material that functions as a donor for electrons and the other a material that functions as an acceptor for electrons. The nano-scale inter-penetrating network concept for “Plastic” solar cells was created at UC Santa Barbara. A simple measure of the impact of this concept can be obtained from a Google search which gives 244,000 “hits” for the Bulk Heterojunction solar cell. Research funded through this program focused on four major areas: 1. Interfacial effects in organic photovoltaics, 2. Charge transfer and photogeneration of mobile charge carriers in organic photovoltaics, 3. Transport and recombination of the photogenerated charge carriers in organic photovoltaics, 4. Synthesis of novel organic semiconducting polymers and semiconducting small molecules, including conjugated polyelectrolytes. Following the discovery of ultrafast charge transfer at UC Santa Barbara in 1992, the nano-organic (Bulk Heterojunction) concept was formulated. The need for a morphology comprising two interpenetrating bicontinuous networks was clear: one network to carry the photogenerated electrons (negative charge) to the cathode and one network to carry the photo-generated holes (positive charge) to the anode. This remarkable self-assembled network morphology has now been established using Transmission electron Microscopy (TEM) either in the Phase Contrast mode or via TEM-Tomography. The steps involved in delivering power from a solar cell to an external circuit are the following: • Photo-excitation of the donor (or the acceptor). • Charge transfer with holes in the donor domain and electrons in the acceptor domain. • Sweep-out to electrodes prior to recombination by the internal electric field. • Energy delivered to the external circuit. Each of these four steps was studied in detail using a wide variety of organic semiconductors with different molecular structures. This UC Santa Barbara group was the first to clarify the origin and the mechanism involved in the ultrafast charge transfer process. The ultrafast charge transfer (time scale approximately 100 times faster than the first step in the photo-synthesis of green plants) is the fundamental reason for the potential for high power conversion efficiency of sunlight to electricity from plastic solar cells. The UCSB group was the first to emphasize, clarify and demonstrate the need for sweep-out to electrodes prior to recombination by the internal electric field. The UCSB group was the first to synthesize small molecule organic semiconductors capable of high power conversion efficiencies. The results of this research were published in high impact peer-reviewed journals. Our published papers (40 in number) provide answers to fundamental questions that have been heavily discussed and debated in the field of Bulk Heterojunction Solar Cells; scientific questions that must be resolved before this technology can be ready for commercialization in large scale for production of renewable energy. Of the forty publications listed, nineteen were co-authored by two or more of the PIs, consistent with the multi-investigator approach described in the original proposal. The specific advantages of this “plastic” solar cell technology are the following: a. Manufacturing by low-cost printing technology using soluble organic semiconductors; this approach can be implemented in large scale by roll-to-roll printing on plastic substrates. b. Low energy cost in manufacturing; all steps carried out at room temperature (approx. a factor of ten less than the use of Silicon which requires high temperature processing). c. Low carbon footprint d. Lightweight, flexible and rugged Because of the resolution of many scientific issues, a significant fraction of which were addressed in the research results of DE-FG02-08ER46535, the power conversion efficiencies are improving at an ever increasing rate. During the funding period of DE-FG02-08ER46535, the power conversion efficiencies of plastic solar cells improved from just a few per cent to values greater than 11% with contributions from our group and from researchers all over the world.

14 SOLAR ENERGY↗
DOE OSTI2016

Moist Process Biases in Simulations of the Madden–Julian Oscillation Episodes Observed during the AMIE/DYNAMO Field Campaign

Two Madden–Julian oscillation (MJO) episodes observed during the 2011 Atmospheric Radiation Measurement Program MJO Investigation Experiment (AMIE)/DYNAMO field campaign are simulated using a regional model with various cumulus parameterizations, a regional cloud-permitting model, and a global variable-resolution model with a high-resolution region centered over the tropical Indian Ocean. Model biases in relationships relevant to existing instability theories of MJO are examined and their relative contributions to the overall model errors are quantified using a linear statistical model. The model simulations capture the observed approximately log-linear relationship between moisture saturation fraction and precipitation, but precipitation associated with the given saturation fraction is overestimated especially at low saturation fraction values. This bias is a major contributor to the excessive precipitation during the suppressed phase of MJO. After accounting for this bias using a linear statistical model, the spatial and temporal structures of the model-simulated MJO episodes are much improved, and what remains of the biases is strongly correlated with biases in saturation fraction. The excess precipitation bias during the suppressed phase of the MJO episodes is accompanied by excessive column-integrated radiative forcing and surface evaporation. A large portion of the bias in evaporation is related to biases in wind speed, which are correlated with those of precipitation. Furthermore, these findings suggest that the precipitation bias sustains itself at least partly by cloud radiative feedbacks and convection–surface wind interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2010

Dipolar Fermi gases in anisotropic traps

The quest for quantum degenerate Fermi gases interacting through the anisotropic and long-range dipole-dipole interaction is an exciting and fast developing branch within the cold-atom research program. Recent experimental progress in trapping, cooling, and controlling polar molecules with large electric dipole moments has, therefore, motivated much theoretical effort. In a recent letter, we have briefly discussed the application of a variational time-dependent Hartree-Fock approach to study theoretically both the static and the dynamic properties of such a system in a cylinder-symmetric harmonic trap. We focused on the hydrodynamic regime, where collisions ensure the equilibrium locally. Here we present a detailed theory extended to encompass the general case of a harmonic trap geometry without any symmetry. After deriving the equations of motion for the gas, we explore their static solutions to investigate key properties like the aspect ratios in both real and momentum space as well as the stability diagram. We find that, despite the lack of symmetry of the trap, the momentum distribution remains cylinder symmetric. The equations of motion are then used to study the low-lying hydrodynamic excitations, where, apart from the quadrupole and monopole modes, the radial quadrupole mode is also investigated. Furthermore, we study the time-of-flight dynamics as it represents an important diagnostic tool for quantum gases. We find that the real-space aspect ratios are inverted during the expansion, while that in momentum space becomes asymptotically unity. In addition, anisotropic features of the dipole-dipole interaction are discussed in detail. These results could be particularly useful for future investigations of strongly dipolar heteronuclear polar molecules deep in the quantum degenerate regime.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2012

The Jefferson Lab High Power Broadband THz Facility

Jefferson Lab operates a broadband terahertz (THz) facility, whose light is generated by 300 fs FWHM 135 pC bunches of electrons at an energy of 125 MeV as they traverse a magnetic compression chicane prior to entering the optical cavity of the Jefferson Lab Free Electron Laser (FEL)[1-3]. The repetition rate of the electron bunches, and hence that of the light pulses, can be as high as 75 MHz continuous. Approximately 1 microJoule of energy is emitted for each pulse, which implies an average power of up to 75W, a peak power of several MW and electric fields of MV/m. Niche applications of this unique source are pump-probe dynamics, both in and out of equilibrium, and large-scale imaging. The beamline itself is designed to deliver an f/8 beam 2mm vertical × 3mm horizontal in size and 60% vertically polarized from a subtended angle of ~200 milliradians horizontal by ~135 milliradians vertical in the accelerator. The first mirror is copper-braid cooled. The philosophy of the optical design is to use a relay optics configuration to transport the beam via reflection off metal mirrors through a series of focal points, with the final one being at a diamond window. The optical scheme was developed using the Synchrotron Radiation Workshop (SRW)[4]. This code performs a full calculation of the electric field from a relativistic electron. It does not handle multiparticle coherent enhancement[5], but it is significantly different from all the other synchrotron radiation calculations because it retains a term called the Coulomb term in which the electric field falls off as 1/R2. In addition to calculating the electric field and the intensity, the code also allows one to propagate the field through the optical focusing elements and apertures. The beam exiting the final diamond window in the THz user lab passes into a series of vacuum chambers with considerable flexibility for attenuation, collimation, and/or beam splitting for imaging, spectroscopy and 2 photon pump-probe measurements. In this talk we will discuss the details of the philosophy, optical design and operating experience of this facility and also present ideas for discussion. We deeply appreciate the excellent support of our colleagues without which these facilities would not have been built and operated. This work was supported by the U.S. Dept. of Energy under contract DE-AC05-06-OR23177, and the Commonwealth of Virginia. The FEL facility is supported by the Office of Naval Research. 1. S. V. Benson et al, ?The VUV/IR/THz Free Electron Laser Program at Jefferson Lab?, Nucl. Instr. and Meth. A649 9 (2011). 2. G.L. Carr, M.C. Martin, W.R. McKinney, K. Jordan, G.R. Neil and G.P. Williams ?High Power Terahertz Radiation from Relativistic Electrons?, Nature 420 153-156 (2002). 3. J.M. Klopf et al. ?The Jefferson Lab high power THz user facility?, Nucl. Instrum. Methods A582, 114-116 (2007). 4. O. Chubar, P. Elleaume, "Accurate And Efficient Computation Of Synchrotron Radiation In The Near Field Region", proc. of the EPAC98 Conference, 22-26 June 1998, p.1177-1179. 5. C. J. Hirschmugl, M. Sagurton and G. P. Williams, ?Multiparticle coherence calculations for synchrotron radiation emission? Physical Review A44, 1316 (1991).

Klopf, John↗
DOE OSTI2023

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Abstract]

Historically, the free gas approximation has been used to treat the thermal scattering of neutrons with energies below a few electron-volts (eV) in unevaluated materials. However, this method inadequately reproduces neutron scattering at these energies. Until recently, only a limited number of materials had available thermal scattering law (TSL) files/libraries in the ENDF nuclear data libraries in this energy range. With advancements in atomistic modeling techniques, such as molecular dynamics, ab-initio molecular dynamics, and density functional theory, TSL libraries have become available for many more materials. This is particularly relevant due to the rising interest in several advanced reactor systems that require novel moderator and reflector materials. While quasi-integral and integral benchmarks have been designed to validate historically important moderator materials (such as light water and polyethylene), there is currently a lack of standard validation methods for TSLs, especially when multiple conflicting TSLs exist. To address this issue, the Oak Ridge National Lab Nuclear Data group has been working on utilizing inelastic neutron scattering (INS) measurements combined with transmission (i.e., total cross section) measurements to evaluate and validate TSLs for different materials. We plan to demonstrate how this method has worked on materials such as polyethylene, lucite, and polystyrene. In addition, we will compare the newly created libraries to ENDF libraries for these materials and explain why integral benchmarks should not be used for validation when multiple TSLs exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

The Potential Role of Concentrating Solar Power within the Context of DOE's 2030 Solar Cost Targets

For solar electricity generating technologies to be cost competitive at a large scale with conventionally generated electricity, cost reductions are needed for both concentrating solar power (CSP) and solar photovoltaic (PV) systems. In 2011, the U.S. Department of Energy (DOE) established solar cost targets that corresponded to reducing CSP and PV prices by approximately 75% in order to achieve a levelized cost of electricity (LCOE) of $0.06 per kilowatt-hour (kWh) for both utility-scale PV and high-capacity factor CSP-TES systems in 2020. Utility-scale PV achieved its 2020 cost target in 2017, and recent estimates for the LCOE of CSP-TES with a molten-salt power tower system are approximately $0.10/kWh for projects that are expected to come online in 2020, which represents a substantial reduction since 2010 - when the LCOE for CSP-TES was around $0.21/kWh. To continue the momentum for cost reductions in solar technologies, DOE recently established cost targets for 2030 that would make solar one of the lowest-cost sources of new electricity in the United States. For CSP-based systems, the new targets correspond to an LCOE in 2030 of $0.05/kWh for a dispatchable, high-capacity factor CSP-TES plant configuration. This aggressive target would have been unimaginable a decade ago. However, building on the previously described reduction in CSP-TES costs over the past decade, recent announcements suggest the next phase of projects will continue this downward trend through lower installation costs, attractive financing, longer-duration PPAs, and the ability to capitalize on the value that the flexibility of storage brings CSP. This presentation summarizes a recent report that evaluates the potential impacts of simultaneously achieving the 2030 cost targets for PV and CSP-TES, and it includes a detailed evaluation of the role that CSP-TES could play in realizing those impacts. The scenarios in this analysis are designed to isolate and assess the potential impacts of achieving DOE's 2030 cost targets for CSP-TES and PV, and they do not reflect the potential benefits or system impacts associated with success in other DOE research programs. The low-cost solar scenarios include DOE's 2030 solar cost targets--which are represented via a roughly 50% reduction in LCOE by 2030 (from current levels) with additional cost reductions thereafter representing technology learning and/or improvements that could result from innovation - coupled with a variety of market, technology, and demand assumptions. With these assumptions, the evolution of the contiguous U.S. electricity system is evaluated with NREL's Renewable Energy Deployment System (ReEDS) model, which was specifically designed to represent the temporal and locational value of renewable generation technologies in the U.S. power system. ReEDS relies on system-wide least-cost optimization to estimate the type and location of future generation and transmission capacity. In addition, it accounts for the locational and temporal variations in variable renewable technologies, including the need for new transmission, curtailment, dynamic capacity value, and the need to hold operating reserves to account for the uncertainty and variability of these technologies. This presentation summarizes the key findings that arise from a detailed evaluation of the impacts of achieving DOE's 2030 cost targets for CSP-TES and PV systems, noting again the inherent challenges associated with modeling future scenarios of the large, complex system electricity system in the contiguous United States within the context of various market features.

14 SOLAR ENERGY↗
DOE OSTI2021

New Capability at ORNL: High Precision Uranium Titration

Destructive analytical measurements establish a nuclear facility’s nuclear material inventory and inventory differences for Nuclear Material Accountancy and Control. A nuclear laboratory’s ability to perform accurate high-precision analytical measurements is key for tracking large inventories within a facility’s material balance areas and during production to track material movement through dynamic processes. For uranium, these measurements are made by using several established high-precision measurement protocols. These include isotope dilution mass spectrometry, gravimetry, and potentiometric titrations. Whatever measurement technique a nuclear lab chooses to use, reference materials (RMs) with certified attributes and accompanying uncertainties are used to calibrate measurement systems, and are the cornerstone for accurate results. In addition to calibration, RMs provide for metrological traceability, are used for method development and validation, and thus provide critical evaluations of the appropriateness and performance of analytical processes used. Evaluations may include validation that a method is fit-for-purpose, quantification of systematic and random biases, and the evaluation of long-term and short-term performance metrics. High-precision measurement techniques require that measured attributes be certified to a high degree of precision in the RMs used—ultimately, to a higher degree than that of the measurement technique itself. The US authority on the production of special nuclear material Certified Reference Materials is the NBL Program Office (NBLPO), formally known as New Brunswick Laboratory (NBL). The NBLPO is responsible for the sales and distribution of existing NBL certified reference materials (CRMs) and for the production of the next generation of nuclear RMs. To accomplish its mission, NBLPO is establishing key base capabilities within the DOE laboratory complex that formerly existed at the NBL laboratory. The Nuclear Analytical Chemistry (NAC) section within the Chemical Sciences Division (CSD) at Oak Ridge National Laboratory (ORNL) is currently working with NBLPO to set up laboratory and measurement capabilities to provide measurements and capabilities for production and/or recertification of existing and future CRMs for uranium assay. The NBL-developed high-precision titration (HPT) method is a critically-evaluated, extremely precise and accurate primary method utilized for the determination of uranium content in a variety of uranium materials. The HPT method, combined with detailed balance weighing protocols, provides for an analytical methodology that is unsurpassed in precision and one in which all sources of error have been evaluated, a requirement of CRM certification. HPT capability within the United States was lost with the closure of the labs at NBL. The NAC has been collaborating with NBLPO to stand up and demonstrate the capability to perform uranium assay via HPT at ORNL. HPT can produce results with an expanded uncertainty of approximately 0.01% for pure uranium compounds, with typical precisions of <0.006%. The major tasks required to stand up the method at ORNL were the refurbishment of a dedicated lab and equipment setup, procedure development, analyst training, establishing method-specific quality assurance, and qualification of an analyst to perform the method. This report summarizes these tasks, outlines the documents drafted, and gives the outcome of the qualification titrations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2016

WIPP Certification of the Retrieval Box Assay System at the Advanced Mixed Waste Treatment Plant, Idaho - 16241

The Retrieval Box Assay System (RBAS) at the Advanced Mixed Waste Treatment Plant (AMWTP) is an Imaging Passive Active Neutron system supplemented with a High Purity Germanium (HPGe) gamma system for the quantification of Pu-239 Fissile Gram Equivalent (FGE) and transuranic (TRU) alpha activity of waste containers of various sizes up to 249 cm x 149 cm x 198 cm (Length x Width x Height). The system was originally calibrated in 2002 (with a range of surrogate matrices and radioactive sources in a 213 cm x 122 cm x 122 cm FRP box) and has been operational for over a decade. As of August 2015, the RBAS has measured approximately 18,000 boxes (50,000 m3 of waste) which has subsequently been fed into the facility's treatment process. Its success, and the wide envelope of waste that it can characterize, has led AMWTP to extend the RBAS mission to include WIPP certification. As part of this new project, the instrument must be capable of meeting WIPP's Waste Acceptance Criteria (WAC) which include reporting activities and masses of the 10 WIPP-tracked radionuclides. The system will be required to assay sludge and debris waste streams, shredded waste and large metal items. The container size includes Standard Waste Box (SWBs), Fiberglass Reinforced Plywood (FRPs) and large Shredder Boxes. Sludge matrices will be loaded into 210-liter (55-gallon) drums and placed inside an SWB. After assay, boxes will be over-packed into SLB-II's for shipment to WIPP. The target waste streams will contain mixtures of Weapons Grade plutonium, Heat Source plutonium (high in Pu-238 content) and Mixed Oxide (MOX). The activity of the waste is expected to encompass a wide dynamic range of activity from low-level waste (LLW) to TRU with up to several hundred grams of plutonium. Some packages will have high gamma activity associated with fission and activation products such as Cs-137 and Co-60. The existing RBAS hardware (detectors, electronics, assay chamber etc.) was retained for this project. Two lead filters were provided and installed on the existing filter wheels to attenuate the gamma ray flux from the highest activity boxes to avoid overload of the HPGe detectors. The new software suite includes PC-FRAM (Version 5), which is used for the isotopic analysis of the gamma spectra. Integration of FRAM and neutron results is performed using a software package known as RBAS WIPP Analysis Reporting Software (RWARS). This package also performs matrix dependent background adjustments. Other upgrades provided include automated tests that identify when assay parameters fall outside a predetermined range and a new Quality Checking (QC) application. Minimum Detectable Activity (MDA) values were determined for the active and passive neutron modes as a function of measured matrix correction indices. The resulting Minimum Detectable Concentration (MDC) has been demonstrated to be less than 3700 Bq/g alpha activity (100 nCi/g). Modifications were performed to update the system TMU parameters to accommodate assay of non-standard size boxes and a wider range of matrix materials. A series of calibration confirmation, measurements were performed on site in 2014 using a representative range of surrogate boxes and radioactive sources including Pu. The WIPP WAC performance criteria were successfully demonstrated up to 157 g Pu and the system passed its first Performance Demonstration Program (PDP) cycle in the Spring of 2015. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗