A thermal analysis was performed for the Advanced Gas Reactor test experiment (AGR-3/4) with time-varying gas gaps. The experiment was irradiated at the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). Several fuel irradiation experiments are planned for the AGR Fuel Development and Qualification Program that supports the development of the Very-High-Temperature Gas-Cooled Reactor (VHTR) under the Next-Generation Nuclear Plant (NGNP) project. AGR-3/4 combines two tests in a series of planned AGR experiments to test tristructural-isotropic (TRISO)-coated, low-enriched uranium oxy-carbide fuel. Forty-eight TRISO-fueled compacts were inserted into 12 separate capsules for the experiment (four compacts per capsule). The purpose of this analysis was to calculate the temperatures of each compact and graphite layer to obtain daily average temperatures using time (fast neutron fluence)-varying gas gaps and to compare with experimentally measured thermocouple data. A finite-element heat transfer model was created for each capsule using the commercial code ABAQUS. Model results are compared to thermocouple data taken during the experiment.
The AGR-1 and AGR-2 experiments irradiated in the Advanced Test Reactor (ATR) have previously been modeled for daily thermal evaluation in Refs.  and . The author discusses similar topics in these references to this article, such as variable gas gaps, mesh sensitivity, thermal conductivity varying with fast neutron fluence, and daily thermal heat rates imported from physics analysis. These first two experiments were shake-down tests for the TRISO fuel particles in compacts.
The AGR-3/4 experiment is composed of 12 individual capsules, approximately 0.06985 m (2.75 in.) diameter by 0.1016 m (4.00 in.) tall, stacked on top of each other to form the test train. Each capsule contains four TRISO-particle compacts that are approximately 0.0127 m (0.5 in.) diameter by 0.0127 m (0.5 in.) long. The compacts are composed of TRISO fuel particles bound together by a carbon matrix. Each compact contains approximately 4150 fissile particles (36 vol% particle packing fraction). Each capsule is supplied with a flowing helium/neon gas mixture to control the test temperature and sweep any fission gases that are released to the fission product monitoring system. Temperature control is accomplished by adjusting the gas mixture ratio of the two gases (helium and neon) with differing thermal conductivities in the gas gaps.
The AGR-3/4 experiment was placed in the north-east flux trap position in the ATR core as shown in Fig. 1. An axial view of two of the 12 capsules is shown in Fig. 2. Four through tubes carry thermocouples (TCs) and gas lines to each individual capsule. All 12 capsules have their own gas mixture and fission product gas return line. Figure 2 also shows the arrangement of how the capsules are stacked together to form the experiment train. Each capsule has a fuel compact in the center surrounded by three graphitic annuli as shown in Fig. 3. Symmetry may not be assumed as the heat rates vary azimuthally for each graphite layer. The graphite annuli proceeding from the compact out are the matrix ring (inside), graphite sleeve, and graphite sink (outside). The largest gas gaps are shown on the outside and inside of the graphite sink. Gas gaps for the two innermost gaps are very small and not perceptible in Fig. 3. Each of these four components is designed to operate at a specified temperature in all 12 capsules. Through tubes containing gas lines and TC wires are displayed in Fig. 3. Four temperature control gas gaps separate these components for a total of 48 gas gaps in the entire test train. Reactor coolant water flows on the outside of the stainless steel capsule shell.
The commercial finite-element heat transfer code ABAQUS  was used for this analysis. The ABAQUS model has a direct volume-for-volume correlation with the physics model. A similar physics model is discussed in Ref.  for the heating of the compacts (each compact is evenly and axially divided into two equal parts). The goal of these predictions is to be able to adjust the TC set points as the fuel burns during the experiment so as to maintain constant fuel temperature.
Numerical Model and Discussion
A cutaway rendering of a typical capsule is shown on the left side of Fig. 4. One of the main goals of this experiment was to make the bulk of the heat from the compacts flow radially out of the capsule instead of axially out of the ends of the capsule. Zirconia, grafoil, and graphite felt insulators are placed on the top and bottom of each capsule shown in Fig. 4. The finite-element mesh with a cutaway view colored by different materials of the entire model is shown on the right side of Fig. 4. A Cartesian coordinate system is appropriate for this model because of the three-dimensionality of the heat flow. Approximately, 400,000 eight-noded hexahedral brick elements were exclusively used in all 12 capsule models. Several mesh convergence studies [1,2] have been performed on the mesh. Identical agreement for this mesh and a mesh with twice as many elements in each direction was obtained. The exclusive use of eight-noded hexahedral brick elements ensures good numerical accuracy for this type of heat conduction model as discussed in Refs.  and . The finite-element models have a mesh density in the radial direction of approximately 38 finite elements per inch with a characteristic length of 0.965 mm (0.038 in.). The gap conductance model was implemented for the outside three gas gaps, while the innermost gas gap had hexagonal brick elements. This innermost gap was modeled with brick elements as all the capsules had the same gas gap. As only one basic mesh was created and propagated to the other 11 capsules, various gas gap conductivities and gap conductances were implemented by taking into account each individual gap dimension. The top and bottom of each model were assumed to be adiabatic. This implies that we are ignoring radiation heat transfer from the top of one capsule to the bottom of the one above. The gas gap between capsules is more than 0.0127 m (0.5 in.).
The fuel compact thermal conductivity was taken from correlations presented from Gontard and Nabielek in Ref. , which gives correlations for conductivity, taking into account temperature, temperature of heat treatment, neutron fluence, and TRISO-coated particle packing fraction (where packing fraction is defined as the total volume of particles divided by the total volume of the compact).
For fluences greater than (), the conductivity increases as fluence increases for higher temperatures because of the annealing of radiation-induced defects in the material with high temperatures, while the opposite occurs at lower temperatures. The thermal conductivity of the matrix ring was taken from the fuel compacts correlation with a fuel particle packing fraction of zero. This was done since a pure matrix material conductivity was not available. A plot similar in shape to the fuel compacts, with higher conductivity, is shown in Fig. 6.
Figure 8 shows a plot of the helium/neon sweep gas thermal conductivity versus temperature and mole fraction of helium. The thermal conductivity increases as the helium mole fraction increases and as the temperature increases. Heat produced in the fuel compacts and graphite components is transferred through the gas gaps surrounding the compacts and components via a gap conductance model using the gap width and the conductivity of the sweep gas as discussed later. Both radiation and conduction heat transfer were considered across every gap. However, because the thermal capacitance of the sweep gas is very low (), advection is not considered in the sweep gas. The sweep gas is modeled as being stationary. The convective heat transfer from these sweep gases would be less than 0.01% of the heat transfer across the gap because of the low density, low flow rate, and low thermal capacitance. The thermal conductivity of the sweep gas mixture was determined using a set of correlations from Brown University  for mixtures of noble gases.
The neon gas fraction for each day was calculated for each capsule using average daily flow rates for helium and neon through each capsule. Graphite and fuel compact material properties vary with neutron fluence. Fluence was imported from the detailed physics daily as-run calculations. The ABAQUS Field Variable model was implemented where the neon fraction was taken as Field Variable 1, and fast neutron fluence was taken as Field Variable 2. ABAQUS thus provides a method of properties being able to vary with fields other than temperature.
The gamma/neutron heating for the various components (including the fuel compacts) was taken from the as-run physics calculations. Typical heat rates for the fuel compacts and graphite components were nominally 120 and , respectively. The test train heat rates exhibit the typical chopped cosine profile that is distinctive of ATR.
All gas gaps were modeled as changing with fast neutron fluence. This was accomplished by having the gas gap conductivity of each capsule change with neutron fluence. Fluence was set as Field Variable 2 in the ABAQUS model. The original finite-element mesh models created in ABAQUS were done with the as-built dimensions for the gas gaps. The gas gaps were assumed to be the hot gas-gap dimension, the hot gas-gap dimension and room temperature gas-gap dimension being virtually the same. Experimental data  obtained from the AGR-1 experiment were used for the compact dimensional change, while data from the AGC-1 experiment  were used for the graphite shrinkage. Matrix dimensional change data were taken from Refs.  and . Table 1 shows the values obtained from measured values from the AGR-1 compact, matrix ring, and graphite in the ring and sink. The table shows two different types of graphite used in the experiment, PCEA and IG-110. This model assumes that both the inside and outside radius of the annuli shrink inward. This assumption is validated by graphite annuli measured from the Peach Bottom reactor since the 1970s .
|Fast fluence ()||AGR-1 compacts||Matrix||PCEA||IG-110|
|Fast fluence ()||AGR-1 compacts||Matrix||PCEA||IG-110|
Figures 12–18 show the results of a small sampling of the results for the first six ATR cycles 151A, 151B, 152B, 154A, 154B, and 155A. Figures 12 and 13 show temperature contours of the various components, while Fig. 14 shows a temperature history plot compared with actual TC measurements. Figures 15–18 show sensitivity to the variable gap model predictions.
A cutaway view of the temperature contours and mesh is shown for capsule 12 (typical) in Fig. 12. The maximum fuel compact temperature is 887°C at the center. Outside stainless steel capsule temperatures are near the temperature of the ATR primary coolant water temperature of 50°C. Gamma heating in the stainless steel end cap shows a radial temperature gradient. Several insulating materials have been placed in the model to prevent heat from transferring in the axial direction and out of the stainless steel end caps. The majority of the heat for these capsules is deposited in the fuel compacts () and the three graphitic ring layers ().
Figure 13 shows temperature contour plots for the () fuel compacts, () matrix, () graphite ring, and () graphite sink. One goal of this experiment is to have as uniform temperature as possible in the fuel compacts and graphite rings. The majority of the compact is between 820 and 870°C as shown in Fig. 13(a). The very center is the hottest and outside edges are cooler as is typical for a heat-generating cylinder with heat transfer on all sides.
The matrix ring temperature contours are shown in Fig. 13(b). Almost the entire matrix ring is between 765 and 800°C. Similar results are shown in Fig. 13(c) for the graphite sleeve, with the vast majority at . The highest temperature in this component is at the very bottom (inside) (not shown). This occurs since the fuel compacts, matrix ring, and graphite sleeve all sit on a thin layer of grafoil that is fairly conductive, yet nonreactive with the materials contacting it. Coolest temperatures are at the top-outside corner.
Figure 13(d) shows the graphite sink temperature contours without the top and bottom lid. Median temperature is 495°C with minimum and maximum . Hot spots occur at the inside in the four locations where the through tubes prevent the heat from evenly transferring to the outside. Coolest temperatures are on the top-outside edges next to the through tube holes. Gamma heating for all of these annular components was implemented in the 90° segments. It appears that the azimuthal temperature variations are very small.
A temperature history plot is shown in Fig. 14 for capsule 2. The TC location for all capsules is in the north-east center of the graphite sink (fourth layer). The top panel plots the actual TC measurements for TCs 1 and 2. The horizontal axis shows the ATR cycle numbers. Each cycle represents approximately 55 effective full-power days. The vertical axis shows the temperature in °C for the TC measurements in panel 1, whereas the temperature difference (measured minus calculated) is shown in panel 2 for TC1. The higher set of triangular shaped dots in the second panel show the temperature difference for a graphite emissivity of 1.0, while the lower set of square shaped dots show emissivity of 0.9 and the heat rates lowered by a multiplying factor. This multiplying factor takes into account the amount of heat calculated to flow into the ATR coolant water in the north-east lobe. Experience has shown that multiplying by this ratio more accurately depicts temperatures in the TCs. These lower set of square shaped dots show temperature predictions of the TCs vary by less than 20°C difference. The breaks in the plot represent times during the cycles that the reactor was shut down and also at the beginning and end of each cycle. There is no simple explanation as to why some capsules predict higher and some lower. Several different inputs are used to calculate the heat rates with the MCNP code such as the driver core fuel burnup, control rods, and outer shim rotation angles. Gas gaps are changing. Vibrations may cause some graphite layers to move by a very small amount toward one direction or another inside the annulus. An ongoing effort is being made to finalize the modeling technique necessary to minimize the difference between measured and calculated TCs. The temperature difference between the measured minus calculated has a maximum of 20°C as the model is doing a good job of predicting the TCs.
A radial temperature profile of capsule 5 is presented in Fig. 15 varying with neutron fast fluence for constant heat rates from the third cycle, day 8 with a constant neon fraction of 0.5. Nearly a 250°C centerline temperature difference is observed between the start of irradiation and a fast neutron fluence of , . This occurs as the conductivity of the graphite decreases as irradiation progresses and is most notable in the matrix ring (second layer, next to compact). Gap 4 is also growing, causing an increasing temperature drop across the gap. This situation shown in Fig. 15 is unrealistic as the heat rate decreases with irradiation, and the gas mixture would be adjusted so as to maintain a constant TC reading but was included to show the sensitivity to fluence. In order to maintain control of the experiment, the entire north-east lobe power increases.
Referring to Fig. 15, a fluence of 0.0 has the lowest peak centerline temperature, whereas a fluence of 5.0 has the highest due to lower thermal conductivity and larger outer gap (gap 4). One interesting note is that gap 2 closes due to neutron shrinkage differences in the graphite and matrix for a fluence of 5.0. Figures 16–18 show results with a different set of runs performed with varying neon fraction in order to maintain the peak centerline temperature constant at 1150°C. Figure 16 shows gap 2 closing at a fluence of 5.0. Notice the temperature gradient in matrix material varying with fluence.
Figure 17 shows the temperature difference across all four gaps for the scenario with constant peak fuel centerline temperature. Gaps 1, 2, and 3 all have a that decreases with fluence. Gap 3 is the only graphite-graphite material of the four gaps, meaning that the dimensional change is constant for both materials. Gap 4 has increasing as the graphite shrinks away from the stainless steel capsule wall, except for a fluence of zero. The for gap 4 initially decreases due to the drastic change in the graphite conductivity between a fluence of 0 and 1.0.
Figure 18 shows the neon fraction varying with fluence for constant heat rate and constant fuel centerline temperature. As the outer gas gap grows, the neon fraction decreases in order to maintain constant temperature. Again, at a fluence of 5.0, gap 2 closes and causes the necessary neon fraction to be constant between a fluence value of 4.0 and 5.0.
Daily as-run thermal analysis has been performed for the AGR-3/4 fuel experiment for all 12 capsules during the first six ATR irradiation cycles of the experiment. A variable gas-gap model changing with fast neutron fluence was implemented. A 3-D finite-element heat transfer model was created to simulate this experiment in the ATR. Volumetric heat rates and fast neutron fluence were imported from a daily as-run detailed physics analysis. Thermal conductivity of the fuel compacts and graphite holders varied with fluence and temperature. Daily helium–neon gas mixtures were implemented into 12 models. Temperature contours of various components have been presented. Daily history plots of actual TC measurements have been compared to simulated results with these models for all 12 capsules for the first three cycles. The temperature predictions appear to correlate fairly closely with the actual TC measurements. The goal of these predictions is to be able to adjust the TC set points as the fuel burns during the experiment so as to maintain constant fuel temperatures.
This work was supported by the U.S. Department of Energy, NGNP Program, Idaho Operations Office Contract DE-AC07-05ID14517.
Funding also was supported by the U.S. Department of Energy, NGNP Program, Idaho Operations Office Contract DE-AC07-05ID14517.
radiation surface area
displacements per atom
view factor from surface 1 to 2
million electron volts
volumetric heat rate,
- , ,
fast neutron fluence,