Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactor Kinetics”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2025

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2021

A percolation theory for designing corrosion-resistant alloys

Iron-chromium and nickel-chromium binary alloys containing sufficient quantities of chromium serve as the prototypical corrosion-resistant metals owing to the presence of a nanometre-thick protective passive oxide film. Should this film be compromised by a scratch or abrasive wear, it reforms with little accompanying metal dissolution, a key criterion for good passive behaviour. This is a principal reason that stainless steels and other chromium-containing alloys are used in critical applications ranging from biomedical implants to nuclear reactor components. Unravelling the compositional dependence of this electrochemical behaviour is a long-standing unanswered question in corrosion science. Herein, we develop a percolation theory of alloy passivation based on two-dimensional to three-dimensional crossover effects that accounts for selective dissolution and the quantity of metal dissolved during the initial stage of passive film formation. We validate this theory both experimentally and by kinetic Monte Carlo simulation. Our results reveal a path forward for the design of corrosion-resistant metallic alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of iodide ions on the speciation of radiolytic transients in molten LiCl–KCl eutectic salt mixtures

The fate of fission-product iodine is critical for the deployment of next generation molten salt reactor technologies, owing to its volatility and biological impacts if it were to be released into the environment. To date, little is known on how ionizing radiation fields influence the redox chemistry, speciation, and transport of iodine in high temperature molten salts. Here we employ picosecond electron pulse irradiation techniques to elucidate for the first time the impact of iodide ions (I – ) on the speciation and chemical kinetics of the primary radiation-induced transient radicals generated in molten chloride salt mixtures (e S – and Cl 2 ˙– ) as a function of temperature (400–700 °C). In the presence of I – ions (≥ 1 wt% KI in LiCl–KCl eutectic), we find that the transient spectrum following the electron pulse is composed of at least three overlapping species: the e S – and the Cl 2 ˙– and ICl ˙– radical anions, for which a deconvoluted spectrum of the latter is reported here for the first time in molten salts. This new transient spectrum was consistent with gas phase density functional theory calculations. The lifetime of the e S – was unaffected by the addition of I – ions. The newly observed interhalogen radical anion, ICl˙ – , exhibited a lifetime on the order of microseconds over the investigated temperature range. The associated chemical kinetics indicate that the predominate mechanism of ICl˙ – decay is via reaction with the Cl 2 ˙– radical anion. The iodine containing product of this reaction is expected to be ICl 2 – , which will have implications for the transport of fission-product iodine in MSR technologies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Elevated temperature nanoscratch of Inconel 617 Superalloy

Inconel 617 superalloy is a main candidate to be used for mechanical and tribo-components in high temperature helium-cooled reactors. Recent findings show that it grows a unique surface oxide, especially under high temperature helium with distinct wear, friction, and contact properties. Here, this study reports the elevated temperature nanoscratch behavior of Inconel 617 and further utilizes it to understand the effect of temperature on contact friction constituent contributors, adhesion and plowing at small scales. Inconel 617 is aged in high temperature helium, and consequently, the total kinetic friction coefficient of the alloy surface oxide is obtained in temperatures ranging from 25 °C to 400 °C. A finite element model is developed and validated based on the experimental results. The model is then utilized along with previously established techniques to determine the adhesion and plowing components of the friction coefficient. At small scale, the experimental results show that with increasing temperature the friction coefficient increases. It was inferred that this increase is mainly due to the increased contribution of plowing friction at high levels of deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
NASA NTRS2008

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗
DOE OSTI2021

Isotope effects on transport in LHD

Abstract Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time ( τ E ) with operational parameters shows positive mass dependence ( M 0.27 ; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence ( M 0.0 ) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Physics↗
DOE OSTI2020

Multidimensional simulations of Mckenna-driven flow tube configuration: Investigating non-ideality in NO x formation flow tube experiments

Multidimensional simulations have been conducted to simulate atmospheric pressure, flat-flame/McKenna-burner-driven-flow tube experiments targeted to obtain NO x speciation data for predicting/analyzing syngas combustion emissions. In a prior work, we demonstrated the impacts of multidimensional transport on post flame region prediction departures from those assuming unidimensional flow/transport conditions. In this work, we develop and utilize a multidimensional laminar reacting flow solver to simulate the fully coupled flame and post flame regions to further elucidate the impacts of the earlier unidimensional modeling assumptions on interpreting post flame NO x experimental data. The model is used to simulate a lean, premixed syngas/air flame and its associated post flame regions within a cylindrical flow-tube-like arrangement. The combustion process takes place under atmospheric condition with trace amount of NO x seeding fed into the inlet gas stream. The spatial evolution of NO x species (NO and NO 2 ) in the flame and in the post-combustion zone suggests two distinct regions: 1) a region encompassing the flame structure itself; and 2) a post flame region in which the temperature decays due to both axial and radial transport processes. The predictions show that for the conditions studied, a pulsatile flow field exists due to the formation of an expanding and contracting recirculation zone in the outer periphery of the flow tube. By resolving the nature of the flow, the resulting time-averaged temperature and species concentrations show improved agreement with existing experimental measurements. The flow-field interaction results in radial inhomogeneities in the NO 2 profiles with the maximum concentration offset from the flow centerline. The location of the peak in NO 2 is coupled with radial temperature gradients from wall cooling effects and their significant influence on NO/NO 2 interconversion kinetics, producing notable NO 2 accumulation in regions near the wall. Geometrical configurations capable of suppressing/minimizing the pulsatile nature are also investigated and the results are compared. Other experimental configurations could be considered in parametric simulations to determine the optimal configuration that would minimize non-idealities in the observations. The work shows the value in performing such computations in advance of settling on a particular design for flow tube/flow reactor experiments.

42 ENGINEERING↗
DOE OSTI2021

In Situ Grain Growth of Nanograined Magnetite under Ion Irradiation at Room Temperature and 500 ℃

Grain growth studies of nanograined metals under both thermal and ion irradiation conditions has been studied extensively. Oxides, on the other hand, have not been studied as extensively, especially under irradiation. Some studies have focused on nanograined oxide grain growth on ZnO under thermal conditions only and on SnO 2 and CeO 2 using 250 keV Ar 2+ and 3 MeV Au + ions at 300 K, respectively. There does not appear to be any studies on grain growth of nanosized oxides which are commonly found to form under light water reactor conditions. Oxides such as magnetite (Fe 3 O 4 ), hematite (α-Fe 2 O 3 ), maghemite (γ-Fe 2 O 3 ), and wustite (FeO) are examples of common oxides which have not been studied. Here, this study focused on the grain growth of nanograined Fe 3 O 4 under both thermal and ion irradiation conditions in situ inside a Transmission Electron Microscope (TEM) to better understand the kinetics of grain growth of oxides.

36 MATERIALS SCIENCE↗
NASA NTRS1984

GCKP84--General Chemical Kinetics Code for Gas-Phase Flow and Batch Processes Including Heat Transfer Effects

A general chemical kinetics code is described for complex, homogeneous ideal gas reactions in any chemical system. The main features of the GCKP84 code are flexibility, convenience, and speed of computation for many different reaction conditions. The code, which replaces the GCKP code published previously, solves numerically the differential equations for complex reaction in a batch system or one dimensional inviscid flow. It also solves numerically the nonlinear algebraic equations describing the well stirred reactor. A new state of the art numerical integration method is used for greatly increased speed in handling systems of stiff differential equations. The theory and the computer program, including details of input preparation and a guide to using the code are given.

David A Bittker↗
DOE OSTI2024

Computational modeling of microalgal biofilm growth in heterogeneous rotating algal biofilm reactors (RABRs) for wastewater treatment

Rotating algal biofilm reactors (RABRs) are innovative systems designed to cultivate microalgae biofilms efficiently. In this paper, we have developed a novel mathematical model to accurately capture the growth dynamics of algae biofilms within RABR. By considering the spatial heterogeneity of the RABR, we introduce a PDE-based model that addresses the spatial variations across the substratum, enabling a more accurate simulation of biofilm growth in RABRs. The photosynthesis process is modeled through reactive kinetics, driving the growth of the algae biofilm. To analyze the system's behavior, we employ finite difference numerical methods to solve the complex PDE model. We then conduct extensive numerical simulations to understand algae biofilm growth in the RABR environment under various operational factors and environmental conditions. One primary focus in these simulations is to investigate the impact of various harvesting strategies, harvesting frequencies, light intensity, and light exposure on the overall biomass productivity of the algae biofilm. The numerical results provide valuable insights into optimizing algae biofilm growth and designing harvesting techniques in RABR systems. Our proposed novel mathematical model provides an effective platform for the theoretical investigation and design of RABRs for wastewater treatment.

09 BIOMASS FUELS↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

NSTX-U liquid metal core-edge facility (LMCE)

