Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fusion reactors”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Taming harmful bursts and heat flux in high-confinement tokamak plasmas

A major challenge in tokamak fusion research is first-wall erosion caused by steady heat loads and sudden energy bursts known as edge-localized modes. Divertor detachment reduces steady-state heat flux, while resonant magnetic perturbations can suppress these instabilities. However, integrating the two has been difficult because they require conflicting operating conditions. Here we demonstrate simultaneous achievement of resonant magnetic perturbations mitigated small edge-localized modes and impurity seeded partial divertor detachment in plasmas with an ITER-similar shape on the DIII-D tokamak. Experiments and simulations show that resonant magnetic perturbations facilitate detachment by redistributing particles, lowering the core density and increasing the scrape-off layer density, thereby reducing the amount of injected gas required. Cooling-gas injection eliminates the secondary heat-flux peak created by three-dimensional magnetic lobes, while edge cooling weakens the plasma response to the applied magnetic fields. These advances illustrate a viable pathway for integrating edge stability control with power exhaust in future fusion reactors.

magnetically confined plasmas↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing Stainless Steel Compatibility in Flowing FLiBe Salt

To address the considerable interest in LiF-BeF 2 (FLiBe) compatibility for fission and fusion reactor applications, static and flowing compatibility experiments were conducted to assess the compatibility with type 316H stainless steel. In static testing at 550° and 650 °C, small mass changes were measured and posttest characterization of the FLiBe showed increased levels of Fe, Cr, Ni and Mn in the salt. Adding Be in the static salt test reduced the dissolution of Fe and Ni. Further, an initial assessment of mass transfer in flowing FLiBe without a Be addition was conducted using a monometallic 316H thermal convection loop (TCL) operated for 1000 h with a peak temperature of 650 °C. Similar to prior results in flowing FLiNaK salt, the 316H specimens exhibited small mass losses in the hot leg. Posttest characterization of the 316H specimens suggested Cr surface depletion in the hot and cold legs and possibly Fe deposition in the cold leg. To further understand this behavior, Cr and Fe dissolution was measured in static FLiBe at 550–650 °C.

316H stainless steel↗
DOE OSTI2025

Deuterium trapping mechanisms in reduced activation ferritic martensitic steels and their correlation with mechanical strengthening

Development of high-strength materials often involves introduction of additional strengthening microstructures that also serve as tritium trapping sites. Such additions in fusion material development could degrade the fuel efficiency in fusion reactors and raise radiological concerns. The contribution of individual microstructure features in hydrogen trapping must be evaluated to ensure fuel efficiency and radiological safety. This study explores the mechanistic origins of deuterium trapping in reduced-activation ferritic–martensitic steels and its correlation to mechanical strengthening. A series of model alloys and engineering steels were fabricated and subjected to different heat treatments to control deuterium trapping site density. Deuterium retention was evaluated using D 2 gas charging and thermal desorption spectroscopy, focusing on the role of grain boundary, dislocation, M 23 C 6 precipitates, and TiC precipitates. Multiscale microstructure characterization and synchrotron X-ray diffraction were performed to characterize microstructure, which was correlated to the deuterium retention property. Results show that TiC precipitates exhibit the highest deuterium trapping capacity, followed by M 23 C 6 precipitates. Dislocation and grain boundary demonstrate the lowest and similar efficiencies. Furthermore, the relationship of trapping quantity and mechanical strengthening of these microstructure features was quantified, demonstrating that TiC precipitates offer highest deuterium trapping per unit of mechanical strengthening.

Retention↗
DOE OSTI2025

Vacuum vessel design with lofted toroidal surfaces for a QHS configuration

A vacuum vessel in between the hot plasma and the magnetic coils is crucial for a fusion reactor. In this study, a tool is presented to perform the initial design of such a vessel automatically. Starting from the last closed flux surface, an algorithm is developed which lofts this manifold out as much as possible. At the same time, the algorithm makes sure that no intersection with the coils takes place. In the presented work, this algorithm is applied to the HSX stellarator and it has been verified that the resulting vessel wall can be built from an engineering perspective. Using field line tracing, the strikepoints on the vessel wall indicate the initial locations for a divertor.

HSX↗
DOE OSTI2025

Implementation of high-speed data acquisition at DIII-D

Research at the DIII-D National Fusion Facility in San Diego focuses on short pulse plasma discharges that specialize on various shaping profiles. High-speed data collection is a critical component for the operation of many of DIII-D’s diagnostics and is fundamental for capturing high-resolution data used in experimental data analysis. Differing techniques enable the plasma control system (PCS) to perform complex real-time feedback control on microsecond time scales. This work presents a comprehensive overview of data acquisition, focusing on the hardware and software used in reliable data acquisition at DIII-D. The robust nature of the data acquisition system allows for various techniques to coexist seamlessly. However, as modern systems capable of nanosecond resolution become more common, existing architectures will need to be modified. Here, by addressing the key challenges of high-speed data acquisition, DIII-D is able to provide real-time data used in plasma operation and has the ability to acquire high fidelity data needed for future experimental fusion reactors, such as ITER.

