Professor Yoshiyuki Saito
Waseda University
Department of Materials Science and Engineering
3-4-1 Okubo, Shinjuku-ku, Tokyo, 169 Japan
yoshi@dice.cache.waseda.ac.jp
This research was supported by the Nippon Steel Corporation.
Simulation of grain growth kinetics using a Monte Carlo method in three dimensions. This is a new algorithm which has the advantage that it prevents grains of the same orientation from coalescing.
Language: | FORTRAN |
Product form: | Source code |
The program reads inputs from a file, in this case called
Graingrowth2.montecarlo.dat which contains the following:
64 64 64 2000 1749573 32 32768 3.00 1
The first, second and third numbers specify the system size which is the number of cells in the three-dimensional Monte Carlo simulation. The fourth number is the time in Monte Carlo units. The fifth number is the initial random number which should be a prime number. The sixth number sets the total number of orientations of the grains in the simulation. The seventh number is the initial number of grains. The eight number is J/kT where J is proportional to the interfacial energy, k is the Boltzmann constant and T is the absolute temperature. The ninth number is the anisotropy of the grain boundary energy (a value of unity represents isotropy, and values greater than unity represent an orientation dependence of interface energy).
The Monte Carlo simulation method is now widely applied to materials science and engineering. This program deals with grain growth kinetics in three-dimensions, with the simulation of interface motion based on the Potts model. As an initial microstructure, an orientation between 1 to Q is assigned to each lattice site at random, where Q is the total number of orientations. The evolution of the grain structure then occurs during the Monte Carlo iterations. It is possible to obtain both the mean grain size and the grain size distribution.
These are read from the file Graingrowth2.montecarlo.data which has a single row as follows :-
Output is to three ASCII files, fort.15, fort.16 and fort.17. See program results (below) for format.
None.
See [2], which has a comparison between calculations and experimental observations.
See full program
32 32 32 100 1749573 32 32768 3.00 1