На главную | База 1 | База 2 | База 3

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

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

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

РД 50-25645.210-85

Москва
ИЗДАТЕЛЬСТВО СТАНДАРТОВ
1986

ИСПОЛНИТЕЛИ:

В.Б. Хвостов, канд. физ.-мат. наук; О.Е. Адамов, канд. техн. наук; А.И. Григорьев, д-р мед. наук; А.В. Губарев; В.П. Демин, канд. техн. наук; Е.Е. Ковалев, д-р техн. наук; Г.В. Красильников; Е.Н. Лесновский, канд. техн. наук; А.П. Лобаков; И.И. Максимов; В.Н. Никитинский; А.А. Оглоблин, д-р техн. наук; В.А. Панин; С.М. Перфильева; И.Я. Ремизов, канд. техн. наук; В.А. Сакович, канд. техн. наук; В.М. Сахаров, канд. техн. наук; В.И. Степакин, канд. техн. наук; И.Ф. Усольцев, канд. техн. наук.

УТВЕРЖДЕНЫ И ВВЕДЕНЫ В ДЕЙСТВИЕ постановлением Государственного комитета СССР по стандартам от 10.12.85 г. № 3881

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

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

Безопасность радиационная экипажа
космического аппарата в космическом полете.
Методика расчета энергетической плотности
потоков протонов и нейтронов в космическом аппарате

РД
50-25645.210-85

Введены впервые

Утверждены постановлением Госстандарта от 10 декабря 1985 г. № 3881, срок введения установлен

с 01.07.87

Настоящие методические указания устанавливают алгоритмы расчета энергетической плотности потоков первичных протонов, а также вторичных протонов и нейтронов (нуклонов) внутри космического аппарата (КА) в диапазоне энергий от 1 до 1000 МэВ, учитывающие пространственное распределение вещества в КА, средняя толщина защиты которого, выраженная в массовых единицах длины, не превышает 15 г/см2.

1. ОБЩИЕ ПОЛОЖЕНИЯ

1.1. Расчет энергетической плотности потоков протонов и вторичных нейтронов внутри КА осуществляют в приближении однократного ядерного взаимодействия первичных протонов с веществом защиты КА.

1.2. В расчетах прохождения космических лучей через защиту КА ядерные взаимодействия протонов с энергией E < 20 МэВ не учитывают.

1.3. Учет потерь энергии протонов на ионизацию осуществляют в приближении непрерывного замедления с использованием ионизационных потерь энергии протона S(E) и ионизационных пробегов R(E), определяемых по РД 50-25645.206-84.

1.4. Ослабление флюенса нуклонов с энергией Е0 за счет ядерных взаимодействий в однородном веществе на геометрическом расстоянии t, см, определяют функцией

f(E0,t) = exp[-τ(E0,t)],                                                            (1)

где τ(Е0, t) - оптическое расстояние в веществе, т.е. расстояние, выраженное в долях среднего пробега нуклона до ядерного взаимодействия λ(Е0), определяемого по ГОСТ 25645.211-85.

1.4.1. Оптическое расстояние для нейтронов с энергией Е0 вычисляют по формуле

                                                     (2)

где ρ - плотность вещества, г/см3.

1.4.2. Оптическое расстояние для протонов с начальной энергией Е0 вычисляют по формуле

                       (3)

где                                             (4)

Остаточную энергию протона Е после прохождения в веществе расстояния t вычисляют из соотношения

R(E) = R(E0) - tρ.                                                                  (5)

1.4.3. Для объекта, состоящего из зон с различными веществами, оптическое расстояние вычисляют как сумму оптических расстояний, пройденных нуклоном в каждой из зон.

1.5. Энергетическое распределение плотности потоков протонов и нейтронов  вычисляют в фиксированных точках  внутри КА в энергетических узлах Ej (j = 1,..., J), заданных в диапазоне энергий от 1 до 1000 МэВ в порядке возрастания энергии.

Примечание. Выбор узловых значений Ej и количество узлов J должны обеспечить вычисление интегрального по энергии флюенса нуклонов с погрешностью не более 20 %. Для энергетических узлов при J > 20 следует принять логарифмическую сетку.

2. ИСХОДНЫЕ ДАННЫЕ ДЛЯ РАСЧЕТА ЭНЕРГЕТИЧЕСКОЙ ПЛОТНОСТИ ПОТОКОВ ПРОТОНОВ И НЕЙТРОНОВ

2.1. Геометрию и состав вещества КА задают по ГОСТ 25645.204-83, причем в качестве внешней границы объекта принимают описанную сферу радиуса R0. Начало неподвижной системы координат совмещают с центром сферы (см. рисунок).

2.2. Каждой зоне K космического аппарата ставят в соответствие индекс вещества l(k), заполняющего эту зону. Пустым зонам присваивают значение индекса l = 0.

2.3. Для каждого вещества с индексом l(k), отличным от нуля, определяют следующие характеристики:

ионизационные потери энергии протонов Sl(Е) по РД 50-25645.206-84;

ионизационные пробеги протонов Rl(Е) по РД 50-25645.206-84;

средние пробеги нуклонов до ядерного взаимодействия λпрl(Е) и λнl(E) по ГОСТ 25645.211-85;

двойные дифференциальные распределения вторичных нуклонов  (m = пр, н - протон или нейтрон, соответственно) по ГОСТ 25645.211-85.

Примечание. Средние пробеги нейтронов с энергией E < 20 МэВ определяют в соответствии с рекомендуемым приложением.

Геометрия расчета энергетической плотности потоков
протонов и нейтронов в космическом аппарате:

1 - космический аппарат (КА); 2 - сфера радиуса R0;  - фиксированная точка внутри КА;
 - точка на поверхности сферы;  - точка внутри КА

2.4. Источник первичных протонов задают плотностью тока протонов через поверхность сферы радиуса R0 в каждой точке  (см. рисунок), равной

                                                (6)

где Ф(Е0) - изотропная в 4p плотность потока протонов космических лучей;  = {θ2, φ2} - единичный вектор направления в точке .

Примечание. Угол 0 ≤ θ2p/2 отсчитывается от нормали  в точке .

2.5. Для моделирования траекторий протона методом Монте-Карло функцию  задают в виде

                                 (7)

где  - начальный статистический вес истории;

 - функция плотности вероятности для энергетического распределения протонов;

f22) = 2μ2, где μ2 = cosθ2 - функция плотности вероятности распределения Ω2 по азимутальному углу θ2, 0 ≤ μ2 ≤ 1; f32) = 1/2p - функция плотности вероятности распределения  по полярному углу φ2, 0 ≤ φ2 ≤ 2p; Еmin и Emax - минимальное и максимальное значения энергии в спектре протонов космических лучей, равные соответственно 30 и 1000 МэВ.

3. АЛГОРИТМ РАСЧЕТА ЭНЕРГЕТИЧЕСКОЙ ПЛОТНОСТИ
ПОТОКА ПЕРВИЧНЫХ ПРОТОНОВ

3.1. Декартовы координаты точки вылета протона  (см. рисунок) с поверхности сферы радиусом R0 моделируют для каждой i-ой истории (i = 1, ..., I) по следующей схеме.

3.1.1. Направление  моделируют из изотропного в 4p распределения для :

                                                                 (8)

где η' и η" - случайные числа, равномерно распределенные в интервале (0, 1).

3.1.2. Направляющие косинусы в декартовой системе координат  вычисляют по формулам

                                                              (9)

3.1.3. Декартовы координаты точки  вычисляют по формулам

                                                                         (10)

3.2. Расстояние Т10 между фиксированной точкой  и точкой  вычисляют по формуле

                               (11)

3.3. Направляющие косинусы в неподвижной декартовой системе координат для направления  из точки  на точку  вычисляют по формулам

                                                                    (12)

3.4. Вычисляют косинус угла cosθ*2 между внешней нормалью к сферической поверхности в точке   и направлением  по формуле

cosθ*2 =                               (13)

причем

т.е.

                                                                          (14)

3.5. Значения энергии протона в точке  принимают равным Ej (j = 1, ..., J) и вычисляют значения энергии первичного протона  в точке  по следующему алгоритму.

3.5.1. Направляют луч из точки  в направлении  и вычисляют длины отрезков tn(, ), (n = 1, ..., N0 - 1) между последовательными пересечениями зон КА, а также последовательность номеров пересекаемых зон kn(, ), (n = 1, ..., N0 - 1) по ГОСТ 25645.204-83.

Примечание. Далее в формулах аргументы при tn и kn опущены.

3.5.2. Значения энергии протона E'n (n = 1, ..., N0 - 1) на границах пересекаемых зон kn вычисляют из рекурентного соотношения для пробегов протона:

Rl(kn)(E'n) = Rl(kn)(E'n-1) + ρl(kn)tn,                                             (15)

полагая E'0 = Ej.

3.5.3. Если значение энергии протона E'n при nN0 - 1 превышает Emax =1000 МэВ, т.е. попадание в точку  протонов с энергией Еj из источника невозможно, то для всех j' ≥ j дальнейшее вычисление энергетического распределения не производится.

3.5.4. Значение энергии первичного протона  принимают равным значению энергии протона E'N0-1

 = E'N0-1.                                                                           (16)

3.6. Оптическое расстояние для протона с энергией  между точками  и  в направлении  вычисляют по формуле

                                     (17)

где

Значения E'n - по п. 3.5.2.

3.7. Значение энергетической плотности потока первичных протонов в точке  при энергии Ej вычисляют по формуле

    (18)

3.8. Значение  заносят в сумматор :

                                                (19)

3.9. Для получения статистически обеспеченной оценки , (j = 1, ..., F) вычисления по пп. 3.1 - 3.8 следует выполнить I раз (для сложной геометрии космического аппарата I ≥ 104).

3.10. Значение , нормированное на единичную плотность потока первичных протонов , вычисляют по формуле

                                             (20)

4. АЛГОРИТМ РАСЧЕТА ЭНЕРГЕТИЧЕСКОЙ ПЛОТНОСТИ ПОТОКА ВТОРИЧНЫХ НУКЛОНОВ

4.1. Декартовы координаты вылета протона  - по п. 3.1.

4.2. Протонной истории присваивают статистический вес W = W0, где W0 определяют по п. 2.5.

4.3. Направление вылета протона из источника (см. рисунок) моделируют по следующему алгоритму.

4.3.1. Подвижную систему координат помещают в точку  направив ось z1 по внешней к сферической поверхности нормали , направляющие косинусы которой в неподвижной системе координат nx, ny, nz вычисляют по формулам (14).

4.3.2. Направление вылета протона  = {θ2, φ2} в подвижной системе координат моделируют, используя f22), f32) определяемые по формулам п. 2.5, причем

                                                                                                  (21)

где η', η" - случайные числа, равномерно распределенные в интервале (0, 1).

4.3.3. Значения направляющих косинусов вектора Ω2 в неподвижной системе координатвычисляют по формулам

            (22)

4.4. Значение энергии первичного протона E1* моделируют из распределения f1(E0), определяемого по п. 2.5, полагая

                                                                 (23)

где η - случайное число, равномерно распределенное в интервале (0, 1).

4.5. Оптическое расстояние в КА для протона с энергией E1*, выходящего из точки  в направлении  вычисляют по следующему алгоритму.

4.5.1. Направляют луч из точки  в направлении  и вычисляют длины отрезков  (n = 1,..., N0 - 1) между последовательными пересечениями зон КА, а также последовательность номеров пересекаемых лучом зон  (n = 1, ..., N0 - 1) по ГОСТ 25645.204-83.

4.5.2. Если на пути луча встречаются только пустые зоны, т.е. l(kn) = 0 для всех n = 1, ..., N0 - 1, то история прерывается.

4.5.3. Если на пути луча есть заполненные веществом зоны, то вычисляют значения энергии протона Е'п на границах зон из рекурентного соотношения для пробегов протона:

Rl(kn)(E'n) = Rl(kn)(E'n-1) + ρl(kn)tn, (n = 1, ..., N0 - 1),                (24)

полагая Е'0 =Е1*.

4.5.4. Если Е'N0-1 < 20 МэВ, то в последовательности пересекаемых зон (n = 1,..., N0 - 1) определяют номер n0, при котором Е'n0 < 20 МэВ.

Примечание. Данное условие означает, что протон из-за ионизационных потерь достигает энергии 20 МэВ в зоне с номером kn.

4.5.5. Оптическое расстояние  вычисляют по формуле

(25)

где   - геометрическое расстояние от точки  до пересечения лучом внешней границы объекта.

4.6. Оптическое расстояние, пройденное протоном в КА до ядерного взаимодействия в точке , вычисляют по формуле

                  (26)

где  - геометрическое расстояние между точками  и ; η - случайное число, равномерно распределенное в интервале (0, 1).

4.7. Новый статистический вес протонной истории, компенсирующий смещенность оценки пробега протона до ядерного взаимодействия вычисляют по формуле

                                  (27)

4.8. В последовательности пересекаемых зон (n = 1, ..., N0 - 1) вычисляют номер n1, для которого выполняется неравенство

                        (28)

где  - значения оптических расстояний протона в пересекаемых зонах КА, определенные по п. 4.5.5.

Примечание. Точка  находится в зоне с номером kn1.

4.9. Энергию протона Е2* в точке  вычисляют по следующему алгоритму.

4.9.1. Вычисляют оптическое расстояние, пройденное протоном в зоне с номером kn1:

                      (29)

4.9.2. Значение оптического расстояния для протона с энергией E*2 в веществе зоны n1 вычисляют по формуле

                             (30)

4.9.3. Значение энергии протона E*2 вычисляют путем обратной интерполяции из зависимости Ql(E) для значения Ql, определяемого по формуле (30).

4.10. Вычисляют геометрическое расстояние T12 по следующему алгоритму.

4.10.1. Расстояние t*n1, пройденное протоном в зоне с номером kn1 вычисляют по формуле

                                  (31)

4.10.2. Расстояние T12 вычисляют по формуле

                                                                     (32)

4.11. Декартовы координаты точки ядерного взаимодействия протона  вычисляют по формулам

                                                                 (33)

4.12. Расстояние Т20 между точкой  и точкой  вычисляют по п. 3.2.

4.13. Направление  из точки  на точку  вычисляют по п. 3.3.

4.14. Угол между векторами  и  вычисляют по формуле

                              (34)

4.15. Направляют луч из точки  в направлении  и вычисляют длины отрезков  (n = 1, ..., N0 - 1) между последовательными пересечениями зон КА, а также последовательность номеров пересекаемых лучом зон  (n = 1, ..., N0 - 1) по ГОСТ 25645.204-83.

4.16. В последовательности пересекаемых зон (n = 1, ..., N0 - 1) определяют номер n2, для которого справедливо неравенство

                                                                 (35)

Примечание. Точка  расположена в зоне с номером kn2.

4.17. Длину отрезка t*n2, проходимую нуклоном в зоне kn2, вычисляют из соотношения

                                                                    (36)

4.18. Длину отрезка tn2 полагают равной t*n2.

4.19. Оптическое расстояние между точками  и  для нейтронов с энергией EjE*2 (j = 1, ..., J) вычисляют по формуле

                                      (37)

4.20. Значение энергетической плотности потока вторичных нейтронов в точке  при энергии Ej вычисляют по формуле

                  (38)

4.21. Значение  заносят в сумматор

                                                 (39)

4.22. Энергию вторичного протона , с которой он должен вылететь из точки , чтобы прийти в точку  с энергией Ej (j = 1,..., J) вычисляют по следующему алгоритму.

4.22.1. Определяют энергии протонов на границах пересекаемых зон из рекурентного соотношения для ионизационных пробегов протона

Rl(kn)(E'n) = Rl(kn)(E'n-1) + ρl(kn)tn, (n = 1,..., n2),                       (40)

полагая Е'0 = Ej.

4.22.2. Если значение энергии протона E'n превышает значение E*2, т.е. попадание в точку  протонов с энергией Ej невозможно, то для всех j' ≥ j дальнейшее вычисление энергетического распределения не производится.

4.22.3. Значение энергии вторичного протона , образующегося в протон-ядерном взаимодействии в точке  принимают равным

4.23. Оптическое расстояние для протона с энергией  (j = 1,..., J) между точками  и  вычисляют по формуле

           (41)

4.24. Значение энергетической плотности потока протонов в точке  при энергии Ej вычисляют по формуле

  (42)

4.25. Значение  заносят в сумматор :

                                                (43)

4.26. Для получения статистически обеспеченной оценки  и , (j = 1,..., J) вычисления по п. 4.1 - 4.25 следует выполнить I раз (I ≥ 104).

4.27. Значения  и , (j = 1, ..., J), нормированные на единичную плотность потока первичных протонов, вычисляют по формулам

                                             (44)

ПРИЛОЖЕНИЕ
Рекомендуемое

1. Средние пробеги нейтронов до ядерного взаимодействия с энергией E < 20 МэВ вычисляют по формуле.

где σн(Е) - полное микроскопическое сечение ядерного взаимодействия нейтронов, барн.

2. Полные микроскопические сечения ядерного взаимодействия нейтронов с E < 20 МэВ рекомендуется брать из справочника*.

__________

* Медведев Ю.А., Степанов Б.М., Труханов Г.Я. Ядерно-физические константы взаимодействия нейтронов с элементами, входящими в состав атмосферы и земной коры: Справочник. - М.: Энергоиздат, 1981.