Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model data volumes”

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 1,027 records · Page 57

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

A Predictive Model for Offgas Composition in Waste Glass Melters - Spring 2021 SULI General Audience Abstract

Approximately 56 million gallons of nuclear waste are currently stored in 177 underground storage tanks at the Hanford site in Washington state. A portion of the tank waste is targeted for immobilization by converting the nuclear waste to borosilicate glass through a process called vitrification. During vitrification, emissions are produced in addition to the intended glass product. These emissions, contained in a volume above the molten glass called the plenum, are composed of water vapor, air components, aerosols, particulates, and other gaseous species that result from evaporation, reactions, and air influx. The formation of these gasses causes the surface of the glass to foam. To reduce foaming, sucrose is added to the waste stream before the conversion to glass. Sucrose also reduces the formation of hazardous gasses in the plenum including nitrous oxides, NOx, and carboxides, COx. However, if excess sucrose is present, the reactions that produce the gasses in the plenum will be incomplete, and increased levels of NOx and COx can be observed. The Waste Treatment and Immobilization Plant at the Hanford site is equipped with the technology to treat the plenum emissions. Though many tests have been run to determine the presence of particulates in the plenum emissions, there is does not exist an established method to confidently predict the composition of the gaseous species in the plenum prior to treatment without physical testing. A model was created in the process modeling software Aspen Plus to predict the composition of NOx and COx in the plenum. The model employs a constant stirred tank reactor and user-defined reaction chemistry to determine the plenum composition. To verify the model, the offgas predicted by the model were compared to available offgas data for ten different feed compositions. It is demonstrated that the developed Aspen Plus model is capable of predicting the composition of NOx and COx in the plenum.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The shard test and nanoporomechanics reverse classical paradigm of cement hydration being contractive

Le Chatelier in 1887 and Powers in 1947 demonstrated that the volume of nanoscale C-S-H (calcium silicate hydrate) particles formed during hydration is smaller than the combined volume of the reactants—the anhydrous Portland cement and water. Hydration has thus been considered as contractive. An experiment shows that the opposite is true above the nanoscale. The porous skeleton of cement paste expands as the growing C-S-H particles push each other apart, similar to crystal growth pressure. This is significant for high-performance concretes (HPC) with low water-cement ratios (w/c ≤ 0.4), where chemical self-desiccation lowers pore relative humidity by 40%, compared to just 1% in traditional concretes (w/c ≈ 0.5). Standard American Society for Testing and Materials (ASTM) C1608 tests, using 10 mm thick water-immersed specimens, show large shrinkage because the half-time of water ingress is many decades, unable to offset shrinkage-causing self-desiccation. The present experiment, using a laser microscopy-topography technique, proves the opposite—expansion, evidenced by measuring the length changes of water-immersed HPC shards 0.5 mm thick in which the diffusion halftime, only about one hour, allows continuous resaturation of pores, canceling self-desiccation. The faster diffusion (halftime of one hour) enables continuous pore resaturation, preventing shrinkage. When sealed with paraffin oil, the shards self-desiccate and shrink. These findings align with studies since 2015, showing that models excluding hydration expansion cannot fit test data across various specimen sizes and sealing conditions. The results suggest that standardized ASTM tests for the so-called chemical shrinkage in modern concretes with very low water-cement ratios are misleading and need revision.

Science & Technology - Other Topics↗

Tropospheric and stratospheric wildfire smoke profiling with lidar: mass, surface area, CCN, and INP retrieval

We present retrievals of tropospheric and stratospheric height profiles of particle mass, volume, surface area, and number concentrations in the case of wildfire smoke layers as well as estimates of smoke-related cloud condensation nuclei (CCN) and ice-nucleating particle (INP) concentrations from backscatter lidar measurements on the ground and in space. Conversion factors used to convert the optical measurements into microphysical properties play a central role in the data analysis, in addition to estimates of the smoke extinction-to-backscatter ratios required to obtain smoke extinction coefficients. The set of needed conversion parameters for wildfire smoke is derived from AERONET observations of major smoke events, e.g., in western Canada in August 2017, California in September 2020, and southeastern Australia in January–February 2020 as well as from AERONET long-term observations of smoke in the Amazon region, southern Africa, and Southeast Asia. The new smoke analysis scheme is applied to CALIPSO observations of tropospheric smoke plumes over the United States in September 2020 and to ground-based lidar observation in Punta Arenas, in southern Chile, in aged Australian smoke layers in the stratosphere in January 2020. These case studies show the potential of spaceborne and ground-based lidars to document large-scale and long-lasting wildfire smoke events in detail and thus to provide valuable information for climate, cloud, and air chemistry modeling efforts performed to investigate the role of wildfire smoke in the atmospheric system.

