A method for quickly obtaining the steady-state space electric field of the insulating layer of a self-powered detector

Abstract

The application discloses a method for quickly obtaining a steady-state space electric field of an insulation layer of a self-powered detector, and comprises the following steps: obtaining radiation field information and detector information at a position where a self-powered neutron detector is located, and simulating to obtain an electric charge deposition amount generated by irradiation current in different radius intervals per unit time of the detector; according to the law of conservation of electric charge, estimating an electric charge loss amount generated by conduction current in different radius intervals of an insulator per unit time under a steady-state condition; combining Ohm's law to establish an expression of the electric charge loss amount generated by the conduction current about an electric field intensity under the steady-state condition, and combining a boundary condition to obtain an electric field intensity distribution, an electric charge density distribution and a potential distribution. The application can effectively improve the calculation efficiency and realize efficient simulation of the sensitivity of the self-powered neutron detector.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of neutron detection, and particularly relates to a method for quickly obtaining a steady-state space electric field of an insulation layer of a self-powered neutron detector. BACKGROUND

[0002] The self-powered neutron detector has been widely applied to the measurement of neutrons in a reactor due to its simple structure, high temperature and pressure resistance, and strong radiation resistance, and is mainly composed of an emitter, an insulation layer and a collector. The emitter is generally a metal material with a large neutron capture cross section, the middle part is provided with an insulation material, and the collector is provided with an alloy with good conductivity. When the detector is in the reactor, the emitter material can capture neutrons, directly generate electrons through decay or generate secondary electrons through gamma ray emission. The electrons can pass through the insulation layer to reach the collector to generate an electric current, and the electric current size is related to the neutron flow intensity, so that the measurement of the neutrons in the reactor can be realized.

[0003] In the working process of the self-powered neutron detector, electrons may stay in the insulation layer due to insufficient kinetic energy, so as to form a space electric field in the insulation layer. The formed space electric field can affect the transport behavior of subsequent electrons in the insulation layer, and further affect the sensitivity of the detector. Therefore, in order to improve the calculation accuracy of the sensitivity of the detector, the space charge effect (also referred to as electrostatic effect) of the insulation layer needs to be quantitatively analyzed. In addition, when designing the detector, the maximum electric field strength of the insulation layer of the detector must be calculated in combination with the working environment, so as to ensure that the maximum electric field strength does not exceed the breakdown strength of the insulation material, otherwise the insulation layer will fail due to breakdown, and finally the detector will fail.

[0004] At present, there are many studies on the electrostatic effect of the self-powered neutron detector, but these studies still have some deficiencies: 1) in most studies, the charge distribution form of the insulation layer is assumed, which is quite different from the actual situation; 2) some studies simulate the charge deposition by the Monte Carlo method, and add a real electric field in the simulation process, and the cumulative process of the electric field is simulated by iteration, which can reflect the distribution and dynamic change of the electric field, but the calculation efficiency of the simulation of the charge deposition is low, especially the simulation of the interaction between the ray and the material is very complex under the condition of containing the electric field, and the calculation time is long, which is difficult to optimize the design; 3) under the actual commercial reactor working condition, the influence of the electric field on the sensitivity of the detector is small, but when designing the size of the detector and selecting the detector material, the size of the space electric field of the insulation layer of the detector needs to be calculated to evaluate whether breakdown occurs, at this time, the cumulative process of the electric field is simulated by iteration, which is too time-consuming, and it is difficult to meet the design requirements. SUMMARY

[0005] The application aims at solving the problems in the prior art, and provides a method for quickly obtaining a steady-state space electric field of an insulation layer of a self-powered neutron detector, so as to solve the problems of too large calculation amount and too long time consumption in obtaining the space electric field of the insulation layer of the self-powered neutron detector.

[0006] To achieve the above object, the application adopts the following technical scheme:

[0007] The method for quickly obtaining the steady-state space electric field of the insulation layer of the self-powered neutron detector comprises the following steps of:

[0008] obtaining radiation field information and detector information at a position where the self-powered neutron detector is located, and simulating to obtain a charge deposition amount generated by irradiation current in different radius intervals in a unit time of the detector;

