Search NASA⌕ Search

SEARCH · Search NASA

Results for “HEAT FLUX”

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 523 records · Page 29

Short-term impacts of 2017 western North American wildfires on meteorology, the atmosphere’s energy budget, and premature mortality

Abstract Western North American fires have been increasing in magnitude and severity over the last few decades. The complex coupling of fires with the atmospheric energy budget and meteorology creates short-term feedbacks on regional weather altering the amount of pollution to which Americans are exposed. Using a combination of model simulations and observations, this study shows that the severe fires in the summer of 2017 increased atmospheric aerosol concentrations leading to a cooling of the air at the surface, reductions in sensible heat fluxes, and a lowering of the planetary boundary layer height over land. This combination of lower-boundary layer height and increased aerosol pollution from the fires reduces air quality. We estimate that from start of August to end of October 2017, ∼400 premature deaths occurred within the western US as a result of short-term exposure to elevated PM 2.5 from fire smoke. As North America confronts a warming climate with more fires the short-term climate and pollution impacts of increased fire activity should be assessed within policy aimed to minimize impacts of climate change on society.

Bernstein, Diana N. (ORCID:0000000308279336)↗

Simulating the Effects of Surface Energy Partitioning on Convective Organization: Case Study and Observations in the US Southern Great Plains

Realistic cloud-resolving simulations were performed to study the effects of surface energy partitioning (surface sensible and latent heat fluxes) on the organization of isolated convection into larger mesoscale convective systems (MCSs) near the US Southern Great Plains. The role of cold pools in mediating surface-convection interactions was explored. Better organized MCSs tended to occur in the experiments with perturbed wetter soil (and more active vegetation), regardless of the effects of soil moisture on the diurnal timing of convective triggering. Wetter soil led to shallower boundary layers and more convective available potential energy than drier soil. The roles of cold pools on convection are lifecycle-stage dependent: A dry surface allows more numerous colliding cold pools, thereby aiding in convective triggering by reducing entrainment in the early stages, and providing gust-front uplift in later stages. However, horizontal propagation of the cold-pool density current can outrun the convective system, creating a slantwise updraft and thus weakening the gust front uplift in later stages. This effect calls into question previous cold pool parameterizations, in which the gust front uplift is mainly proportional to the negative buoyancy of cold air. Finally, both the model and observation show an enhancement of surface latent heat flux during the passage of a gust front at night, suggesting that a positive feedback between the surface and convection helps MCSs to persist into the nighttime.

54 ENVIRONMENTAL SCIENCES↗

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental evidence of enhanced radial transport in small ELM regimes at DIII-D

Small/type-II edge-localized-modes (ELMs), carrying 1% of the plasma stored energy, are found to deposit only 45 ± 5% of the ELM power near the strike point, and the remaining 55 ± 5% to the far scrape-off-layer (SOL). Small ELMs spread their power over a larger area compared to type-I ELMs, where such a ratio is about 60% and 40% for near- and far-SOL regions, respectively. The larger spread is reflected in the heat flux width (λq) in the SOL for the intra-small ELMs profile of 6.0 mm, almost a factor 2 larger than that of type-I ELMs of 3.15 mm, for similar plasma conditions and magnetic configuration. At the ELM peak, the small ELMs λq is found to be up to 4 times larger than for the type-I ELMs, going from 2 to 7.9 mm, indicating enhanced radial transport in the neon-seeded small ELM scenario. Inter-ELM λqs have been also calculated at the secondary outer divertor in quasi-double-null (QDN) discharges. It is found that, on average, λq is 2.2 times larger in the high-separatrix-density small ELM regime, compared to a reference type-I ELM one. These findings are supported by small ELMs radial velocity profiles, measured at the outer midplane with a fast reciprocating probe, showing a decay length (λvr) in the SOL of 12.8 cm, which is 3.3 times larger than that for the type-I ELMs of 3.9 cm. This analysis shows that small ELMs, although attractive for future machines due to low peak heat flux and large λq, might be of concern for the larger flux to the outer wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗

Compatibility of divertor detachment and ELM suppression in DIII-D high- β p plasmas with ITER-similar shape

