Results for “Dynamic Line Rating” — Page 10

DOE OSTI2015

Spontaneous toroidal flow generation due to negative effective momentum diffusivity

Spontaneous structure formation, and in particular, zonal flows, is observed in a broad range of natural and engineered systems, often arising dynamically as the saturated state of a linear instability. Flows in tokamaks are known to self-organise on small scales, but large scale toroidal flows also arise even when externally applied torques are zero. This has previously been interpreted as the result of small externally imposed breaking of a symmetry. However, we show that for large enough field line pitch, a robust spontaneous symmetry breaking occurs, leading to the generation of strong toroidal flow structures; parameters are typical of Spherical Tokamak discharges with reversed shear profiles. The short wavelength dynamics are qualitatively similar to the growth of poloidal flow structures, and toroidal flow gradients nonlinearly saturate at levels where the shearing rate is comparable to linear growth rate. On long wavelengths, we measure Prandtl numbers of around zero for these systems, in conjunction with the formation of structured toroidal flows, and we show that this is consistent with a model of momentum transport where fluxes act to reinforce small flow gradients: the effective momentum diffusivity is negative. Toroidal flow structures are largely unaffected by collisional damping, so this may allow toroidal bulk flows of order the ion thermal velocity to be maintained with zero momentum input. This phenomenon also provides a mechanism for the generation of localised meso-scale structures like transport barriers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

INFRARED VARIABILITY OF EVOLVED PROTOPLANETARY DISKS: EVIDENCE FOR SCALE HEIGHT VARIATIONS IN THE INNER DISK

We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks in the IC 348 cluster that show significant infrared variability. Using 3-8 {mu}m and 24 {mu}m photometry along with 5-40 {mu}m spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5 {mu}m spectroscopy, optical spectroscopy, and near-infrared photometry, covering timescales of days to years, we examine the variability in the disk, stellar, and accretion flux. We find substantial variations (10%-60%) at all infrared wavelengths on timescales of weeks to months for all of these young stellar objects. This behavior is not unique when compared to other cluster members and is consistent with changes in the structure of the inner disk, most likely scale height fluctuations on a dynamical timescale. Previous observations, along with our near-infrared photometry, indicate that the stellar fluxes are relatively constant; stellar variability does not appear to drive the large changes in the infrared fluxes. Based on our near-infrared spectroscopy of the Pa{beta} and Br{gamma} lines we find that the accretion rates are variable in most of the evolved disks but the overall rates are probably too small to cause the infrared variability. We discuss other possible physical causes for the variability, including the influence of a companion, magnetic fields threading the disk, and X-ray flares.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Evaporation of Water Nanodroplets on Heated Surfaces: Does Nano Matter?

While experiments and continuum models have provided a relatively good understanding of the evaporation of macroscopic water droplets, elucidating how sessile nanodroplets evaporate is an open question critical for advancing nanotechnological applications where nanodroplets can play an essential role. In this study, using molecular dynamics simulations, we find that evaporating nanodroplets, in contrast to their macroscopic counterparts, are not always in thermal equilibrium with the substrate and that the vapor concentration on the nanodroplet surface does not reach a steady state. As a result, the evaporative behavior of nanodroplets is significantly different. Regardless of hydrophobicity, nanodroplets do not follow conventional evaporation modes but instead exhibit dynamic wetting behavior characterized by huge, non-equilibrium, isovolumetric fluctuations in the contact angle and contact radius. For hydrophilic nanodroplets, the evaporation rate, controlled by the vapor concentration, decays exponentially over time. Hydrophobic nanodroplets follow stretched exponential kinetics arising from the slower thermalization with the substrate. The evaporative half-lifetime of the nanodroplets is directly related to the thermalization time scale and therefore increases monotonically with the hydrophobicity of the substrate. Finally, the evaporative flux profile along the nanodroplet surface is highly nonuniform but does not diverge at the contact line as the macroscopic continuum models predict.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2011

CHARACTERIZING EXOPLANETARY ATMOSPHERES THROUGH INFRARED POLARIMETRY

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties, and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g., due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1%, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line of sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Pyruvate, phosphate dikinase regulatory protein impacts light response of C4 photosynthesis in Setaria viridis

