skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Analysis of simulation methodology for calculation of the heat of transport for vacancy thermodiffusion

Abstract

Computation of the heat of transport Q{sub a}{sup *} in monatomic crystalline solids is investigated using the methodology first developed by Gillan [J. Phys. C: Solid State Phys. 11, 4469 (1978)] and further developed by Grout and coworkers [Philos. Mag. Lett. 74, 217 (1996)], referred to as the Grout-Gillan method. In the case of pair potentials, the hopping of a vacancy results in a heat wave that persists for up to 10 ps, consistent with previous studies. This leads to generally positive values for Q{sub a}{sup *} which can be quite large and are strongly dependent on the specific details of the pair potential. By contrast, when the interactions are described using the embedded atom model, there is no evidence of a heat wave, and Q{sub a}{sup *} is found to be negative. This demonstrates that the dynamics of vacancy hopping depends strongly on the details of the empirical potential. However, the results obtained here are in strong disagreement with experiment. Arguments are presented which demonstrate that there is a fundamental error made in the Grout-Gillan method due to the fact that the ensemble of states only includes successful atom hops and hence does not represent an equilibrium ensemble. This placesmore » the interpretation of the quantity computed in the Grout-Gillan method as the heat of transport in doubt. It is demonstrated that trajectories which do not yield hopping events are nevertheless relevant to computation of the heat of transport Q{sub a}{sup *}.« less

Authors:
;  [1]
  1. Advanced Material Processing and Analysis Center and Department of Physics, University of Central Florida, 4000 Central Florida Blvd., Orlando, Florida 32816 (United States)
Publication Date:
OSTI Identifier:
22308966
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Applied Physics; Journal Volume: 116; Journal Issue: 2; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; ATOMS; CALCULATION METHODS; EQUILIBRIUM; HEAT; HEAT TRANSFER; INTERACTIONS; SIMULATION; SOLIDS; THERMAL DIFFUSION; TRAJECTORIES; VACANCIES

Citation Formats

