Труды КНЦ вып.5. ГЕЛИОГЕОФИЗИКА. 8/2019(10)

В представленной работе предлагается использовать численные методы Монте-Карло и современные модели взаимодействия частиц, реализованные в пакете для разработки программ GEANT4 [8] для решения задачи определения параметров потоков вторичных частиц и количественной оценки скорости ионизации, вызванной прохождением протонов ГКЛ через атмосферу Земли во время минимума солнечной активности. 1. Параметризация модели и входные данные для источника первичных частиц Для параметризации атмосферы Земли при расчете прохождения протонов ГКЛ используется эмпирическая модель NRLMSISE-00 [9]. Для каждого значения широты и долготы рассчитываются параметры температуры, плотности и процентного соотношения химических элементов в интервале высот от 0 км до 80 км с шагом 0,1 км. Рис. 2. Блок-схема алгоритма расчета значений высоты и плотности для одного слоя атмосферы в модели прохождения протонов КЛ. Этот процесс повторяется N-1 раз, a H a N расчет производится посредством простого усреднения оставшихся значений в таблице данных NRLMSISE. Следует заметить, что показатели процентного содержания веществ (водород, гелий, кислород, азот и аргон) вычисляются также, как и в случае с плотностью, по уже известным данным Fig. 2. The flowchart of the algorithm for calculating altitude and density values for one atmosphere layer in the model of the CR proton transport. This process is repeated N-1 times, and N is calculated by simply averaging the remaining values in the NRLMSISE data table. It should be noted that the percentages of substances (hydrogen, helium, oxygen, nitrogen and argon) are calculated as well as in the case of density, according to already known data 242