Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive simulation”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Kinematically Hot, Extremely Metal-Poor C-19 Stellar Stream in DESI DR2

Stellar streams are the result of a host galaxy's gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a globular cluster. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicites for over 10 million stars reaching significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$ km s$^{-1}$ and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel 'spur' feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a 'hot', metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Mohammed, Nasser [Toronto U.] (ORCID:0009000812240↗
DOE OSTI2026

The Kinematically Hot, Extremely Metal-Poor C-19 Stellar Stream in DESI DR2

Stellar streams are the result of a host galaxy's gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a globular cluster. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicites for over 10 million stars reaching significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$ km s$^{-1}$ and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel 'spur' feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a 'hot', metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Mohammed, Nasser [Toronto U.] (ORCID:0009000812240↗
DOE OSTI2026

Peering down the barrel with DESI DR2: 10 000+ inflows at $z$ < 0.6 reveal how galaxies accrete cold gas

Direct observational constraints on how galaxies acquire their gas remain remarkably limited, hindering our understanding of the baryon cycle. We present a search for down-the-barrel NaI D absorption towards 15.6 million galaxies at $z < 0.6$ in DESI Data Release 2. We use Bayesian evidence ratios to assess whether the absorption requires additional components tracing interstellar gas distinct from the systemic component of the galaxy. We construct a catalogue of 50 088 (27 420) galaxies with moderate (strong) evidence for down-the-barrel absorption. The inferred absorption components are broadly distributed in velocity, with approximately 50% at $v_{\rm flow} < -50$ km/s, 30% within 50 km/s of the systemic velocity and the remaining 20% at $v_{\rm flow} > 50$ km/s. We find strong evidence for a large population of low-velocity, infalling absorbers with velocities $\sim$20 km/s in edge-on galaxies, consistent with radial inflows predicted in simulations. The stronger correlation in early-type galaxies between inflow velocity and stellar velocity dispersion, compared to that with stellar mass, suggests that a portion of these inflows may be associated with accreting satellites. These results reveal the multiple pathways in which galaxies accrete gas at redshift $z < 0.6$ for the first time in a statistically significant sample.

Weng, S. [Marseille, Lab. Astrophys.]↗
DOE OSTI2026

Cosmic Pairs: A DESI Census of Dual and Offset AGN as Precursors to Massive Black Hole Binaries

We present a systematic census of dual and offset active galactic nuclei (AGN) using spectroscopic data from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). After correcting for observational systematics, our final sample contains $>7,000$ dual AGN and 27,000 galaxy pairs containing one AGN over the redshift range $0 \lesssim z \lesssim 3.6$. This sample expands the known dual AGN sample by $\sim 1-2$ orders of magnitude at $0.2 \lesssim z \lesssim 0.4$, includes $\sim 50$ dwarf dual AGN candidates in a regime where only a handful were previously known, and triples the census at $z>2$. Dual AGN are preferentially found at small separations, consistent with merger-driven triggering of AGN activity. The two members of a pair differ in their star formation response: the more massive (primary) host changes little with separation, while the less massive (secondary) lies $\sim 0.3$ dex above matched inactive and one-AGN companions at the same projected separation in main-sequence offset. Using ASTRID simulations, we predict that the fraction of DESI dual AGN whose central black holes will merge by $z \sim 0$ increases with redshift, reaching $\sim 76\%$ by $z \sim 2$, while the fraction producing LISA-detectable mergers peaks at $\sim 37\%$ near $z \sim 0.9$. These results provide the largest uniformly selected spectroscopic sample of kpc-scale dual and offset AGN candidates from a single survey, connecting their host-galaxy and AGN demographics to the progenitor population of massive black hole mergers detectable by LISA.

