Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric instabilities”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Interpreting and Stabilizing Machine-Learning Parametrizations of Convection

Neural networks are a promising technique for parameterizing subgrid-scale physics (e.g., moist atmospheric convection) in coarse-resolution climate models, but their lack of interpretability and reliability prevents widespread adoption. For instance, it is not fully understood why neural network parameterizations often cause dramatic instability when coupled to atmospheric fluid dynamics. This paper introduces tools for interpreting their behavior that are customized to the parameterization task. First, we assess the nonlinear sensitivity of a neural network to lower-tropospheric stability and the midtropospheric moisture, two widely studied controls of moist convection. Second, we couple the linearized response functions of these neural networks to simplified gravity wave dynamics, and analytically diagnose the corresponding phase speeds, growth rates, wavelengths, and spatial structures. To demonstrate their versatility, these techniques are tested on two sets of neural networks, one trained with a super-parameterized version of the Community Atmosphere Model (SPCAM) and the second with a near-global cloud-resolving model (GCRM). Additionally, even though the SPCAM simulation has a warmer climate than the cloud-resolving model, both neural networks predict stronger heating/drying in moist and unstable environments, which is consistent with observations. Moreover, the spectral analysis can predict that instability occurs when GCMs are coupled to networks that support gravity waves that are unstable and have phase speeds larger than 5 m s -1 . In contrast, standing unstable modes do not cause catastrophic instability. Using these tools, differences between the SPCAM-trained versus GCRM-trained neural networks are analyzed, and strategies to incrementally improve both of their coupled online performance unveiled.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spectral and temporal shaping of spectrally incoherent pulses in the infrared and ultraviolet

Laser–plasma instabilities (LPIs) hinder the interaction of high-energy laser pulses with targets. Simulations show that broadband spectrally incoherent pulses can mitigate these instabilities. Optimizing laser operation and target interaction requires controlling the properties of these optical pulses. We demonstrate closed-loop control of the spectral density and pulse shape of nanosecond spectrally incoherent pulses after optical parametric amplification in the infrared (∼1053 nm) and sum–frequency generation to the ultraviolet (∼351 nm) using spectral and temporal modulation in the fiber front end. The high versatility of the demonstrated approaches can support the generation of high-energy, spectrally incoherent pulses by future laser facilities for improved LPI mitigation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2020

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planetary waves can activate resonant drag instabilities in 3D dusty gaseous discs

Resonant drag instabilities (RDIs) in protoplanetary discs are driven by the aerodynamic back-reaction of dust on gas and occur when the relative dust-gas motion resonate with a wave mode intrinsic to the gas fluid. Axisymmetric models indicate that the RDI generates filamentary perturbations, leading to grain clumping and planetesimal formation. Motivated by these findings, we investigate the dust–gas interaction in a non-axisymmetric inviscid protoplanetary disc with an embedded low-mass planet ($M_{\mathrm{p}}\in [0.3, 3] \, \rm M_{\oplus }$, here $\, \rm M_{\oplus }$ is the Earth mass). We conduct global 3D high-resolution two-fluid simulations, with the dust being parametrized by the Stokes number $\mathrm{St}\in [0.01,0.5]$. We find that planetary waves (PWs; also known as Rossby waves), which propagate along the downstream separatrices of the horseshoe region, resonate with the streaming motion and trigger the RDI. The consequent development of a global-scale filamentary dust distribution does not sensitively depend on the Stokes number, nor does it depend on the fast dust settling that takes place in an inviscid disc. The rapid onset of this instability, which is comparable to the dynamical orbital time-scale, suppresses the formation of asymmetric structures in the dust in the vicinity of the planet (such as dust voids and filaments). Additionally, we find that the dust feedback enables buoyancy resonances in an otherwise non-buoyant (globally isothermal) disc. Therefore, our results provide the first numerical evidence of RDIs generation driven by planetary waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tunable Laser-Plasma Amplifier (Final Report)

