Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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

Reaction-space analysis of homogeneous charge compression ignition combustion with varying levels of fuel stratification under positive and negative valve overlap conditions

Full-cycle computational fluid dynamics simulations with gasoline chemical kinetics were performed to determine the impact of breathing and fuel injection strategies on thermal and compositional stratification, combustion and emissions during homogeneous charge compression ignition combustion. The simulations examined positive valve overlap and negative valve overlap strategies, along with fueling by port fuel injection and direct injection. The resulting charge mass distributions were analyzed prior to ignition using ignition delay as a reactivity metric. The reactivity stratification arising from differences in the distributions of fuel–oxygen equivalence ratio (Φ FO ), oxygen molar fraction (χ O2 ) and temperature (T) was determined for three parametric studies. In the first study, the reactivity stratification and burn duration for positive valve overlap valve events with port fuel injection and early direct injection were nearly identical and were dominated by wall-driven thermal stratification. nitrogen oxide (NO) and carbon monoxide (CO) emissions were negligible for both injection strategies. In the second study, which examined negative valve overlap valve events with direct injection and port fuel injection, reactivity stratification increased for direct injection as the Φ FO and T distributions associated with direct fuel injection into the hot residual gas were positively correlated; however, the latent heat absorbed from the hot residual gas by the evaporating direct injection fuel jet reduced the overall thermal and reactivity stratification. These stratification effects were offsetting, resulting in similar reactivity stratification and burn durations for the two injection strategies. The higher local burned gas temperatures with direct injection resulted in an order of magnitude increase in NO, while incomplete combustion of locally over-lean regions led to a sevenfold increase in CO emissions compared to port fuel injection. The final study evaluated positive valve overlap and negative valve overlap valve events with direct injection. Furthermore, relative to positive valve overlap, the negative valve overlap condition had a wider reactivity stratification, a longer burn duration and higher NO and CO emissions associated with reduced fuel–air mixing.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

High-Entropy Alloys for Accelerator Beam Window Applications

Development of novel high-entropy alloys (HEAs) is currently underway for potential use as beam windows in future multi-megawatt target systems at Fermilab. HEAs encompass a new class of materials with a vast design space allowing for material properties to be tailored for particular applications and to potentially offer improved resistance to beam-induced radiation damage and thermal shock effects. The alloy systems being studied consist of several compositions of AlCoCrMnTiV with 4 6 component elements. These alloys are all predicted by CALPHAD simulation to have a single-phase BCC crystal structure and low density, with some compositions displaying ordered, nanoscale precipitates. This presentation will briefly discuss alloy design and synthesis before giving a detailed description of the characterization studies of these HEAs in both the pristine state and post-irradiation by low-energy heavy ions to high damage levels. Electron microscopy techniques to quantify elemental homogeneity and composition, determine grain size, shape, and orientation, and quantify lattice parameters, defect structures and precipitate phases are all being used to study alloy microstructures. Mechanical properties of the alloys at the microscale will be reported. The evolution of these properties as a function of radiation damage will also be described. Bulk thermal characteristics of these HEAs have been tested to measure specific heat capacity and coefficient of thermal expansion as a function of temperature. To determine bulk tensile properties a miniature tensile testing apparatus is under development; it s commissioning will be covered briefly. The talk will conclude with our plans for alloy down-selection.

Burleigh, A. [Fermilab]↗
DOE OSTI2014

Non-destructive quantification of alignment of nanorods embedded in uniaxially stretched polymer films

Among several methods developed for uniaxial alignment of metallic nanorods for optical applications, alignment by film stretching consists in embedding the rods in a transparent thin film of thermoplastic polymer, followed by simultaneous heating and uniaxial stretching of the composite film. As to the quantification of the resulting alignment, it has been limited to statistical calculations based on microscopic examination, which is incomplete, subject to errors due to geometric distortions of the scanning electron microscope images and destructive, since it involves cutting of samples. In contrast, we present in this paper a non-destructive quantification of the average orientation of the rods, based on a probabilistic approach combined with numerical simulations of absorbance spectra and spectrometric characterization of the composite film. Assuming electromagnetically non-interacting rods, we consider the longitudinal absorbance peak of their ensemble to consist of the superposition of their individual spectra that we obtain by numerical simulation using the size and shape adapted dielectric function of the metal and the finite difference time domain method. The accuracy of the solution depends on the number of discretization intervals, the accuracy of the numerical simulations, and the accurate knowledge of the polydispersity of the rods. For the sake of concreteness, we used nanorods to describe the quantification steps but the method is equally valid for any dichroic particles.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Mathematical simulation of the amplification of 1790-nm laser radiation in a nuclear-excited He – Ar plasma containing nanoclusters of uranium compounds

