SEARCH · Search NASA
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.
Presentation to introduce a number of university-led NEUP projects on the subject of thermal-hydraulic behavior of accident tolerant fuel materials for the 2020 AFC Annual Review Meeting.
Abstract We present the results of GENE gyrokinetic calculations based on a series of JET–ITER-like-wall (ILW) type I ELMy H-mode discharges operating with similar experimental inputs but at different levels of power and gas fuelling. We show that turbulence due to electron-temperature-gradient (ETGs) modes produces a significant amount of heat flux in four JET–ILW discharges, and, when combined with neoclassical simulations, is able to reproduce the experimental heat flux for the two low gas pulses. The simulations plausibly reproduce the high-gas heat fluxes as well, although power balance analysis is complicated by short ELM cycles. By independently varying the normalised temperature gradients ( ω T e ) and normalised density gradients ( ω n e ) around their experimental values, we demonstrate that it is the ratio of these two quantities η e = ω T e / ω n e that determines the location of the peak in the ETG growth rate and heat flux spectra. The heat flux increases rapidly as η e increases above the experimental point, suggesting that ETGs limit the temperature gradient in these pulses. When quantities are normalised using the minor radius, only increases in ω T e produce appreciable increases in the ETG growth rates, as well as the largest increases in turbulent heat flux which follow scalings similar to that of critical balance theory. However, when the heat flux is normalised to the electron gyro-Bohm heat flux using the temperature gradient scale length L T e , it follows a linear trend in correspondence with previous work by different authors.
Accurate monitoring of boiling heat transfer is critical for safeguarding high-power systems operating in environments where conventional optical diagnostics are hindered by radiation fields or restricted visual accessibility. This study presents a non-intrusive framework that integrates hydroacoustic sensing with deep learning to infer near-wall boiling characteristics and enable predictive thermal assessment without visual access. In a prototypical subcooled flow-boiling facility representative of the Isotope Production Facility (IPF) at Los Alamos, hydrophones capture boiling-induced acoustic emissions that are transformed into background-removed Short-Time Fourier Transform (STFT) spectrograms. A convolutional neural network (CNN) then regresses heat flux, wall superheat, and key bubble parameters directly from these spectrograms. The CNN achieved predictive accuracy under nominal conditions and demonstrated robustness and generalization under acoustic noise for Signal-to-Noise Ratios (SNRs) down to approximately 0 dB. When integrated into an ANSYS CFX wall-boiling model, the acoustically inferred parameters reproduced boiling curve and critical heat flux (CHF) values consistent with image-based benchmarks. Furthermore, the model retained reliable performance under moderate variations in bulk temperature, flow rate, and hydrophone placement, confirming its generalizability across practical boundary conditions. These results demonstrate the feasibility of hydroacoustic-based deep learning as a viable path toward real-time, radiation-tolerant boiling diagnostics and predictive thermal safety assessment in inaccessible systems such as the IPF.
Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,\textrm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16 $\,\pm\,7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e = 5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream n e and T e profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,\textrm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.
Condensation tests were performed using a newly developed test facility for scaling the passive containment cooling system (PCCS) to a small modular reactor (SMR). The PCCS of the SMR plays a pivotal role in ensuring greater safety, reliability, and compactness than what is afforded by traditional reactors. Therefore, a well-designed PCCS is essential to SMRs. However, previous studies and test data were unsuitable for scaling, due to high variation in the test geometry and operating conditions. This study intends to close this research gap by using a novel designed scaled test facility consisting of vertical condensing test sections featuring 1-, 2-, and 4-inch-diameter condensing tubes with annular water cooling, and by applying superheated and saturated steam with different steam mass flow ranges of 5–25 g/s. Further, the primary test data, including axial temperatures, mass flow rates, and pressures, were used in conjunction with a standard data reduction method to estimate critical parameters such as heat fluxes, heat transfer coefficients, and condensation rates. These scaled test data would support improving empirical correlations and validating condensation models to identify scaling distortion for SMR PCCSs.
The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. Here, the heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.
The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.
Dissipating large heat fluxes from a surface is critically important in numerous industrial and natural applications. Boiling based spray cooling and surface texturing are two of the most promising methods being investigated to address this problem. Although our understanding on these topics has significantly improved over past decades, critical gaps remain in the knowledgebase stymieing the realization of their full potential. As an example, while bubble growth in pool boiling have been investigated in detail, comparatively little is known about how the bubbles evolve inside boiling drops. In the present work, we have investigated for the first time, the microbubble dynamics inside water droplets boiling on superhydrophilic textured substrates using high-speed X-ray phase contrast imaging (XRPCI). Our observations show that the transient bubble density variation follows similar characteristics irrespective of the texture spacing at a given surface temperature. For an example microstructure, we found that the number of discrete bubbles on the surface decreases as temperature is increased although their growth rate increases. Here, we observe that bubble growth is highly non-uniform during the lifetime of a drop on the surface. Initially, bubbles grow under diffusion-limited regime, but at later times they grow as ~t 1.45 due to combined effects of coalescence and evaporation. In some conditions, we found that bubbles shrink dramatically af- ter the initial growth spurt presumably due to severe quenching of the surface, and migration of bubbles on the surface. Using the bubble sizes, for the first time we analyzed the heat flux removed by a single bubble and also by all the bubbles at a given time. We find that the highest dissipation through latent- heat component (~600 W/cm 2 ) occurs just in the beginning and thereafter it decreases. We expect that our findings and the analysis would guide further work on the topic and will aid in the overarching goal of engineering surfaces that are more efficient in boiling heat transfer.
Heat flux sensors compatible with hot environments are critical to advance aerospace, materials, and energy generation technologies that cope with extreme thermal conditions. In this work, we report on the development and characterization of a high-temperature heat flux sensor using the transverse Seebeck effect in rhenium single crystals. The sensor leverages refractory alloys and ceramics compatible with temperatures exceeding 1000 °C. The heat flux sensor was characterized from room temperature to 500 °C using a temperature-controlled calibration facility. At constant temperature, the sensor’s voltage output is linear with respect to the absorbed heat flux. The responsivity of the sensor varies with temperature, from 1.3 μV/(W/cm 2 ) at room temperature to −3.2 μV/(W/cm 2 ) at 500 °C, increasing monotonically in magnitude after changing sign from positive to negative at approximately 300 °C. The experimental results are in good agreement with analytical predictions of the sensor’s temperature-dependent responsivity, which suggest a further increase in magnitude up to −7.4 μV/(W/cm 2 ) at 1000 °C. These results highlight the unique characteristics of rhenium as a TSE transducer. The design offers compatibility with a wide range of operating temperatures and yields a measurement sensitivity that increases as the environmental conditions become more challenging.
In future tokamak reactors like ITER and the Fusion Pilot Plant (FPP), real-time feedback control of heat flux to the plasma-facing components (PFC) will be critical for steady-state operation. This work presents the first experimental demonstration of real-time divertor heat flux estimation with infrared thermography and feedback control with impurity seeding on the DIII-D tokamak. The flexible infrastructure of the Plasma Control System (PCS) on DIII-D makes this new capability possible. The PCS software runs on a gateway computer system, and five real-time compute nodes. An array of low latency streaming digitizers from D-TACQ Solutions connects to these real-time computers to collect and process data, and send commands to actuators during plasma discharges. This system handles the signal IO from the tokamak and allows the PCS to utilize the diagnostic data necessary to perform control in real-time. Feedback control on heat flux was accomplished by feeding infrared camera data from the “Infrared TV” (IRTV) camera to a custom-developed User Datagram Protocol (UDP) server. This server transmits infrared data to a newly developed PCS algorithm that estimates the heat flux to PFC. Here, a proportional integral derivative (PID) controller minimizes the error between a heat flux reference and the real-time estimate by injecting nitrogen gas into the divertor.
Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.
Solutions to the inverse heat conduction problem (IHCP) are methods that can be used to quantify surface heat flux in multi-layer materials for components in which there are limited subsurface (internal) temperature measurements, such as coated components. A critical consideration is to capture high frequency fluctuations using a practical heat flux sensor. To that end, this paper highlights key parameters for calculating accurate surface heat transfer. Specifically, this research extends the available solutions to the IHCP for multi-substrate structures through an impulse response methodology. The sensitivity of the impulse method was quantified with respect to practical measurements. Further, when compared to the inverse case, the impulse method resulted in lower errors when calculating surface heat flux over a range of conditions. Overall, this work provides a foundation for deducing heat flux from a subsurface heat flux sensor while maintaining a high-frequency response.
Near-resonant energy transfer to large-scale stable modes is shown to reduce transport above the linear critical gradient, contributing to the onset of transport at higher gradients. This is demonstrated for a threshold fluid theory of ion temperature gradient turbulence based on zonal-flow-catalyzed transfer. Furthermore, the heat flux is suppressed above the critical gradient by resonance in the triplet correlation time, a condition enforced by the wave numbers of the interaction of the unstable mode, zonal flow, and stable mode.
Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.
In this project, robust heat-flux measurement systems were developed. The heat-flux sensors utilize thermoelectric effects to directly transduce the heat-flux inputs to analog electrical voltage signals. They were constructed from dedicated materials that can withstand temperatures of at least 1000°C and maintain adequate performance at these conditions for prolonged periods of time. The proposed approaches took into account numerous considerations, including system cost, sensor head resilience, sensor footprint, data accuracy, response time, and maintenance requirements. Through modern thermoelectric materials design, methodical materials selection and rigorous testing in materials characterization labs and medium-scale fire research facilities, we have demonstrated functioning laboratory prototypes, upon which one could base industrial heat-flux sensing platforms capable of operating in the challenging high-temperature, corrosive environments of the boilers of coal-fired power plants. A distributed sensor array for heat-flux measurements throughout the furnace water-wall, the superheater area and the economizer coils can provide critical data for the power plant control systems to increase efficiency, improve safety and reduce down times. For example, the combined heat-flux sensor/control systems can contribute to the optimization of burner and boiler operations under flexible loads, the optimization of heat-exchange conditions and overall reduction of heat rate and emissions, the prediction of imminent overheating conditions, and the optimization of the soot-blowing protocols.
We report simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.
Abstract In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.
For the reliable long-term operation of fusion power plants, it is crucial to understand and predict the lifetime of materials in use. These materials include all structural and functional materials utilized at the first wall, blanket, magnets, and shielding. The key challenge is, that the harsh environment including high heat fluxes, high thermal stress and stress cycling, neutron irradiation, and sputtering on such materials should not be viewed separately. Currently, the synergistic loads cannot be evaluated experimentally because of the lack of adequate facilities. The purpose of that work is to design a synergetic Neutron Irradiation and Thermomechanical Experiment (NITE) for fusion materials. This design will leverage the existing Advanced-Test-Reactor (ATR), a fission reactor at the Idaho National Laboratory. We also acknowledge that with existing fission reactors the exact fusion condition cannot be created, and the limitations are critically discussed. The combination of neutron irradiation with a high heat flux is the focus. This is realized with an irradiation capsule design that includes a TRISO fueled region inside the capsule to enable a steady-state heat flux on one side of the specimen. In conclusion, the experimental design modeling showed that steady-state heat fluxes of 2.4 MW/m 2 with a thermal gradient of above 250°C can be achieved in a 5 mm thick specimen.