Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A Parallel Computing Framework for Finding the Supported Solutions to a Biobjective Network Optimization Problem

Abstract Solving multi‐objective combinatorial optimization problems can quickly become computationally challenging when applied to large networks generated from big data. Present‐day first‐rate Mixed Integer Programming (MIP) solvers have parallelism built‐in to take advantage of multicore architectures; but specialized network optimization algorithms that can often solve graph problems more efficiently than a general MIP solver are typically programmed serially. Thus, these network algorithm implementations do not take full advantage of modern multicore computer capabilities. Many parallel computing languages include a fork/join framework into their concurrency packages that allow for simple conversion of recursive serial algorithms to operate in parallel without having to use complicated low‐level message‐passing interfaces. Fork/join is particularly well suited to divide and conquer algorithms such as Non‐Inferior Set Estimation (NISE) method for computing the supported Pareto front of a multi‐objective optimization problem. This work develops a simple general parallel NISE (pNISE) implementation using the Java 7 concurrency tools, and then offers results from a specific implementation of solving a bi‐objective shortest path problem. The results indicate that this method is capable of sizeable speed‐ups on both small‐scale personal computers as well as large‐scale shared memory supercomputers. Finally, other network problems suitable to this method are discussed, including multi‐objective variants of the minimum spanning tree problem, transportation problem, assignment problem, maximum flow problem, and the minimum‐cost flow problem. Copyright © 2015 John Wiley & Sons, Ltd.

Medrano, Fernando Antonio↗
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

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

Flexible nuclear screening approximation to the two-electron spin–orbit coupling based on ab initio parameterization

The derivation, implementation, and validation of a new approximation to the two-electron spin–orbit coupling (SOC) terms is reported. The approximation, referred to as flexible nuclear screening spin–orbit, is based on the effective one-electron spin–orbit operator and accounts for two-electron SOC effects by screening nuclear charges. A highly flexible scheme for the nuclear screening is developed, mainly using parameterization based on ab initio atomic SOC calculations. Tabulated screening parameters are provided for contracted and primitive Gaussian-type basis functions of the ANO-RCC basis set for elements from H to Cm. The strategy for their adaptation to any other Gaussian basis set is presented and validated. A model to correct for the effect of splitting of transition metal d orbitals on their SOC matrix elements is introduced. The method is applied to a representative set of molecules, and compared to exact treatment and other approximative approaches at the same level of relativistic theory. The calculated SOC matrix elements are in very good agreement with their “exact” values; deviation below 1% is observed on average. The presented approximation is considered to be generally applicable, simple to implement, highly efficient, and accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Microscopic study of the coupled-wire construction and plausible realization in spin-dependent optical lattices

Coupled-wire constructions offer particularly simple and powerful models to capture the essence of strongly correlated topological phases of matter. They often rely on effective theories valid in the low-energy and strong-coupling limits, which impose severe constraints on the physical systems where they could be realized. In this work, we investigate the microscopic relevance of a class of coupled-wire models and their possible experimental realization in cold-atom experiments. We connect with earlier results and prove the emergence of fractional quantum Hall states in the limit of strong interwire tunneling. Contrary to previous studies relying on renormalization group arguments, our microscopic approach exposes the connection between coupled-wire constructions and model wave functions in continuum Landau levels. Then, we use exact-diagonalization methods to investigate the appearance of these fractional quantum Hall states in more realistic settings. We examine the parameter regimes where these strongly correlated phases arise, and provide a way to detect their appearance in cold-atom experiments through standard time-of-flight measurements. Motivated by this experimental probe, we finally propose a realization of our model with cold atoms in spin-dependent optical lattices. Our estimates show that the previous fractional quantum Hall phases lie within experimentally accessible parameter regimes, giving a viable route toward their experimental study.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Application of a Triple Bubbler Inside Molten Salt for Pyroprocessing Safeguards Purposes

