The Role of Surface Tension and Mobility Model in Simulations of Grain Growth
Pith reviewed 2026-05-24 15:33 UTC · model grok-4.3
The pith
In three dimensions grain growth simulations produce substantially different grain size distributions for isotropic versus anisotropic reduced mobilities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The results indicate that in two dimensions, the different surface tension and mobility models do not play a significant role in the stationary grain size distribution. However, in three dimensions, there is a substantial difference between the distributions obtained from the same three models, depending on whether the reduced mobilities are isotropic or anisotropic. Additional results show that in three dimensions, the misorientation distribution function of a grain network with random orientation texture returns to the close vicinity of the Mackenzie distribution even if started very far from it.
What carries the argument
Threshold dynamics algorithms that capture the Herring angle condition at junctions and automatically handle topological transitions.
If this is right
- In two dimensions the stationary grain size distribution is insensitive to the choice of surface tension or mobility model.
- In three dimensions isotropic and anisotropic reduced mobilities yield substantially different distributions.
- The misorientation distribution function approaches the Mackenzie distribution from arbitrary initial conditions in three dimensions.
- Large scale simulations can be performed without manual handling of junction conditions or topology changes.
Where Pith is reading between the lines
- The robustness of the misorientation distribution suggests that random texture evolution is largely independent of starting misorientation statistics.
- Three-dimensional simulations may need anisotropic mobility models to accurately predict experimental grain size distributions in polycrystals.
- Similar threshold dynamics approaches could be applied to study grain growth under additional driving forces such as curvature or stored energy.
Load-bearing premise
The threshold dynamics algorithms correctly enforce the Herring angle condition at triple junctions and manage topological changes without distorting the long-time distributions.
What would settle it
Comparison of the three-dimensional stationary distributions obtained here against those from an independent numerical method such as the phase field approach to see if the isotropic-anisotropic difference persists.
read the original abstract
We explore the effects of surface tension and mobility models in simulations of grain growth using threshold dynamics algorithms that allow performing large scale simulations, while naturally capturing the Herring angle condition at junctions and automatically handling topological transitions. The results indicate that in two dimensions, the different surface tension / mobility models considered do not play a significant role in the stationary grain size distribution. However, in three dimensions, there is a substantial difference between the distributions obtained from the same three models, depending on whether the reduced mobilities are isotropic or anisotropic. Additional results show that in three dimensions, the misorientation distribution function of a grain network with random orientation texture returns to the close vicinity of the Mackenzie distribution even if started very far from it.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses threshold-dynamics algorithms to simulate grain growth and examines the influence of three surface-tension/mobility models. It reports that in 2D the stationary grain-size distributions are insensitive to model choice, whereas in 3D the distributions differ substantially when reduced mobilities are taken to be isotropic versus anisotropic. A second result is that, in 3D, the misorientation distribution function of a random-orientation network relaxes to the vicinity of the Mackenzie distribution regardless of the initial texture.
Significance. If the numerical fidelity of the method is confirmed, the 3D model-sensitivity result would indicate that anisotropy in reduced mobility must be retained for quantitative microstructure predictions, while the relaxation to the Mackenzie distribution supplies a falsifiable statement about texture evolution under curvature-driven motion.
major comments (2)
- [Abstract and §3] Abstract and §3 (Methods): the statement that the threshold-dynamics scheme “naturally captures the Herring angle condition at junctions and automatically handling topological transitions” is presented without any quantitative verification (analytic triple-junction test, comparison with phase-field or front-tracking codes, or measured angle error statistics). Because the headline 3D differences are attributed to the mobility models, any systematic discrepancy between the discrete rule and the continuum variational problem could produce model-dependent artifacts rather than physical effects.
- [§4] §4 (Results): the claim of a “substantial difference” between isotropic and anisotropic reduced-mobility cases in 3D is stated without reported error bars, Kolmogorov-Smirnov statistics, or convergence checks with respect to grid size or time-step. The absence of these controls leaves open the possibility that the observed spread is within numerical uncertainty.
minor comments (2)
- [Figures] Figure captions should explicitly state the number of grains, the number of independent runs, and the binning procedure used for the histograms.
- [§2] Notation for reduced mobility (isotropic vs. anisotropic) should be introduced once in §2 and used consistently thereafter.
Simulated Author's Rebuttal
We thank the referee for the careful review and constructive comments. We address each major comment below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [Abstract and §3] Abstract and §3 (Methods): the statement that the threshold-dynamics scheme “naturally captures the Herring angle condition at junctions and automatically handling topological transitions” is presented without any quantitative verification (analytic triple-junction test, comparison with phase-field or front-tracking codes, or measured angle error statistics). Because the headline 3D differences are attributed to the mobility models, any systematic discrepancy between the discrete rule and the continuum variational problem could produce model-dependent artifacts rather than physical effects.
Authors: We agree that the manuscript would be strengthened by explicit quantitative verification of the junction conditions. The threshold-dynamics algorithm is derived from a variational formulation whose continuum limit recovers the Herring condition and handles topology changes; this property is established in the foundational references for the method. Nevertheless, to address the concern directly, we will add a short verification subsection (including a triple-junction angle test and comparison to the expected Herring angles) to the Methods section of the revised manuscript. revision: yes
-
Referee: [§4] §4 (Results): the claim of a “substantial difference” between isotropic and anisotropic reduced-mobility cases in 3D is stated without reported error bars, Kolmogorov-Smirnov statistics, or convergence checks with respect to grid size or time-step. The absence of these controls leaves open the possibility that the observed spread is within numerical uncertainty.
Authors: We acknowledge that §4 does not report error bars, statistical tests, or explicit convergence data. The observed differences between isotropic and anisotropic cases are visually large and appear consistently across independent runs, but we agree that quantitative controls are required to rule out numerical uncertainty. In the revision we will add (i) standard-error estimates obtained from multiple independent realizations, (ii) Kolmogorov-Smirnov statistics comparing the distributions, and (iii) a brief grid-convergence study for the 3D cases. revision: yes
Circularity Check
No circularity: results are direct numerical outputs from simulations
full rationale
The paper reports grain-size and misorientation distributions obtained by running threshold-dynamics simulations under three surface-tension/mobility models. These quantities are generated by the algorithm itself and compared to the externally known Mackenzie distribution; no fitted parameters, self-definitional relations, or load-bearing self-citations reduce any reported stationary distribution to an input quantity by construction. The derivation chain consists of standard numerical evolution followed by post-processing statistics, which remains self-contained against external benchmarks.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.