No Arabic abstract
The changing thermal conductivity of an irradiated material is among the principal design considerations for any nuclear reactor, but at present few models are capable of predicting these changes starting from an arbitrary atomistic model. Here we present a simple model for computing the thermal diffusivity of tungsten, based on the conductivity of the perfect crystal and resistivity per Frenkel pair, and dividing a simulation into perfect and athermal regions statistically. This is applied to highly irradiated microstructures simulated with Molecular Dynamics. A comparison to experiment shows that simulations closely track observed thermal diffusivity over a range of doses from the dilute limit of a few Frenkel pairs to the high dose saturation limit at 3 displacements per atom (dpa).
Ion-irradiated FeCr alloys are useful for understanding and predicting neutron-damage in the structural steels of future nuclear reactors. Previous studies have largely focused on the structure of irradiation-induced defects, probed by transmission electron microscopy (TEM), as well as changes in mechanical properties. Across these studies, a wide range of irradiation conditions has been employed on samples with different processing histories, which complicates the analysis of the relationship between defect structures and material properties. Furthermore, key properties, such as irradiation-induced changes in thermal transport and lattice strain, are little explored. Here we present a systematic study of Fe3Cr, Fe5Cr and Fe10Cr binary alloys implanted with 20 MeV Fe$^{3+}$ ions to nominal doses of 0.01 dpa and 0.1 dpa at room temperature. Nanoindentation, transient grating spectroscopy (TGS) and X-ray micro-beam Laue diffraction were used to study the changes in hardness, thermal diffusivity and strain in the material as a function of damage and Cr content. Our results suggest that Cr leads to an increased retention of irradiation-induced defects, causing substantial changes in hardness and lattice strain. However, thermal diffusivity varies little with increasing damage and instead degrades significantly with increasing Cr content in the material. We find significant lattice strains even in samples exposed to a nominal displacement damage of 0.01 dpa. The defect density predicted from the lattice strain measurements is significantly higher than that observed in previous TEM studies, suggesting that TEM may not fully capture the irradiation-induced defect population.
We report a molecular dynamics simulation of melting of tungsten (W) nanoparticles. The modified embedded atom method (MEAM) interatomic potentials are used to describe the interaction between tungsten atoms. The melting temperature of unsupported tungsten nanoparticles of different sizes are found to decrease as the size of the particles decreases. The melting temperature obtained in the present study is approximately a decreasing function of inverse radius, in a good agreement with the predictions of thermodynamic models. We also observed that the melting of a W nanoparticle is preceded by the premelting of its outer skin at a temperature lower than its melting temperature.
The heat flux autocorrelation functions of carbon nanotubes (CNTs) with different radius and lengths is calculated using equilibrium molecular dynamics. The thermal conductance of CNTs is also calculated using the Green-Kubo relation from the linear response theory. By pointing out the ambiguity in the cross section definition of single wall CNTs, we use the thermal conductance instead of conductivity in calculations and discussions. We find that the thermal conductance of CNTs diverges with the CNT length. After the analysis of vibrational density of states, it can be concluded that more low frequency vibration modes exist in longer CNTs, and they effectively contribute to the divergence of thermal conductance.
Hydrogen isotopes are retained in materials for fusion power applications, changing both hydrogen embrittlement and tritium inventory as the microstructure undergoes irradiation damage. But modelling of highly damaged materials - exposed to over 0.1 displacements per atom (dpa) - where asymptotic saturation is observed, for example tungsten facing the plasma in a fusion tokamak reactor, is difficult because a highly damaged microstructure cannot be treated as weakly interacting isolated defect traps. In this paper we develop computational techniques to find the defect content in highly irradiated materials without adjustable parameters. First we show how to generate converged high dose (>1 dpa) microstructures using a combination of the creation-relaxation algorithm and molecular dynamics simulations of collision cascades. Then we make robust estimates of point defects and void regions with simple developments of the Wigner-Seitz decomposition of lattice sites. We use our estimates of the void surface area to predict the deuterium retention capacity of tungsten as a function of dose. This is then compared to 3He nuclear reaction analysis (NRA) measurements of tungsten samples self-irradiated at 290 K to different damage doses and exposed to deuterium plasma at low energy at 370 K. We show that our simulated microstructures give an excellent match to the experimental data, with both model and experiment showing 1.5-2.0 at.% deuterium retained in tungsten in the limit of high dose.
Measurements of thermal diffusivity in several insulators have been shown to reach a Planckian bound on thermal transport that can be thought of as the limit of validity of semiclassical phonon scattering. Beyond this regime, the heat transport must be understood in terms of incoherent motion of the atoms under strongly anharmonic interactions. In this work, we propose a model for heat transport in a strongly anharmonic system where the thermal diffusivity can be lower than the Planckian thermal diffusivity bound. Similar to the materials which exhibit thermal diffusivity close to this bound, our scenario involves complex unit cell with incoherent intra-cell dynamics. We derive a general formalism to compute thermal conductivity in such cases with anharmonic intra-cell dynamics coupled to nearly harmonic inter-cell coupling. Through direct numerical simulation of the non-linear unit cell motion, we explicitly show that our model allows sub-Planckian thermal diffusivity. We find that the propagator of the acoustic phonons becomes incoherent throughout most of the Brillouin zone in this limit. We expect these features to apply to more realistic models of complex insulators showing sub-Planckian thermal diffusivity, suggesting a multi-species generalization of the thermal diffusivity bound that is similar to the viscosity bound in fluids.