The application of a triple bubbler system within molten salt environments, specifically a eutectic mixture of LiCl-KCl, is explored for nuclear safeguards purposes, focusing on nuclear material accountancy (NMA) and holdup accountancy in pyroprocessing techniques. Traditional holdup accountancy methods, which involve periodic sampling and lengthy lead times, are being supplemented by the development of actinide sensors using electrochemistry. However, these methods alone are insufficient to meet regulatory requirements set by agencies like NRC and IAEA. The triple bubbler device measures salt level, density, and surface tension simultaneously using three dip tubes of varying diameters and heights. This study investigates the accuracy and robustness of the triple bubbler in molten salt conditions at temperatures ranging from 450-550°C. The updated bubbler design, utilizing tantalum metal and larger tube geometry, demonstrated minimal plugging and maintained measurement accuracy within acceptable uncertainties: 0.611% for depth, 1.34% for density, and approximately 40% for surface tension. Results indicate that the triple bubbler system is effective for real-time, accurate measurements in pyroprocessing environments, making it a promising tool for ensuring compliance with nuclear safeguards and supporting nuclear non-proliferation efforts.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Cone beam x-ray luminescence computed tomography: A feasibility study

Purpose: The appearance of x-ray luminescence computed tomography (XLCT) opens new possibilities to perform molecular imaging by x ray. In the previous XLCT system, the sample was irradiated by a sequence of narrow x-ray beams and the x-ray luminescence was measured by a highly sensitive charge coupled device (CCD) camera. This resulted in a relatively long sampling time and relatively low utilization of the x-ray beam. In this paper, a novel cone beam x-ray luminescence computed tomography strategy is proposed, which can fully utilize the x-ray dose and shorten the scanning time. The imaging model and reconstruction method are described. The validity of the imaging strategy has been studied in this paper. Methods: In the cone beam XLCT system, the cone beam x ray was adopted to illuminate the sample and a highly sensitive CCD camera was utilized to acquire luminescent photons emitted from the sample. Photons scattering in biological tissues makes it an ill-posed problem to reconstruct the 3D distribution of the x-ray luminescent sample in the cone beam XLCT. In order to overcome this issue, the authors used the diffusion approximation model to describe the photon propagation in tissues, and employed the sparse regularization method for reconstruction. An incomplete variables truncated conjugate gradient method and permissible region strategy were used for reconstruction. Meanwhile, traditional x-ray CT imaging could also be performed in this system. The x-ray attenuation effect has been considered in their imaging model, which is helpful in improving the reconstruction accuracy. Results: First, simulation experiments with cylinder phantoms were carried out to illustrate the validity of the proposed compensated method. The experimental results showed that the location error of the compensated algorithm was smaller than that of the uncompensated method. The permissible region strategy was applied and reduced the reconstruction error to less than 2 mm. The robustness and stability were then evaluated from different view numbers, different regularization parameters, different measurement noise levels, and optical parameters mismatch. The reconstruction results showed that the settings had a small effect on the reconstruction. The nonhomogeneous phantom simulation was also carried out to simulate a more complex experimental situation and evaluated their proposed method. Second, the physical cylinder phantom experiments further showed similar results in their prototype XLCT system. With the discussion of the above experiments, it was shown that the proposed method is feasible to the general case and actual experiments. Conclusions: Utilizing numerical simulation and physical experiments, the authors demonstrated the validity of the new cone beam XLCT method. Furthermore, compared with the previous narrow beam XLCT, the cone beam XLCT could more fully utilize the x-ray dose and the scanning time would be shortened greatly. The study of both simulation experiments and physical phantom experiments indicated that the proposed method was feasible to the general case and actual experiments.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Innovative approaches to neutron beam optimization: A case study of FP14 at LANSCE

This paper describes an experimental method to acquire high resolution energy- and spatially- resolved neutron beam spots using the time-gated neutron imaging system with Teledyne Pi-MAX4 camera. These experimental data offer a unique opportunity for benchmarking beam spot simulations. High-quality simulations depend significantly on a high-fidelity geometry model, which can be challenging for legacy facilities. We informed our MCNPX geometry model by latest metrology survey employing Leica laser tracker ATS600. It gave us a high fidelity description of our facility geometry. Such a robust integration of novel tools and methods yields a previously unattainable level of accuracy in both predicting and capturing neutron beam spots, marking a notable advancement over traditional methods reliant on static image plates. Here, to demonstrate the practical application of these tools, we are showing a non-uniform beam spot challenge at our Flight Path 14 (FP14) at the Los Alamos Neutron Science Center (LANSCE). Our precise MCNPX prediction of beam spot shifting as function of neutron energy was confirmed by experimental beam spot measured with extremely high level of detail. Results of this research demonstrate a significant leap in neutron beam optimization at LANSCE and set a new benchmark in beam spot characterization. The advanced methods presented here have potential for adoption at similar research facilities worldwide, aiming at substantial improvement in neutron beam delivery for experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2004

