CN110579790B - Method for evaluating electrostatic effect of self-powered neutron detector - Google Patents
Method for evaluating electrostatic effect of self-powered neutron detector Download PDFInfo
- Publication number
- CN110579790B CN110579790B CN201910982838.8A CN201910982838A CN110579790B CN 110579790 B CN110579790 B CN 110579790B CN 201910982838 A CN201910982838 A CN 201910982838A CN 110579790 B CN110579790 B CN 110579790B
- Authority
- CN
- China
- Prior art keywords
- self
- insulating layer
- neutron detector
- charge
- powered neutron
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/005—Details of radiation-measuring instruments calibration techniques
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
A method for evaluating electrostatic effect of a self-powered neutron detector comprises the following steps: (1) obtaining the charge deposition rate V of the insulating layer at the zero moment by Monte Carlo program simulationq,dep,0(r, t), and at this time the charge loss rate Vq,out,0(r, t) is zero; (3) from Vq,dep,0(r, t) and Vq,out,0(r, t) obtaining a net deposition rate V of the charge at time zeroq,0(r, t); (4) setting the time step as delta t, and taking the charge deposition rate and the charge loss rate of the previous time point as initial values of the current time point to obtain the charge distribution Q (r, t) of the insulating layer at the current time; (5) obtaining the charge density rho of the insulating layer according to Q (r, t) and the volume of the insulating layer0(r, t) and fitting to obtain an analytic expression rho (r, t) thereof; (6) substituting rho (r, t) into a Poisson equation and solving to obtain the potential and the electric field of the insulating layer; (7) adding the electric field into the Monte Carlo simulation process to obtain new net deposition rate of electric charge at the current time point, and repeating the steps (4) - (7) until stable net deposition rate of electric charge V is obtainedq(r, t); (8) if VqAnd (r, t) outputting the electric field and the electric potential if the convergence requirement is met, otherwise, repeating the steps.
Description
Technical Field
The invention belongs to the technical field of neutron detection, and particularly relates to an evaluation method for electrostatic effect of a self-powered neutron detector.
Background
A self-powered neutron detector generally consists of an emitter, an insulating layer, and a collector. The emitter is generally made of a material with a large neutron capture cross section, and secondary electrons generated in the emitter or an insulating layer by rays generated by reaction of neutrons and the emitter material reach the collector to generate current, so that neutron flux is measured. The material of the insulating layer generally has good electrical insulating properties, and is resistant to high temperature and high pressure and strong irradiation. The collector is used for receiving electrons emitted from the emitter or the insulating layer, and the material has the characteristics of high temperature and high pressure resistance and corrosion resistance, and has small deformation at high temperature so as to ensure that the self-powered neutron detector has a stable mechanical structure.
Electrons moving to the insulating layer may stay in the insulating layer due to insufficient energy or the electrons in the insulating layer are ionized to leave the original position and leave positive charges, and the charges cannot flow away in time due to poor conductivity of the insulating layer and are accumulated, so that the generated electric field influences the transportation of subsequent electrons in the insulating layer, and further influences the detection signal, namely the electrostatic effect of the self-powered neutron detector. Therefore, the electrostatic effect should be taken into account when building a calculation model for the response of the self-powered neutron detector.
At present, scholars at home and abroad make a lot of theoretical researches and experimental researches on the response of a self-powered neutron detector, so that a lot of achievements are obtained, and meanwhile, the scholars have some defects: 1) the electrostatic effect is considered in most theoretical model researches, but the distribution of the electrostatic charge in the insulating layer comes from the assumption, the resistivity of the material of the insulating layer and the real ray flux density near the detector in the pile are not considered, and therefore, the electrostatic effect is not considered in the assumed form to have any practical significance; 2) the response of a self-powered neutron detector can be well simulated by establishing a numerical calculation model, but the relevant literature does not consider the electrostatic effect in the detector numerical model, so that the calculation model is incomplete.
Disclosure of Invention
In order to solve the problems in the prior art, an object of the present invention is to provide a method for evaluating an electrostatic effect of a self-powered neutron detector, which quantitatively describes an electrostatic field and an electrostatic potential existing in an insulating layer due to electrostatic charge deposition, and can quantitatively study the influence of the electrostatic effect on the performance of the self-powered neutron detector, so as to better guide the application of the self-powered neutron detector in practice.
In order to achieve the above purpose, the invention adopts the following technical scheme:
a method for evaluating electrostatic effect of a self-powered detector comprises the following steps:
step 1: providing the neutron-photon source distribution near a self-powered neutron detector in an actual reactor as the input of a Monte Carlo program; by the Monte Carlo programSequentially simulating the particle transport process in the self-powered neutron detector to obtain the charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of an emitter of the self-powered neutron detector in unit timeq，dep(r，t)；
Step 2: average neutron flux density phi of emitter of neutron detector with actual self-power in reactorn，emiAverage photon flux density phiγ，emiAnd step 1, obtaining charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of the emitter of the self-powered neutron detector in unit timeq，dep(r, t) is substituted into the formula (1) to obtain the charge deposition rate distribution V corresponding to different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，dep，0(r，t)；
Vq，dep，0(r，t)＝(φn，emi+φγ，emi)·vq，dep(r，t) (1)
And step 3: the electric charge quantity accumulated in the insulating layer of the self-powered neutron detector at the zero moment is zero, and no electric field exists, so that the charge loss rate distribution V caused by the existence of the electric field at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，out，0(r, t) is zero;
and 4, step 4: v obtained in step 2 and step 3q，dep，0(r, t) and Vq，out，0(r, t) is substituted into the formula (2) to obtain the net deposition rate distribution V of charges at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，0(r，t)；
Vq，0(r，t)＝Vq，dep，0(r，t)-Vq，out，0(r，t) (2)
Vq，cur(r，t)＝Vq，dep，pre(r，t)-Vq，out，pre(r，t) (2-1)
Vq(r，t)＝Vq，dep(r，t)-Vq，out(r，t) (2-2)
And 5: setting the time step length as delta t, and obtaining the unitThe charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the previous time point of each time stepq，dep，pre(r, t) and a charge loss rate distribution Vq，out，pre(r, t) is taken as an initial value of the current time point, linear change of the total charge deposition rate and the total charge loss rate in a single time step is assumed, charge self-balancing iterative calculation in the single time step is carried out, and stable charge deposition rate distribution V at different radius positions in an insulating layer of the self-powered neutron detector at the current time point is obtainedq，dep(r, t) and a charge loss rate distribution Vq，out(r, t) by the following steps: (a) will Vq，dep，pre(r, t) and Vq，out，pre(r, t) is substituted into the formula (2-1) to obtain the charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq，our(r, t) and distributing the net deposition rate V of the electric charges at different radius positions in the insulating layer of the energy-saving neutron detector at a time point before the current time stepq(r, t- Δ t), charge distribution Q (r, t- Δ t), and Vq，cur(r, t) and the time step delta t are substituted into a formula (3) to obtain the charge distribution Q (r, t) at different radius positions in the insulating layer of the self-powered neutron detector at the current time point; (b) obtaining the charge density distribution rho at different radius positions in the insulating layer of the self-powered neutron detector at the current time point according to the Q (r, t) obtained in the step (a) and the volume of the insulating layer of the self-powered neutron detector0(r, t) and fitting to obtain an analytic expression rho (r, t) thereof; (c) substituting rho (r, t) obtained in the step (b) into a Poisson equation (4) and solving to obtain the potential distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point
And obtaining the electric field distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point by the formula (5)
(d) Subjecting the product obtained in step (c)
Adding the new charge deposition rate distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point by adopting the methods of the step 1 and the step 2 in the Monte Carlo simulation process; (e) due to the existence of the electric field, current can be formed in the insulating layer of the self-powered neutron detector, and the current density distribution at different radius positions
Obtained by the formula (6); (f) charge loss rate distribution V caused by electric field at different radius positions in insulating layer of self-powered neutron detector at current time pointq，out，cur(r, t) is obtained from formula (7); (g) if the charge loss rate distribution V obtained by two times of iterative calculationq，out(r，t，i-1)、Vq，out(r, t, i) satisfying the condition (8), converging the self-balancing iterative computation, otherwise, repeating the steps (a) - (f);
Q(r，t)＝Q(r，t-Δt)+0.5·Δt·[Vq(r，t)+Vq(r，t-Δt)] (3)
wherein, the dielectric constant is the relative dielectric constant of the insulating layer material of the self-powered neutron detector;
wherein sigma is the conductivity of the insulating layer material of the self-powered neutron detector at the corresponding temperature;
in the formula (I), the compound is shown in the specification,
the current density distribution at a radius position behind different radius intervals in an insulating layer of the self-powered neutron detector is realized; s (r) -2 pi rL is the axial sectional area of the insulating layer of the self-powered neutron detector at different radius positions, and L is the length of the self-powered neutron detector;
|Vq，out(t，i)-Vq，out(t，i-1)|＜0 (8)
in the formula, Vq，out(t，i)＝∑Vq，out(r, t, i) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the ith iterative calculation of charge self-balancing calculation in a single time step; vq，out(t，i-1)＝∑Vq，out(r, t, i-1) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the i-1 th iteration calculation of charge self-balancing calculation in a single time step;0calculating a convergence residual for self-balancing iteration;
step 6: the stable charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time point obtained in the step 5q，dep(r, t) and a charge loss rate distribution Vq，out(r, t) is substituted into the formula (2-2) to obtain stable charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq(r，t)。
If V is obtained from step 6q(r, t) satisfies the condition (9), which indicates that the charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is stable, the global balance iterative computation is converged, and the potential distribution at different radius positions of the insulating layer of the self-powered neutron detector is output
And electric field distribution
Otherwise, repeating the step 5 and the step 6Continuously carrying out global balance iterative calculation;
|Vq(t)|＜1 (9)
in the formula, Vq(t)＝∑Vq(r, t) is the stable net deposition rate of total charge in the insulating layer of the self-powered neutron detector at the current point in time;1a convergence residual is calculated for the global balanced iteration.
Compared with the traditional method, the method has the following advantages: (1) the charge accumulation process in the insulating layer of the self-powered neutron detector can be simulated, and stable charge distribution in the insulating layer is obtained finally without approximating the charge distribution in an assumed mode; (2) the real ray flux density in the pile and the resistivity of the insulating layer material at a specific temperature can be considered, so that the calculation result is more accurate.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a self-balancing charge calculation flow diagram for an insulating layer of a self-powered neutron detector in a single time step.
FIG. 3 is a graphical representation of net deposition rate of charge in an insulating layer of a self-powered neutron detector as a function of time.
FIG. 4 is a schematic diagram of a self-powered neutron detector.
FIG. 5 is a graph of charge deposition rate distribution at different radii of an insulating layer per average neutron flux density of an emitter of a self-powered neutron detector at time zero per unit time.
FIG. 6 is a graph of charge deposition rate distribution at different radial locations in an insulating layer of an energy-self neutron detector at time zero.
FIG. 7 is a graph of net deposition rate distribution of charge at different radial locations in an insulating layer of an energy-self neutron detector at time zero.
FIG. 8 is a graph of net deposition rate distribution of charge at different radial locations in an insulating layer of an energy-efficient neutron detector at a first point in time.
FIG. 9 is a charge distribution within an insulating layer of a self-powered neutron detector at a first point in time.
FIG. 10 is a graph of electric field distribution at different radial locations in an insulating layer of an energy-efficient neutron detector at a first point in time.
FIG. 11 is a graph of potential distribution at different radial locations in an insulating layer of an energy-efficient neutron detector at a first point in time.
FIG. 12 is a graph of net deposition rate of charge in an insulating layer of an energy-self-powered neutron detector over time.
FIG. 13 is a steady state charge distribution at different radial locations in an insulating layer of an energy-efficient neutron detector.
FIG. 14 is a steady state electric field distribution at different radial locations in an insulating layer of an energy-efficient neutron detector.
FIG. 15 is a steady state potential distribution at different radial positions in an insulating layer of an energy-efficient neutron detector.
Detailed Description
The invention is described in further detail below with reference to the following figures and specific embodiments.
The invention discloses an evaluation method of electrostatic effect of a self-powered detector, which comprises the following steps: the specific process is shown in figure 1.
Step 1: the neutron-photon source distribution in the actual reactor near the self-powered neutron detector is given as input to the monte carlo procedure. Simulating a particle transport process in the self-powered neutron detector by a Monte Carlo program to obtain charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of an emitter of the self-powered neutron detector in unit timeq，dep(r，t)；
Step 2: average neutron flux density phi of emitter of neutron detector with actual self-power in reactorn，emiAverage photon flux density phiγ，emiAnd step 1, obtaining charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of the emitter of the self-powered neutron detector in unit timeq，dep(r, t) is substituted into the formula (1) to obtain the charge deposition rate distribution V corresponding to different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，dep，0(r，t)；
Vq，dep，0(r，t)＝(φn，emi+φγ，emi)·vq，dep(r，t) (1)
And step 3: the electric charge quantity accumulated in the insulating layer of the self-powered neutron detector at the zero moment is zero, and no electric field exists, so that the charge loss rate distribution V caused by the existence of the electric field at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，out，0(r, t) is zero;
and 4, step 4: v obtained in step 2 and step 3q，dep，0(r, t) and Vq，out，0(r, t) is substituted into the formula (2) to obtain the net deposition rate distribution V of charges at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，0(r，t)；
Vq，0(r，t)＝Vq，dep，0(r，t)-Vq，out，0(r，t) (2)
Vq，cur(r，t)＝Vq，dep，pre(r，t)-Vq，out，pre(r，t) (2-1)
Vq(r，t)＝Vq，dep(r，t)-Vq，out(r，t) (2-2)
And 5: as shown in FIG. 2, the time step is set to be delta t, and the charge deposition rate distribution V at different radius positions in the insulating layer of the energy-saving neutron detector at the time point before the single time stepq，dep，pre(r, t) and a charge loss rate distribution Vq，out，pre(r, t) is taken as an initial value of the current time point, linear change of the total charge deposition rate and the total charge loss rate in a single time step is assumed, charge self-balancing iterative calculation in the single time step is carried out, and stable charge deposition rate distribution V at different radius positions in an insulating layer of the self-powered neutron detector at the current time point is obtainedq，dep(r, t) and a charge loss rate distribution Vq，out(r, t) by the following steps: (a) will Vq，dep，pre(r, t) and Vq，out，pre(r, t) is substituted into the formula (2-1) to obtain the charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq，cur(r, t) and applying a time point before the current time step to the energized neutronsNet deposition rate distribution of charge V at different radial positions in the detector's insulating layerq(r, t- Δ t), charge distribution Q (r, t- Δ t), and Vq，cur(r, t) and the time step delta t are substituted into a formula (3) to obtain the charge distribution Q (r, t) at different radius positions in the insulating layer of the self-powered neutron detector at the current time point; (b) obtaining the charge density distribution rho at different radius positions in the insulating layer of the self-powered neutron detector at the current time point according to the Q (r, t) obtained in the step (a) and the volume of the insulating layer of the self-powered neutron detector0(r, t) and fitting the equation by adopting a polynomial (3-1) to obtain an analytic expression rho (r, t); (c) substituting rho (r, t) obtained in the step (b) into a Poisson equation (4) and solving to obtain the potential distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time
And obtaining the electric field distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time according to the formula (5)
(d) Subjecting the product obtained in step (c)
Adding the new charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time point obtained in the step 1 and the step 2 in the Monte Carlo simulation processq，dep(r, t); (e) due to the existence of the electric field, current can be formed in the insulating layer of the self-powered neutron detector, and the current density distribution at different radius positions
Obtained by the formula (6); (f) charge loss rate distribution V caused by electric field at different radius positions in insulating layer of self-powered neutron detector at current time pointq，out(r, t) is obtained from formula (7); (g) if V is obtained by two times of iterative calculationq，out(r, t) satisfies the condition (8), the self-balancing iterative computation converges, otherwise steps (a) - (f) are repeated.
Q(r，t)＝Q(r，t-Δt)+0.5·Δt·[Vq(r，t)+Vq(r，t-Δt)] (3)
In the formula, aiAnd biIs a polynomial coefficient; n and m are the maximum order of the polynomial.
Wherein the dielectric constant of the insulating layer material of the self-powered neutron detector is shown in the formula.
In the formula (I), the compound is shown in the specification,
the current density distribution at different radius positions in an insulating layer of the self-powered neutron detector is realized; sigma is the conductivity of the insulating layer material of the self-powered neutron detector at the corresponding temperature;
the distribution of electric fields at different radius positions in an insulating layer of the neutron detector is self-powered.
In the formula (I), the compound is shown in the specification,
for self-powered neutronsCurrent density distribution at a radius position behind different radius intervals in the insulating layer of the detector; and S (r) -2 pi rL is the axial sectional area of the insulating layer of the self-powered neutron detector at different radius positions, and L is the length of the self-powered neutron detector.
|Vq，out(t，i)-Vq，out(t，i-1)|＜0 (8)
In the formula, Vq，out(t，i)＝∑Vq，out(r, t, i) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the ith iterative calculation of charge self-balancing calculation in a single time step; vq，out(t，i-1)＝∑Vq，out(r, t, i-1) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the i-1 th iteration calculation of charge self-balancing calculation in a single time step;0the convergence residual is calculated for self-balancing iterations.
The specific process for solving the poisson equation is as follows: since the neutron detector is in a cylindrical geometry, the potential in the insulating layer
Related only to the radius r, and thus Poisson's equation (4) can be written as equation (8-1)
Substituting formula (3-1) for formula (8-1) to obtain formula (8-2)
In the formula, A and B are constant coefficients.
To obtain finally
General solution of (1):
from the boundary conditions:
in the formula: r isminIs the inner radius of the insulating layer, rmaxThe outer radius of the insulating layer.
Obtaining:
wherein:
solving to obtain the coefficient:
so the solution is determined as:
the electric field strength is therefore:
step 6: the stable charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time point obtained in the step 5q，dep(r, t) and a charge loss rate distribution Vq，out(r, t) is substituted into the formula (2-2) to obtain stable charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq(r，t)。
As indicated by t in FIG. 35Shown if V is obtained from step 6q(r, t) satisfies the condition (9), which indicates that the charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is stable, the global balance iterative computation is converged, and the potential distribution at different radius positions of the insulating layer of the self-powered neutron detector is output
And electric field distribution
Otherwise, repeating the step 5 and the step 6 to continue the global balance iterative computation.
|Vq(t)|＜1 (9)
In the formula, Vq(t)＝∑Vq(r, t) is the stable net deposition rate of total charge in the insulating layer of the self-powered neutron detector at the current point in time;1a convergence residual is calculated for the global balanced iteration.
Examples of the applications
Taking a vanadium self-powered neutron detector as an example, the structure and the size of the detector are shown in fig. 4. The emitter material is natural vanadium and has a density of 6.11g/cm3. The insulating layer is made of alumina and has a density of 1.9g/cm3. The collector material is Inconel-600, and the density is 8.51g/cm3. Resistivity of the insulation layer the resistivity of the insulation layer was determined according to the literature [ Zareen Khan Abdul Jalil Khan, Mohd Idris Taib, Izhar Abu Husin,&Nurfarhana Ayuni(2010).Comparison study on in-core neutron detector for online neutron flux mapping of research and power reactor.RnD Seminar 2010：Research and Development Seminar2010，Malaysia.http：//inis.iaea.org/search/search.aspx？orig_q＝RN：43056490]is set to 6 x 1010Ω·m。
The neutron-photon source distribution near the self-powered neutron detector in the actual reactor adopts that the energy of the outer surface of the collector is 0.0253A monoenergetic isotropic neutron surface source of eV as input to the monte carlo program. Simulating the particle transport process in the self-powered neutron detector by a Monte Carlo program to obtain the charge deposition rate distribution v at different radiuses of the insulating layer under the unit average neutron flux density of the emitter of the self-powered neutron detector in unit timeq，dep(r, t) as shown in FIG. 5.
Neutron flux density phi of emitter of actual self-powered neutron detector in reactorn，emiGet 1013n·cm-2·s-1. By phin，emiAnd vq，dep(r, t) obtaining the charge deposition rate distribution V corresponding to different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，dep，0(r, t) as shown in FIG. 6.
The electric charge quantity accumulated in the insulating layer of the self-powered neutron detector at the zero moment is zero, and no electric field exists, so that the charge loss rate distribution V caused by the existence of the electric field at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，out，0(r, t) is zero.
From Vq，dep，0(r, t) and Vq，out，0(r, t) obtaining the net deposition rate distribution V of electric charges at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，0(r, t) as shown in FIG. 7.
Setting the time step length to be 0.02s, and distributing the charge deposition rate V at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，dep，0(r, t) and a charge loss rate distribution Vq，out，0(r, t) is taken as an initial value of the current time point, linear change of the total charge deposition rate and the total charge loss rate in a single time step is assumed, charge self-balancing iterative calculation in the single time step is carried out, and stable charge deposition rate distribution V at different radius positions in an insulating layer of the self-powered neutron detector at the current time point is obtainedq，dep(r, t) and a charge loss rate distribution Vq，out(r, t). Charge density distribution rho (r, t) at different radius positions in an insulating layer for a self-powered neutron detector in charge self-balancing iterative computation within a single time stepThe parameters n and m in the polynomial (3-1) in this example take 5 and 1, respectively.
Stable charge deposition rate distribution V at different radius positions in insulating layer of neutron detector powered by current time pointq，dep(r, t) and a charge loss rate distribution Vq，out(r, t) obtaining a stable net deposition rate profile V of chargeq(r, t) as shown in FIG. 8. At this time, the charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is shown in fig. 9, and the corresponding electric field distribution and electric potential distribution are shown in fig. 10 and 11, respectively. It can be seen that the charge density in the insulating layer of the self-powered neutron detector gradually decreases with the increase of the radius, and the electric field and the electric potential are smaller. The potential has a maximum of about 1.086kV, so electrons from the emitter or collector must have an energy greater than 1.086keV to be able to pass through the insulating layer.
Stable net deposition rate distribution V of electric charge at different radius positions in the insulating layer of the neutron detector from power if the current time pointq(r, t) satisfies the condition (9), which indicates that the charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is stable, the global balance iterative computation is converged, otherwise, the steps are repeated, and the global balance iterative computation is continued.
As shown in fig. 12, through the global balance iterative calculation, the net deposition rate of charge in the insulating layer of the self-powered neutron detector finally approaches zero, which indicates that the charge distribution tends to be stable. The steady-state charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is finally obtained as shown in fig. 13, and the corresponding electric field distribution and electric potential distribution are respectively shown in fig. 14 and fig. 15. At the moment, the charge density in the insulating layer of the self-powered neutron detector is increased, the insulating layer presents a distribution with a high middle and two low sides, and the electric field and the electric potential are increased accordingly. The potential has a maximum of about 486.816kV and therefore electrons from the emitter or collector must have an energy of more than 486.816keV to be able to penetrate the insulating layer. It follows that electrostatic effects do affect the transport of electrons in insulators, and that accurate potential and electric field distributions must be obtained by iterative calculations.
While the invention has been described with reference to embodiments, it is understood that the terminology used is intended to be in the nature of words of description and illustration, rather than of limitation. As the present invention may be embodied in several forms without departing from the spirit or essential characteristics thereof, it should also be understood that the above-described embodiments are not limited by any of the details of the foregoing description, but rather should be construed broadly within its spirit and scope as defined in the appended claims, and therefore all changes and modifications that fall within the meets and bounds of the claims, or equivalences of such meets and bounds are therefore intended to be embraced by the appended claims.
Claims (2)
1. A self-powered neutron detector electrostatic effect assessment method is characterized by comprising the following steps: the method comprises the following steps:
step 1: providing the neutron-photon source distribution near a self-powered neutron detector in an actual reactor as the input of a Monte Carlo program; simulating a particle transport process in the self-powered neutron detector by a Monte Carlo program to obtain charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of an emitter of the self-powered neutron detector in unit timeq，dep(r，t)；
Step 2: average neutron flux density phi of emitter of neutron detector with actual self-power in reactorn，emiAverage photon flux density phiγ，emiAnd step 1, obtaining charge deposition rate distribution v at different radius positions of an insulating layer of the self-powered neutron detector under the unit average neutron and photon flux density of the emitter of the self-powered neutron detector in unit timeq，dep(r, t) is substituted into the formula (1) to obtain the charge deposition rate distribution V corresponding to different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，dep，0(r，t)；
Vq，dep，0(r，t)＝(φn，emi+φγ，emi)·vq，dep(r，t) (1)
And step 3: the electric charge quantity accumulated in the insulating layer of the self-powered neutron detector at the zero moment is zero, and no electric field exists, so that the charge loss rate distribution V caused by the existence of the electric field at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，out，0(r, t) is zero;
and 4, step 4: v obtained in step 2 and step 3q，dep，0(r, t) and Vq，out，0(r, t) is substituted into the formula (2) to obtain the net deposition rate distribution V of charges at different radius positions in the insulating layer of the self-powered neutron detector at the zero momentq，0(r，t)；
Vq，0(r，t)＝Vq，dep，0(r，t)-Vq，out，0(r，t) (2)
Vq，cur(r，t)＝Vq，dep，pre(r，t)-Vq，out，pre(r，t) (2-1)
Vq(r，t)＝Vq，dep(r，t)-Vq，out(r，t) (2-2)
And 5: setting the time step as delta t, and setting the charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the time point before the single time stepq，dep，pre(r, t) and a charge loss rate distribution Vq，out，pre(r, t) is taken as an initial value of the current time point, linear change of the total charge deposition rate and the total charge loss rate in a single time step is assumed, charge self-balancing iterative calculation in the single time step is carried out, and stable charge deposition rate distribution V at different radius positions in an insulating layer of the self-powered neutron detector at the current time point is obtainedq，dep(r, t) and a charge loss rate distribution Vq，out(r, t) by the following steps: (a) will Vq，dep，pre(r, t) and Vq，out，pre(r, t) is substituted into the formula (2-1) to obtain the charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq，cur(r, t) and distributing the net deposition rate V of the electric charges at different radius positions in the insulating layer of the energy-saving neutron detector at a time point before the current time stepq(r, t- Δ t), charge distribution Q (r, t- Δ t), and Vq，cur(r, t) and the time step delta t are substituted into a formula (3) to obtain the charge distribution Q (r, t) at different radius positions in the insulating layer of the self-powered neutron detector at the current time point; (b) q (r, t) obtained in step (a) and an insulating layer body of a self-powered neutron detectorObtaining charge density distribution rho of different radius positions in the insulating layer of the self-powered neutron detector at the current time point through product0(r, t) and fitting to obtain an analytic expression rho (r, t) thereof; (c) substituting rho (r, t) obtained in the step (b) into a Poisson equation (4) and solving to obtain the potential distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point
And obtaining the electric field distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point by the formula (5)
(d) Subjecting the product obtained in step (c)
Adding the new charge deposition rate distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point by adopting the methods of the step 1 and the step 2 in the Monte Carlo simulation process; (e) due to the existence of the electric field, current can be formed in the insulating layer of the self-powered neutron detector, and the current density distribution at different radius positions
Obtained by the formula (6); (f) charge loss rate distribution V caused by electric field at different radius positions in insulating layer of self-powered neutron detector at current time pointq，out，cur(r, t) is obtained from formula (7); (g) if the charge loss rate distribution V obtained by two times of iterative calculationq，out，cur(r，t，i-1)、Vq，out，cur(r, t, i) satisfying the condition (8), converging the self-balancing iterative computation, otherwise, repeating the steps (a) - (f);
Q(r，t)＝Q(r，t-Δt)+0.5·Δt·[Vq，cur(r，t)+Vq(r，t-Δt)] (3)
wherein, the dielectric constant is the relative dielectric constant of the insulating layer material of the self-powered neutron detector;
wherein sigma is the conductivity of the insulating layer material of the self-powered neutron detector at the corresponding temperature;
in the formula (I), the compound is shown in the specification,
the current density distribution at a radius position behind different radius intervals in an insulating layer of the self-powered neutron detector is realized; s (r) -2 pi rL is the axial sectional area of the insulating layer of the self-powered neutron detector at different radius positions, and L is the length of the self-powered neutron detector;
|Vq，out，cur(t，i)-Vq，out，cur(t，i-1)|＜0 (8)
in the formula, Vq，out，cur(t，i)＝∑Vq，out，cur(r, t, i) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the ith iterative calculation of charge self-balancing calculation in a single time step; vq，out，cur(t，i-1)＝∑Vq，out，cur(r, t, i-1) is the total charge loss rate in the insulating layer of the self-powered neutron detector at the current time point obtained by the i-1 th iteration calculation of charge self-balancing calculation in a single time step;0calculating a convergence residual for self-balancing iteration;
step 6: obtained in step 5Stable charge deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq，dep(r, t) and a charge loss rate distribution Vq，out(r, t) is substituted into the formula (2-2) to obtain stable charge net deposition rate distribution V at different radius positions in the insulating layer of the self-powered neutron detector at the current time pointq(r, t); if V is obtained from step 6q(r, t) satisfies the condition (9), which indicates that the charge distribution Q (r, t) in the insulating layer of the self-powered neutron detector is stable, the global balance iterative computation is converged, and the potential distribution at different radius positions of the insulating layer of the self-powered neutron detector is output
And electric field distribution
Otherwise, repeating the step 5 and the step 6 to continue the global balance iterative computation;
|Vq(t)|＜1 (9)
in the formula, Vq(t)＝∑Vq(r, t) is the stable net deposition rate of total charge in the insulating layer of the self-powered neutron detector at the current point in time;1a convergence residual is calculated for the global balanced iteration.
2. The method for evaluating electrostatic effects of a self-powered neutron detector of claim 1, wherein: the stable net charge deposition rate distribution at different radius positions in the insulating layer of the self-powered neutron detector at the current time point is obtained by subtracting the charge loss rate distribution caused by the existence of the electric field from the charge deposition rate distribution caused by the rays at the different radius positions in the insulating layer of the self-powered neutron detector at the current time point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910982838.8A CN110579790B (en) | 2019-10-16 | 2019-10-16 | Method for evaluating electrostatic effect of self-powered neutron detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910982838.8A CN110579790B (en) | 2019-10-16 | 2019-10-16 | Method for evaluating electrostatic effect of self-powered neutron detector |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110579790A CN110579790A (en) | 2019-12-17 |
| CN110579790B true CN110579790B (en) | 2020-12-08 |
Family
ID=68814829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910982838.8A Active CN110579790B (en) | 2019-10-16 | 2019-10-16 | Method for evaluating electrostatic effect of self-powered neutron detector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110579790B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003139861A (en) * | 2001-11-06 | 2003-05-14 | Toshiba Corp | Detector for monitoring environment at local region |
| CN104900280A (en) * | 2015-04-09 | 2015-09-09 | 中国核动力研究设计院 | Signal delay elimination method for silver self-powered detector based on H infinity filtering |
| CN105373698A (en) * | 2015-10-14 | 2016-03-02 | 中科华核电技术研究院有限公司 | Method for eliminating signal time delay of self-powered neutron detector |
| CN108828650A (en) * | 2018-08-08 | 2018-11-16 | 中国原子能科学研究院 | A kind of design method of the shielding thickness of neutron measurement device |
| CN109671513A (en) * | 2018-12-27 | 2019-04-23 | 中国原子能科学研究院 | The continuously adjustable neutron howitzer structure of the angle of divergence and its scaling method |
| CN109903867A (en) * | 2019-02-28 | 2019-06-18 | 西安交通大学 | A kind of method that circuit composition and its parameter are eliminated in determining self-power neutron detector delay |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004073326A2 (en) * | 2003-02-09 | 2004-08-26 | Structured Materials Industries, Inc. | Smart portable detector and microelectronic radiation detector |
-
2019
- 2019-10-16 CN CN201910982838.8A patent/CN110579790B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003139861A (en) * | 2001-11-06 | 2003-05-14 | Toshiba Corp | Detector for monitoring environment at local region |
| CN104900280A (en) * | 2015-04-09 | 2015-09-09 | 中国核动力研究设计院 | Signal delay elimination method for silver self-powered detector based on H infinity filtering |
| CN105373698A (en) * | 2015-10-14 | 2016-03-02 | 中科华核电技术研究院有限公司 | Method for eliminating signal time delay of self-powered neutron detector |
| CN108828650A (en) * | 2018-08-08 | 2018-11-16 | 中国原子能科学研究院 | A kind of design method of the shielding thickness of neutron measurement device |
| CN109671513A (en) * | 2018-12-27 | 2019-04-23 | 中国原子能科学研究院 | The continuously adjustable neutron howitzer structure of the angle of divergence and its scaling method |
| CN109903867A (en) * | 2019-02-28 | 2019-06-18 | 西安交通大学 | A kind of method that circuit composition and its parameter are eliminated in determining self-power neutron detector delay |
Non-Patent Citations (1)
| Title |
|---|
| 铑自给能探测器电流信号计算方法及影响因素;李树成等;《核电子学与探测技术》;20110531;565-567 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110579790A (en) | 2019-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Highly stable upwind FEM for solving ionized field of HVDC transmission line | |
| Liu et al. | An efficient model to simulate stable glow corona discharges and their transition into streamers | |
| CN107729608B (en) | Short-gap gas discharge numerical simulation method based on time domain spectral element method | |
| Lu et al. | Characteristics of Trichel pulse parameters in negative corona discharge | |
| Zhou et al. | Review of charge accumulation on spacers of gas insulated equipment at DC stress | |
| Jaworski et al. | Modification of the electron energy distribution function during lithium experiments on the National Spherical Torus Experiment | |
| CN103823114A (en) | Apparatus and method for measuring radiation induced conductivity of medium material at different temperatures | |
| Rodríguez‐Serna et al. | Finite‐element‐analysis models for numerical simulation of partial discharges in spherical cavities within solid dielectrics: a review and a novel method | |
| CN110765618A (en) | Method for calculating response current of self-powered neutron detector in pressurized water reactor | |
| Lu et al. | The effect of environmental temperature on negative corona discharge under the action of photoionization | |
| He et al. | Simulation of positive streamer propagation in an air gap with a GFRP composite barrier | |
| CN110579790B (en) | Method for evaluating electrostatic effect of self-powered neutron detector | |
| Liu et al. | Research on the simulation method for HVDC continuous positive corona discharge | |
| Khaled et al. | Experimental and analytical study on the performance of novel design of efficient two‐stage electrostatic precipitator | |
| WO2024099010A1 (en) | Method and system for assessing risk of internal charging of dielectrics of spacecraft on synchronous orbit, and terminal | |
| CN113933610A (en) | Cable insulation medium space charge dynamic distribution calculation method | |
| Bian et al. | Comparative study on corona-generated audible noise and radio noise of ac long-term operating conductors with two bundle types | |
| Pan et al. | Investigation of cavity PD physical processes at DC voltage by simulation | |
| Bian et al. | Negative DC corona discharge effects on the ultraviolet photon count in rod to plane air gaps | |
| Li et al. | Dividing and setting strategy of improving calculation efficiency for needle electrode corona discharge with a large‐scale space | |
| Juda et al. | 3D-simulation of ionisation gauges and comparison with measurements | |
| Du et al. | Computation of total electric field considering natural wind under high-altitude UHVDC transmission lines | |
| Hruby et al. | Three-dimensional Hybrid computer modeling of langmuir probes of finite dimensions in medium pressure plasmas | |
| He et al. | Charging models for airborne suspended particles around HVDC lines | |
| Babich et al. | Initiation of positive streamers near uncharged ice hydrometeors in the thundercloud field |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |