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 739 records · Page 41

Engineering scale demonstration of a prospective Cast Stone process

This report documents an engineering-scale demonstration with non-radioactive simulants that was performed at SRNL using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. Over three days, the SCPF was used to fill a 1600 gallon container, staged outside the facility, with simulated Cast Stone grout. The container, staged outside the building approximately 60 ft from the SCPF, was instrumented with x-, y-, and z-axis thermocouples to monitor curing temperature. The container was also fitted with two formed core sampling vials. For the operation, the targeted grout production rate was 1.5 gpm. This required a salt solution flow rate of approximately 1 gpm and a premix feed rate of approximately 580 lb/h. During the final day of operation, the dry feed rate was increased to evaluate the ability of the system to handle increased throughput. Although non-steady state operational periods created free surface liquids, no bleed water was observed either before or after operations. The final surface slope at a fill height of 39.5 inches was 1-1.5 inches across the 8.5 foot diameter container, highest at the final fill point and lowest diametrically opposed to the fill point. During processing, grout was collected in cylindrical containers from both the mixer discharge and the discharge into the container. These samples were stored in a humid environment either in a closed box proximal to the container or inside the laboratory. Additional samples collected at these sampling points were analyzed for rheological properties and density. Both the rheological properties (plastic viscosity and yield strength) and density were consistent with previous and later SCPF runs.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗

Effect of spatial density variation and O+ concentration on the growth and evolution of electromagnetic ion cyclotron waves

We simulate electromagnetic ion cyclotron (EMIC) wave growth and evolution within three regions, the plasmasphere (or plasmaspheric plume), the plasmapause, and the low-density plasmatrough outside the plasmapause. First, we use a ring current simulation with a plasmasphere model to model the particle populations that give rise to the instability for conditions observed on 9 June 2001. Then, using two different models for the cold ion composition, we do a full scale hybrid code simulation in dipole coordinates of the EMIC waves on a meridional plane at MLT = 18 and at 1900 UT within a range of L shell from L = 4.9 to 6.7. EMIC waves were observed during June 9, 2001 by Geostationary Operational Environmental Satellite (GOES) spacecraft. While an exact comparison between observed and simulated spectra is not possible here, we do find significant similarities between the two, at least at one location within the region of largest wave growth. We find that the plasmapause is not a preferred region for EMIC wave growth, though waves can grow in that region. The density gradient within the plasmapause does, however, affect the orientation of wave fronts and wave vector both within the plasmapause and in adjacent regions. There is a preference for EMIC waves to be driven in the He+ band (frequencies between the O+ and He+ gyrofrequencies) within the plasmasphere, although they can also grow in the plasmatrough. If present, H+ band waves are more likely to grow in the plasmatrough. This fact, plus L dependence of the frequency and possible time evolution toward lower frequency waves, can be explained by a simple model. Large O+ concentration limits the frequency range of or even totally quenches EMIC waves. This is more likely to occur in the plasmatrough at solar maximum. Such large O+ concentration significantly affects the H+ cutoff frequency and hence the width in frequency of the stop band above the He+ gyrofrequency. EMIC wave surfaces predicted by cold plasma theory are altered by the finite temperature of the ring current H+.

58 GEOSCIENCES↗

Analysis of Monolith Cores from an Engineering Scale Demonstration of a Prospective Cast Stone Process

