
In addition to the seismic configuration, another in-
put to the Parmec model is configuration of cracked fuel
bricks within the graphite core. Due to years of expo-
sure to high temperatures and irradiation, some of the
fuel bricks within the reactor are cracking, causing them
to break into two pieces. The presence and configura-
tion of these cracks has an impact on the reaction of the
core to seismic loading. It is possible that up to 40% of
the fuel bricks will eventually crack, although it is diffi-
cult to determine or predict where and when cracks will
occur.
The relationship between crack configuration and
seismic response of core components is complex, hence
the Parmec model consists of many thousands of pa-
rameters and equations. In addition, there are over
102500 possible permutations of crack configuration, as-
suming 40% cracking. With each configuration requir-
ing around 2 hours to compute the seismic response
via Parmec, it is clearly impractical to generate data
for even a small percentage of them. Instead, industry
practice is to generate random configurations of cracks,
passing each through Parmec in order to build up a
stochastic distribution of the seismic response.
2.2. Previous Machine Learning Surrogate Model of
Parmec
In previous machine learning assessments of AGR
graphite core seismic analysis [Jones et al. (2022)], each
crack configuration is considered an individual data in-
stance, with the encoding of cracked bricks being the
input features and the response of core components to
the earthquake being the output labels. The Parmec
software generates a time-history of the earthquake re-
sponse for all of the thousands of components within
the core. For the sake of simplicity and focus, the MLS
model was trained to predict the earthquake response
for a single interstitial brick at a single time frame - see
Figure 2.
To summarise the features of the MLS, each instance
has an input size of 1988 with this being the number
of fuel bricks within the AGR graphite. This input was
arranged into a 3D tensor which retains physical posi-
tional relationships within the actual AGR graphite core
(Figure 3). Each element is either a 1, -1 or 0 repre-
senting a cracked brick, uncracked brick or ‘empty’ po-
sition. The 3-dimensional encoding of the input fea-
tures also allows the dataset to be used with a convolu-
tional neural network [Albawi et al. (2017)] which was
found to be the best performing type of machine learn-
ing model.
For the aforementioned study, a dataset of approx-
imately 8300 instances was created using the random
crack pattern generator and the Parmec software. Out
of these instances, 6300 (75%) were used for training
with the remaining 2000 samples retained for testing.
Figure 2: A Top Down Diagram of the AGR Graphite Core Parmec
Model. Bricks are arranged into channels of two different types: fuel
(blue) and interstitial (grey). Both types of channel are the same
height, with fuel bricks being stacked seven high and the shorter in-
terstitial bricks being stacked 12 high. The cracking status of all
1988 fuel bricks is included in the input features (whether the brick
is cracked or not) of the surrogate machine learning model. For the
output labels, only the earthquake response of the upper most inter-
stitial brick (orange) is predicted by the surrogate machine learning
model.
Figure 3: Visualisation of a 3-dimensional Feature Encoding. This
example represents a single instance with each data-point representing
a fuel brick. Yellow and black data points represent uncracked and
cracked bricks, respectively.
2.3. Data Augmentation
Data augmentation is frequently employed in classi-
fication problems within the field of machine learning
[Shorten and Khoshgoftaar (2019)], where the model
predicts a discrete category for each dataset instance.
A classic example of classification is in computer vi-
sion, where a 2D or 3D tensor representing an image
4