Dadiani, Ekaterine [Carnegie Mellon U.] (ORCID:000↗
DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Model Form Error Correction for a Black-Box Thermal Battery Heat Transfer Simulation

Thermal batteries are crucial for supplying power to high-consequence engineering applications such as rockets. Computational simulations have been developed to predict thermal battery behavior, but these simulations often suffer from modeling errors, including model form uncertainty. Addressing this uncertainty can be achieved by quantifying either the model discrepancy in the output or the model form error (MFE) in the governing equation. MFE is particularly valuable as it can be better extrapolated beyond observed outputs, which is essential for predictions involving changes in external system loading, system configuration and geometry, or output quantities. This paper employs a state estimation approach to estimate MFE using experimental data and then utilizes machine learning (ML) to model its relationship with state variables. A nonintrusive technique is used to estimate MFE in a black-box thermal battery heat transfer simulation. The trained machine learning model for MFE is then applied to correct simulation predictions under extrapolated initial conditions and battery configurations. In conclusion, the methodology's performance is evaluated using additional experimental data, demonstrating its effectiveness in improving prediction accuracy.

Batteries↗
DOE OSTI2024

Effects of Size on Characteristics of Electroslag Remelted (ESR) HAYNES®282® alloy ingots

The effects of varying ingot diameters on the characteristics of round ESR ingots of 282 alloy are evaluated through a combined experimental and modeling study. Ingots of four different diameters --100 mm, 150 mm, 200 mm and 500 mm, were sectioned and analyzed using metallography to assess spatial variations in dendritic structure (primary / secondary dendrite arm spacings – PDAS/SDAS) across these ingots. Melt parameters of the processes were input into the simulation code, MeltFlow-ESR™ to theoretically predict PDAS/SDAS in the ingots. The values for PDAS/SDAS predicted by the simulations are compared with those measured by experiments and quantitative differences in these values are discussed in the context of the capability and limitations of the model and measurement methods. Both the predicted and measured values show that the PDAS/SDAS values increase with increasing ingot size. In addition, the relation with other ingot characteristics, such as local solidification time as well as the likelihood for freckle formation are discussed. The implication of these results for down-stream ingot processing and for ingot soundness are discussed.

Krishnamurthy, Ramanathan↗
DOE OSTI2023

Electron–Ion Temperature Ratio in Astrophysical Shocks

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T e /T p inferred from the Hα profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T e /T p for supernova remnant shocks using Hα profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the Hα determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks, and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence, and varying preshock conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improving deep learning performance for predicting large-scale geological ${{CO}_{2}}$ sequestration modeling through feature coarsening

Physics-based reservoir simulation for fluid flow in porous media is a numerical simulation method to predict the temporal-spatial patterns of state variables (e.g. pressure p) in porous media, and usually requires prohibitively high computational expense due to its non-linearity and the large number of degrees of freedom (DoF). This work describes a deep learning (DL) workflow to predict the pressure evolution as fluid flows in large-scale 3-dimensional(3D) heterogeneous porous media. In particular, we develop an efficient feature coarsening technique to extract the most representative information and perform the training and prediction of DL at the coarse scale, and further recover the resolution at the fine scale by spatial interpolation. We validate the DL approach to predict pressure field against physics-based simulation data for a field-scale 3D geologic CO 2 sequestration reservoir model. We evaluate the impact of feature coarsening on DL performance, and observe that the feature coarsening not only decreases the training time by >74% and reduces the memory consumption by >75%, but also maintains temporal error 0.63% on average. Besides, the DL workflow provides predictive efficiency with 1406 times speedup compared to physics-based numerical simulation. The key findings from this research significantly improve the training and prediction efficiency of deep learning model to deal with large-scale heterogeneous reservoir models, and thus it can also be further applied to accelerate workflows of history matching and reservoir optimization for close-loop reservoir management.

58 GEOSCIENCES↗
DOE OSTI2022

Atmospheric teleconnection associated with the Atlantic multidecadal variability in summer: assessment of the CESM1 model

Observational analysis shows that the Atlantic multidecadal variability (AMV) is associated with climate variability in the Northern Hemisphere through a zonal atmospheric teleconnection extending from the North Atlantic Ocean and propagating eastward around the Northern Hemisphere. Here, we studied the fidelity of model simulations in reproducing the observed summer AMV and the associated impacts on the mid-latitude climate by analysing simulations using the National Centre for Atmospheric Research Community Earth System Model Version 1 (CESM1), including CESM1 North Atlantic idealized and pacemaker simulations, CESM1 large ensemble twentieth century uninitialized simulations and large ensemble initialized CESM1 decadal predictions. To further compare the fidelity of CESM1, we also analysed large ensemble simulations from three other models. Our results suggest that the uninitialized large ensemble simulations from all models can produce an AMV time evolution and its regional climate impacts similar to the observations to certain degree. By initializing the observed oceanic condition in decadal prediction simulations, the simulated AMV and its regional impacts are closer to the observed ones than those in uninitialized ensemble simulations. In addition, the pacemaker simulations that nudged the time-evolving observed North Atlantic sea surface temperature anomalies produce spatiotemporal characteristics of the AMV and AMV climate impacts closer to the observed ones than the uninitialized simulations. We conclude that although coupled models can produce AMV and its regional impacts similar to observed, proper initialization and bias correction of the sea surface temperature spatial and temporal structure can improve this capability.

Initialization↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
DOE OSTI2022

Dynamic force and stability prediction for milling using feed rate scheduling software and time-domain simulation

This paper describes: 1) the use of feed rate scheduling software to predict the radial depth of cut variation for three-axis milling toolpaths and; 2) the use of the radial depth profile in a time-domain simulation to predict dynamic cutting forces. The time-domain simulation, which also includes the tool tip frequency response functions and force model (which relates the cutting force components to the chip geometry) as inputs, enables dynamic force profiles to be predicted and parameter combinations that cause chatter to be identified. A ramp geometry is selected that provides constantly varying radial depth and force predictions are completed at multiple axial depths for comparison to measured forces. Both stable and unstable (chatter) milling conditions were observed with good agreement between time-domain simulation and measurement results. The value of combining the feed rate scheduling software and time-domain simulation is demonstrated.

42 ENGINEERING↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