120 likes | 651 Views
A Kinetic Monte Carlo Study Of Ordering in a Binary Alloy. Group 3: Tim Drews (ChE) Dan Finkenstadt (Physics) Xuemin Gu (MSE) CSE 373/MatSE 385/Physics 363 Final Project University of Illinois at Urbana-Champaign December 14, 2000. Code Introduction.
E N D
A Kinetic Monte Carlo Study Of Ordering in a Binary Alloy Group 3: Tim Drews (ChE) Dan Finkenstadt (Physics) Xuemin Gu (MSE) CSE 373/MatSE 385/Physics 363 Final Project University of Illinois at Urbana-Champaign December 14, 2000
Code Introduction • Main part of project: development of Kinetic Monte Carlo (KMC) code to simulate ordering in a binary superalloy • Equilibrium simulations that compute a vacancy pathway through the bulk alloy • Simulated alloy: similar to Fe-Al • Exhibits three phase behavior: • B2 ordered phase • A2 disordered phase • A2+B2 mixed phase • Compared the results to phase diagrams from the literature • Developed data analysis code to compute various order parameters • Developed many visualization techniques to determine the phase of the alloy
Move 1 2 3 4 5 6 7 8 0 1 where Nearest Neighbor Sites Kinetic Monte Carlo Method Jump Frequency and Potential Function
Phase Diagrams Theoretical Simulated with Grand Canonical MC Simulations Reproduced from A Monte-Carlo Study of B2 Ordering and Precipitation Via Vacancy Mechanism in B.C.C. Lattices. Athenes, M., Bellon, P., Martin, G., and Haider, F. Acta Metallurgica. Vol. 44, No. 12, pp. 4739-4748, 1996.
BCC Lattice Two Simple Cubic Lattices Bulk BCC Lattice • sublattice (arbitrary) sublattice (arbitrary)
2 plot sublattice sublattice 64 Cell Cubic Lattice - Disordered A2 Phase T = 1000 K u1 = -0.04 1*107 MC Steps cB = 0.25
32 Cell Cubic Lattice - High and Low Order B2 Phase T = 700 K MC steps = 1*107 u1 = -0.04 cB = 0.25 2 plot sublattice sublattice T = 700 K MC steps = 1*107 u1 = -0.04 cB = 0.50
2 plot sublattice sublattice 64 Cell Cubic Lattice - High Order B2 Phase T = 700 K u1 = -0.04 1*107 MC Steps cB = 0.45
16, 32, and 64 Cell Cubic Lattice - High Order B2 Phase T = 700 K u1 = -0.04 cB = 0.45 1*106 MC steps sublattice sublattice 1*107 MC steps sublattice sublattice sublattice sublattice
32 Cell Cubic Lattice - Mixed A2+B2 Phase T = 200 K u1 = -0.04 1*107 MC Steps cB = 0.25 2 plot sublattice sublattice
Conclusions • KMC code can generate expected equilibrium phase behavior, for a given temperature and concentration, provided enough MC steps are taken for the given simulation cell size • Observed dynamic growth and ageing • Did not distinguish phase transitions • Could do this by computing total free energy and looking for kinks in the free energy diagram • Could also look for critical slowing down or finite-size scaling phenomena • Did not vary the activation energy • Done by Athenes, et al., and this energy was found to have a significant effect • Can be readily done with this code