Abstract Integration of transient and steady-state divertor heat fluxes control with a high-performance core is necessary for future fusion reactors. In recent DIII-D high- β p experiments, divertor detachment and simultaneous edge localized mode (ELM) suppression are demonstrated while the plasma confinement quality is maintained high in ITER-similar shape. By optimizing the neon injection in high- β p scenario with ITER-similar shape, deep detachment and ELM suppression are achieved with a high-performance core ( β N ∼ 2.8, β p ∼ 2.3) at q 95 ∼ 7.5. Partial divertor detachment and suppression of large ELMs are achieved at q 95 ∼ 6. The stability analyses suggest that with low neon injection, the density pedestal becomes higher and steeper and the T i profile also increases, therefore the increased edge pressure and higher current density destabilize the Peeling-Ballooning mode (PBM), which would lead to a large ELM collapse. With strong neon gas puffing, the significantly reduced pedestal pressure and current density, due to the degraded T e pedestal, lead to the stabilization of PBM and ELMs are suppressed. For both cases, the coupling between the large radius internal transport barrier (ITB) and edge pedestal is the key reason for maintaining high global performance. The formation of large radius ITB compensates for pedestal degradation. Such results could provide an attractive scenario to well control the transient and steady-state heat flux onto the divertor plates while maintaining good plasma performance, which is an important step toward the steady-state operation of future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope effects on transport in LHD

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

Physics↗

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantifying the impacts of atmospheric rivers on the surface energy budget of the Arctic based on reanalysis

We present a comprehensive analysis of Arctic surface energy budget (SEB) components during atmospheric river (AR) events identified by integrated water vapor transport exceeding the monthly 85th-percentile climatological threshold in 3-hourly ERA5 reanalysis data from January 1980 to December 2019. Analysis of average anomalies in SEB components, net SEB, and the overall AR contribution to both the seasonal SEB components and net SEB climatology reveals clear seasonality and distinct land–sea–sea ice contrast patterns. Over the sea-ice-covered central Arctic Ocean, ARs significantly impact net SEB, inducing substantial surface warming in fall, winter, and spring. This warming is primarily driven by large anomalies in surface downward longwave radiation (LWD), which average 29–45 W m −2 during the cold seasons. In contrast, AR-related LWD anomalies are smaller in summer, averaging around 15 W m −2 , indicating a reduced impact during this season. Over sub-polar oceans, ARs have the most substantial positive impact on net SEB in cold seasons, mainly attributed to significant positive turbulent heat flux anomalies. AR-related turbulent heat anomalies reduce the upward turbulent flux, contributing up to −11 % relative to its seasonal climatology. In summer, ARs induce negative impacts on net SEB, primarily due to reduced shortwave radiation from increased cloud cover during AR events. Over continents, ARs generate smaller absolute impacts on net SEB because the large LWD anomalies are largely offset by corresponding increases in upward longwave radiation, particularly during cold seasons. Additionally, the seemingly large relative contributions of ARs to the net SEB over land primarily reflects the small magnitude of the climatological net SEB over continents. Greenland, especially western Greenland, exhibits significant downward longwave radiation anomalies associated with ARs, which drive large net SEB anomalies and contribute >54 % to mean SEB and induce amplified surface warming year-round. This holds significance for melt events, particularly during summer. Additionally, results of AR-related SEB impacts strongly depend on detection methods, as restrictive AR detection algorithms that emphasize extreme AR events, with large AR-related anomalies, do not necessarily indicate a large overall contribution to the SEB climatology due to the low occurrence frequency of these events. This study quantifies the role of ARs in the surface energy budget, contributing to our understanding of the Arctic warming and sea ice decline in ongoing Arctic amplification.

Arctic Sea ice change↗

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continuous observations of the surface energy budget and meteorology over the Arctic sea ice during MOSAiC

Abstract The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) was a yearlong expedition supported by the icebreaker R/V Polarstern , following the Transpolar Drift from October 2019 to October 2020. The campaign documented an annual cycle of physical, biological, and chemical processes impacting the atmosphere-ice-ocean system. Of central importance were measurements of the thermodynamic and dynamic evolution of the sea ice. A multi-agency international team led by the University of Colorado/CIRES and NOAA-PSL observed meteorology and surface-atmosphere energy exchanges, including radiation; turbulent momentum flux; turbulent latent and sensible heat flux; and snow conductive flux. There were four stations on the ice, a 10 m micrometeorological tower paired with a 23/30 m mast and radiation station and three autonomous Atmospheric Surface Flux Stations. Collectively, the four stations acquired ~928 days of data. This manuscript documents the acquisition and post-processing of those measurements and provides a guide for researchers to access and use the data products.

