Application of Adaptive Neuro-Fuzzy Inference ...

6 downloads 16474 Views 835KB Size Report
your own website. You may .... are the fuzzy sets and pi, qi and ri are the design parameters. (i = 1, 2). ... modifiable parameters to make the ANFIS output match.
Application of Adaptive Neuro-Fuzzy Inference System for Prediction of Neutron Yield of IR-IECF Facility in High Voltages A. Adineh-Vand, M. Torabi, G. H. Roshani, M. Taghipour, S. A. H. Feghhi, M. Rezaei & S. M. Sadati Journal of Fusion Energy ISSN 0164-0313 J Fusion Energ DOI 10.1007/s10894-013-9631-z

1 23

Your article is protected by copyright and all rights are held exclusively by Springer Science +Business Media New York. This e-offprint is for personal use only and shall not be selfarchived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com”.

1 23

Author's personal copy J Fusion Energ DOI 10.1007/s10894-013-9631-z

ORIGINAL RESEARCH

Application of Adaptive Neuro-Fuzzy Inference System for Prediction of Neutron Yield of IR-IECF Facility in High Voltages A. Adineh-Vand • M. Torabi • G. H. Roshani M. Taghipour • S. A. H. Feghhi • M. Rezaei • S. M. Sadati



Ó Springer Science+Business Media New York 2013

Abstract This paper presents a soft computing based artificial intelligent technique, adaptive neuro-fuzzy inference system (ANFIS) to predict the neutron production rate (NPR) of IR-IECF device in wide discharge current and voltage ranges. A hybrid learning algorithm consists of back-propagation and least-squares estimation is used for training the ANFIS model. The performance of the proposed ANFIS model is tested using the experimental data using four performance measures: correlation coefficient, mean absolute error, mean relative error percentage (MRE%) and root mean square error. The obtained results show that the proposed ANFIS model has achieved good agreement with the experimental results. In comparison to the experimental data the proposed ANFIS model has MRE% \1.53 and 2.85 % for training and testing data respectively. Therefore, this model can be used as an efficient tool to predict the NPR in the IR-IECF device. Keywords ANFIS  Prediction  NPR  Inertial electrostatic confinement fusion  IR-IECF device

A. Adineh-Vand Computer Department, Engineering Faculty, Islamic Azad University, Kermanshah, Iran M. Torabi  G. H. Roshani  S. A. H. Feghhi  M. Rezaei  S. M. Sadati Shahid Beheshti University, Tehran, Iran G. H. Roshani (&) Energy Faculty, Kermanshah University of Technology, Kermanshah, Iran e-mail: [email protected] M. Taghipour Department of Biomedical Engineering, Faculty of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran

Introduction Inertial electrostatic confinement (IEC) devices are compact neutron/proton generators that can confine highenergy ions in an electrostatic potential well. Spherical configuration is the most common one that consists of two-concentrate electrodes inside a vacuum chamber. Usually the gridded central electrode connects to a negative high voltage and the outer one grounded. Ions produced by a glow discharge in the inter electrode space gain energy from the applied eclectic field and accelerated to the center inside the hollow cathode and therefore give raise a relatively hot and dense plasma inside the cathode. Because of spherically converging ions, a potential well forms in the plasma core within the cathode, which plays an important role in fusion production [1, 2]. The ions which do not loss (colliding to the grid is the most important loss agent) have a chance to recalculate in the potential well and fuse together. In this case, depending on the working condition beam–beam interaction and beam to target are most dominant process in fusion reaction [3, 4]. As shown in Eqs. (1)–(3), if the IECF device runs by deuterium gas a high flux of neutron and proton can achieve. D þ D ! He3 ð0:82 MeVÞ þ nð2:45 MeVÞ

ð1Þ

D þ D ! Tð1:01 MeVÞ þ pð3:03 MeVÞ

ð2Þ

D þ He3 ! að3:52 MeVÞ þ Pð14:7 MeVÞ

ð3Þ