[0009] based on the charge deposition amount generated by the irradiation current in the different radius intervals in the unit time of the detector, estimating a charge loss amount generated by conduction current in different radius intervals of the insulator in a unit time under a steady-state condition according to a charge conservation law;

[0010] combining Ohm's law, establishing an expression of the charge loss amount generated by the conduction current about an electric field intensity under the steady-state condition, and combining a boundary condition to obtain a distribution of the electric field intensity, a charge density distribution and a potential distribution.

[0011] Further, the charge deposition amount generated by the irradiation current in the different radius intervals in the unit time of the detector is simulated by Monte Carlo simulation software.

[0012] The radiation field information and the detector information at the position where the self-powered neutron detector is located are input into a Monte Carlo program.

[0013] The interaction between rays and the detector is simulated by the Monte Carlo program, and the charge deposition amount generated by the irradiation current in different radius intervals of the detector is counted.

[0014] Further, the radiation field information comprises particle types, particle energy and momentum distribution and radiation field flux density, and the detector information comprises a detector geometric structure and a material.

[0015] Further, the calculation process of the charge loss amount generated by the conduction current in different radius intervals of the insulator in a unit time under the steady-state condition is as follows:

[0016] Under the steady-state condition, the free charge density in the insulator is unchanged, that is,

[0017]

[0018] It is indicated that the charge deposition amount generated by the irradiation current and the charge loss amount generated by the conduction current in a unit time are the same in value and opposite in sign, that is,

[0019] Q esp (r~r+Δr)=-Q dep (r~r+Δr)

[0020] Substitute the conduction current density, we get

[0021] Q esp (r~r+Δr)=-J cond (r)S(r)+J cond (r+Δr)S(r+Δr)

[0022] =-2πL[J cond (r)*r-J cond (r+Δr)*(r+Δr)]

[0023] Where, Q dep (r~r+Δr) represents the amount of charge deposited per unit time generated by the irradiation current, Q esp (r~r+Δr) represents the amount of charge loss per unit time generated by the conduction current, L represents the length of the detector, J cond (r) represents the conduction current density at radius r in steady state.

[0024] Further, in the steady state, the amount of charge loss generated by the conduction current with respect to the electric field intensity is expressed as:

[0025]

[0026] Where, σ represents the conductivity of the detector insulating layer, L represents the length of the detector, r k represents the inner radius of the kth cylindrical shell layer when the insulating layer of the detector is evenly divided into N cylindrical shell layers along the radius, r k+1 represents the outer radius, k = 1, 2, 3, …, N; Q k represents the amount of charge loss per unit time generated by the conduction current in the interval of radius (r k , r k+1 ), E k represents the electric field intensity at radius r k .

[0027] Further, in the steady state, the calculation process of the electric field intensity distribution, charge density distribution and potential distribution is:

[0028] Combined with the boundary conditions, a linear equation set of electric field intensity and the amount of charge loss per unit time generated by the conduction current is established, and is arranged in matrix form;

[0029] According to the relationship between potential and electric field, the potential difference ΔV of the inner and outer conductors is represented as

[0030]

[0031] The electric field is integrated to calculate the electric field intensity distribution, the charge density distribution and the potential distribution.

[0032] Further, the method for integrating the electric field is a numerical integration method or a mirror method.

[0033] A system for quickly obtaining a steady-state space electric field of an insulation layer of a self-powered detector, comprising:

[0034] An acquisition module is configured to acquire radiation field information and detector information at a position where the self-powered neutron detector is located, and simulate charge deposition amounts generated by irradiation current in different radius intervals per unit time of the detector;

[0035] An estimation module is configured to estimate, according to the charge conservation law, charge loss amounts generated by conduction current in different radius intervals of the insulator per unit time in a steady state based on the charge deposition amounts generated by the irradiation current in the different radius intervals per unit time of the detector;

[0036] A calculation module is configured to establish an expression of the charge loss amounts generated by the conduction current with respect to the electric field intensity in a steady state in combination with Ohm's law, and acquire the electric field intensity distribution, the charge density distribution and the potential distribution in combination with a boundary condition.

[0037] A terminal device comprises a memory, a processor and a computer program stored on the memory and executable on the processor, and the processor implements the steps of the method when executing the computer program.

[0038] A computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the steps of the method.

[0039] Compared with the prior art, the present application has the following beneficial effects:

[0040] The present application provides a method for quickly obtaining a steady-state space electric field of an insulation layer of a self-powered detector, which estimates charge loss amounts generated by conduction current in different radius intervals of an insulator in a steady state according to the charge conservation law and the charge deposition amounts generated by irradiation current, and acquires the electric field intensity distribution, the charge density distribution and the potential distribution according to Ohm's law. BRIEF DESCRIPTION OF DRAWINGS

[0041] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiments. It should be understood that the following drawings only show some of the embodiments of the present application, and therefore should not be considered as limiting the scope. For those skilled in the art, other related drawings can also be obtained without creative labor.

[0042] Figure 1 The flow chart of the method for quickly obtaining the steady-state space electric field of the self-powered probe insulating layer of the present application.

[0043] Figure 2 The charge density distribution diagram of the self-powered probe insulating layer in the steady-state condition of the present application.

[0044] Figure 3 The electric field intensity distribution diagram of the self-powered probe insulating layer in the steady-state condition of the present application.

[0045] Figure 4 The potential distribution diagram of the self-powered probe insulating layer in the steady-state condition of the present application.

[0046] Figure 5 The system structure schematic diagram of the preferred embodiment of the present application for quickly obtaining the steady-state space electric field of the self-powered probe insulating layer.

[0047] Figure 6 The electronic device structure schematic diagram of the preferred embodiment of the present application. DETAILED DESCRIPTION

[0048] The exemplary embodiments of the present application will be described below with reference to the accompanying drawings, which include various details of the embodiments of the present application to help the understanding, and should be considered as merely exemplary. Therefore, those skilled in the art should recognize that various changes and modifications can be made to the embodiments described herein without departing from the scope and spirit of the present application. Also, in order to be clear and concise, the description below omits the description of the well-known functions and structures.

[0049] Obviously, the described embodiments are part of the embodiments of the present application, not all. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0050] It should be noted that the terminal involved in the embodiments of the present application can include, but is not limited to, a mobile phone, a personal digital assistant (PDA), a wireless handheld device, a tablet computer, a personal computer (PC), an MP3 player, an MP4 player, a wearable device (for example, smart glasses, a smart watch, a smart bracelet, etc.), a smart home device, and the like.

[0051] In addition, the term "and / or" in this paper is only a description of the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B, which can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in this paper generally represents that the front and rear associated objects are in an "or" relationship.

[0052] The application will be described in further detail below with reference to the drawings:

[0053] Referring to Figure 1 The application provides a method for quickly obtaining a steady-state space electric field of an insulating layer of a self-sustaining detector, comprising the following steps:

[0054] Step 1: obtaining the charge deposition amount Q generated by the irradiation current in different radius intervals per unit time of the detector dep (r~r+Δr).

[0055] Wherein, the charge deposition amount Q generated by the irradiation current per unit time dep (r~r+Δr) is calculated by Monte Carlo simulation software, that is, the radiation field information (including particle types, particle energy and momentum distribution, radiation field flux density, etc.) at the position where the self-sustaining neutron detector is located and the calculation detector information (geometric structure, material, etc. of the detector) are given as the input of the Monte Carlo program; the interaction of rays and the detector is simulated by the Monte Carlo program, and the charge deposition amount Q dep (r~r+Δr) generated by the irradiation current in different radius intervals of the detector is counted.

[0056] The radiation field information (including particle types, particle energy and momentum distribution, radiation field flux density, etc.) at the position where the self-sustaining neutron detector is located and the required parameters (geometric structure, material, etc. of the detector) can also be given as the input of the Monte Carlo program; the interaction of rays and the detector is simulated by the Monte Carlo program, and the irradiation current density at different positions of the detector is counted, and the charge deposition amount Q dep (r~r+Δr) generated by the irradiation current in different radius intervals of the detector is calculated by the following formula.

[0057] Q dep (r~r+Δr) = -2πL[J rad (r)*r - J rad (r+Δr)*(r+Δr)]

[0058] where L is the length of the detector;

