Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Line Rating”

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 217 records · Page 12

Dynamic modeling of a single-stage downward firing, entrained flow gasifier

The gasifier is the heart of the integrated gasification combined cycle (IGCC), a technology that has emerged as an attractive alternative to conventional coal-fired power plant technology due to its higher efficiency and cleaner environmental performance especially with the option of CO{sub 2} capture and sequestration. Understanding the optimal performance of the gasifier is therefore paramount for the efficient operation of IGCC power plants. Numerous gasifier models of varying complexity have been developed to study the various aspects of gasifier performance. These range from simple one-dimensional (1D) process-type models to rigorous higher order 2-3D models based on computational fluid dynamics (CFD). Whereas high-fidelity CFD models can accurately predict most key aspects of gasifier performance, they are computationally expensive and typically take hours to days to execute on high-performance computers. Therefore, faster 1D partial differential equation (PDE)-based models are required for use in dynamic simulation studies, control system analysis, and training applications. A number of 1D gasifier models can be found in the literature, but most are steady-state models that have limited application in the practical operation of the gasifier. As a result, 1D PDE-based dynamic models are needed to further study and predict gasifier performance under a wide variety of process conditions and disturbances. In the current study, a 1D transient model of a single-stage downward-fired GE/Texaco-type entrained-flow gasifier has been developed. The model comprises mass, momentum and energy balances for the gas and solid phases. The model considers the initial gasification processes of water evaporation and coal devolatilization. In addition, the key heterogeneous and homogeneous chemical reactions have been modeled. The resulting time-dependent PDE model is solved using the well-known method of lines approach in Aspen Custom Modeler®, whereby the PDEs in the spatial domain are discretized and the resulting differential algebraic equations (DAEs) are then integrated over time using a dynamic integrator. The dynamic response results of the gasifier performance parameters to certain disturbances commonly encountered during practical operation are presented. These disturbances include ramp and step changes to input variables such as coal flow rate, oxygen-to-coal ratio and water-to-coal ratio among others. Comparison of model predictions to available dynamic data will also be discussed.

01 COAL, LIGNITE, AND PEAT↗

MODEL FOR ALFVEN WAVE TURBULENCE IN SOLAR CORONAL LOOPS: HEATING RATE PROFILES AND TEMPERATURE FLUCTUATIONS

It has been suggested that the solar corona may be heated by dissipation of Alfven waves that propagate up from the solar photosphere. According to this theory, counterpropagating Alfven waves are subject to nonlinear interactions that lead to turbulent decay of the waves and heating of the chromospheric and coronal plasma. To test this theory, better models for the dynamics of Alfven waves in coronal loops are required. In this paper, we consider wave heating in an active region observed with the Solar Dynamics Observatory in 2010 May. First a three-dimensional (3D) magnetic model of the region is constructed, and ten magnetic field lines that match observed coronal loops are selected. For each loop we construct a 3D magnetohydrodynamic model of the Alfven waves near the selected field line. The waves are assumed to be generated by footpoint motions inside the kilogauss magnetic flux elements at the two ends of the loop. Based on such models, we predict the spatial and temporal profiles of the heating along the selected loops. We also estimate the temperature fluctuations resulting from such heating. We find that the Alfven wave turbulence model can reproduce the observed characteristics of the hotter loops in the active region core, but the loops at the periphery of the region have large expansion factors and are predicted to be thermally unstable.

79 ASTRONOMY AND ASTROPHYSICS↗

EXPLORING MAGNETIC FIELD STRUCTURE IN STAR-FORMING CORES WITH POLARIZATION OF THERMAL DUST EMISSION

The configuration and evolution of the magnetic field in star-forming cores are investigated in order to directly compare simulations and observations. We prepare four different initial clouds having different magnetic field strengths and rotation rates, in which magnetic field lines are aligned/misaligned with the rotation axis. First, we calculate the evolution of such clouds from the prestellar stage until long after protostar formation. Then, we calculate the polarization of thermal dust emission expected from the simulation data. We create polarization maps with arbitrary viewing angles and compare them with observations. Using this procedure, we confirmed that the polarization distribution projected on the celestial plane strongly depends on the viewing angle of the cloud. Thus, by comparing the observations with the polarization map predicted by the simulations, we can roughly determine the angle between the direction of the global magnetic field and the line of sight. The configuration of the polarization vectors also depends on the viewing angle. We find that an hourglass configuration of magnetic field lines is not always realized in a collapsing cloud when the global magnetic field is misaligned with the cloud rotation axis. Depending on the viewing angle, an S-shaped configuration of the magnetic field (or the polarization vectors) appears early in the protostellar accretion phase. This indicates that not only the magnetic field but also the cloud rotation affects the dynamical evolution of such a cloud. In addition, by comparing the simulated polarization with actual observations, we can estimate properties of the host cloud such as the evolutionary stage, magnetic field strength, and rotation rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗

Application of the Generic Electronics for Time Projection Chamber ($\mathrm{GET}$) readout system for heavy Radioactive isotope collision experiments

Here we have implemented the Generic Electronics for Time Projection Chamber (GET) in a SAMURAI Pion Reconstruction and Ion-Tracker (SπRIT) readout system for heavy radioactive ion collision experiments at RIKEN-RIBF. The SπRIT experiment is designed for heavy ion collision experiments with radioactive ion beams, where a Time Projection Chamber (TPC) with 12096 pixelized readout pads is employed as the main device. Since the TPC is located on the beam line, the readout electronics must handle small signals from pions as well as very large signals from beam or large fragment particles. Operation of the GET electronics during experiment functioned well using 270 time-bucket readout with 25 MHz sampling at an event Data acquisition (DAQ) rate of 60 Hz. Using the slope information of acquired signals it is possible to extend the dynamic range of dE/dx information compared to using the peak height information. However, huge signals arising from energetic δ-rays produced by un-interacted projectiles induce dead channels, which can be recovered after 70 μs on average.

47 OTHER INSTRUMENTATION↗

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗

Measurements and Prediction of Ash Deposition in a Cyclone-Fired Boiler Operating under Variable Load Conditions

Measurements of ash deposition rates were made between the secondary superheater and reheater sections of a 450 MW cyclone-fired lignite boiler as the operational load varied from 33 to 100%. Significant reductions in deposition rates with a decrease in operational load were observed. To uncover the causative mechanisms behind these observations, operational data from the power plant were used to carry out computational fluid dynamic (CFD) simulations of the boiler. After ascertaining that the gas temperatures and velocities at various sections within the boiler were being represented adequately, decoupled simulations of the ash deposition process on the deposit probe were carried out using a finely resolved boundary layer mesh. Fly ash particle size distribution (PSD) and its concentration for the decoupled calculations were determined from stand-alone cyclone barrel simulations. The ash partitioning (mass %) between the fly ash and slag was found to be ~50:50, which was in line with previous field observations, and it did not vary significantly across different cyclone loads. The predicted PSD of the deposit ash was concentrated in the size range 10–30 microns, which was in agreement with cross-sectional images of the deposit obtained from the measurements. At lower loads, sharp variations in the deposition rates were predicted in the gas temperature range 950–1150 K. The particle kinetic energy—particle viscosity-based capture methodology utilized in this study in conjunction with appropriate ash compositions, ash viscosity models and gas temperature estimates can help estimate slagging propensities at different loads reasonably well in these systems.

Krishnamoorthy, Gautham (ORCID:0000000255205092)↗

Study of the Power Supply Ripple Effect on teh Dynamics at SPEAR

For long term stability analysis time variation of tunes is important. We have proposed and tested a technique for measuring the magnitude of this variation. This was made possible by using tune extraction algorithms that require small number of turns thus giving an instantaneous tune of the machine. In this paper we demonstrate the measured effect of the tune modulation with 60 Hz power supplies ripple, power line interference from SLAC linac operating at 30 Hz repetition rate, and nonperiodic variation.

43 PARTICLE ACCELERATORS↗

Considerations about future hard x-ray area detectors

X-ray sources continue to advance in both intensity and temporal domains, thereby opening new ways to analyze the structure and properties of matter, provided that the resultant x-ray images can be efficiently and quantitatively recorded. In this perspective we focus on specific limitations of pixel area x-ray detectors. Although pixel area x-ray detectors have also advanced in recent years, many experiments are still detector limited. Specifically, there is need for detectors that can acquire successive images at GHz rates; detectors that can accurately measure both single photon and millions of photons per pixel in the same image at frame rates of hundreds of kHz; and detectors that efficiently capture images of very hard x-rays (20 keV to several hundred keV). The data volumes and data rates of state-of-the-art detection exceeds most practical data storage options and readout bandwidths, thereby necessitating on-line processing of data prior to, or in lieu of full frame readouts.

36 MATERIALS SCIENCE↗

RADIO OBSERVATIONS REVEAL UNUSUAL CIRCUMSTELLAR ENVIRONMENTS FOR SOME TYPE Ibc SUPERNOVA PROGENITORS

