Uniformity Assessment of TRISO Fuel Particle Distribution in Spherical HTGR Fuel Element Using Voronoi Tessellation and Delaunay Triangulation.
The spherical fuel element (SFE) is the key component in the Chinese high-temperature gas-cooled reactor (HTGR) project. Usually, an SFE is stuffed with about 12,000 tri-structural-isotropic (TRISO) fuel particles [1, 2]. Ideally, the TRISO particles are uniformly distributed, giving a good temperature distribution, as well as a homogeneous temperature gradient field [3, 4]. An illustration of the SFE for HTGR is shown in Figure 1. The diameter of the fuel zone is 50 mm encased with a 5 mm fuel-free outer graphite matrix. The nonuniform distribution of the TRISO particles in the SFE may lead to a nonuniform thermal expansion and excessive thermal gradients. Close clusters of particles may form local hotspots, increasing the stress on TRISO particle layers as the stress is directly dependent on the temperature in pressure vessel model of the particle . Consequently, the failure probability of the TRISO particle would increase, leading to the release of fission products. Thus, it is necessary to measure and evaluate the TRISO particle distribution to ensure they are relatively uniform before operation. Monte Carlo method was used to simulate TRISO particle distributions and the distance of nearest neighbors was calculated by Nabielek . In our previous research, 3D visualization of particles in an SFE was implemented and 3D volume rendering was obtained for subjective evaluation, as shown in Figure 2. The integrity of SFE, i.e., whether there were fabrication defects, was discussed in . In order to assess the spatial uniformity of the particles, we propose in this study to process the TRISO particles in a 3-dimensional (3D) coordinate system and assess the uniformity of the TRISO particle distribution using the geometric metrics extracted from the Voronoi tessellation and Delaunay triangulation. The number of particles and their locations were obtained by the use of an image-processing method, after a cone-beam computed tomography (CBCT) was used for data acquisition. Then, the Voronoi tessellation and Delaunay triangulation were constructed for the assessment of the distribution uniformity of the TRISO particles.
A Voronoi tessellation, or Voronoi diagram, partitions a space into a set of mutually exclusive regions, with a number of given seeds. Each seed locates at the center of one region, also known as a Voronoi cell. The region is composed of the set of all points that are closer to that seed than any others. A Delaunay triangulation is simply the dual graph for a Voronoi tessellation [8,9]. The Voronoi diagram was successfully used in cylindrical particle packing by Rycroft . Herein, nearest neighbor distance and Voronoi cell volume were investigated to describe the characteristics of the spatial distribution of TRISO particles.
2.1. Nondestructive Data Acquisition and 3D Image Reconstruction. X-ray inspection is a commonly used qualitative evaluation technique for SFEs . In this paper, CBCT was employed as a powerful and quantitative tool to evaluate the TRISO particle distribution nondestructively [12-14].
The sketch and a photo of the CBCT system are shown in Figures 3(a) and 3(b), respectively. A 160 kVp micro focus X-ray generator (Philips HOMX-161) was used as the X-ray source. The tube voltage and the current were set to be 140 kV and 0.2 mA, respectively. The focal size was adjusted to 20-50 [micro]m. A flat panel (Varian PaxScan 2520) was used as the X-ray photon detector. The effective spatial resolution (pixel size) was 80 [micro]m . As shown in Figures 3(a) and 3(b), the SFE is placed between the X-ray source and the detector. X-ray was emitted from the source and penetrated the object of interest before being detected by the flat-panel detector. Computed tomography (CT) images, or cross-section images, can be reconstructed from the projection data (attenuation of the X-ray by the object of interest). To acquire projection data at different angles, the SFE was rotated. The source-to-isocenter (rotation center) distance was 25 cm and the source-to-detector distance was 75 cm. Projection data were obtained at one-degree interval for 360[degrees]. Thus 360 sets of projection data were collected at different angles. Tomographic images were reconstructed using the fast FDK algorithm as shown in Figure 3(c) and a 3D image was generated nondestructively.
2.2. Construction of the Voronoi Tessellation and Delaunay Triangulation. We developed an image segmentation method to identify the centroids of TRISO particles in the 3D reconstructed SFE tomographic image. With the coordinates of the TRISO particle centroids, the 3D Voronoi tessellation and Delaunay triangulation were constructed [8, 9]. Note that the size of the TRISO particles was ignored hereafter to simplify the analysis. In the following steps, all TRISO particles are represented by their centroids.
We take a 2-dimensional (2D) case as an example to demonstrate how the Voronoi tessellation and Delaunay triangulation may be used to evaluate the uniformity of a group of particles. Randomly distributed particles and the corresponding Voronoi tessellation and Delaunay triangulation are shown in Figure 4. Figure 4(c) shows the dual relationship clearly. It can be imagined that the Voronoi cells and the edges of the Delaunay triangles tend to be similar if the particles are uniformly distributed.
3. Metrics for Assessment of Uniformity of TRISO Particles
To investigate the relation between the local density of particles and geometric characteristics of the Voronoi tessellation, five metrics were extracted, including the volume of Voronoi cell, number of neighbors, area of Voronoi facet, volume of Delaunay simplex, and length of simplex edges . Intuitively, one expects the spread of the distributions of the metrics to be narrow if particles are uniformly distributed. In other words, the uniformity of the particle is closely related to the width of the distributions. We may thus quantify the uniformity of the particles by characterizing the distributions of the metrics.
Every TRISO particle has its own Voronoi cell during the tessellation construction. The number of neighbors indicates how many TRISO particles are related/connected in the specified field. The nearest neighbor indicates how close a TRISO particle is to its neighbor. The volume of Voronoi cells indicates the mean size of a region TRISO particles control. The volume of Voronoi cells is smaller in areas where particles are more densely packed.
We use 2D examples to demonstrate how the aforementioned metrics may help us assess the uniformity of the particle distribution. As shown in Figure 5, for instance, the number of neighbors is equal to 6 for uniformly distributed particles, while the same metric tends to be different in the random example.
In Figure 6, 121 particles randomly distributed in a 2D space are shown. The distributions of the five metrics are also sketched in the same figure.
4. Experimental Results
An SFE was processed using the above methods. The number of the measured TRISO particles in the SFE was 11,603, after image segmentation and region labeling. Given that the nominal number of TRISO particles is 12,000, the relative error is approximately 3% (397/12000). Table 1 is the xyz position of each TRISO particle in orthogonal coordinate system (Cartesian). The location unit is pixel which is equal to about 80 [micro]m. Voronoi volume distribution of all particles is shown in Figure 7, close to a log-normal distribution for randomly uniformly distributed particles in theory.
The distribution of the distance to the nearest neighbor of TRISO particles is shown in Figure 8, which is also close to a log-normal distribution for randomly uniformly distributed particles in theory .
The experimental results demonstrated that the distribution of the TRISO particles is highly similar to the lognormal distribution, which represents randomly uniformly distributed particles . The majority of the particles had 5 or 6 neighbors, the most probable volumes of Voronoi cells were between 0.003 and 0.015, the most probable areas of Voronoi facets were between 0.2 and 0.6, the most probable Delaunay simplex volumes were less than 0.01, and the lengths of simplex edges peaked at 0.01 and roughly symmetrically distributed between 0 and 0.2.
In this study, uniformity assessment method based on the Voronoi tessellation and the Delaunay triangulation gives quantitative results in a 3D volume rendering for TRISO particles in an SFE. The metrics extracted from the Voronoi tessellation and Delaunay triangulation demonstrated that, in a real SFE, the locations of the TRISO particles closely followed a randomly uniform distribution. A more rigorous algorithm is under development, in which statistical test will be used to reach a more convincing conclusion in terms of the uniformity of the TRISO particle distribution inside a 3D spherical space.
The figures and experimental setup data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
This work was jointly supported by the National Natural Science Foundation of China (No. 61571262, No. 11575095, and No. 61171115) and National Key Research and Development Program (No. 2016YFC010 5406). Y. P. was also supported in part by the Fundamental Research Funds for the Central Universities through Beijing Jiaotong University (W15JB0030). The authors would like to thank Professor Bing Liu at Institute of Nuclear and New Energy Technology at Tsinghua University for providing the surrogate fuel pebbles and Mr. Stephen Conway for grammar correction.
 Z. Zhang, Z. Wu, D. Wang et al., "Current status and technical description of Chinese 2 x 250 MWffc HTR-PM demonstration plant," Nuclear Engineering and Design, vol. 239, no. 7, pp. 1212-1219, 2009.
 M. A. Futterer, L. Fu, C. Sink et al., "Status of the very high temperature reactor system," Progress in Nuclear Energy, vol. 77, pp. 266-281, 2014.
 "IAEA-TECDOC 1674: Advances in High Temperature Gas Cooled Reactor Fuel Technology," Vienna, 2012.
 K. Verfondern, H. Nabielek, and J. M. Kendall, "Coated particle fuel for high temperature gas cooled reactors," Nuclear Engineering and Technology, vol. 39, no. 5, pp. 603-616, 2007.
 K. Verfondern, "Triso fuel performance modeling and simulation," Comprehensive Nuclear Materials, vol. 3, pp. 755-788, 2012.
 H. Nabielek, TRISO Particles in The Fuel Zone: What Is the Limiting Packing Fraction? IAEA, Vienna, Austria, 2013.
 G. Yu, Y. Du, X. Xiang, Y. Liu, Z. Li, and X. Wang, "3D nondestructive visualization and evaluation of TRISO particles distribution in HTGR fuel pebbles using cone-beam computed tomography," Science and Technology of Nuclear Installations, vol. 2017, Article ID 3857075, 6 pages, 2017.
 Y. H. Peng, X. G. Wang, Y. Du, and H. J. Chen, "Using voronoi tessellation and Delaunay triangulation to evaluate spatial uniformity of particle distribution," Applied Mechanics and Materials, vol. 614, pp. 413-416, 2014.
 R. Y. Yang, R. P. Zou, and A. B. Yu, "Voronoi tessellation of the packing of fine uniform spheres," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 65, no. 4, 2002.
 C. H. Rycroft, "VORO++: a three-dimensional voronoi cell library in C++," Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 19, no. 4, 2009.
 D. Yi, X. Wang, X. Xiang, and L. Bing, "Automatic X-ray inspection for the HTR-PM spherical fuel elements," Nuclear Engineering and Design, vol. 280, pp. 144-149, 2014.
 S. Knol, P. R. Hania, B. Janssen, M. Heijna, and S. De Groot, "PYCASSO: X-ray tomography on HTR coated particles," Nuclear Engineering and Design, vol. 271, pp. 206-208, 2014.
 P. R. Hania, B. Janssen, A. V. Fedorov et al., "Qualification of HTR pebbles by X-ray tomography and thermal analysis," Nuclear Engineering and Design, vol. 251, pp. 47-52, 2012.
 G. Lecomte, D. Tisseur, J. M. Letang, J. Banchet, and M. P. Vitali, "Quantitative homogeneity and in-contact particles of High Temperature Reactors (HTR) compacts determination via X-ray tomography," in Proceedings of the 34th Annual Review of Progress in Quantitative Nondestructive Evaluation,pp. 1646-1653, 2008.
 X. G. Wang, Y. Ye, and Z. P. Guo, "Two typical 3D computed tomography systems and application," CT Theory and Applications, vol. 18, no. 4, pp. 26-33, 2009.