The possibility of applying dusty active media in nuclearpumped lasers has been considered. The amplification of 1790-nm radiation in a nuclear-excited dusty He – Ar plasma is studied by mathematical simulation. The influence of nanoclusters on the component composition of the medium and the kinetics of the processes occurring in it is analysed using a specially developed kinetic model, including 72 components and more than 400 reactions. An analysis of the results indicates that amplification can in principle be implemented in an active laser He – Ar medium containing 10-nm nanoclusters of metallic uranium and uranium dioxide. (lasers)

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Peptide-Directed PdAu Nanoscale Surface Segregation: Toward Controlled Bimetallic Architecture for Catalytic Materials

Bimetallic nanoparticles are of immense scientific and technological interest given the synergistic properties observed when mixing two different metallic species at the nanoscale. This is particularly prevalent in catalysis, where bimetallic nanoparticles often exhibit improved catalytic activity and durability over their monometallic counterparts. Yet despite intense research efforts, little is understood regarding how to optimize bimetallic surface composition and structure synthetically using rational design principles. Recently, it has been demonstrated that peptide-enabled routes for nanoparticle synthesis result in materials with sequence-dependent catalytic properties, providing an opportunity for rational design through sequence manipulation. In this study, bimetallic PdAu nanoparticles are synthesized with a small set of peptides containing known Pd and Au binding motifs. The resulting nanoparticles were extensively characterized using high-resolution scanning transmission electron microscopy, X-ray absorption spectroscopy and high-energy X-ray diffraction coupled to atomic pair distribution function analysis. Structural information obtained from synchrotron radiation methods were then used to generate model nanoparticle configurations using reverse Monte Carlo simulations, which illustrate sequence-dependence in both surface structure and surface composition. Replica exchange solute tempering molecular dynamic simulations were also used to predict the modes of peptide binding on monometallic surfaces, indicating that different sequences bind to the metal interfaces via different mechanisms. As a testbed reaction, electrocatalytic methanol oxidation experiments were performed, wherein differences in catalytic activity are clearly observed in materials with identical bimetallic composition. Finally, taken together, this study indicates that peptides could be used to arrive at bimetallic surfaces with enhanced catalytic properties, which could be leveraged for rational bimetallic nanoparticle design using peptide-enabled approaches.

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

Dose response of selected solid state detectors in applied homogeneous transverse and longitudinal magnetic fields

Purpose: MR-Linac devices under development worldwide will require standard calibration, commissioning, and quality assurance. Solid state radiation detectors are often used for dose profiles and percent depth dose measurements. The dose response of selected solid state detectors is therefore evaluated in varying transverse and longitudinal magnetic fields for this purpose. Methods: The Monte Carlo code PENELOPE was used to model irradiation of a PTW 60003 diamond detector and IBA PFD diode detector in the presence of a magnetic field. The field itself was varied in strength, and oriented both transversely and longitudinally with respect to the incident photon beam. The long axis of the detectors was oriented either parallel or perpendicular to the photon beam. The dose to the active volume of each detector in air was scored, and its ratio to dose with zero magnetic field strength was determined as the “dose response” in magnetic field. Measurements at low fields for both detectors in transverse magnetic fields were taken to evaluate the accuracy of the simulations. Additional simulations were performed in a water phantom to obtain few representative points for beam profile and percent depth dose measurements. Results: Simulations show significant dose response as a function of magnetic field in transverse field geometries. This response can be near 20% at 1.5 T, and it is highly dependent on the detectors’ relative orientation to the magnetic field, the energy of the photon beam, and detector composition. Measurements at low transverse magnetic fields verify the simulations for both detectors in their relative orientations to radiation beam. Longitudinal magnetic fields, in contrast, show little dose response, rising slowly with magnetic field, and reaching 0.5%–1% at 1.5 T regardless of detector orientation. Water tank and in air simulation results were the same within simulation uncertainty where lateral electronic equilibrium is present and expectedly differed at the beam edge in transverse field orientations only. Due to the difference in design, the two detectors behaved differently. Conclusions: When transverse magnetic fields are present, great care must be taken when using diamond or diode detectors. Dose response varies with relative detector orientation, magnetic field strength, and between detectors. This response can be considerable (∼20% for both detectors). Both detectors in longitudinal fields exhibit little to no dose response as a function of magnetic field. Water tank simulations seem to suggest that the diode detector is better suited to general beam commissioning, and each detector must be investigated separately.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2015