Control↗
DOE OSTI2026

Corroborating VNA and thermal measurements of transmission loss on the DIII-D ECH waveguide system

Electron cyclotron heating (ECH) and current drive (ECCD) will play a large role in tokamak-based fusion reactors. At the DIII-D tokamak, 110 GHz microwaves injected into the plasma can provide core heating and current drive as well as impurity control, neoclassical tearing mode mitigation, and breakdown assistance. Understanding the physics of these processes relies on accurate estimates of injected ECH power. DIII-D’s ECH system consists of six MW-class Microwave Power Products (MPP) gyrotron microwave sources. Operating the gyrotrons far from the tokamak removes them from magnetic field interference, so 31.75 mm inner-diameter corrugated waveguides transmit the microwave power the 80 m from the gyrotrons to steerable launchers in the tokamak chamber. Estimates of injected power rely on knowing the generated power at the source and then subtracting transmission loss. Conventional transmission loss measurements based on calorimetric dummy loads are onerous and only possible during extended maintenance periods. This work examines two tools that provide more flexibility for the transmission loss measurements. Furthermore, a resistive temperature detector (RTD) array installed along a waveguide measures heat lost to the transmission line, and low power time domain reflectometry (TDR) measurements with a vector network analyzer (VNA) allows loss measurements without burdensome hardware modifications.

ECH↗
DOE OSTI2024

Characteristics of oxide-dispersion strengthened alloys produced by high-temperature severe deformation

This study is to explore an economically attractive and technically feasible processing method for oxide-nanoparticle strengthened alloys for fusion reactor application. Despite many scientific merits of the advanced oxide-dispersion strengthened (ODS) alloys, such as the nanostructured ferritic alloy (NFA) 14YWT, the only viable production path for a high-quality NFA is the high-power mechanical alloying process. This process is often a multi-day high-speed ball milling of alloy powder with a small quantity of yttria (Y 2 O 3 ) powder, followed by the milled-powder consolidation using extrusion or other methods and additional thermomechanical processing (TMP) for property control. This complex production path, including the low-temperature mechanical alloying in particular, has limited technical advancement toward the cost-effective and scale-up production of ODS alloy components. To overcome such a practical limitation, we proposed to explore alternative low-cost processing routes using traditional thermomechanical processing (TMP) method only. Further, a series of continuous TMP cycles, which were designed to impose high-temperature severe plastic deformation (HT-SPD) conditions to the consolidated powder mixtures, were applied to achieve the effective distribution of oxide particles in nanograin structure and thus desirable mechanical properties. Since the reduced-activation ferritic-martensitic (RAFM) alloy powders (Fe-10Cr and Fe-14Cr alloys with various Y contents) are available in our inventory, we focused to utilize the new solid-state synthesis approach for controlling oxide (oxygen source) dissolution and nanoscale clustering in nanograin structure in those alloys. A combination of powder consolidation at 900 °C and continuous thermomechanical activation at 600 °C yielded two essential ODS alloy microstructure contents–nanograin structure and nanoparticle distribution–and thus demonstrated a good combination of strength and ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The interaction of edge dislocations with hydrogen-helium bubbles in tungsten

In fusion reactors, plasmas-facing components undergo degradation due to displacement damage and injection of gas impurities. Significant amounts of hydrogen (H) and helium (He) impurities can be introduced into materials through plasma exposure and nuclear transmutation, and their synergistic effects lead to the formation of mixed He/H bubbles under irradiation, changing their hardening mechanisms. In this study, molecular dynamics (MD) models were developed to study the interaction between ½<111> edge dislocations with mixed He/H-bubbles in W. The critical resolved shear stress (CRSS) that is required for dislocation breakaway increases with He/H-to-vacancy ratio at 600 K, indicating strong pinning effects of mixed He/H bubbles. However, as the temperature increases up to 1400 K, the mixed He/H bubbles become unstable as the H atoms are increasingly emitted from the mixed He/H bubbles and migrate into W matrix, which lowers bubble pressure and CRSS. Overall, these results highlight the synergistic effects of He and H on dislocation behavior in irradiated W during deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2024

ToF-SIMS in material research: A view from nanoscale hydrogen detection

Hydrogen in materials has attracted tremendous interest as its incorporation leads to significant alterations in nanoscale structure, composition, and chemistry, impacting functional properties. It has also been integral to nuclear fusion reactors and is considered a future clean energy source. However, nanoscale characterization and manipulation of hydrogen in materials are challenging as only a selected few analytical techniques can readily detect hydrogen, among which time-of-flight secondary ion mass spectrometry (ToF-SIMS) is a unique and powerful one due to its excellent detection limit along with decent depth and lateral resolutions. Here, in this review, we discuss, using selected examples, how to detect and quantify hydrogen in materials by ToF-SIMS and its impact on revealing the hydrogenation/protonation-induced novel functional states in different classes of materials. In addition, we present our protocols on sample preparation and experimental conditions optimization, allowing us to achieve the best possible results. Finally, we highlight future research directions that can lead to the discovery of novel functional states and ultimately provide a deeper understanding of scientific questions in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2026