In C4 plants, the pyruvate (Pyr), phosphate dikinase regulatory protein (PDRP) regulates the activity of the C4 pathway enzyme Pyr, phosphate dikinase (PPDK) in a light-/dark-dependent manner. The importance of this regulatory action to C4 pathway function and overall C4 photosynthesis is unknown. To resolve this question, we assessed in vivo PPDK phospho-regulation and whole leaf photophysiology in a CRISPR–Cas9 PDRP knockout (KO) mutant of the NADP-ME C4 grass green millet (Setaria viridis). PDRP enzyme activity was undetectable in leaf extracts from PDRP KO lines. Likewise, PPDK phosphorylated at the PDRP-regulatory Thr residue was immunologically undetectable in leaf extracts. PPDK enzyme activity in rapid leaf extracts was constitutively high in the PDRP KO lines, irrespective of light or dark pretreatment of leaves. Gas exchange analysis of net CO2 assimilation revealed PDRP KO leaves had markedly slower light induction kinetics when leaves transition from dark to high-light or low-light to high-light. In the initial 30 min of the light induction phase, KO leaves had an ~15% lower net CO2 assimilation rate versus the wild-type (WT). Despite the impaired slower induction kinetics, we found growth and vigor of the KO lines to be visibly indistinguishable from the WT when grown in normal air and under standard growth chamber conditions. However, the PDRP KO plants grown under a fluctuating light regime exhibited a gradual multi-day decline in Fv/Fm, indicative of progressive photosystem II damage due to the absence of PDRP. Collectively, our results demonstrate that one of PDRP’s functions in C4 photosynthesis is to ensure optimal photosynthetic light induction kinetics during dynamic changes in incident light.

Plant Sciences↗

DOE OSTI2026

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on Banshee Distribution Network

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on the Banshee Distribution Network: Preprint

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2014

Heating-induced glass-glass and glass-liquid transformations in computer simulations of water

Water exists in at least two families of glassy states, broadly categorized as the low-density (LDA) and high-density amorphous ice (HDA). Remarkably, LDA and HDA can be reversibly interconverted via appropriate thermodynamic paths, such as isothermal compression and isobaric heating, exhibiting first-order-like phase transitions. We perform out-of-equilibrium molecular dynamics simulations of glassy water using the ST2 model to study the evolution of LDA and HDA upon isobaric heating. Depending on pressure, glass-to-glass, glass-to-crystal, glass-to-vapor, as well as glass-to-liquid transformations are found. Specifically, heating LDA results in the following transformations, with increasing heating pressures: (i) LDA-to-vapor (sublimation), (ii) LDA-to-liquid (glass transition), (iii) LDA-to-HDA-to-liquid, (iv) LDA-to-HDA-to-liquid-to-crystal, and (v) LDA-to-HDA-to-crystal. Similarly, heating HDA results in the following transformations, with decreasing heating pressures: (a) HDA-to-crystal, (b) HDA-to-liquid-to-crystal, (c) HDA-to-liquid (glass transition), (d) HDA-to-LDA-to-liquid, and (e) HDA-to-LDA-to-vapor. A more complex sequence may be possible using lower heating rates. For each of these transformations, we determine the corresponding transformation temperature as function of pressure, and provide a P-T “phase diagram” for glassy water based on isobaric heating. Our results for isobaric heating dovetail with the LDA-HDA transformations reported for ST2 glassy water based on isothermal compression/decompression processes [Chiu et al., J. Chem. Phys. 139, 184504 (2013)]. The resulting phase diagram is consistent with the liquid-liquid phase transition hypothesis. At the same time, the glass phase diagram is sensitive to sample preparation, such as heating or compression rates. Interestingly, at least for the rates explored, our results suggest that the LDA-to-liquid (HDA-to-liquid) and LDA-to-HDA (HDA-to-LDA) transformation lines on heating are related, both being associated with the limit of kinetic stability of LDA (HDA)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

The Kinetics of the Liquid Phase Nucleation in a Stretched FCC Crystal: A Molecular Dynamics Simulation