The primary disposition path of Low Activity Waste (LAW) at the DOE Hanford Site is vitrification. A cementitious waste form is one of the alternatives being considered for the supplemental immobilization of the LAW that will not be treated by the primary vitrification facility. Washington River Protection Solutions (WRPS) has been directed to generate and collect data on cementitious or pozzolanic waste forms such as Cast Stone. This report documents the coring and leach testing of monolithic samples cored from an engineering-scale demonstration (ES Demo) with non-radioactive simulants. The ES Demo was performed at SRNL in October of 2013 using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft. diameter x 3.25 ft. high container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average LAW composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. In 2014 core samples originally obtained approximately six months after filling the ES Demo were tested along with bench scale molded samples that were collected during the original pour. A latter set of core samples were obtained in late March of 2015, eighteen months after completion of the original ES Demo. Core samples were obtained using a 2” diameter x 11” long coring bit. The ES Demo was sampled in three different regions consisting of an outer ring, a middle ring and an inner core zone. Cores from these three lateral zones were further segregated into upper, middle and lower vertical segments. Monolithic core samples were tested using the Environmental Protection Agency (EPA) Method 1315, which is designed to provide mass transfer rates (release rates) of inorganic analytes contained in monolithic material under diffusion controlled release conditions as a function of leaching time. Compressive strength measurements and drying tests were also performed on the 2015 samples. Leachability indices reported are based on analyte concentrations determined from dissolution of the dried samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Molecular Modeling of Surfactant Interaction on Phospholipid Bilayers Mimicking Corneal Epithelium

Surfactants found in consumer products can compromise eye corneal membrane integrity upon accidental exposure. Traditional in vitro and in vivo approaches to evaluate membrane–surfactant interaction pose experimental limitations such as species variability, reproducibility, and most often do not provide the overall picture. These limitations motivate the use of in silico models to study phenomena like cellular disruption assays caused by surfactants at the molecular scale. In this work, coarse-grained molecular dynamics simulations have been employed to investigate how nonionic alcohol ethoxylate (AE) and anionic surfactant alcohol ethoxy sulfate (AES) interact with lipid bilayer liposomes that mimic corneal epithelial cell membranes. The spherical liposome is composed of 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phospho-l-serine (DOPS), and cholesterol, resembling the composition of the corneal epithelial cells’ membrane bilayer. The simulation consisted of varying degrees of representative surfactant compositions and two initial types of surfactant configurations within or outside the liposome. Our results reveal that both surfactants induce outer leaflet bulging, agreeing with membrane solubilization models. The more highly ethoxylated surfactant, AE, caused more consistent inner leaflet disruption than AES, resulting in significantly more water permeation and membrane thinning. In addition, both surfactants increase the lateral diffusion of lipids within the membrane layers, with higher ethoxylated AE showing a stronger effect than AES. This study demonstrates how surfactant structure and localization influence bilayer membrane integrity, offering mechanistic insights into the irritation potential, thus guiding the rational design of effective surfactant-based formulations.

Lipids↗

SU-E-T-380: Particle Microdosimetry Study Based On 3D-Cylindrical Silicon Radiaton Detectors

Purpose: A new design of a solid-state-microdetector based on silicon 3D microfabrication and its performance to characterise Lineal energy, Specific Energy, dose, LET and other microdosimetric variables required for modelling particle relative biological effectiveness (RBE) is presented. Methods: A microdosimeter formed by a matrix of independent sensors with well-defined micrometric cylindrical shape and with a volume similar to those of cellular dimensions is used to measure microdosimetric variables. Each sensor measures the radiation deposited energy which, divided by the mean cord length of the sensors, provides us with the Linear Energy (y) of the radiation as well as its energy distribution, and frequencymean. Starting from the these distributions in different points of a proton beam, we generate biophysical data (e.g. Linear Energy Transfer (LET), Specific Energy (z), etc…) needed for relative biological effectiveness (RBE) calculations radiation effect models used in particle radiotherapy treatment planning. In addition, a Tissue Equivalent Proportional Counter (TEPC) will be used as baseline to calibrate the “y” magnitude of the microdosimeter unit-cells. Results: The experimental measurements will soon be carried out at the Perelman Center for Advanced Medicine (University of Pennsylvania), which provides proton beam for clinical research proposals. The results of distributions measured of the microdosimetric variables from the first tests developed in the proton facility will be presented and compared with Monte Carlo simulations using the Geant4 code. Conclusion: The use of 3D microdosimeters such as the one presented here will enhance the accuracy of RBE calculations normally affected by the inherent uncertainty of monte carlo simulations due to the approximation of material composition and energy dependent physical laws involved in such calculations. The effect of such approximations will be quatified by comparison with absolute measurement of radiation quality parameters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗

SENSITIVITY ANALYSIS OF FUEL PHYSICAL PROPERTY EFFECTS ON SPARK IGNITION ENGINE PERFORMANCE

Alternative fuels are of interest to automakers and regulators due to their potential to reduce net greenhouse gas emissions from transportation sources. Alternative fuels also have fuel properties which may enable advanced combustion modes with higher engine thermal efficiencies. There has been previous work to identify the relationship between various fuel properties and engine performance, but most of this work has been experiments or simulations where the change in properties was obtained through changing the fuel composition, making isolating the effects of individual fuel properties difficult. In this study, numerical simulations have been used to investigate the effects of individual fuel physical properties such as viscosity or heat of vaporization (HoV) on engine performance. Simulations have been performed of two different engine platforms, the first an optical, single-cylinder research engine and the second a multi-cylinder production engine. Both engines are direct-injection spark-ignition engines with pent-roof heads and are designed for automotive applications. Each engine was run at a different operating condition, one stable and one knock-limited. Different base fuels provided a variety of simulated conditions. Up to six different fuel properties were varied as part of Global Sensitivity Analyses performed for each of the engines with multiple performance targets including thermal efficiency, combustion efficiency and combustion phasing. Results show trends that are largely consistent with previous experimental findings using multiple fuels. The engine thermal efficiency was primarily sensitive to the fuel's HoV, with other fuel physical properties having smaller effects. For optical engine results, the magnitude of the effect was greater in this study than expected based on previous experimental results were many fuel physical and chemical properties were varied simultaneously. However, for the multi-cylinder production engine, the relationship between thermal efficiency and HoV was slightly smaller.

engine simulation↗

Quantitative interpretation of molecular dynamics simulations for X-ray photoelectron spectroscopy of aqueous solutions

Over the past decade, energy-dependent ambient pressure X-ray photoelectron spectroscopy (XPS) has emerged as a powerful analytical probe of the ion spatial distributions at the vapor (vacuum)-aqueous electrolyte interface. These experiments are often paired with complementary molecular dynamics (MD) simulations in an attempt to provide a complete description of the liquid interface. There is, however, no systematic protocol that permits a straightforward comparison of the two sets of results. XPS is an integrated technique that averages signals from multiple layers in a solution even at the lowest photoelectron kinetic energies routinely employed, whereas MD simulations provide a microscopic layer-by-layer description of the solution composition near the interface. Here, we use the National Institute of Standards and Technology database for the Simulation of Electron Spectra for Surface Analysis (SESSA) to quantitatively interpret atom-density profiles from MD simulations for XPS signal intensities using sodium and potassium iodide solutions as examples. We show that electron inelastic mean free paths calculated from a semi-empirical formula depend strongly on solution composition, varying by up to 30% between pure water and concentrated NaI. The XPS signal thus arises from different information depths in different solutions for a fixed photoelectron kinetic energy. XPS signal intensities are calculated using SESSA as a function of photoelectron kinetic energy (probe depth) and compared with a widely employed ad hoc method. SESSA simulations illustrate the importance of accounting for elastic-scattering events at low photoelectron kinetic energies (<300 eV) where the ad hoc method systematically underestimates the preferential enhancement of anions over cations. Finally, some technical aspects of applying SESSA to liquid interfaces are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High‐Performance and Stable Semi‐Transparent Perovskite Solar Cells through Composition Engineering