54 ENVIRONMENTAL SCIENCES↗

Using Robotic Manipulators for Radioactive Waste Inspection - 20090

The global nuclear industry has a growing volume of nuclear waste which needs to be scanned, sorted according to its activity and material type, then processed into the correct waste packages for long term storage and disposal. It is vital that there is a detailed understanding of the waste inventory stored in long term waste containers, as knowledge of their contents could predict or prevent any adverse effects in storage. The numerous 'scan and sort' tables which are currently used at many different facilities around the world to sort waste into their correct containers are human operated and require very slow gamma scanning procedures combined with educated guesswork to manually sort the waste. This often leads to excessive conservatisms, with placement of lower activity wastes in higher activity containers, which in turn costs significantly more to store. In the United Kingdom it costs UK Pounds 46 k per cubic meter to store intermediate level waste compared to just UK Pound 2.9 k per cubic meter to store low level waste according to a 2008 Department of Energy and Climate Change report in the UK. A proposed solution to this problem, is the use of a robotic manipulator to automatically inspect the 'scan and sort table' in order to produce an accurate 3D model of the table's waste contents and attach an overlaid radiation map. The radiation map contains spectrometry data and can in consequence be used to distinguish and locate specific radioisotopes. The 3D model should be as accurate as possible in order to allow for a second robot arm with an attached gripper to grasp the objects and place them into their designated long-term storage container. Various scanning procedures are explored in this study including basic raster scanning, adaptive raster scanning and point sampling. The optimal solution will in practice be defined by the required application and activity level of the wastes being inspected. The results presented in this study indicate that it is possible to produce a centimeter accurate 3D model of a mixed assortment of components on a nuclear waste 'scan and sort' table. In addition, it was shown that the waste objects emitting radiation could be accurately identified and located, with an overlaid radiation map. This study is applicable across the nuclear waste management sector. Many of the ideas and concepts developed in this study are applicable in other decommissioning settings for example, dismantling of legacy gloveboxes or routine inspection of nuclear waste packages in storage. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017) (and subsequently titled) Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020) (Final Technical Report)

This is the Final Technical Report for the DoE Measurement Innovation (previously Diagnostic Development) grant DE-SC0006077 titled HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017), and subsequently titled Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 May 2011–17 July 2020. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that will improve understanding of transport, magnetic equilibrium, and electric fields in magnetically confined plasmas. We have developed novel hardware and advanced techniques to expand beam-based diagnostic capabilities, and extend measurements to new operating regimes. Some have been realized using a Heavy Ion Beam Probe (HIBP) having traditional features, while other applications have joined our innovations with subsets of advantageous HIBP features. This work has established smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key Measurement Innovations include: Computer Models and Simulations of HIBP Applications-Computer modeling and simulating of HIBP applications is central to predicting interaction of beam particles with plasmas, designing diagnostic systems, and anticipating measurement characteristics. We have extended the capability and improved measurement fidelity of the HIBP diagnostic on the Madison Symmetric Torus (MST) reversed field pinch (RFP); simulated feasibility of HIBP operation in the Helically Symmetric eXperiment (HSX) optimized stellarator; and investigated the plausibility of HIBP measurements in the ASDEX Upgrade tokamak. We have also developed a new technique that uses ion optics to model a finite phase-space beam. It yields more realistic (than traditional method) estimates of sample volume characteristics, which influence measurement resolution and sensitivity. Extension of Diagnostic Capabilities and Measurements on the MST RFP-Capabilities have been extended using the first and only HIBP installation on an RFP. Challenges associated with operation of the diagnostic on MST include the three-dimensional nature of particle trajectories, temporal topology and amplitude changes in the plasma equilibrium, and strong particle and radiative emission from the plasma. We have addressed these challenges through development of detailed and higher precision diagnostic simulations, calibrated subtraction of noise to resolve secondary ion signals, and implementation of hardware that enables higher fidelity measurement of plasma fluctuations. Development of New Hardware and Measurement Techniques - We have advanced beam-based measurement capabilities via development of new detectors and techniques to determine the poloidal magnetic flux ψ in the plasma, and enable placement of detectors close to the plasma and thus allow a smaller and less expensive diagnostic. We have demonstrated measurement of ψ, an ability that was made possible by development of highly effective noise reduction techniques. We also modeled the effect of non-ideal and instrumentation effects by developing a virtual (simulated) diagnostic. Application of beam-based diagnostics on various devices and multiple magnetic configurations (e.g. stellarator, tokamak, and reversed field pinch) enables investigation of critical physics issues and measurements resulting in the broad parameter space data needed to test and validate theory & modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anomaly Detection and Approximate Similarity Searches of Transients in Real-time Data Streams