[0059] Since irradiation can produce unstable nuclei, which due to long half-life, cannot be completely decayed in the current time step, in order to quickly calculate the steady state, it is necessary to count the charge deposition amount generated by the complete decay of these unstable nuclei.

[0060] Step 2: According to the law of conservation of charge, estimate the charge loss amount Q generated by the conduction current in the different radius intervals of the insulator per unit time in the steady state esp (r~r+Δr).

[0061] Since the self-powered detector structure is coaxial cylinder, and the length L of the detector is much larger than the radius R, it can be considered that the charge and current distribution is only related to r.

[0062] In the steady state, the free charge density in the insulator is constant, that is

[0063]

[0064] The above formula shows that the charge deposition amount generated by the irradiation current per unit time is the same in value and opposite in sign as the charge loss amount generated by the conduction current, that is

[0065] Q esp (r~r+Δr) = -Q dep (r~r+Δr)

[0066] Substituting the conduction current density, we can get

[0067] Q esp (r~r+Δr) = -J cond (r) S(r) + J cond (r+Δr) S(r+Δr)

[0068] = -2πL[J cond (r)*r - J cond (r+Δr)*(r+Δr)]

[0069] Step 3: Combine Ohm's law to establish the expression of the charge loss amount generated by the conduction current about the electric field intensity in the steady state.

[0070] According to Ohm's law, we have

[0071] J cond (r) = σE(r)

[0072] Thus, we get

[0073] Q esp (r) = -2πσL [E(r) * r - E(r + Δr) * (r + Δr)]

[0074] Where, J cond (r) is the conduction current density at radius r in steady state, E(r) is the electric field intensity at radius r in steady state, σ is the conductivity of the detector insulating layer, L is the length of the detector.

[0075] Divide the detector insulating layer along the radius into N cylindrical shells, the inner radius of the kth layer is denoted as r k , and the outer radius is denoted as r k+1 , then the inner radius of the detector insulating layer is r1, and the outer radius is r N+1 , where k = 1, 2, 3, …, N; the radius (r k , r k+1 ) interval is the amount of charge loss per unit time generated by the conduction current, denoted as Q k ; the electric field intensity at radius r k is denoted as E k , we get

[0076]

[0077] Step 4: Combine the boundary conditions to establish a linear equation set of the electric field intensity and the amount of charge loss per unit time generated by the conduction current, and arrange it in matrix form. Through matrix inversion method, the electric field intensity distribution, charge density distribution and potential distribution are obtained, and the three physical quantities can be determined by one of them. According to the relationship between potential and electric field, the potential difference ΔV of the inner and outer conductors can be expressed as

[0078]

[0079] The emitter and collector of the detector are both conductors, and their potentials are considered to be 0, i.e. ΔV = 0.

[0080] For the integral of the electric field, it can be expressed in two ways:

[0081] a. Direct numerical integration (Gaussian integration, rectangular integration, trapezoidal integration, etc.), and arrange the superposition coefficients of the electric field intensity at different radii as a supplementary equation. Here, take the trapezoidal integration as an example to illustrate.

[0082]

[0083] If the insulating layer is uniformly divided, it can be arranged as

[0084]

[0085] Combine the following formula

[0086]

[0087] Rearranging into matrix form, we get

[0088] AE = Q

[0089]

[0090] E = [E1, E2, ..., E N+1 ] T

[0091]

[0092] Therefore, the electric field intensity distribution can be obtained by inverting the matrix, i.e.

[0093] E = A -1 Q

[0094] According to Maxwell's equations, the relationship between the electric field strength and the free charge density within an insulator is:

[0095]

[0096] In the formula, ε represents the dielectric constant of the insulator.

[0097] Therefore, the free charge density distribution within the insulator can be obtained from the following equation.

[0098]

