High-precision and wide-range real-time neutron flux monitor system through multipoint linear calibration

High-precision and wide-range real-time neutron flux monitor system through multipoint linear calibration 1 2 first-author https://orcid.org/0000-0002-4221-0298 1 2 https://orcid.org/0000-0003-4582-8314 1 2 https://orcid.org/0000-0002-8809-1503 1 2 4 1 2 5 1 2 corresp Corresponding author, xuxf@mail.ustc.edu.cn Corresponding author, xuxf@mail.ustc.edu.cn xuxf@mail.ustc.edu.cn 3 3 3 1 2 State Key Laboratory of Particle Detection and Electronics, University of Science and Technology of China, Hefei 230026, China State Key Laboratory of Particle Detection and Electronics, University of Science and Technology of China, Hefei 230026, China Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China Center for Fusion Science of Southwestern Institute of Physics, Chengdu 610041, China Center for Fusion Science of Southwestern Institute of Physics, Chengdu 610041, China National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China College of Electronic Engineering, National University of Defense Technology, Hefei 230037, China College of Electronic Engineering, National University of Defense Technology, Hefei 230037, China The neutron flux monitor (NFM) system is an important diagnostic subsystem introduced by large nuclear fusion devices such as international thermonuclear experimental reactor (ITER), Japan torus-60, tokamak fusion test reactor, and HL-2A. Neutron fluxes can provide real-time parameters for nuclear fusion, including neutron source intensity and fusion power. Corresponding to different nuclear reaction periods, neutron fluxes span over seven decades, thereby requiring electronic devices to operate in counting and Campbelling modes simultaneously. Therefore, it is crucial to design a real-time NFM system to encompass such a wide dynamic range. In this study, a high-precision NFM system with a wide measurement range of neutron flux is implemented using real-time multipoint linear calibration. It can automatically switch between counting and Campbelling modes with variations in the neutron flux. We established a testing platform to verify the feasibility of the NFM system, which can output the simulated neutron signal using an arbitrary waveform generator. Meanwhile, the accurate calibration interval of the Campbelling mode is defined well. Based on the abovementioned design, the system satisfies the requirements, offering a dynamic range of 10 8 cps, temporal resolution of 1 ms, and maximal relative error of 4% measured at the signal–noise ratio of 15.8 dB. Additionally, the NFM system is verified in a field experiment involving HL-2A, and the measured neutron flux is consistent with the results. Multipoint linear calibration Real time Simulated neutron signal Neutron flux monitoring ITER HL-2A I. Introduction 1 ITER is a core experimental platform constructed to investigate the controlled nuclear fusion reaction and provide guidance for the optimized design and operation of fusion power plants [ 1 1 ]. To comprehensively monitor the entire nuclear fusion reaction process, it is vital to develop reliable diagnostic systems to collect all the necessary plasma parameters to protect the machine, improve plasma performance, and promote physics studies [ 2 2 – 7 7 2–7 ]. Among them, the NFM system is one of the major diagnostic systems used for measuring the total neutron flux, and it has been used since the initial design phase of the ITER [ 8 8 – 10 10 8–10 ]. Neutron flux measurement is performed to evaluate neutron emissivity and is central for understanding the behavior of nuclear fusion reactions[ 11 11 ]. Generally, the NFM system comprises two parts, i.e., fission chamber detectors and an electronic processing module [ 12 12 ]. Because many gamma signals are generated during the reaction, it is important to select fission chamber detectors with an high neutron gama (n/γ) suppression ratio. Meanwhile, the traditional NFM system is based on analog circuits with simple discrimination algorithms [ 13 13 , 14 14 ]. Consequently, complex offline analysis must be performed to ensure the accuracy of neutron multiplicity counting, and the raw data used to require a considerable amount of storage space [ 15 15 - 17 17 15-17 ]. Meanwhile, the neutron flux of the ITER exhibits a wide range from 10 14 to approximately 10 21 n/s [ 18 18 ], rendering it difficult to satisfy such a wide requirement range through only the single counting or Campbelling mode in electronics. In this study, we designed a high-precision and wide-range real-time NFM system through multipoint linear calibration. The digital pulse shape discrimination (DPSD) module, which is based on a high-speed analog–to-digital converter (ADC) and field-programmable gate array (FPGA), can be applied to organic scintillators and liquid scintillation detectors in counting mode. Owing to the high degree integration of digital circuits, we can easily implement real-time discriminating algorithms. To encompass the wide dynamic range of neutron fluxes, the counting and Campbelling modes (known as "fluctuation mode" or "mean square voltage mode") were integrated into the system. Real-time multipoint linear calibration can be used to solve the crossover between counting and Campbelling modes and improve the accuracy of neutron fluxes in the Campbelling mode. Using multipoint linear calibration in the overlapping interval of the Campbelling and counting modes, we can obtain the exact conversion coefficients between the Campbelling value and the real neutron flux. Subsequently, the measurement dynamic range of the neutron flux can be extended to 10 8 counts per second (cps). Compared with the DSP unit prototype of Toshiba [ 19 19 ], the NFM system has a faster time response of 0.2 ms and a temporal resolution of 1 ms. In particular, the simulated neutron signal processed by an arbitrary waveform generator, which is consistent with the characteristics of the output waveform of the fission chamber detector, was performed to measure the accuracy of the system. Using simulated neutron signals overcomes the shortage of actual neutron sources. Furthermore, the simulated neutron signal has universal applicability; hence, it can be used widely in the evaluation of neutron flux measurement systems. Based on the design above, the noise margin of the system was obtained, and the relative error of the neutron counting rate did not exceed 4% because the accurate calibration intervals of the Campbelling mode was analyzed well; the measured counting rate was much higher than the 10% accuracy required for the neutron flux measurement [ 20 20 ]. II. Design Principles of NFM System 1 The NFM system primarily comprises four modules, i.e., fission chambers (FCs), charge sensitive preamplifiers (PAs), main electronics (MA), and monitoring software of the upper computer. The overall block diagram is shown in Fig. 1 Fig. 1 , where PXI represents peripheral component interconnection extensions for instrumentation. A. Fission Chamber 2 The fission chamber is a gas counter containing 235 U [ 12 12 , 21 21 - 23 23 21-23 ] mainly composed of fissile material coating and a gas ionization chamber. When the neutron enters the ionization chamber, a nuclear reaction occurs with fissile material UO 2 on the cathode, and hundreds of MeV of energy is released. Subsequently, the generated fission fragments further react with the gas when passing through the ionization chamber. Finally, the generated ions are collected by the electric fields of the positive and negative electrodes to form a charge pulse signal. Correspondingly, the energy released by the neutron are orders of magnitude larger than the γ-ray and hard X-ray; therefore, pulse shape discrimination is unnecessary for 235 U FCs. In field experiments involving the HL-2A Tokamak device, the amplitudes of the gamma and noise were approximately tens of millivolts, but the amplitudes of the neutron pulse signal exceeded 200 mV. Therefore, we set approximate discrimination threshold, and only signals with a pulse amplitude greater than the threshold were recognized as neutron pulses. B. Preamplifier Design 2 The charge signal output from the fission chamber is extremely weak and disappears easily owing to the high output impedance. Meanwhile, some electromagnetic interference and noise effects were inevitable in the experimental environment. The pulse signal cannot be distinguished from the noise if the long-line transmission (150 m from FC to MA) is performed directly. Therefore, the charge pulse signal output from the fission chamber must be amplified and processed with a signal conditioning circuit to improve the anti-interference ability. The charge-sensitive preamplifier is mainly composed of four-stage amplifiers (see Fig. 2 Fig. 2 ). The negative pulse output from the FC is first input to the first stage (A 1 ) for integral amplification, which involves a charge-sensitive amplifier with high input impedance. Subsequently, the signal passes through the second and third stages (A 2 and A 3 , respectively) of the high-gain amplifier. Finally, the signal passes through the fourth stage (A 4 ) of the low output impedance inverted amplifier and then outputs a positive pulse. Considering that the range of input signals for the MA should not exceed 5 V, the amplitude of a single neutron pulse should be limited to less than 300 mV after being processed by the preamplifiers. C. High-Speed MA Design 2 In recent years, with the increase in ADC sampling rate, DPSD technology has piqued significant research interest owing to its flexible processing and strong scalability for digital signals. In this study, we designed MA capable of extracting valuable information from analog signals. The ADC sampling rate was as high as 250 million samples per second with a 12-bit width. The designed MA is crucial for processing the digital signal of the entire system and controlling the operation of each module. The system can support three boards to collect data in parallel and then upload the data to the host computer through the PXI bus interface. The overall design diagram is shown in Fig. 2 Fig. 2 , where ACC represents alternating current coupling, PGA represents programmable gain amplifier, and DCI represents digital Campbelling integrator. An FPGA was adopted to calculate the Campbelling integral value of the sampling data from the ADC, and the calculation was based on the Campbelling theorem [ 24 24 ]. The output voltage of the detector can be described based on the current impact Qδ ( t ), where Q is the strength of the current impact, i.e., the amount of charge carried. Considering that the detector is a linear time-invariant system that can be depicted by an impulse response function h ( t ), the impulse in the interval from τ to τ+ d τ will contribute to the output voltage at time t . Here, λ represents the average counting rates of the impulse signal, and the events of nuclear radiation induced by the detector obey the principle of Poisson distribution over time; therefore, the average counting rate is λ d τ for the interval from τ to τ+ d τ . The average contribution of the output voltage at time t is as follows: <math display="block" id="M1"><mrow><mover accent="true"><mrow><mtext>d</mtext><mi>v</mi></mrow><mo stretchy="true">¯</mo></mover><mo>=</mo><mi>λ</mi><mtext>d</mtext><mi>τ</mi><mi>Q</mi><mi>h</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><mi>u</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>,</mo></mrow></math> where u ( t ) is the unit step function. The root mean square deviation of the output voltage at time t is attributed entirely to the randomness of the number of pulses n . According to the Poisson mean square distribution equation, the root mean square deviation can be obtained, as follows: <math display="block" id="M2"><mrow><mi>σ</mi><msub><mrow><mo stretchy="false">(</mo><mtext>d</mtext><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><msqrt><mrow><mi>λ</mi><mtext>d</mtext><mi>τ</mi></mrow></msqrt><mi>Q</mi><mi>h</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><mi>u</mi><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>.</mo></mrow></math> As each small interval d τ is independent of each other, the mean square deviation of the output voltage generated by the entire impulse sequence at time t is as follows: <math display="block" id="M3"><mrow><msup><mi>σ</mi><mn>2</mn></msup><msub><mrow><mo stretchy="false">(</mo><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><mi>λ</mi><msup><mi>Q</mi><mn>2</mn></msup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mrow><msup><mi>h</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><msup><mi>u</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo>−</mo><mi>τ</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>τ</mi></mrow></mrow></mstyle><mo>.</mo></mrow></math> Subsequently, <math display="block" id="M4"><mrow><msup><mi>σ</mi><mn>2</mn></msup><msub><mrow><mo stretchy="false">(</mo><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><mi>λ</mi><msup><mi>Q</mi><mn>2</mn></msup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup><mrow><msup><mi>h</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>t</mi></mrow></mrow></mstyle><mo>.</mo></mrow></math> In fact, the charge amount of each current pulse output by the detector exhibits a specific fluctuation that is related to the incident particle’s energy spectrum and the statistical fluctuation of the detector. Herein, <math id="M5"><mrow><mover accent="true"><mrow><msup><mi>Q</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover></mrow></math> is used to represent the average squared amount of charge. Similarly, we obtained the expression of the Campbell theorem as follows: <math display="block" id="M6"><mrow><msup><mi>σ</mi><mn>2</mn></msup><msub><mrow><mo stretchy="false">(</mo><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><mi>λ</mi><mover accent="true"><mrow><msup><mi>Q</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup><mrow><msup><mi>h</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>t</mi></mrow></mrow></mstyle><mo>.</mo></mrow></math> The variance is calculated as <math display="block" id="M7"><mrow><msup><mi>σ</mi><mn>2</mn></msup><msub><mrow><mo stretchy="false">(</mo><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><mover accent="true"><mrow><msup><mi>v</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover><mo>−</mo><msup><mover accent="true"><mi>v</mi><mo>¯</mo></mover><mn>2</mn></msup></mrow></math> As the system uses AC coupling, <math id="M8"><mrow><mover accent="true"><mi>v</mi><mo>¯</mo></mover><mo>=</mo><mn>0</mn></mrow></math> . Therefore, the variance evolves into <math display="block" id="M9"><mrow><mover accent="true"><mrow><msup><mi>v</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover><mo>=</mo><msup><mi>σ</mi><mn>2</mn></msup><msub><mrow><mo stretchy="false">(</mo><mi>v</mi><mo stretchy="false">)</mo></mrow><mi>n</mi></msub><mo>=</mo><mi>λ</mi><mover accent="true"><mrow><msup><mi>Q</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup><mrow><msup><mi>h</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>t</mi></mrow></mrow></mstyle><mo>.</mo></mrow></math> According to formula (7), <math id="M10"><mrow><mover accent="true"><mrow><msup><mi>v</mi><mn>2</mn></msup></mrow><mo stretchy="true">¯</mo></mover></mrow></math> is proportional to the pulse count rate. Therefore, by calculating the square of the voltage value sampled by the high-speed ADC and performing a linear calibration [ 25 25 , 26 26 ], the measurement range of the neutron flux can be widened to 10 8 cps. Considering the high-speed parallel operation ability of the FPGA, we integrated the main algorithms to the FPGA to process the signal in real time. The output signal of the preamplifier was transmitted to the analog signal processing module of the MA for conditioning. This module is mainly used for 1) AC coupling, which can remove the DC signal (required by Campbell’s theorem); 2) gain controlling, which enables the signal to reach the ADC sampling requirement via FPGA control; 3) translating the single-ended signal into a differential signal, which enables the ADC to perform sampling. Next, the signal was input to the ADC digital sampling module of the MA; hence, the analog signal was converted into a digital signal (250M). Finally, the digital signal was processed by the FPGA module of the MA, which exhibits high-speed computing capabilities owing to the parallel computing afforded by the FPGA. Limited by the uploading speed of the PXI bus, we ultimately focused on only the neutron flux for the experimental results. Therefore, the raw data sampled by the ADC was first processed using the FPGA; subsequently, the as-obtained result of counting rate and Campbelling value without calibrating were uploaded to the host PC through the PXI bus. This module is crucial for 1) performing n/γ discrimination according to the pulse amplitude, pulse width, and other parameters; 2) calculating counting values and the integral value of Campbelling; 3) controlling the operation of the system; and 4) packaging the calculation result of the FPGA and uploading them to the PXI bus. D. Software Design 2 The data processed by the FPGA can only be stored in RAM resources; hence, the contents will disappear once the power is turned off. Therefore, it is necessary to upload the calculated data to the host computer through the PXI bus. Subsequently, we calibrated the Campbelling value in real time based on the multipoint linear strategy. The overall design of the host computer software is shown in Fig. 3 Fig. 3 . When the neutron flux was in the overlapping range of 10 4 to 10 5 cps [ 25 25 , 27 27 ], the counting and Campbelling values were both accurate. Furthermore, the rising process of the neutron flux must be sufficiently wide to contain the calibration interval in the actual fusion reaction such that the Campbelling values can be calibrated using the data of the overlapping region. Subsequently, the correlation between the actual neutron flux and the Campbelling values can be obtained. Finally, the actual neutron flux can be calculated directly from the calibration result. As shown in Fig. 3 Fig. 3 , the data uploaded by the FPGA is mainly the counting and Campbelling values per time interval. The calibration module continuously monitors the counting values uploaded from the FPGA. To reduce the error caused by the calibration, the accurate calibration intervals of Campbelling and counting were divided equally into 10 subintervals. Only a set of data containing a counting value and a Campbelling value for each subinterval was selected for calibration in the same experiment. The 10 as-obtained sets of data were linearly fitted via the least-squares method to obtain values of linear correlation coefficients, intercepts, and slopes. Generally, there exists only one rising edge containing calibration interval in the same experiment. Therefore, calibration was usually performed only once in each discharge experiment; subsequently, the neutron flux was obtained using the calibration coefficients and Campbell values. III. Simulated Signals of Neutron Pulse 1 Owing to the harsh conditions of using the neutron source and the complex operating environment, it was challenging to directly obtain the actual neutron signal pulses that represent the actual neutron flux of fusion reaction in the laboratory. Therefore, we developed a system for simulating the output signal of the fission chamber. First, the data for simulating a neutron signal was generated using MATLAB; subsequently, a waveform was output using an arbitrary waveform generator (Keysight 81160A). A. Principle 2 The single neutron signal of the fission chamber satisfies formula (8), <math display="block" id="M11"><mrow><msub><mi>V</mi><mrow><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mrow></msub><mo>=</mo><mfrac><mrow><mi>Q</mi><mi>R</mi></mrow><mi>τ</mi></mfrac><mo stretchy="false">(</mo><mo stretchy="false">(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>w</mi><mi>t</mi></mrow></msup><mo stretchy="false">)</mo><mi>u</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>-</mo><mo stretchy="false">(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>w</mi><mo stretchy="false">(</mo><mi>t</mi><mo>-</mo><mi>τ</mi><mo stretchy="false">)</mo></mrow></msup><mo stretchy="false">)</mo><mi>u</mi><mo stretchy="false">(</mo><mi>t</mi><mo>-</mo><mi>τ</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>,</mo></mrow></math> where Q is the amount of charge released from the fission chamber [ 28 28 ] (set to Q = 0.5 pC), R is the equivalent input resistance of the preamplifier (set to R = 5 kΩ), τ is the electronic collection time of a fission reaction in the fission chamber (set to τ = 400 ns), w = 1 / RC , and C = 200 pF. The resulting single pulse neutron waveform was similar to the waveform shown in Fig. 5 Fig. 5 . Meanwhile, the waveform of simulated single neutron pulse is based on the measured single neutron pulse on HL-2A ( Fig. 8 Fig. 8 ), and the waveforms of the two agreed well with each other. The statistical characteristics of the fission chamber’s output signal satisfy the following: 1) In the nuclear physics experiment, the number of particles recorded by the detector within a certain period of time is completely independent. Therefore, the number of particles recorded by the detector in a certain period of time can be regarded as a Bernoulli experiment. When the probability of event A is sufficiently small ( p<< 1) and the number of experiments n is large and np tends to a finite value λ , then the binomial distribution is approximately a Poisson distribution. 2) The amplitudes of the fission chamber pulses are proportional to the energy absorbed by the incident particles, and random fluctuations exist in the absorbed energy, which can be represent by the Gaussian distribution in a general simulation. The standard deviation of the distribution σ= R × A/ 2.33, where R represents the detector’s intrinsic energy resolution, and A represents the average value of the pulse amplitude. 3) Owing to the interference of environmental factors such as detectors and electronic devices, noise simulated by the Gaussian distribution noise will be introduced during the measurement of the neutron signal. 4) The PA is treated as a linear time-invariant system that satisfies the principle of superposition. The superposition of neutron waveforms was processed by superimposing the amplitudes linearly. According to the principles above, the simulation of the neutron pulse waveform is formed by the superposition of three parts, i.e., a single neutron pulse waveform whose amplitude satisfies the Gaussian distribution, a Poisson distribution sequence, and a noise signal satisfying a random distribution ( Fig. 4 Fig. 4 ). B. Results 2 Based on the method above, the simulated neutron signals with neutron fluxes ranging from 10 5 to 10 8 cps can be obtained ( Fig. 5 Fig. 5 ). In MATLAB, we designed a graphical user interface that can adjust the parameters to simulate the waveforms for all neutron flux ranges. IV. Results and Discussion 1 Under laboratory conditions, it is difficult to obtain the neutron source that satisfies the counting rate requirement for the acquisition system. Hence, we simulated neutron signals to verify the feasibility of the NFM system in advance. Simulated waveforms with counting rates of 10 3 –10 8 cps were generated using an arbitrary waveform generator (Keysight 81160A) and then input to the MA for neutron flux measurement. We obtained the relative error by comparing the measured and theoretical results. A. Union Test Using Simulated Signals 2 Based on the characteristics of the neutron pulse, a pulse waveform with counting rates of 10 3 to 10 8 cps was generated using a Keysight 81160A arbitrary waveform generator. The width of a single neutron pulse was from 200 to 400 ns, the pulse amplitude A satisfied the Gaussian distribution A ~ N (μ, σ 2 ), where ( μ = 250 mV, σ = 10 mV), and the Gaussian white noise V pp = 150 mV (99.7%). The counting and Campbelling values were uploaded to the PC and compared with the theoretical values. The obtained relative error is shown in Fig. 6 Fig. 6 . As depicted in Fig. 6 Fig. 6 , when the counting rate of the neutron was less than 4.04 × 10 5 cps, the relative error of the counting mode was less than 2.23%. When the counting rate of the neutron was greater than 2.2 × 10 4 cps, the relative error of Campbelling mode was less than 3.14%. Correspondingly, the relative error was independent of the width of the single neutron signal. In particular, the relative error with the neutron counting rate from 2 × 10 4 to 4 × 10 5 cps was less than 2.28% in the Campbelling mode and 2.23% in the counting mode; therefore, this overlapping interval was regarded as the calibration interval for the Campbelling mode in the following experiments, and the counting rate of 4 × 10 5 cps was set as the switching point for the counting and Campbelling modes. Furthermore, the simulation verified that the relative error of the NFM system was less than 4%, indicating that the NFM system significantly improved the accuracy of the neutron flux measurement. B. Experimental Results on HL-2A 2 A field experiment was performed on this system using the HL-2A Tokamak device of the Southwest Institute of Physics, where high-purity 235 U fission chambers were used. The plateau length was from 300 to 800 V, slope of the plateau was ≤ 0.9%, sensitivity was 1±0.15 cps/nυ, and insulation resistance was 2 × 10 9 Ω [ 29 29 ]. The completed system is shown in Fig. 7 Fig. 7 , including the HL-2A Tokamak device and the NFM system. In addition, a hardware trigger module was introduced to the system to compare with other measurement results of HL-2A; this trigger module can respond to the trigger signal of HL-2A and then triggers the board of the MA to begin operation through the PXI bus (clock zero). When operating in the counting mode, the noise threshold must be set appropriately to remove the effects of noise and γ rays. First, a graph with a neutron pulse was captured in the oscilloscope mode (another function of the system), as shown in Fig. 8 Fig. 8 ; the noise amplitude of the system was ±50 mV. To improve the system accuracy, we set the noise threshold above 50 mV to remove the effect of noise on the counting rate of the neutron. Furthermore, the figure that the rising time of the single neutron pulse was approximately 80 ns, the pulse width was approximately 300 ns, and the amplitude was approximately 300 mV, indicating a typical neutron pulse signal. In shot 31231, the MA measured the number of neutron pulse signals output from the fission chamber; the specific results are shown in Fig. 8 Fig. 8 . The counting rate of the neutron pulse signals measured by this system increased from (502 ms, 0 cps) to (511 ms, 40000 cps) and then stabilized for approximately 489 ms before finally decreasing from (1000 ms, 50000 cps) to (1037 ms, 0 cps). In general, the curve of the counting rate of the neutron pulse signals shows the same trend as that of the power of the neutral beam injection (PNBI). In detail, the PNBI began increasing at 499.14 ms and then reduced to 0.57 kW at 1003.41 ms. The duration of the neutron pulse signals measured by the NFM system was approximately 535 ms, in accordance with the duration of neutral beam injection (NBI) for 504.27 ms as the neutrons were only generated during NBI. Furthermore, the measured duration of the neutron pulse signals was longer owing to the relaxation time of the nuclear reaction. However, the neutron flux was insufficient to reach the switching point of the Campbelling and counting modes. Based on the experimental results above, we confirmed that the system is reliable and applicable for measuring neutron fluxes. V. Conclusion 1 In this study, we designed a high-precision and wide-range real-time NFM system and proposed a multipoint linear calibration algorithm. Both the counting and Campbelling modes were employed for the MA, and the measured dynamic range of the neutron flux was extended to 10 8 cps with the maximal relative error of 4%. In addition, the characteristics of the neutron pulse signals were described, and neutron signals superimposed with noise were simulated. Using the simulated neutron signals solved the deficiency of the actual neutron source, and the simulation afforded universal applicability; hence, it can be used widely to detect neutron fluxes. Additionally, the NFM system was validated through a field experiment involving HL-2A, and the measured neutron flux was consistent with the results. Despite the achievements above, the system requires further verifications based on field experiments involving the ITER. Furthermore, it is still challenging to resist noise and suppress interference under strong electromagnetic fields when the cable between the FC and PA is increased to 50 m long. References J.W. Yang , Q.W. Yang , G.S. Xiao et al ., Fusion neutron flux monitor for ITER . Plasma Sci. Technol . 10 : 141 ( 2008 ). https://doi.org/10.1088/1009-0630/10/2/01 https://doi.org/10.1088/1009-0630/10/2/01 ; Fusion neutron flux monitor for ITER A. J. H. Donné , A. E. Costley , R. Barnsley et al ., Chapter 7: Diagnostics . Nucl. Fusion 47 : S337 ( 2007 ). https://doi.org/10.1088/0029-5515/47/6/S07 https://doi.org/10.1088/0029-5515/47/6/S07 ; Chapter 7: Diagnostics L. Bertalot , V. Krasilnikov , L. Core et al ., Present Status of ITER Neutron Diagnostics Development . J. Fusion Energy ( 2019 ). https://doi.org/10.1007/s10894-019-00220-w https://doi.org/10.1007/s10894-019-00220-w ; Present Status of ITER Neutron Diagnostics Development G. Vayakis , S. Arshad , D. Delhom et al ., Development of the ITER magnetic diagnostic set and specification . Rev. Sci. Instrum . 83 : ( 2012 ). https://doi.org/10.1063/1.4732077 https://doi.org/10.1063/1.4732077 ; Development of the ITER magnetic diagnostic set and specification C.J. Tang , C.J. Xiao , Y.G. Yao et al ., Development of neutron depth profiling system at CARR . Nucl. Tech . 42 : 040402 ( 2019 ). https://doi.org/10.11889/j.0253-3219.2019.hjs.42.040402 https://doi.org/10.11889/j.0253-3219.2019.hjs.42.040402 (in Chinese) Development of neutron depth profiling system at CARR Z.J. Liu , Y.W. Yang , L. Zheng et al ., Measurement and analysis of the 232Th(n,2n) reaction rate in a polyethylene shell with DT neutrons . Nucl. Tech . 41 : 060502 ( 2018 ). https://doi.org/10.11889/j.0253-3219.2018.hjs.41.060502 https://doi.org/10.11889/j.0253-3219.2018.hjs.41.060502 (in Chinese) Measurement and analysis of the 232Th(n,2n) reaction rate in a polyethylene shell with DT neutrons B. W. Zheng , C. Y. Jiang , Z. H. Liu et al ., Correction and verification of HL-2A Tokamak Bonner sphere spectrometer in monoenergetic neutron fields from 100 keV to 5 MeV . Nucl. Sci. Tech . 30 : 159 ( 2019 ). https://doi.org/10.1007/s41365-019-0689-9 https://doi.org/10.1007/s41365-019-0689-9 Correction and verification of HL-2A Tokamak Bonner sphere spectrometer in monoenergetic neutron fields from 100 keV to 5 MeV M. Sasao , A. V. Krasilnikov , T. Nishitani et al ., Overview of neutron and confined/escaping alpha diagnostics planned for ITER . Plasma Phys. Control. Fusion 46 : S107 ( 2004 ). https://doi.org/10.1088/0741-3335/46/7/S08 https://doi.org/10.1088/0741-3335/46/7/S08 ; Overview of neutron and confined/escaping alpha diagnostics planned for ITER C. W. Barnes and A. L. Roquemore ., Neutron source strength monitors for ITER . Rev. Sci. Instrum . 68 : 573 ( 1997 ). https://doi.org/10.1063/1.1147656 https://doi.org/10.1063/1.1147656 ; Neutron source strength monitors for ITER H.K. Zhang , Y.T. Zhao , A. Liu et al ., On-line monitoring of the thermal neutron sensitivity of Co-SPND detector . Nucl. Tech . 43 : 040003 ( 2020 ). https://doi.org/10.11889/j.0253-3219.2020.hjs.43.040003 https://doi.org/10.11889/j.0253-3219.2020.hjs.43.040003 (in Chinese) On-line monitoring of the thermal neutron sensitivity of Co-SPND detector B. W. Zheng , W. Zhang , T. Y. Wu et al ., Development of the real-time double-ring fusion neutron time-of-flight spectrometer system at HL-2M . Nucl. Sci. Tech . 30 : 175 ( 2019 ). https://doi.org/10.1007/s41365-019-0702-3 https://doi.org/10.1007/s41365-019-0702-3 Development of the real-time double-ring fusion neutron time-of-flight spectrometer system at HL-2M G. L. Yuan , Q. W. Yang , J. W. Yang et al ., Fusion Neutron Flux Detector for the ITER . Plasma Sci. Technol . 16 : 168 ( 2014 ). https://doi.org/10.1088/1009-0630/16/2/14 https://doi.org/10.1088/1009-0630/16/2/14 ; Fusion Neutron Flux Detector for the ITER M. J. Joyce , M. D. Aspinall , F. D. Cave et al ., The design, build and test of a digital analyzer for mixed radiation fields . IEEE Trans. Nucl. Sci . 57 : 2625 ( 2010 ). https://doi.org/10.1109/TNS.2010.2044245 https://doi.org/10.1109/TNS.2010.2044245 ; The design, build and test of a digital analyzer for mixed radiation fields K. A. A. Gamage , M. J. Joyce , and J. C. Adams ., Combined digital imaging of mixed-field radioactivity with a single detector . Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip . 635 : 74 ( 2011 ). https://doi.org/10.1016/J.NIMA.2011.01.033 https://doi.org/10.1016/J.NIMA.2011.01.033 ; Combined digital imaging of mixed-field radioactivity with a single detector S. Normand , B. Mouanda , S. Haan et al ., Discrimination methods between neutron and gamma rays for boron loaded plastic scintillators . Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip . 484 : 342 ( 2002 ). https://doi.org/10.1016/S0168-9002(01)02016-2 https://doi.org/10.1016/S0168-9002(01)02016-2 ; Discrimination methods between neutron and gamma rays for boron loaded plastic scintillators Y.-K. Kim , S. K. Lee , B. H. Kang et al ., Performance improvement of neutron flux monitor at KSTAR . J. Instrum . 7 : C06013 ( 2012 ). https://doi.org/10.1088/1748-0221/7/06/C06013 https://doi.org/10.1088/1748-0221/7/06/C06013 Performance improvement of neutron flux monitor at KSTAR J. Wu , X. Zhou , C. Yuan et al ., A real-time online data acquisition system for Dragon-I linear induction accelerator . Nucl. Sci. Tech . 28 : 1 ( 2017 ). https://doi.org/10.1007/s41365-017-0182-2 https://doi.org/10.1007/s41365-017-0182-2 A real-time online data acquisition system for Dragon-I linear induction accelerator L. Bertalot , R. Barnsley , M. F. Direz et al ., Fusion neutron diagnostics on ITER tokamak . J. Instrum . 7 : C04012 ( 2012 ). https://doi.org/10.1088/1748-0221/7/04/C04012 https://doi.org/10.1088/1748-0221/7/04/C04012 ; Fusion neutron diagnostics on ITER tokamak M. Isobe , K. Ogawa , H. Miyake et al ., Wide dynamic range neutron flux monitor having fast time response for the Large Helical Device . Rev. Sci. Instrum . 85 : 11E114 ( 2014 ). https://doi.org/10.1063/1.4891049 https://doi.org/10.1063/1.4891049 ; Wide dynamic range neutron flux monitor having fast time response for the Large Helical Device Y. Kaschuck , A. Krasilnikov , A. Alekseyev et al ., Neutron flux monitoring system for ITER-FEAT (abstract) . Rev. Sci. Instrum . 72 : 823 ( 2001 ). https://doi.org/10.1063/1.1323487 https://doi.org/10.1063/1.1323487 ; Neutron flux monitoring system for ITER-FEAT (abstract) T. Nishitani , S. Kasai , L. C. Johnson et al ., Neutron monitor using microfission chambers for the International Thermonuclear Experimental Reactor . Rev. Sci. Instrum . 70 : 1141 ( 1999 ). https://doi.org/10.1063/1.1149297 https://doi.org/10.1063/1.1149297 ; Neutron monitor using microfission chambers for the International Thermonuclear Experimental Reactor M. Ishikawa , T. Kondoh , T. Nishitani et al ., Design of microfission chamber for ITER operations . Rev. Sci. Instrum . 79 : 1 ( 2008 ). https://doi.org/10.1063/1.2969286 https://doi.org/10.1063/1.2969286 ; Design of microfission chamber for ITER operations W. Di Wang , H. R. Cao , J. Cao et al ., A study of beryllium moderator thickness for a fission chamber with fast neutron measurements . Nucl. Sci. Tech . 28 : 1 ( 2017 ). https://doi.org/10.1007/s41365-017-0283-y https://doi.org/10.1007/s41365-017-0283-y A study of beryllium moderator thickness for a fission chamber with fast neutron measurements R. A. Dubridge ., Campbell theorem - system concepts and results . IEEE Trans. Nucl. Sci. NS - 14 : 241 ( 1967 ). https://doi.org/10.1109/TNS.1967.4324422 https://doi.org/10.1109/TNS.1967.4324422 ; Campbell theorem - system concepts and results S.P. Li , X.F. Xu , H.R. Cao , et al ., Dynamic linear calibration method for a wide range neutron flux monitor system in ITER . Nucl. Sci. Tech . 24 : 8 ( 2013 ). https://doi.org/10.13538/j.1001-8042/nst.2013.04.002 https://doi.org/10.13538/j.1001-8042/nst.2013.04.002 ; Dynamic linear calibration method for a wide range neutron flux monitor system in ITER H.R. Cao , S.P. Li , X.F. Xu , et al ., An automatic gain adjustment Campbell integrator for neutron flux detection in ITER . Nucl. Sci. Tech . 23 : 114 ( 2012 ). An automatic gain adjustment Campbell integrator for neutron flux detection in ITER Y. Endo , T. Ito , and E. Seki ., A counting-campbelling neutron measurement system and its experimental results by test reactor . IEEE Trans. Nucl. Sci . 29 : 714 ( 1982 ). https://doi.org/10.1109/TNS.1982.4335943 https://doi.org/10.1109/TNS.1982.4335943 ; A counting-campbelling neutron measurement system and its experimental results by test reactor N. Qiao , D. Li , and H. S. Xiong ., Research on fission chamber signal simulation of wide-range nuclear instrument system . At. Energy Sci. Technol . 47 : 1892 ( 2013 ). (in Chinese) Research on fission chamber signal simulation of wide-range nuclear instrument system J. W. Yang , X. Y. Song , W. Zhang et al ., Development of prototype neutron flux monitor for ITER . Plasma Sci. Technol . 7 : 2860 ( 2005 ). https://doi.org/10.1088/1009-0630/7/3/018 https://doi.org/10.1088/1009-0630/7/3/018 ; Development of prototype neutron flux monitor for ITER 10.1007/s41365-020-00798-3 This work was supported by the National Natural Science Foundation of China (Nos.11475131, 11975307, and 11575184) and the National Magnetic Confinement Fusion Energy Development Research (2013GB104003). 2019-10-24 2020-07-13 2020-07-17 2020-09-03 2020-09-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00798-3 http://dx.doi.org/10.1007/s41365-020-00798-3 . 2020 Nuclear Science and Techniques 9 31