We present extensive radio observations of the nearby Type Ibc supernovae (SNe Ibc) 2004cc, 2004dk, and 2004gq spanning {Delta}t Almost-Equal-To 8-1900 days after explosion. Using a dynamical model developed for synchrotron emission from a slightly decelerated shock wave, we estimate the velocity and energy of the fastest ejecta and the density profile of the circumstellar medium. The shock waves of all three supernovae are characterized by non-relativistic velocities of v-bar approx. (0.1-25)c and associated energies of E Almost-Equal-To (2-10) Multiplication-Sign 10{sup 47} erg, in line with the expectations for a typical homologous explosion. Smooth circumstellar density profiles are indicated by the early radio data and we estimate the progenitor mass-loss rates to be M-dot approx. (0.6-13) x 10{sup -5} M{sub Sun} yr{sup -1} (wind velocity, v{sub w} = 10{sup 3} km s{sup -1}). These estimates approach the saturation limit ( M-dot {approx}10{sup -4} M{sub Sun} yr{sup -1}) for line-driven winds from Wolf-Rayet stars, the favored progenitors of SNe Ibc including those associated with long-duration gamma-ray bursts. Intriguingly, at later epochs all three supernovae show evidence for abrupt radio variability that we attribute to large density modulations (factor of {approx}3-6) at circumstellar radii of r Almost-Equal-To (1-50) Multiplication-Sign 10{sup 16} cm. If due to variable mass loss, these modulations are associated with progenitor activity on a timescale of {approx}10-100 years before explosion. We consider these results in the context of variable mass-loss mechanisms including wind clumping, metallicity-independent continuum-driven ejections, and binary-induced modulations. It may also be possible that the SN shock waves are dynamically interacting with wind termination shocks; however, this requires the environment to be highly pressurized and/or the progenitor to be rapidly rotating prior to explosion. The proximity of the density modulations to the explosion sites may suggest a synchronization between unusual progenitor mass loss and the SN explosion, reminiscent of Type IIn supernovae. This study underscores the utility of radio observations for tracing the final evolutionary stage(s) of SN progenitor systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Megahertz Rate Optical Diagnostics of Explosively Generated Soot

Detonation of a solid explosive produces a turbulent and luminous post‐detonation fireball containing condensed carbon soot. Diagnostics of soot dynamics are needed for model validation and to interpret emission signals. Diffuse back‐illumination extinction imaging (DBI‐EI) and laser‐induced incandescence (LII) are two common optical diagnostics for flame soot. This work extends both to measure time‐resolved soot dynamics from a 12 mm HMX hemisphere. DBI‐EI measured line‐of‐sight optical extinction exceeds 99% in some regions. Results are tomographically inverted to obtain a signal proportional to the soot volume fraction. Pulse‐burst LII measures are presented at 1 MHz. For the first time, LII results are combined with DBI‐EI extinction measurements to quantify and correct for signal trapping. Following this, spatially and temporally resolved DBI‐EI and LII measures are shown to be in reasonable agreement. Finally, experimental results are compared against recently developed simulations. Quantitative differences in the soot mixing dynamics are resolved. These findings motivate future model improvements and demonstrate ongoing needs for diagnostic advancements for heavily sooting environments.

diffuse back-illumination extinction imaging↗

Correlation between the gas-phase metallicity and ionization parameter in extragalactic H II regions

The variations of the metallicity and ionization parameter in H II regions are usually thought to be the dominant factors that produce the variations we see in the observed emission line spectra. There is an increasing amount of evidence that these two quantities are physically correlated, although the exact form of this correlation is debatable in the literature. Simulated emission line spectra from photoionized clouds provide important clues about the physical conditions of H II regions and are frequently used for deriving metallicities and ionization parameters. Through a systematic investigation on the assumptions and methodology used in applying photoionization models, we find that the derived correlation has a strong dependence on the choice of model parameters. On the one hand, models that give consistent predictions over multiple emission-line ratios yield a positive correlation between the metallicity and ionization parameter for the general population of H II regions or star-forming galaxies. On the other hand, models that are inconsistent with the data locus in high-dimensional line ratio space yield discrepant correlations when different subsets of line ratios are used in the derivation. The correlation between the metallicity and ionization parameter has a secondary dependence on the surface density of the star formation rate (SFR), with the higher SFR regions showing a higher ionization parameter but weaker correlations. The existence of the positive correlation contradicts the analytical wind-driven bubble model for H II regions. We explore assumptions in both dynamical models and photoionization models, and conclude that there is a potential bias associated with the geometry. However, this is still insufficient to explain the correlation. Mechanisms that suppress the dynamical influence of stellar winds in realistic H II regions might be the key to solving this puzzle, though more sophisticated combinations of dynamical models and photoionization models to test are required.

79 ASTRONOMY AND ASTROPHYSICS↗

Low temperature spectral dynamics of single molecules in ultrathin polymer films

We studied the spectral dynamics of single fluorescent dye molecules embedded in ultrathin films (5 – 100 nm) of the amorphous polymer polyisobutylene at cryogenic temperatures and its variation with film thickness. Noticeable portion of molecules in the ensemble shows a behavior which is inconsistent with the standard tunneling model: Their spectral lines are subject to irreversible spectral jumps, continuous shifting, and abrupt chaotic changes of the linewidth or jumping rate. In films thinner than 100 nm, the occurrence of “non-standard” spectral behavior increases with decreasing sample thickness at fixed excitation intensity. In addition, it also increases with laser intensity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

Ultrahigh-resolution spectroscopy with atomic or molecular dark resonances: Exact steady-state line shapes and asymptotic profiles in the adiabatic pulsed regime

Exact and asymptotic line shape expressions are derived from the semiclassical density matrix representation describing a set of closed three-level {Lambda} atomic or molecular states including decoherences, relaxation rates, and light shifts. An accurate analysis of the exact steady-state dark-resonance profile describing the Autler-Townes doublet, the electromagnetically induced transparency or coherent population trapping resonance, and the Fano-Feshbach line shape leads to the linewidth expression of the two-photon Raman transition and frequency shifts associated to the clock transition. From an adiabatic analysis of the dynamical optical Bloch equations in the weak field limit, a pumping time required to efficiently trap a large number of atoms into a coherent superposition of long-lived states is established. For a highly asymmetrical configuration with different decay channels, a strong two-photon resonance based on a lower states population inversion is established when the driving continuous-wave laser fields are greatly unbalanced. When time separated resonant two-photon pulses are applied in the adiabatic pulsed regime for atomic or molecular clock engineering, where the first pulse is long enough to reach a coherent steady-state preparation and the second pulse is very short to avoid repumping into a new dark state, dark-resonance fringes mixing continuous-wave line shape properties and coherent Ramsey oscillations are created. Those fringes allow interrogation schemes bypassing the power broadening effect. Frequency shifts affecting the central clock fringe computed from asymptotic profiles and related to the Raman decoherence process exhibit nonlinear shapes with the three-level observable used for quantum measurement. We point out that different observables experience different shifts on the lower-state clock transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global Carbon Budget 2024

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC) are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The global net uptake of CO2 by the ocean (SOCEAN, called the ocean sink) is estimated with global ocean biogeochemistry models and observation-based fCO2 products (fCO2 is the fugacity of CO2). The global net uptake of CO2 by the land (SLAND, called the land sink) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2023, EFOS increased by 1.3 % relative to 2022, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (10.3 ± 0.5 GtC yr−1 when the cement carbonation sink is not included), and ELUC was 1.0 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 ± 0.9 GtC yr−1 (40.6 ± 3.2 GtCO2 yr−1). Also, for 2023, GATM was 5.9 ± 0.2 GtC yr−1 (2.79 ± 0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 2.3 ± 1.0 GtC yr−1, with a near-zero BIM (−0.02 GtC yr−1). The global atmospheric CO2 concentration averaged over 2023 reached 419.31 ± 0.1 ppm. Preliminary data for 2024 suggest an increase in EFOS relative to 2023 of +0.8 % (−0.2 % to 1.7 %) globally and an atmospheric CO2 concentration increase by 2.87 ppm, reaching 422.45 ppm, 52 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2023, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink. This living-data update documents changes in methods and datasets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2024 (Friedlingstein et al., 2024).

Friedlingstein, Pierre (ORCID:0000000333094739)↗

Dynamic Hydride Vapor Phase Epitaxy (D-HVPE) as a Route to Inexpensive, High-Performance III-V Materials and Devices

This work describes the recent development of Dynamic Hydride Vapor Phase Epitaxy (D-HVPE), which extends the standard HVPE process to allow for the deposition of relatively complex heterostructures and devices at growth rates in excess of 300 um/h. This process uses low-cost elemental metal precursors, and uses all of the reactants efficiently to create high-performance devices. D-HVPE also lends itself to potential in-line growth, which when combined with the low-cost precursors and extremely fast growth rates for high-quality III-V materials, can lead to significantly lower epitaxial growth costs. Here, we focus on the growth of high-efficiency photovoltaics as an example device to show the capabilities of the D-HVPE growth process.

14 SOLAR ENERGY↗

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