Abstract We present Lightcurve Anomaly Identification and Similarity Search ( LAISS ), an automated pipeline to detect anomalous astrophysical transients in real-time data streams. We deploy our anomaly detection model on the nightly Zwicky Transient Facility (ZTF) Alert Stream via the ANTARES broker, identifying a manageable ∼1–5 candidates per night for expert vetting and coordinating follow-up observations. Our method leverages statistical light-curve and contextual host galaxy features within a random forest classifier, tagging transients of rare classes ( spectroscopic anomalies), of uncommon host galaxy environments ( contextual anomalies), and of peculiar or interaction-powered phenomena ( behavioral anomalies). Moreover, we demonstrate the power of a low-latency (∼ms) approximate similarity search method to find transient analogs with similar light-curve evolution and host galaxy environments. We use analogs for data-driven discovery, characterization, (re)classification, and imputation in retrospective and real-time searches. To date, we have identified ∼50 previously known and previously missed rare transients from real-time and retrospective searches, including but not limited to superluminous supernovae (SLSNe), tidal disruption events, SNe IIn, SNe IIb, SNe I-CSM, SNe Ia-91bg-like, SNe Ib, SNe Ic, SNe Ic-BL, and M31 novae. Lastly, we report the discovery of 325 total transients, all observed between 2018 and 2021 and absent from public catalogs (∼1% of all ZTF Astronomical Transient reports to the Transient Name Server through 2021). These methods enable a systematic approach to finding the “needle in the haystack” in large-volume data streams. Because of its integration with the ANTARES broker, LAISS is built to detect exciting transients in Rubin data.

79 ASTRONOMY AND ASTROPHYSICS↗

Hybrid Large Eddy Simulation/Reynolds-Averaged Navier–Stokes Analysis of a Premixed Ethylene-Fueled Dual-Mode Scramjet Combustor

Hydrocarbon fuels offer optimal high energy per volume for scramjet applications for sustained hypersonic flight but require additional residence time due to slower ignition delays (compared with hydrogen fuel). The injection of ethylene at the start of the isolator of a dual-mode scramjet combustor, operating in ramjet mode, allows sufficient mixing to achieve efficient premixed turbulent combustion. A cavity flameholder anchors the flame, supplying sufficient radicals to sustain a stable flame in the high-speed environment. This work investigates flame structure and stabilization limits of a new configuration with a scaled-down cavity embedded in the flow path of the combustor with a strut and insert. The cavity is reduced in size by one-third to enable eventual direct numerical simulations of the flame stabilization process. This work, however, focuses on modeling the full isolator/combustor geometry using a hybrid large eddy simulation/Reynolds-averaged Navier–Stokes simulation strategy. Additionally, particle image velocimetry, planar laser-induced fluorescence, coherent anti-Stokes Raman spectroscopy, and pressure measurements are compared with numerical predictions to analyze and characterize the conditions within the combustor, including flame structure, flow velocities, species composition, and wall pressure. An adjustable air throttle in the extender is able to control the placement of the shock train in the isolator and maintain a stable flame at various equivalence ratios. The simulations show reasonably good agreement with the experimental scalar and velocity data and predict a flame angle consistent with premixed turbulent flame-speed correlations.

42 ENGINEERING↗

Report of METL Mutual Inductance Level Sensor Development for Use in Liquid Metals – FY2023

