Program Library

- Provenance of code.
- Purpose of code.
- Specification.
- Description of subroutine's operation.
- References.
- Parameter descriptions.
- Error indicators.
- Accuracy estimate.
- Any additional information.
- Example of code
- Auxiliary subroutines required.
- Keywords.
- Download source code.
- Links.

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.

- Y. Saito, ISIJ International,
__38__(1998) 559-566. (Key paper entitled "Monte Carlo Simulation of Grain Growth in Three Dimensions".) - Y. Saito, Memoirs of the School of Science and Engineering, 59 (1995) 75-99. (This review paper can be obtained via an e-mail request to the author).
- Y. Saito and M. Enomoto, ISIJ International,
__32__(1992) 267-274. - Y. Saito and M. Enomoto, Computer Innovation of New Materials, 1991, 809-814, published by Elsevier Science, North-Holland.
- Y. Saito and H. Harada, Materials Science and Engineering,
__A223__(1997) 1-9. - Y. Saito, Materials Science and Engineering,
__A223__(1997) 10-16. - Y. Saito, Materials Science and Engineering,
__A223__(1997) 114-124. - Y. Saito, Materials Science and Engineering,
__A223__(1997) 125-133.

These are read from the file `Graingrowth2.montecarlo.data` which has a single row as follows :-

**I0**- integer- Number of cells in the X direction.
**J0**- integer- Number of cells in the Y direction.
**K0**- integer- Number of cells in the Y direction.
**IMCSTP**- integer- Time in Monte Carlo units.
**IDUMMY**- integer- Initial random number which is a prime number.
**IQ**- integer- Number of different grain orientations.
**NOGR**- integer- Number of grains.
**AJ**- real- Factor proportional to interfacial energy normalised by kT
**IANIS**- real- Represents anisotropy of interfacial energy.

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

There are four output files created,