# Prediction of Neutron Yield of IR-IECF Facility in High Voltages Using Artificial Neural Network

Authors: A. Sadighzadeh; A. Salehizadeh; M. Mohammadzadeh; F. Shama; S. Setayeshi; S. A. H. Feghhi; S. M. Sadati; M. Rezaei; E. Haji Ebrahimi; G. H. Roshani
Journal: Journal of Engineering
 (2014)

---

## Abstract

Artificial neural network (ANN) is applied to predict the number of produced neutrons from IR-IECF device in wide discharge current and voltage ranges. Experimentally, discharge current from 20 to 100 mA had been tuned by deuterium gas pressure and cathode voltage had been changed from −20 to −82 kV (maximum voltage of the used supply). The maximum neutron production rate (NPR) of 1.46 × 107 n/s had occurred when the voltage was −82 kV and the discharge current was 48 mA. The back-propagation algorithm is used for training of the proposed multilayer perceptron (MLP) neural network structure. The obtained results show that the proposed ANN model has achieved good agreement with the experimental data. Results show that NPR of 1.855 × 108 n/s can be achieved in voltage and current of 125 kV and 45 mA, respectively. This prediction shows 52% increment in maximum voltage of power supply. Also, the optimum discharge current can increase 1270% NPR.


---

## Body

## 1. Introduction

Inertial electrostatic confinement fusion (IECF) device is a compact and simple structure for nuclear fusion researches by electrical discharge, which can operate in the pulsed or continuous mode. It consists of two concentric (or coaxial) electrodes in which usually the central one (cathode) is negatively high-voltage-biased and the outward electrode (anode) is grounded. In this configuration, strong electric fields between electrodes lead to iodinating the filling gas and then accelerate the created ions toward the center, where the electrons are placed in the opposite direction. As a result, rather hot and dense plasma is formed in the center of cathode. In this situation, the continuous nuclear fusion reactions occur, which are the results of the beam-target and beam-beam interactions depending on the working conditions, although the beam-beam interactions will be negligible compared to beam-target reactions in the plasma conditions in this study. Therefore IECF is considered as a source of hot and dense plasma, highly energetic ions, and fast neutrons (when using deuterium or mixture of deuterium-tritium gas). IECF device is an excellent apparatus because of its ability to generate fast neutrons with high-flux from a small source for many applications, such as medicine (e.g., boron neutron capture cancer therapy) [1–3], radiography or tomography of thick materials, space propulsion system [4], inspection system and explosive landmine detection [5], neutron activation analysis, mine and petroleum exploration, and security screening. Therefore neutron production rate optimization in this device will increase its efficiency for the above applications.IEC researches have been performed since the 1950s decade. Recently, over the last 10 years, many researches and developments on IECF neutron source have been done [2, 4, 6–8]. Most of these developments are created by the usage of the different ion sources and enhancement of voltage and current [2, 4, 6]. IR-IECF is the prototype IEC device that is structured at the Atomic Energy Organization of Iran (AEOI) and the preliminary experiments on neutron yield were performed [7]. D-D neutrons have lower energy (2.5 MeV) and yield compared with D-T neutrons (14 MeV) and therefore they can be moderated within a thinner moderation layer, and thermal neutron flux can be obtained with D-D neutrons. The mentioned benefit and other advantages of D-D neutrons caused that deuterium gas to be used in IECF device.Reported results in previous works suggest that using the D-T gas mixture and increasing the ion current via ion source would give a neutron rate of 7.3 × 1012 n/sec at 1.2 mTorr, 75 kV, and 1.5 A ion current [6, 8, 9]. Researchers in Wisconsin University, in order to increase and optimize the neutron production rate in the UW IECF device, studied effect of cathode’s size (diameter), geometry, and material composition on NPR [10]. Dietrich in MIT University investigated multigrid IEC device theoretically and experimentally for improving particle confinement [11]. In this study, for the first time, the ability of artificial neural network for modeling the IECF device in order to optimize the neutron production rate was investigated and the maximum neutron yield was obtained. In the recent work [12], dependence of neutron production rate (NPR) on cathode voltage and discharge current in IR-IECF device for deuterium gas has been studied.Objective of this work is to train an ANN based on the data presented in [12] to predict NPR at high voltages that cannot be achieved experimentally by the current facilities. The use of ANN is discussed in detail in the following sections to predict NPR for voltages up to 130 kV.


## 2. Experimental Setup