A robust mutual inductance level sensor (MILS) for use in high temperature liquid metals has been developed at Argonne National Laboratory (ANL). The current design utilizes mineral insulated cables with a 300-series stainless steel sheath, magnesium oxide insulation, and a single copper conductor. Two coils are wrapped on a common core made of a 300-series stainless steel tube. One coil is energized with an alternating current (AC) source, and this electromagnetically couples with the second coil to generate an induced voltage. The voltage is measured to determine the mutual induction between the coils, and the mutual inductance measurement can be used to determine the level of a nearby liquid metal volume. The MILS can be located in an isolating thimble that is fully sealed to a vessel containing liquid metal for ease of maintenance and replacement. When the primary coil is energized with the AC source, the coil not only electromagnetic (EM) couples with the second coil, but also the surrounding liquid metal. The EM coupling with the surrounding liquid metal reduces the EM coupling with the second coil, producing an inverse-linear relationship between liquid metal level and secondary coil voltage. This has all been demonstrated at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at ANL. The most current iteration of the MILS system is the MILS-MK-II. This sensor system has been commissioned in a non-sodium test stand where calibrations were performed using a sodium analog. The calibrations proved to be highly linear and repeatable. The MILS-MK-II has been installed in the METL expansion tank where high temperature sodium tests have been performed. The MILS-MK-II has been calibrated against a known standard at temperature of 300°C, and the calibrations have proved to be highly linear and repeatable. The calibrated MILS-MK-II has been in operation in the METL facility for several thousands of hours at temperatures around 300°C. The experimental data has been used to develop and validate electromagnetic finite element models in COMSOL Multiphysics, and now these models can be used to advance the development of the sensor system. Next steps will include temperature compensation to allow for operation at various temperatures. Additionally, efforts to incorporate a self-calibration methodology are underway.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗

Preliminary IHLW Formulation Algorithm Description

This report documents the initial algorithm that could be used by the Waste Treatment and Immobilization Plant (WTP) in batching high-level waste (HLW) and glass-forming chemicals (GFCs) in the HLW melter feed preparation vessel (MFPV) (HFP-VSL-00001 and -00005). Not all Hanford tank waste can be accommodated by the models developed for this report and significant expansion of the model boundaries could be achievable to reduce the WTP mission life and total canister production count. The immobilized HLW (IHLW) must meet a series of constraints to be acceptable for disposal in the Monitored Geologic Repository, which are contained in the Specification 1 of the Contract (DOE 2000), the Waste Acceptance Product Specifications (WAPS, DOE 1996), and the Waste Acceptance System Requirements Document (WASRD, DOE 2007). The IHLW Waste Form Compliance Plan (WCP, 24590-HLW-PL-RT-07-0001, Rev 3) specifies that the formulation algorithm will be developed and used to comply with the constraints associated with glass composition and properties. This report is not an engineering calculation, does not provide design input, and is not an engineering study. Algorithm inputs include the chemical analyses of the blended HLW in the HLW blend vessel (HBV) (HLP-VSL-00028, the volume and composition of the MFPV heel, the volume and composition of the MFPV after waste addition, the volume and composition of MFPV batch after GFC addition, the compositions of individual GFCs, and the mass of glass in each canister. In addition to these inputs, uncertainties in the HLW composition and processing parameters are included in the algorithm. Using the above inputs, the algorithm calculates the following outputs: 1) the volume of HLW to be transferred from the HBV to the MFPV, 2) the mass of each GFC for addition to the MFPV, 3) the composition of the glass that will be produced along with uncertainties, and 4) the predicted properties, with associated uncertainties, of the resulting IHLW. The algorithm uses the property-composition models to calculate properties with associated uncertainties and compares them with various constraints to ensure that a processable feed is formulated and a compliant IHLW is produced. The GFC additions are determined using an optimization approach to provide high confidence that the HLW glass will meet all product quality requirements and key processing constraints. For most HLW batches there are many possible glass compositions that meet all constraints. In these cases, the glass composition is optimized for a series of target component concentrations and target property values. The algorithm also incorporates process measurement and product quality uncertainties, based on the work of Piepel et al. (2005). Estimates of the various process and measurement uncertainties that affect glass compositions and predicted glass properties have been previously reported (Piepel et al. 2005, 2006) and the impacts of these estimated uncertainties on the IHLW composition envelope that meets product quality and processing-related properties with sufficient confidence were evaluated. The details of work performed to date to develop this initial GFC addition and batching algorithm are summarized in Sections 4 and 5. An example data set is used to illustrate the calculations of the algorithm summarized in Section 6. Finally, in Section 7, there is a statement of the required work to achieve a final operational IHLW formulation control algorithm. This report is not an engineering calculation, does not provide design input, and is not an engineering study.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Isotopic and Fuel Lattice Parameter Trends in Extended Enrichment and Higher Burnup LWR Fuel, Volume II: BWR