Exploration at the laser intensity frontier has always offered new avenues for physics and reaching beyond this frontier is a grand challenge. Present-day petawatt-class lasers provide focused intensities on target of 10 22 W/cm 2 , corresponding to electric fields of 200 TV/m, while laser-plasma amplification opens a route for focused intensities well above 10 23 W/cm 2 . Intensities in this range provide the ability to test quantum electrodynamics in the unexplored low-energy, strong-field regime where signatures for new physics may arise. The behavior of matter under such extraordinary conditions is a rich and fascinating subject, not only in its own right in fundamental plasma physics, but for the many potential applications that promise to enrich the natural sciences in the future, including compact electron-beam, ion-beam particle accelerators, and ultra-bright X-ray sources. Although laser systems are now under construction internationally to access intensities of 10 23 W/cm 2 , the current technologies used appear to be fundamentally limited to these intensities. The realization of intensities beyond 10 23 W/cm 2 using parametric amplification in plasmas promises a breakthrough in high-energy density physics. Parametric amplification using Raman scattering in a plasma could provide the enabling technology for the generation of ultra-high-power laser pulses, but a more complete understanding of the nonlinear optics of plasmas is required. There is a significant gap in well-diagnosed laser-plasma instability studies of nonlinear plasma-wave phenomena, which are critical to understand for future laser-plasma devices. To achieve an efficient laser-plasma amplifier, plasma waves must be driven to large amplitude where significant energy can be rapidly transferred from the pump to the seed over the pulse duration of the seed. Simulations suggest that this nonlinear pump depletion regime can be achieved in the “pi-pulse” amplification regime. While simulations show this optimal regime with efficient amplification, it has remained elusive in experiments and there is a growing consensus within the community that thermal effects and pump beam limitations prevent laser-plasma amplifiers from progressing through the linear regime into the nonlinear pump depletion regime. Previous experiments have been significantly limited by the laser power available at the necessary wavelengths for the seed laser; therefore, the amplification is required to start in the linear regime where it is sensitive to many deleterious effects. The enabling technology (currently unique to plasma-wave amplification in the world) at the University of Rochester is the ability to provide a seed pulse with sufficient power (4 mJ/100 fs seed) to immediately drive nonlinear plasma waves into the pi-pulse regime and to tune its wavelength to optimize the efficiency of energy transfer. This in combination with the state-of-the-art OMEGA heater beams providing multiple kilojoules in a nanosecond to sufficiently heat the plasma make this system distinct from previous studies. These heater beams will provide, for the first time in Raman amplification studies, a homogeneous electron temperature high enough to prevent pump beam propagation issues that have plagued previous experiments. These systems will provide a platform for driving electron-plasma waves into the nonlinear regime where pump depletion and pulse shortening are predicted to lead to high amplification efficiencies (>30%). The Team has made significant progress through prior support from DOE Fusion Energy Sciences [DOE Office of Science Award Number DE-SC0016253 (2016-2022)]. This includes twenty-two peer-reviewed manuscripts, one patent, ten contributed talks presented at international conferences, and research that was highlighted as invited talks at fifteen international conferences. The broader impacts of this research are evident in the support of early career scientists, two Ph.D. theses, four current graduate students, a Masters Project, two undergraduate researchers, and an underrepresented minority student hired through the California Alliance for Minority Participation who now works as a Research Engineer in the group. This research met all of the funded research objectives and the highlights from primary Raman amplification thrust of this work are discussed below and form the foundation for the proposed research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Grey Galaxies’ as an endpoint of the Kerr-AdS superradiant instability