A Web Centric Architecture for Deploying Multi-Disciplinary Engineering Design Processes

There are continuous needs for engineering organizations to improve their design process. Current state of the art techniques use computational simulations to predict design performance, and optimize it through advanced design methods. These tools have been used mostly by individual engineers. This paper presents an architecture for achieving results at an organization level beyond individual level. The next set of gains in process improvement will come from improving the effective use of computers and software within a whole organization, not just for an individual. The architecture takes advantage of state of the art capabilities to produce a Web based system to carry engineering design into the future. To illustrate deployment of the architecture, a case study for implementing advanced multidisciplinary design optimization processes such as Bi-Level Integrated System Synthesis is discussed. Another example for rolling-out a design process for Design for Six Sigma is also described. Each example explains how an organization can effectively infuse engineering practice with new design methods and retain the knowledge over time.

Woyak, Scott↗
DOE OSTI2019

Direct evaluation of dynamical large-deviation rate functions using a variational ansatz

We describe a simple form of importance sampling designed to bound and compute large-deviation rate functions for time-extensive dynamical observables in continuous-time Markov chains. We start with a model, defined by a set of rates, and a time-extensive dynamical observable. We construct a reference model, a variational ansatz for the behavior of the original model conditioned on atypical values of the observable. Direct simulation of the reference model provides an upper bound on the large-deviation rate function associated with the original model, an estimate of the tightness of the bound, and, if the ansatz is chosen well, the exact rate function. The exact rare behavior of the original model does not need to be known in advance. We use this method to calculate rate functions for currents and counting observables in a set of network- and lattice models taken from the literature. Straightforward ansätze yield bounds that are tighter than bounds obtained from Level 2.5 of large deviations via approximations that involve uniform scalings of rates. Here, we show how to correct these bounds in order to recover the rate functions exactly. Our approach is complementary to more specialized methods and offers a physically transparent framework for approximating and calculating the likelihood of dynamical large deviations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Confidence-Based Feature Acquisition

Confidence-based Feature Acquisition (CFA) is a novel, supervised learning method for acquiring missing feature values when there is missing data at both training (learning) and test (deployment) time. To train a machine learning classifier, data is encoded with a series of input features describing each item. In some applications, the training data may have missing values for some of the features, which can be acquired at a given cost. A relevant JPL example is that of the Mars rover exploration in which the features are obtained from a variety of different instruments, with different power consumption and integration time costs. The challenge is to decide which features will lead to increased classification performance and are therefore worth acquiring (paying the cost). To solve this problem, CFA, which is made up of two algorithms (CFA-train and CFA-predict), has been designed to greedily minimize total acquisition cost (during training and testing) while aiming for a specific accuracy level (specified as a confidence threshold). With this method, it is assumed that there is a nonempty subset of features that are free; that is, every instance in the data set includes these features initially for zero cost. It is also assumed that the feature acquisition (FA) cost associated with each feature is known in advance, and that the FA cost for a given feature is the same for all instances. Finally, CFA requires that the base-level classifiers produce not only a classification, but also a confidence (or posterior probability).

Wagstaff, Kiri L.↗
DOE OSTI2020

Risk Model for EM Decision Support Toolsets

The Government Office of Accountability (GAO) has published several reports identifying the need for the Department of Energy (DOE) Office of Environmental Management (EM) to address mounting costs for DoE's cleanup program. DOEEM could greatly benefit from independent decision tool-sets/models that would allow them to evaluate options and inform business decision at the enterprise level considering site-specific life cycle costs and system plans. Program goals: Develop a tool-set that enables EM to independently evaluate alternatives, assess outcomes from different contracting strategies, and inform critical decisions for the enterprise. Project goals: Adaption of a Risk Model for integration with complimentary tool-sets for project level decision making. Operational Events: Discrete event model that evaluates operational variables (e.g. capacity, throughput, maintenance constraints) and identify bottlenecks. Identifying and Bounding Project Risk: Identify risks that impact confidence in meeting goals (e.g. production). Life cycle Cost: Evaluate impacts of staffing levels, inventory, and capital investments on life cycle costs. Methods and approach: Monte Carlo Analysis is being executed to generate results: Input derives from Risk Register Data; Tied to Projected and Target Schedules; Incorporates float duration within the model. Assumes associated risk mitigation actions being completed within a timeline of five fiscal years. Metrics include: Confidence in Meeting Production Goal; Confidence in Safety Standards; Confidence in Continuous Operation; Other Metrics can be added as appropriate regarding specific site needs. The adapted risk model can be used as a stand alone decision tool or can be integrated complementary tool-sets (i.e. process and cost models) for project-specific decisions. These support tool-sets can then be integrated with others for site- and complex- level evaluations. Future work includes designing an adaptable and modular framework that would allow integration of multiple tool-sets for holistic and/or targeted evaluation of alternative strategies for decision making that could lead to risk and cost reduction across the enterprise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Comparison of Shape Optimization Methods for Heat Exchanger Fins Using Computational Fluid Dynamics