54 ENVIRONMENTAL SCIENCES↗

Thermal shock resistance of multilayer silicon carbide receiver tubes for 800°C molten salt concentrating solar power application

CSP power tower receiver systems during rapid transient weather periods can be vulnerable to thermal shock conditions from rain that which can facilitate the onset of leaks and failures that can have catastrophic consequences. Silicon carbide (SiC) materials have attractive receiver application characteristics for being light weight, having high-strength and excellent thermal shock resistance performance which make them a particularly good fit for receiver absorber materials in CSP. In this investigation, the performance characteristics of Ceramic Tubular Products (CTP) SiC ceramic matrix composite (CMC), multilayered tubes were explored with respect to thermal shock performance for solar receiver applications in next generation CSP plants. Here, thermal shock testing was performed at the Sandia National Laboratories (SNL) Solar Furnace facility using a dynamic stage and thermal shock tube test setup. The tubes tested under incident solar heat flux of 100 W/cm2 were heated with inner tube temperatures reaching approximately 800 °C, with outer temperatures exceeding or just reaching 1000 El for the multilayer and monolithic SiC tubes respectively. The tubes were then quenched with simulated rain. The tubes were then cooled and subjected to hoop stress analysis using an Instron device to assess their subsequent mechanical strength. Furthermore, the on-sun study experimental results indicate an average of 24.2% and 97% higher hoop strength for the CMC tubes than those composed of monolithic SiC and aluminum oxide (A1203) respectively.

14 SOLAR ENERGY↗

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗

Gyrokinetic prediction of core tungsten peaking in a WEST plasma with nitrogen impurities

Tungsten peaking is predicted in the core of a WEST plasma with total-f gyrokinetic simulations, including both collisional and turbulent transport. This prediction is validated with a synthetic diagnostic of the bolometry. Although nitrogen impurities are shown to reduce the neoclassical peaking of tungsten on-axis, the overall tungsten peaking increases when nitrogen impurities are present, as they reduce the turbulence screening off-axis. This finding is important for the plasma current ramp-up phase of ITER, where light impurities seeding will be desirable to achieve low temperatures at the plasma-facing components and reduce tungsten sputtering. It provides further argument for applying early ECRH heating to maintain margins on the core power balance. The neoclassical peaking factor is cross-verified between XGC and FACIT. The heat flux at separatrix and the heat load width are modeled by XGC and compared to WEST data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of the KSTAR experiments and future plan

The Korea Superconducting Tokamak Advanced Research (KSTAR) device has recently entered a new phase of operation following the installation of a tungsten mono-block divertor, providing a unique platform to investigate high-performance, long-pulse plasma scenarios for future reactors. The 2023–2024 experimental campaigns demonstrated significant progress in commissioning the tungsten divertor, developing plasma operation scenarios, and achieving real-time control. Comparative studies revealed increased core radiation losses and associated performance degradation relative to the carbon divertor environment. Alleviation strategies, such as optimized neutral beam injection timing, boron powder seeding, and impurity control through gas puffing, were shown to reduce tungsten accumulation and extend high-performance phases, including record H-mode operation exceeding 100 s. International collaboration enabled the first demonstration of high poloidal beta scenarios with an internal transport barrier on KSTAR. Advances in control included upgraded real-time resonant magnetic perturbation schemes and disruption forecasting. Furthermore, divertor detachment control using a new surrogate-model-based approach with real-time radiation imaging demonstrated active impurity and heat flux management. Supported by enhanced heating, current drive, and diagnostic systems, these achievements offer valuable insights into tungsten impurity behavior, transport physics, and control strategies for future reactors including ITER. Planned upgrades toward a full-tungsten wall and ITER-aligned real-time plasma control will further establish KSTAR as a leading experimental platform for developing operational scenarios required to achieve steady-state, burning plasma conditions.

KSTAR↗