Интернет - Библиотека
Нормативная правовая документация

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

БЕЗОПАСНОСТЬ РАДИАЦИОННАЯ ЭКИПАЖА
КОСМИЧЕСКОГО АППАРАТА
В КОСМИЧЕСКОМ ПОЛЕТЕ

МЕТОДИКА РАСЧЕТА ЭФФЕКТИВНОСТИ
ЭЛЕКТРОСТАТИЧЕСКОЙ ЗАЩИТЫ ПРИ РАДИАЦИОННОМ
ВОЗДЕЙСТВИИ

 

РД 50-25645.221-90

 

 

 

ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО УПРАВЛЕНИЮ
КАЧЕСТВОМ ПРОДУКЦИИ И СТАНДАРТАМ

Москва 1991

 

РУКОВОДЯЩИЙ НОРМАТИВНЫЙ ДОКУМЕНТ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

БЕЗОПАСНОСТЬ РАДИАЦИОННАЯ ЭКИПАЖА
КОСМИЧЕСКОГО АППАРАТА В КОСМИЧЕСКОМ
ПОЛЕТЕ

Методика расчета эффективности электростатической
защиты при радиационном воздействии

РД
50-25645.221-90

Дата введения 01.01.92

Настоящие методические указания устанавливают алгоритм расчета эффективности электростатической защиты космического аппарата (КА) при воздействии на него электронов радиационных поясов Земли (РПЗ).

Методика предназначена для расчетов дозовых нагрузок на биологические и технические объекты за электростатической защитой (ЭСЗ) и определения показателя ее эффективности.

1. ОБЩИЕ ПОЛОЖЕНИЯ. ВХОДНАЯ ИНФОРМАЦИЯ ДЛЯ ПРОВЕДЕНИЯ РАСЧЕТОВ

1.1. Методика оценки эффективности ЭСЗ при радиационном воздействии основана на учете влияния защитного электростатического поля на динамику движения электронов и моделировании методом Монте-Карло переноса в веществе конструкции ЭСЗ электронов остаточного потока и генерируемого им тормозного излучения.

1.2. Эффективность ЭСЗ (μэфф ), см2/г, характеризуют показателем, представляющим собой аналог эффективного коэффициента ослабления

                                                                   (1)

где H - мощность дозы, создаваемой электронами или тормозным излучением в расчетной точке поглотителя, расположенной непосредственно за ЭСЗ, Зв/с;

δЭСЗ - приведенная толщина конструкций ЭСЗ, г/см2.

1.3. Под ЭСЗ в методических указаниях понимают высоковольтное электрофизическое устройство в виде набора идентичных секций, каждая из которых представляет собой электростатический конденсатор, образованный двумя тонкими плоскопараллельными электродами, разделенными вакуумной изоляцией. Конструктивно набор секций закреплен на опорных изоляторах цилиндрической формы.

1.4. При описании геометрических характеристик модели ЭСЗ и при моделировании движения электрона в поле защиты используют декартову систему координат, привязанную к точке пересечения оси симметрии модели ЭСЗ с внешним (верхним) электродом (стартовая система координат), а при моделировании прохождения электрона и тормозного γ-кванта в веществе конструкций ЭСЗ - локальную декартову систему координат, привязанную к координатам точки их входа в вещество.

1.5. Геометрические параметры модели ЭСЗ задают в соответствии с ГОСТ 25645.204 (разд. 2), описывая набор высоковольтных электродов и изоляторов совокупностью зон с постоянными физическими свойствами в пределах зон.

1.5.1. Зоны, соответствующие высоковольтным электродам, задают плоскостями, параллельными плоскости XY в декартовой системе координат.

1.5.2. Зоны, соответствующие изоляторам, задают цилиндрическими поверхностями, параллельными оси Z.

1.6. Атомный номер ядер вещества (z) высоковольтных электродов задают в диапазоне 6 - 13, а атомный номер ядер вещества изоляторов в диапазоне 13 - 18.

1.7. Распределение потенциала (U), В, в модели ЭСЗ описывают функцией

                                      (2)

где E = ΔU/dc - напряженность электрического поля в секции модуля ЭСЗ, В/м;

ΔU - разность потенциала между двумя соседними электродами секции, В;