Inverse design techniques are one way to leverage advances in 3D printing, artificial intelligence, and computational resources to achieve increased performance of heat exchangers. Two optimization techniques (genetic algorithm and particle swarm) and three geometry representations (binary level set, composite Bézier, and free form deformation) are used to increase the heat transfer and reduce the pressure drop of a heat exchanger fin. After running 210,810 OpenFOAM simulations, results indicate that a significant performance increase of the fin can be realized in less than 48 hrs, allowing for such a process to be integrated in traditional design processes. The best design increased the performance of the objective function, compared to the baseline rectangular geometry by 75%. A custom distributed infrastructure was built, allowing for all methods to reach 95% of the final objective values in a little over 4 hrs, handling 1674 OpenFoam simulations per hour.

42 ENGINEERING↗
DOE OSTI2014

An ab initio investigation of the ground and low-lying singlet and triplet electronic states of XNO{sub 2} and XONO (X = Cl, Br, and I)

A systematic ab initio treatment of the nitryl halides (XNO{sub 2}) and the cis- and trans- conformers of the halide nitrites (XONO), where X = Cl, Br, and I, have been carried out using highly correlated methods with sequences of correlation consistent basis sets. Equilibrium geometries and harmonic frequencies have been accurately calculated in all cases at the explicitly correlated CCSD(T)-F12b level of theory, including the effects of core-valence correlation for the former. Where experimental values are available for the equilibrium structures (ClNO{sub 2} and BrNO{sub 2}), the present calculations are in excellent agreement; however, the X-O distances are slightly too long by about 0.01 Å due to missing multireference effects. Accurate predictions for the iodine species are made for the first time. The vertical electronic excitation spectra have been calculated using equation-of-motion coupled cluster methods for the low-lying singlet states and multireference configuration interaction for both singlet and triplet states. The latter also included the effects of spin-orbit coupling to provide oscillator strengths for the ground state singlet to excited triplet transitions. While for ClNO{sub 2} the transitions to excited singlet states all occur at wavelengths shorter than 310 nm, there is one longer wavelength singlet transition in BrNO{sub 2} and two in the case of INO{sub 2}. The long wavelength tail in the XNO{sub 2} species is predicted to be dominated by transitions to triplet states. In addition to red-shifting from X = Cl to I, the triplet transitions also increase in oscillator strength, becoming comparable to many of the singlet transitions in the case of INO{sub 2}. Hence in particular, the latter species should be very photolabile. Similar trends are observed and reported for the halogen nitrites, many of which for the first time.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRSDate not supplied

Generating Essential Climate Variables from Multiple Satellite Hyperspectral Remote Sensors

Hyperspectral observations from satellite-based sensors provide high information content for the Earth’s atmospheric and surface properties. Traditionally, long-term climate products are derived by performing spatial and temporal averaging of level-2 satellite products. There are two shortcomings of this approach. First, it is a time-consuming process to generate level-2 data products since modern hyperspectral satellite sensors have millions of observations each day with thousands of spectral channels for each observation. Secondly, differences in level-2 retrieval algorithms can lead to errors in the fused multi-satellite data. We have developed a radiometrically consistent spectral fingerprinting method, which overcomes the above-mentioned shortcomings, to derive climate change signals from multiple satellite sensors using spatiotemporally averaged level-1 data. We have applied this method to data collected from Atmospheric Infrared Sounder (AIRS) on Aqua satellite and Cross-track Infrared Sounder (CrIS) on SNPP and NOAA20 and generated decade-long climate data records for atmospheric temperature, water vapor, cloud, trace gases, and surface skin temperature. A key component to this work is a set of observational-based radiative kernels produced from CrIS level-1 data using a single field of view (SFOV) optimal estimation retrieval algorithm. Only limited CrIS level-1 data (e.g., 1-2 years of data) are needed to the derive radiative kernels. Our Principal Component-based Radiative Model (PCRTM) enables us to perform SFOV retrievals under all sky conditions and provides radiative kernels (including those for clouds) needed by the spectral fingerprinting method. In this presentation, we will describe the basic methodology, the details of the algorithm, and results from NASA Aqua AIRS and Suomi-NPP CrIS data. The method can be applied to study future hyperspectral remote sensors such as CLARREO (Climate Absolute Radiance and Refractivity Observatory) Pathfinder (CPF), Tropospheric Emissions: Monitoring of Pollution (TEMPO), Surface Biology and Geology (SBG), Atmosphere Observing System (AOS).

