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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Criticality Calculations with MCNP6 - Practical Lectures

These slides are used to teach MCNP (Monte Carlo N-Particle) usage to nuclear criticality safety analysts. The following are the lecture topics: course information, introduction, MCNP basics, criticality calculations, advanced geometry, tallies, adjoint-weighted tallies and sensitivities, physics and nuclear data, parameter studies, NCS validation I, NCS validation II, NCS validation III, case study 1 - solution tanks, case study 2 - fuel vault, case study 3 - B&W core, case study 4 - simple TRIGA, case study 5 - fissile mat. vault, criticality accident alarm systems. After completion of this course, you should be able to: Develop an input model for MCNP; Describe how cross section data impact Monte Carlo and deterministic codes; Describe the importance of validation of computer codes and how it is accomplished; Describe the methodology supporting Monte Carlo codes and deterministic codes; Describe pitfalls of Monte Carlo calculations; Discuss the strengths and weaknesses of Monte Carlo and Discrete Ordinants codes; The diffusion theory model is not strictly valid for treating fissile systems in which neutron absorption, voids, and/or material boundaries are present. In the context of these limitations, identify a fissile system for which a diffusion theory solution would be adequate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Determining acceptable limits of fast-electron preheat in direct-drive-ignition–scale target designs

In direct-drive–ignition designs, preheat by fast electrons created by the two-plasmon–decay or stimulated Raman instabilities can increase the adiabat in the fuel layer and reduce compression and neutron yield. Since eliminating the preheat entirely is a major challenge, it is necessary to understand the levels of preheat that preclude ignition in a direct-drive target. Two 1-D ignition-scale target designs serve as the basis for examining the effects of synthetically increasing the levels of fast electrons using the 1-D radiation–hydrodynamic code LILAC , which includes two models of fast-electron transport. The first is an ignition design adapted from a 2-D polar-direct-drive (PDD) design for the National Ignition Facility (NIF). The second is derived from the first but uses a laser pulse that generates stronger shocks and a higher fuel adiabat. This more stable “alpha-burner” design approaches ignition and achieves yield multiplication as a result of alpha heating. The designs are then re optimized to recover performance. The igniting design, when fast-electron transport was modeled with diffusion, was found to tolerate 50% more fast-electron preheat of the cold (sub-50 eV) DT ice layer when the laser pulse was optimized using the optimizer Telios . When a straight-line fast-electron transport model was used, the effects of optimization were negligible. For the alpha-burner design, an increase of over a factor of at least 3 in the tolerable level of fast-electron preheat was obtained for both transport models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Methane production and emissions in trees and forests

Forest ecosystem methane (CH 4 ) research has focused on soils, but trees are also important sources and sinks in forest CH 4 budgets. Living and dead trees transport and emit CH 4 produced in soils; living trees and dead wood emit CH 4 produced inside trees by microorganisms; and trees produce CH 4 through an abiotic photochemical process. Here, we review the state of the science on the production, consumption, transport, and emission of CH 4 by living and dead trees, and the spatial and temporal dynamics of these processes across hydrologic gradients inclusive of wetland and upland ecosystems. Emerging research demonstrates that tree CH 4 emissions can significantly increase the source strength of wetland forests, and modestly decrease the sink strength of upland forests. Scaling from stem or leaf measurements to trees or forests is limited by knowledge of the mechanisms by which trees transport soil-produced CH 4 , microbial processes produce and oxidize CH 4 inside trees, a lack of mechanistic models, the diffuse nature of forest CH 4 fluxes, complex overlap between sources and sinks, and extreme variation across individuals. Understanding the complex processes that regulate CH 4 source–sink dynamics in trees and forests requires cross-disciplinary research and new conceptual models that transcend the traditional binary classification of wetland vs upland forest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Effects of Moisture on Diffusion in Unmodified Wood Cell Walls: A Phenomenological Polymer Science Approach

Despite the importance of cell wall diffusion to nearly all aspects of wood utilization, diffusion mechanisms and the detailed effects of moisture remain poorly understood. In this perspective, we introduce and employ approaches established in polymer science to develop a phenomenological framework for understanding the effects of moisture on diffusion in unmodified wood cell walls. The premise for applying this polymer-science-based approach to wood is that wood polymers (cellulose, hemicelluloses, and lignin) behave like typical solid polymers. Therefore, the movement of chemicals through wood cell walls is a diffusion process through a solid polymer, which is in contrast to previous assertions that transport of some chemicals occurs via aqueous pathways in the cell wall layers. Diffusion in polymers depends on the interrelations between free volume in the polymer matrix, molecular motions of the polymer, diffusant dimensions, and solubility of the diffusant in the polymer matrix. Because diffusion strongly depends on whether a polymer is in a rigid glassy state or soft rubbery state, it is important to understand glass transitions in the amorphous wood polymers. Through a review and analysis of available literature, we conclude that in wood both lignin and the amorphous polysaccharides very likely have glass transitions. After developing and presenting this polymer-science-based perspective of diffusion through unmodified wood cell walls, suggested directions for future research are discussed. A key consideration is that a large difference between diffusion through wood polymers and typical polymers is the high swelling pressures that can develop in unmodified wood cell walls. This pressure likely arises from the hierarchical structure of wood and should be taken into consideration in the development of predictive models for diffusion in unmodified wood cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Model-Independent Characterization of Charge Diffusion in Thick Fully Depleted CCDs

