CN105572718A - Flight space atmosphere neutron fluence rate calculating method and system thereof - Google Patents
Flight space atmosphere neutron fluence rate calculating method and system thereof Download PDFInfo
- Publication number
- CN105572718A CN105572718A CN201410550332.7A CN201410550332A CN105572718A CN 105572718 A CN105572718 A CN 105572718A CN 201410550332 A CN201410550332 A CN 201410550332A CN 105572718 A CN105572718 A CN 105572718A
- Authority
- CN
- China
- Prior art keywords
- fluence rate
- target location
- atmospheric neutron
- neutron fluence
- longitude
- 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.)
- Granted
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
The invention relates to a flight space atmosphere neutron fluence rate calculating method and a system thereof. The method comprises the following steps of acquiring a cut-off rigidity; acquiring a first multiple corresponding to the cut-off rigidity; acquiring a second multiple corresponding to a target height of a target position and a third multiple corresponding to a standard height; according to the first multiple, the second multiple and the third multiple, calculating a fourth multiple; according to a standard fluence rate of a standard position in a preset energy scope and the fourth multiple, calculating a first fluence rate to be converted; acquiring an energy threshold value and a fifth multiple corresponding to the energy threshold value; according to the first fluence rate to be converted and the fifth multiple, calculating an atmosphere neutron fluence rate of target equipment at the target position. Through the technical scheme of the invention, the atmosphere neutron fluence rate can be calculated in a height scope including an area near ground; a longitude influence is considered; and calculating accuracy of the atmosphere neutron fluence rate is increased so that accurate and wide assessment can be performed on an atmosphere-neutron single particle effect.
Description
Technical field
The present invention relates to neutron fluence rate technical field, in particular to a kind of flying area atmospheric neutron fluence rate computing method and a kind of flying area atmospheric neutron fluence rate computing system.
Background technology
Air environment with the complicated microelectronic component of storage organization will inevitably meet with the high energy atmospheric neutron of about every square centimeter of 300 ~ 18000 1MeV ~ 1000MeV per hour in the place environment of flying height (3000 ~ 20000 meters).These high-energy neutrons can penetrate cabin skin, beat on the kernel instruction control module or critical data storage unit of air environment, produce soft error and hard fault, cause navigating (navigation neceiver), radar-probing system (Active Phased Array Radar), data network (the AFDX network switch), communication (optical fiber/bus), high speed computer system, avionic device, engine (FADEC), telex system, automatic Pilot technology, flight alarm, display screen, there is blank screen in other flight system etc. containing electron device, crash, reset, restart, loss of data, the safety hazards such as order loss.
Up to now, domestic to air environment meet with in place environment high-energy neutron do not carry out specific aim protection with evaluate.In order to set up protection and appraisement system, quantitatively characterizing must be carried out to the harm of atmospheric neutron single particle effect.And atmospheric neutron fluence rate is the important parameter of harm quantitatively characterizing.
Mainly contain two methods in the world and calculate atmospheric neutron fluence rate, the broadcast model of Ge Shi Boeing exploitation, another is NASA-Langley model, is called for short NASA model.
Based on the value actually of the 1-10MeV atmospheric neutron fluence rate carried out the sixties in the world, Boeing develops initial simplification Boeing model.This model hypothesis 1-10MeV atmospheric neutron fluence rate can be decomposed into three factors, and a factor changes along with the change of height, and a factor changes along with the change of latitude, and a factor considers the energy of neutron.The people such as A.Taber, by carrying out matching to a large amount of detection data, have shown that 1-10MeVBoeing empirical model is used for predicting the atmospheric neutron flux of differing heights and latitude.Its concrete formula is as follows:
dN/dE＝0.3459E
-0.9219×exp[-0.01522·(lnE)
2]
N(E)dE＝26E
-1.16±0.2×exp[-(0.0069x)]·dE
φ
1-10(ω
Lat)＝0.6252exp{-0.461[cos(2×ω
Lat)]
2-0.94cos(2×ω
Lat)+0.252}
In formula: N, φ are neutron flux;
E is neutron energy;
X is atmospheric depth, g/cm2;
ω Lat is latitude.
U.S. NASA is also by having researched and proposed NASA-Langley model to atmospheric neutron distribution situation.NASA model is called as AIR model, is to utilize the exploitation of the flight survey data of the sixties to the seventies to form.In NASA model, atmospheric neutron flux affects by three major parameters, is respectively: atmospheric density (g/cm
2), vertical cut rigidity (GV) and sun environmental baseline.
NASA model is prediction 1-10MeV atmospheric neutron flux empirical model more accurately, and this model basic calculating formula is as follows:
Wherein:
F(R,C)＝(Λ/λ)·f(R,C)·exp(x
m/Λ-x
m/λ)
f(R,C)＝exp(250/λ)φ(250,R,C)
λ＝165+2R
x
m＝50+ln{2000+exp[-2(C-100)]}
In formula, the atmospheric neutron flux that φ 1-10 (x, R, C) is 1-10MeV;
X is atmospheric density, g/cm
2; R is cut-off rigidity;
C is solar activity constant;
λ, xm and Λ are the intermediate parameters calculated.
Wherein, air is thick relevant to height A, and A unit is feet:
The factor wherein affecting atmospheric neutron fluence rate mainly contains 4: atmospheric density, cut-off rigidity, energy range, solar activity.Wherein atmospheric density is by highly determining, cut-off rigidity is then determined by longitude and latitude.
Boeing model of the prior art mainly considers height and the impact of latitude on atmospheric neutron fluence rate, have ignored the impact of longitude, therefore applies the atmospheric neutron fluence rate that Boeing model obtains accurate not.Compared with Boeing model, although NASA model is comparatively accurate.But NASA model is only applicable to lower than 20km height.And due to the impact of parameter lambda=165+2R, require that atmospheric density is lower than 250g/cm
2, inapplicable and floor level.
Summary of the invention
Technical matters to be solved by this invention is, how to design one and is applicable to various height, and the algorithm of the high atmospheric neutron fluence rate of counting accuracy or system, to carrying out extensively atmospheric neutron single particle effect and assessing accurately.
For this purpose, the present invention proposes a kind of flying area atmospheric neutron fluence rate computing method, comprising: inquiry longitude, latitude is corresponding with cut-off rigidity shows, the cut-off rigidity that the longitude of acquisition and target location, latitude are corresponding; The multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtains first multiple corresponding with described cut-off rigidity; The multiple proportion table of inquiry height, latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtains second multiple corresponding with the object height of described target location, obtains the triple corresponding with described calibrated altitude; According to described first multiple, the second multiple and triple, calculate the quadruple of described target location atmospheric neutron fluence rate relative to described normal place atmospheric neutron fluence rate in the horizontal plane; According to the standard fluence rate of described normal place within the scope of preset energy and described quadruple, calculate described target location and first treat conversion fluence rate within the scope of described preset energy; Acquisition can cause the energy threshold of the target device generation single particle effect being positioned at described target location, the multiple proportion table of the fluence rate within the scope of the inquiry fluence rate of energy threshold and described preset energy, obtains five multiple corresponding with described energy threshold; The atmospheric neutron fluence rate of described target device when described target location is calculated according to described first fluence rate to be converted and described 5th multiple.
Preferably, the step calculating described quadruple specifically comprises: the first ratio calculating described second multiple and described triple, calculates the product of described first multiple and described first ratio as described quadruple.
Preferably, inquiry described longitude, latitude is corresponding with cut-off rigidity show before also comprise: the longitude judging whether to get described target location, if the longitude of described target location can be got, then inquire about described longitude, latitude is corresponding with cut-off rigidity shows, obtain the cut-off rigidity corresponding with the longitude of target location, latitude, if the longitude of described target location can not be got, then inquire about described height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of described target location and target latitude, according to described standard fluence rate and described 6th multiple, calculate described target location within the scope of described preset energy second treats conversion fluence rate, obtain described 5th multiple, the atmospheric neutron fluence rate of described target device when described target location is calculated according to described second fluence rate to be converted and described 5th multiple.
Preferably, also comprise: assess according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
Preferably, described preset energy scope is that 10MeV is to just infinite.
The invention allows for a kind of flying area atmospheric neutron fluence rate computing system, comprise: query unit, for inquiring about longitude, latitude is corresponding with cut-off rigidity to be shown, obtain the longitude with target location, the cut-off rigidity that latitude is corresponding, the multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtain first multiple corresponding with described cut-off rigidity, inquiry height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain second multiple corresponding with the object height of described target location, obtain the triple corresponding with described calibrated altitude, acquisition can cause the energy threshold of the target device generation single particle effect being positioned at described target location, the multiple proportion table of the fluence rate within the scope of the inquiry fluence rate of energy threshold and described preset energy, obtain five multiple corresponding with described energy threshold, computing unit, for according to described first multiple, the second multiple and triple, calculate the quadruple of described target location atmospheric neutron fluence rate relative to described normal place atmospheric neutron fluence rate in the horizontal plane, according to the standard fluence rate of described normal place within the scope of preset energy and described quadruple, calculate described target location and first treat conversion fluence rate within the scope of described preset energy, calculate the atmospheric neutron fluence rate of described target device when described target location according to described first fluence rate to be converted and described 5th multiple.
Preferably, described computing unit, for calculating the first ratio of described second multiple and described triple, calculates the product of described first multiple and described first ratio as described quadruple.
Preferably, also comprise: judging unit, for judging whether the longitude that can get described target location, wherein, described query unit judges the Longitude In System that can get described target location at described judging unit, inquire about described longitude, latitude is corresponding with cut-off rigidity to be shown, obtain the longitude with target location, the cut-off rigidity that latitude is corresponding, the Longitude In System that can not get described target location is judged at described judging unit, inquire about described height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of described target location and target latitude, and described computing unit is according to described standard fluence rate and described 6th multiple, calculate described target location within the scope of described preset energy second treats conversion fluence rate, obtain described 5th multiple, the atmospheric neutron fluence rate of described target device when described target location is calculated according to described second fluence rate to be converted and described 5th multiple.
Preferably, also comprise: assessment unit, for assessing according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
Preferably, described preset energy scope is that 10MeV is to just infinite.。
Pass through technique scheme, atmospheric neutron fluence rate can be calculated in the altitude range comprising Near Ground, and consider the impact of longitude, improve the counting accuracy of atmospheric neutron fluence rate, to carry out more accurately for atmospheric neutron single particle effect and to assess widely.
Accompanying drawing explanation
Can understanding the features and advantages of the present invention clearly by reference to accompanying drawing, accompanying drawing is schematic and should not be construed as and carry out any restriction to the present invention, in the accompanying drawings:
Fig. 1 shows the schematic flow diagram of flying area atmospheric neutron fluence rate computing method according to an embodiment of the invention;
Fig. 2 shows the schematic flow diagram of flying area atmospheric neutron fluence rate computing method in accordance with another embodiment of the present invention;
Fig. 3 shows the schematic block diagram of flying area atmospheric neutron fluence rate computing system according to an embodiment of the invention.
Embodiment
Can more clearly understand above-mentioned purpose of the present invention, feature and advantage, below in conjunction with the drawings and specific embodiments, the present invention is further described in detail.It should be noted that, when not conflicting, the feature in the embodiment of the application and embodiment can combine mutually.
Set forth a lot of detail in the following description so that fully understand the present invention; but; the present invention can also adopt other to be different from other modes described here and implement, and therefore, protection scope of the present invention is not by the restriction of following public specific embodiment.
As shown in Figure 1, flying area atmospheric neutron fluence rate computing method, comprising: S1 according to an embodiment of the invention, and inquiry longitude, latitude is corresponding with cut-off rigidity shows, the cut-off rigidity that the longitude of acquisition and target location, latitude are corresponding; Wherein, longitude, latitude be corresponding with cut-off rigidity to be shown as shown in table 1, and wherein north latitude is just, south latitude is negative, and longitude is represented by east longitude 0-360 °, is 360-x ° of W.As 30S, 120W, be-30N, 240E.
Table 1
Particularly, such as target location be in north latitude 75 °, east longitude 30 °, be highly 6100 meters, so can inquire north latitude 75 ° according to table 1, the cut-off rigidity R=0.15 that east longitude 30 ° is corresponding.
S2, the multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtains first multiple corresponding with cut-off rigidity; Wherein under cut-off rigidity and calibrated altitude, the multiple proportion table of atmospheric neutron fluence rate is as shown in table 2,
| Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple |
| 0.00 | 561.70 | 2.20 | 506.89 | 8.38 | 193.2 | 13.43 | 111.25 |
| 0.01 | 561.70 | 2.30 | 503.85 | 8.99 | 174.36 | 13.50 | 110.34 |
| 0.02 | 561.70 | 2.32 | 497.65 | 9.23 | 167.52 | 13.57 | 109.77 |
| 0.03 | 561.70 | 2.40 | 493.72 | 9.56 | 163.33 | 13.62 | 109.49 |
| 0.05 | 561.70 | 2.56 | 483.35 | 9.57 | 162.04 | 13.76 | 108.20 |
| 0.06 | 561.70 | 2.73 | 463.07 | 9.66 | 161.49 | 13.77 | 108.10 |
| 0.08 | 561.70 | 2.75 | 463.07 | 10.17 | 153.78 | 13.80 | 107.74 |
| 0.09 | 561.70 | 2.85 | 454.95 | 10.55 | 146.52 | 13.88 | 107.01 |
| 0.13 | 561.70 | 3.01 | 447.68 | 10.81 | 142.74 | 13.92 | 106.75 |
| 0.14 | 561.70 | 3.05 | 441.25 | 10.97 | 140.69 | 14.06 | 105.50 |
| 0.17 | 561.69 | 3.20 | 429.35 | 10.98 | 140.12 | 14.10 | 105.15 |
| 0.19 | 561.70 | 3.30 | 419.22 | 11.10 | 138.42 | 14.19 | 104.37 |
| 0.20 | 561.70 | 3.56 | 399.60 | 11.23 | 132.08 | 14.20 | 104.20 |
| 0.24 | 561.70 | 3.68 | 385.15 | 11.25 | 139.55 | 14.35 | 103.09 |
| 0.25 | 561.70 | 3.71 | 387.28 | 11.32 | 134.72 | 14.37 | 102.84 |
| 0.27 | 561.70 | 4.18 | 356.77 | 11.59 | 131.82 | 14.46 | 102.01 |
| 0.28 | 561.70 | 4.28 | 346.05 | 11.61 | 131.57 | 14.64 | 100.62 |
| 0.36 | 561.70 | 4.30 | 360.00 | 11.66 | 130.16 | 14.65 | 100.54 |
| 0.38 | 561.70 | 4.33 | 347.28 | 11.67 | 130.16 | 14.71 | 100.05 |
| 0.42 | 561.70 | 4.55 | 333.40 | 11.72 | 131.18 | 14.94 | 98.24 |
| 0.52 | 561.70 | 4.71 | 329.90 | 12.22 | 123.64 | 14.95 | 98.16 |
| 0.55 | 561.70 | 4.90 | 313.20 | 12.35 | 122.38 | 15.16 | 96.56 |
| 0.68 | 561.69 | 5.02 | 311.59 | 12.36 | 122.50 | 15.37 | 95.00 |
| 0.74 | 561.67 | 5.37 | 290.39 | 12.39 | 121.59 | 16.00 | 90.21 |
| 0.78 | 561.63 | 5.43 | 287.54 | 12.67 | 119.06 | 17.00 | 82.41 |
| 0.79 | 561.67 | 5.58 | 280.60 | 12.72 | 119.06 | ||
| 0.85 | 561.50 | 5.83 | 272.62 | 12.73 | 118.84 | ||
| 0.89 | 561.48 | 6.04 | 260.48 | 12.87 | 116.93 | ||
| 0.93 | 561.39 | 6.10 | 262.51 | 12.88 | 116.82 | ||
| 0.99 | 560.73 | 6.63 | 239.35 | 12.91 | 116.30 | ||
| 1.08 | 560.16 | 6.78 | 233.48 | 12.99 | 115.68 |
| Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple | Cut-off rigidity | Multiple |
| 1.10 | 559.93 | 6.86 | 232.13 | 13.02 | 115.37 | ||
| 1.15 | 560.04 | 7.13 | 222.10 | 13.06 | 114.96 | ||
| 1.19 | 558.02 | 7.17 | 230.46 | 13.07 | 114.66 | ||
| 1.44 | 551.49 | 7.41 | 213.65 | 13.20 | 113.35 | ||
| 1.69 | 540.63 | 7.57 | 211.90 | 13.22 | 113.35 | ||
| 1.85 | 530.70 | 7.73 | 204.62 | 13.26 | 112.96 | ||
| 1.86 | 534.99 | 7.89 | 196.71 | 13.31 | 112.66 | ||
| 2.19 | 514.32 | 8.03 | 195.44 | 13.38 | 111.74 |
Table 2
Particularly, calibrated altitude can be set to 40000 feet, approximates 12190 meters, and so when cut-off rigidity R=0.15, the first multiple is 561.7.
S3, the multiple proportion table of inquiry height, latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtains second multiple corresponding with the object height of target location, obtains the triple corresponding with calibrated altitude; Wherein, highly, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate is in the horizontal plane as shown in table 3,
Table 3
Because latitude and longitude factor were considered in step s 2, so only need inquire about and highly corresponding multiple in table 3, second multiple corresponding with object height 6100 meters is 0.13, and the triple corresponding with calibrated altitude 12190 meters is 0.88.
S4, according to the first multiple, the second multiple and triple, calculated target positions atmospheric neutron fluence rate is relative to the quadruple of normal place atmospheric neutron fluence rate in the horizontal plane;
Particularly, the ratio that quadruple equals the second multiple and triple is multiplied by the first multiple, and namely 561.7 × 0.13/088, approximate 83.
S5, according to the standard fluence rate of normal place within the scope of preset energy and quadruple, calculated target positions first treats conversion fluence rate within the scope of preset energy;
It is just infinite that preset energy scope can be that 10MeV arrives, and so first treats that conversion fluence rate equals standard fluence rate and quadruple is long-pending, namely 14.4 × 83, approximate 1195/cm
2h.
S6, obtains and can cause the energy threshold of the target device generation single particle effect being positioned at target location, and the multiple proportion table of the fluence rate within the scope of the fluence rate of inquiry energy threshold and preset energy, obtains five multiple corresponding with energy threshold; The multiple proportion table of the fluence rate within the scope of the fluence rate of energy threshold and preset energy is as shown in table 4,
| Sequence number | Energy (MeV) | Neutron fluence rate (#/cm2h) | With > 10MeV multiple proportion |
| 1 | >1 | 9200 | 1.53 |
| 2 | >2 | 8300 | 1.38 |
| 3 | >3 | 7700 | 1.28 |
| 4 | >4 | 7300 | 1.22 |
| 5 | >5 | 7000 | 1.17 |
| 6 | >6 | 6800 | 1.13 |
| 7 | >7 | 6500 | 1.08 |
| 8 | >8 | 6400 | 1.07 |
| 9 | >9 | 6200 | 1.03 |
| 10 | >10 | 6000 | 1.00 |
| 11 | 1-10 | 3200 | 0.50 |
Table 4
The energy threshold such as got is 2MeV, and the 5th so corresponding with energy threshold multiple is 1.38, certainly, also can directly inquire about the neutron fluence rate corresponding with energy threshold, and being divided by with the neutron fluence rate in preset range obtains the 5th multiple.
S7, calculates the atmospheric neutron fluence rate of target device when target location according to the first fluence rate to be converted and the 5th multiple.
Particularly, atmospheric neutron fluence rate equals the product that first treats conversion fluence rate and the 5th multiple, and namely 1195 × 1.38, approximate 1653/cm
2h, thus draw, energy threshold be the target device of 2MeV north latitude 75 °, east longitude 30 °, during the target location of height 6100 meters, its atmospheric neutron fluence rate is 1653/cm
2h.And then can assess according to obtaining the NEUTRON PROTECTION of atmospheric neutron fluence rate to target device.
Preferably, step S4 comprises: the first ratio calculating the second multiple and triple, calculates the product of the first multiple and the first ratio as quadruple.
As shown in Figure 2, preferably, also comprised before step S1: the longitude judging whether to get target location, if can get the longitude of target location, then enters step S1; If the longitude of target location can not be got, then inquire about the multiple proportion table of height, latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of target location and target latitude.
Due in certain type situation down to height and the dimension data that can get target location, so then directly can inquire about the 6th multiple corresponding to this height and dimension from table 3.The dimension of such as target location is 76 °, is highly 6100 meters, and so can inquire corresponding 6th multiple according to table 3 is 90.5.
According to standard fluence rate and the 6th multiple, within the scope of preset energy second of calculated target positions treats conversion fluence rate.Standard fluence rate is constant, is still 14.4/cm
2h, so second treats that conversion fluence rate equals standard fluence rate and the 6th multiple is long-pending, namely 14.4 × 90.5, approximate 1303/cm
2h.
When the energy threshold of target device is 2MeV, from table 4, obtain the 5th multiple, calculate the atmospheric neutron fluence rate of target device when target location according to the second fluence rate to be converted and the 5th multiple.Namely 1303 × 1.38, approximate 1803/cm
2h.Namely energy threshold be the target device of 2MeV north latitude 75 °, during the height target location of 6100 meters, its atmospheric neutron fluence rate is 1803/cm
2h.
Above-mentioned normal place can elect New York as.
Preferably, also comprise: assess according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
Preferably, preset energy scope is that 10MeV is to just infinite.
As shown in Figure 3, flying area atmospheric neutron fluence rate computing system 20 comprises according to an embodiment of the invention: query unit 21, for inquiring about longitude, latitude is corresponding with cut-off rigidity to be shown, obtain the longitude with target location, the cut-off rigidity that latitude is corresponding, the multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtain first multiple corresponding with cut-off rigidity, inquiry height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain second multiple corresponding with the object height of target location, obtain the triple corresponding with calibrated altitude, acquisition can cause the energy threshold of the target device generation single particle effect being positioned at target location, the multiple proportion table of the fluence rate within the scope of the inquiry fluence rate of energy threshold and preset energy, obtain five multiple corresponding with energy threshold, computing unit 22, for according to the first multiple, the second multiple and triple, calculated target positions atmospheric neutron fluence rate is relative to the quadruple of normal place atmospheric neutron fluence rate in the horizontal plane, according to the standard fluence rate of normal place within the scope of preset energy and quadruple, calculated target positions first treats conversion fluence rate within the scope of preset energy, calculates the atmospheric neutron fluence rate of target device when target location according to the first fluence rate to be converted and the 5th multiple.
Preferably, computing unit 22, for calculating the first ratio of the second multiple and triple, calculates the product of the first multiple and the first ratio as quadruple.
Preferably, also comprise: judging unit 23, for judging whether the longitude that can get target location, wherein, query unit 21 judges the Longitude In System that can get target location at judging unit 23, inquiry longitude, latitude is corresponding with cut-off rigidity to be shown, obtain the longitude with target location, the cut-off rigidity that latitude is corresponding, the Longitude In System that can not get target location is judged at judging unit 23, inquiry height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of target location and target latitude, and computing unit 22 is according to standard fluence rate and the 6th multiple, within the scope of preset energy second of calculated target positions treats conversion fluence rate, obtain the 5th multiple, the atmospheric neutron fluence rate of target device when target location is calculated according to the second fluence rate to be converted and the 5th multiple.
Preferably, also comprise: assessment unit 24, for assessing according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
Preferably, preset energy scope is that 10MeV is to just infinite.
According to the embodiment of the present invention, additionally providing a kind of non-volatile machine readable media, storing the program product for calculating atmospheric neutron fluence rate, program product comprises above-mentioned two kinds of program products.
According to the embodiment of the present invention, additionally provide a kind of machine readable program, program makes machine perform atmospheric neutron fluence rate computing method arbitrary in as above technical scheme.
According to the embodiment of the present invention, additionally provide a kind of storage medium storing machine readable program, wherein, machine readable program makes machine perform atmospheric neutron fluence rate computing method arbitrary in as above technical scheme.
More than be described with reference to the accompanying drawings technical scheme of the present invention, consider in correlation technique, for the calculating of atmospheric neutron fluence rate, what have does not consider longitude factor, thus cause accuracy in computation lower, assessment accuracy for NEUTRON PROTECTION is reduced, and the calculating being difficult to carry out atmospheric neutron fluence rate in Near Ground altitude range had, makes the scope of assessment for NEUTRON PROTECTION limited.By the technical scheme of the application, atmospheric neutron fluence rate can be calculated in the altitude range comprising Near Ground, and consider the impact of longitude, improve the counting accuracy of atmospheric neutron fluence rate, to carry out more accurately for atmospheric neutron single particle effect and to assess widely.
In the present invention, term " first ", " second ", " the 3rd ", " the 4th ", " the 5th ", " the 6th " only for describing object, and can not be interpreted as instruction or hint relative importance.Term " multiple " refers to two or more, unless otherwise clear and definite restriction.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. flying area atmospheric neutron fluence rate computing method, is characterized in that, comprising:
Inquiry longitude, latitude is corresponding with cut-off rigidity shows, the cut-off rigidity that the longitude of acquisition and target location, latitude are corresponding;
The multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtains first multiple corresponding with described cut-off rigidity;
The multiple proportion table of inquiry height, latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtains second multiple corresponding with the object height of described target location, obtains the triple corresponding with described calibrated altitude;
According to described first multiple, the second multiple and triple, calculate the quadruple of described target location atmospheric neutron fluence rate relative to described normal place atmospheric neutron fluence rate in the horizontal plane;
According to the standard fluence rate of described normal place within the scope of preset energy and described quadruple, calculate described target location and first treat conversion fluence rate within the scope of described preset energy;
Acquisition can cause the energy threshold of the target device generation single particle effect being positioned at described target location, the multiple proportion table of the fluence rate within the scope of the inquiry fluence rate of energy threshold and described preset energy, obtains five multiple corresponding with described energy threshold;
The atmospheric neutron fluence rate of described target device when described target location is calculated according to described first fluence rate to be converted and described 5th multiple.
2. flying area atmospheric neutron fluence rate computing method according to claim 1, it is characterized in that, the step calculating described quadruple specifically comprises:
Calculate the first ratio of described second multiple and described triple, calculate the product of described first multiple and described first ratio as described quadruple.
3. flying area atmospheric neutron fluence rate computing method according to claim 1, is characterized in that, the described longitude of inquiry, latitude is corresponding with cut-off rigidity show before also comprise:
Judge whether the longitude that can get described target location, if the longitude of described target location can be got, then inquire about described longitude, latitude is corresponding with cut-off rigidity shows, obtain the cut-off rigidity corresponding with the longitude of target location, latitude;
If the longitude of described target location can not be got, then inquire about described height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of described target location and target latitude, according to described standard fluence rate and described 6th multiple, calculate described target location within the scope of described preset energy second treats conversion fluence rate, obtain described 5th multiple, the atmospheric neutron fluence rate of described target device when described target location is calculated according to described second fluence rate to be converted and described 5th multiple.
4. flying area atmospheric neutron fluence rate computing method according to any one of claims 1 to 3, is characterized in that, also comprise:
Assess according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
5. flying area atmospheric neutron fluence rate computing method according to any one of claims 1 to 3, is characterized in that, described preset energy scope is that 10MeV is to just infinite.
6. a flying area atmospheric neutron fluence rate computing system, is characterized in that, comprising:
Query unit, for inquiring about longitude, latitude is corresponding with cut-off rigidity shows, and obtains the cut-off rigidity corresponding with the longitude of target location, latitude,
The multiple proportion table of atmospheric neutron fluence rate under inquiry cut-off rigidity and calibrated altitude, obtains first multiple corresponding with described cut-off rigidity,
The multiple proportion table of inquiry height, latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtains second multiple corresponding with the object height of described target location, obtains the triple corresponding with described calibrated altitude,
Acquisition can cause the energy threshold of the target device generation single particle effect being positioned at described target location, the multiple proportion table of the fluence rate within the scope of the inquiry fluence rate of energy threshold and described preset energy, obtains five multiple corresponding with described energy threshold;
Computing unit, for according to described first multiple, the second multiple and triple, calculates the quadruple of described target location atmospheric neutron fluence rate relative to described normal place atmospheric neutron fluence rate in the horizontal plane,
According to the standard fluence rate of described normal place within the scope of preset energy and described quadruple, calculate described target location and first treat conversion fluence rate within the scope of described preset energy,
The atmospheric neutron fluence rate of described target device when described target location is calculated according to described first fluence rate to be converted and described 5th multiple.
7. flying area atmospheric neutron fluence rate computing system according to claim 5, it is characterized in that, described computing unit, for calculating the first ratio of described second multiple and described triple, calculates the product of described first multiple and described first ratio as described quadruple.
8. flying area atmospheric neutron fluence rate computing system according to claim 5, is characterized in that, also comprise:
Judging unit, for judging whether the longitude that can get described target location,
Wherein, described query unit judges the Longitude In System that can get described target location at described judging unit, inquire about described longitude, latitude is corresponding with cut-off rigidity to be shown, obtain the longitude with target location, the cut-off rigidity that latitude is corresponding, the Longitude In System that can not get described target location is judged at described judging unit, inquire about described height, the multiple proportion table of latitude and normal place atmospheric neutron fluence rate in the horizontal plane, obtain six multiple corresponding with the object height of described target location and target latitude, and described computing unit is according to described standard fluence rate and described 6th multiple, calculate described target location within the scope of described preset energy second treats conversion fluence rate, obtain described 5th multiple, the atmospheric neutron fluence rate of described target device when described target location is calculated according to described second fluence rate to be converted and described 5th multiple.
9. flying area atmospheric neutron fluence rate computing system according to any one of claim 5 to 8, is characterized in that, also comprise:
Assessment unit, for assessing according to the NEUTRON PROTECTION of atmospheric neutron fluence rate to described target device.
10. flying area atmospheric neutron fluence rate computing system according to any one of claim 5 to 8, is characterized in that, described preset energy scope is that 10MeV is to just infinite.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410550332.7A CN105572718B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410550332.7A CN105572718B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate calculation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105572718A true CN105572718A (en) | 2016-05-11 |
| CN105572718B CN105572718B (en) | 2019-02-12 |
Family
ID=55883059
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410550332.7A Active CN105572718B (en) | 2014-10-16 | 2014-10-16 | Flying area atmospheric neutron fluence rate calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572718B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014136990A1 (en) * | 2013-03-08 | 2014-09-12 | 国立大学法人九州大学 | High-sensitivity neutron detection method using self-radioactivation of scintillator |
-
2014
- 2014-10-16 CN CN201410550332.7A patent/CN105572718B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014136990A1 (en) * | 2013-03-08 | 2014-09-12 | 国立大学法人九州大学 | High-sensitivity neutron detection method using self-radioactivation of scintillator |
Non-Patent Citations (4)
| Title |
|---|
| E. NORMAND ET AL: "Altitude and Latitude Variations in Avionics SEU and Atmospheric Neutron Flux", 《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》 * |
| G. HUBERT ET AL: "A Neutron Spectrometer for Avionic Environment Investigations", 《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》 * |
| MUNEHIKO KOWATARI ET AL: "Evaluation of Geomagnetic Latitude Dependence of the Cosmic-ray Induced Environmental Neutrons in Japan", 《JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY》 * |
| 张紫霞等: "空间环境中子效应及测量技术", 《装备环境工程》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105572718B (en) | 2019-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| De Winter et al. | Wind extremes in the North Sea Basin under climate change: An ensemble study of 12 CMIP5 GCMs | |
| Dubyagin et al. | Can flow bursts penetrate into the inner magnetosphere? | |
| CN104007479B (en) | A kind of Ionospheric Tomography based on multiple dimensioned subdivision and Ionospheric delay correcting method | |
| WO2008117787A1 (en) | Map matching system, map matching method and program | |
| Hsu et al. | Assessment of the impact of FORMOSAT‐7/COSMIC‐2 GNSS RO observations on midlatitude and low‐latitude ionosphere specification: Observing system simulation experiments using ensemble square root filter | |
| CN105974905A (en) | Simulated test system and method for single-particle overturn fault of aviation data bus | |
| CN107315182A (en) | Passive location excitation equivalence aptitude tests method based on semi-physical simulation | |
| Ratnam et al. | Performance evaluation of selected ionospheric delay models during geomagnetic storm conditions in low-latitude region | |
| Brecht et al. | The role of the Martian crustal magnetic fields in controlling ionospheric loss | |
| Lin et al. | Atmospheric gravity waves in the ionosphere and thermosphere during the 2017 solar eclipse | |
| Hisadomi et al. | Multiple gamma‐ray glows and a downward TGF observed from nearby thunderclouds | |
| Ribstein et al. | The interaction between gravity waves and solar tides: Results from 4‐D ray tracing coupled to a linear tidal model | |
| Settimi et al. | The IONORT-ISP-WC system: inclusion of an electron collision frequency model for the D-layer | |
| Guo et al. | Propagation of Large‐Scale Solar Wind Events in the Outer Heliosphere From a Numerical MHD Simulation | |
| CN104933316A (en) | Ionized layer obscuration retrieval method based on two-parameter mixture regularization | |
| CN105572717B (en) | Avionic device atmospheric neutron fluence rate calculation method and system | |
| CN105572718A (en) | Flight space atmosphere neutron fluence rate calculating method and system thereof | |
| Lora-Clavijo et al. | PBH mass growth through radial accretion during the radiation dominated era | |
| Wang et al. | Eddy‐induced acceleration of Argo floats | |
| CN105718714A (en) | Method and system for determining atmospheric neutron single-particle upset rate of microcircuit | |
| CN103675845B (en) | A kind of single carrier multiple antennas occultation signal generating system | |
| CN103837130B (en) | For data processing method and the device of airborne lidar system | |
| CN105676016B (en) | The method and device in neutron single-particle effect device sensitivity section is obtained using BGR | |
| CN103605100A (en) | Positioning error simulation method for lightning detection system | |
| CN109470947A (en) | A method of atmospheric neutron single particle effect section is determined using test data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |