SEARCH · Engineering Papers

Results for “density limit”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Disentangling core and edge mechanisms of the density limit in DIII-D negative triangularity plasmas

The density limit is investigated in the DIII-D negative triangularity plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans show distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered by the onset of a non-disruptive, high-field side radiation front and the associated cooling, which clamps the edge density below . In contrast, the core density continues to rise until it saturates, a state characterized by enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis ≈ 6), increased Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative cooling and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

density limit↗

DOE OSTI2025

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit↗

DOE OSTI2026

Real-time avoidance of the L-mode and H-mode density limit via machine-learned stability metrics

Reliable operation of burning plasma tokamaks will require robust control strategies to avoid macroscopic instability limits such as the L-mode and H-mode density limits (LDL, HDL). In this work, we explore closed-loop avoidance of these phenomena at DIII-D using machine-learned risk metrics. Feedback control is implemented via the ‘DL Supervisor’ scheme, which regulates the chosen risk metric by reducing the density target or increasing NBI heating in real-time. Using the LDL 25 risk metric, the LDL is reproducibly suppressed. We also introduce an HDL risk metric in this study, HDL 25 , which reduces the False Positive Rate by 2x compared to the Greenwald fraction. Applying this scaling to a plasma current ramp-down, we successfully avoid an HDL-driven H/L back-transition. These experiments constitute the first demonstration of real-time DL avoidance using machine-learned risk metrics. These instability metrics outline a path to safer high-density operation, more reliable ramp-down scenarios, and improved off-normal control for next-step devices such as ITER and SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗

DOE OSTI2026

Local Spin Density Approximation Strongly Improved by a Better-Informed Local Scaling of Its Self-Interaction Correction

The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.

Approximation↗

DOE OSTI2025

Turbulence suppression at extreme plasma densities on DIII-D and EAST

Recent high-poloidal-beta (high-βP) experiments on DIII-D and EAST have made coordinated breakthroughs for high confinement quality at high density near the Greenwald limit. Density gradient amplification of turbulence suppression at high βP can explain both of these achievements. Experiments on DIII-D have achieved Greenwald fraction (fGr = line-averaged density/Greenwald density) above 1 simultaneously with normalized energy confinement (H98y2) around 1.5, as required in fusion reactor designs but never before verified in tokamak experiments with the divertor configuration. A synergy between increased H98y2 and fGr is observed with strong gas puffing, due to the build-up of an internal transport barrier at large radius in the temperature and density channels. Transport simulations reveal that the favorable trend of reduced turbulent energy transport at higher density is only expected when increasing the density gradient at high local safety factor and high β, thus at high βP to ensure strong α-stabilization. These conditions are crucial to many conceptual designs for steady-state reactors. New experiments on EAST have nearly doubled the ion temperature at fGr ∼ 0.9, consistent with predict-first modeling results based on the same physics revealed from the DIII-D analysis. All previous EAST long-pulse H-modes have Ti ≪ Te near plasma axis. Transport modeling indicates that the profiles are limited by ion-temperature-gradient modes at mid-radius. The modeling also suggested potential solutions, including reducing magnetic shear, enhancing density gradients, and higher impurity concentration. Following this guidance, EAST experiments directly show a strong enhancement of Ti achieved with a combination of a second plasma current ramp-up, a density gradient increase, and a Zeff perturbation by a short pulse (100 ms) of impurity injection, as predicted by the earlier modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗

DOE OSTI2026

Binder-Free Graphite Anodes for Next-Generation High-Performance Lithium-Ion Batteries

High-energy density anodes are crucial for next-generation lithium-ion batteries (LIBs) particularly for electric vehicle (EV) applications. Sluggish lithium-diffusion kinetics coupled with conventional anode fabrication processes containing polymeric binders hinder fast-charging capabilities and high-energy density of graphite. Herein, we introduce a binder-free graphite anode fabrication strategy using the electrospinning technique that contains ~2.41% carbon nanotubes (CNTs). Our strategy relies on the formation of an interconnecting conductive CNT network coupled with an ultrathin N-doped carbon coating on graphite particles from sacrificial binders. This combination enhances both structural integrity and electrical conductivity and, in turn, improves fast-charging capabilities and high energy density of LIBs. The binder-free graphite anode achieves ~335.0 mAh g–1 capacity at C/3 rate over 400 cycles with capacity retention of >95% and average Coulombic efficiencies >99.95%. These promising results suggest that the binder-free anode fabrication with a multifunctional design approach could elevate the energy-density limits of the graphite anodes, solving high-energy density requirements of EVs, and potentially provides a path forward for the development of economically feasible energy storage systems for various applications.

Ozcan, Muca [ORNL] (ORCID:0000000320020474)↗

DOE OSTI2024

Characterizing the negative triangularity reactor core operating space with integrated modeling