We present a new method to measure charge diffusion in charge-coupled devices (CCDs). The method is based on a statistical characterization of the shapes of charge clouds produced by low-energy X-rays using known properties of the two-dimensional Gaussian point-spread function (PSF). The algorithm produces reliable upper and lower bounds on the size of the PSF for photons converting near the entrance window of a device. It is optimally suited to the case of thick back-illuminated CCDs where the X-ray absorption length is smaller than the silicon thickness and the diffusion scale is of the same order as the pixel size. The only assumptions are that the charge cloud width is a monotonically increasing function of distance from the conversion point to the buried channel, and that the conversion probability is peaked at the surface. Otherwise, no physical models of carrier transport or knowledge of the electric field profile in the CCD are needed. In suboptimal conditions, the upper bound increases and the lower bound is unaffected, so confidence in the correctness of results is retained. The new method has been benchmarked against Monte Carlo simulations and tested on X-ray images measured on thick high-resistivity prototype CCDs for the Large Synoptic Survey Telescope. In Monte Carlo simulations of noiseless images having the optimal diffusion scale, the upper bound approximated the true PSF within 5%, increasing to 10% in simulations with added noise. Even with severely undersampled or truncated PSFs, the method brackets the true value to within 25%. Our method is accurate and computationally efficient and offers a fast and simple experimental setup.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Modeling the Ionic Strength Effect on Diffusion in Clay. The DR-A Experiment at Mont Terri

Solute diffusion in compacted clays depends on ionic strength through its control on the thickness of the electrical double layer (EDL) on the charged clay surfaces. In the DR-A field experiment (Mont Terri, Switzerland), synthetic porewater was circulated through a borehole for 189 days, leading to the out-diffusion of a variety of tracers into the Opalinus Clay. The borehole solution was then replaced with a higher-salinity solution for an additional 540 days, leading to the diffusion of Cs + , Ca 2+ , Mg 2+ , and Sr 2+ back into the borehole and to an increase in the out-diffusion of anions (I - , Br - ) and 3 H. The experimental results were interpreted using the CrunchClay code, which includes a mean electrostatic potential model for the EDL. The EDL corresponds to a second continuum in addition to bulk electrically neutral porewater. Species-specific diffusion (Nernst-Planck equation) occurs through both domains. A 1D radial model considered a single pore diffusion coefficient (Dp = 10-9 m2/s) for cations and 3H in the bulk porosity, and a smaller D p (3 × 10 -10 m 2 /s) for anions. Dp values in the EDL were smaller (10 -11 m 2 /s), except for Cs + and K + (5 × 10 -10 and 2 × 10 -10 m 2 /s, respectively). The model reproduced well the experimental results and showed the capability to consider temporal changes in geochemical conditions affecting the transport and retention of potentially important radionuclide contaminants (e.g., 137 Cs + , 90 Sr 2+ , 129 I - ) in underground geological nuclear waste repositories. Finally, coupled multicomponent diffusion together with the electrostatic properties of the charged surfaces are essential in the development of predictive models for ion transport in clays.

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

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2015

Impact of multi-component diffusion in turbulent combustion using direct numerical simulations

This study presents the results of DNS of a partially premixed turbulent syngas/air flame at atmospheric pressure. The objective was to assess the importance and possible effects of molecular transport on flame behavior and structure. To this purpose DNS were performed at with two proprietary DNS codes and with three different molecular diffusion transport models: fully multi-component, mixture averaged, and imposing the Lewis number of all species to be unity.

10 SYNTHETIC FUELS↗
DOE OSTI2014

Interactions of mobile helium clusters with surfaces and grain boundaries of plasma-exposed tungsten

We report results of atomistic computations for the interactions of small mobile helium clusters (He{sub n}) with free surfaces and grain boundaries (GBs) in tungsten toward development of continuum drift-diffusion-reaction models for the dynamics of mobile helium clusters in plasma-exposed tungsten. Molecular-statics (MS) simulations based on reliable many-body interatomic potentials are carried out for He{sub n} (1 ≤ n ≤ 7) clusters near sinks to obtain the potential energy profiles of the He{sub n} clusters as a function of the clusters' center-of-mass distance from a sink. Sinks investigated include surfaces, GBs, and regions in the vicinity of junctions where GBs intersect free surfaces. Elastic interaction potentials based on elastic inclusion theory provide an excellent description of the MS results for the cluster-sink interactions. The key parameter in the elastic models is the sink segregation strength, which is found to increase with increasing cluster size. Such cluster-sink interactions are responsible for the migration of small helium clusters by drift and for helium segregation on surfaces and grain boundaries in tungsten. Such helium segregation on sinks is observed in large-scale molecular-dynamics simulations of helium aggregation in model polycrystalline tungsten at 933 K upon helium implantation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Effect of burst and recombination models for Monte Carlo transport of interacting carriers in a-Se x-ray detectors on Swank noise

Purpose: The authors describe the modification to a previously developed Monte Carlo model of semiconductor direct x-ray detector required for studying the effect of burst and recombination algorithms on detector performance. This work provides insight into the effect of different charge generation models for a-Se detectors on Swank noise and recombination fraction. Methods: The proposed burst and recombination models are implemented in the Monte Carlo simulation package, ARTEMIS, developed byFang et al. [“Spatiotemporal Monte Carlo transport methods in x-ray semiconductor detectors: Application to pulse-height spectroscopy in a-Se,” Med. Phys. 39(1), 308–319 (2012)]. The burst model generates a cloud of electron-hole pairs based on electron velocity, energy deposition, and material parameters distributed within a spherical uniform volume (SUV) or on a spherical surface area (SSA). A simple first-hit (FH) and a more detailed but computationally expensive nearest-neighbor (NN) recombination algorithms are also described and compared. Results: Simulated recombination fractions for a single electron-hole pair show good agreement with Onsager model for a wide range of electric field, thermalization distance, and temperature. The recombination fraction and Swank noise exhibit a dependence on the burst model for generation of many electron-hole pairs from a single x ray. The Swank noise decreased for the SSA compared to the SUV model at 4 V/μm, while the recombination fraction decreased for SSA compared to the SUV model at 30 V/μm. The NN and FH recombination results were comparable. Conclusions: Results obtained with the ARTEMIS Monte Carlo transport model incorporating drift and diffusion are validated with the Onsager model for a single electron-hole pair as a function of electric field, thermalization distance, and temperature. For x-ray interactions, the authors demonstrate that the choice of burst model can affect the simulation results for the generation of many electron-hole pairs. The SSA model is more sensitive to the effect of electric field compared to the SUV model and that the NN and FH recombination algorithms did not significantly affect simulation results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Charge decay in the spatial afterglow of plasmas and its impact on diffusion regimes

Abstract The spatial afterglow is a region at the boundary of a non-equilibrium plasma where charged species relax into ambient equilibrium. In many applications, the spatial afterglow is the part of the plasma that interacts with surfaces, such as suspended particles or a material substrate. However, compared to the bulk plasma, there has been little effort devoted to studying the properties of the spatial afterglow, and a fundamental analysis has not yet been developed. Here, we apply double Langmuir probe measurements and develop an advection-diffusion-recombination model to provide a detailed description of charged species in the spatial afterglow over a wide range of pressures, temperatures, plasma dimensions, and flow rates. We find that the density of charged species in the spatial afterglow decays by orders of magnitude, which leads to a transition from ambipolar to free diffusion. These insights can be used to explain or predict experimental observations of phenomena, such as the charging of dust grains and the dose of charged species to a biomaterial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Efficient Reformulation of Solid Phase Diffusion in Electrochemical-Mechanical Coupled Models for Lithium-Ion Batteries: Effect of Intercalation Induced Stresses

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms with an appropriate discretization or approximation for the solid phase diffusion within the electrode particle. One of the major difficulties in simulating Li-ion battery models is the need for simulating solid-phase diffusion in the second radial dimension r within the particle. It increases the complexity of the model as well as the computation time/cost to a great extent. This is Particularly true for the inclusion of pressure induced diffusion inside particles experiencing volume change. A computationally efficient representation for solid-phase diffusion is discussed in this paper. The operating condition has a significant effect on the validity, accuracy, and efficiency of various approximations for the solid-phase transport governed by pressure induced diffusion. This paper introduces efficient methods for solid phase reformulation - (1) parabolic profile approach and (2) a mixed order finite difference method for approximating/representing solid-phase concentration variations within the active materials of porous electrodes for macroscopic models for lithium-ion batteries. (C) 2013 The Electrochemical Society. All rights reserved.

De, S↗
DOE OSTI2021

Cherenkov Telescope Array sensitivity to the putative millisecond pulsar population responsible for the Galactic Centre excess

