Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Flow”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Origins of heat transport anisotropy in MoTe2 and other bulk van der Waals materials

Knowledge of how heat flows anisotropically in van der Waals (vdW) materials is crucial for thermal management of emerging 2D materials devices and design of novel anisotropic thermoelectric materials. Despite the importance, anisotropic heat transport in vdW materials is yet to be systematically studied and is often presumably attributed to anisotropic speeds of sound in vdW materials due to soft interlayer bonding relative to covalent in-plane networks of atoms. In this work, we investigate the origins of the anisotropic heat transport in vdW materials, through time-domain thermoreflectance (TDTR) measurements and first-principles calculations of anisotropic thermal conductivity of three different phases of MoTe 2 . MoTe 2 is ideal for the study due to its weak anisotropy in the speeds of sound. We find that even when the speeds of sound are roughly isotropic, the measured thermal conductivity of MoTe 2 along the c-axis is 5–8 times lower than that along the in-plane axes. We derive meaningful characteristic heat capacity, phonon group velocity, and relaxation times from our first principles calculations for selected vdW materials (MoTe 2 , BP, h-BN, and MoS 2 ), to assess the contributions of these factors to the anisotropic heat transport. Interestingly, we find that the main contributor to the heat transport anisotropy in vdW materials is anisotropy in heat capacity of the dominant heat-carrying phonon modes in different directions, which originates from anisotropic optical phonon dispersion and disparity in the frequency of heat-carrying phonons in different directions. Further, the discrepancy in frequency of the heat-carrying phonons also leads to ~2 times larger average relaxation times in the cross-plane direction, and partially explains the apparent dependence of the anisotropic heat transport on the anisotropic speeds of sound. This work provides insight into understanding of the anisotropic heat transport in vdW materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

200 h of discharge cycling with an all-aqueous copper thermally regenerative ammonia battery

Thermally regenerative ammonia batteries (TRABs) offer an approach to energy storage and electricity generation by harnessing low-grade heat (T < 150 °C). A TRAB discharge resembles that of a flow battery, where electrochemical reactions produce electrical power from energy stored within aqueous electrolytes. While there are many types of TRABs, the all-aqueous copper TRAB (Cu aq -TRAB) has produced the largest power and energy storage densities. Despite many improvements to TRAB performance, most tests have only lasted a few hours, which are not representative of operation times expected of these devices. Herein, we operated a Cu aq -TRAB for 200 h of constant current discharging to assess battery performance and component stability. After 200 h of testing, the average power density increased slightly to 7.2 mW cm -2 , which was within 0.01 % of starting conditions. Likewise, the average energy density for the final cycle was only 0.03 % lower than the initial cycle. The overall insensitivity of the power cell to cycling represents a major milestone in the advancement of TRAB systems. In conclusion, energy dispersive X-ray spectroscopy provides evidence to suggest that these small changes in power and energy density are likely to the membrane acclimating to TRAB electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2025

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗
DOE OSTI2023

Analysis techniques for blob properties from gas puff imaging data

Filamentary structures, also known as blobs, are a prominent feature of turbulence and transport at the edge of magnetically confined plasmas. They cause cross-field particle and energy transport and are, therefore, of interest in tokamak physics and, more generally, nuclear fusion research. Several experimental techniques have been developed to study their properties. Among these, measurements are routinely performed with stationary probes, passive imaging, and, in more recent years, Gas Puff Imaging (GPI). In this work, we present different analysis techniques developed and used on 2D data from the suite of GPI diagnostics in the Tokamak à Configuration Variable, featuring different temporal and spatial resolutions. Although specifically developed to be used on GPI data, these techniques can be employed to analyze 2D turbulence data presenting intermittent, coherent structures. We focus on size, velocity, and appearance frequency evaluation with, among other methods, conditional averaging sampling, individual structure tracking, and a recently developed machine learning algorithm. We describe in detail the implementation of these techniques, compare them against each other, and comment on the scenarios to which these techniques are best applied and on the requirements that the data must fulfill in order to yield meaningful results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Westinghouse Accident Tolerant Fuel Phase 2B with Higher Enriched and Higher Burnup Add-On Project Final Technical Report Deliverable Volume 1

The Westinghouse Electric Company LLC (Westinghouse) accident tolerant fuel (ATF) program funded by Westinghouse and the US Department of Energy (DOE) utilized chromium (Cr) coated zirconium alloy cladding with doped UO 2 (ADOPT(TM)) and uranium silicide (U 3 Si 2 ) high density/high thermal conductivity fuel for its lead test rod (LTR) program with irradiation beginning in 2019. ADOPT is a Cr 2 O 3 +Al 2 O 3 doped UO 2 pellet with increased oxidation resistance and density, and increased resistance to fission gas release. Due to the issues with U 3 Si 2 fuel reaction in pressurized water reactor (PWR) environments, uranium nitride (U 15 N) has been substituted for U 3 Si 2 as a long-term fuel option. Cr coated cladding with ADOPT fuel is the near-term Westinghouse EnCore fuel product. The long-term EnCore ® fuel product is SiGA ® SiC/SiC composite cladding with high density/high thermal conductivity UN fuel. In 2020, the higher burnup and higher enriched and (HBHE) program was integrated into the ATF program. The objective of this expanded program is to extend ATF burnups to at least 75 MWd/kgU to economically facilitate 24-month cycles in PWRs. The ATF program now includes ADOPT fuel with enrichments >5% 235 U. Over the past several years, Westinghouse has tested the Cr coated zirconium (Zr) and silicon carbide (SiC) claddings in Westinghouse Churchill autoclaves and the Massachusetts Institute of Technology (MIT) reactor. Additionally, the Cr coated cladding has been tested in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and in LTR programs at the Doel 4 and Bryon 2 commercial reactors. High temperature tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried out to determine the time and temperature limits for the Cr coated zirconium claddings. These tests indicate that Cr coated Zr cladding can take temperatures up to about 1500°C for short periods of time without becoming totally oxidized. The main issue is the formation of the Cr-Zr eutectic at 1333°C resulting in the migration of this eutectic inward with the formation of ZrCr 2 which rapidly oxidizes to ZrO 2 on the outside of the tube. The Cr coated Zr cladding also delays the bursting of the tube, reduces the burst area, and reduces hydriding and embrittlement of the Zr which is a major benefit for reducing fuel fragmentation, redistribution, and dispersal (FFRD), the major licensing issue faced by HBHE. The manufacturing parameters for the SiC have been found to have a significant effect on the corrosion rate of the SiC in light water reactor (LWR) conditions and significant improvements have been made in the yield, quality, and oxidation resistance of the SiC cladding by identifying and tightening manufacturing process parameters. Current autoclave results for SiC composite claddings indicate that a corrosion rate of fewer than 2 micrometers per year can be achieved which meets corrosion requirements under normal operating conditions. High temperature steam oxidation tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried and determined that SiC cladding can withstand temperatures up to about 1800°C to 1900°C before excessive corrosion reactions begin, and that the SiC will not balloon and burst. ADOPT fuel pellets were incorporated into the Westinghouse ATF program because of their increased density and oxidation resistance. A topical report supporting the use of ADOPT pellets has been submitted to the NRC and will be approved shortly. ADOPT pellets have been a product in Europe for over 15 years. The additives work to make larger grain sizes which appear to reduce fission gas release during transients as well as increase the density of the pellet. Fuel rod and assembly design in preparation for the lead test rod (LTR) and lead test assembly (LTA) programs is underway as well as licensing efforts with the Nuclear Regulatory Commission (NRC). Finally, accident analyses coupled with economic evaluations for both operating savings as well as fuel savings have been initiated. Modular Accident Analysis Program 5 (MAAP5) calculations indicate that both Cr coated cladding and SiC composite cladding can increase the time to fuel rod loss of geometry for up to two hours longer than current Zr based cladding in a station blackout scenario. This additional time is due to the much higher oxidation resistance of the SiC and Cr coated cladding. These two hours can be used to implement additional response options instituted through the Diverse and Flexible Mitigation Capability (FLEX) program by reactor operators. Both ATF cladding options also reduce hydrogen production which dramatically reduces primary system and containment pressure and the risk of fission product release beyond containment in the unlikely event of an accident. The lower pressure in the system allows more time to feed cooling water to the core, resulting in the avoidance of fuel melting. The coping time is extended indefinitely if the modest water flow provided by FLEX continues. Experimental work on methods to rapidly fabricate SiC composite structures with high density and reduce the fabrication price of SiC fibers while maintaining a high level of performance is needed. Methods for mitigating or stopping the Cr-Zr eutectic are also needed to further improve the oxidation resistance of Cr coated Zr. Minimal (<1%) swelling of U 3 Si 2 and subsequent fission gas release has been demonstrated up to 20 MWd/kgU. Irradiation experiments with U 3 Si 2 fuel in ATR to determine these properties at 40 MWd/kgU were completed but post irradiation examinations were not done since U 3 Si 2 has been replaced by UN. UN has three issues that are being addressed. The first is increasing the oxidation resistance so that there is not excessive reaction up to ~1500°C. While UN increases the pellet density and additives to the UN postponed the temperature at which rapid oxidation occurred by ~100°C to 200°C, this is not enough, and pellet coating options are now being pursued. The second is developing a method to manufacture that does not require the multi-step process of UF 6 to UO 2 to UC to UN. Efforts are underway looking at UF 6 to UN 2 to UN and UF 6 to UF 4 to UN 2 to UN reactions. Finally, an economically acceptable process for potentially enriching the 15 N content of the nitrogen used to make UN to >95% 15N was identified though not experimentally validated. In addition to the work supported by the DOE, research and testing activities are being carried on in a world-wide effort funded by many countries such as Sweden, United Kingdom, Belgium, Netherlands, Spain, Germany, Japan, and France. This work is being facilitated through the Westinghouse led Collaboration for Advanced Research on Accident Tolerant Fuel (CARAT) program. Annual meetings were organized in 2018 and 2019 by Westinghouse as a venue for presentation of this work and to provide for the cross-fertilization of ideas among the many outstanding researchers in the ATF area. No meetings were held in 2020 and 2021 due to the COVID-19 related restrictions on traveling and gatherings. Since SiC, coated cladding, and high-density fuel options are not currently used in the nuclear industry; support from the government and industry members is needed to further the significant effort of setting new standards. The same is true for the Nuclear Regulatory Commission (NRC) which must review and approve the commercial use of these new fuels since all current regulations are oriented toward Zr/UO 2 fuel. Several meetings have been held with the NRC to generate a fast-track approach to licensing ATF using a combination of atomic modeling, in-rod sensors, and in-reactor testing. This approach was memorialized in an “Accelerated Fuel Qualification White Paper” prepared by the Accelerated Fuel Qualification Working Group lead by General Atomics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Westinghouse Accident Tolerant Fuel Phase 2B with Higher Enriched and Higher Burnup Add-On Project Final Technical Report Deliverable Volume 2

