some paraview post-processing scripts to highlight different orientations in colour
- what
- Minimising surface energy by surface diffusion (surface divergence of the mean curvature) [Mullins1957]. This is not equivalent to evolving shape by surface tension (mean curvature). The regularisation method used is well-documented in [Torabi2009] and the gamma construction method is due to [Siem2005].
- how
- The surfaces of different families are identified and the surface energy associated with the surface are specified.
- the alpha values are normalised according to:
- gamma-plot is obtained by fine tuning of the parameters alpha and w to ensure no overlapping between the energy minima. The units of xyz are dimensionless as gamma is normalised by gamma_0
- similarly xi vector is obtained by differentiating
$hat{gamma}=rgamma(vec{n})$ in$r$ space where$vec{r}$ is the position vector
- missing orientations are found by means of the inverse-gamma-plot [Herring1951] and the well-known criteria [Sekerka2005]
- A-TiO2 gamma-plot inner convex hull can be readily extracted from above:
- of course the equilirium shape is only bounded by 8 {101} surfaces - a result that can be easily deduced from Wulff construction. Here we demonstrate that the convexified xi-plot gives the equivalent result:
- However the time it requires to reach this state follows a power law according to Mullins and it's heavily dependent on the particle size. Since the transient shape is of interest here we only allow it to evolve for certain time. The result is shown below. Finally colour the surfaces by orientations. This is done in paraview (see 001.py, 112.py etc.)
ToDo
- kinetics of deposition by ballistic transport
Reference
[Mullins1957] | W.W. Mullins, Theory of Thermal Grooving, Journal of Applied Physics 28, 333 (1957) |
[Sekerka2005] | R.F. Sekerka, Analytical criteria for missing orientations on three-dimensional equilibrium shapes, Journal of Crystal Growth, Volume 275, Issues 1–2 (2005) |
[Herring1951] | C.Herring, `Surface Tension as a Motivation for Sintering, The Physics of Powder Metallurgy, p. 143, McGraw-Hill, New York (1951) |
[Torabi2009] | S.Torabi, J. Lowengrub, A. Voigt, S. Wise, A new phase-field model for strongly anisotropic systems, Proc. R. Soc., 465 1337-1359 (2009) |
[Siem2005] | E.J. Siem, W. C. Carter, Orientation-dependent surface tension functions for surface energy minimizing calculations, Journal of Materials Science, Volume 40, Issue 12, pp 3107–3113 (2005) |