IR-IECF device consists of a spherical grounded anode, which serves as a chamber, and a grid cathode that are set concentrically as shown in Figure1. This device consists of 13.5 cm diameter stainless steel cathode, 41 cm diameter anode with a 60 cm diameter, and 60 cm height vacuum chamber [7].Figure 1
Schematic diagram of IR-IECF device and parameters used in this study.The IR-IECF facility accompanied with the cathode and the anode geometry is shown in Figure2. In the experiments, facility is adjusted in a desired gas pressure and then by changing the cathode negative voltage glow plasma inside the cathode grid is produced. Produced ions are accelerated toward the cathode, and most of them penetrate the hollow cathode spherical wire net. In this situation, fusion reactions through beam-beam and beam-target collisions occurred. By increasing the cathode voltage and decreasing gas pressure, desired dense plasma is obtained.(a) IR-IECF device, (b) cathode (13.5 cm diameter), and (c) anode (41 cm diameter).
(a)
(b)
(c)In order to measure the neutron production rate (NPR), a calibrated neutron probe LB6411, which consists of a cylindrical3He proportional counter tube, was placed in 60 cm distance from cathode center of IR-IECF device [12]. The cathode voltage and discharge current have a significant impact on the NPR. The experiments had been carried out with the tunable 82 kV DC high voltage. Maximum NPR was 1.46 × 107 n/s and had occurred when the voltage was 82 kV and the discharge current was 48 mA. In the next sections using artificial neural network and these experimental data the NPR has been predicted for these ranges: current discharge = 40 mA–120 mA and voltage = 50 kV–130 kV. The extrapolation beyond the range of experimental data is one of the most important applications of using ANNs to model the nonlinear unknown-response systems such as the studied IECF. This method has been used in many published works such as [13–15].


## 3. Neural Network Computations

Artificial neural networks (ANNs) are mathematical models inspired in the human brain. Ability of learning to use experimental data is the main characteristic of this technique. This model can predict behaviors and patterns from a finite set of experimental data, called the “training set” of the ANN [16]. In general, the structure of ANN typically is comprised of three layers including input, hidden, and output layers. First primary data is collected in input layer and then sent to different processing units (neuron), which constitute the hidden layer of the networks [17]. The structure of a neuron consists of two major terms of weight and activation function as shown in Figure 3.Figure 3
The neuron mathematical model.Activation function receives an argumentn and generates an output a. The input of network for each neuron is sum of all input values in which each is multiplied by its weight and a bias term. The output value can be calculated by feeding the network input into the activation function of the neuron. Many activation functions are applied in ANN model such as triangular basis, pure line, soft max, log sigmoid, tan sigmoid, and hard limit. To achieve an acceptable prediction, the ANN is trained to minimize the error between ANN output and experimental data.In this research, multilayer perceptron (MLP) neural networks were used. MLP is a feed-forward artificial neural network model that maps input data sets onto a set of appropriate outputs. What makes a multilayer perceptron a different tool to predict the nonlinear unknown-response systems is that each neuron uses a nonlinear activation function (such as tangent sigmoid), which was developed to model the frequency of action potentials, or firing, of biological neurons in the brain [18, 19]. The proposed MLP model is shown in Figure 4, where the inputs are discharge current (mA) and voltage (kV) and the output is NPR (n/s).Figure 4
Architecture for the proposed MLP model.The input to the nodet in the first hidden layer is given by [20, 21]
(1)
η
t
=
∑
u
=
1
2
X
u
W
u
t
+
θ
t
,
t
=
1,2
,
…
,
10
.
The output from the tth neuron of the second hidden layer is given by
(2)
O
t
=
f
η
t
,
where X is the input variables, θ is the bias term, W is the weighting factor, and f is the activation function of the hidden layer. The output of the neuron in the output layer is given by
(3)
Y
=
∑
u
=
1
6
O
u
W
u
+
b
.
Required data set for training the network is obtained by use of IR-IECF experimental data. The presented MLP network is trained by Levenberg-Marquardt (LM) algorithm [22]. In this method, first derivative and second derivative (Hessian) are used for network weight correction [23]. The numbers of samples for training and testing data are 36 (about 70%) and 15 (about 30%), respectively. In this study, for optimizing the ANN configuration many different structures with one hidden layer and two, three, and four hidden layers were tested by applying different number of neurons in each layer and different epochs. MATLAB 7.0.4 software was used for training the ANN model. Table 1 shows the specification of the proposed ANN model being used in this study. In this table, an epoch represents completion of processing on the entire training data set [23].Table 1
Specification of the proposed ANN model.
Neural network
MLP
Number of neurons in the input layer
2
Number of neurons in the first hidden layer
10
Number of neurons in the second hidden layer
6
Number of neurons in the output layer
1
Number of epochs
1200
Activation function
tansig


## 4. Results and Dissection

Table2 shows the obtained errors for the proposed ANN model, where the mean relative error percentage (MRE%) and the roots mean square error (RMSE) are calculated by
(4)
RMSE
=
∑
m
=
1
N
X
m
Exp
-
X
m
Pred
2
N
0.5
,
MRE
%
=
100
×
1
N
∑
m
=
1
N
X
m
(
Exp
)
-
X
m
(
Pred
)
X
m
(
Exp
)
,
where N is the number of data and “X
(
Exp
)” and “X
(
Pred
)” stand for experimental and predicted (ANN) values, respectively.Table 2
Obtained errors for training and testing results of the proposed ANN model.
Error
Train
Test
MRE%
0.0072
0.6563
RMSE
0.0445
0.0719Figure5 shows the comparison of the experimental and predicted results using the proposed ANN model for training and testing data. The comparisons between experimental and predicted (ANN) results for training and testing data are listed in Tables 2 and 3, respectively.Table 3
The data that were used for training the network and predicted number of produced neutrons by ANN.
Discharge current (mA)
Voltage (kV)
Produced neutron (×106)
Predicted neutron (×106)
Difference
20
30
0.055164
0.0552
-
3.648
E
-
05
20
40
0.364079
0.3641
-
2.08
E
-
05
20
60
2.846438
2.8464
3.8
E
-
05
20
70
4.468245
4.4682
4.5
E
-
05
30
40
0.358563
0.3667
-
0.0081371
30
50
1.152918
1.2269
-
0.073982
30
60
3.033994
3.0024
0.031594
40
30
0.052405
0.0584
-
0.0059947
40
50
1.130852
1.1437
-
0.012848
40
60
3.072608
3.0582
0.014408
50
30
0.042476
0.0469
-
0.0044241
50
50
1.125336
1.0738
0.051536
50
70
5.95766
6.1068
-
0.14914
60
40
0.292367
0.2832
0.0091667
60
50
1.10327
1.0192
0.08407
60
70
6.729949
6.5857
0.144249
70
30
0.027582
0.029
-
0.0014182
70
50
0.926747
0.9802
-
0.0534529
70
60
3.199484
3.155
0.044484
80
30
0.017652
0.0214
-
0.0037477
80
40
0.248236
0.2421
0.0061358
80
60
3.089157
3.1426
-
0.053443
80
70
7.060931
7.0959
-
0.034969
90
40
0.24272
0.2172
0.0255195
90
50
0.910198
0.9312
-
0.0210019
90
60
3.056059
3.0759
-
0.019841
100
30
0.012136
0.0121
3.597
E
-
05
100
40
0.171007
0.1823
-
0.0112931
100
60
2.923667
2.8965
0.027167
100
70
6.840276
6.8374
0.002876
50
75
7.99871
7.988
0.01071
96
77
11.80499
11.8062
-
0.001207
55
77
11.25336
11.2547
-
0.001342
26
81
7.722893
7.7191
0.003793
75
82
13.51506
13.513
0.002062
48
82
14.61833
14.6195
-
0.001167Comparison of experimental and predicted results for (a) training data and (b) testing data.
(a)
(b)For training results, variance and standard deviation of differences are obtained as1.98
E
-
03 and 0.0445, respectively.For testing result, the variance and standard deviation of differences are obtained as 0.005088 and 0.07134, respectively. From Tables3 and 4 and Figure 5, clearly the predicted produced neutrons by ANN model are close to the experimental results. These results show the applicability of ANN as an accurate and reliable model for the prediction of NPR in IECF device according to the discharge current and voltage. Figure 6 shows the prediction of NPR in the range of current discharge = 40 mA–120 mA and voltage = 50 kV–130 kV.Table 4
The data that were used for testing the network and predicted number of produced neutrons by ANN.
Discharge current (mA)
Voltage (kV)
Produced neutron (×106)
Predicted neutron (×106)
Difference
20
50
1.130852
1.2495
-
0.11865
30
30
0.052957
0.0723
-
0.01934
30
70
4.964717
5.131
-
0.16628
40
40
0.34753
0.3329
0.01463
40
70
5.681843
5.6062
0.075643
50
60
3.166386
3.1047
0.061686
50
40
0.342014
0.3058
0.036214
60
30
0.035856
0.0373
-
0.00144
60
60
3.210517
3.1387
0.071817
70
40
0.28685
0.2629
0.02395
70
70
7.060931
6.955
0.105931
80
50
0.915714
0.9536
-
0.03789
90
30
0.014894
0.0133
0.001594
90
70
7.060931
6.9842
0.076731
100
50
0.904682
0.8952
0.009482Figure 6
Prediction of NPR for current discharge = 40 mA–120 mA and voltage = 50 kV–130 kV.In this range, maximum number of produced neutrons happened inI
=
45 mA and V
=
125 kV. In this situation, IECF can produce 1.855 × 108 neutrons per second. In another optimum point where I
=
72.5 mA and V
=
115 kV, the IECF can produce 1.854 × 108 neutrons per second. These results show significant NPR increase in the IECF device according to increasing voltage. 52% increment in maximum voltage of power supply (82 kV to 125 kV) was found and, in optimum operational discharge current (in this case I
=
45 mA), the neutron produced rate can be increased to 1270% NPR (1.46 × 107 n/s to 1.855 × 108 n/s).


## 5. Conclusion

In this paper, the development of an ANN model for prediction of the neutron yield in IECF devices is presented. The comparison between experimental and predicted results by ANN model shows that there is a good agreement between them with mean relative error (MRE) less than 0.66%. Therefore, the MLP network can be used as an efficient tool to predict the NPR in the IECF device. For developing the model, the input parameters are voltage (kV) and discharge current (mA) and the output is neutron yield (n/s). Results show that in optimum discharge current when maximum voltage of power supply increases 1.5 times, NPR can increase more than 12 times.



---

*Source: 798160-2014-12-16.xml*