Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Joint EM-NE-International Study of Glass Behavior over Geologic Time Scales - 12303

Vitrification has been chosen as the best demonstrated available technology for waste immobilization worldwide. To date, the contributions of physical and chemical processes controlling the long-term glass dissolution rate in geologic disposal remain uncertain, leading to a lack of international consensus on a rate law for glass corrosion. Existing rate laws have overcome uncertainty through conservatism, but a thorough mechanistic understanding of waste form durability in geologic environments would improve public and regulator confidence. If it is possible to take credit for the true durability of the waste form in repository system evaluations, then it is possible to design the repository with less conservatism with concomitant cost savings. To gain a fundamental understanding of the dissolution rate law, six nations have joined together to formulate a joint plan for collaborative research into the mechanisms controlling the long-term corrosion of glass. This report highlights the technical program plan behind the US portion of this effort, with an emphasis on the current understanding (and limitations) of several mechanistic theories for glass corrosion. Some recent results are presented to provide an example of the ongoing research. Atom probe tomography has been used to provide a high-resolution analysis of elemental concentration gradients present at the hydrated glass / pristine glass interface in SON68 after 25.75 years of corrosion in a simulated granitic groundwater at 90 deg. C. The most valuable result of these initial studies is the success of the technique. Characterization by APT had never been previously demonstrated for glass corrosion layers. The resolution of APT is a powerful addition to the tools with which we can investigate the mechanisms dominating glass corrosion. Some other key results of this study include the observation that the elemental interfacial width between the hydrated glass and pristine glass appears to be much sharper (∼2 nm for B, Na and Al) than had been previously measured using nanoSIMS (∼240 nm). It is not clear whether the APT analysis and nanoSIMS characterizations were possibly performed on topographically unique regions, or whether nanoSIMS overestimated the elemental width. However, the APT data seems very convincing that the elemental width can be much sharper than was previously thought. This result calls into question some of the assumptions made for the diffusion-control models of glass dissolution, since such a sharp profile would not match the diffusion coefficients used to date. Other results, such as the observation of apparently layered concentration profiles, show that gel evolution is not as simple as is currently assumed in nearly every model. This task is a good example of the collaborative nature of the I-TEAM effort. Based on experimental needs and differences in expertise, scientists from DOE and CEA worked together to change the level of understanding in the field. These types of interactions are nearly ubiquitous among the tasks in the technical program plan. With the excellence of the team in place and the willingness of the participants to work together for a common understanding, the stated goal of consensus on the mechanistic basis for radionuclide release from glass is well within reach. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Dual-pace transient heat conduction in vertically aligned carbon nanotube arrays induced by structure separation

Vertically aligned arrays of carbon nanotubes (CNTs) are attractive for a wide range of macroscopic applications which can exploit the remarkable properties of individual nanotubes. In this work, an abnormal behavior of CNT bundles in heat conduction is discovered under the transient electro-thermal characterization. The measured voltage change over the sample shows a dual-pace thermal response (DTR), which could not be fitted using a single thermal diffusivity heat transfer model. Instead, two thermal diffusivities are determined from the DTR phenomenon. Three rounds of cryogenic experiments are conducted to investigate the physics behind DTR phenomenon. It starts to emerge when the temperature is reduced to a certain level. After two rounds of cryogenic experiments, the DTR phenomenon becomes permanent from 295 K to 10 K. The nano-scale structure separation induced by the increased thermal strain leads to macroscale structure separation, which results in two parallel heat conduction paths responsible for the DTR phenomenon. By building a new parallel heat transfer model taking both the transient and steady-state electrical and thermal response into consideration, the area ratio of separated CNTs is determined. Furthermore, this result uncovers the existence of coiled morphology and nano-structure evolution of CNTs under cryogenic state and its effect on thermal and electrical transport, especially on their transient behaviors. Knowledge of the thermal transport in these CNTs arrays is important considering the fact that transient thermal response strongly affects their mechanical, optical, and electrical behaviors.

42 ENGINEERING↗
DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Theoretical limit of spatial resolution in diffuse optical tomography using a perturbation model

We have assessed the limit of spatial resolution of timedomain diffuse optical tomography (DOT) based on a perturbation reconstruction model. From the viewpoint of the structure reconstruction accuracy, three different approaches to solving the inverse DOT problem are compared. The first approach involves reconstruction of diffuse tomograms from straight lines, the second – from average curvilinear trajectories of photons and the third – from total banana-shaped distributions of photon trajectories. In order to obtain estimates of resolution, we have derived analytical expressions for the point spread function and modulation transfer function, as well as have performed a numerical experiment on reconstruction of rectangular scattering objects with circular absorbing inhomogeneities. It is shown that in passing from reconstruction from straight lines to reconstruction using distributions of photon trajectories we can improve resolution by almost an order of magnitude and exceed the accuracy of reconstruction of multi-step algorithms used in DOT. (optical tomography)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

FLARES IN THE CRAB NEBULA DRIVEN BY UNTWISTING MAGNETIC FIELDS

The recent discovery of PeV electrons from the Crab Nebula, produced on rapid timescales of one day or less with a sharply peaked gamma-ray spectrum without hard X-rays, challenges traditional models of diffusive shock acceleration followed by synchrotron radiation. Here, we outline an acceleration model involving a DC electric field parallel to the magnetic field in a twisted toroidal field around the pulsar. Sudden developments of resistivity in localized regions of the twisted field are thought to drive the particle acceleration, up to PeV energies, resulting in flares. This model can reproduce the observed timescales of T Almost-Equal-To 1 day, the peak photon energies of U{sub {Phi},rr} Almost-Equal-To 1 MeV, maximum electron energies of U{sub e,rr} Almost-Equal-To 1 PeV, and luminosities of L Almost-Equal-To 10{sup 36} erg s{sup -1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Diffuse-Interface Capturing Methods for Compressible Two-Phase Flows

Simulation of compressible flows became a routine activity with the appearance of shock-/contact-capturing methods. These methods can determine all waves, particularly discontinuous ones. However, additional difficulties may appear in two-phase and multimaterial flows due to the abrupt variation of thermodynamic properties across the interfacial region, with discontinuous thermodynamical representations at the interfaces. To overcome this difficulty, researchers have developed augmented systems of governing equations to extend the capturing strategy. These extended systems, reviewed here, are termed diffuse-interface models, because they are designed to compute flow variables correctly in numerically diffused zones surrounding interfaces. In particular, they facilitate coupling the dynamics on both sides of the (diffuse) interfaces and tend to the proper pure fluid–governing equations far from the interfaces. This strategy has become efficient for contact interfaces separating fluids that are governed by different equations of state, in the presence or absence of capillary effects, and with phase change. As a result, more sophisticated materials than fluids (e.g., elastic–plastic materials) have been considered as well.

42 ENGINEERING↗
DOE OSTI2015

Evaluation of Technetium Getters to Improve the Performance of Cast Stone

Cast Stone has been selected as the preferred waste form for solidification of aqueous secondary liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) process condensates and low-activity waste (LAW) melter off-gas caustic scrubber effluents. Cast Stone is also being evaluated as a supplemental immobilization technology to provide the necessary LAW treatment capacity to complete the Hanford tank waste cleanup mission in a timely and cost effective manner. One of the major radionuclides that Cast Stone has the potential to immobilize is technetium (Tc). The mechanism for immobilization is through the reduction of the highly mobile Tc(VII) species to the less mobile Tc(IV) species by the blast furnace slag (BFS) used in the Cast Stone formulation. Technetium immobilization through this method would be beneficial because Tc is one of the most difficult contaminants to address at the U.S. Department of Energy (DOE) Hanford Site due to its complex chemical behavior in tank waste, limited incorporation in mid- to high-temperature immobilization processes (vitrification, steam reformation, etc.), and high mobility in subsurface environments. In fact, the Tank Closure and Waste Management Environmental Impact Statement for the Hanford Site, Richland, Washington (TC&WM EIS) identifies technetium-99 ( 99 Tc) as one of the radioactive tank waste components contributing the most to the environmental impact associated with the cleanup of the Hanford Site. The TC&WM EIS, along with an earlier supplemental waste-form risk assessment, used a diffusion-limited release model to estimate the release of different contaminants from the WTP process waste forms. In both of these predictive modeling exercises, where effective diffusivities based on grout performance data available at the time, groundwater at the 100-m down-gradient well exceeded the allowable maximum permissible concentrations for 99 Tc. (900 pCi/L). Recent relatively short-term (63 day) leach tests conducted on both LAW and secondary waste Cast Stone monoliths indicated that 99 Tc diffusivities were at or near diffusivities where the groundwater at the 100-m down-gradient well would exceed the allowable maximum permissible 99 Tc concentrations. There is, therefore, a need and an opportunity to improve the retention of Tc in the Cast Stone waste form. One method to improve the performance of the Cast Stone waste form is through the addition of “getters” that selectively sequester Tc inside Cast Stone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Characterization and Modeling of Grain Boundary Chemistry Evolution in Ferritic Steels under Irradiation

