Additive Manufacturing Powder Particle Size Distributions: Comparison of Histogram Binning Methods

— Additive manufacturing powders require a well-defined Particle Size Distribution (PSD) and spherical morphology to ensure good flowability. To simplify characterisation, powders can be prepared using standard metallurgical techniques followed by optical imaging of the cross-sectioned particles. Measured PSDs of particle sections are typically underestimates of the true PSD; hence, stereological corrections are required. Variations arise in the histogram binning methods (central binning versus upper limit binning) of commonly used stereological corrections. Although the results show some sensitivity to the binning method used, the GCO method seemed reasonably robust to changes in the binning method. However, authors are encouraged to follow the method as it is intended within the literature, which was found to be especially true when using Saltykov’s method.

show some sensitivity to the binning method used, the GCO method seemed reasonably robust to changes in the binning method.However, authors are encouraged to follow the method as it is intended within the literature, which was found to be especially true when using Saltykov's method.

I. INTRODUCTION
A well-graded Particle Size Distribution (PSD) of spherical particles is necessary for good powder flowability and particle packing in powder-based additive-manufacturing processes [1,2].Hence, users of powder-based additive-manufacturing processes must regularly assess the PSD to ensure optimal part production.However, many users are Small and Medium-sized Enterprises (SMEs) without direct access to dedicated particle analysis facilities such as Laser Size Diffraction (LSD) or Scanning Electron Microscopy.This introduces the need to outsource analysis to third-party laboratories, increasing costs and lead times.
Preparing powders for optical analysis using standard metallurgical procedures of hot-mounting, grinding, and polishing provides a low-cost procedure of analysis with a quick turnaround, thereby reducing the frequency of thirdparty outsourcing.However, as shown in Fig. 1, the diameters of the particle sections observed on the polished plane are typically an underestimate of the true particle size.Hence, stereological corrections must be applied to the apparent PSD of the sections to estimate the true PSD.
Stereology involves using probability, statistical measures, and integral geometry to relate 2D measurements to parameters defining 3D structures.In this case, a histogram of measured particle section profiles is used to derive a histogram of spheres sizes which could have produced the profile histogram via a reconstruction of sphere size distribution approach.However, a question that arises is which binning method to use.Some techniques use central binning whilst others use upper limit binning.
The first reconstructive procedure of this kind was proposed by Wicksell [3] who considered the probability of finding an apparent section diameter, , centered within a bin and derived a transformation matrix for 15 bins in which the bin limits were ( − 0.5)Δ and ( + 0.5)Δ where  represents the bin number, and Δ represents the bin size.Wicksell also had a 1/4ℎ bin with a range of 0 − 0.5Δ.Saltykov [4] proposed a successive subtraction method based on a method developed by Scheil [5] which had been modified by Schwartz [6].Saltykov constructed a histogram using bin limits of ( − 1)Δ and Δ.The upper limits of the bins represented the profile and particle size for each class.Goldsmith [7] and Cruz-Orive [8] each independently improved these previous methods, resulting in a method that used the same bin limits as Saltykov but used class midpoints as a measure of the particle size, similar to Wicksell's method.
In some cases, binning methods are used interchangeably without explanation or critical evaluation.Applications of Saltykov's method (an upper limit method) have been found in literature using the central binning method [9].Differences in outcomes have not been investigated to determine whether it is an acceptable practice to change binning methods.Changes to binning is therefore investigated in this work.

II. MATERIALS AND METHODS
Plasma-atomised Grade 23 Ti-6Al-4V titanium alloy powder was used in this study.The sample preparation, digital imaging and image analysis were conducted as described in refs.[10,11].LSD analysis was performed to provide a ground truth for comparison purposes.The apparent PSD was generated via automated image analysis.The Saltykov (SSS) and Goldsmith-Cruz-Orive (GCO) methods were compared using both central and upper limit binning.Mean Absolute Errors (MAE) between each stereological correction output and ground truth (LSD) data were calculated.

III. RESULTS
The resulting cumulative size distributions from the SSS and GCO methods, plotted using upper and central binning methods, are shown in Fig. 2 with the cumulative size distributions of the particle sections (  ) and the LSD data.IV.DISCUSSION Although literature recommends using the coefficients of the GCO method rather than those of the SSS method [12], in this work, the SSSUpper method performs slightly better than the GCOCenter method.Both methods produce the same mean particle size (26.5 μm), skewness (1.0) and kurtosis (4.5) values, but the SSSUpper method results in a lower standard deviation (8.6 μm vs 9.0 μm) and lower MAE (2.3% vs 2.4%) when compared to the GCOCenter method.However, if the coefficients of the GCO method are used but the resulting   distribution placed at the upper bin limits (GCOUpper), the MAE value decreases from 2.4% to 1.3% which becomes the smallest MAE value, indicating the best fit to the ground truth (LSD) data, in this study.Furthermore, the GCOUpper method results in a mean particle size of 28.3 μm which, although larger than the mean from LSD, results in a smaller percentage difference of 1.4% (0.4 μm) compared to 5.0% for both the GCOCenter and SSSUpper methods (1.4 μm), and 11.8% for the SSSCenter method (3.3 μm).Using the coefficients of Saltykov's method and plotting at the central bin limits is not advised in this case as it produces a PSD with greater MAE than even the apparent PSD.
V. CONCLUSIONS Although results show limited sensitivity in the GCO method to changes in the binning method used (central binning or upper limit binning), seemingly, no significant differences could be established.However, authors are encouraged to follow any method as it is described in the literature.This is particularly recommended for the application of Saltykov's method where upper limit binning must be used as originally intended.As shown here, Saltykov method with central binning gave an erroneous result.

Fig. 1
Fig. 1 Diagrams in (a) 3D and (b) 2D showing how particle section diameter, , is typically an underestimate of true particle diameter, .

Fig. 2
Fig. 2 Cumulative fraction line graphs comparing the particle size distributions obtained from 2D apparent section sizes, LSD analysis, and stereological corrections via the SSS and GCO methods.Measures of central tendency and dispersion for each PSD are provided in Table I along with the calculated MAE.

manufacturing powders require a well- defined Particle Size Distribution (PSD) and spherical morphology to ensure good flowability. To simplify characterisation, powders can be prepared using standard metallurgical techniques followed by optical imaging of the cross-sectioned particles. Measured PSDs of particle sections are typically underestimates of the true PSD; hence, stereological corrections are required. Variations arise in the histogram binning methods (central binning versus upper limit binning) of commonly used stereological corrections. Although the results
Additive Manufacturing Powder Particle Size Distributions: Comparison of Histogram Binning Methods C. Gallagher 1, * , E. Kerr 1 , and S. McFadden 1 1 Faculty of Computing, Engineering and the Built Environment, Ulster University, Derry/Londonderry, Northern Ireland *Contact: Gallagher-c64@ulster.ac.uk, phone +44-28 7167 5522 keywords: Additive Manufacturing, Stereology, Particle Size Distribution

Table I
Measures of central tendency and dispersion with mean absolute error calculation for each size distribution.