Ethane-xenon mixtures under shock conditions

Mixtures of light elements with heavy elements are important in inertial confinement fusion. We explore the physics of molecular scale mixing through a validation study of equation of state (EOS) properties. Density functional theory molecular dynamics (DFT-MD) at elevated temperature and pressure is used to obtain the thermodynamic state properties of pure xenon, ethane, and various compressed mixture compositions along their principal Hugoniots. In order to validate these simulations, we have performed shock compression experiments using the Sandia Z-Machine. A bond tracking analysis correlates the sharp rise in the Hugoniot curve with the completion of dissociation in ethane. Furthermore, the DFT-based simulation results compare well with the experimental data along the principal Hugoniots and are used to provide insight into the dissociation and temperature along the Hugoniots as a function of mixture composition. Interestingly, we find that the compression ratio for complete dissociation is similar for several compositions suggesting a limiting compression for C-C bonded systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Hydrogen Isotope Separation from Aqueous Streams using Zeolite Membranes

Efficient separation of heavy hydrogen isotopes {sup 3}H (or T, tritium) and {sup 2}H (or D, deuterium) from {sup 1}H has been a long-standing endeavor in nuclear fuel cycle research and development. For example, tritium (half-life 12.3 years) is a fission product from nuclear plants and is found as tritiated water (HTO) in effluent streams from nuclear installations. Tritium must be concentrated and removed from such aqueous streams for waste storage. The separation of deuterium from deuterated water (HDO) has applications in heavy water production, and can also be used as a nonradioactive simulant in HTO separation research. Several methods have been explored such as ion exchange/adsorption and membrane permeation. Adsorption by molecular sieving materials such as zeolites can separate tritiated water from water with moderate selectivity (HTO/H{sub 2}O separation factor∼1.2-1.4) due to a subtle difference in adsorption characteristics between the two species. However, there remains a requirement for high-selectivity continuous separation of heavy hydrogen isotopes from water. Here we demonstrate new results on the use of zeolite membranes as a potential high-performance, low-cost, radiation-stable candidate for such separations. In zeolite membranes, separation of HTO (or HDO) from water would take place by a combination of adsorption and diffusion characteristics. Here we describe the hydrogen isotope separation properties of silicoaluminophosphate (SAPO-34) membranes synthesized on porous α-alumina supports. D{sub 2}O/H{sub 2}O feed mixtures of realistic compositions are used as a non-radioactive simulant for tritiated water feeds. Our results show promising applications of zeolite membranes for hydrogen isotope separation and recovery. In summary, SAPO-34 membranes were fabricated and investigated to demonstrate their pervaporative separation properties for D{sub 2}O-H{sub 2}O binary systems. These membrane exhibit high selectivity (up to ∼19) for D{sub 2}O and good permeance (> 600 GPU) even at low concentrations of D{sub 2}O in the feed, making them highly interesting candidates for capture of heavy hydrogen isotopes (D and T) from aqueous streams encountered in the nuclear fuel cycle. The mechanism of separation appears to involve both diffusivity and adsorption differences as expected in nanoporous materials, even though detailed identification of their respective contributions is a subject of future work. We also hypothesize that tritiated (T{sub 2}O/H{sub 2}O) mixtures will also show similar performance in relation to the present deuterated mixtures. This is planned to be further investigated through pervaporation measurements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Investigation of variable compositions on the removal of technetium from Hanford Waste Treatment Plant low activity waste melter off-gas condensate simulant