Ferritic/martensitic (FM) steels such as HT-9, T-91 and NF12 with chromium concentrations in the range of 9-12 at.% Cr and high Cr ferritic steels (oxide dispersion strengthened steels with 12-18% Cr) are receiving increasing attention for advanced nuclear applications, e.g. cladding and duct materials for sodium fast reactors, pressure vessels in Generation IV reactors and first wall structures in fusion reactors, thanks to their advantages over austenitic alloys. Predicting the behavior of these alloys under radiation is an essential step towards the use of these alloys. Several radiation-induced phenomena need to be taken into account, including phase separation, solute clustering, and radiation-induced segregation or depletion (RIS) to point defect sinks. RIS at grain boundaries has raised significant interest because of its role in irradiation assisted stress corrosion cracking (IASCC) and corrosion of structural materials. Numerous observations of RIS have been reported on austenitic stainless steels where it is generally found that Cr depletes at grain boundaries, consistently with Cr atoms being oversized in the fcc Fe matrix. While FM and ferritic steels are also subject to RIS at grain boundaries, unlike austenitic steels, the behavior of Cr is less clear with significant scatter and no clear dependency on irradiation condition or alloy type. In addition to the lack of conclusive experimental evidence regarding RIS in F-M alloys, there have been relatively few efforts at modeling RIS behavior in these alloys. The need for predictability of materials behavior and mitigation routes for IASCC requires elucidating the origin of the variable Cr behavior. A systematic detailed high-resolution structural and chemical characterization approach was applied to ion-implanted and neutron-irradiated model Fe-Cr alloys containing from 3 to 18 at.% Cr. Atom probe tomography analyses of the microstructures revealed slight Cr clustering and segregation to dislocations and grain boundaries in the ion-irradiated alloys. More significant segregation was observed in the neutron irradiated alloys. For the more concentrated alloys, irradiation did not affect existing Cr segregation to grain boundaries and segregation to dislocation loops was not observed perhaps due to a change in the dislocation loop structure with increasing Cr concentration. Precipitation of α’ was observed in the neutron irradiated alloys containing over 9 at.% Cr. However ion irradiation appears to suppress the precipitation process. Initial low dose ion irradiation experiments strongly suggest a cascade recoil effect. The systematic analysis of grain boundary orientation on Cr segregation was significantly challenged by carbon contamination during ion irradiation or by existing carbon and therefore carbide formation at grain boundaries (sensitization). The combination of the proposed systematic experimental approach with atomistic modeling of diffusion processes directly addresses the programmatic need for developing and benchmarking predictive models for material degradation taking into account atomistic kinetics parameters

36 MATERIALS SCIENCE↗
DOE OSTI2012

Multispecies Diffusion Capability For The AMP Nuclear Fuel Performance Code (LANL Milestone M31MS060301 Final Report)