The kinetics of a spontaneous formation of liquid phase in a stretched (superheated) Lennard-Jones fcc crystal is studied. Molecular dynamics experiments are used to determine the main parameters of the nucleation process: nucleation frequency J, diffusion coefficient of nuclei D{sub *}, nonequilibrium Zel’dovich factor Z, and critical nucleus size R{sub *}. The calculations are performed at negative pressures from the endpoint of the melting line and at positive pressures that are higher by a factor of eight than the critical pressure. The simulation results are compared to the classical homogeneous nucleation theory. It is found that the theory qualitatively correctly reproduces the dynamics of developing the process. The theory and the simulation demonstrate good quantitative agreement for the transition rate of the liquid phase nucleus through the critical size, but there is large difference in the numbers of critical nuclei in the unit volume of the metastable phase. In the case of significant superheatings and negative pressures, the contribution of the energy of elastic stresses to the moving force of the phase transformation is small and it can be neglected in a first approximation. The mismatch between the theory and the simulation results can be eliminated taking that the surface free energy of a curved “crystal–liquid droplet” interface is smaller than that of a plane interface by 30–35%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Formation of Power-law Electron Energy Spectra in Three-dimensional Low- β Magnetic Reconnection

While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of nonrelativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in nonrelativistic low-β reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index p ~ 4 in the 3D simulation, which persist throughout the nonlinear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant acceleration rate for electrons at different energies. Here, we construct a model that explains the observed power-law spectral index in terms of the dynamical balance between particle acceleration and escape from main acceleration regions, which are defined based upon a threshold for the curvature drift acceleration term. This result could be important for explaining the formation of power-law energy spectrum in solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Effect of surface free energies on the heterogeneous nucleation of water droplet: A molecular dynamics simulation approach

Heterogeneous nucleation of water droplet on surfaces with different solid-liquid interaction intensities is investigated by molecular dynamics simulation. The interaction potentials between surface atoms and vapor molecules are adjusted to obtain various surface free energies, and the nucleation process and wetting state of nuclei on surfaces are investigated. The results indicate that near-constant contact angles are already established for nano-scale nuclei on various surfaces, with the contact angle decreasing with solid-liquid interaction intensities linearly. Meanwhile, noticeable fluctuation of vapor-liquid interfaces can be observed for the nuclei that deposited on surfaces, which is caused by the asymmetric forces from vapor molecules. The formation and growth rate of nuclei are increasing with the solid-liquid interaction intensities. For low energy surface, the attraction of surface atoms to water molecules is comparably weak, and the pre-existing clusters can depart from the surface and enter into the bulk vapor phase. The distribution of clusters within the bulk vapor phase becomes competitive as compared with that absorbed on surface. For moderate energy surfaces, heterogeneous nucleation predominates and the formation of clusters within bulk vapor phase is suppressed. The effect of high energy particles that embedded in low energy surface is also discussed under the same simulation system. The nucleation preferably initiates on the high energy particles, and the clusters that formed on the heterogeneous particles are trapped around their original positions instead of migrating around as that observed on smooth surfaces. This feature makes it possible for the heterogeneous particles to act as fixed nucleation sites, and simulation results also suggest that the number of nuclei increases monotonously with the number of high energy particles. The growth of nuclei on high energy particles can be divided into three sub-stages, beginning with the formation of a wet-spot, increase of contact angle with near-constant contact line, and finally growth with constant contact angle. The growth rate of nuclei also increases with the size of high energy particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

High Speed In-situ X-ray Imaging of 3D Freeze Printing of Aerogels

3D freeze printing (3DFP) combines drop-on-demand (DOD) inkjet printing with freeze casting to fabricate lightweight and multifunctional aerogels with customized geometries. Freeze casting is an efficient and easily implemented method capable of fabricating porous, sponge-like structures for many different applications. This process enables tailoring the microstructure of the final product (i.e., pore morphology, alignment, average size distribution, etc.) by controlling the fabrication conditions and freezing kinetics. Furthermore, its combination with DOD printing provides the capability of engineering the macrostructure without relying on a mold as reported for 3D freeze-printed aerogels made from graphene, silver nanowires, and other nanocomposites. In this paper, we performed in-situ X-ray imaging to understand the inside process dynamics in 3DFP using a commercially available colloidal silica ink. We investigated the 3DFP process with the following hierarchy: first, single droplets; then, uniform lines obtained from coalescence of droplets; and finally, three consecutive lines deposited layer by layer. With the help of X-ray imaging, the importance of the balance between material deposition and freezing rates was shown in-situ by the observation inside of the freeze front following the tip of the printed line. The effects of the substrate temperature on the elimination of undesired interfacial boundaries were also shown by the observed ice crystals penetrating from lower to upper layer.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Role of Metal Cocatalysts in the Photocatalytic Production of Hydrogen from Water Revisited