Commercial light water reactor (LWR) operators and fuel vendors in the United States are pursuing changes to nuclear fuel that include extended enrichment (EE) and accident-tolerant fuel (ATF) designs to further improve reactor safety and plant economics. Extended fuel enrichments above 5% 235 U pin enrichment and up to 10% 235 U are a subset of high assay low-enriched uranium (HALEU) that may be deployed in commercial US LWRs in the near term. ATF features, such as cladding coatings or alternative cladding materials, are designed to improve fuel system performance under accident conditions. One goal of EE is to improve fuel cycle economy by enabling fuel to be depleted to higher burnup than the typical current limits (62 gigawatt-days per metric ton of uranium [GWd/MTU]). Adoption of EE, ATF, and high burnup (HBU) fuels in the US commercial fleet requires a clear understanding of the effects on core physics parameters and used fuel isotopic content, as well as confidence in the accuracy of computer code predictions over an expanded range of materials, enrichment, and burnup. A thorough understanding of the applicability and adequacy of benchmark data (e.g., criticality, decay heat, isotopic content) for computer code validation is necessary to ensure that appropriate safety margins are maintained. As part of the US Nuclear Regulatory Commission (NRC) agreement number 31310019N0008, “SCALE Code Development, Assessment and Maintenance,” the effects of EE, ATF, and HBU are being assessed for selected representative LWR fuel designs. The project is divided into phases, and this report summarizes the findings of the Phase 1 work, which focuses on lattice physics parameter and used fuel isotopic changes for a conventional GE14 10 x 10 boiling water reactor (BWR) design with GNF-2 part length rod patterns to model a modern BWR assembly design.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nickel-rich, volatile depleted iron meteorites: Relationships and formation processes

Ungrouped iron meteorites Tishomingo, Willow Grove, and Chinga, and group IVB iron meteorites, are Ni-rich. Similarities include enrichments of 10–100 × CI for some refractory siderophile elements, and equivalent depletions in more volatile siderophile elements. Superimposed on the overall enrichment/depletion trend, certain siderophile elements (P, W, Fe, Mo) are depleted relative to elements of similar volatility. All three ungrouped irons derive from parent bodies formed in the early Solar System. Willow Grove and Chinga are characterized by cosmic ray exposure corrected 182 W/ 184 W consistent with metal-silicate segregation on their parent bodies within 1–3 Myr of Solar System formation, within the age range determined for segregation of magmatic iron meteorite parent bodies, including group IVB irons. Tishomingo is characterized by a younger model age 4–5 Myr subsequent to Solar System formation, reflecting either late stage melting resulting from 26 Al decay, or an impact resetting. The discovery of stishovite in Tishomingo, indicating exposure to a minimum shock pressure of 8–9 GPa, is consistent with the latter. Stishovite in Tishomingo and chromite included in troilite-daubréelite in IVB irons allows oxygen isotopic composition comparison between these meteorites. Different mass independent oxygen isotopic compositions of IVB irons and Tishomingo indicate genetically distinct parent bodies. By contrast, mass independent Mo isotopic compositions overlap within analytical uncertainties, indicating a similar, carbonaceous chondrite (CC) type genetic heritage. Molybdenum and 183 W isotopic data for Chinga and Willow Grove indicate derivation from CC type parent bodies. Willow Grove shares Mo and 183 W isotopic compositions with the Ni-rich South Byron Trio (SBT) grouplet and the Milton pallasite. These Ni-rich meteorites likely formed in the same general nebular environment as other CC planetesimals, likely the outer Solar System. Highly siderophile element (HSE) abundances of Willow Grove and Tishomingo are similar to some IVB meteorites, consistent with formation by moderate degrees of fractional crystallization from initial metallic melts with low S and P, and modestly fractionated HSE. The comparable HSE abundances of Tishomingo and Willow Grove to some IVB irons, yet substantially higher Ni concentrations, indicate formation on parent bodies with lower bulk HSE abundances or HSE concentration in proportionally smaller volumes of metal. HSE abundances in Chinga are considerably lower than in IVB irons, highly fractionated, and processes responsible for these remain elusive. For IVB irons and these ungrouped irons, high temperature condensation likely dominated the enrichment and depletion of the refractory and volatile siderophile elements, respectively. Parent body degassing may have also played a role. Relative depletion of volatile siderophile elements is not, however, a universal feature of high-Ni meteorites. The SBT and Milton pallasite are Ni-rich, but less depleted in the more volatile siderophile elements. Nickel enrichment was likely driven by oxidation of Fe metal during parent body accretion or core segregation. Oxidation of the Tishomingo and Willow Grove parent bodies may have occurred at ~ IW + 1, indicated by relative Mo and W depletions due to metal/water reaction during differentiation. Late-stage reduction, indicated by the presence of Cr-bearing sulfides in Tishomingo and IVB irons, may have resulted from exhaustion of the oxidant.

