Abstract
In the complex situation in which we live, preventive measures are increasingly being taken against the threat of nuclear attack and ionizing radiation, nuclear shelters are being built for private and public purposes. The purpose of the authors is to provide a practical approach to quickly calculate the thickness of the concrete shelter wall for such cases. The shelter wall thickness is assessed on the basis of published data on radiation dose calculations, corroborated by field data from nuclear tests, and information or calculations of the dose attenuation in the shelter wall. Doses from nuclear tests or simulations are summarized in tables in relation to bomb yield and distance from the explosion point. A dose criterion is chosen as the acceptable dose behind the shield. Neutron dose transmission factors, defined as ratio of the dose criterion to the dose without shielding are calculated. The concrete wall thickness corresponding to the neutron dose transmission factor was taken from reference sources. For gamma rays, the shield thickness is calculated based on an analytical relationship between dose without shield, the dose criterion, and the shield thickness. Data tables are provided with the estimated concrete wall thickness in relation to bomb yield and distance.
1 Introduction
The authors' first goal was to find a reliable source for a selection of doses caused by fission neutrons, fission gamma rays, and secondary gamma rays from nuclear fission bombs. Such information was found in Ref. [1]. The second goal was to estimate the dose transmission factor (DTF) for neutrons as a ratio between the dose behind the shield and the dose without shield and to find a source where information is given on the relationship between DTF and the thickness of the shelter wall. Information was found in Ref. [2]. The third goal was to propose a formula to estimate the shelter wall thickness against gamma rays, based on the relation between the dose behind the shield, dose without shield and wall thickness. Appropriate commonly used materials such as different concrete types are proposed for both neutron and gamma ray protection. Fallout gamma radiation is taken into account to estimate the roof wall thickness. A comparison is made with the results of applying different approaches to estimate the thickness of protective walls. The most protective place for a nuclear shelter is below the ground level. In the case of underground location of the shelter, the soil shield should be taken into account in the protection of the shelter construction.
1.1 Diagram of the Shelter Wall Thickness Evaluation Process
The recommended shelter wall thickness is the highest estimated thickness value.
2 Nuclear Explosion Above Ground Level. Above Ground Shelter
2.1 Shielding Against Neutrons.
The energies of the neutrons received at some distance from nuclear explosion cover a very wide range, from several millions down to a fraction of an eV (1 eV = 1.602 × 10−19J). With spectra for the finite number of energy groups results of neutron dose calculations with a neutron transport computer codes corroborated by test data, are given in Fig. 1 developed in accordance with Fig. 8.123a from Ref. [1].
2.1.1 Determination of the Initial Neutron Dose.
From Fig. 1, we read the initial neutron dose per kiloton total yield as a function of slant range from fission weapon air bursts at elevation about of 90 m, based on 0.9 normal sea-level air density, based on Ref. [1].
The doses in the figure are valid for air-burst above 90 m, 0.9 normal sea-level air density [1]. In cases of contact surface burst, the doses from Fig. 1 shall be multiplied by factor of 0.5 [1]. For bursts between the surface and 90 m an approximate value of a dose can be calculated by linear interpolation between a surface burst and one at 90 m or above [1].
2.1.2 Presentation of the Expected Neutron Doses From a Fission Bomb as a Function of Bomb Yield and Distance From the Nuclear Explosion Point.
The effective dose can be calculated for fission neutrons also based on neutron fluence to effective dose conversion coefficients. Neutron fluence per kiloton energy yield can be found in Ref. [1] in Fig. 8.117a. Values of neutron fluence to effective dose conversion coefficients depending on neutron energy can be found in Ref. [3] in Fig. 5 of the cited document. The energy spectrum of neutrons leaking from a weapon is primarily a fission spectrum [4] with average energy of about 2 MeV [5]. For 2 MeV neutrons, we obtain similar doses to those shown in Table 1, but the data in Table 1 are based on more direct method.
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 30 | 60 | 90 | 120 | 150 | 300 | 1500 | 3000 |
1500 | 1371.6 | 1.5 | 3 | 4.5 | 6 | 7.5 | 15 | 75 | 150 |
2000 | 1828.8 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 1 | 5 | 10 |
2500 | 2286 | 0.01 | 0.02 | 0.03 | 0.04 | 0.05 | 0.1 | 0.5 | 1 |
3000 | 2743.2 | 0.001 | 0.002 | 0.003 | 0.004 | 0.005 | 0.01 | 0.05 | 0.1 |
3500 | 3200.4 | 0.00009 | 0.00018 | 0.00027 | 0.00036 | 0.00045 | 0.0009 | 0.0045 | 0.009 |
4000 | 3857.6 | 0.00001 | 0.00002 | 0.00003 | 0.00004 | 0.00005 | 0.0001 | 0.0005 | 0.001 |
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 30 | 60 | 90 | 120 | 150 | 300 | 1500 | 3000 |
1500 | 1371.6 | 1.5 | 3 | 4.5 | 6 | 7.5 | 15 | 75 | 150 |
2000 | 1828.8 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 1 | 5 | 10 |
2500 | 2286 | 0.01 | 0.02 | 0.03 | 0.04 | 0.05 | 0.1 | 0.5 | 1 |
3000 | 2743.2 | 0.001 | 0.002 | 0.003 | 0.004 | 0.005 | 0.01 | 0.05 | 0.1 |
3500 | 3200.4 | 0.00009 | 0.00018 | 0.00027 | 0.00036 | 0.00045 | 0.0009 | 0.0045 | 0.009 |
4000 | 3857.6 | 0.00001 | 0.00002 | 0.00003 | 0.00004 | 0.00005 | 0.0001 | 0.0005 | 0.001 |
where No is the number of the emitted fast neutrons; ri is the distance from the source; L(=constant for given distance and energy) is the relaxation length; f(=constant for given energy) is the change in the relaxation length [6,7].
It can be seen that the exponent in formula 3 can be omitted only if L≫ ri as also stated in [6]. L varies between 250 m and 850 m for neutrons in air, increasing with the neutron energy [7] and is smaller than the distances in Table 1. Calculations of the inverse square of the distance are not acceptable in this case. Based on formulas 3 and 2, we have calculated results, which are similar to those given in Table 1.
2.1.3 Selection of the Dose Criterion (the Acceptable Dose Behind the Shield).
The dose criterion here is defined as DC = 0.0001 Sv dose equivalent, equal to the usual permitted daily dose from external exposure of the operator in a controlled area in nuclear facilities without special permit. Different DCs can be defined, for example, based on the accident reference levels.
2.1.4 Dose Transmission Factor (DTF) Estimation.
DTF is defined as the ratio between the personal dose equivalent of photons or neutrons behind the shield and the personal dose equivalent in front of the shield.
Dij is the neutron personal dose equivalent in front of the shielding (here taken from Table 1).
2.1.5 Application of Appropriate Relations Between DTF and Shield Thickness.
Based on Fig. 1 in Ref. [2], calculated with the Monte Carlo N-particle transport code (MCNP code), Fig. 2 was developed, which plots the DTF curves for various materials versus shield thickness for fission neutrons. Instead of percentages of the DTF on the original figure [2] we use the corresponding numbers. It is known that the energy spectrum of neutrons leaking from a weapon is predominantly a fission spectrum [4]. The interesting material for the shelter wall is concrete. It can be seen from Fig. 2 that heavy concrete attenuates fission neutrons more effectively than concrete of lower density. From Fig. 2 we read the shelter wall thickness corresponding to the DTF of interest for the desired material, here heavy concrete with a density of 3200 kg/m3. For example, the heavy concrete shield thickness for DTF = 10−4 for fission neutrons according Fig. 2 is about 0.78 m. To achieve lower value of the DTF than those shown on the Fig. 2 (<10−5) we can multiply two or more DTFs from Fig. 2 and the resulting shield thickness will be the sum of the corresponding to every DTF wall thicknesses. (The relationship between DTF and x on the logarithmic scale in Fig. 2 is a straight line). To achieve a DTF of 10−7, we can sum 1 m, which reduces the dose with a factor 10−5, and 0.38 m, which reduces the dose with additional factor of 10−2. The resulting heavy concrete wall thickness is 1.38 m. For our purposes, we can similarly consider the curve of another material.
2.1.6 Calculation of Heavy Concrete Shield Thickness.
The wall thicknesses of the heavy concrete according to Fig. 2 based on Ref. [2], as a function of the calculated DTFs in Table 2 are listed in Table 3. Table 4 is for the case if one wants to calculate from Fig. 2 shield thickness for material other than concrete, e.g. water.
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | m | ||||||||
1000 | 914.4 | 3.33 × 10−6 | 1.67 × 10−6 | 1.11 × 10−6 | 8.33 × 10−7 | 6.67 × 10−7 | 3.33 × 10−7 | 6.67 × 10−8 | 3.33 × 10−8 |
1500 | 1371.6 | 6.67 × 10−5 | 3.33 × 10−5 | 2.22 × 10−5 | 1.67 × 10−5 | 1.33 × 10−5 | 6.67 × 10−6 | 1.33 × 10−6 | 6.67 × 10−7 |
2000 | 1828.8 | 1.00 × 10−3 | 5.00 × 10−4 | 3.33 × 10−4 | 2.50 × 10−4 | 2.00 × 10−4 | 1.00 × 10−4 | 2.00 × 10−5 | 1.00 × 10−5 |
2500 | 2286 | 1.00 × 10−2 | 5.00 × 10−3 | 3.33 × 10−3 | 2.50 × 10−3 | 2.00 × 10−3 | 1.00 × 10−3 | 2.00 × 10−4 | 1.00 × 10−4 |
3000 | 2743.2 | 1.00 × 10−1 | 5.00 × 10−2 | 3.33 × 10−2 | 2.50 × 10−2 | 2.00 × 10−2 | 1.00 × 10−2 | 2.00 × 10−3 | 1.00 × 10−3 |
3500 | 3200.4 | >1 | 5.56 × 10−1 | 3.70 × 10−1 | 2.78 × 10−1 | 2.22 × 10−1 | 1.11 × 10−1 | 2.22 × 10−2 | 1.11 × 10−2 |
4000 | 3857.6 | >1 | >1 | >1 | >1 | >1 | 1 | 2.00 × 10−1 | 1.00 × 10−1 |
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | m | ||||||||
1000 | 914.4 | 3.33 × 10−6 | 1.67 × 10−6 | 1.11 × 10−6 | 8.33 × 10−7 | 6.67 × 10−7 | 3.33 × 10−7 | 6.67 × 10−8 | 3.33 × 10−8 |
1500 | 1371.6 | 6.67 × 10−5 | 3.33 × 10−5 | 2.22 × 10−5 | 1.67 × 10−5 | 1.33 × 10−5 | 6.67 × 10−6 | 1.33 × 10−6 | 6.67 × 10−7 |
2000 | 1828.8 | 1.00 × 10−3 | 5.00 × 10−4 | 3.33 × 10−4 | 2.50 × 10−4 | 2.00 × 10−4 | 1.00 × 10−4 | 2.00 × 10−5 | 1.00 × 10−5 |
2500 | 2286 | 1.00 × 10−2 | 5.00 × 10−3 | 3.33 × 10−3 | 2.50 × 10−3 | 2.00 × 10−3 | 1.00 × 10−3 | 2.00 × 10−4 | 1.00 × 10−4 |
3000 | 2743.2 | 1.00 × 10−1 | 5.00 × 10−2 | 3.33 × 10−2 | 2.50 × 10−2 | 2.00 × 10−2 | 1.00 × 10−2 | 2.00 × 10−3 | 1.00 × 10−3 |
3500 | 3200.4 | >1 | 5.56 × 10−1 | 3.70 × 10−1 | 2.78 × 10−1 | 2.22 × 10−1 | 1.11 × 10−1 | 2.22 × 10−2 | 1.11 × 10−2 |
4000 | 3857.6 | >1 | >1 | >1 | >1 | >1 | 1 | 2.00 × 10−1 | 1.00 × 10−1 |
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | Heavy concrete (density 3200 kg/m3) shield thickness x, m | ||||||||
1000.00 | 914.4 | 1.06 | 1.10 | 1.16 | 1.19 | 1.21 | 1.26 | 1.41 | 1.46 |
1500.00 | 1371.6 | 0.83 | 0.89 | 0.92 | 0.95 | 0.99 | 1.01 | 1.15 | 1.21 |
2000.00 | 1828.8 | 0.58 | 0.65 | 0.68 | 0.72 | 0.75 | 0.78 | 0.93 | 1.00 |
2500.00 | 2286 | 0.38 | 0.45 | 0.48 | 0.52 | 0.54 | 0.58 | 0.73 | 0.78 |
3000.00 | 2743.2 | 0.18 | 0.25 | 0.27 | 0.30 | 0.33 | 0.38 | 0.53 | 0.58 |
3500.00 | 3200.4 | 0 | 0.05 | 0.06 | 0.08 | 0.12 | 0.18 | 0.33 | 0.38 |
4000.00 | 3857.6 | 0 | 0 | 0 | 0 | 0 | 0 | 0.14 | 0.18 |
kT | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|
yd. | Heavy concrete (density 3200 kg/m3) shield thickness x, m | ||||||||
1000.00 | 914.4 | 1.06 | 1.10 | 1.16 | 1.19 | 1.21 | 1.26 | 1.41 | 1.46 |
1500.00 | 1371.6 | 0.83 | 0.89 | 0.92 | 0.95 | 0.99 | 1.01 | 1.15 | 1.21 |
2000.00 | 1828.8 | 0.58 | 0.65 | 0.68 | 0.72 | 0.75 | 0.78 | 0.93 | 1.00 |
2500.00 | 2286 | 0.38 | 0.45 | 0.48 | 0.52 | 0.54 | 0.58 | 0.73 | 0.78 |
3000.00 | 2743.2 | 0.18 | 0.25 | 0.27 | 0.30 | 0.33 | 0.38 | 0.53 | 0.58 |
3500.00 | 3200.4 | 0 | 0.05 | 0.06 | 0.08 | 0.12 | 0.18 | 0.33 | 0.38 |
4000.00 | 3857.6 | 0 | 0 | 0 | 0 | 0 | 0 | 0.14 | 0.18 |
2.1.7 Alternative Calculations of the Concrete Shelter Wall Thickness for Neutron Releases.
We can use also the neutron DTFs from Ref. [1].
i.e., to achieve DTF1 from some DTF0 shown in Table 5, we need to multiply DTF0Z times and add the corresponding wall thickness Z times. It is necessary to multiply 0.229 m with Z to calculate the required wall thickness of ordinary concrete. There is a small difference in the concrete density (2200 kg/m3) in Fig. 2 and Table 5 (2340 kg/m3) as used elsewhere in Ref. [1], but the uncertainty in reading from graphs and in the data in Table 5 is dominating. The assessed wall thicknesses from application of Fig. 2 are lower compared to the application of the DTFs from Table 5, but are more accurate, because the curves in Fig. 2 were calculated with MCNP according to Ref. [2]. For neutron shielding, we recommend the wall thickness of the shelter according to Table 3 up to 1.50 m of heavy concrete with density of 3200 kg/m3 (2.27 m for ordinary concrete with density of approximately 2200 kg/m3), calculated for fission bomb yield of 1000 kT at 914.4 m distance from the explosion point. 1.50 m concrete shelter wall thickness is also recommended in Ref. [9].
Structure | Neutrons |
---|---|
Concrete blockhouse shelter with density of 2340 kg/m3 | |
0.229 m (9 in.) walls | 0.3–0.5 |
0.305 m (12 in.) walls | 0.2–0.4 |
Shelter, partly above grade | |
With 0.61 m (2 ft.) earth cover | 0.02–0.08 |
Structure | Neutrons |
---|---|
Concrete blockhouse shelter with density of 2340 kg/m3 | |
0.229 m (9 in.) walls | 0.3–0.5 |
0.305 m (12 in.) walls | 0.2–0.4 |
Shelter, partly above grade | |
With 0.61 m (2 ft.) earth cover | 0.02–0.08 |
2.2 Concrete Wall Thickness Estimation for Initial (Fission and Secondary) Gamma Ray Shielding.
From a 1-kiloton weapon the estimated production of fission gamma rays from 7.5 × 1022 fissions occurring in the last generation would produce a total of 6 × 1023 fission gamma rays having a spectrum of energies between about 0.01 and 10 MeV. Most of the important gamma rays are in the energy region from 1 to 4.5 MeV [4]. Secondary (nitrogen capture) gamma rays range in energy from about 1.5 MeV to about 11 MeV, with those between 4 and 7.5 MeV, together with a 10.83-MeV gamma ray, being the most prominent [4]. The fission product gamma-ray and the secondary gamma-ray contribute approximately equal dose at slant ranges up to about 2743 m (3000 yd.). For explosions of higher yield, however, hydrodynamic enhancement may cause the fission product gamma-ray dose to exceed the secondary gamma-ray dose, particularly for larger ranges [1].
2.2.1 Selection of Formulas and Parameters, Describing the Relation Between Shield Thickness and Dose (Rate).
D is the dose in front the shield – from Table 6 for fission gamma rays and from Table 7 for secondary gamma rays; x is the shield thickness; μ is the linear attenuation coefficient of the shielding material (Table 8), B is the build-up factor (Table 9), DC is the dose criterion.
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 0.7 | 12.6 | 28 | 36.4 | 70 | 105 | 210 | 1260 | 2100 |
2000 | 1828.8 | 0.01 | 0.18 | 0.45 | 0.7 | 1 | 1.8 | 4 | 60 | 180 |
3000 | 2743.2 | 0.0002 | 0.0036 | 0.009 | 0.016 | 0.02 | 0.036 | 0.1 | 1.4 | 7 |
3500 | 3200.4 | 0.00002 | 0.0002 | 0.0004 | 0.0006 | 0.0008 | 0.001 | 0.002 | 0.01 | 0.02 |
4000 | 3857.6 | 0.0000035 | 0.000035 | 0.00007 | 0.000105 | 0.00014 | 0.000175 | 0.00035 | 0.00175 | 0.0035 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 0.7 | 12.6 | 28 | 36.4 | 70 | 105 | 210 | 1260 | 2100 |
2000 | 1828.8 | 0.01 | 0.18 | 0.45 | 0.7 | 1 | 1.8 | 4 | 60 | 180 |
3000 | 2743.2 | 0.0002 | 0.0036 | 0.009 | 0.016 | 0.02 | 0.036 | 0.1 | 1.4 | 7 |
3500 | 3200.4 | 0.00002 | 0.0002 | 0.0004 | 0.0006 | 0.0008 | 0.001 | 0.002 | 0.01 | 0.02 |
4000 | 3857.6 | 0.0000035 | 0.000035 | 0.00007 | 0.000105 | 0.00014 | 0.000175 | 0.00035 | 0.00175 | 0.0035 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 1.7 | 17 | 34 | 51 | 68 | 85 | 170 | 850 | 1700 |
1500 | 1371.6 | 0.17 | 1.7 | 3.4 | 5.1 | 6.8 | 8.5 | 17 | 85 | 170 |
2000 | 1828.8 | 0.03 | 0.3 | 0.6 | 0.9 | 1.2 | 1.5 | 3 | 15 | 30 |
2500 | 2286 | 0.004 | 0.04 | 0.08 | 0.12 | 0.16 | 0.2 | 0.4 | 2 | 4 |
3000 | 2743.2 | 0.0007 | 0.007 | 0.014 | 0.021 | 0.028 | 0.035 | 0.07 | 0.35 | 0.7 |
3500 | 3200.4 | 0.00012 | 0.0012 | 0.0024 | 0.0036 | 0.0048 | 0.006 | 0.012 | 0.06 | 0.12 |
4000 | 3857.6 | 0.000025 | 0.00025 | 0.0005 | 0.00075 | 0.001 | 0.00125 | 0.0025 | 0.0125 | 0.025 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv | Sv |
1000 | 914.4 | 1.7 | 17 | 34 | 51 | 68 | 85 | 170 | 850 | 1700 |
1500 | 1371.6 | 0.17 | 1.7 | 3.4 | 5.1 | 6.8 | 8.5 | 17 | 85 | 170 |
2000 | 1828.8 | 0.03 | 0.3 | 0.6 | 0.9 | 1.2 | 1.5 | 3 | 15 | 30 |
2500 | 2286 | 0.004 | 0.04 | 0.08 | 0.12 | 0.16 | 0.2 | 0.4 | 2 | 4 |
3000 | 2743.2 | 0.0007 | 0.007 | 0.014 | 0.021 | 0.028 | 0.035 | 0.07 | 0.35 | 0.7 |
3500 | 3200.4 | 0.00012 | 0.0012 | 0.0024 | 0.0036 | 0.0048 | 0.006 | 0.012 | 0.06 | 0.12 |
4000 | 3857.6 | 0.000025 | 0.00025 | 0.0005 | 0.00075 | 0.001 | 0.00125 | 0.0025 | 0.0125 | 0.025 |
The linear attenuation coefficient μ is a measure of how much the incident energy beam is weakened by the material it is passing through.
Mean-free-path (mfp) is the average distance traveled by a moving particle in a target medium between interactions with the target material [2].
Thickness (mfp) =μ–1, m | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Energy (MeV) | 0.005 | 0.01 | 0.02 | 0.03 | 0.04 | 0.05 | 0.06 | 0.07 | 0.08 | 0.09 |
0.5 | 1.48 | 2.33 | 4.17 | 6.43 | 9.12 | 12.4 | 16 | 20.41 | 25.42 | 37.14 |
1 | 1.58 | 2.05 | 3.31 | 4.88 | 6.59 | 8.54 | 10.22 | 12.18 | 15.46 | 21.09 |
2 | 1.44 | 1.86 | 2.83 | 3.72 | 4.77 | 5.81 | 7.12 | 8.24 | 9.36 | 11.88 |
4 | 1.37 | 1.77 | 2.23 | 3.01 | 3.45 | 4.16 | 4.88 | 5.56 | 6.11 | 7.32 |
6 | 1.31 | 1.54 | 2.04 | 2.46 | 2.89 | 3.37 | 3.81 | 4.24 | 4.87 | 5.78 |
8 | 1.27 | 1.52 | 1.88 | 2.23 | 2.64 | 2.92 | 3.27 | 3.76 | 4.13 | 4.63 |
10 | 1.22 | 1.41 | 1.73 | 2.02 | 2.38 | 2.62 | 3.11 | 3.22 | 3.54 | 4.23 |
Thickness (mfp) =μ–1, m | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Energy (MeV) | 0.005 | 0.01 | 0.02 | 0.03 | 0.04 | 0.05 | 0.06 | 0.07 | 0.08 | 0.09 |
0.5 | 1.48 | 2.33 | 4.17 | 6.43 | 9.12 | 12.4 | 16 | 20.41 | 25.42 | 37.14 |
1 | 1.58 | 2.05 | 3.31 | 4.88 | 6.59 | 8.54 | 10.22 | 12.18 | 15.46 | 21.09 |
2 | 1.44 | 1.86 | 2.83 | 3.72 | 4.77 | 5.81 | 7.12 | 8.24 | 9.36 | 11.88 |
4 | 1.37 | 1.77 | 2.23 | 3.01 | 3.45 | 4.16 | 4.88 | 5.56 | 6.11 | 7.32 |
6 | 1.31 | 1.54 | 2.04 | 2.46 | 2.89 | 3.37 | 3.81 | 4.24 | 4.87 | 5.78 |
8 | 1.27 | 1.52 | 1.88 | 2.23 | 2.64 | 2.92 | 3.27 | 3.76 | 4.13 | 4.63 |
10 | 1.22 | 1.41 | 1.73 | 2.02 | 2.38 | 2.62 | 3.11 | 3.22 | 3.54 | 4.23 |
Gamma buildup factor for concrete depends on the material and the energy of the photons.
2.2.2 Gamma-Ray Energy Range Selection.
For fission product gamma rays we chose the important energy range around 2 MeV, for secondary gamma rays is selected the energy range around 5 MeV [4].
2.2.3 Mean-Free-Path Calculation.
The thickness of the media ranges from 0.5 to 10 mfp and gamma energy ranges from 0.5 to 10 MeV. The maximum relative error is 2% [11].
2.2.4 Information About Gamma Buildup Factors for Concrete (Table 9).
2.2.5 Selection of Barite Concrete Radiological Parameters for Gamma Rays (Table 10).
2.2.6 Fission Product Gamma-Ray.
2.2.6.1 Fission Product Gamma-Ray Component of the Initial Nuclear Radiation Dose.
The necessary correction for the differences in cloud rise velocity and hydrodynamic enhancement, each of which is a function of total energy yield, is provided by the use of effective bomb yield instead of actual bomb yield [1].
Figure 3, based on Fig. 8.130b in Ref. [1], shows the effective yield as a function of actual yield. From the plots in Fig. 3, we find the effective yield for the actual yield of interest at a given distance and multiply the registered dose per kT from Fig. 4 at that distance with the effective yield.
Here we consider air burst at high above 90 m.
2.2.6.2 Evaluation of Barite Concrete Shield Thickness for Fission Product Gamma Rays.
The shield thickness is estimated by solving Eq. (7) using the parameters given in Table 10 for barite concrete (ρ = 3490 kg/m3) and fission gamma-ray energy of 2 MeV. The results are listed in Table 11.
Gamma energy, MeV | μ, m−1 | mfp = 1/μ, m | Buildup factor B |
---|---|---|---|
5 | 18.76 | 0.0533 | ∼4 |
2 | 11.21 | 0.0892 | ∼10 |
Gamma energy, MeV | μ, m−1 | mfp = 1/μ, m | Buildup factor B |
---|---|---|---|
5 | 18.76 | 0.0533 | ∼4 |
2 | 11.21 | 0.0892 | ∼10 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m |
1000 | 914.4 | 0.9952 | 1.2530 | 1.3243 | 1.3477 | 1.4060 | 1.4422 | 1.5040 | 1.6639 | 1.7094 |
2000 | 1828.8 | 0.6162 | 0.8741 | 0.9558 | 0.9952 | 1.0270 | 1.0795 | 1.1507 | 1.3923 | 1.4903 |
3000 | 2743.2 | 0.2672 | 0.5251 | 0.6068 | 0.6581 | 0.6780 | 0.7305 | 0.8216 | 1.0570 | 1.2006 |
3500 | 3200.4 | 0.0618 | 0.2672 | 0.3291 | 0.3652 | 0.3909 | 0.4108 | 0.4726 | 0.6162 | 0.6780 |
4000 | 3857.6 | 0 | 0.1118 | 0.1736 | 0.2098 | 0.2354 | 0.2553 | 0.3172 | 0.4607 | 0.5226 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m |
1000 | 914.4 | 0.9952 | 1.2530 | 1.3243 | 1.3477 | 1.4060 | 1.4422 | 1.5040 | 1.6639 | 1.7094 |
2000 | 1828.8 | 0.6162 | 0.8741 | 0.9558 | 0.9952 | 1.0270 | 1.0795 | 1.1507 | 1.3923 | 1.4903 |
3000 | 2743.2 | 0.2672 | 0.5251 | 0.6068 | 0.6581 | 0.6780 | 0.7305 | 0.8216 | 1.0570 | 1.2006 |
3500 | 3200.4 | 0.0618 | 0.2672 | 0.3291 | 0.3652 | 0.3909 | 0.4108 | 0.4726 | 0.6162 | 0.6780 |
4000 | 3857.6 | 0 | 0.1118 | 0.1736 | 0.2098 | 0.2354 | 0.2553 | 0.3172 | 0.4607 | 0.5226 |
2.2.7 Secondary Gamma-Ray
2.2.7.1 Assessment of Secondary Gamma-Ray Component of the Initial Nuclear Radiation Dose.
As pointed out in Ref. [4] in some radiation transport calculations for “free-field” radiation doses produced by predominantly fission weapons, gamma rays produced by neutron capture in the nitrogen of the atmosphere, secondary gamma rays, have been the only gamma rays considered other than fission-product gamma rays.
Air-secondary gamma-ray component of the initial nuclear radiation dose per kiloton yield as a function of slant range from fission weapon airbursts is given in Fig. 5, developed based on Fig. 8. 127a, defense curve, from Ref. [1]. Corrections of the actual doses in this case are not necessary.
The data from the defense curve is selected as more conservative.
2.2.7.2 Selection of Important Energy Diapason.
For secondary gamma rays, the energy diapason around 5 MeV.
2.2.7.3 Shelter Wall Thickness Estimation for Secondary Gamma-Ray Component.
The shield thickness x is estimated by solving Eq. (7) using the parameters given in Table 10 for a secondary gamma-ray energy of 5 MeV. The results are listed in Table 12.
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m |
1000 | 914.4 | 0.5931 | 0.7159 | 0.7528 | 0.7744 | 0.7898 | 0.8017 | 0.8386 | 0.9244 | 0.9614 |
1500 | 1371.6 | 0.4704 | 0.5931 | 0.6301 | 0.6517 | 0.6670 | 0.6789 | 0.7159 | 0.8017 | 0.8386 |
2000 | 1828.8 | 0.3779 | 0.5007 | 0.5376 | 0.5592 | 0.5746 | 0.5865 | 0.6234 | 0.7092 | 0.7462 |
2500 | 2286 | 0.2705 | 0.3933 | 0.4302 | 0.4518 | 0.4672 | 0.4791 | 0.5160 | 0.6018 | 0.6387 |
3000 | 2743.2 | 0.1776 | 0.3004 | 0.3373 | 0.3589 | 0.3743 | 0.3862 | 0.4231 | 0.5089 | 0.5458 |
3500 | 3200.4 | 0.0836 | 0.2064 | 0.2433 | 0.2649 | 0.2803 | 0.2921 | 0.3291 | 0.4149 | 0.4518 |
4000 | 3857.6 | 0.00 | 0.1227 | 0.1597 | 0.1813 | 0.1966 | 0.2085 | 0.2455 | 0.3313 | 0.3682 |
kT | 1 | 10 | 20 | 30 | 40 | 50 | 100 | 500 | 1000 | |
---|---|---|---|---|---|---|---|---|---|---|
yd. | m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m | x, m |
1000 | 914.4 | 0.5931 | 0.7159 | 0.7528 | 0.7744 | 0.7898 | 0.8017 | 0.8386 | 0.9244 | 0.9614 |
1500 | 1371.6 | 0.4704 | 0.5931 | 0.6301 | 0.6517 | 0.6670 | 0.6789 | 0.7159 | 0.8017 | 0.8386 |
2000 | 1828.8 | 0.3779 | 0.5007 | 0.5376 | 0.5592 | 0.5746 | 0.5865 | 0.6234 | 0.7092 | 0.7462 |
2500 | 2286 | 0.2705 | 0.3933 | 0.4302 | 0.4518 | 0.4672 | 0.4791 | 0.5160 | 0.6018 | 0.6387 |
3000 | 2743.2 | 0.1776 | 0.3004 | 0.3373 | 0.3589 | 0.3743 | 0.3862 | 0.4231 | 0.5089 | 0.5458 |
3500 | 3200.4 | 0.0836 | 0.2064 | 0.2433 | 0.2649 | 0.2803 | 0.2921 | 0.3291 | 0.4149 | 0.4518 |
4000 | 3857.6 | 0.00 | 0.1227 | 0.1597 | 0.1813 | 0.1966 | 0.2085 | 0.2455 | 0.3313 | 0.3682 |
To determine the effectiveness of a given shield as a protection against the initial radiation from a nuclear explosion, it is recommended to use the higher of the values for fission and secondary gamma rays [1].
2.2.8 Derivation of Shelter Wall Thickness for Different Materials, Based on Results for Concrete.
2.2.9 Alternative Calculations of the Concrete Shelter Wall Thickness for Initial Gamma Rays.
2.2.9.1 Tenth-Value Thickness.
A tenth-value thickness is defined as the thickness of the specified material which reduces the radiation dose or dose rate to one tenth of that falling upon it [1]. The effective tenth-value thickness of some materials of interest for broad beams of gamma rays emitted by the fission products in the first minute after the detonation and for secondary gamma rays are given in Ref. [1]. For example, the approximate tenth-value shield thickness for concrete (density 2340 kg/m3) is 0.2794 m (11 in.) [1]. We can compare the calculation results in Table 11 with the results from tenth-value thickness application using Eq. (5) with DTF0 equal to 0.1 and make the conclusion that both results are similar and expected given the differences in the concrete density.
2.2.9.2 Estimates of the Initial Gamma-Ray Shielding Afforded in Terms of a “Dose Transmission Factor”.
Concrete shield thickness estimates (density 2340 kg/m3) can be found using the DTF between 0.1 and 0.2 for a 0.2286 m (9-in.) shelter wall thickness [1]. For a shelter partially above grade, we use a DTF of 0.03-0.07, corresponding to a wall thickness of 0.6096 m (2 ft.) [1]. We can compare the calculation results in Table 11 with results from application of the DTF between 0.1 and 0.2 using Eq. (5). Both results are similar and expected given the differences in the concrete density. Underground shelter with an earth cover will have dose reduction factor orders of magnitude lower than the concrete blockhouse shelter.
2.3 Fallout Radiation.
Radioactive fallout contaminates large areas and is an immediate and extreme biological hazard. Radioactive fallout emits alpha, beta, and gamma radiation. Estimates of attenuation of residual radiation in various structures have been made, based partly on calculations and partly on measurements with simulated fallout [1]. According to Ref. [1] blockhouse shelter wall thickness of 0.2286 m (9 in.) for protection from fallout radiation corresponds to DTF of 0.007–0.09. Total-dose contours from early fallout at 1, 6, and 18 h after a surface burst with a total yield of 2 megatons and 1-megaton fission yield (24.14 km/h (15 mph) effective wind speed) are given in Fig. 9.86b in Ref. [1]. At a distance r of about 3000–4000 m from the bomb burst according to Ref. [1] the expected total dose from fallout in the first hour is 10 Sv. To estimate the thickness of the concrete shield to protect against exposure to 10 Sv from fallout radiation, we use the dose criterion DC = 0.0001 Sv and calculate DTF1 = 0.0001/10 = 10−5. We conservatively select DTF0 = 0.09 with corresponding concrete wall thickness of 0.2286 m. According to Eq. (5) the number Z to multiply the wall thickness of 0.2286 m (9 in.) is equal to 4.78. The necessary concrete roof and wall thickness x for protection from fallout gamma-ray radiation from 1 MT fission yield is 1.093 m.
2.3.1 Rule of Thumb.
The total dose increases with time over a given distance, due to continued deposition of radioactive substances, but dose rate decreases with time due to the radioactive decay and the corresponding exposure rate is expected to decrease as well. Based on empirical measurements of fallout radiation with dose-rate meters the 7:10 Rule of thumb has been derived [13]. The 7:10 Rule of Thumb states that for every 7-fold increase in time after detonation, there is a 10-fold decrease in the exposure rate, where the rate is the same unit as the time increase; For 7 × 7 × 7 = 343 h (two weeks) the dose rate will decrease 10 × 10 × 10 = 1000 times. For 7 × 7 = 49 h the dose rate will decrease 10 × 10 = 100 times [1,9]. The rule assumes that the time of detonation is known and that fallout from only one detonation is present in relatively significant quantities. For accuracy and reliability, nothing can replace a direct instrument reading. However, depending on circumstances, it may become necessary to apply this general rule to predict when conditions may allow short trips outside a shelter [1].
3 Conclusions
This practical guideline proposes a simplified approach to estimate the required thickness of heavy concrete wall (density 3200 kg/m3) for shielding against fission neutrons and initial gamma rays. The expected neutron and initial gamma ray doses per kT yield are taken from reference source [1] at different distances and multiplied by the bomb yield. The required dose behind the shield DC is defined as 0.0001 Sv = 0.1 mSv. The required DTFs are calculated as the ratio of dose behind the shield and dose without shielding.
3.1 Neutron Shielding.
Heavy concrete shield thicknesses calculated with MCNP corresponding to a specific DTF were taken from reference source [2]. The results are tabulated in relation to the distance and bomb yield. A comparison is made based on DTFs from a reference source [1] and both methods were found to give similar results given the differences in the concrete density and accuracy of the methods. For neutron shielding is recommended a wall about 1.50 m thick made of heavier concrete with density ≥ 3200 kg/m3 is assessed for 1000 kT fission bomb yield at a distance of 914.4 m (1000 yd.).
3.2 Shielding of Initial (Fission and Secondary) Gamma Rays.
Fission product gamma ray doses are taken from the reference source [1] and they are corrected with the effective yield as a function of actual yield [1]. For secondary gamma rays no correction is necessary according to Ref. [1]. Relation between the doses in-front and behind the shield and shield thickness is used according to Eq. (6). Parameters of the function are chosen for barite concrete and the corresponding gamma energy. The results were compared with the application of tenth-value thickness and DTFs for fission and secondary gamma rays found in Ref. [1]. The results of applying the analytical function, tenth-value thickness, and empirical DTFs are similar. Based on all methods applied, a 1.70 m thick protective wall made of barite concrete (density 3420 kg/m3) assessed for 1000 kT fission bomb yield at a distance of 914.4 m (1000 yd.) is recommended for protection against fission gamma rays.
3.3 Fallout Radiation Shielding.
From the reference source [1] are taken DTFs and the corresponding concrete wall thicknesses. The calculated thickness of the walls for fallout radiation protection is about 1.10 m of ordinary concrete with density of 2340 kg/m3.
3.4 Recommended Thickness of the Wall of a Nuclear Shelter.
Based on the results of the evaluation of the shield thickness against fission neutrons, initial (fission and secondary) gamma rays, and fallout radiation, the most conservative value of 1.70 m of barite concrete with a density of 3420 kg/m3 is recommended, calculated for fission gamma rays for a 1000 kT fission bomb yield at a distance of 914.4 m (1000 yd.).
Manually operated shelter ventilation systems must be installed to protect against internal exposure from fallout radiation. Shelter wall openings must be accordingly shielded, for example, with armored doors, labyrinths, movable screens, etc. The provisions, water, air, etc. shall be ensured at least for the first two weeks for the time until the fallout radiation is reduced naturally or by organizational actions.
Nomenclature
- B =
buildup factor
- Bconcr =
buildup factor for concrete
- Bsoil =
buildup factor for soil
- D =
dose in front of the shield, Sv
- DC =
dose criterion, Sv
- Dconcr =
dose behind the concrete shield, Sv
- Dsoil =
dose behind the soil shield, Sv
- DTF =
dose transmission factor
- DTF0 =
empirical DTF
- DTF1 =
requested DTF
- E =
neutron or photon energy, J
- f =
factor of change in the relaxation length
- L =
relaxation length, m
- mfp =
mean-free-path, m
- N0 =
number of the emitted fast neutrons
- r =
distance, m
- t =
time, s
- x =
shield thickness, m
- xconcr =
thickness of the concrete shield, m
- xsoil =
thickness of the soil shield, m
- Y =
bomb yield, kg
- Z =
represents the relationship between two DTfs