The Westinghouse Electric Company LLC (Westinghouse) accident tolerant fuel (ATF) program funded by Westinghouse and the US Department of Energy (DOE) utilized chromium (Cr) coated zirconium alloy cladding with doped UO 2 (ADOPT(TM)) and uranium silicide (U 3 Si 2 ) high density/high thermal conductivity fuel for its lead test rod (LTR) program with irradiation beginning in 2019. ADOPT is a Cr 2 O 3 +Al 2 O 3 doped UO 2 pellet with increased oxidation resistance and density, and increased resistance to fission gas release. Due to the issues with U 3 Si 2 fuel reaction in pressurized water reactor (PWR) environments, uranium nitride (U 15 N) has been substituted for U 3 Si 2 as a long-term fuel option. Cr coated cladding with ADOPT fuel is the near-term Westinghouse EnCore ® fuel product. The long-term EnCore fuel product is SiGA ® SiC/SiC composite cladding with high density/high thermal conductivity UN fuel. In 2020, the higher burnup and higher enriched and (HBHE) program was integrated into the ATF program. The objective of this expanded program is to extend ATF burnups to at least 75 MWd/kgU to economically facilitate 24-month cycles in PWRs. The ATF program now includes ADOPT fuel with enrichments >5% 235 U. Over the past several years, Westinghouse has tested the Cr coated zirconium (Zr) and silicon carbide (SiC) claddings in Westinghouse Churchill autoclaves and the Massachusetts Institute of Technology (MIT) reactor. Additionally, the Cr coated cladding has been tested in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and in LTR programs at the Doel 4 and Bryon 2 commercial reactors. High temperature tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried out to determine the time and temperature limits for the Cr coated zirconium claddings. These tests indicate that Cr coated Zr cladding can take temperatures up to about 1500°C for short periods of time without becoming totally oxidized. The main issue is the formation of the Cr-Zr eutectic at 1333°C resulting in the migration of this eutectic inward with the formation of ZrCr 2 which rapidly oxidizes to ZrO 2 on the outside of the tube. The Cr coated Zr cladding also delays the bursting of the tube, reduces the burst area, and reduces hydriding and embrittlement of the Zr which is a major benefit for reducing fuel fragmentation, redistribution, and dispersal (FFRD), the major licensing issue faced by HBHE. The manufacturing parameters for the SiC have been found to have a significant effect on the corrosion rate of the SiC in light water reactor (LWR) conditions and significant improvements have been made in the yield, quality, and oxidation resistance of the SiC cladding by identifying and tightening manufacturing process parameters. Current autoclave results for SiC composite claddings indicate that a corrosion rate of fewer than 2 micrometers per year can be achieved which meets corrosion requirements under normal operating conditions. High temperature steam oxidation tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried and determined that SiC cladding can withstand temperatures up to about 1800°C to 1900°C before excessive corrosion reactions begin, and that the SiC will not balloon and burst. ADOPT fuel pellets were incorporated into the Westinghouse ATF program because of their increased density and oxidation resistance. A topical report supporting the use of ADOPT pellets has been submitted to the NRC and will be approved shortly. ADOPT pellets have been a product in Europe for over 15 years. The additives work to make larger grain sizes which appear to reduce fission gas release during transients as well as increase the density of the pellet. Fuel rod and assembly design in preparation for the lead test rod (LTR) and lead test assembly (LTA) programs is underway as well as licensing efforts with the Nuclear Regulatory Commission (NRC). Finally, accident analyses coupled with economic evaluations for both operating savings as well as fuel savings have been initiated. Modular Accident Analysis Program 5 (MAAP5) calculations indicate that both Cr coated cladding and SiC composite cladding can increase the time to fuel rod loss of geometry for up to two hours longer than current Zr based cladding in a station blackout scenario. This additional time is due to the much higher oxidation resistance of the SiC and Cr coated cladding. These two hours can be used to implement additional response options instituted through the Diverse and Flexible Mitigation Capability (FLEX) program by reactor operators. Both ATF cladding options also reduce hydrogen production which dramatically reduces primary system and containment pressure and the risk of fission product release beyond containment in the unlikely event of an accident. The lower pressure in the system allows more time to feed cooling water to the core, resulting in the avoidance of fuel melting. The coping time is extended indefinitely if the modest water flow provided by FLEX continues. Experimental work on methods to rapidly fabricate SiC composite structures with high density and reduce the fabrication price of SiC fibers while maintaining a high level of performance is needed. Methods for mitigating or stopping the Cr-Zr eutectic are also needed to further improve the oxidation resistance of Cr coated Zr. Minimal (<1%) swelling of U 3 Si 2 and subsequent fission gas release has been demonstrated up to 20 MWd/kgU. Irradiation experiments with U 3 Si 2 fuel in ATR to determine these properties at 40 MWd/kgU were completed but post irradiation examinations were not done since U 3 Si 2 has been replaced by UN. UN has three issues that are being addressed. The first is increasing the oxidation resistance so that there is not excessive reaction up to ~1500°C. While UN increases the pellet density and additives to the UN postponed the temperature at which rapid oxidation occurred by ~100°C to 200°C, this is not enough, and pellet coating options are now being pursued. The second is developing a method to manufacture that does not require the multi-step process of UF 6 to UO 2 to UC to UN. Efforts are underway looking at UF 6 to UN 2 to UN and UF 6 to UF 4 to UN 2 to UN reactions. Finally, an economically acceptable process for potentially enriching the 15 N content of the nitrogen used to make UN to >95% 15 N was identified though not experimentally validated. In addition to the work supported by the DOE, research and testing activities are being carried on in a world-wide effort funded by many countries such as Sweden, United Kingdom, Belgium, Netherlands, Spain, Germany, Japan, and France. This work is being facilitated through the Westinghouse led Collaboration for Advanced Research on Accident Tolerant Fuel (CARAT) program. Annual meetings were organized in 2018 and 2019 by Westinghouse as a venue for presentation of this work and to provide for the cross-fertilization of ideas among the many outstanding researchers in the ATF area. No meetings were held in 2020 and 2021 due to the COVID-19 related restrictions on traveling and gatherings. Since SiC, coated cladding, and high-density fuel options are not currently used in the nuclear industry; support from the government and industry members is needed to further the significant effort of setting new standards. The same is true for the Nuclear Regulatory Commission (NRC) which must review and approve the commercial use of these new fuels since all current regulations are oriented toward Zr/UO 2 fuel. Several meetings have been held with the NRC to generate a fast-track approach to licensing ATF using a combination of atomic modeling, in-rod sensors, and in-reactor testing. This approach was memorialized in an “Accelerated Fuel Qualification White Paper” prepared by the Accelerated Fuel Qualification Working Group lead by General Atomics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Influences of Mg Ions and Al Ions on Deposition Behavior of Supersaturated Silicic Acid - 20261