Depending on the fusion reaction rate and operation mode (continues or pulse mode) IECF device could have variety of application from teaching and material process [5, 6] to high flux neutron sources. Among all of these applications, compact neutron/proton source are most required for cancer detection, cancer treatment boron

123

Author's personal copy J Fusion Energ

neutron capture therapy (BNCT) [7–9] and personal landmine identifications [10, 11]. Research and development on IECF devices, as a neutron/proton source, have been done extendedly and significant results have been obtained [4, 12, 13]. From economical and research point of view, simulation and theoretical method that can predict neutron production rate (NPR) could be very useful for optimized design and performance of this facility. In this study, we predict the NPR of IR-IECF device [14] by ANFIS at high voltage and discharge current.

Fig. 1 Schematic diagram of IR-IECF device

Fig. 2 IR-IEC device with electrodes in the right [14]

123

Experimental Set Up The IR-IECF is a glow discharge based device with a spherical configuration [14, 15]. As shown schematically in the Fig. 1, it consists of a cylindrical vacuum chamber with 60 cm height, 13.5 cm diameter stainless still cathode and 41 cm diameter anode. Figure 2 shows a whole device with electrodes in the right. In the experiments; the device has been adjusted in a given deuterium gas pressure then by raising the cathode voltage, desired plasma has obtained. Finally, for different values of fixed cathode voltage from

Author's personal copy J Fusion Energ

20 to 70 kV the mutual discharge current from 20 to 100 mA has been fixed by controlling the gas pressure. For measuring the NPR a calibrated neutron probe LB6411 has been used.

[18, 19]. ANFIS is based on the first-order Sugeno fuzzy model. An ANFIS consisting of a set of TSK-type fuzzy IF–THEN rules is used to map the system inputs to outputs. To understand the working of ANFIS in simple meaning, for a first order Sugeno fuzzy model, two fuzzy if–then rules are considered:

Modeling Based on ANFIS Rule 1 : If x is A1 and y is B1 ;

then

f1 ¼ p1 x þ q1 y þ r1

ANFIS Architecture

Rule 2 : If x is A2 and y is B2 ;

then

f2 ¼ p2 x þ q2 y þ r2

The adaptive neuro-fuzzy inference system (ANFIS), is one of the examples of neuro-fuzzy systems in which a fuzzy inference system (FIS) is implemented in the framework of adaptive networks. ANFIS is a FIS trained by an artificial neural network (ANN)-learning algorithm [16, 17]. The most important reason for combining FISs with ANNs is to use the learning capability of ANNs

where x and y are the inputs, fi are the outputs, Ai and Bi are the fuzzy sets and pi, qi and ri are the design parameters (i = 1, 2). Figure 3 illustrates the inference method for this Sugeno model. Figure 4 shows a typical ANFIS structure [20, 21]. In this figure, there are two adaptive layers (layers 1 and 4). Layer 1 has modifiable parameters related to the input

Fig. 3 The inference method of the sugeno model

Fig. 4 ANFIS structure

123

Author's personal copy J Fusion Energ

membership function. The parameters in this layer are called premise parameters. The node functions of this layer are given by Jang [18, 19] and Jang et al. [19]: Ai ðxÞ ¼ lAi ðxÞ;

i ¼ 1; 2

ð4Þ

Bi ðyÞ ¼ lBi ðyÞ;

i ¼ 1; 2

ð5Þ

where x and y are the inputs of the model and lAi(x) and lBi(y) are the fuzzy membership functions. A typical membership function is Dsigmf (built-in membership function composed of difference between two sigmoidal membership functions) with four parameters a1, b1, a2 and b2 given by the following equation: Dsigmf ðxÞ ¼

1 1 þ expða1 ðx  b1 ÞÞ 1  1 þ expða2 ðx  b2 ÞÞ

ð6Þ

Every node in layer 2 represents the firing strength of the each rule. The fixed node in layer 3 calculates the ratio of the ith rule’s firing strengths to the sum of all rule’s firing strengths. The output of layer 3 is called normalized firing strengths. The output of layer 2 and layer 3 are given by Eqs. (7) and (8), respectively [20, 21]. Wi ¼ Ai ðxÞ  Bi ðyÞ; i ¼ 1; 2 wi wi ¼ ; i ¼ 1; 2 w1 þ w2

method and least squares method to update parameters. This algorithm consists of two stages: forward pass and backward pass. In the forward pass, the least-squares method is used to identify consequent linear parameters (the parameters in the layer 4) and in the backward pass, the gradient descent method is employed to tune premise nonlinear parameter (the parameters in the layer 1). Detailed information of ANFIS can be found in Jang [18]. Proposed ANFIS Model The purpose of this paper is to investigate the accuracy of an ANFIS model to predict the NPR of IR-IECF device in wide discharge current and voltage ranges. In the developed ANFIS, two variables consisting of voltage and discharge current are selected as input variables to predict the NPR, which is the target variable. Figure 5 shows a simplified overview of the proposed ANFIS model. ANFIS modeling process starts by obtaining a data set (input– output data pairs) and dividing it into training (to find the

ð7Þ ð8Þ

The outputs of the adaptive nodes in layer 4 are given by: O4;i ¼ wi fi ¼ wi ðpi x þ qi y þ ri Þ;

i ¼ 1; 2

ð9Þ

where pi, qi and ri are called consequent parameters. Layer 5 computes the weighted average of the output signals of the layer 4 by the following equation: O5;1 ¼

w 1 f1 þ w 2 f2 w1 þ w2

ð10Þ

Fig. 6 Comparison of the experimental and predicted results for the training data

The training of the ANFIS model is to tune all the modifiable parameters to make the ANFIS output match the training data. ANFIS applies a hybrid-learning algorithm [19, 22], and combines the gradient descent

Fig. 5 The proposed ANFIS model

123

Fig. 7 Comparison of the experimental and predicted results for the testing data

Author's personal copy J Fusion Energ Table 1 The obtained NPR for training data using the proposed ANFIS model in comparison with the experimental Voltage (-kV)

Discharge current (mA)

NPR (n/s) (9106) Experimental

ANFIS

Table 2 The obtained NPR for testing data using the proposed ANFIS model in comparison with the experimental Voltage (-kV)

NPR (n/s) (9106)

Discharge current (mA)

Experimental

ANFIS

30

20

0.055164

0.058982895

50

20

1.130852

1.126840829

40

20

0.364079

0.346367534

30

30

0.052957

0.053360385

60 70

20 20

2.846438 4.468245

2.886784383 4.46734229

70 40

30 40

4.964717 0.34753

5.069799452 0.332687046

40

30

0.358563

0.367689799

70

40

5.681843

5.618292704

50

30

1.152918

1.131691919

60

50

3.166386

3.122481544

60

30

3.033994

2.984878986

40

50

0.342014

0.318160515

30

40

0.052405

0.049922121

30

60

0.035856

0.034549608

50

40

1.130852

1.150745065

60

60

3.210517

3.153693021

60

40

3.072608

3.085758073

40

70

0.28685

0.287435393

30

50

0.042476

0.043210529

70

70

7.060931

7.524078956

50

50

1.125336

1.149567037

50

80

0.915714

0.964636975

70

50

5.95766

5.950154471

30

90

0.014894

0.015490663

40

60

0.292367

0.295518196

70

90

7.060931

6.981909746

50

60

1.10327

1.04687175

50

100

0.904682

0.877110299

70

60

6.729949

6.737933147

30

70

0.027582

0.026099747

50 60

70 70

0.926747 3.199484

0.957147715 3.19684675

Table 3 The obtained errors for the proposed ANFIS model

30

80

0.017652

0.017609357

Error

40

80

0.248236

0.250854839

60

80

3.089157

3.097387806

70

80

7.060931

7.059458937

40

90

0.24272

0.236276293

50

90

0.910198

0.908886184

60

90

3.056059

3.041974107

30

100

0.012136

0.011904104

40

100

0.171007

0.175240295

60

100

2.923667

2.930299327

70

100

6.840276

6.840539717

75

50

7.99871

77

96

11.80499

11.80512942

77

55

11.25336

11.24778307

81 82

26 75

7.722893 13.51506

7.7227766 13.51536227

82

48

14.61833

14.61929399

8.003800703

initial premise parameters for membership functions) and testing data (to validate the trained ANFIS structure) sets. In this study, training and testing data has been obtained by use of IR-IECF experimental data. Total data divided into two sets: 36 samples (about 70 %) are used for training and 15 samples (about 30 %) are used for testing the trained ANFIS model. Testing and training data must be different and are selected randomly from the original data set. MATLAB 7.0.4 software is used for training the ANFIS

Data Train

Test

MAE

0.0102

0.0622

RMSE

0.0172

0.1277

MRE%

1.522

2.845

model. We have tested many different ANFIS architecture with different type of membership functions, number of membership functions and number of epochs to obtain the best ANFIS configuration.

Results and Discussion The optimized ANFIS structure was selected based on the minimal residual error and maximum number of 150 epochs. The proposed ANFIS model being used in this study utilizes Sugeno type FISs. It has 3 Dsigmf membership functions for each input and 9 linear membership functions for the output. The number of MFs assigned to each input variable is chosen by trial and error. It employs 36 linear parameters, 24 nonlinear parameters, and 9 fuzzy rules to predict the output. Figures 6 and 7, Tables 1 and 2 show the comparisons between experimental and predicted results using the proposed ANFIS model, where the correlation coefficient (CC) is calculated by:

123

Author's personal copy J Fusion Energ

"P CC ¼ 1 

N i¼1

ðXi ðExpÞ  Xi ðPredÞÞ2 PN 2 i¼1 ðXi ðExpÞÞ

# ð11Þ

Table 3 shows the obtained errors for the proposed ANFIS model, where the mean relative error percentage (MRE%), the mean absolute error (MAE), and the root mean square error (RMSE) are calculated by: MRE% ¼ 100 

MAE ¼

 N   1X Xi ðExpÞ  Xi ðPredÞ   N i¼1 Xi ðExpÞ

Z 1X jXi ðExpÞ  Xi ðPredÞj N i¼1

"P RMSE ¼

N i¼1

ðXi ðExpÞ  Xi ðPredÞÞ2 N

Fig. 8 The obtained NPR using the proposed ANFIS model (in the range of the training and testing data)

Fig. 9 The obtained NPR using the proposed ANFIS model (outside the range of the training and testing data)

123

ð12Þ

ð13Þ #0:5 ð14Þ

In Eqs. (11)–(14), N is the number of data and ‘X(Exp)’ and ‘X(Pred)’ stand for experimental and predicted (ANFIS) values respectively. According to the Figs. 6 and 7, Tables 1, 2 and 3, it is clear that the predicted NPR using the proposed ANFIS model is close to the experimental. These results show the applicability of ANFIS as an accurate and reliable model for the prediction of NPR of IR-IECF device according to the discharge current and voltage. The obtained NPR using the proposed ANFIS model versus voltage and discharge current in the range of the training and testing data is shown Fig. 8. In this range, the maximum NPR was happened in current = 69.5 mA and voltage = -82 kV. In this situation IECF can produce 20.161 9 106 neutrons/s. Also the minimum NPR was happened in current = 100 mA and voltage = -30 kV. In this situation IECF can produce 0.0119 9 106 neutrons/s. The obtained NPR using the proposed ANFIS model versus voltage and

Author's personal copy J Fusion Energ

discharge current outside the range of the training and testing data is shown Fig. 9. In this range, the maximum NPR was happened in current = 68 mA and voltage = -130 kV. In this situation IECF can produce 43.01 9 106 neutrons/s.

Conclusion In this paper, an ANFIS model is used for prediction of the neutron yield in IECF devices. For developing of the model, the input parameters are voltage and discharge current and the output is NPR. The comparison between the experimental and predicted results using the proposed ANFIS model shows that there is a good agreement between them with MRE% \1.53 and 2.85 % for training and testing data respectively. Therefore, the proposed ANFIS model can be used as an efficient tool to predict the NPR in the IECF device.

References 1. T.A. Thorson, R.D. Durst, R.J. Fonck, L.P. Wainwrite, Phys. Plasmas 4, 4–16 (1997) 2. R.L. Hirsh, Phys. Fluids 11, 2486 (1968) 3. H. George et al., IEEE Trans. Plasma Sci. 25(4), 733–739 (1997) 4. T. Takamatsu et al., Nucl. Fusion 46, 142–148 (2006) 5. S. Semsari et al., The effect of high temperature He? implantation on polycrystalline tungsten in IR-IECF. J. Fusion Energ. 31(4), 389–395 (2012) 6. S. Semsari et al., Comparison of high energy He? and D? irradiation impact on tungsten surface in the IR-IECF device. J. Fusion Energ. 32(1), 142–149 (2013) 7. G.L. Kulcinski, Non-electric applications of fusion energy—an important precursor to commercial electric power. Fusion Technol. 34, 477 (1998)

8. R.P. Ashley, The UW IEC Team, in Experimental Progress in 2003 of the UW IEC Facility. 6th US–Japan IEC Workshop, Tokyo, Japan, 20–21 October (2003) 9. R.P. Ashley, G.L. Kulcinski, J.F. Santarius, S. KrupakarMurali, G. Piefer, R. Radel, Steady state D He3 proton production in an IEC fusion device. Fusion Technol. 39, 546 (2001) 10. Application of Nuclear Techniques to Anti-Personnel Landmines Identification, IAEA/PS/RC-799, IAEA (2000) 11. K. Yoshikawa et al., in Research and Development on a Compact Discharge-Driven D–D Fusion Neutron Source for Explosive Detection. The 2nd Joint International Conference on Sustainable Energy and Environment, Bangkok, Thailand, 21–23 November (2006) 12. K. Nozak et al., in International Symposium on Discharges and Electrical Insulation in Vacuum, pp. 650–653 (2006) 13. J.F. Santarius et al., Fusion Sci. Technol. 47, 1238–1243 (2005) 14. V. Damideh et al., Experimental study of the Iranian inertial electrostatic confinement fusion device as a continuous neutron generator. J. Fusion Energ. 31(2), 109–111 (2012) 15. E. Haji Ebrahimi et al. The influence of cathode voltage and discharge current on neutron production rate of inertial electrostatic confinement fusion (IR-IECF). J Fusion Energ. 32(1), 62–65 (2012) 16. S. Shao, Fuzzy self-organizing controller and its application for dynamic processes. Fuzzy Sets Syst. 26, 151–164 (1998) 17. H. Takagi, in Application of Neural Networks and Fuzzy Logic to Consumer Products. Proceedings of the International Conference on Industrial Fuzzy Electronics, Control, Instrumentation, and Automation, vol. 3, San Diego, November 2000, pp. 1629–1639 18. J.S.R. Jang, ANFIS: adaptive-network-based fuzzy inference systems. IEEE Trans. Syst. Man Cybern. 23, 665–685 (1993) 19. J.S.R. Jang, C.T. Sun, E. Mizutani, A Computational Approach to Learning and Machine Intelligence Neuro-Fuzzy and Soft Computing (Prentice Hall, US edition, 1997) 20. S. Shahbudin, A. Hussain, A. El-Shafie, N.M. Tahir, S.A Samad, Adaptive-neuro fuzzy inference system for human posture classification using a simplified shock graph. Vis. Inform. Bridging Res. Pract. 5857, 585–595 (2009) 21. T. Culliere, A. Titli, J. Corrieu, in Neuro-Fuzzy Modeling of Nonlinear Systems for Control Purposes. Proceedings of the IEEE International Conference on Fuzzy Systems, Yokohama, pp. 2009–2016 (1995) 22. J.S.R. Jang, C.T. Sun, Neuro-fuzzy modeling and control. IEEE Spec. Issue Fuzzy Logic Eng. Appl. 83, 378–406 (1995)

123