Настоящие методические указания устанавливают методики расчета спектров линейной энергии для тяжелых заряженных частиц (далее - ТЗЧ) с зарядом от 1 до 32 единиц абсолютной величины заряда электрона и энергией на нуклон от 0,1 до 104 МэВ в тканеэквивалентном веществе при размерах шарового микрообъема от 0,1 до 20 мкм.
Методические указания предназначены для расчетов микродозиметрических характеристик полей ионизирующих излучений (далее - микродозиметрических характеристик), воздействующих на биологические объекты в космических полетах.
Пояснения терминов, применяемых в методических указаниях, приведены в приложении 1.
1.1. Под линейной энергией y, кэВ/мкм, понимают относящуюся к событию поглощения стохастическую величину, равную частному от деления фактически поглощенной в микрообъеме энергии
(1)Для шарового микрообъема диаметром l0:
.Различают частотный и дозовый спектры линейной энергии (далее - СЛЭ). Под частотным СЛЭ f(y), мкм/кэВ, понимают относительное число событий поглощения в микрообъеме, приходящихся на элементарный интервал значений линейной энергии в окрестности y. Под дозовым СЛЭ d(y), мкм/кэВ, понимают относительную долю поглощенной дозы, приходящейся на элементарный интервал значений линейной энергии в окрестности y. Частотный и дозовый СЛЭ нормированы на единицу.
1.2. В качестве основных микродозиметрических характеристик выбирают частотный f(y) и дозовый d(y) СЛЭ, а также определяемые по ним частотное
Другие микродозиметрические характеристики вычисляют по f(y), d(y),
1.3. Методики, представленные в разд. 2 и 3, основаны на предположениях, что при расчете СЛЭ пренебрегают:
- кривизной траекторий ТЗЧ вблизи и внутри микрообъема;
- дополнительными событиями поглощения в микрообъеме и изменением энергий заряженных частиц, обусловленными ядерными взаимодействиями и радиационными потерями вблизи и внутри микрообъема.
1.4. СЛЭ для ТЗЧ, рассчитываемые по методикам разд. 2 и 3, относятся к шаровому микрообъему, выделенному в однородном тканеэквивалентном веществе.
1.5. Методика расчета СЛЭ по методу Монте-Карло (разд. 2) установлена для случая, когда относительная погрешность вычисления
2.1. Настоящая методика основана на моделировании методом Монте-Карло прохождения ТЗЧ, а также возникающих под их действием дельта-электронов вблизи и внутри микрообъема и вычислении поглощенных энергий в нем, соответствующих прохождениям отдельных ТЗЧ. Считают, что центр микрообъема помещен в начало декартовой системы координат {0x1, 0x2, 0x3}, а ТЗЧ движутся в направлении оси 0x3. Изменением энергии ТЗЧ в слое вещества толщиной, достаточной для установления электронного равновесия, пренебрегают.
2.2. В качестве исходных данных для расчета частотного f(y) и дозового d(y) СЛЭ выбирают:
- энергию E, МэВ, атомный номер Z и массовое число A ТЗЧ;
- диаметр l0, мкм, шарового микрообъема;
- значения линейной энергии yj, кэВ/мкм, определяющие интервалы
, для усреднения значений СЛЭ. Полагают y0 = 0.Примечание. Следует применять логарифмическую сетку, удовлетворяющую требованию, чтобы дополнительная погрешность вычисления дозового среднего значения линейной энергии по формуле
(4)связанная с выбором значений линейной энергии yj, не превышала 2%.
- максимальную энергию Tmax, кэВ, дельта-электрона по формулам:
(5) (6)EA = 1,066·10-3 E/A; (7)
- линейную передачу энергии
(8)где dE/dx - ионизационные потери ТЗЧ в тканеэквивалентном веществе, кэВ/мкм;
T0 - минимальная энергия моделируемых дельта-электронов, определяемая из соотношения: T0 = max{0,1; 0,007l0}, кэВ;
I - средний эффективный потенциал ионизации вещества по РД 50-25645.206, кэВ;
(9)- граничную энергию Tm, кэВ, дельта-электронов, до которой учитывается пространственная корреляция траекторий дельта-электронов с траекторией ТЗЧ, по формуле
(11)где Rp(Tmax) - практический пробег, мкм, электрона с энергией Tmax;
T(10l0) - кинетическая энергия электрона, практический пробег которого равен 10l0, кэВ;
- долю ионизационных потерь ТЗЧ,
(12)- радиус am, мкм, сечения области моделирования плоскостью, перпендикулярной траектории ТЗЧ, по формуле
(13)где Rp(Tm) - практический пробег электрона с энергией Tm.
Примечания:
1. Значения dE/dx вычисляют при E/A >= 2 МэВ/нуклон по РД 50-25645.206, а при E/A < 2 МэВ/нуклон - по данным табл. 1 и формуле (124) приложения 3. Практические пробеги электронов определяют по данным табл. 2 приложения 4.
2. Область моделирования - микрообъем и прилегающие к нему слои вещества, в пределах которых производится моделирование прохождения ТЗЧ и возникающих дельта-электронов с учетом пространственной корреляции их траекторий.
2.4. Область моделирования определяют неравенствами
|x1| <= am; (14)
|x2| <= am; (15)
(16)2.5. Алгоритм расчета СЛЭ для ТЗЧ, пересекающих область моделирования, состоит в следующем.
2.5.1. По очередному случайному числу
(17)x2,0 = 0; (18)
(19)2.5.2. Вычисляют параметры n-й истории:
- статистический вес координаты x1,0 по формуле
WR = x1,0; (20)
- координату x3,вых, соответствующую точке выхода ТЗЧ из области моделирования, по формуле
(21)- угол
и пересекающимися по линии, совпадающей с траекторией ТЗЧ, по формуле (22)где l - свободный параметр (l0 < l <= 10l0), мкм, выбираемый из условия, чтобы вероятностью попадания в микрообъем дельта-электронов с энергией менее Tm, вылетающих из точки возникновения с x1,0 >= l в противоположном по отношению к нему направлении, можно было бы пренебречь.
2.5.3. Среднюю энергию
(23)Истинное значение энергии
(24)где
(25)Значения
При
полагают (26)где m - целое число, удовлетворяющее, при очередном
![]() При
полагают![]() где
При
полагают .2.5.5. По очередному
(27)где
Значения
(28)где
(29)а Tr - энергия дельта-электрона, имеющего практический пробег
![]() 2.5.6. При x3,i >= x3,вых полагают
и переходят к вычислениям п. 2.5.11.При x3,i < x3,вых по очередным
(30)значения величин
(31) (32)и статистический вес WT,i энергии Ti дельта-электрона
(33)2.5.7. При
к текущему значению При
переходят к вычислениям п. 2.6.2.5.8. Вычисляют значения D по формуле
(34)При D <= 0 повторяют вычисления с п. 2.5.5.
2.5.9. Расстояние S0 от точки испускания i-го дельта-электрона до ближайшей точки пересечения луча в направлении движения дельта-электрона с поверхностью
вычисляют по формуле (35)Далее переходят к п. 2.6.
2.5.10. К текущему значению
Значение линейной энергии yn и статистический вес Wn для n-й истории вычисляют по формулам:
(36)где i0 - номер дельта-электрона последнего перед выходом ТЗЧ из области моделирования.
Находят наименьшее значение индекса jm, при котором yn <= yj, где yj - выбранные узлы разбиения шкалы линейной энергии (j = 1, 2, ..., J). (Далее j = jm).
2.5.12. В сумматоры y*F, y*D, f*(yj) и d*(yj) заносят вклады от n-й истории, равные, соответственно, ynWn,
(38) (39) (40) (41)где
(42)2.5.13. В сумматор числа событий поглощения N заносят единицу. При N не кратном 20 повторяют вычисления по п. 2.5.1.
2.5.14. Вычисляют и запоминают оценки частотного и дозового средних значений линейной энергии для очередной серии из 20 событий поглощения:
(43) (44)где m - индекс, означающий, что помеченная им величина относится к m-й серии, а также текущие значения
(45) (46)2.5.15. При выполнении условия (для m > 10)
моделирование траекторий ТЗЧ прекращают, переходя к п. 2.5.16, если иначе, то продолжают расчет, начиная с п. 2.5.1.
2.5.16. Рассчитывают окончательные оценки частотного
(48) (49) j = 1, 2,..., J, (50) j = 1, 2,..., J, (51)где
(52)Далее переходят к вычислению п. 2.7.
2.6. Траектории дельта-электронов (далее - электронов) моделируют с учетом их кривизны и возможности рождения вторичных, третичных и т.д. поколений электронов. Процедура вычисления энергии
2.6.1. Присваивают исходные значения сумматору поглощенных энергий
- направляющим косинусам единичного вектора
, задающего начальное направление движения электрона в системе координат {0x1, 0x2, 0x3}, (53) (54) (55)- координатам радиуса-вектора
, задающего точку начала моделируемой траектории (j = 1, 2, 3): (56)- энергии электрона
(57)где Rp - практический пробег электрона с энергией Ti;
2.6.2. Ограниченные линейные передачи энергии
и упругих взаимодействий, полное 2.6.3. Длину пути
(58)где
Координаты радиуса-вектора
точки взаимодействия вычисляют по формуле (j = 1, 2, 3): (59)2.6.4. При
2.6.5. Длине S части отрезка
, принадлежащей микрообъему, присваивают в зависимости от знака параметра (60)следующие значения:
при D <= 0 полагают S = 0;
при D > 0 полагают
(61)и
(62)если
, то взаимодействие неупругое. Для его моделирования переходят к п. 2.6.10.2.6.7. Энергию электрона в точке
2.6.8. Элемент, на котором произошло упругое рассеяние, определяют по очередному
(63)2.6.9. Величины
и , определяющие направление вылета электрона из точки упругого взаимодействия, вычисляют по формулам (64) (65)где
F(
Далее выполняют вычисления, начиная с п. 2.6.14.
2.6.10. Потерю энергии
в точке неупругого взаимодействия вычисляют по очередному Энергии электронов, покидающих точку
(66) (67) (68) (69) (70) (71)где
2.6.11. При
к текущему значению (72)2.6.12. При
запоминают радиус-вектор (74) (75)и переходят к п. 2.6.14.
(77) (78)где
или
моделирование траектории электрона данного поколения прекращают. При
к текущему значению В случае невыполнения условий (79) и (80) переходят к следующему пункту.
2.6.15. Направляющие косинусы вектора
, задающего направление движения рассматриваемого электрона из точки (81) (82) (83)2.6.16. Координатам вектора
(84) (85)и повторяют расчеты начиная с п. 2.6.2.
2.7. При
расчет завершают, полагая искомые вычисляют их по формулам: (86) (87) (88) (89)где помеченные индексом (e) величины относятся к событиям поглощения, формируемым дельта-электронами с энергиями свыше Tm. Эти величины, одинаковые для всех ТЗЧ одной скорости, но разных зарядов, рассчитывают согласно п. 2.8.
2.8. Методика расчета
, концентричной рассматриваемому шаровому микрообъему диаметром l0, и моделировании прохождения электронов внутри этой сферы методом Монте-Карло состоит в следующем.2.8.1. Нормированный на единицу интегральный спектр флюенса Ф(T) у поверхности сферы диаметром l вычисляют по формулам:
(90) (91)где Le(T) - линейная передача энергии, кэВ/мкм, для электрона с энергией T в тканеэквивалентном веществе, определяемая по данным табл. 4 приложения 4.
2.8.2. Для точки вылета электрона в n-й истории принимают x1,0 = 0, x2,0 = 0, x3,i = -l/2,
(92) (93)где
- функция, обратная Ф(T);2.8.3. Выполняют вычисления согласно п. 2.6 с той лишь разницей, что вместо предусмотренного в пп. 2.6.4, 2.6.6 и 2.6.14 перехода в п. 2.5.10, переходят в п. 2.8.4.
2.8.4. В случае
2.8.5. Окончательные оценки искомых величин для событий поглощения, формируемых дельта-электронами с энергией более Tm, рассчитывают по формулам:
(94) (95) j = 1, 2, ..., J; (96) j = 1, 2, ..., J, (97)где
(98)3.1. Методика основана на предположении, что разбросом энергетических потерь ТЗЧ на отрезках траекторий внутри микрообъема можно пренебречь и что вся потерянная ТЗЧ энергия поглощается в точках их взаимодействий с веществом.
3.2. В качестве исходных данных для расчетов частотного f(y) и дозового d(y) СЛЭ выбирают:
- энергию E, МэВ, атомный номер Z и атомную массу A ТЗЧ;
- диаметр l0, мкм, микрообъема.
3.3. Линейные передачи энергии L = -dE/dx, кэВ/мкм, и ионизационные пробеги R, мкм, ТЗЧ в тканеэквивалентном веществе, используемые в расчетах, вычисляют при E/A >= 2 МэВ/нуклон по РД 50-25645.206, а при меньших энергиях - по формулам и данным приложения 3.
3.4. В случае, когда соблюдается условие
применяют следующие формулы для СЛЭ, частотного
(100) (101) (102) (103)3.5. В случае, когда условие (99) не соблюдается, частотный f(y) и дозовый d(y) СЛЭ представляют в виде:
(104) (105)где c - постоянная величина, определяемая из условия нормировки на единицу
(106)Физический смысл и формулы для расчета каждого слагаемого при равномерно распределенных в среде источниках ТЗЧ приведены в пп. 3.5.1 - 3.5.4.
3.5.1. Слагаемое N1(y) определяет вклад в СЛЭ от ТЗЧ, треки которых полностью принадлежат микрообъему. Значения N1(y) рассчитывают по формуле
(107)где
дельта-функция, а (108)Здесь и далее Ey = 2·10-3l0y/3.
3.5.2. Слагаемое N2(y) определяет вклад в СЛЭ от ТЗЧ, треки которых начинаются внутри микрообъема, но заканчиваются вне его. Значения N2(y) рассчитывают по формуле
(109)где x = R(E) - R(E - Ey),
(110)3.5.3. Слагаемое N3(y) определяет вклад в СЛЭ от ТЗЧ, треки которых начинаются вне микрообъема, но заканчиваются внутри него. Значения N3(y) рассчитывают по формуле
(111)где
.3.5.4. Слагаемое N4(y) определяет вклад в СЛЭ от ТЗЧ, пронизывающих микрообъем. Значения N4(y) рассчитывают по формуле
(112)где x = R(E) - R(E - Ey),
(113)3.5.5. Частотное
Справочное
К ТЕРМИНАМ, ПРИМЕНЯЕМЫМ В МЕТОДИЧЕСКИХ УКАЗАНИЯХ
Справочное
СО СПЕКТРАМИ ЛИНЕЙНОЙ ЭНЕРГИИ (СЛЭ)
В микродозиметрии и ее приложениях, помимо линейной энергии y, частотного f(y) и дозового d(y) СЛЭ, частотного
Под удельной энергией z, Гр, понимают стохастическую величину, равную частному от деления фактически поглощенной в микрообъеме энергии
(114)При y, кэВ/мкм, z, Гр, и диаметре шарового микрообъема l0, мкм, справедливы следующие соотношения:
(115) (116) (117) (118) (119)Для вычисления плотности распределения удельной энергии при заданной поглощенной дозе следует использовать формулу
(120)где f(k)(z) - k-кратная свертка от f1(z), определяемая с помощью рекуррентного соотношения
(121)При достаточно больших
и достаточно малых поглощенных дозах fD(z) вычисляют по формулам: (122)и
(123)соответственно, где
Справочное
(ТЗЧ) С ЭНЕРГИЯМИ НА УКЛОН МЕНЕЕ 2 МэВ
Ионизационные потери dE/dx, кэВ/мкм, в тканеэквивалентном веществе для ТЗЧ при энергиях на нуклон менее 2 МэВ следует рассчитывать по ионизационным потерям (dE/dx)p для протонов, представленным в таблице, согласно формуле:
где Zэфф(Z) - эффективный заряд ТЗЧ с зарядом ядра Z, определяемый формулой (10) разд. 2.
Ионизационные пробеги R(E), мкм, ТЗЧ, включая протоны, следует определять по ее ионизационным потерям согласно формуле
(125)где E - энергия ТЗЧ, МэВ.
Таблица 1
в тканеэквивалентном веществе
Примечание. Для получения dE/dx при промежуточных значениях E следует применять линейную интерполяцию в двойном логарифмическом масштабе, а при E < 0,001 МэВ/нуклон - линейную экстраполяцию в обычном масштабе.
Справочное
ЭЛЕКТРОНОВ В ТКАНЕЭКВИВАЛЕНТНОМ ВЕЩЕСТВЕ
В табл. 2 представлены значения практического пробега Rp электрона в тканеэквивалентном веществе при энергиях T в диапазоне от 0,1 до 200 кэВ. Значения Rp для энергий этого диапазона, не представленных в табл. 2, следует вычислять методом линейной интерполяции в двойном логарифмическом масштабе.
Таблица 2
в тканеэквивалентном веществе
В табл. 3 представлены значения функции
, определяющей зависимость энергии электрона от глубины его проникновения S, мкм, в тканеэквивалентное вещество (126)где T(0) - начальная энергия электрона, кэВ;
Rp - практический пробег, мкм, электрона с энергией T(0).
Значения Rp при промежуточных T следует вычислять методом линейной интерполяции в обычном масштабе.
Таблица 3
В табл. 4 и 5 представлены значения ионизационных пробегов Re, полных Le и ограниченных
и полного упругих взаимодействий Значения представленных в табл. 4 и 5 величин при промежуточных T следует получать методом линейной интерполяции в двойном логарифмическом масштабе.
Таблица 4
для электронов в тканеэквивалентном веществе
Таблица 5
для электронов в тканеэквивалентном веществе
Параметр экранирования ядра атомными электронами
рассчитывают по формуле (127)где
Zk - атомный номер k-го элемента;
Нумерация и значения параметров Zk,
Таблица 6
тканеэквивалентного вещества
Справочное
В НЕУПРУГИХ ВЗАИМОДЕЙСТВИЯХ
Потери энергии электронов в неупругих взаимодействиях с атомными электронами в тканеэквивалентном веществе, превышающие заданный порог
(128)c1 = 1 - 2c0; (129)
(130)a1 = c1/c2 + c2ln(2c0); (131)
a2 = c1/(1 - c0) - c2ln[2(1 - c0)]; (132)
(133) (134)2. Для очередного случайного числа
(135) (136)4. Для очередного
(137) (138)Если оно выполнено, то переходят в п. 6, в противном случае вычисляют
(139)и переходят в п. 8.
(140)7. Для очередного
(141)Если оно выполнено, то возвращаются в п. 6, в противном случае переходят к следующему пункту.
(142)Вернуться в "Каталог нормативных документов"
Источник информации: https://internet-law.ru/documents/prod/metodika/7/rd_35579.html
На правах рекламы:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||