Comparative Study on Snowflake Dendrite Solidification Modeling Using a Phase-Field Model and by Cellular Automaton

Dendrite is among the most frequently observed structures during the solidification process. Different dendrite morphologies caused by environmental conditions can affect the physical properties of materials. The formation of snowflakes can generate various morphologies under different conditions, and is used in this work as an example. Simulation technologies provide insight into the correlation between a resulting morphology and its impact parameter, including the phase-field method (PF) and cellular automaton (CA). The PF method is derived from thermodynamic functions and kinetic equations, while the CA model is established by interaction rules between subsystems. It is difficult to solve the PF method due to the coupled differential equations, wherein the actual physical parameters are included. The CA model is conceptually simple and computationally efficient; however, the physical meaning of the parameters is absent. In this work, an example of snowflake formation is considered by PF with all the impact factors defined first, and then parameters in CA are searched by iterations to approximate the result, i.e., latent heat and the anisotropic coefficient in the PF method correspond to the initial distribution and the environmental effect in the CA model. In addition, the discrete time of each iteration in the CA model is identified according to the dendritic growth speed of these two models. A systematic identification process for the CA parameters' physical meaning is demonstrated by the comparison with the PF method, and an approximate simulation of the PF method can be obtained simply by the CA model. The combination of the PF method and the CA model can be used to investigate the influence of environmental factors on dendritic morphology.

Automated summary

This paper investigates the influence of environmental conditions on dendritic morphology and explores the application of the phase-field method and cellular automaton method. It is found that the result of the phase-field method can be approximated by adjusting the parameters of the cellular automaton model.