36 MATERIALS SCIENCE↗

Chemometric modeling of power conversion efficiency of organic dyes in dye sensitized solar cells for the future renewable energy

The present study reports chemometric modeling of power conversion efficiency (PCE) of dye sensitized solar cells (DSSCs) using the biggest available data set till date which comprises around 1200 dyes covering 7 chemical classes. To extract the best structural features required for higher PCE, we have developed multiple partial least squares (PLS) quantitative structure-property relationship (QSPR) models for the Triphenylamine, Phenothiazine, Indoline, Porphyrin, Coumarin, Carbazole and Diphenylamine chemical classes using descriptors derived from the best subset selection method followed by selection of best five models in each dataset based on the Mean Absolute Error (MAE) values. The models were validated both internally and externally followed by the consensus predictions employing “Intelligent Consensus Predictor” tool to examine whether the quality of predictions can be improved with the “intelligent” selection of multiple PLS models. The quality of predictions for the respective external sets showed that the consensus models (CM) are better than the individual models (IM) in most of the cases. From the insights of the developed models, we concluded that attributes like a packed structure toward higher conductivity of electrons, auxiliary donor fragment of aromatic tertiary amines, number of thiophenes inducing the bathochromic shift and augmenting the absorption, presence of additional electron donors, enhancement of electron-donating abilities, number of non-aromatic conjugated C(sp 2 ) which helps as conjugation extension units to broaden the absorption and highly conjugated $π$-systems exert positive contributions to the PCE. On the contrary, features negatively contributing to PCE are the followings: fragments which lower the tendency of localized $π–π*$ transition, fragments related to larger volume and surface area of dyes along with hydrophobicity resulting in poor adhesion, fragment RC = N causing dye hydrolysis, steric hindrance for π electronic mobility, fragments enhancing polarity, etc. Finally, the identified features from the best QSPR model of the coumarin dataset was employed in designing of ten more efficient coumarin dyes (predicted %PCE ranging from 8.93 to 10.62) than the existing ones.

14 SOLAR ENERGY↗

Large Area Airborne Contamination Monitoring