[0099] The electric field intensity distribution expression can be fitted and then directly differentiated; alternatively, the charge density distribution can be obtained through numerical differentiation. For example, using an intermediate difference scheme, it is assumed that (r k ,r k+1 The charge density over the radius interval is a constant ρ k ,So

[0100]

[0101]

[0102]

[0103] Similarly, the expression for electric potential is:

[0104]

[0105] The electric field intensity distribution expression can be fitted and then directly integrated; alternatively, the potential distribution can be obtained through numerical integration, for example, by using the trapezoidal rule.

[0106]

[0107]

[0108] b. According to the mirror method, the electric field can be obtained from the charge density distribution in the insulating layer, the steady-state charge density in the interval of radius (r k ,r k+1 ) is denoted as ρ k , and the electric field intensity E(r) in the steady-state case is

[0109]

[0110]

[0111]

[0112]

[0113] Since the electric field intensity E(r) is obtained by linear superposition of the charge density distribution ρ m (m = 1, 2, …, N) in the insulating layer, it can be arranged as:

[0114]

[0115]

[0116]

[0117] On this basis,

[0118]

[0119] Let

[0120] a k,m = 2πσL(e m+1,k r m+1 -e m,k r m )

[0121] After arrangement,

[0122] AX = Q

[0123] A = [a k,m ] N*N

[0124] X = [ρ m ] N*1

[0125] Q = [Q k ] N*1

[0126] The steady-state charge density distribution can be obtained by matrix inversion

[0127] X = A -1 Q

[0128] Based on the charge density distribution, the electric field strength distribution can be directly calculated, i.e.

[0129]

[0130] Based on the charge density distribution under the steady-state condition, the steady-state potential distribution calculation process is as follows:

[0131]

[0132]

[0133]

[0134]

[0135] The sensitivity of the detector under the steady-state electric field condition is counted:

[0136] The sensitivity of the detector under the steady-state electric field condition is calculated by the Monte Carlo simulation software, i.e. the radiation field information (including particle species, particle energy and momentum distribution, radiation field flux density, etc.) at the position where the self-powered neutron detector is located and the required parameters (detector geometry, material, calculated electric field distribution, etc.) are given as the input of the Monte Carlo program; the interaction between the rays and the detector under the condition of containing the electric field is simulated by the Monte Carlo program, and the sensitivity of the detector is counted.

[0137] The counting method of the sensitivity of the detector is to determine the position where the electric field of the insulating layer of the detector is 0 according to the calculated electric field distribution, and count the irradiation current at the position.

[0138] Since the total current of the detector is equal to the sum of the conduction current and the irradiation current, i.e. I = I R + I C

[0139] In the formula, I represents the total current of the detector, I R represents the irradiation current, and I C represents the conduction current.

[0140] For the position where the electric field is 0, the conduction current is 0, and the total current I is equal to the irradiation current I R .

[0141] The sensitivity S of the detector is equal to the ratio of the total current I to the neutron flux at the position of the detector, i.e.

[0142]

[0143] It should be noted that the electric field calculated by the method can be used as the initial value for iterative calculation of the electric field. Since the result of the fast calculation is close to the real electric field distribution, the calculation result of the method as the initial value input to the iterative calculation can further improve the accuracy of the calculation result with less calculation time.

[0144] Example 1:

[0145] This embodiment adopts a typical rhodium self-powered neutron detector with a length of 5 cm, and the detector size and material are shown in the following table.

[0146] Table 1 Detector size and material

[0147] Geometry Outer radius (cm) Material Density (g / cm 3 )]]> Emitting electrode 0.025 103 Rh]]> ​ 12.4 Insulating layer 0.04 Al2O3 2.9 Collecting electrode 0.065 Inconel 8.4

[0148] The radioactive source is a cylindrical isotropic radioactive source located at the surface of the emitter. A thermal neutron source at 600K is used, with a flux density of 1E13 n·cm -2 s -1 The resistivity of the insulating layer is set to 5E12 Ω·m. Since the rhodium detector is a delayed SPND, the charge deposition amount generated by the decay of all unstable nuclei needs to be counted during simulation.

[0149] The calculation results of the steady-state charge density distribution are shown in Figure 2 The fast calculation result only calculates once (10 6 instances, taking 2.3 minutes), so the statistical fluctuation is large, but the calculation time is extremely short; the iterative calculation result is iterated 43 times to reach the steady state (10 6 instances each time, total time 2.1 hours), the statistical error is low, but the calculation cost is high.