Kerr-AdS d+1 black holes for d ≥ 3 suffer from classical superradiant instabilities over a range of masses above extremality. We conjecture that these instabilities settle down into Grey Galaxies (GGs) — a new class of coarse-grained solutions to Einstein’s equations which we construct in d = 3. Grey Galaxies are made up of a black hole with critical angular velocity ω = 1 in the ‘centre’ of AdS, surrounded by a large flat disk of thermal bulk gas that revolves around the centre of AdS at the speed of light. The gas carries a finite fraction of the total energy, as its parametrically low energy density and large radius are inversely related. GGs exist at masses that extend all the way down to the unitarity bound. Their thermodynamics is that of a weakly interacting mix of Kerr-AdS black holes and the bulk gas. Their boundary stress tensor is the sum of a smooth ‘black hole’ contribution and a peaked gas contribution that is delta function localized around the equator of the boundary sphere in the large N limit. We also construct another class of solutions with the same charges; ‘Revolving Black Holes (RBHs)’. RBHs are macroscopically charged SO(d, 2) descendants of AdS-Kerr solutions, and consist of ω = 1 black holes revolving around the centre of AdS at a fixed radial location but in a quantum wave function in the angular directions. RBH solutions are marginally entropically subdominant to GG solutions and do not constitute the endpoint of the superradiant instability. Nonetheless, we argue that supersymmetric versions of these solutions have interesting implications for the spectrum of supersymmetric states in, e.g. N = 4 Yang-Mills theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗
DOE OSTI2021

Long-Term Creep-Rupture Behavior of Alloy Inconel 740/740H

To explore potential application of Ni-based alloys for power generation at the higher temperatures and pressures needed to achieve high thermal to electrical-energy conversion efficiency, an extensive creep-rupture dataset covering up to 875 °C and almost 70,000 h for Inconel 740/740H was analyzed using Larson–Miller parameter and Wilshire approaches. The results were used to assess the relative effectiveness of the two analytical methods, both in describing the experimental data and, because of the breadth of the dataset, using analyses of its shorter-time data to make creep lifetime predictions for much more extended times, which were then directly compared to the measured rupture times. The respective methods were also used to predict creep-limited lifetimes relevant to power production (that is, 100,000 h or greater). Despite the complexity of the precipitation-strengthened Inconel 740/740H alloy and the generalized parametric approach of these methodologies, the predictions based on such were reasonably accurate when the entire dataset was analyzed. However, when the analysis was confined to only data for conditions yielding creep-rupture times < 5,000 h( about 65 % of the entire dataset), the Wilshire correlation yielded better prediction for longer time lifetimes due to the inherent instability of the specific Larson–Miller formalism used in this analysis when extrapolated significantly outside its analysis range.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗
DOE OSTI2022

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Implementation and Exploration of Parameterizations of Large-Scale Dynamics in NCAR's Single Column Atmosphere Model SCAM6

A single column model with parameterized large-scale (LS) dynamics is used to better understand the response of steady-state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large-scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR's Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies -Weak Temperature Gradient (GASS-WTG) Intercomparison project. Using the same parameter settings and simulation configuration as in the GASS-WTG Intercomparison project, SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the large-scale circulation's damping of tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS-WTG Intercomparison project. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid-scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column's reference radiative-convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Control Parameter Sensitivity Study for Inverter-Based-Resource Dominated Grids: A Small Signal Stability Approach and Framework

The growing adoption of renewable energy is driving the prevalence of inverter-based resources (IBRs) within power grids. Future power grids will integrate both grid-following IBRs (GFM-IBRs) and grid-forming IBRs (GFL-IBRs) alongside synchronous generators. Therefore, it is crucial to perform stability studies that account for all components and especially control interactions related to IBRs. Extensive research has performed to study the IBR-related stability, however, the sensitivity study of IBRs' control parameters on system stability has not been adequately studied yet, especially from a systematic way. Therefore, this paper conducts a small signal stability analysis for a generic grid with multiple types of resources and develops an analytical framework for assessing the sensitivity of control parameters affecting stability margins. To achieve that, the non-autonomous reduced-order non-linear dynamic model is developed for a generic power system with multiple synchronous generator-based resources (SGBRs), GFM-IBRs, and GFL-IBRs. Based on the analytic model, a systematic framework for parametric sensitivity on systems' asymptotic stability is developed. A parameter sensitivity analysis based on eigenvalue methods is proposed. The impact of the droop controllers of GFM-IBRs, PQ-dispatch and the PLL controller of GFL-IBR on the system asymptotic stability is discussed. This sensitivity study is aiming to provide deep insights on control parameters' impact on system stability, and gives direction for parameter tuning in case of instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗
DOE OSTI2025

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