The Low Activity Waste (LAW) vitrification facility at the Hanford Waste Treatment and Immobilization Plant (WTP) will generate an aqueous condensate recycle stream (LAW Off-Gas Condensate) from the offgas system. The plan for disposition of this stream during baseline operations is to send it to the WTP Pretreatment Facility, where it will be blended with LAW, concentrated by evaporation and recycled to the LAW vitrification facility again. The primary reason to recycle this stream is so that the semi-volatile 99 Tc isotope eventually becomes incorporated into the glass. This stream also contains non-radioactive salt components that are problematic in the melter, so diversion of this stream to another process would eliminate recycling of these salts and would enable simplified operation of the LAW melter and the Pretreatment Facilities. This diversion from recycling this stream within WTP would have the effect of decreasing the LAW vitrification mission duration and quantity of glass waste. The concept being tested here involves removing the 99 Tc so that the decontaminated aqueous stream, with the problematic salts, can be disposed elsewhere.

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

U-net architected deep material network training with microstructure local field information

The Deep Material Network (DMN) has recently emerged as a powerful reduced-order modeling framework for simulating the mechanical response of heterogeneous materials such as composites. Unlike most data-driven approaches that directly learn a material’s response under prescribed loading, the DMN acts as a homogenization operator, learning the kinematic constraints and mechanical interactions of the underlying microstructure. However, traditional DMN training relies exclusively on homogenized effective properties derived from Direct Numerical Simulations (DNS), discarding the rich local field data that govern microstructural interactions. In this work, we extend the DMN framework to incorporate such local field information into the offline training process. Utilizing a U-Net architecture, we augment the DMN training objective to include the first and second statistical moments of the local stress fields obtained from linear DNS. This ensures that the learned network topology not only fits the effective stiffness but also accurately reflects the internal local stress and strain partitioning of the microstructure. The results confirm that supervising the localization process during training yields a superior surrogate model, reducing local prediction errors by an order of magnitude and significantly improving generalization to unseen nonlinear constitutive behaviors compared to traditional DMNs.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Cinema image-based in situ analysis and visualization of MPAS-ocean simulations

We created an in situ exploration visualization of an MPAS-Ocean simulation.We leveraged compositing in Cinema to provide interactive exploration.We decreased the storage footprint of the analysis and visualization results. Due to power and I/O constraints associated with extreme scale scientific simulations, in situ analysis and visualization will become a critical component to scientific exploration and discovery. Current analysis and visualization options at extreme scale are presented in opposition: write files to disk for interactive, exploratory analysis, or perform in situ analysis to save data products about phenomena that a scientists knows about in advance. In this study, we demonstrate extreme scale visualization of MPAS-Ocean simulations leveraging a third option based on Cinema, which is a novel framework for highly interactive, image-based in situ analysis and visualization that promotes exploration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2014

Angular dependence of the response of the nanoDot OSLD system for measurements at depth in clinical megavoltage beams