Silicic acid can change into various species such as deposit, colloidal, and soluble form under certain conditions: for example, gradual pH changes around geological disposal repositories. In particular, silicic acid deposits narrow groundwater flow paths, thereby restricting the flow of groundwater. Such effects may contribute to a delay in radionuclide migration. Previous studies have reported that silicic acid supersaturated in sodium and calcium ions can narrow the flow paths. This study examined the effect of magnesium (Mg) and aluminum (Al) ions on the deposition of supersaturated silicic acid on the solid phase. The deposition behavior of supersaturated silicic acid on the solid phase was investigated in a batch system while changing the initial supersaturation concentration, the coexisting ion concentration, and the amount of added solid phase. In the presence of Mg ions, silicic acid precipitation increased with increasing Mg ion concentration due to the influence of the electrolyte, and crosslinking reactions. The coexistence of Al ions significantly suppressed silicic acid deposition due to the formation of a covering on the solid surface during Al ion hydrolysis under conditions of high Al concentration (1.4 mM). On the basis of these experimental results, the apparent deposition rate constant of silicic acid was estimated to be in the range from 10{sup -11} to 10{sup -10} m/s. The Damkoehler number (the dimensionless number expressed by the ratio of the characteristic time of mass transfer and the chemical reaction constant) was sufficient to cause the clogging of the host rock by deposition. This suggests that supersaturated silicic acid deposition, regardless of the presence of Mg ions (and Al ions), contributes to improving barrier performance against radionuclide migration in the host rock. (authors)

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