Abstract Semi‐transparent perovskite solar cells (ST‐PeSCs) have tremendous potential as solar windows owing to their higher efficiency and visible transmittance. However, studies toward this application are still nascent, particularly in unraveling the interplay between how the perovskite composition impacts the achievable device performance and stability. Here, the role of A‐ and X‐site modification in APbX 3 perovskites is studied to understand their influence on these factors. Through detailed experimental and simulation work, it is found that a perovskite composition consisting of cesium (Cs) and formamidinium (FA) at the A‐site delivers the best device performance over a range of band gaps, which are tuned by changes to the X‐site anion. Using this optimized perovskite composition, power conversion efficiencies of 15.5% and 4.1% are achieved for ST‐PeSCs with average visible transmittance values between 20.7% and 52.4%, respectively. Furthermore, the CsFA‐based ST‐PeSCs show excellent long‐term stability under continuous illumination and heating. The stability of the precursor solutions across each of the studied compositions has also been considered, showing dramatic differences in the structural properties of the perovskites and their device performance for all mixed A‐site compositions possessing the archetypal methyl ammonium species, while also confirming the superior stability of the CsFA precursor solutions.

14 SOLAR ENERGY↗

Epithermal Neutron Resonance Analysis Using a Compact D-T Generator

Nuclear material identification and isotopic composition analysis are important capabilities for safeguards and the verification of nuclear arms control and nonproliferation treaties. Epithermal neutron resonance imaging has been proposed for spent fuel assay measurements, and the development of a mobile neutron imaging system could provide such a tool for on-site applications. This technique is highly sensitive to isotopes of interest for safeguards applications (e.g. 235,238 U, 238,239,240 Pu) and yields unique transmission spectra for a target containing actinides and some mid- and high-Z elements. Paired with a position-sensitive neutron detector, such a system could provide spatial information of the target isotopic composition. Recent progress in the development of high intensity neutron sources could make mobile epithermal neutron imaging designs increasingly feasible to serve as a tool for nuclear inspectors. In this work, Monte Carlo (MC) simulations were performed to further optimize a previous design for an epithermal neutron resonance-based analysis apparatus using a D-T portable neutron generator. Additionally, experimental work is being performed to test the feasibility of this method. The MC simulations of moderation time distributions for different moderator compositions and thicknesses were measured to provide a multi-group correction for neutron time-of-flight energy reconstruction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A composite electrodynamic mechanism to reconcile spatiotemporally resolved exciton transport in quantum dot superlattices

Quantum dot (QD) solids are promising optoelectronic materials; further advancing their device functionality requires understanding their energy transport mechanisms. The commonly invoked near-field Förster resonance energy transfer (FRET) theory often underestimates the exciton hopping rate in QD solids, yet no consensus exists on the underlying cause. In response, we use time-resolved ultrafast stimulated emission depletion (STED) microscopy, an ultrafast transformation of STED to spatiotemporally resolve exciton diffusion in tellurium-doped cadmium selenide–core/cadmium sulfide–shell QD superlattices. We measure the concomitant time-resolved exciton energy decay due to excitons sampling a heterogeneous energetic landscape within the superlattice. The heterogeneity is quantified by single-particle emission spectroscopy. This powerful multimodal set of observables provides sufficient constraints on a kinetic Monte Carlo simulation of exciton transport to elucidate a composite transport mechanism that includes both near-field FRET and previously neglected far-field emission/reabsorption contributions. Uncovering this mechanism offers a much-needed unified framework in which to characterize transport in QD solids and additional principles for device design.

36 MATERIALS SCIENCE↗

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗

Phase-field modeling of the interactions between an edge dislocation and an array of obstacles

Obstacles, such as voids and precipitates, are prevalent in crystalline materials. They strengthen crystals by serving as barriers to dislocation glide. Here in this work, we develop a phase-field dislocation dynamics (PFDD) technique for investigating the interactions between dislocations and second-phase obstacles, which can be either voids or precipitates. The PFDD technique is constructed to account for elastic heterogeneity, elastic anisotropy, dissociation of the dislocation, and dislocation transmission across bicrystalline interfaces. Within the framework, we present a model for “pseudo-voids”, which are voids shearable by dislocations, in contrast to unphysical, unshearable voids in conventional phase-field dislocation formulations. We employ the PFDD technique to investigate the in-plane interactions between an edge dislocation and an array of nano-scale obstacles with different spacings. In this application, the interactions take place in glide planes of either a face-centered cubic (FCC) Cu or a body-centered cubic (BCC) Nb matrix, while the precipitates have a Cu 1-x Nb x composition, with x varying from 0.1 to 0.9. Our atomistic simulations find that the alloy precipitates can have an FCC, an amorphous, or a BCC phase, depending on the compositional ratio between Cu and Nb, i.e., value of x. Among all types of obstacles, the critical stresses for dislocation bypass are the highest for unshearable amorphous precipitates, followed by shearable crystalline precipitates, and then the pseudo-voids.