NSTX-U/LMCE will provide a unique and world-leading research facility to address the primary challenge to delivering economic and timely magnetic fusion energy, namely the need to develop a power and particle exhaust and first-wall system that can withstand very high edge heat fluxes, maximize energy confinement, and avoid the production of large masses of solid eroded first-wall material. The NSTX-U/LMCE facility will assess the ability of liquid metals (LMs) – especially liquid lithium – to provide a new boundary condition for magnetic fusion systems, to extend the lifetime of the plasma facing components (PFCs) and improve core plasma confinement. Such capability is needed to establish the basis for next-step fusion facilities including fusion pilot plants, and to maintain U.S. world leadership in core-edge integration research. NSTX-U/LMCE will leverage the ability to generate very high divertor perpendicular heat flux q⊥ ~ 100MW/m 2 , extensive diagnostics, and liquid-metal-applicable infrastructure of NSTX-U. NSTX-U/LMCE will provide access to a high-confinement plasma core with majority self-driven plasma current, the flexibility to test a range of liquid metal divertor concepts, access to a range of separatrix collisionalities (from high to very low), and the ability to controllably vary the first-wall temperature to vary the plasma- wall interaction physics on liquid lithium components. Further, NSTX-U/LMCE will utilize more reactor-relevant high-Z refractory-metal PFC substrates. With these capabilities the NSTX-U/LMCE facility will explore the full continuum of core-edge solutions ranging from high core radiated power, to conditions with radiative losses concentrated in the scrape-off layer (SOL), and ultimately low recycling conditions. The low collisionality SOL that may be accessible in the low recycling regime is relatively unexplored and will require a kinetic treatment of the edge, which can be addressed theoretically, and with experiments in LTX-β. Additional smaller-scale preparatory R&D facilities will be required to reduce the risk of premature technical/engineering failure of liquid metal systems implemented in NSTX-U. The NSTX-U/LMCE facility aligns very well with recommendations in the FESAC Long-Range Plan and NASEM Pilot Plant reports and the Bold Decadal Vision, will be unique in the world program throughout the next decade, and is garnering private company interest in utilizing NSTX-U/LMCE for development of LM PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Combined Nitrogen and Phosphorous Recovery via Electrochemical Technology Integration into Municipal Wastewater Treatment Plants

This project developed electroN-P, an electrochemical technology for recovering nitrogen and phosphorus as a fertilizer product from anaerobic digester centrate, targeting an advancement from TRL 4 to 6. The reactor was scaled 220× from a 250 mL batch cell to a continuous 4-channel system treating 55 L of wastewater, recovering >80% of phosphorus. Energy consumption under constant-voltage operation was lower than the embedded energy of conventional fertilizers, while constant-current operation produced a cost-competitive product with a smaller reactor footprint. TEA and LCA outcomes were highly sensitive to the magnesium source; alternative magnesium salt configurations projected reduced costs and embedded energy compared to sacrificial rods, but at the expense of significantly longer operation times. The recovered fertilizer performed comparably to diammonium phosphate and triple superphosphate in soil and plant trials. Integrating the technology into a whole-plant model reduced aeration energy by nearly 50%, lowering the levelized cost of water treatment from $0.668/m 3 to $0.648/m 3 . Results support targeting commercialization at smaller (≤1 MGD) facilities with high-strength digester streams. Further investigation into the transport and corrosion kinetics governing sacrificial anode wear is recommended to support continuous, longer-term operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Methanol oxidation up to 100 atm in a supercritical pressure jet-stirred reactor

Methanol (CH 3 OH) has attracted considerable attention as a renewable fuel or fuel additive with low greenhouse gas emissions. Methanol oxidation was studied using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at pressures of 10 and 100 atm, at temperatures from 550 to 950 K, and at equivalence ratios of 0.1, 1.0, and 9.0 in experiments and simulations. Here, the experimental results show that the onset temperature of CH 3 OH oxidation at 100 atm is around 700 K, which is more than 100 K lower than the onset at 10 atm and this trend cannot be predicted by the existing kinetics models. Furthermore, a negative temperature coefficient (NTC) behavior was clearly observed at 100 atm at fuel rich conditions for methanol for the first time. To understand the observed temperature shift in the reactivity and the NTC effect, we updated some key elementary reaction rates of relevance to high pressure CH 3 OH oxidation from the literature and added some new low-temperature reaction pathways such as CH 2 O + HO 2 = HOCH 2 O 2 (RO 2 ), RO 2 + RO 2 = HOCH 2 O (RO) + HOCH 2 O (RO) + O 2 , and CH 3 OH + RO 2 = CH 2 OH + HOCH 2 O 2 H (ROOH). Although the model with these updates improves the prediction somewhat for the experimental data at 100 atm and reproduces well high-temperature ignition delay times and laminar flame speed data in the literature, discrepancies still exist for some aspects of the 100 atm low-temperature oxidation data. In addition, it was found that the pressure-dependent HO 2 chemistry shifts to lower temperature as the pressure increases such that the NTC effect at fuel-lean conditions is suppressed. Therefore, as shown in the experiments, the NTC phenomenon was only observed at the fuel-rich condition where fuel radicals are abundant and the HO 2 chemistry at high pressure is weakened by the lack of oxygen resulting in comparatively little HO 2 formation.

33 ADVANCED PROPULSION SYSTEMS↗