Distributed fiber sensing systems for 3D combustion temperature field monitoring in coal-fired boilers using optically generated acoustic waves (Final Report)

In this project, we have developed and tested three kinds of fiber optic sensing systems for real time monitoring of temperature variations within an industrial scale boiler furnace. The fiber optic sensing systems target spatial and temporal distributions of high temperature profiles in a boiler furnace in fossil power plants. The reconstructed temperature profile will provide critical input for the control mechanisms to optimize the combustion process. This temperature profile will address the essential problem for fossil power plants in achieving higher efficiency and fewer pollutant emissions. Acoustic pyrometer systems have been used to reconstruct temperature field of power plant boilers based on measuring TOF (times-of-flight) of sound waves along some straight paths in a 2D cross-section of the boiler. In this project, optically generated acoustic signals from a fiber optic sensing system have replaced the acoustic signals generated from an electrical transducer. A 3D reconstruction algorithm replaced the previous 2D model. In this project, three kinds of fiber optic sensing systems have been developed and tested. They are fiber optic sensing system I, fiber optic sensing system II (Distributed Sensing System I) and fiber optic sensing system III (Distributed Sensing System II). For fiber optic sensing system I, the fiber optic ultrasound generator acts as a signal generator. A microphone, hydrophone or other electronic devices serve as a signal receiver. In this system, there are one generator and one receiver. Distance test, water temperature test, air temperature test, air temperature reconstruction, and GE ISBF pilot test were performed by Fiber optic sensing system I. The fiber optic sensing system I successfully detected temperature in all these tests. We got 2D temperature reconstruction results by using the fiber optic sensing system I and it matched the reference data. The fiber optic sensing system I successfully survived in GE ISBF boiler environment (480 °F). For fiber optic sensing system II (Distributed Sensing System I), it is an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Fiber Bragg Grating (FBG) and Fabry-Perot (FP) sensor act as a signal receiver. In this system, there is one generator and one receiver. Aluminum plate temperature test, furnace high temperature test, and GE ISBF pilot test were performed by the fiber optic sensing system II. Fiber optic sensing system II successfully detected the temperature in all these tests. The fiber optic sensing system II successfully survived at up to 700 °C furnace environment and 320 °C GE ISBF boiler environment. For fiber optic sensing system III (Distributed Sensing System II), it is also an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Multiple FP fiber sensors act as signal receivers. In this system, there is one generator and three receivers. Three GE ISBF pilot tests were performed by fiber optic sensing system III. The fiber optic sensing system III survived in the cold flow tests in GE’s ISBF pilot test facility. However, we didn’t get high temperature data by using this system since the nanosecond laser issues. During the period of the project, test trials, data simulation and algorithm optimization was performed successfully. For real time temperature field construction, the sampling rate must be fast enough to capture the field variations. The technology of Code-division multiple access (CDMA) is well studied which could allow parallel multiplexing, even if signals overlap in time or frequencies. Moreover, it has been known that extending the length of signal significantly improves SNR. For acoustic signals, these multiplexing techniques have also been widely used, mainly for sonar and acoustic communications. The CDMA modulation technique has been proposed and studied to guarantee high network throughput, low channel access delay and low energy consumption. We have studied the temperature field reconstruction using Gaussian Radial Basis Functions (GRBF)-based approximation approach. Reconstruction of 3D temperature field using Neural Networks with measured TOF and known propagation paths is feasible. 2D and 3D temperature field reconstruction simulation results are achieved. The milestone status is shown in Table 1. We finished milestone 1-8 and milestone 10. For milestone 9, we did three pilot tests by using the fiber optic sensing system III (Distributed Sensing System II) at GE Power. However, due to the failure of the ns laser, we did not get the temperature results. We conducted some additional tasks that were not originally proposed: 1) We fabricated a fiber optic sensing system I and did a pilot test based on this system. 2) In the proposal, we proposed two pilot tests at GE Power. In reality, we finished at least seven pilot tests at GE Power. GE Power has made a lot of efforts for supporting the pilot tests. 3) We got a simulation results based on CDMA. In summary, most of the tasks have been accomplished. The outcome of this project removed a few barriers that hinder the achievement of the final product of the distributed sensing systems. With the successful accomplishment of this project, a prototype of the fiber optic sensing system can be fabricated to attract more interests from companies and other funding agencies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