Purpose: The purpose of this investigation was to assess the angular dependence of a commercial optically stimulated luminescence dosimeter (OSLD) dosimetry system in MV x-ray beams at depths beyondd{sub max} and to find ways to mitigate this dependence for measurements in phantoms. Methods: Two special holders were designed which allow a dosimeter to be rotated around the center of its sensitive volume. The dosimeter's sensitive volume is a disk, 5 mm in diameter and 0.2 mm thick. The first holder rotates the disk in the traditional way. It positions the disk perpendicular to the beam (gantry pointing to the floor) in the initial position (0°). When the holder is rotated the angle of the disk towards the beam increases until the disk is parallel with the beam (“edge on,” 90°). This is referred to as Setup 1. The second holder offers a new, alternative measurement position. It positions the disk parallel to the beam for all angles while rotating around its center (Setup 2). Measurements with five to ten dosimeters per point were carried out for 6 MV at 3 and 10 cm depth. Monte Carlo simulations using GEANT4 were performed to simulate the response of the active detector material for several angles. Detector and housing were simulated in detail based on microCT data and communications with the manufacturer. Various material compositions and an all-water geometry were considered. Results: For the traditional Setup 1 the response of the OSLD dropped on average by 1.4% ± 0.7% (measurement) and 2.1% ± 0.3% (Monte Carlo simulation) for the 90° orientation compared to 0°. Monte Carlo simulations also showed a strong dependence of the effect on the composition of the sensitive layer. Assuming the layer to completely consist of the active material (Al{sub 2}O{sub 3}) results in a 7% drop in response for 90° compared to 0°. Assuming the layer to be completely water, results in a flat response within the simulation uncertainty of about 1%. For the new Setup 2, measurements and Monte Carlo simulations found the angular dependence of the dosimeter to be below 1% and within the measurement uncertainty. Conclusions: The dosimeter system exhibits a small angular dependence of approximately 2% which needs to be considered for measurements involving other than normal incident beams angles. This applies in particular to clinicalin vivo measurements where the orientation of the dosimeter is dictated by clinical circumstances and cannot be optimized as otherwise suggested here. When measuring in a phantom, the proposed new setup should be considered. It changes the orientation of the dosimeter so that a coplanar beam arrangement always hits the disk shaped detector material from the thin side and thereby reduces the angular dependence of the response to within the measurement uncertainty of about 1%. This improvement makes the dosimeter more attractive for clinical measurements with multiple coplanar beams in phantoms, as the overall measurement uncertainty is reduced. Similarly, phantom based postal audits can transition from the traditional TLD to the more accurate and convenient OSLD.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Insights from the Lattice-Strain Evolution on Deformation Mechanisms in Metallic-Glass-Matrix Composites

In this paper, in situ high-energy synchrotron X-ray diffraction experiments and micromechanics-based finite element simulations have been conducted to examine the lattice-strain evolution in metallic-glass-matrix composites (MGMCs) with dendritic crystalline phases dispersed in the metallic-glass matrix. Significant plastic deformation can be observed prior to failure from the macroscopic stress–strain curves in these MGMCs. The entire lattice-strain evolution curves can be divided into elastic–elastic (denoting deformation behavior of matrix and inclusion, respectively), elastic–plastic, and plastic–plastic stages. Characteristics of these three stages are governed by the constitutive laws of the two phases (modeled by free-volume theory and crystal plasticity) and geometric information (crystalline phase morphology and distribution). The load-partitioning mechanisms have been revealed among various crystalline orientations and between the two phases, as determined by slip strain fields in crystalline phase and by strain localizations in matrix. Finally, implications on ductility enhancement of MGMCs are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗
DOE OSTI2016

Calcium Phosphate Constraints in HLW Glasses: Phosphate Liquidus Model

