Three-dimensional normal grain growth was appropriately simulated using a Potts model Monte Carlo algorithm. The quasi-stationary grain size distribution obtained from simulation agreed well with the experimental result of pure iron. The Weibull function with a parameter β=2.77 and the Yu-Liu function with a parameter ν =2.71 fit the quasi-stationary grain size distribution well. The grain volume distribution is a function that decreased exponentially with increasing grain volume. The distribution of boundary area of grains has a peak at S/ Three-dimensional normal grain growth was appropriate simulated using a Potts model Monte Carlo algorithm. The quasi-stationary grain size distribution obtained from simulation agreed well with the experimental result of pure iron. The Weibull function with a parameter β = 2.77 and the Yu- Liu function with a parameter ν = 2.71 fit the quasi-stationary grain size distribution well. The grain volume distribution is a function that decreased exponentially with increasing grain volume. The distribution of boundary area of ​​grains has a peak at S /