Комплексные исследования полярной ионосферы / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Ленинград : Наука, 1970. – 184 с.

углов рассеяния велика, и обратный поток фотонов будет давать значительный вклад при регистрации ненаправленным де­ тектором. Точный аналитический расчет потока и энергетического спектра фотонов и их распределения по глубине проникновения затруд­ нителен. Поэтому задача решается путем моделирования данного процесса на ЦВМ методом Монте-Карло. Сущность расчета сводится к тому, что, задаваясь случайными числами в соответствии с теми законами распределения, которым подчиняются те или иные вычисляемые параметры, можно просле­ дить за движением фотона до его поглощения или до выхода из среды в обратном направлении (альбедо). Основным элементом модели является определение длины свободного пробега фотона до взаимодействия, типа взаимодействия (фотоэффект или эффект Комптона), новых энергий и направления движения фотона в случае комптоновского рассеяния. Длина пробега фотона в воздухе определяется по формуле г = ~ г т 1пт . (4) с(а) £ * ' ' где I — случайное число, распределенное по равномерному закону в интервале (0.1); о (а) — полный массовый коэффициент погло­ щения; I — длина свободного пробега. Ряд значений -g задан таблично, а промежуточные значе­ ния можно вычислить с помощью интерполяционной формулы 1 1 I а ~~ а п+і Г 1_____________ 1 1 ( 2 \ ° и к+і) а « - а п +1 L °п ю °я+і к+і) J ’ где ая< а < а я+1. Глубина проникновения фотона после і-го акта рассеяния определяется по формуле * h i = 2 ^ COS Ѳ., J=0 J J где lj — длина свободного пробега фотона после /-го акта рассея­ ния; Ѳ^. — угол между вертикалью и направлением движения фотона. После расчета пробега, если фотон не вышел из среды (h > 0), определяется тип взаимодействия. Коэффициент поглощения о складывается из двух величин: 0 = °ф+ ак. (4) где — коэффициент поглощения для фотоэффекта; ож— коэф­ фициент поглощения для комптоновского рассеяния. JG0