In the paper, we focus on atom diffusion behavior in Ni-based superalloys, which have important applications in the aero-industry. Specifically, the expressions of the key physical parameter – transition rate (jump rate) in the diffusion can be given from the diffusion theory in solids and the kinetic Monte Carlo (KMC) method, respectively. The transition rate controls the diffusion process and is directly related to the energy of vacancy formation and the energy of migration of atom from density functional theory (DFT). Moreover, from the KMC calculations, the diffusion coefficients for Ni and Al atoms in the γ phase (Ni matrix) and the γʹ phase (intermetallic compound Ni3Al) of the superalloy have been obtained. We propose a strategy of time stepping to deal with the multi-time scale issues. In addition, the influence of temperature and Al concentration on diffusion in dilute alloys is also reported.

We present results of simulations of the energetics and dynamics involved in the realization of the NV (nitrogen-vacancy) center in diamond. We use the self-consistent charge-density functional tight-binding approximation and show that when the nitrogen resides on a single substitutional site, it fails to attract a vacancy, hence no NV center can be formed. However, if it occupies a split interstitial site and two vacancies reside on the second or third neighbor sites, an NV center will form following annealing at temperatures as low as 300°C and 650°C, respectively. These results provide guidelines to experimentalists on how to increase the efficiency of NV formation in diamond.

We employ the lattice Boltzmann method and random walk particle tracking to simulate the time evolution of hydrodynamic dispersion in bulk, random, monodisperse, hard-sphere packings with bed porosities (interparticle void volume fractions) between the random-close and the random-loose packing limit. Using Jodrey-Tory and Monte Carlo-based algorithms and a systematic variation of the packing protocols we generate a portfolio of packings, whose microstructures differ in their degree of heterogeneity (DoH). Because the DoH quantifies the heterogeneity of the void space distribution in a packing, the asymptotic longitudinal dispersion coefficient calculated for the packings increases with the packings’ DoH. We investigate the influence of packing length (up to 150 dp, where dp is the sphere diameter) and grid resolution (up to 90 nodes per dp) on the simulated hydrodynamic dispersion coefficient, and demonstrate that the chosen packing dimensions of 10 dpx 10 dpx 70 dp and the employed grid resolution of 60 nodes per dp are sufficient to observe asymptotic behavior of the dispersion coefficient and to minimize finite size effects. Asymptotic values of the dispersion coefficients calculated for the generated packings are compared with simulated as well as experimental data from the literature and yield good to excellent agreement.