ABSTRACT The leading explanation of the Fermi Galactic Centre γ-ray excess is the extended emission from an unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt γ-rays, also inject large quantities of electrons/positrons (e±) into the interstellar medium. These e± could potentially inverse-Compton (IC) scatter ambient photons into γ-rays that fall within the sensitivity range of the upcoming Cherenkov Telescope Array (CTA). In this article, we examine the detection potential of CTA to this signature by making a realistic estimation of the systematic uncertainties on the Galactic diffuse emission model at TeV-scale γ-ray energies. We forecast that, in the event that e± injection spectra are harder than E−2, CTA has the potential to robustly discover the IC signature of a putative Galactic bulge MSP population sufficient to explain the Galactic Centre excess for e± injection efficiencies in the range of ≈2.9–74.1 per cent, or higher, depending on the level of mismodelling of the Galactic diffuse emission components. On the other hand, for spectra softer than E−2.5, a reliable CTA detection would require an unphysically large e± injection efficiency of ${\gtrsim} 158{{\ \rm per\ cent}}$. However, even this pessimistic conclusion may be avoided in the plausible event that MSP observational and/or modelling uncertainties can be reduced. We further find that, in the event that an IC signal were detected, CTA can successfully discriminate between an MSP and a dark matter origin for the radiating e±.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Diffusion-driven growth of nanowires by low-temperature molecular beam epitaxy

With ZnTe as an example, we use two different methods to unravel the characteristics of the growth of nanowires (NWs) by gold-catalyzed molecular beam epitaxy at low temperature. In the first approach, CdTe insertions have been used as markers, and the nanowires have been characterized by scanning transmission electron microscopy, including geometrical phase analysis and energy dispersive electron spectrometry; the second approach uses scanning electron microscopy and the statistics of the relationship between the length of the tapered nanowires and their base diameter. Axial and radial growth are quantified using a diffusion-limited model adapted to the growth conditions; analytical expressions describe well the relationship between the NW length and the total molecular flux (taking into account the orientation of the effusion cells), and the catalyst-nanowire contact area. A long incubation time is observed. This analysis allows us to assess the evolution of the diffusion lengths on the substrate and along the nanowire sidewalls, as a function of temperature and deviation from stoichiometric flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Diffusive mass transport in agglomerated glassy fallout from a near-surface nuclear test

Aerodynamically-shaped glassy fallout is formed when vapor phase constituents from the nuclear device are incorporated into molten carriers (i.e. fallout precursor materials derived from soil or other near-field environmental debris). The effects of speciation and diffusive transport of condensing constituents are not well defined in models of fallout formation. Previously we reported observations of diffuse micrometer scale layers enriched in Na, Fe, Ca, and 235 U, and depleted in Al and Ti, at the interfaces of agglomerated fallout objects. Here in this paper, we derive the timescales of uranium mass transport in such fallout as it cools from 2500 K to 1500 K by applying a 1-dimensional planar diffusion model to the observed 235 U/ 30 Si variation at the interfaces. By modeling the thermal transport between the fireball and the carrier materials, the time of mass transport is calculated to be <0.6 s, <1 s, <2 s, and <3.5 s for fireball yields of 0.1 kt, 1 kt, 10 kt, and 100 kt respectively. Based on the calculated times of mass transport, a maximum temperature of deposition of uranium onto the carrier material of ~2200 K is inferred (1σ uncertainty of ~200 K). We also determine that the occurrence of micrometer scale layers of material enriched in relatively volatile Na-species as well as more refractory Ca-species provides evidence for an oxygen-rich fireball based on the vapor pressure of the two species under oxidizing conditions. These results represent the first application of diffusion-based modeling to derive material transport, thermal environments, and oxidation-speciation in near-surface nuclear detonation environments.

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

Spatially dependent cluster dynamics model of He plasma surface interaction in tungsten for fusion relevant conditions

In fusion reactors, plasma facing components (PFC) and, in particular, the divertor will be irradiated with high fluxes of low-energy (similar to 100 eV) helium and hydrogen ions. Tungsten is one of the leading candidate divertor materials for ITER and DEMO fusion reactors. However, the behaviour of tungsten under high dose, coupled helium/hydrogen exposure remains to be fully understood. The PFC response and performance changes are intimately related to microstructural changes, such as the formation of point defect clusters, helium and hydrogen bubbles or dislocation loops. Computational materials' modelling results are described here that investigate the mechanisms controlling microstructural evolution in tungsten. The aim of this study is to understand and predict sub-surface helium bubble growth under high flux helium ion implantation (similar to 10 22 m -2 s -1 ) at high temperatures (>1000 K). Here, we report results from a spatially dependent cluster dynamics model based on reaction-diffusion rate theory to describe the evolution of the microstructure under these conditions. The key input parameters to the model (diffusion coefficients, migration and binding energies, initial defect production) are determined from a combination of atomistic modelling and available experimental data. The results are in good agreement with results of an analytical model that is presented in a separate paper. In particular, it is found that the sub-surface evolution with respect to bubble size and concentration of the helium bubbles strongly depends on the flux and temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