doc_act

РД 50-25645.210-85 МУ. Безопасность радиационная экипажа. Методика расчета энергетической плотности потоков протонов

Реклама

  Скачать документ



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

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

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



Реклама

РД 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 МэВ определяют в соответствии с рекомендуемым приложением.

РД 50-25645.210-85 МУ. Безопасность радиационная экипажа. Методика расчета энергетической плотности потоков протонов

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

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

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

(6)

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

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

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

(7)

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

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

f2(?2) = 2?2, где ?2 = cos?2 - функция плотности вероятности распределения ?2 по азимутальному углу ?2, 0 ? ?2 ? 1; f3(?2) = 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. Направляющие косинусы в неподвижной декартовой системе координат для направления из точки на точку вычисляют по формулам

РД 50-25645.210-85 МУ. Безопасность радиационная экипажа. Методика расчета энергетической плотности потоков протонов (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 при n ? N0 - 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} в подвижной системе координат моделируют, используя f2(?2), f3(?2) определяемые по формулам п. 2.5, причем

(21)

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

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

РД 50-25645.210-85 МУ. Безопасность радиационная экипажа. Методика расчета энергетической плотности потоков протонов (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. Оптическое расстояние между точками и для нейтронов с энергией Ej ? E*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.

СОДЕРЖАНИЕ

1. Общие положения. 1

2. Исходные данные для расчета энергетической плотности потоков протонов и нейтронов. 2

3. Алгоритм расчета энергетической плотности потока первичных протонов. 4

4. Алгоритм расчета энергетической плотности потока вторичных нуклонов. 6

Приложение. 10




Реклама: ;


Самые популярные документы раздела



Рейтинг@Mail.ru Яндекс.Метрика