The use of photocatalysts to promote the production of molecular hydrogen from water, following the so-called water splitting reaction, continues to be a promising route for the green production of fuels. The molecular basis of this photocatalysis is the photoexcitation of electrons from the valence band of semiconductors to their conduction band, from which they can be transferred to chemical reactants, protons in the case of water, to promote a reduction reaction. The mechanism by which such a process takes place has been studied extensively using titanium oxide, a simple material that fulfills most requirements for water splitting. However, photocatalysis with TiO 2 tends to be highly inefficient; a cocatalyst, commonly a late transition metal (Au, Pt) in nanoparticle form, needs to be added to facilitate the production of H 2 . The metal is widely believed to help with the scavenging of the excited electrons from the conduction band of the semiconductor in order to prevent their recombination with the accompanying hole formed in the valence band, a step that cancels the initial photon absorption and competes with the photolytic chemical reduction. Here we review and analyze the molecular basis for that mechanism and argue for an alternative explanation, that the role of the metal is to help with the recombination of the atomic hydrogen atoms produced by proton reduction on the semiconductor surface instead. First, we summarize what is known about the electronic structure of these photocatalysts and how the electronic levels need to line up for the reduction of protons in water to be feasible. Next, we review the current understanding of the dynamics of the steps associated with the absorption of photons, the de-excitation via electron–hole pair recombination and fluorescence decay, and the electronic transitions that lead to proton reduction, and contrast those with the rates of the chemical steps required to produce molecular hydrogen. The following section addresses the changes introduced by the addition of the metal cocatalyst, comparatively evaluating its role as either an electron scavenger or a promoter of the recombination of hydrogen atoms. A discussion of the viable chemical mechanisms for the latter pathway is included. Finally, we briefly mention other associated aspects of this photocatalysis, including the possible promotion of H 2 production with visible light via resonant excitation of the surface plasmon of Au nanoparticles, the use of single-metal (Au, Pt) atom catalysts and of yolk–shell nanostructures, and the reduction of organic molecules. We end with a brief personal perspective on the possible generality of the concepts introduced in this Critical Review.

08 HYDROGEN↗

DOE OSTI2013

Electrode-plasma-driven radiation cutoff in long-pulse, high-power microwave devices

The impact of electrode plasma dynamics on the radiation production in a high power microwave device is examined using particle-in-cell simulations. Using the design of a compact 2.4 GHz magnetically insulated line oscillator (MILO) as the basis for numerical simulations, we characterize the time-dependent device power and radiation output over a range of cathode plasma formation rates. These numerical simulations can self-consistently produce radiation characteristics that are similar to measured signals in long pulse duration MILOs. This modeling capability should result in improved assessment of existing high-power microwave devices and lead to new designs for increased radiation pulse durations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Bubbling Behavior in a Waste Glass Melter

Legacy tank waste from five decades of plutonium production at the Hanford Site during the Cold War will be vitrified into a stable waste form for disposal. A Waste Treatment Plant, currently under construction, houses Joule heated ceramic-lined waste glass melters wherein radioactive waste will be combined with glass formers to produce borosilicate glass. Through testing conducted in pilot-scale melters with surrogate waste, as well as operational experience at other sites, it has been demonstrated that the use of bubblers can dramatically improve the glass production rate. However, knowledge of the bubble dynamics in the molten glass is necessary for prediction and design of improved melting processes. A computational fluid dynamics model of the melter has been developed to investigate the processes affecting feed melt rate. The model shows that the convection currents within the melt pool due to Joule heating supplied by the electrodes are enhanced by forced air bubbling. In this paper, parameters influencing the formation and rise of a bubble stream in a viscous liquid are investigated. Experimental data and observations from laboratory and pilot-scale testing are integrated with first principles to develop a set of equations that form a mechanistic bubbling model to be incorporated into the melter model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