In the Mix : A Workshop Merging Computational Chemistry and Electrochemistry Alongside Data Science

As chemistry expands to more complex and interdisciplinary areas, a new generation of diverse researchers must engage with science and learn effective cross-disciplinary collaboration and communication. To these ends, we designed and implemented In the Mix, a graduate student-led, two-day workshop for undergraduate students promoting collaborative science in the context of energy storage innovations. Here, the interactive workshop was designed for future and emerging researchers to gain hands-on experience with data science, computational chemistry, and electrochemistry techniques that are critical for developing materials for battery technologies. Participants also visited commercial renewable energy facilities to help them connect discovery-based research with industry and broader societal considerations. The workshop content and structure ensured that participants experienced the interrelatedness of the fields and understood the importance of collaborative research to yield scientific advances with real-world applications. An external team evaluated the workshop and participants’ perceptions of their experiences. While our research context was energy storage, the workshop goals and outcomes are applicable to other contexts. Interdisciplinary, experiential workshops are a key avenue to broadening participation in science and research, and the ideas presented here can be readily modified for other scientific contexts and/or incorporated as broader impact activities.

25 ENERGY STORAGE↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

In situ estimation of ice crystal properties at the South Pole using LED calibration data from the IceCube Neutrino Observatory

The IceCube Neutrino Observatory instruments about 1 km 3 of deep, glacial ice at the geographic South Pole. It uses 5160 photomultipliers to detect Cherenkov light emitted by charged relativistic particles. An unexpected light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow direction of the ice. We examine birefringent light propagation through the polycrystalline ice microstructure as a possible explanation for this effect. The predictions of a first-principles model developed for this purpose, in particular curved light trajectories resulting from asymmetric diffusion, provide a qualitatively good match to the main features of the data. This in turn allows us to deduce ice crystal properties. Since the wavelength of the detected light is short compared to the crystal size, these crystal properties include not only the crystal orientation fabric, but also the average crystal size and shape, as a function of depth. By adding small empirical corrections to this first-principles model, a quantitatively accurate description of the optical properties of the IceCube glacial ice is obtained. In this paper, we present the experimental signature of ice optical anisotropy observed in IceCube light-emitting diode (LED) calibration data, the theory and parameterization of the birefringence effect, the fitting procedures of these parameterizations to experimental data, and the inferred crystal properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