Negative triangularity (NT) has received renewed interest as a fusion reactor regime due to its beneficial power-handling properties, including low scrape-off layer power and a larger divertor wetted area that facilitates simple divertor integration. NT experiments have also demonstrated core performance on par with positive triangularity (PT) high confinement mode (H-mode) without edge-localized modes (ELMs), encouraging further study of an NT reactor core. In this work, we use integrated modeling to scope the operating space around two NT reactor strategies. The first is the high-field, compact fusion pilot plant concept Modular, Adjustable, NT ARC (MANTA) and the second is a low field, high aspect ratio concept based on work by Medvedev et al. By integrating equilibrium, core transport, and edge ballooning instability models, we establish a range of operating points with less than 50 MW scrape-off layer power and fusion power comparable to PT H-mode reactor concepts. Heating and seeded impurities are leveraged to accomplish the same fusion performance and scrape-off layer exhaust power for various pressure edge boundary conditions. Scans over these pressure edge conditions accommodate any current uncertainty of the properties of the NT edge and show that the performance of an NT reactor will be extremely dependent on the edge pressure. The high-field case is found to enable lower scrape-off layer power because it is capable of reaching high fusion powers at a relatively compact size, which allows increased separatrix density without exceeding the Greenwald density limit. Adjustments in NT shaping exhibit small changes in fusion power, with an increase in fusion power density seen at weaker NT. Infinite-n ballooning instability models indicate that an NT reactor core can reach fusion powers comparable to leading PT H-mode reactor concepts while remaining ballooning-stable. Seeded krypton is leveraged to further lower scrape-off layer power since NT does not have a requirement to remain in H-mode while still maintaining high confinement. We contextualize the NT reactor operating space by comparing to popular PT H-mode reactor concepts, and find that NT exhibits competitive ELM-free performance with these concepts for a variety of edge conditions while maintaining relatively low scrape-off layer power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗

DOE OSTI2026

Controlled expansion for pairing in a polarized band with strong repulsion

Can strong repulsive interactions be shown to give rise to pairing in a controlled way? We find that for a single flavor polarized band, there is a small expansion parameter in the low density limit, once the Bloch wavefunction form factor is taken into account. A perturbative expansion is possible, even if the interaction is much stronger than the Fermi energy ε F . As a matter of principle, our work shows analytically how strong pairing can emerge from strong repulsion. Here, we illustrate our method with two examples: a 2D Dirac model and a 1D tight binding model with two orbitals. In the latter case, using density matrix renormalization group, we show that the analytical theory indeed guided us to discover the parameter regime where p-wave pairing with order-1 strength is dominant.

Ferromagnetic superconductors↗

DOE OSTI2025

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗

DOE OSTI2025

Next-generation tunnel FETs: exploring material perspectives and areal tunneling configurations

The end of Dennard scaling, which facilitated proportional increases in computing power without added energy costs until the mid-2000s, has underscored the urgent need for innovative semiconductor devices that can enhance energy efficiency. Tunnel field-effect transistors (TFETs) have emerged as promising candidates to surpass the energy efficiency of conventional metal oxide semiconductor field-effect transistors (MOSFETs). Unlike MOSFETs, which rely on thermionic emission to overcome the source-channel potential barrier, TFETs operate through quantum tunneling, potentially enabling sub-60 mV dec −1 subthreshold swing (SS) for low-voltage operation. However, lateral TFETs have faced challenges in achieving adequate on-state current (I ON ) and a broad SS operation window, limiting their practical utility. This review article advocates for areal TFETs, which utilize face-to-face tunnel junctions that ideally offer step-function current turn-on characteristics and allow I ON to scale with device area rather than width. We highlight recent advancements in integrating 2D materials into tunneling structures, which could facilitate efficient band-to-band tunneling through atomically thin layers, while addressing challenges of gate field screening. We then discuss the nearer-term prospects of epitaxial areal TFETs comprising III–V compound semiconductors and group-IV semiconductors based on recent experimental progress. The review examines both quantum mechanical and semiclassical modeling approaches for TFETs, including techniques to reduce the computational complexity. The article delves into ongoing challenges in material synthesis, interface engineering, device fabrication, and integration pathways, concluding with recommendations for future research directions to overcome the fundamental power density limitations of conventional transistor technology.

2D materials↗

NASA NTRS2026

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry↗

DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗

DOE OSTI2025

Elucidating Electric Field-Induced Rate Promotion of Brønsted Acid-Catalyzed Alcohol Dehydration

Applied potentials have been demonstrated as a powerful tool to promote heterogeneous Brønsted acid catalysis by orders of magnitude, leveraging interfacial electric fields to stabilize protonated intermediates. However, the use of flat two-dimensional electrodes with inherently low active site densities limits the application of conventional thermochemical characterization techniques that can probe the nature of catalytic active sites. Here, we use kinetic analyses with an electrostatics-based model to elucidate the intricacies of potential-induced rate promotion, employing liquid-phase dehydration of 1-methylcyclopentanol catalyzed by carboxylic acid groups on carbon nanotubes as a probe system. By using a basket electrode to directly polarize catalyst powder, we demonstrate that thermocatalytic reaction rates can be promoted by 100,000-fold, exhibiting a log–linear dependence on applied potential with rate-potential scalings as high as 125 ± 4 mV per 10-fold rate increase. In agreement with model predictions, we show that lower ionic strengths attenuate potential sensitivity, resulting from a weakening of the interfacial electric field that interacts with the acidic proton. Furthermore, we experimentally confirm the model-predicted “isokinetic potential” (at ∼0.6 V vs Ag/AgCl)─the potential at which all rate scaling lines at various ionic strengths intersect, making the rate independent of ionic strength. Base titrations reveal that only ∼8% of the carboxylic acid sites are catalytically active, yet these same active sites are operational at the highest and lowest potentials. Collectively, our results provide a key methodology for modeling catalytic effects of electric fields, quantifying active sites under applied potential, and demonstrating fundamental principles of electric field-induced rate promotion.

Catalysts↗