[0150] The calculation results of the steady-state electric field intensity distribution are shown in Figure 3 The fast calculation result is slightly larger than the iterative calculation result, but the overall trend is consistent. The maximum field strength of the iterative calculation steady state is 3.52×10 8 V / m, the maximum field strength of the fast calculation steady state is 4.05×10 8 V / m, and the relative error of the maximum field strength is 15.18%; the calculation results of the steady-state potential distribution are shown in Figure 4As shown, the fast calculation result is slightly larger than the iterative calculation result, but the overall trend is consistent. The steady-state maximum potential of the iterative calculation is 10.24 kV, the steady-state maximum potential of the fast calculation is 11.70 kV, and the relative error of the maximum potential is 16.68%. Therefore, by using the fast calculation method provided in the application, the charge, electric field and potential distribution of the insulating layer can be quickly calculated, and whether the electric field strength of the insulating layer exceeds the breakdown field strength can be estimated.

[0151] Embodiment 2

[0152] The embodiment 2 provided by the application is an embodiment of the system for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector provided by the application, as shown in the embodiment 2, Figure 5 The embodiment of the system comprises an acquisition module, an estimation module and a calculation module.

[0153] The acquisition module is used to acquire the radiation field information and the detector information at the position where the self-sustaining energy neutron detector is located, and simulate the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector;

[0154] The estimation module is used to estimate the charge loss amount generated by the conduction current in different radius intervals of the insulator per unit time under the steady-state condition according to the charge conservation law based on the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector;

[0155] The calculation module is used to establish an expression of the charge loss amount generated by the conduction current with respect to the electric field strength under the steady-state condition in combination with the Ohm's law, and acquire the electric field strength distribution, the charge density distribution and the potential distribution in combination with the boundary conditions.

[0156] It can be understood that the system for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector provided by the application corresponds to the method for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector provided by the above-mentioned embodiments, and the related technical features of the system for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector can refer to the related technical features of the method for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector, which will not be described here.

[0157] As shown in the embodiment 2, Figure 6 Another object of the application is to provide an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the steps of the method for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector.

[0158] The method for quickly obtaining the steady-state space electric field of the insulating layer of the self-sustaining energy detector comprises the following steps:

[0159] Obtain the radiation field information and the detector information at the position where the self-powered neutron detector is located, and simulate the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector;

[0160] Based on the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector, the charge loss amount generated by the conduction current in different radius intervals per unit time of the insulator in a steady state is estimated according to the charge conservation law;

[0161] The expression of the charge loss amount generated by the conduction current with respect to the electric field intensity in a steady state is established in combination with the Ohm's law, and the electric field intensity distribution, the charge density distribution and the potential distribution are obtained in combination with the boundary conditions.

[0162] The fourth object of the present application is to provide a computer readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method for quickly obtaining the steady-state space electric field of the insulating layer of the self-powered detector.

[0163] The method for quickly obtaining the steady-state space electric field of the insulating layer of the self-powered detector comprises the following steps:

[0164] Obtain the radiation field information and the detector information at the position where the self-powered neutron detector is located, and simulate the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector;

[0165] Based on the charge deposition amount generated by the irradiation current in different radius intervals per unit time of the detector, the charge loss amount generated by the conduction current in different radius intervals per unit time of the insulator in a steady state is estimated according to the charge conservation law;

[0166] The expression of the charge loss amount generated by the conduction current with respect to the electric field intensity in a steady state is established in combination with the Ohm's law, and the electric field intensity distribution, the charge density distribution and the potential distribution are obtained in combination with the boundary conditions.

[0167] Those skilled in the art will understand that embodiments of the present application can be provided as methods, systems or computer program products. Therefore, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage etc.) containing computer usable program code.

[0168] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flow or blocks. Figure 1 one or more flow or blocks.

[0169] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. Figure 1 one or more flow or blocks. Figure 1 one or more flow or blocks.

[0170] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flow or blocks. Figure 1 one or more flow or blocks.

[0171] Finally, it should be noted that the above-mentioned embodiments are merely used to illustrate the technical solutions of the present application, but not limit the technical solutions of the present application. Although the present application has been described in detail with reference to the above-mentioned embodiments, those skilled in the art should understand that the technical solutions of the present application can be modified or equivalent replaced without departing from the spirit and scope of the present application, and any modification or equivalent replacement should be covered in the protection scope of the claims of the present application.

Claims