Current D&D operations at Hanford have demonstrated a flaw in the current state of the art capability of defining airborne contamination boundaries. Airborne particulate emissions of Pu-239 from CM2H operations on the Hanford Pu Finishing Plant were detected well beyond areas controlled for airborne plutonium, putting numerous workers at risk for radiological assimilations. Forty-two Hanford workers had positive bioassay results for plutonium uptake. Work areas determined to be free of airborne contamination as measured by state-of-the-art continuous alpha monitoring (CAM) systems were later discovered not to be. Airborne contamination was scattered to a such degree that CAM systems were surveying insufficient volumes of air to accurately predict where respiratory protection was required. A $\$ 19,000$ state of the art portable CAM system used at both Hanford and at SRS screens air for alpha activity at the rate of 1 CFM. This research project will explore the feasibility of using commercially available air purifying systems to screen significantly larger volumes of air for airborne actinide contamination. Commercial household HEPA based air purifying units screen large volumes of air, a simple $\$ 300$ unit being looked at for this study could screen a volume of air that would require over 300 of the $\$ 19,000$ CAMs. The goal of the project was to determine the feasibility of building a detection system capable of screening the air purifiers pre-filters and HEPA filters for any actinide bearing hot particles that have been captured. The simplest technologies available to Field-Radiological-Protection personnel are hand-held alpha field survey meters. With the size and the corrugated structure of these HEPA filters, hand screening is not practical. This project researched alternatives capable of automating this function. Two potential technologies were examined. One technology evaluated was the measurement of xray emissions from the actinides, taking advantage of the higher emission rates of the L x-rays compared to the gamma emissions of isotopes such as Pu-238 or Pu-239. High resolution to low resolution systems were evaluated. Branching ratio data on the actinides in the literature is incomplete, and a study was initiated to measure branching ratios of the various L emissions of a number of the plutonium isotopes. A Monte Carlo simulation was also run to calculate the effect the attenuating HEPA media would have on x-ray measurements on spectrometers centered in each quadrant above the surface of the HEPA. The second technology evaluated was measuring the light from a scintillating media applied to the surface of the HEPA filter using a ThorLabs Scientific-Grade Digital Camera Model 8051M. A light tight enclosure capable of holding the 2’x2’x1” HEPA filter was fabricated. Numerous scintillation media were evaluated for both the sensitivity and for resolution of measuring alpha activity. Figure 1 shows images of a 4”x4” section of HEPA filter with a 3 x 3 grid of Pu-238 contamination loaded inside the crevices of the HEPA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Wellbore Fracture Imaging Using Inflow Detection Measurements

One of the most striking measurements taken during DOE’s EGS Collab project at the 4850-foot depth location was the so-called ‘sewer cam’, which enabled direct visualization of the flow of water into the production well through fractures during the stimulation. The ability to see directly which fractures were flowing and (roughly) how much was a breakthrough in understanding the topology of the created fracture network. Achieving this kind of fracture flow imaging at FORGE would be more challenging because of the 225°C temperature, but equally or even more valuable if it could be achieved. In 2017, a joint project between Sandia and Stanford developed a downhole tool concept to measure the enthalpy of multiphase fluid entering a geothermal well from individual fractures (Gao et al., 2017). For the FORGE project, measuring enthalpy is of less interest because the fluid is expected to be single-phase liquid water. However, the foundation of the device was the measurement of chloride ion concentration, which could form the basis for a direct measurement of inflow from fractures. During the 2017 project, this novel chloride sensing system was implemented into a laboratory test instrument, and we confirmed the capability of the system to measure the ion concentration of fluid entering a model wellbore through a small entry port. The wellbore was a 6-inch diameter model well, and the port was approximately 0.08 inch (2mm) in diameter. The device could measure the chloride concentration accurately even when the well was flowing in a bubbly flow. Given its accuracy, the tool should be able to identify locations of water entering the wellbore even if the ion concentration differs only slightly from that of the water in the well. It is likely that different fractures may flow slightly different chloride concentrations, which would make it feasible to detect individual fractures as well as to estimate the volume of their flow. Ultimately, we could also recognize different fractures flowing back significantly different ion concentrations after fracturing in the FORGE wells. This could be realized by adding different ions in the fracturing fluids in different fractures created at different stages of stimulation (and modifying the tool to include different ion specificity). Sandia’s tool was shown during the study to have the capability to withstand the 225°C temperature, and the electrochemical sensing elements were tested in the laboratory to 225°C at 1500 psia for 24 hours. An early implementation of the fully integrated downhole electrochemical tool, including high-temperature electronics, robust housing, and wireline truck interface, had previously been constructed and tested successfully at Sandia; thus, hardware development tasks focused on advancing the technology readiness level (TRL) of this promising technology for FORGE deployment, rather than on developing a new scientific basis for its operation. The data collection electronics in this tool allowed for several other sensors (pressure, temperature, flow spinner) to be implemented in parallel as well. The research was a new collaboration between Stanford and Sandia to modify and refine the tool for FORGE deployment, to make the downhole measurements, and to characterize the evolving fractures.

15 GEOTHERMAL ENERGY↗