dc - расстояние между двумя соседними электродами секции, м;

dм - высота модуля, м.

Число секций модуля (Nс) определяют по формуле

Nc = 2Uраб/ΔU,                                                                     (3)

где Uраб - рабочее напряжение ЭСЗ, максимальное значение которого не должно превышать 1/e · Tmax, В;

e - заряд электрона, Кл;

Tmax - максимальная энергия падающих электронов, Дж.

1.8. Дифференциальное по энергии распределение плотности потока электронов, падающих на внешний электрод модели ЭСЗ, задают в соответствии с действующими методиками.

1.9. Моделирование процесса переноса электронного излучения в веществе, моделирование процессов рождения тормозных γ-квантов, а также переноса γ-квантов осуществляют с использованием предварительно подготовленных массивов и параметров, алгоритмы вычисления которых даны в приложении.

1.9.1. При моделировании отрезка траектории электрона (толщины слоя Мольер) tМ используют массив параметров С, элементы которого определены для тех значений энергии TМ, для которых справедливо соотношение (28), выполняемое в интервале энергий (TМ(i-1), TМ(i)), где 0,06 МэВ < TМ(i) <10 МэВ.

1.9.2. При моделировании процесса углового рассеяния по теории Мольер используют массив интегральных распределений Pυ для розыгрыша параметра υ в теории Мольер, а также массив Bм вспомогательного параметра B, входящего в выражение для угла рассеяния (ΘМ) электрона после прохождения слоя tМ. Массивы Pυ и Bм рассчитывают для значений TМ из упомянутого выше интервала энергий.

1.9.3. При моделировании ионизационных потерь энергии ΔT электроном в слое Мольер используют распределение вероятности Pb, b2) для розыгрыша безразмерного параметра b2, через который по теории Блунка и Лейзенганга определяют значение ΔT. Массив Рb рассчитывают для нескольких значений параметра b2, который является функцией средних ионизационных потерь, толщины слоя и атомного номера вещества. Значения параметров b2 для диапазона значений энергии электрона TМ определяют из массива данных .

1.9.4. При моделировании процесса образования γ-кванта используют массив данных по макроскопическим сечениям Mγ для расчета вероятности образования γ-кванта при прохождении электроном слоя tМ. Для розыгрыша энергии γ-кванта используют распределение вероятности Pк(K, T), рассчитанное для пяти значений энергий электронов.

2. МЕТОДИКА РАСЧЕТА ЭФФЕКТИВНОСТИ ЭЛЕКТРОСТАТИЧЕСКОЙ ЗАЩИТЫ ПРИ РАДИАЦИОННОМ ВОЗДЕЙСТВИИ

2.1. Учет влияния отклоняющего электрического поля на движение электрона в единичной секции модуля ЭСЗ

2.1.1. Проверяют условия принадлежности траектории движения электрона конусу «разрешенных» направлений движения,

                           (4)

то электрон отклоняется полем (вне конуса);

если

                           (5)

то электрон достигает высоковольтного электрода (внутри конуса),

где Θe - угол вхождения электрона в поле относительно нормали к электроду, рад;

e - заряд электрона, Кл;

Tн - начальная энергия электрона, Дж;

mc2 - энергия покоя электрона, Дж.

2.1.2. Определяют координаты вхождения электрона в высоковольтный электрод после прохождения им защитного электрического поля в плоскости движения электрона.

                              (6)

                         (7)

где c - скорость света, м/с;

- начальный импульс электрона, кг · м · с-1;

P0Z = P0 · cosΘe;

P0y = P0sinΘe;

t - время движения электрона между двумя электродами, с;

2.1.3. Определяют угол вхождения электрона (Θвх), рад, в высоковольтный электрод относительно нормали к электроду по формуле

                        (8)

2.1.4. Определяют энергию электрона (Tп), Дж, после прохождения области электрического поля по формуле

Tп = Tн + |e|Edc.                                                                    (9)

2.1.5. Определяют координаты вхождения электрона в электрод в системе координат, привязанной к стартовой системе координат по формулам

Xгл = Xн + Xвх · sinφ,                                                            

Yгл = Yн + Yвх · cosφ,                                                             (10)

Zгл = Zн + Zвх,                                                                       