1. A method for rapidly obtaining the steady-state spatial electric field of the insulating layer of a self-powered detector, characterized in that, include: By acquiring radiation field information and detector information at the location of the self-powered neutron detector, the amount of charge deposition generated by the irradiation current in different radius intervals of the detector per unit time can be simulated. Based on the amount of charge deposited by the irradiation current in different radius intervals of the detector per unit time, and according to the law of charge conservation, the amount of charge lost by the conduction current in different radius intervals of the insulator per unit time under steady state is estimated. By combining Ohm's law, we establish an expression for the charge loss caused by the conduction current in steady state with respect to the electric field strength. By combining the boundary conditions, we obtain the electric field strength distribution, charge density distribution, and electric potential distribution. Under the aforementioned steady-state condition, the calculation process for the charge loss caused by the conduction current within different radius intervals of the insulator per unit time is as follows: Under steady-state conditions, the free charge density within the insulator remains constant, i.e. This indicates that, per unit time, the amount of charge deposited by irradiation current and the amount of charge lost by conduction current are the same in value but opposite in sign. Substituting the conduction current density, we obtain in, This represents the amount of charge deposited per unit time generated by the irradiation current. It represents the amount of charge lost per unit time due to conduction current. Indicates the length of the detector. This represents the conduction current density at radius r in steady state; Under the aforementioned steady-state condition, the expression for the charge loss caused by the conduction current with respect to the electric field strength is as follows: in, Indicates the conductivity of the detector's insulating layer. Indicates the length of the detector. This represents the inner radius of the k-th cylindrical shell layer when the detector's insulating layer is uniformly divided into N cylindrical shell layers along the radius. Represents the outer radius, k = 1, 2, 3, ..., N; Represents radius The amount of charge loss per unit time due to conduction current within the interval. Represents radius The electric field strength at that location; The calculation process for the electric field intensity distribution, charge density distribution, and electric potential distribution under the aforementioned steady-state condition is as follows: By combining the boundary conditions, a set of linear equations is established for the electric field strength and the charge loss per unit time caused by the conduction current, and then rearranged into matrix form; Based on the relationship between electric potential and electric field, the potential difference between the inner and outer conductors... Represented as Both the detector's emitter and collector are conductors, and their potentials are assumed to be 0. ; By integrating the electric field, the electric field intensity distribution, charge density distribution, and electric potential distribution can be calculated.

2. The method for rapidly obtaining the steady-state spatial electric field of the insulating layer of a self-powered detector according to claim 1, characterized in that, The amount of charge deposition generated by the irradiation current within different radius intervals per unit time of the detector was obtained by Monte Carlo simulation software: Input the radiation field information and detector information at the location of the self-powered neutron detector into the Monte Carlo program; The interaction between the radiation and the detector was simulated using a Monte Carlo program, and the amount of charge deposited by the irradiation current in different radius ranges of the detector was statistically analyzed.

3. The method for rapidly obtaining the steady-state spatial electric field of the insulating layer of a self-powered detector according to claim 2, characterized in that, The radiation field information includes particle type, particle energy and momentum distribution, and radiation field flux density; the detector information includes detector geometry and materials.

4. The method for rapidly obtaining the steady-state spatial electric field of the insulating layer of a self-powered detector according to claim 1, characterized in that, The method for integrating the electric field is either numerical integration or the method of images.

5. A system for rapidly acquiring the steady-state spatial electric field of the insulating layer of a self-powered detector, characterized in that, The steps for implementing the method of claim 1 include: The acquisition module is used to acquire radiation field information and detector information at the location of the self-powered neutron detector, and to simulate the amount of charge deposition generated by the irradiation current in different radius intervals of the detector per unit time. The estimation module is used to estimate the amount of charge loss caused by conduction current in different radius intervals of the insulator under steady-state conditions, based on the amount of charge deposition generated by the irradiation current in different radius intervals of the detector per unit time and according to the law of conservation of charge. The calculation module is used to establish an expression for the charge loss generated by the conduction current with respect to the electric field strength under steady-state conditions, based on Ohm's law, and to obtain the electric field strength distribution, charge density distribution, and electric potential distribution by combining boundary conditions.

6. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 4.

7. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 4.

Images