42 ENGINEERING↗

Simulation of Spent Nuclear Fuel loading into a Final Waste Repository

Previous work towards studying repository loading have used: - spent nuclear fuel (SNF) with an average burnup composition; - a lumped capacitance thermal model for calculating temperature in a Cyclus repository model. The goal of this work is to improve on the repository models and use U.S. historical SNF inventory data in simulations to more accurately study the loading of a waste repository. These goals will be achieved by: - using UNF-ST-and-DARDS Unified Database (UDB) that has historic assembly-specific data (e.g, isotopic composition, heat) in Cyclus simulations; - implementing a more accurate thermal model within a Cyclus repository model. Objectives: - Create a Cyclus spent fuel conditioning model that packages spent fuel bundles into packages which have user-defined properties. - Create a Cyclus medium-fidelity repository model that gives accurate time and spatial dependent temperature values and loads the repository based on a user-selected loading strategy. Cyclus is an agent-based extensible framework for modeling flow of material through user-defined nuclear fuel cycles. In Cyclus, each facility in the fuel cycle is modeled individually and the facilities interact with one another as independent agents. The spent fuel conditioning model accepts spent fuel bundles and puts them into a cylindrical waste package. In the spent fuel conditioning model, the user can define variables: For each layer, radius, thermal conductivity, thermal diffusivity; For each package, Number of spent fuel bundles, Radius and height. The waste repository model accepts waste packages and emplaces them into available positions within the waste repository based on a thermal criteria. The thermal criteria is a temperature limit at the interface between the waste package surface and the host geology, that is set based on the repository's host geology. In the waste repository model, the user can define the variables: Capacity, Distance between waste packages, Distance between drifts, Repository host geology, Loading Strategy. After the addition of new waste packages at each time step, the waste repository model recalculates the temperature at each location in the repository. If the addition of this new package causes its temperature to exceed the thermal limit, it will be placed back into the buffer. A thermal model that relies on a transient 'outside' model and quasi-steady-state 'inside' model is used to accurately determine the temperature in the repository. The 'outside' model assumes a homogenous medium with the Engineered Barrier System (EBS) replaced by the geologic medium. Temperature solutions for the central waste package, adjacent point and line sources are superimposed to calculate the temperature at specific points in the repository. The central drift consists of one finite line source which represents the central waste package. The central drift also consists of point sources that represent neighboring waste packages in the central drift. The neighboring drifts are represented by infinite line sources. The 'inside' model is considered to be at a quasi-steady-state condition because EBS has a relatively low thermal mass compared to the infinite geologic medium. The steady state calculation is performed at each time step with the heat source and interface temperature as boundary conditions. Future Work: Run Cyclus simulations with U.S. historical SNF inventory data, the spent fuel conditioning and repository models to study how waste package acceptance strategies impact repository loading.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-D-19A-05: The Dosimetric Impact of Using Xoft Axxent® Electronic Brachytherapy Source TG-43 Dosimetry Parameters for Treatment with the Xoft 30 Mm Diameter Vaginal Applicator