In situ studies on heavy ion irradiation and partial oxidation induced stacking faults and nanograins in tungsten nanolaminates

Understanding the microstructural evolution of tungsten (W) under extreme irradiation environments is critical for its application as a plasma-facing material in future fusion reactors. Heavy ion irradiation study offers expedited damage accumulation and thus allows irradiation tolerance property prediction in a short period of time. In this study we explore in situ Kr ion irradiation of W nanolaminates at 800 °C and uncover a complex interplay of irradiation-induced transformations, including the emergence of stacking faults, nanograin formation, and generation of thickened grain boundaries defined as GB regions with the thickness of up to 15 nm. High-resolution transmission electron microscopy studies reveal the formation of extended planar faults and a metastable hexagonal close-packed (hcp) phase within the body-centered cubic (bcc) matrix. These structural transitions are facilitated by irradiation-induced shear. They are further stabilized by the presence of oxygen, which promotes stacking fault formation and vacancy trapping. Elevated temperatures enhance defect mobility, enabling dynamic recrystallization into ultrafine grains and the thickening of GBs due to defect absorption and impurity segregation. These findings unravel non-equilibrium phase transformation pathways in heavy ion irradiated W and highlight the critical role of impurity-mediated defect dynamics in governing its radiation tolerance properties.

Wazeer, Adil [Purdue University]↗
DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗
DOE OSTI2026

A dual dynamic shutter system for accelerating ion irradiation sample throughput via lateral gas implantation gradients

Ion irradiation for material performance testing is limited due to its serial nature, which allows for only one value of the implantation (appm) versus dose (dpa) parameter space to be explored for each ion and experiment at a time. While ion irradiation can accelerate the process by up to three orders of magnitude compared to neutron irradiation experiments, the sample throughput for ion irradiation remains relatively low. To address these limitations, a novel capability has been developed at the Michigan Ion Beam Laboratory (MIBL), enabling for the creation of single- and two-dimensional lateral ion implantation gradients using recently installed motorized-controlled ion-beam shutters. This advancement can generate a wide scope of the two-dimensional (H+, He2+) implantation parameter space within a single sample. Integration of this new capability now allows for dual- and triple-ion beam experiments to be performed with full user control over not only the ion implantation depth, but also laterally across the sample by imposing ion implantation concentration gradients, thus providing researchers with a high-throughput means for material testing under various irradiation conditions. Furthermore, recent improvements in MIBL's microbeam ion-beam analysis (IBA) target station now allow for probing these concentration gradients in irradiated alloys with exceptional spatial resolution, down to 10 µm. These two approaches promise to significantly improve ion irradiation capabilities and increase the sample throughput by several orders of magnitude. The application of the shutter technique plus the subsequent microbeam characterization of the imposed implantation gradients are showcased by two proof-of-principle ion-irradiated experiments, one performed on single-crystal Si and the other on the fusion-candidate alloy F82H-IEA. These advancements mark a substantial leap in ion-beam technology, offering researchers a robust, high-throughput method to efficiently investigate candidate alloys with high technological readiness for both advanced fission and fusion reactor applications, in a time- and cost-effective manner.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Surface normal emissivity measurements in uncoated and oxidized B-coated plasma facing component materials in the long wavelength infrared spectrum

Accurate surface emissivity measurements are integral to the fidelity of the infrared thermography diagnostic evaluation in fusion reactors. The emissivities of ATJ TM graphite, Sigrafine® (R6510) graphite, and TZM alloy were measured as a function of temperature in accordance with the Contact Thermometer Method in the ASTM E1933 standard. Samples were heated resistively in a high-vacuum chamber, and the surface temperature was monitored using a surface thermocouple and a Telops long-wavelength infrared camera. The surfaces of Sigrafine® graphite and TZM alloy samples were coated with 10 nm and 20 nm layers of oxidized boron and were also measured in the uncoated condition. The thickness was assessed from measurement by a quartz crystal microbalance (QCM). ATJ TM graphite was uncoated. The measured emissivity of ATJ TM Graphite ranged from 0.82 to 0.83, uncoated Sigrafine® graphite ranged from 0.76 to 0.98, and uncoated TZM alloy ranged from 0.11 to 0.13. Oxidized boron coatings increased the emissivity of the TZM alloy to 0.18–0.23 but decreased the emissivity of Sigrafine® graphite to 0.56–0.66 for a 10 nm coating and to 0.65–0.74 for a 20 nm coating. In contrast to the expected blackbody radiance, the emissivity of uncoated Sigrafine® graphite and TZM alloy did not monotonically increase with temperature from 100 to 500 ℃.

Emissivity↗