I. Введение
Методы расчета прохождения ЗЧ КП через материалы основаны на двух подходах: метода статистических испытаний (метод Монте-Карло) и метода численного решения кинетического уравнения Больцмана для плотности потока частиц в рамках различных приближений (детерминистский подход). У обоих методов имеются как преимущества, так и недостатки в сравнении друг с другом.
Недостаток детерминистского подхода заключается в том, что при вычислении плотности потока частиц в среде координаты энергии и углы, которые характеризуют состояние частиц, изменяются на дискретной сетке, а не непрерывно. Поэтому проведение расчетов для гетерогенных структур, корректный учет граничных условий и учет процессов генерации вторичных частиц затруднено. Преимуществом перед методом Монте-Карло (МК) являются существенно более низкие затраты времени счета. Метод МК представляет собой прямой численный эксперимент, основанный на прослеживания большого числа траекторий частиц в рассматриваемом объеме. Метод МК можно применять без ограничений для объектов произвольной геометрии и материального состава. Но для снижения статистических погрешностей необходимо увеличение времени счета. Данный недостаток компенсирует появление быстродействующих ЭВМ. Поэтому классический метод МК и его модификации в настоящее время являются основным инструментом при решении различных задач в области радиационной физики, в том числе и в космических приложениях[8,9].
По результатам расчетов, когда набор статистики необходимый для расчетов закончен, определяются энергетические групповые спектры протонов, электронов и вторичного тормозного излучения (ТИ) за защитой. Данные спектры, если необходимо, будет нетрудно преобразовать в дифференциальные. В практических приложениях, применительно к радиационным условиям на борту КА, поглощенная доза за защитой рассчитывается для кремния, хотя, путем суммирования вкладов от выходящих из барьера частиц, может быть определена для любого материала.