где Xн, Yн, Zн - начальные (стартовые) координаты электрона, м;

φ - азимутальный угол в декартовой системе координат, рад.

2.2. Определение толщины слоя Мольер (tм), г/см2, проводят по формуле

tМ = Cм(T)(1,13 + 3,76γ2)/γ2,                                               (11)

где Cм(T) определяют из массива данных Cм путем интерполирования при заданном значении энергии T;

β2 = (2TE0 + T2)/(T + E0)2,

E0 = mc2.

2.3. Определение угла рассеяния после прохождения слоя Мольер

2.3.1. Определяют значение вспомогательных параметров B из массива Bм путем интерполирования по T и χc, рад, по формуле

                       (12)

где A - атомный вес материала электрода.

2.3.2. Для текущего значения случайного числа из датчика случайных чисел, равномерно распределенных на интервале (0,1), методом кусочно-линейной интерполяции определяют υ из массива данных Pυ.

2.3.3. Определяют угол рассеяния по формуле

2.4. Определяют направление движения электрона  относительно нормали  к слою толщиной tМ после i-го шага моделирования переноса электрона по соотношениям:

cosφi = (c1sinφi-1 + c2cosφi-1)/sinΘi;                                     (13)

sinφi = (-c1cosφi-1 + c2sinφi-1)/sinΘi,

где c1 = -sinΘМ · sinφ'i; c2 = sinΘi-1cosΘМ + cosΘi-1sinΘМcosφ'.

Считают, что значения угла φ' распределены равномерно в интервале (0 - 2p). Определение значения случайного угла проводят с помощью датчика случайных чисел x по формуле φ' = 2px.

2.5. Определение ионизационных потерь энергии электроном в слое толщиной tМ

2.5.1. Определяют значение вспомогательного параметра b2 по формуле

                                              (14)

где  - средние потери энергии на ионизацию на отрезке tМ, МэВ.

2.5.2. Из массива данных  находят два ближайших значения b12 и b22 (b12b2b22).

2.5.3. Из массива распределения вероятностей Pв по случайному числу находят методом кусочно-линейной интерполяции параметры λb1, λb2, соответствующие b12 и b22, и затем методом интерполяции по b2 определяют λb.

2.5.4. Потерю энергии электрона с начальной энергией T на ионизацию при прохождении им слоя tМ определяют по формуле

                            (15)

2.6. Определение вероятности образования и энергии тормозного γ-кванта

2.6.1. Вероятность образования тормозного γ-кванта при прохождении электроном с энергией T слоя вещества толщиной tМ определяют соотношением

P = 1 - exp(-Θγ),                                                                   (16)

где Θγ = μγ(T) · tМ.

Значение μγ(T) определяют из массива μγ линейной интерполяцией по T.

2.6.2. Значение энергии γ-кванта определяют методом двойной линейной интерполяции по случайному числу x и по энергии электрона Т из массива Pк.

Образовавшемуся γ-кванту с энергией Еγ приписывают статистический вес Wγ, равный

μγ(T) · tМ.                                                                             (17)

2.6.3. Угол, определяющий направление движения рожденного γ-кванта, относительно направления движения электрона перед прохождением слоя Мольер принимают равным углу рассеяния электрона ΘМ.

2.7. Определение параметров γ-кванта при его прохождении в веществе

2.7.1. Пробег γ-кванта до взаимодействия определяют по формуле

                                                                           (18)

где x - случайное число из равномерного распределения чисел на промежутке (0 - 1); μ(Eγ) - коэффициент ослабления γ-квантов, см2/2.

2.7.2. Вид взаимодействия в конце пробега γ-кванта определяют из условий:

при  - имеет место комптон-эффект;

при  - имеет место фотоэффект,                                                  (19)

где σк(Еγ) и σф(Eγ) - макроскопические сечения комптоновского взаимодействия и фотоэффекта соответственно.

Далее прохождение γ-кванта рассматривают только при комптоновском взаимодействии.

2.7.3. Координаты точки взаимодействия γ-кванта с веществом определяют по формулам:

Xγ(b) = Xγ - λcosφγsinΘγ,

Yγ (b) = Yγ - λsinφγsinΘγ,                                                         (20)

Zγ(b) = Zγ + λcosΘγ,

