Contrast sensitivity in 14MeV fast neutron radiography

Contrast sensitivity in 14MeV fast neutron radiography 1 2 first-author 1 corresp Corresponding.E-mail address: baojie@ciae.ac.cn Corresponding.E-mail address: baojie@ciae.ac.cn baojie@ciae.ac.cn 3 2 1 1 1 Key Laboratory of Nuclear Data, China Institute of Atomic Energy, Beijing 102413, China Key Laboratory of Nuclear Data, China Institute of Atomic Energy, Beijing 102413, China Xi’an Institute of High Technology, Xi’an 710025, China Xi’an Institute of High Technology, Xi’an 710025, China Science and Technology Information Center, China Academy of Engineering Physics, Mian-yang 621999, China Science and Technology Information Center, China Academy of Engineering Physics, Mian-yang 621999, China Fast neutron radiography (FNR) is an effective non-destructive testing technique. Due to the scattering effect and low detection efficiency, the detection limit of FNR under certain conditions cannot be determined. In order to obtain the minimum detectable thickness by FNR, we studied the contrast sensitivity of FNR lead samples, both theoretically and experimentally. We then clarified the relationship between pixel value and irradiation time, and sample materials and thickness. Our experiment, using a 4-cm-thick lead sample, verified our theoretical expression of FNR contrast sensitivity. Fast neutron radiography Contrast sensitivity Experimental research 1. Introduction 1 Both fast neutron radiography (FNR) and thermal neutron radiography (TNR) are important non-destructive testing techniques. [ 1 1 - 3 3 1-3 ] These 2 methods have similar mechanisms and complementary effects; however, compared to TNR, the neutron in the FNR method has higher penetrability, and can better detect light materials in thick heavy metals, which extends the application of neutron radiography. Presently, there exists a widely recognized international non-destructive testing standard, ASTME545-86, for TNR; however, no quantitative analytical standard has been established for the detection ability of FNR due to its scattering effect and low detection efficiency. [ 4 4 - 6 6 4-6 ] Generally speaking, the FNR technique is not as mature or developed as the TNR technique, which makes it crucial to develop new analysis approaches and criterion for FNR testing. In recent years, international research has mainly focused on the construction and optimization of large-scale imaging systems, [ 7 7 - 9 9 7-9 ] the development of conversion screens, [ 10 10 - 12 12 10-12 ] parameter simulations, [ 13 13 - 15 15 13-15 ] and preliminary applications. [ 16 16 - 18 18 16-18 ] In terms of spatial resolution, FNR uses the KLASENS method to determine the modulation transfer function (MTF) of the system. In terms of contrast sensitivity, FNR uses tiered samples (made from the same materials, but with differing levels of thickness) for preliminary quantitative analysis. [ 19 19 - 21 21 19-21 ] Researchers have not yet solved an inherent problem in FNR's contrast sensitivity, i.e., the longitudinal detection limit of FNR when a fixed-thickness sample made from a certain material is irradiated by neutrons with the same energy and fixed intensity. As such, in this study, we conducted an experiment on the contrast sensitivity of 14 MeV FNR generated by the 600kV Cockcroft-Walton accelerator in the Key Laboratory of Nuclear Data of the China Institute of Atomic Energy. 2. Modeling and Analysis 1 The contrast sensitivity of an FNR system is defined as the minimum discernible variance in the thickness of a fixed-thickness sample along the incidence direction of the detection beam. Contrast sensitivity serves as an index for judging a detection system’s capacity for identifying the longitudinal thickness of a sample. If the thickness variance of a sample is less than the contrast sensitivity of the system, the sample is deemed undetectable. If the variance is equivalent to the contrast sensitivity, it is necessary to consider the comprehensive influence of irradiation time and scattering on the statistical fluctuation of the image contrast. Contrast sensitivity is primarily a TNR technical index, where its expression can be modified according to the characteristics of FNR applications. Figure 1 Figure 1 shows the schematic of FNR data conversion. In Figure 1 Figure 1 , I 0 represents the intensity of incident collimated neutron beams, x is the thickness of the sample, I is the intensity of the neutron beam transmitted through the sample, the yellow rectangle represents a thin scintillator as a convertor, Q is the number intensity of fluorescence photons converted from the neutron beam by the convertor, and P stands for the pixel value of the image on the camera. When passing through the sample, the intensity of the neutron beam I has 2 parts: directly-attenuated neutrons and scattered neutrons, where the latter depends on the geometry and material function of the sample. Therefore, I can be expressed as: <math display="block" id="M1"><mrow><mi>I</mi><mo>=</mo><msub><mi>I</mi><mn>0</mn></msub><msup><mtext>e</mtext><mrow><mo>−</mo><mi>μ</mi><mi>x</mi></mrow></msup><mo>+</mo><msub><mi>I</mi><mn>0</mn></msub><mi>S</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math> where µ is the macroscopic cross-section, µ s is the material cross-section, and S ( x , y , z , µ s ) is the geometry and material function of the sample. The fast neutrons continue to react with the converter. In the experiment, we used a thin BC400 plastic scintillator plate (consisting of carbon and hydrogen) as the convertor. The reactions between fast neutrons and carbon (or hydrogen nuclei) generate recoil nucleons and protons. According to energy and momentum conservation principles, the equation of recoil nucleus energy E A can be expressed as: <math display="block" id="M2"><mrow><msub><mi>E</mi><mi>A</mi></msub><mo>=</mo><msub><mi>K</mi><mi>b</mi></msub><mi>E</mi><mi>c</mi><mi>o</mi><msup><mi>s</mi><mn>2</mn></msup><mi>θ</mi><mo>,</mo></mrow></math> where E represents the incident neutron energy, θ isthe angle between the recoil nucleon and the incident neutron, A is the mass number of the recoil nucleon, and <math id="M3"><mrow><msub><mi>K</mi><mtext>b</mtext></msub><mo>=</mo><mfrac><mrow><mn>4</mn><mi>A</mi></mrow><mrow><mi>（</mi><mn>1</mn><mo>+</mo><mi>A</mi><msup><mi>）</mi><mn>2</mn></msup></mrow></mfrac></mrow></math> . E A represents the energy of a recoil carbon nucleon if A=12, and E A is a recoil proton if A=1. The ratio of K b between a recoil carbon nucleon and a recoil proton is 1:3.57. Therefore, the energy of a recoil proton is much larger than that of a recoil carbon nucleon under the same initial conditions. Even for the same energy, a recoil proton generates more photon yield than a recoil carbon nucleon because of different reaction cross-sections and ionization energy. Accordingly, the reaction of fast neutrons with a scintillator BC400 can be simplified as an interaction between neutrons and protons. During the reaction process, the recoil protons impact the BC400 scintillator and generate fluorescence photons, whose intensity Q is: <math display="block" id="M4"><mrow><mi>Q</mi><mo>∝</mo><mi>ε</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mi>I</mi><msub><mi>σ</mi><mrow><mi>B</mi><mi>C</mi></mrow></msub><mi>m</mi><mi>n</mi><mo> </mo><mo>,</mo></mrow></math> Where ε( E ) is the energy deposition function of the recoil proton energy ( E ) in the BC400 scintillator, σ BC is the reaction cross-section of the fast neutrons with hydrogen in the BC400 scintillator, m is the inherent nucleon density of the BC400 scintillator, n is the coefficient used to compensate for the recoil effect of carbon nucleons on photon intensity Q . During the conversion process, we regarded each spot on the convertor that generated photons as a point of light source that illuminated at a full 4π solid angle. Blocked by the converter and influenced by the transmission path, the camera only recorded photons entering the camera lens at a very small solid angle. Ultimately, the fluorescence photons passed through the air, were reflected by the mirror, and transmitted through the lens before being recorded by the CCD, which inevitably led to some attenuation. Considering all of these factors, we chose a parameter L to represent the equivalent distance from the light source to the camera, and coefficients K 1 , K 2 ,and K 3 to take into account airborne transmission efficiency, the reflection efficiency of the mirror, and the transmittance capacity of the lens, respectively. The unit-time pixel value P / T detected by the CDD camera can be expressed as: <math display="block" id="M5"><mrow><mfrac><mi>P</mi><mi>T</mi></mfrac><mo>=</mo><mfrac><mrow><mi>ϑ</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow><mi>h</mi><mi>v</mi><mi>Q</mi></mrow><mrow><mn>4</mn><mi>π</mi><msup><mi>L</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>K</mi><mn>1</mn></msub><msub><mi>K</mi><mn>2</mn></msub><msub><mi>K</mi><mn>3</mn></msub><mo>,</mo></mrow></math> where <math id="M6"><mrow><mi>ϑ</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></mrow></math> represents a positive correlation function between the pixel value P and the photon energy e , h stands for the Planck constant, v stands for the photon frequency, and T is the irradiation time. If the imaging system is free of stray light and the camera works stably under electrical refrigeration, the background image B(w,v) with no neutron beam depends solely on the pixel coordinates w and v , and is constant with respect to time. By combining Equations (2) and (4), the pixel value P is: <math display="block" id="M7"><mrow><mi>P</mi><msub><mrow><mtext>=λTI</mtext></mrow><mn>0</mn></msub><msup><mtext>e</mtext><mrow><mo>−</mo><mtext>μx</mtext></mrow></msup><mo>+</mo><msub><mrow><mtext>λTI</mtext></mrow><mn>0</mn></msub><mtext>S</mtext><mrow><mo>(</mo><mrow><mtext>x</mtext><mo>,</mo><mtext>y</mtext><mo>,</mo><mtext>z</mtext><mo>,</mo><msub><mtext>μ</mtext><mtext>s</mtext></msub></mrow><mo>)</mo></mrow><mo>+</mo><mtext>B</mtext><mrow><mo>(</mo><mrow><mtext>w</mtext><mo>,</mo><mtext>v</mtext></mrow><mo>)</mo></mrow><mtext>T,</mtext></mrow></math> where <math id="M8"><mrow><mi>λ</mi><mo>=</mo><mfrac><mrow><mi>ε</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mi>ϑ</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow><mi>h</mi><mi>v</mi><msub><mi>σ</mi><mrow><mi>B</mi><mi>C</mi></mrow></msub><mi>m</mi><mi>n</mi></mrow><mrow><mn>4</mn><mi>π</mi><msup><mi>L</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>K</mi><mn>1</mn></msub><msub><mi>K</mi><mn>2</mn></msub><msub><mi>K</mi><mn>3</mn></msub></mrow></math> is a parameter related to the entire testing and imaging system. In our experiment, we only considered variances in thickness Δ x along the direction of the neutron beam. By differentiating both sides of Equation (5), we obtained <math display="block" id="M9"><mrow><mtext>Δ</mtext><mi>P</mi><msub><mrow><mtext>=-μλTI</mtext></mrow><mtext>0</mtext></msub><msup><mtext>e</mtext><mrow><mtext>-μx</mtext></mrow></msup><msub><mrow><mtext>Δx+λTI</mtext></mrow><mtext>0</mtext></msub><mover accent="true"><mi>S</mi><mo>˙</mo></mover><msub><mrow><mtext>(x,y,z,μ</mtext></mrow><mtext>s</mtext></msub><mtext>)Δx</mtext><mtext>.</mtext></mrow></math> By setting Δ P as the minimum signal in the pixel value that can be detected by the computer system after digital statistics (which is related to the fluctuation of the system’s counting statistics) we were able to calculate the minimum variance in thickness Δ x as: <math display="block" id="M10"><mrow><mi>Δ</mi><mi>x</mi><mo>=</mo><mi>Δ</mi><mi>P</mi><mo>/</mo><mi>λ</mi><mi>T</mi><msub><mi>I</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mover accent="true"><mi>S</mi><mo>˙</mo></mover><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>-</mo><mi>μ</mi><msup><mi>e</mi><mrow><mo>−</mo><mi>μ</mi><mi>x</mi></mrow></msup><mo stretchy="false">)</mo><mo>.</mo></mrow></math> If the scattering term <math id="M11"><mrow><mover accent="true"><mi>S</mi><mo>˙</mo></mover><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow></math> had a minor effect on the pixel value, then we were able to simplify Equation (7) as: <math display="block" id="M12"><mrow><mo>|</mo><mi>Δ</mi><mi>x</mi><mo>|</mo><mo>≈</mo><mo>|</mo><mo stretchy="false">(</mo><mi>Δ</mi><mi>P</mi><msup><mi>e</mi><mrow><mi>μ</mi><mi>x</mi></mrow></msup><mo>/</mo><mo stretchy="false">(</mo><mi>λ</mi><mi>T</mi><msub><mi>I</mi><mn>0</mn></msub><mi>μ</mi><mo stretchy="false">)</mo><mo>|</mo><mo>.</mo></mrow></math> As shown in Equation (8), when the influence of scattering could be ignored (i.e., if the scattering was too insignificant to be considered), the contrast sensitivity of the system, (the minimum discernible thickness (Δ x )) was inversely proportional to the irradiation time T , the macroscopic cross-section μ of the sample, and the intensity of the neutron source I 0 . Moreover, it had a positive correlation to the minimum discernible signal of the pixel value Δ P . 3. Experimental Study 1 Figure 2(a) Figure 2(a) depicts the schematic diagram of the relative positions of the sample and experimental facility. We generated a 14.1 MeV fast neutron source with a flux of ~1.5×10 10 n/s and an energy of 14.1 MeV using a 600kV high voltage multiplier. The neutron source emitted into the collimator from the large inlet at a 4πsolid angle, and incidents perpendicularly on the sample in a plane wave after passing through the collimator. We used a BC400 plastic scintillator plate with a size of 200 mm×200 mm×10 mm, a density of 1.032 g/cm3, and an H:C~1.103:1as the convertor. The collimator we used in the experiment was a composite structure consisting of lead, polyethylene, stainless steel, and red copper. It had a length of L =147 cm, with a maximum diameter at the outlet of D =8 cm, and a minimum diameter inside of d =5.14 cm. Based on the results of the experiment and the MC simulation, the flux ratios between the neutrons and the γ-ray near the collimator inlet (on the neutron source side) and the outlet were ~1:2.81*10 -3 and ~1:5.68*10 -4 , respectively. The ratio between the neutrons and the γ flux at the collimator outlet was approximately 25:1. Based on the experimental and simulation data, we determined that the X-ray and the γ-ray contribution to the dose rate on the BC400 scintillator was relatively low when compared to that of the neutrons. Futhermore, based on the simulation, the average conversion efficiency of the X-ray and the γ-ray with a 1 cm thickness on the BC400 scintillator was only 0.2%; whereas the neutron efficiency was approximately 2.2%, which was much higher than the gamma efficiency. The peak center of the BC400 scintillation light spectrum was located in the UV range (maximum wavelength 423 nm), which could barely be recorded by the CCD camera. In short, the experimental images’ contribution to the pixel values from the X-ray and the γ-ray were relatively small compared to that of the neutrons. The camera was a 1024×1024-pixel high-sensitivity scientific CCD. Each pixel corresponded to a 0.084 mm×0.084 mm area on the converter. The sample was located 2 cm from the outlet of the collimator. The system collimation ratio was 294, and the geometric dullness of the system U g was ~0.0029 cm. We processed all of the images in the experiment after dark field deduction and flat field analysis had occurred. The testing times of the dark field and flat field were equivalent to the exposure time of the experiment images. Based on the relationship between various parameters in the sample system and the pixel value of the image as stated in Equation (5), we used 2 types of samples consisting of different materials and varied thickness. We then used different irradiation times for the 2 types of samples in the experiments. We chose lead and polyethylene as the materials for the 2 samples, because the former (lead) has a weak scattering effect of fast neutrons, while the scattering effect of the latter (polyethylene) is strong. Both samples were arranged in the shape of a triangular prism with a 5 cm×10 cm×5√5 cm right triangle as the base, and a thickness of 5 cm, as shown in Figures 2(b) and 2(c) Figures 2(b) and 2(c) , which also illustrate the direction of the transmitted neutrons. Our analysis consisted of 3 steps: (1) First, we compared the influence of the different sample materials on the image pixel values; (2) Second, we derived the relationship between the irradiation time and the image pixel values, and (3) Third, we determined how the variations of sample thickness affected the pixel values for a certain material. Once we determined all 3 relationships, we were able to establish the corresponding discriminant of the system contrast sensitivity. 3.1 Experimental Study on Influential Factors in Pixel Values 2 Figure 3(a)-(c) Figure 3(a)-(c) shows the fluorescence images of the lead sample’s irradiation times at 1200 s and 1800 s, and the irradiation time of the polyethylene sample at 1800 s, respectively. The change of thickness in the images is presented in a horizontal direction. We subtracted the background of all of the fluorescence images to ensure equal irradiation times. We also filtered the medians before saving the data in ASCII format. Figure 3(d)-(f) Figure 3(d)-(f) shows the experimental data and fitting curves of the pixel values versus the sample thickness for the corresponding configurations. By visually comparing the images in Figure 3(a)-(c) Figure 3(a)-(c) , one can see that Figure 3(c) Figure 3(c) has the lightest coloring and superior gray scale layering, while Figure 3(a) Figure 3(a) has the darkest coloring and poorest gray scale layering. When analyzing the pixel values of the images, we chose a specific area marked as a rectangle with red edges, as shown in Figure 3(a)-(c) Figure 3(a)-(c) . By averaging the pixel values along the width of the rectangle, we obtained the experimental data of the pixel values with respect to changes in thickness (along the length of the rectangular area) as shown in Figure 3(d)-(f) Figure 3(d)-(f) . In order to choose an appropriate function to fit the experiment data of the pixel value referred to in Equation (5), we eliminated the term B(w,v) T because we had already subtracted the background images from the data. We then ignored the scattering term S ( x , y , z , µ s ) because it did not have an explicit expression. Next, we fitted the function using Y = KAe -µx , and evaluated the reliability of the fitted curve based on the goodness-of-fit R 2 . In Y = KAe -µx , Y , K , A , µ , x , and x are related to the pixel value P ,system parameter λ, neutron accumulation TI 0 , material of the sample, and the thickness of the sample, respectively. Table 1 Table 1 shows the fitting functions of the experimental data in the rectangular areas in Figure 3(a)-(c) Figure 3(a)-(c) and the goodness-of-fit of each curve. 1 1 Name 1 1 Sample 1 1 Fitting Function 1 1 Goodness-of-Fit 1 1 a 1 1 1200 s Lead Sample 1 1 Y = 840.3130e - 0.0131 x 1 1 R² = 0.9925 1 1 b 1 1 1800 s Lead Sample 1 1 Y= 1211.3447e - 0.0131 x 1 1 R² = 0.9943 1 1 c 1 1 1800 s Polyethylene Sample 1 1 Y= 1330.3464e - 0.0087 x 1 1 R² = 0.9902 Judging by the goodness-of-fit data, the adoption of function Y = KAe -µx provided consistently high degrees of fitting, which indicates that choosing the exponential function Y as the fitting function was appropriate. However, compared to the lead sample of 1200 s of irradiation time ( Figure 3(a) Figure 3(a) ) and 1800 s of irradiation time ( Figure 3(b) Figure 3(b) ), the goodness-of-fit of the 1800-s-irradiation-time polyethylene sample ( Figure 3(c) Figure 3(c) ) was the lowest. This may be the result of polyethylene having a stronger scattering effect on fast neutrons (compared to lead samples), where the scattering term S ( x , y , z , µ s ) may have had a relatively larger contribution to the pixel values of the images, and should be treated carefully rather than being ignored. The lead samples with 1200-s-irradiation-times and 1800-s-irradiation-times had the same fitting value of µ (0.0131), while the 1800-s-irradiation-time polyethylene sample had a fitting value of µ (0.0087). This demonstrates that the parameter µ depends on the materials of the sample, and is free of the constraints of time. When considering the value of parameter KA , we depended on the system parameter λ , the irradiation time T , and neutron intensity I 0. We determined that if all of the tests were achieved in one system and had stable neutron sources and equivalent irradiation times, the KA values were the same for all samples, regardless of material type. However, the fitted data of the 1800-s-irradiation-time lead sample and the 1800-s-irradiation-time polyethylene sample had KA values of 1211.34 and 1330.34, respectively. The disparity in KA values between the lead and polyethylene samples with equal irradiation times implies that the scattering term S ( x , y , z , µ s ) exerts influence on the pixel values, or that the drift of the neutron beam intensity I 0 for a relatively long time (thousands of seconds) has an impact on KA values. Overall, for Figure 3(a)-(c) Figure 3(a)-(c) , the pixel value undergoes a negative exponential function with respect to approximate thickness x . Our analysis showed that Y = KAe -µx can satisfy the fitting of pixel value curves under certain limitations, such as the scattering effect. Moreover, the results of our experiment agree well with the theoretical derivation of the pixel values. 3.2 Establishing the Contrast Sensitivity Discriminants 2 We used the quantitative relationships between pixel values and irradiation time, and sample thickness and sample material to establish a contrast sensitivity discrimination standard. After obtaining the function Y = KAe -µx of the pixel values, followed by changes in thickness, the pixel value Y 0 at thickness X 0 can be calculated as: <math display="block" id="M13"><mrow><msub><mi>Y</mi><mn>0</mn></msub><mo>=</mo><mi>K</mi><mi>A</mi><msup><mi>e</mi><mrow><mo>−</mo><mo>µ</mo><msub><mi>x</mi><mrow><msub><mrow/><mn>0</mn></msub></mrow></msub></mrow></msup></mrow></math> According to the fluctuation theory of counting statistics, when the count is Y 0 , the range of the statistical fluctuation is: <math display="block" id="M14"><mrow><msub><mi>Y</mi><mn>0</mn></msub><mo>−</mo><msqrt><mrow><msub><mi>Y</mi><mn>0</mn></msub></mrow></msqrt><mo><</mo><msub><mi>σ</mi><mn>0</mn></msub><mo><</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>+</mo><msqrt><mrow><msub><mi>Y</mi><mn>0</mn></msub></mrow></msqrt></mrow></math> When the pixel value Y falls within the range of σ 0 , it is deemed undetectable under condition Y 0 ; when Y is out of the range of σ 0 , it is assumed to be detectable. Based on Equations (9) and (10), when the sample thickness was X 0 and the irradiation intensity was A , we calculated the minimum detectable concave depth of the sample X c and the minimum detectable convex height X v as follows: <math display="block" id="M15"><mrow><msub><mi>X</mi><mtext>c</mtext></msub><mo>=</mo><msub><mi>X</mi><mn>0</mn></msub><mo>−</mo><mfrac><mn>1</mn><mi>μ</mi></mfrac><mi>ln</mi><mfrac><mrow><mi>K</mi><mi>A</mi></mrow><mrow><msub><mi>Y</mi><mn>0</mn></msub><mo>+</mo><msqrt><mrow><msub><mi>Y</mi><mn>0</mn></msub></mrow></msqrt></mrow></mfrac><mo>,</mo></mrow></math> <math display="block" id="M16"><mrow><msub><mi>X</mi><mtext>v</mtext></msub><mo>=</mo><mfrac><mn>1</mn><mi>μ</mi></mfrac><mi>ln</mi><mfrac><mrow><mi>K</mi><mi>A</mi></mrow><mrow><msub><mi>Y</mi><mn>0</mn></msub><mo>−</mo><msqrt><mrow><msub><mi>Y</mi><mn>0</mn></msub></mrow></msqrt></mrow></mfrac><mo>−</mo><msub><mi>X</mi><mn>0</mn></msub></mrow></math> Equations (11) and (12) are the discriminants of FNR contrast sensitivity, which can be used to determine the detection limit of FNR under the irradiation of fast neutrons with specific energy, fixed intensity, and for a sample of specific materials and thickness. The capacity of an FNR system to identify material defects depends on its spatial resolution and ability to identify variations in the thickness in the longitudinal direction. Take groove defects for an example: when the width of a groove is much larger than the spatial resolution of the imaging system, the system’s ability to identify the groove depends solely on whether the depth of the groove is larger than the minimum discernible thickness, i.e., the contrast sensitivity of the system. Otherwise, the system’s ability to identify defects is affected by the system’s spatial resolution on the edge extension. Spatial resolution and scattering terms are not considered in Equations (11) and (12);therefore, Equations (11) and (12) are applicable only in cases where (1) the dimension of a defect is much larger than the system’s spatial resolution, (2) the edge that broadens due to the spatial resolution limit has less of an effect on the pixel values than the effect of sample thickness variation (i.e., the effects of sample thickness variation must exert greater influence on the pixel values), and (3) the scattering effect is weak. 3.3 Experimental Verification of the Reliability of Contrast Sensitivity Discriminants 2 In order to verify the reliability of the contrast sensitivity discriminants, we conducted an analysis based on the fitted curve Y =1211.3447 e -0.0131 x with high goodness-of-fit for the lead sample with an 1800 s irradiation time. When the sample’s thickness X 0 was 40 mm, the corresponding pixel value Y 0 was 717.3 under the conditions of the experiment. Based on Equation (10), pixel values fall into the range of 717.3±26.8, and were considered undetectable. The minimum detectable concave depth X c and convex height X v were 2.8 mm and 2.9 mm, respectively, using Equations (11) and (12). Considering systematic errors, we aimed at the minimum discernible thickness of 3 mm for a 40mm thick lead sample in the FNR system. For the experiment, we designed 2 cuboid lead samples with geometrical proportions of 60 mm×40 mm×20 mm. One of the samples had 5 grooves with different depths on the 60 mm×40 mm surface (as shown in Figure 4(a) Figure 4(a) ). The depth and width for each slot were the same: 1 mm, 2 mm, 3 mm, 4 mm, and 5 mm for the 5 grooves, respectively. During the experiment, the 2 samples were stacked to form a block measuring 60 mm×40 mm×40 mm. The experimental conditions included the neutron beam intensity and irradiation times being consistent with those in Equations (11) and (12). Figure 4(b) Figure 4(b) shows the fluorescence image of the lead sample in which the grooves (easily distinguishable by eye) have depths of 3 mm, 4 mm, and 5 mm, respectively. We vertically added the pixel values in the red rectangle in Figure 4(b) Figure 4(b) using the software Andor Solis for Imaging . We then plotted a statistical diagram of the pixel values in the horizontal axis, from which one can see 3 peaks which correspond to the groove depths of 3 mm, 4 mm, and 5 mm, respectively. Since the pixel values associated with the 1 mm and 2 mm deep grooves were equivalent to the statistical fluctuations, they were undetectable. We only analyzed data from the central area of the images because they were more creditable than the areas near the boundary. In conclusion, the minimum detectable thickness of the 40mm thick lead sample was ~3 mm, which demonstrated the reliability of Equations (11) and (12). Moreover, the contrast sensitivity discriminant theory stated above was applicable for evaluating the minimum detectable thickness of a fix-thickness low-scattering sample in an FNR testing system. When applying such a method, attention must be paid to the applicable range and conditions of the criterion. 4. Conclusion 1 We conducted a study on the contrast sensitivity of a 14 MeV FNR system using the 600kV Cockcroft-Walton accelerator in the China Institute of Atomic Energy. We derived an equation for pixel values with respect to system sample parameters and irradiation time, and determined the theoretical expression of contrast sensitivity. We selected lead and polyethylene samples to derive the relationship between pixel values and irradiation times, and sample materials and thickness.; We then established the discriminants of contrast sensitivity and the reliability of the fitted curves. Afterwards, we verified the discriminants using a 40 mm thick lead sample. References M. Segawa , T. Kai , T. Sakai et al ., Development of a high-speed camera system for neutron imaging at a pulsed neutron source . Nucl. Instrum. Methods A 697 , 77 - 83 ( 2013 ). doi: 10.1016/j.nima.2012.08.079 http://doi.org/10.1016/j.nima.2012.08.079 Development of a high-speed camera system for neutron imaging at a pulsed neutron source B. Tang , C. G. Zhou , H. Y. Huo et al ., The pilot experimental study of 14 MeV fast neutron digital radiography . Sci. China Ser. G . 52 , 1090 - 1096 ( 2009 ). (in Chinese) doi: 10.1007/s11433-009-0184-0 http://doi.org/10.1007/s11433-009-0184-0 The pilot experimental study of 14 MeV fast neutron digital radiography S. Fujinei , K. Yoneda , K. Yoshii et al ., Development of imaging techniques for fast neutron radiography in Japan . Nucl. Instrum. Methods A 424 , 190 - 199 ( 1999 ). doi: 10.1016/S0168-9002(98)01326-6 http://doi.org/10.1016/S0168-9002(98)01326-6 Development of imaging techniques for fast neutron radiography in Japan H. Li , Y. Wu , C. Cao et al ., Design optimization, manufacture and response measurements for fast-neutron radiography converters made of scintillator and wavelength-shifting fibers . Nucl. Instrum. Methods A 762 , 64 - 69 ( 2014 ). doi: 10.1016/j.nima.2014.05.079 http://doi.org/10.1016/j.nima.2014.05.079 Design optimization, manufacture and response measurements for fast-neutron radiography converters made of scintillator and wavelength-shifting fibers Y. B. Zou , L. A. Guo , Z. Y. Guo et al ., Development of a converter made of scintillator and wavelength-shifting fibers for fast neutron radiography . Nucl. Instrum. Methods A 651 , 62 - 66 ( 2011 ). doi: 10.1016/j.nima.2009.01.130 http://doi.org/10.1016/j.nima.2009.01.130 Development of a converter made of scintillator and wavelength-shifting fibers for fast neutron radiography R. Zboray , I. Mor , V. Dangendorf et al ., High-frame rate imaging of two-phase flow in a thin rectangular channel using fast neutrons . Appl. Radiat. Isot . 90 , 122 - 131 ( 2014 ). doi: 10.1016/j.apradiso.2014.03.023 http://doi.org/10.1016/j.apradiso.2014.03.023 High-frame rate imaging of two-phase flow in a thin rectangular channel using fast neutrons B. Tang , Y. Wu , H. Li et al ., The physics analysis and expriment study of zinc sulphide scintillator for fast neutron radiography . Nucl. Instrum. Methods A 729 , 327 - 333 ( 2013 ). doi: 10.1016/j.nima.2013.07.063 http://doi.org/10.1016/j.nima.2013.07.063 The physics analysis and expriment study of zinc sulphide scintillator for fast neutron radiography H. Wang , V. Tang , M. J. Carrick et al ., Reconstruction algorithm for point source neutron imaging through finite thickness scintillator . Nucl. Instrum. Methods A 693 , 294 - 301 ( 2012 ). doi: 10.1016/j.nima.2012.07.018 http://doi.org/10.1016/j.nima.2012.07.018 Reconstruction algorithm for point source neutron imaging through finite thickness scintillator S. Bishnoi , R.G. Thomas , P. S. Sarkar et al ., Simulation study of fast neutron radiography using GEANT4 . JINST . 10 , P02002 ( 2015 ). doi: 10.1088/1748-0221/10/02/P02002 http://doi.org/10.1088/1748-0221/10/02/P02002 Simulation study of fast neutron radiography using GEANT4 P. Andersson , J. B. Valldor , E. S. Andersson et al ., Design and initial 1 D radiography tests of the FANTOM mobile fast neutron radiography and tomography system . Nucl. Instrum. Methods A 756 , 82 - 93 ( 2014 ). doi: 10.1016/j.nima.2014.04.052 http://doi.org/10.1016/j.nima.2014.04.052 Design and initial 1 D radiography tests of the FANTOM mobile fast neutron radiography and tomography system J. M. Qi , F. Q. Zhang , J. C. Chen et al ., Experimental research on light output characteristic of scintillator applied in pulsed y-ray and fast neutron imaging systems . Atom. Energ. Sci. Technol . 48 , 715 - 720 ( 2014 ). doi: 10.7538/yzk.2014.48.04.0715 http://doi.org/10.7538/yzk.2014.48.04.0715 Experimental research on light output characteristic of scintillator applied in pulsed y-ray and fast neutron imaging systems C. B. Lu , P. Xu , J Bao et al ., Simulation analysis and experimental verification of fast neutron radiography . Acta. Phys. Sin . 64 , 198702 ( 2014 ). (in Chinese) . doi: 10.7498/aps.64.198702 http://doi.org/10.7498/aps.64.198702 Simulation analysis and experimental verification of fast neutron radiography G. Nicholas , Y. Liu , R. Y. James et al ., Development and commercialization of a fast-neutron/X-ray cargo scanner . Proc. Scie. Eng . 978 , 330 - 336 ( 2010 ). Development and commercialization of a fast-neutron/X-ray cargo scanner H. Li , Y. Zou , S. Wang et al ., Preliminary study on neutron radiography with several hundred keV fast neutrons . Chin. Phys. Lett . 29 , 1 - 4 ( 2012 ). doi: 10.1088/0256-307X/29/11/118701 http://doi.org/10.1088/0256-307X/29/11/118701 Preliminary study on neutron radiography with several hundred keV fast neutrons C. B. Lu , P. Xu , J. Bao et al ., Preliminary experimental research of detail sensitivity in fast neutron radiography . Nucl. Tech . 38 , 080201 ( 2015 ). (in Chinese) . doi: 10.11889/j.0253-3219.2015.hjs.38.080201 http://doi.org/10.11889/j.0253-3219.2015.hjs.38.080201 Preliminary experimental research of detail sensitivity in fast neutron radiography L. F. He , S. B. Han , H. L. Wang et al ., Design of real-time neutron radiography at China advanced research reactor . Phys. Proc . 43 , 48 - 53 ( 2013 ). Design of real-time neutron radiography at China advanced research reactor I. S. Anderson , A. J. Hurd , R. L. McGreevy et al ., Neutron imaging and applications . Springer-Verlag New York Inc , 2009 . Neutron imaging and applications Z. Q. Cui , Z. J. Chen , X. F. Xie , et al ., Design of a magnetic shielding system for the time of flight enhanced diagnostics neutron spectrometer at Experimental Advanced Superconducting Tokamak . Rev. Sci. Instrum . 85 , 11D829 ( 2014 ). doi: 10.1063/1.4890536 http://doi.org/10.1063/1.4890536 Design of a magnetic shielding system for the time of flight enhanced diagnostics neutron spectrometer at Experimental Advanced Superconducting Tokamak L. F. He , S. B. Han , Y. Wang et al ., Calculation and analysis of the neutron radiography spatial resolution . Nucl. Tech . 37 , 040502 ( 2014 ). (in Chinese) Calculation and analysis of the neutron radiography spatial resolution X. H. Zheng , Y. B. Nie , X. C. Ruan et al ., Experimental study of fast neutron radiography based on Cockcroft Walton Generator . Atom. Energ. Sci. Technol . 49 , 1117 - 1121 ( 2015 ). doi: 10.7538/yzk.2015.49.06.1117 http://doi.org/10.7538/yzk.2015.49.06.1117 Experimental study of fast neutron radiography based on Cockcroft Walton Generator C. Cao , S. Wang , K. Tang et al ., Comparison of two approaches to magnetic field quantification by polarized neutron images . Acta. Phys. Sin . 63 , 182801 ( 2014 ). (in Chinese) . doi: 10.7498/aps.63.182801 http://doi.org/10.7498/aps.63.182801 Comparison of two approaches to magnetic field quantification by polarized neutron images We would like to express our heartfelt gratitude to Professor Zhou Zu-Ying and Professor Tang Hong-Qing (of the China Institute of Atomic Energy) for their guidance. We would also like to thank the accelerator team for providing high-quality neutron beams. 10.1007/s41365-017-0228-5 2016-04-19 2016-12-01 2016-12-05 2017-04-28 2017-06-01 2021-04-11 © 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. 2017 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-017-0228-5 http://dx.doi.org/10.1007/s41365-017-0228-5 . 2017 Nuclear Science and Techniques 06 28