Xu Liu↗
DOE OSTI2018

Advanced Laboratory and Field Arrays: Novel Design and Assessment Methodologies for Wave Energy Converter Design (Wave-SPARC) (Task 8)

This task focused on advancing the mission of the Wave-SPARC project by (1) establishing a user-centered Set-Based Design methodology to enable industrial partners to design higher performance WECs, (2) assessing the means through which WEC concepts and WEC design processes are evaluated by employing and potentially improving Wave-SPARC’s Technological Performance Level (TPL) metric, and (3) clarifying how the Set-Based Design approach to WEC design would be applicable for devices of different scales and/or maritime markets.

13 HYDRO ENERGY↗
DOE OSTI2024

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗
DOE OSTI2024

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗
DOE OSTI2025

Corrosion Testing of Refractory inContact with Molten Glasses Designed for Waste Vitrification - VSL Touchpoint Matrix Glasses

It is known that the predictive life of the refractory ceramic liner of nuclear waste glass melters is conservative, as demonstrated by performance of these materials such as in the Defense Waste Processing Facility (DWPF). The motivation for this task is to maximize the useful life of the melters that will be operated at the Waste Treatment and Immobilization Plant (WTP), which will in turn minimize procurement and disposal costs and melter outage times, as well as to identify maximum loadings in the waste glass of those species that corrode melter components. This task was initiated jointly with Pacific Northwest National Laboratory (PNNL) with the objective to develop a methodology and model to enable more accurate prediction of refractory service life under prototypic conditions from laboratory-scale material corrosion tests. Refractory corrosion is generally reported as physical material loss, measured in units of distance (e.g., inch) or as physical material loss rate, measured in units of distance per time (e.g., inch/day). In post-operational melters, the refractory corrosion is measured directly, sometimes reported as corrosion depth. Crucible tests are used in the laboratory to accelerate the refractory corrosion to facilitate a meaningful measurement in a commensurate amount of time. Crucible tests are particularly useful in understanding refractory corrosion across a large glass composition space, where operational testing would be prohibitive. Some of the critical parameters that are known to influence refractory corrosion by molten glass in a crucible test are temperature, system redox, molten salt phases, glass chemistry, and test duration. The majority of data collected for Monofrax® K-3 (hereafter referred to as K-3) corrosion is from crucible tests, but a small amount comes directly from scaled and production melters. Crucible test data has been collected under varying conditions, whereas data collected from operational melters is relatively fewer and represents conditions specific to the melter campaign. The result is that the published data can be grouped and analyzed in multiple ways, not all of which are readily comparable. The Standard Test Method for Isothermal Corrosion Resistance of Refractories to Molten Glass (ASTM C621) outlines the general guidelines used across industry. That method describes a sealed, static test in which the surface area of the refractory coupon and the volume of glass are fixed. A significant portion of the crucible data pertaining to nuclear waste glasses has been collected in a modified configuration; the most notable differences being the surface area of the refractory coupon to volume of the glass and use of a method for bubbling the melt. To our knowledge, the influence of those parameters on the refractory corrosion has not been quantified. In this work, it was determined that static tests and bubbled tests would be performed. Savannah River National Laboratory (SRNL) was tasked with setting up and performing static testing while PNNL was tasked with setting up and performing bubbled testing. Initial activities were performed to establish laboratory methods that reproduce data comparable to existing data sets of K-3 refractory corrosion by low activity waste (LAW) and high-level waste (HLW) glass compositions. Later activities were focused on refining the test parameters to establish a standard test practice to be used between Laboratories and collecting additional data to be used in the enhanced waste glass model development. This document serves primarily to convey the refractory loss measurement results from corrosion testing of K-3 refractory with waste glass compositions developed for use in the WTP melters. The data will be used in the enhanced property/composition models being developed for waste glass vitrification and melter operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