This work addresses only diffusion. The contact solver in AMP was not sufficiently developed this year to attempt treatment of species contact. A cylindrical tensor diffusion coefficient model was added to the AMP code, with the KHHS model [1] implemented into the AMP material library as a specific example. A cylindrical tensor diffusion operator manufactured solution verification example was coded. Before meeting the full text of the milestone task, it remains to: (1) code and run a cylindrical tensor diffusion solver manufactured solution (2) code and run the validation example of [1] (3) document results. These are dependent on developing new capabilities for the AMP code requiring close collaboration with the AMP team at ORNL. The model implemented provides a good intermediate first step toward a general multi-species solver. The multi-species capability of the AMP nuclear fuel code [2] is intended to allow the modeling of radiation-driven redistribution of various elements through solid metal nuclear reactor fuels. The initial model AMP provides for U-Pu-Zr fuels is based on the analysis of the Integral Fast Reactor (IFR) fuel development program experiment X419 post-irradiation data described in [1], referred to here as the KHHS model. This model may be specific to that experiment, but it was thought to provide a good start for the AMP code, because it (1) is formulated at the engineering scale, (2) decouples the species from each other, (3) predetermines the phase boundaries so that reference to a phase diagram is not needed, and (4) one of the authors (Hayes) was the NEAMS Fuels IPSC manager for FY11. The KHHS model is formulated for radial fluxes as little axial redistribution is seen experimentally. As U-Pu-Zr fuel is irradiated, the constituents migrate to form three annular regions. The center region is Zr-enriched and U-depleted, the middle region is Zr-depleted and U-enriched, and the outer region is Zr-enriched and U-depleted. The Pu concentration stays roughly constant throughout with slight enrichment in the center and depletion near the surface. Pu acts as a solvent for the mixture. The experiment was only run to 1.9% burnup, so the model is not at this time applicable to the high-burnup scenarios that the AMP code is intended to eventually model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Communication: Modeling of concentration dependent water diffusivity in ionic solutions: Role of intermolecular charge transfer

The translational diffusivity of water in solutions of alkali halide salts depends on the identity of ions, exhibiting dramatically different behavior even in solutions of similar salts of NaCl and KCl. The water diffusion coefficient decreases as the salt concentration increases in NaCl. Yet, in KCl solution, it slightly increases and remains above bulk value as salt concentration increases. Previous classical molecular dynamics simulations have failed to describe this important behavior even when polarizable models were used. Here, we show that inclusion of dynamical charge transfer among water molecules produces results in a quantitative agreement with experiments. Our results indicate that the concentration-dependent diffusivity reflects the importance of many-body effects among the water molecules in aqueous ionic solutions. Comparison with quantum mechanical calculations shows that a heterogeneous and extended distribution of charges on water molecules around the ions due to ion-water and also water-water charge transfer plays a very important role in controlling water diffusivity. Explicit inclusion of the charge transfer allows us to model accurately the difference in the concentration-dependent water diffusivity between Na{sup +} and K{sup +} ions in simulations, and it is likely to impact modeling of a wide range of systems for medical and technological applications.

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

Erosion and re-deposition of lithium and boron coatings under high-flux plasma bombardment

Lithium and boron coatings are applied to the walls of many tokamaks to enhance performance and protect the underlying substrates. Li and B-coated high-Z substrates are planned for use in NSTX-U and are a candidate plasma-facing component (PFC) for DEMO. However, previous measurements of Li evaporation and thermal sputtering on low-flux devices indicate that the Li temperature permitted on such devices may be unacceptably low. Thus it is crucial to characterize gross and net Li erosion rates under high-flux plasma bombardment. Additionally, no quantitative measurements have been performed of the erosion rate of a boron-coated PFC during plasma bombardment. A realistic model for the compositional evolution of a Li layer under D bombardment was developed that incorporates adsorption, implantation, and diffusion. A model was developed for temperature-dependent mixed-material Li-D erosion that includes evaporation, physical sputtering, chemical sputtering, preferential sputtering, and thermal sputtering. The re-deposition fraction of a Li coating intersecting a linear plasma column was predicted using atomic physics information and by solving the Li continuity equation. These models were tested in the Magnum-PSI linear plasma device at ion fluxes of 10^23-10^24 m^-2 s^-1 and Li surface temperatures less than 800 degrees C. Li erosion was measured during bombardment with a neon plasma that will not chemically react with Li and the results agreed well with the erosion model. Next the ratio of the total D fluence to the areal density of the Li coating was varied to quantify differences in Li erosion under D plasma bombardment as a function of the D concentration. The ratio of D/Li atoms was calculated using the results of MD simulations and good agreement is observed between measurements and the predictions of the mixed-material erosion model. Li coatings are observed to disappear from graphite much faster than from TZM Mo, indicating that fast Li diffusion into the bulk graphite substrate occurred, as predicted. Li re-deposition fractions very close to unity are observed in Magnum-PSI, as predicted by modeling. Finally, predictions of Li coating lifetimes in the NSTX-U divertor are calculated. The gross erosion rate of boron coatings was also measured for the first time in a high-flux plasma device.

liquid metals↗