Libing Zhu, (1) Xincheng Xiang, (1) Yi Du, (1) Gongyi Yu, (1) Ziqiang Li, (1) Yahui Peng (ID), (2) and Xiangang Wang (ID) (1)
(1) Institute of Nuclear and New Energy Technology, Tsinghua University, Beijing 100084, China
(2) School of Electronic and Information Engineering, Beijing Jiaotong University, Beijing 100044, China
Correspondence should be addressed to Yahui Peng; email@example.com and Xiangang Wang; firstname.lastname@example.org
Received 20 June 2018; Accepted 6 September 2018; Published 1 October 2018
Academic Editor: Eugenijus Uspuras
Caption: Figure 1: Sketch of a spherical fuel element used in HTGR stuffed with TRISO fuel particles.
Caption: Figure 2: 3D volume rendering of TRISO particles in a spherical fuel element.
Caption: Figure 3: (a) Scheme of 3D cone-beam CT, (b) facility, and (c) a tomographic image of a spherical fuel element.
Caption: Figure 4: An example of the Voronoi tessellation (cyan solid lines, (a)), and its dual graph, the Delaunay triangulation (red dashed lines, (b)) for the same set of seed points (red solid dots). The overlap of both aforementioned graphs to illustrate the duality (c).
Caption: Figure 5: Examples of a uniformly distributed seed pattern (a) and a randomly distributed seed pattern (b) with both Voronoi tessellation and Delaunay triangulation.
Caption: Figure 6: The distributions of the five metrics (from top to bottom, right to left): number of neighbors, volume of Voronoi cells, area of Voronoi facets, volume of Delaunay simplices, and length of simplex edges for randomly distributed particles in a 2-dimensional space.
Caption: Figure 7: Voronoi volume distribution (Influential Zone).
Caption: Figure 8: Nearest neighbor distances.
Table 1: Mass center (XYZ) of every TRISO. Index X Y Z 1 925.500 817.800 6.300 2 874.204 844.408 5.633 3 934.897 850.759 2.897 4 905.167 862.333 2.000 5 935.704 885.481 2.074 6 793.854 899.902 9.073 7 865.333 922.667 5.429 8 1020.160 927.520 3.360 9 770.222 937.778 6.000 ...... ...... ...... ...... 11603 1034.489 981.244 1503.200
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Zhu, Libing; Xiang, Xincheng; Du, Yi; Yu, Gongyi; Li, Ziqiang; Peng, Yahui; Wang, Xiangang|
|Publication:||Science and Technology of Nuclear Installations|
|Date:||Jan 1, 2018|
|Previous Article:||Modelling of Severe Accident and In-Vessel Melt Retention Possibilities in BWR Type Reactor.|
|Next Article:||The R&D of HTR-STAC Program Package: Source Term Analysis Codes for Pebble-Bed High-Temperature Gas-Cooled Reactor.|