About 50 million gallons of high-level mixed waste is currently stored in underground tanks at the United States Department of Energy’s (DOE’s) Hanford site in the State of Washington. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) will provide DOE’s Office of River Protection (ORP) with a means of treating this waste by vitrification for subsequent disposal. The tank waste will be separated into low- and high-activity waste fractions, which will then be vitrified respectively into Immobilized Low Activity Waste (ILAW) and Immobilized High Level Waste (IHLW) products. The ILAW product will be disposed in an engineered facility on the Hanford site while the IHLW product is designed for acceptance into a national geological disposal facility for high-level nuclear waste. The ILAW and IHLW products must meet a variety of requirements before they can be accepted for disposal. Glass property-composition models will be used to control the vitrification process at the WTP to ensure that the glass meets product quality and processability requirements. The quality and the validity ranges of the property-composition models and the supporting database will greatly influence waste treatment operations, including the waste loading that is achievable in the glass for each waste composition. The current baseline models for operation of the WTP were developed under the Bechtel National, Inc. (BNI) WTP contract. These models were based on LAW and HLW glass property data over a range of compositions that was guided primarily by the WTP contract and the set of tanks specified therein. The HLW testing was directed towards waste streams from the then-planned early feed tanks for the WTP (i.e., AZ-101, AZ-102, C-106/AY-102, and C-104/AY-101). These wastes are high in iron (AZ-101, AZ-102, and C-106/AY-102) or thorium (C-104/AY-101). While the baseline models may be adequate for commissioning of the facility and processing of the specified tanks, refinements and extensions or augmentations of these models will be required to process the balance of the tank waste inventory. The overall mission of the WTP will require the immobilization of tank waste compositions that are variously dominated by mixtures of aluminum, chromium, bismuth, iron, phosphorous, zirconium, and sulfur compounds as waste-loading-limiting components. Further, the WTP contract includes minimum waste loading requirements that are generally rather modest while loadings significantly higher than those in the current baseline are likely achievable. In order to achieve these higher waste loadings and to expand the applicable composition regions, new models will be needed to support the optimized operation of the WTP. Since the development of the baseline WTP models, ORP directed a program of work at the Vitreous State Laboratory (VSL) to examine various approaches for optimization of the performance of the WTP, including increasing the waste loadings in the LAW and HLW glass products. That work involved HLW rich in Al, Al+Na, Cr, Bi, S, Zr, Sr, P and Fe, in contrast to old baseline models that were based mostly on HLW tanks that were high in Fe, with one other that was also high in Th. Glass property data collected under the ORP program were compiled and combined with data collected to develop the WTP models. The combined databases have been used to update pre-existing HLW composition models or to develop new models for properties such as Product Consistency Test (PCT) responses, spinel crystallization (one-percent crystal fraction temperature, T 1% ), melt viscosity, melt electrical conductivity, nepheline formation upon slow cooling (i.e., Canister Centerline Cooling [CCC]), and Toxicity Characteristic Leach Procedure (TCLP) Cd response. While spinel is typically the most commonly occurring phase in HLW glasses, other phases can be limiting depending on the waste and glass composition. For example, various phosphate phases can become important as the phosphorus content increases. Lithium phosphate crystallization can impact PCT releases of HLW glasses, while calcium phosphate phases forming in the cold-cap can deleteriously impact feed processing in the melter. At present, however, no models exist that relate glass composition and crystallization of phosphate phases in general and calcium phosphate in particular. The present work was therefore directed towards the initial development of such models. Future work may be directed at studying amorphous phase separation in high phosphate HLW glasses. Within the River Protection Project, feed compositions are projected using the Hanford Tank Waste Operations Simulator (HTWOS), a dynamic computer model that tracks the waste through storage, retrieval, feed staging, and treatment processes. A review of the HTWOS data shows that WTP will process 4106 batches of HLW glass for the overall mission (System Plan Rev. 6, ORP-11242), each with its own waste feed and glass composition. Some of the projected HTWOS feed batches are expected to have relatively high concentrations of phosphorous compared to typical HLW streams. In the 4106 batches of HLW for the overall WTP mission P 2 O 5 varies from 0.59 wt% to 13.00 wt%, with 95 batches containing at least 10 wt% P 2 O 5 . Occurrences of phosphate crystals have commonly been noted in glasses with high bismuth and phosphorus contents. In addition, recently collected data from a set of statistically designed HLW matrix glasses prepared and characterized to support enhancement of the HLW glass property-composition models identified apatite in a number of glass samples. However, the current HLW model for crystallization (one-percent crystal fraction temperature) is applicable only to spinel crystallization in glasses. Previous projections of waste loadings and glass volume (such as those in System Plan Rev. 6, ORP-11242) have used preliminary constraints for calcium phosphate crystallization that are based on the limited available data. Additional data and improved constraints or, preferably, models that relate glass composition to crystallization of phosphate-containing phases, are therefore needed to support HLW vitrification at the WTP. The present study focuses on the development of such a model based on data from previously characterized glasses and a set of newly fabricated glasses.

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

Full-wave simulation of a three-dimensional metamaterial prism

In our article, a negative-index metamaterial prism based on a composite unit cell containing a split-ring resonator and a z-dipole is designed and simulated. The design approach combines simulations of a single unit cell to identify the appropriate cell design (yielding the desired negative-index behavior) together with subcell modeling (which simplifies the mesh representation of the resonator geometry and allows for a larger number of resonator cells to be handled). Furthermore, to describe the methodology used in designing a n = -1 refractive index prism, our results include the effective-medium parameters, the far-field scattered patterns, and the near-zone field distributions corresponding to a normally incident plane-wave excitation of the prism.

42 ENGINEERING↗