Purpose: A full TG-43 dosimetric characterization has not been performed for the Xoft Axxent ® electronic brachytherapy source (Xoft, a subsidiary of iCAD, San Jose, CA) within the Xoft 30 mm diameter vaginal applicator. Currently, dose calculations are performed using the bare-source TG-43 parameters and do not account for the presence of the applicator. This work focuses on determining the difference between the bare-source and sourcein- applicator TG-43 parameters. Both the radial dose function (RDF) and polar anisotropy function (PAF) were computationally determined for the source-in-applicator and bare-source models to determine the impact of using the bare-source dosimetry data. Methods: MCNP5 was used to model the source and the Xoft 30 mm diameter vaginal applicator. All simulations were performed using 0.84p and 0.03e cross section libraries. All models were developed based on specifications provided by Xoft. The applicator is made of a proprietary polymer material and simulations were performed using the most conservative chemical composition. An F6 collision-kerma tally was used to determine the RDF and PAF values in water at various dwell positions. The RDF values were normalized to 2.0 cm from the source to accommodate the applicator radius. Source-in-applicator results were compared with bare-source results from this work as well as published baresource results. Results: For a 0 mm source pullback distance, the updated bare-source model and source-in-applicator RDF values differ by 2% at 3 cm and 4% at 5 cm. The largest PAF disagreements were observed at the distal end of the source and applicator with up to 17% disagreement at 2 cm and 8% at 8 cm. The bare-source model had RDF values within 2.6% of the published TG-43 data and PAF results within 7.2% at 2 cm. Conclusion: Results indicate that notable differences exist between the bare-source and source-in-applicator TG-43 simulated parameters. Xoft Inc. provided partial funding for this work.

60 APPLIED LIFE SCIENCES↗

Oxidation and Condensation of Zinc Fume From Zn-CO 2 -CO-H 2 O Streams Relevant to Steelmaking Off-Gas Systems

Here the objective of this research was to study the condensation of zinc vapor to metallic zinc and zinc oxide solid under varying environments to investigate the feasibility of in-process separation of zinc from steelmaking off-gas dusts. Water vapor content, temperature, degree of cooling, gas composition, and initial zinc partial pressure were varied to simulate the possible conditions that can occur within steelmaking off-gas systems, limited to Zn-CO 2 -CO-H 2 O gas compositions. The temperature of deposition and the effect of rapidly quenching the gas were specifically studied. A homogeneous nucleation model for applicable experiments was applied to the analysis of the experimental data. It was determined that under the experimental conditions, oxidation of zinc vapor by H 2 O or CO 2 does not occur above 1108 K (835 °C) even for highly oxidizing streams (CO 2 /CO = 40/7). Rate expressions that correlate CO 2 and H 2 O oxidation rates to gas composition, partial pressure of water vapor, temperature, and zinc partial pressure were determined to be as follows: Rate $ \left(\frac{mol}{m^2s}\right) $ = 406 exp $ \left(\frac{-50.2 kJ/mol}{RT}\right) $ (pZnpCO 2 $-$ PCO/K eq CO 2 ) $\frac{mol}{m^2 x s}$ Rate $ \left(\frac{mol}{m^2s}\right) $ = 32.9 exp $ \left(\frac{-13.7 kJ/mol}{RT}\right) $ (pZnPH 2 O $-$ PH 2 /K eq H 2 O) $\frac{mol}{m^2 x s}$. It was proven that a rapid cooling rate (500 K/s) significantly increases the ratio of metallic zinc to zinc oxide as opposed to a slow cooling rate (250 K/s). SEM analysis found evidence of heterogeneous growth of ZnO as well as of homogeneous formation of metallic zinc. The homogeneous nucleation model fit well with experiments where only metallic zinc deposited. An expanded model with rates of oxidation by CO 2 and H 2 O as shown was combined with the homogenous nucleation model and then compared with experimental data. The calculated results based on the model gave a reasonable fit to the measured data. For the conditions used in this study, the rate equations for the oxidation of zinc by carbon dioxide and water vapor as well as the homogeneous nucleation model of metallic zinc were applicable for various temperatures, zinc partial pressures, CO 2 :CO ratios, and H 2 O partial pressures.

36 MATERIALS SCIENCE↗