где Θγ - угол, образуемый γ-квантом с осью Z, рад;

φγ - азимутальный угол, рад.

2.8. Определение характеристик комптоновского взаимодействия

2.8.1. Энергию γ-кванта после комптоновского взаимодействия определяют по следующим соотношениям:

(21)

где K0 и K - энергии первичного и рассеянного γ-кванта в единицах m0c2;

S = K0/(1 + 0,5625K0);

x - случайное число.

2.8.2. Угол рассеянного γ-кванта (после комптоновского рассеяния) Θ'γ относительно направления движения γ-кванта с энергией K0 определяют по формуле

                                                            (22)

2.8.3. Угол вылета образовавшегося электрона Θе и его энергию Tе определяют по формулам:

                                                             (23)

Te = (K0 - K)mc2,                                                                   (24)

где

2.9. Определение функций распределения электронов и тормозных γ-квантов за ЭСЗ

2.9.1. Функции распределения ту и nγ строят с использованием дифференциальных токовых величин , , получающихся в результате моделирования траектории электронов и тормозных γ-квантов. Связь между функцией распределения и дифференциальной плотностью тока через поверхность (разделяющую ЭСЗ и тканеэквивалентный поглотитель, например в точке А) определяют по соотношениям:

 - для электронов,      (25)

 - для γ-квантов.       (26)

2.9.2. Построение функций ne, nγ сводят к процедуре построения двумерных гистограмм с интервалом по энергии ΔTe = Tmax/Ne и ΔTγ = Tmax/Nγ. Здесь Tmax - максимальная энергия электрона РПЗ. Ne выбирают из интервала чисел (25 - 50), а Nγ - из интервала (100 - 200).

2.10. Мощность дозы электронного  и тормозного  излучений в заданной точке  за ЭСЗ определяют с использованием функций распределения  и  по формулам:

                           (27)

где  - средние массовые потери энергии на ионизацию в веществе поглотителя, МэВ · см2/г;

μ(Eγ) - массовый коэффициент передачи энергии для γ-квантов с энергией Еγ, см2/г;

kn - переходной коэффициент от поглощенной энергии к дозе, равный 1,6 · 10-10;

ρ - плотность поглотителя, г/см3.

2.11. Определение показателя эффективности ЭСЗ

Серией расчетов суммарных значений мощностей дозы электронного и тормозного излучений  за ЭСЗ, проводимых при различных значениях Uраб(δЭСЗ), т.е. при различном приведенном весе (ЭСЗ), устанавливают кривую ослабления дозы - зависимость  от веса ЭСЗ, т.е.

Значение показателя эффективности в соответствии с п. 1.2 определяют по формуле (1)

ПРИЛОЖЕНИЕ

Обязательное

ОПРЕДЕЛЕНИЕ ВСПОМОГАТЕЛЬНЫХ МАССИВОВ Cм, Bм, Pυ, Mb2, Pb, Mγ и Pк

1. Подготовка данных для моделирования углового рассеяния

Устанавливают интервалы энергий, в пределах которых можно считать параметр B постоянным. Это выполнимо для тех интервалов энергии T, в пределах которых выполняются условия:

                            (28)

где tМ - толщина слоя Мольер;

β2 = (2TE0 + T2)/(T + E0)2;

Е0 = mc2 - энергия покоя электрона.

Совокупность констант См(Т) составляет массив См.

Значение B находят по формуле

B = ln(1, 1min(1, 4m)).                                                         (30)

Здесь m - среднее число актов рассеяния при прохождении электроном слоя вещества tМ. Значение m находят по формуле

m = 6680z1/3(z + 1)tМ/[β2A(1 + 3,35γ2)],                              (31)

где A - атомный номер;

γ2 = z2/(137β)2.

Совокупность значений B составляет массив данных Bм, а соответствующие им величины T составляют массив данных Tм.

Для установленных интервалов энергий строят распределения

                                                          (33)

где fM = F0(υ) + B-1F1(υ) + B-2F2(υ) - табулированные функции, представленные в табл. 1.

Значение υmax выбирают из интервала 7 - 10.

Совокупность значений Pυ(υ) составляет массив данных Pυ.

Таблица 1

Константы F0, F1, F2

υ

F0

F1

F2

0

2,00 E00

8,456 E-1

2,492 E00

0,2

1,92 E00

7,032 E-1

2,069 E00

0,4

1,72 E00

3,437 E-1

1,048 E00

0,6

1,40 E00

-7,770 E-2

-4,400 E-3

0,8

1,05 E00

-3,981 E-1

-6,068 E-1

1,0

7,34 E-1

-5,285 E-1

-6,359 E-1

1,2

4,74 E-1

-4,770 E-1

-3,086 E-1

1,4

2,82 E-1

-3,183 E-1

+5,250 E-2

1,6

1,54 E-1

-1,396 E-1

+2,423 E-1

1,8

7,83 E-2

-6,000 E-4

+2,386 E-1

2,0

3,66 E-2

+7,820 E-2

+1,316 E-1

2,2

1,58 E-3

+1,054 E-1

+1,960 E-2

2,4

6,30 E-3

+1,008 E-1

-4,670 E-2

2,6

2,32 E-3

+8,262 E-2

-6,460 E-2

2,8

7,90 E-4

+6,247 E-2

-5,460 E-2

3,0

2,50 E-4

+4,550 E-2

-3,568 E-2

3,2

7,30 E-5

+3,288 E-2

-1,923 E-2

3,4

1,90 E-5

+2,402 E-2

-8,470 E-3

3,6

4,70 E-5

+1,791 E-2

-2,640 E-3

3,8

1,10 E-6

+1,366 E-2

+5,000 E-5

4,0

2,30 E-7

+1,064 E-2

+1,074 E-3

4,5

3,00 E-9

+6,140 E-3

+1,229 E00

5,0

2,00 E-11

+3,831 E-3

+8,326 E-1

5,5

2,00 E-13

+2,527 E-3

+5,368 E-1

6,0

5,00 E-16

+1,739 E-3

+3,495 E-1

7,0

1,00 E-21

+9,080 E-4

+1,584 E-1

2. Подготовка данных для моделирования ионизационных потерь

Для интервалов энергий, установленных в п. 1 приложения, вычисляют значения параметра b2 по формуле

                                             (34)

где  - средние потери энергии на ионизацию, на пути tМ, МэВ.

Совокупность параметров b2 составляет массив данных Mb2.

Для рассчитанных параметров b2 строят распределения

                                                      (35)

где f(λ) - функция Блунка и Лайзенганга, описывающая разброс ионизационных потерь

                                  (36)

где ci, γi, λi - постоянные, представленные в табл. 2.

Таблица 2

Значение постоянных ci, γi, λi (i = 1, 2, 3, 4)

i

1

2

3

4

ci

0,174

6,058

0,019

0,007

γi

1,8

2,0

3,0

5,0

λi

0,0

3,0

6,5

11,0

Совокупность значений Pb2(λ, b2) составляет массив данных Pb.

3. Подготовка данных для моделирования рождения тормозных γ-квантов

3.1. Значения параметра (μγ(T)), см2/2, используемые при определении вероятности образования тормозного γ-кванта, определяют из соотношения

  (37)

где NA - число Авогадро;

z, A - атомный номер и атомная масса рассеивающих атомов;

R, L, Н - постоянные, зависящие от энергии электрона Т и атомного номера, табулированные в табл. 3;

Kmin, Kmax - минимальная и максимальная энергии γ-кванта, Kmax = 0,98Т.

Совокупность значений μi составляет массив данных Мγ.

3.2. Для моделирования энергии γ-кванта подготавливают распределения

                                                    (38)

где в качестве сечения тормозного излучения используют выражение

где R, L, Н - постоянные, зависящие от T и z, приведенные в табл. 3, T = 0,06, 0,4, 1 и 10 МэВ.

Совокупность значений Pк(K, T) составляет массив данных Pк.

Таблица 3

Значения коэффициентов в формуле для dσ/dK

z

T, МэВ

R

L

H

4

0,025

8,885

0,74

11,291

 

0,050

10,230

0,70

11,800

 

0,075

10,833

0,63

12,278

 

0,100

11,159

0,60

12,396

 

0,200

11,600

0,53

12,515

 

0,300

11,900

0,49

12,600

 

0,400

12,300

0,48

12,800

 

0,500

12,500

0,47

13,000

 

0,600

12,600

0,47

13,000

 

0,800

12,900

0,45

13,300

 

1,000

13,230

0,47

13,600

 

1,200

13,400

0,48

13,790

 

1,400

13,450

0,49

13,870

 

1,600

13,600

0,49

13,030

 

1,800

13,650

0,51

14,120

 

2,000

13,700

0,52

14,230

 

5,000

14,150

0,54

15,110

 

10,000

12,694

0,67

15,084

 

15,000

12,456

0,74

14,993

 

20,000

12,000

0,79

14,924

 

25,000

11,600

0,83

14,850

 

30,000

11,400

0,67

14,850

13

0,025

4,933

0,82

9,503

 

0,050

6,861

0,68

10,488

 

0,075

7,931

0,65

10,916

 

0,100

8,639

0,62

11,186

 

0,200

10,053

0,57

11,636

 

0,300

10,746

0,54

11,879

 

0,400

11,122

0,53

12,012

 

0,500

11,441

0,52

12,179

 

0,600

11,550

0,50

12,220

 

0,800

11,900

0,50

12,490

 

1,000

12,000

0,51

12,600

 

1,200

12,100

0,51

12,700

 

1,400

12,300

0,54

12,880

 

1,600

12,500

0,54

13,090

 

1,800

12,500

0,56

13,090

 

2,000

12,550

0,56

13,210

 

5,000

12,800

0,65

13,948

 

10,000

11,900

0,79

13,700

 

15,000

11,400

0,83

13,700

 

20,000

11,000

0,97

13,000

 

25,000

10,700

1,04

13,000

 

30,000

10,500

1,06

13,000

26

0,025

2,933

1,04

8,178

 

0,050

4,663

0,74

9,374

 

0,075

5,891

0,68

10,027

 

0,100

6,743

0,64

10,413

 

0,200

8,734

0,58

11,201

 

0,300

9,725

0,55

11,554

 

0,400

10,139

0,54

11,659

 

0,500

10,562

0,54

11,796

 

0,600

10,866

0,54

11,961

 

0,800

11,277

0,53

12,209

 

1,000

11,547

0,54

12,398

 

1,200

11,608

0,55

12,429

 

1,400

11,740

0,56

12,552

 

1,600

11,832

0,57

12,653

 

1,800

11,920

0,58

12,758

 

2,000

12,042

0,58

12,911

 

5,000

12,100

0,73

13,144

 

10,000

11,107

0,79

12,979

 

15,000

10,900

1,00

12,900

 

20,000

10,100

1,02

12,800

 

25,000

9,700

1,02

12,800

 

30,000

9,600

1,03

12,800

 

ИНФОРМАЦИОННЫЕ ДАННЫЕ

1. РАЗРАБОТАН И ВНЕСЕН Министерством здравоохранения СССР

РАЗРАБОТЧИКИ

А.А. Волобуев; В.А. Гончарова; А.И. Григорьев, д-р мед. наук; Е.Е. Ковалев, д-р. техн. наук; Т.Я. Королькова; Е.Н. Лесновский, канд. техн. наук; В.А. Панин; Е.В. Пашков, канд. техн. наук; В.М. Петров, канд. физ.-мат. наук; Ю.В. Потапов, канд. физ.-мат. наук; Т.Я. Рябова, канд. физ.-мат. наук; В.Н. Шлапак, канд. физ.-мат. наук

2. УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по управлению качеством продукции и стандартам от 29.12.90 № 3760

3. СРОК ПРОВЕРКИ - 1996 г., периодичность проверок - 5 лет

4. ВВЕДЕН ВПЕРВЫЕ

5. ССЫЛОЧНЫЕ НОРМАТИВНО ТЕХНИЧЕСКИЕ ДОКУМЕНТЫ

Обозначение НТД, на который дана ссылка

Номер пункта

ГОСТ 25645.204-83

1.5

 

СОДЕРЖАНИЕ

1. Общие положения. Входная информация для проведения расчетов. 1

2. Методика расчета эффективности электростатической защиты при радиационном воздействии. 3

Приложение Определение вспомогательных массивов Cм, Bм, Pυ, Mb2, Pb, Mγ и Pк 7