Tucker, William C., and Schelling, Patrick K., E-mail: patrick.schelling@ucf.edu. Analysis of simulation methodology for calculation of the heat of transport for vacancy thermodiffusion. United States: N. p., 2014. Web. doi:10.1063/1.4887121.
Tucker, William C., & Schelling, Patrick K., E-mail: patrick.schelling@ucf.edu. Analysis of simulation methodology for calculation of the heat of transport for vacancy thermodiffusion. United States. doi:10.1063/1.4887121.
Tucker, William C., and Schelling, Patrick K., E-mail: patrick.schelling@ucf.edu. Mon . "Analysis of simulation methodology for calculation of the heat of transport for vacancy thermodiffusion". United States. doi:10.1063/1.4887121.
@article{osti_22308966,
title = {Analysis of simulation methodology for calculation of the heat of transport for vacancy thermodiffusion},
author = {Tucker, William C. and Schelling, Patrick K., E-mail: patrick.schelling@ucf.edu},
abstractNote = {Computation of the heat of transport Q{sub a}{sup *} in monatomic crystalline solids is investigated using the methodology first developed by Gillan [J. Phys. C: Solid State Phys. 11, 4469 (1978)] and further developed by Grout and coworkers [Philos. Mag. Lett. 74, 217 (1996)], referred to as the Grout-Gillan method. In the case of pair potentials, the hopping of a vacancy results in a heat wave that persists for up to 10 ps, consistent with previous studies. This leads to generally positive values for Q{sub a}{sup *} which can be quite large and are strongly dependent on the specific details of the pair potential. By contrast, when the interactions are described using the embedded atom model, there is no evidence of a heat wave, and Q{sub a}{sup *} is found to be negative. This demonstrates that the dynamics of vacancy hopping depends strongly on the details of the empirical potential. However, the results obtained here are in strong disagreement with experiment. Arguments are presented which demonstrate that there is a fundamental error made in the Grout-Gillan method due to the fact that the ensemble of states only includes successful atom hops and hence does not represent an equilibrium ensemble. This places the interpretation of the quantity computed in the Grout-Gillan method as the heat of transport in doubt. It is demonstrated that trajectories which do not yield hopping events are nevertheless relevant to computation of the heat of transport Q{sub a}{sup *}.},
doi = {10.1063/1.4887121},
journal = {Journal of Applied Physics},
number = 2,
volume = 116,
place = {United States},
year = {Mon Jul 14 00:00:00 EDT 2014},
month = {Mon Jul 14 00:00:00 EDT 2014}
}
  • We apply the recently developed constrained-dynamics method to elucidate the thermodiffusion of vacancies in a single-component material. The derivation and assumptions used in the method are clearly explained. Next, the method is applied to compute the reduced heat of transport Q{sub v}{sup *}−h{sub fv} for vacancies in a single-component material. Results from simulations using three different Morse potentials, with one providing an approximate description of Au, and an embedded-atom model potential for Ni are presented. It is found that the reduced heat of transport Q{sub v}{sup *}−h{sub fv} may take either positive or negative values depending on the potential parametersmore » and exhibits some dependence on temperature. It is also found that Q{sub v}{sup *}−h{sub fv} may be correlated with the activation entropy. The results are discussed in comparison with experimental and previous simulation results.« less
  • Thermal diffusion causes a gradient in water vapor pressure in air-vapor mixtures in a temperature gradient. Estimates of the thermal diffusion constant for air-vapor mixtures permit the magnitude of the thermal diffusion effect to be estimated using conventional Onsager relations. When expected thermal gradients in an air gap adjacent to a high level waste container in a hypothetical geological repository with optimized geometries are considered, the thermal convection of the gas phase (thermogravitation) can result in a significant vertical separation of water vapor and air. For geometries including a vertical container surface and an air gap on the order ofmore » 2 cm wide and 600 cm tall, the relative humidity at the top of the convective system is up to 30% higher than when the gas mixture between a waste container and nearby backfill is assumed to be homogeneous. Such a difference may have important consequences, in particular, for the rate of corrosion of the waste container. However, high thermogravitational separation in the vicinity of the waste container is restricted to particular geometries, and waste repository can be designed to minimize the effect. Geometries that preserve an air gap of more than 5 cm or geometries prohibiting free convection are expected to eliminate thermogravitational separation of water vapor and air.« less
  • A fundamental understanding of spatial and temporal thermal distributions is crucial for predicting solidification and solid-state microstructural development in parts made by additive manufacturing. While sophisticated numerical techniques that are based on finite element or finite volume methods are useful for gaining insight into these phenomena at the length scale of the melt pool (100 - 500 µm), they are ill-suited for predicting engineering trends over full part cross-sections (> 10 x 10 cm) or many layers over long process times (> many days) due to the necessity of fully resolving the heat source characteristics. On the other hand, itmore » is extremely difficult to resolve the highly dynamic nature of the process using purely in-situ characterization techniques. This article proposes a pragmatic alternative based on a semi-analytical approach to predicting the transient heat conduction during powder bed metal additive manufacturing process. The model calculations were theoretically verified for selective laser melting of AlSi10Mg and electron beam melting of IN718 powders for simple cross-sectional geometries and the transient results are compared to steady state predictions from the Rosenthal equation. It is shown that the transient effects of the scan strategy create significant variations in the melt pool geometry and solid-liquid interface velocity, especially as the thermal diffusivity of the material decreases and the pre-heat of the process increases. With positive verification of the strategy, the model was then experimentally validated to simulate two point-melt scan strategies during electron beam melting of IN718, one intended to produce a columnar and one an equiaxed grain structure. Lastly, through comparison of the solidification conditions (i.e. transient and spatial variations of thermal gradient and liquid-solid interface velocity) predicted by the model to phenomenological CET theory, the model accurately predicted the experimental grain structures.« less
  • Abstract not provided.
  • Ferrofluids consist of magnetic nanoparticles dispersed in a carrier liquid. Their strong thermodiffusive behaviour, characterised by the Soret coefficient, coupled with the dependency of the fluid's parameters on magnetic fields is dealt with in this work. It is known from former experimental investigations on the one hand that the Soret coefficient itself is magnetic field dependent and on the other hand that the accuracy of the coefficient's experimental determination highly depends on the volume concentration of the fluid. The thermally driven separation of particles and carrier liquid is carried out with a concentrated ferrofluid (φ = 0.087) in a horizontalmore » thermodiffusion cell and is compared to equally detected former measurement data. The temperature gradient (1 K/mm) is applied perpendicular to the separation layer. The magnetic field is either applied parallel or perpendicular to the temperature difference. For three different magnetic field strengths (40 kA/m, 100 kA/m, 320 kA/m) the diffusive separation is detected. It reveals a sign change of the Soret coefficient with rising field strength for both field directions which stands for a change in the direction of motion of the particles. This behaviour contradicts former experimental results with a dilute magnetic fluid, in which a change in the coefficient's sign could only be detected for the parallel setup. An anisotropic behaviour in the current data is measured referring to the intensity of the separation being more intense in the perpendicular position of the magnetic field: S{sub T‖} = −0.152 K{sup −1} and S{sub T⊥} = −0.257 K{sup −1} at H = 320 kA/m. The ferrofluiddynamics-theory (FFD-theory) describes the thermodiffusive processes thermodynamically and a numerical simulation of the fluid's separation depending on the two transport parameters ξ{sub ‖} and ξ{sub ⊥} used within the FFD-theory can be implemented. In the case of a parallel aligned magnetic field, the parameter can be determined to ξ{sub ‖} = (2.8; 9.1; 11.2) × 10{sup −11} · D{sub ‖} kg/(A{sup 2}m) for the different field strengths and in dependence on the magnetic diffusion coefficient D{sub ‖}. An adequate fit in the perpendicular case is not possible, by ξ{sub ⊥} = 1 × 10{sup −17} kg/(Am{sup 2}) a rather good agreement between numerical and experimental data can be found for a field strength of 40 kA/m, a change in the coefficient's sign in the perpendicular setup is not numerically determinable via this theory. The FFD-theory is only partly applicable to calculate the concentration profile in concentrated magnetic fluids established due to a temperature gradient and magnetic field applied.« less