WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

 

На правах рукописи

Косарев Николай Иванович

РАДИАЦИОННЫЕ РЕЗОНАНСНЫЕ ПРОЦЕССЫ В ОПТИЧЕСКИ ПЛОТНЫХ СРЕДАХ

Специальность 01.04.05 – Оптика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

доктора физико-математических наук

Красноярск - 2010

Работа выполнена в Институте вычислительного моделирования СО РАН и Федеральном государственном образовательном учреждении высшего профессионального образования «Сибирский юридический институт Министерства внутренних дел Российской Федерации (г.Красноярск)»

Научный консультант        доктор физико-математических наук

       Шапарев Николай Якимович

Официальные оппоненты:        доктор физико-математических наук

       Карпов Сергей Васильевич;

       доктор физико-математических наук

       Слабко Виталий Васильевич;

       доктор физико-математических наук

       Солдатов Анатолий Николаевич

Ведущая организация:        Институт оптики атмосферы СО РАН, г. Томск

Защита состоится «16» апреля 2010 г. в 14 часов на заседании диссертационного совета Д 003.055.01 при Институте физике им. Л.В. Киренского СО РАН по адресу: 660036, г.Красноярск 36, Академгородок, 50, стр. 38. Институт физики им. Л.В. Киренского СО РАН.

С диссертацией можно ознакомиться в библиотеке Института физики им. Л.В. Киренского СО РАН.

Автореферат разослан        «__» __________ 2010 г.

Ученый секретарь

диссертационного совета,

доктор физико-математических наук        Втюрин А.Н.

Общая характеристика работы

В работе обсуждаются результаты математического моделирования процессов воздействия резонансного излучения на динамические оптически плотные газовые и плазменные среды. Построенные модели и алгоритмы численного решения систем интегродифференциальных уравнений применены для исследования задач: лазерной резонансной флуоресценции паров металлов и ультрахолодной плазмы; резонансного оптического разряда; ионизации и свечения искусственных бариевых облаков (БО) под действием солнечного света.

Актуальность исследований. Явление резонансной флуоресценции широко используется в спектроскопии для диагностики светящегося газа и плазмы. По спектральным данным, полученным от оптически тонких сред, однозначно определяют параметры излучающего газа. В плотных газовых и плазменных средах их физическое состояние завуалировано эффектами многократного переизлучения (переноса) фотонов атомами и ионами. Для описания процесса блуждания фотонов в плотных средах Холстейн1 предложил ввести такую оптическую характеристику, как эффективное время высвечивания. Он и независимо от него Биберман2

также выдвинули предположение о том, что в процессе рассеяния фотонов атомами реализуется модель полного перераспределения по частотам (ППЧ). В итоге, задача о высвечивании оптически плотной среды в приближении Бибермана-Холстейна сводилась к интегральному в стационарном, а к интегродифференциальному в нестационарном случае уравнениям, которые позволяли находить концентрации частиц в возбужденном состоянии для двухуровневой модели атома. Одновременно с выводом данных уравнений авторы предложили и приближенные методы их решения, которые до сих пор часто используются для интерпретации спектроскопических данных по свечению лабораторной и космической плазмы и газов.

Дальнейшее развитие теории переноса излучения было стимулировано задачами астрофизики, связанными с исследованием спектров свечения различных космических объектов, звездных и планетных атмосфер. Их решение позволило сформировать и построить классическую теорию радиационного переноса. Пионерскими в данной области следует считать монографии таких авторов, как Соболев В.В., Иванов В.В., Михалас Д., Нагирнер Д.И., и др.

В настоящее время большой класс задач спектроскопии касается воздействия лазерного излучения на плотные газовые среды и ультрахолодную плазму (УП). Следует отметить, что впервые идею охлаждения плазмы резонансным излучением и ее удержание в магнитооптической ловушке предложили авторы3. Сегодня УП является объектом повышенного интереса у теоретиков и экспериментаторов, что вызвано их желанием исследования физических условий процесса кристаллизации плазмы.

Оптический канал получения информации нашел большое применение и в связи с развитием активных ракетных экспериментов в верхней атмосфере Земли и открытом Космосе. Для исследования физических свойств космического пространства с начала 80 – х годов на высотах от 160 км и выше инжектируют искусственные облака (ИО). Находясь в поле солнечного света ИО рассеивают излучение в резонансных атомных и ионных линиях. Спектроскопические данные фиксируются наземной теле- фотоаппаратурой, на основе чего определяется массоперенос, направления и величины высокоширотных ветров, геомагнитных и электрических полей.

Особенности моделей, описывающих взаимодействие излучения с веществом в перечисленных выше постановках задач таковы, что возникает необходимость учета целого ряда сопутствующих факторов: трехмерную геометрию газовой среды; большое количество линий атомов и ионов; различные механизмы формирования спектральной линии и модели перераспределения фотонов по частотам; макроскопический разлет вещества и неоднородное пространственное распределение его параметров (концентрацию, температуру и вектор направления скорости разлета частиц); локальное возбуждение малого объема среды и наличие возбужденных атомов/ионов в среде, формирующих функцию источников; динамику параметров среды и воздействующего излучения. Перенос излучения в таких динамических многоуровневых системах изучен крайне слабо4. Таким образом, решение описанного круга задач является возможным только с привлечением численных алгоритмов и развитием методов математического моделирования в теории переноса излучения.

Цель работы состояла в исследовании задач резонансного воздействия солнечного и лазерного излучения на оптически плотные газы и плазму, на основе построенных моделей, численных алгоритмов и пакетов прикладных программ.

В программу исследований входило:

  • построение физико-математических моделей, описывающих взаимодействие солнечного и лазерного излучения с оптически плотными средами;
  • разработка численных алгоритмов для решения полученных систем интегродифференциальных уравнений, учитывающих многоуровневость атомов (ионов), конечный замкнутый объем газовой среды, эффекты полного и частичного перераспределения фотонов по частотам, макроскопический разлет вещества, неоднородное пространственное распределение плотности и температуры вещества, совместное действие различных механизмов уширения контура спектральной линии;
  • проведение численных расчетов по лазерной резонансной флуоресценции и ионизации паров металлов; поглощению лазерного излучения, рассеянию и испусканию света ультрахолодной плазмой; ионизации и свечению БО широкополосным солнечным светом;
  • анализ расчетных данных и сравнение с натурными и лабораторными экспериментами.

Сложность решения задач о переносе резонансного излучения в многоуровневых динамических средах вызвана также тем, что реализация численных алгоритмов на ЭВМ значительно усложняется необходимостью учета сферической или цилиндрической геометрии газа, локального фотовозбуждения малого объема среды лазерным пучком, частичного перераспределения по частотам, макроскопического движения вещества и др. Одновременный учет хотя бы некоторых из этих физических факторов в одной численной модели сопряжен с серьезными трудностями в вычислительном плане и требует привлечения и дальнейшее развитие численных методов в теории переноса излучения и в вычислительной математике.

Выполненные исследования вносят значительный вклад в развитие научного направления «Перенос резонансного излучения в оптически плотных, динамических, газовых и плазменных средах». Результаты численного моделирования задач: о фотоионизации и свечении БО; лазерной резонансной ионизации и флуоресценции паров щелочных металлов; поглощении, испускании и рассеянии света в расширяющейся ультрахолодной плазме, являются приоритетными в области математического моделирования в лазерной физике, в оптике и спектроскопии, а также в физике космической плазмы и газа.

Методы исследований. Для выполнения поставленных целей автором использовались известные методики расчета сечений и скоростей радиационно-столкновительных процессов, физические справочники оптических характеристик химических элементов, модели и методы физики плазмы и теории переноса излучения, численные методы решения систем обыкновенных дифференциальных уравнений, дискретно-разностные схемы и методы, квадратурные формулы при вычислении соответствующих интегралов.

Достоверность и обоснованность подтверждается качественным и количественным совпадением результатов с экспериментальными данными натурных ионосферных и лабораторных исследований. Численные модели основаны на уравнениях баланса населенностей, скоростные коэффициенты которых рассчитывались по известным методикам, основанных на квантовомеханических вычислениях. Построенные численные алгоритмы протестированы на аналитических решениях.

На защиту выносятся:

1. Математическая модель процесса распространения солнечного и лазерного излучения в многоуровневых средах сложных геометрических конфигураций.

2. Численные алгоритмы решения систем интегродифференциальных уравнений для сред в форме плоскопараллельного слоя, сферы и цилиндра при полном и частичном перераспределении по частотам, макроскопическом разлете вещества и совместном влиянии различных типов уширения спектральной линии.

3. Результаты численного исследования лазерной флуоресценции и эффективного времени высвечивания паров натрия и лития для различных функций перераспределения фотонов по частотам и геометрий газовых сред, в условии фотовозбуждения всей среды и при локальном облучении малого её объема лазерным пучком.

4. Результаты численных экспериментов по формированию резонансного оптического разряда, расчет уширения ионизованного канала в условии переноса излучения из зоны насыщения на периферию среды и скорости волны ионизации, метод определения характерного времени ионизации паров.

5. Модель, алгоритмы и результаты моделирования процесса испускания света, поглощения и рассеяние лазерного излучения в ультрахолодной плазме, метод определения скорости разлета плазмы.

6. Численное моделирование фотоионизации, фотовозбуждения и свечения бариевого облака под действием солнечного света, расчет спектрально-яркостных характеристик свечения и характерного времени ионизации бария.

Научная новизна. Исследованы динамика переноса излучения в многоуровневых средах, которые имеют форму плоскопараллельного слоя, цилиндра и сферы, характеристики радиационных полей и распределение плотности частиц по объему. Проведен анализ роли возбужденных атомов (ионов), формирующих функцию источников среды, в задачах: о лазерно-индуцированной флуоресценции паров щелочных металлов; резонансной лазерной ионизации паров натрия; испускании света, поглощении и рассеянии лазерного излучения ультрахолодной плазмой; фотоионизации и свечении искусственных бариевых облаков под действием солнечного света.

К числу приоритетных исследований автор относит следующие:

  • физико-математическую модель и алгоритмы численного решения задачи о переносе широкополосного излучения в оптически плотных многоуровневых средах сложных геометрических конфигураций, которые учитывают полное и частичное перераспределение по частотам, радиальный разлет вещества. Благодаря использованию осевой симметрии поставленных задач удалость более чем в 10 раз сократить размерность системы интегродифференциальных уравнений;
  • радиационно-столкновительную модель, алгоритмы и результаты численного исследования процесса формирования резонансного лазерного разряда в парах натрия;
  • исследование влияния расширения зоны насыщения на «аномальную» флуоресценцию паров щелочных металлов, обсуждаемую в некоторых экспериментах, а также роли радиационного переноса в расширении плазменного канала, созданного лазерным излучением;
  • исследование эффективного времени жизни возбужденных атомов натрия и лития, описываемые различными моделями функций перераспределения фотонов по частотам и геометриями газовой среды, для объяснения и интерпретации экспериментальных данных;
  • модель, численные алгоритмы и результаты моделирования флуоресценции УП в поле лазерного излучения, методику оценки скорости радиального разлета кальциевой плазмы, основанную на измерении проходящего лазерного излучения в крыле линии;
  • результаты моделирования фотоионизации и расчет спектрально-яркостных характеристик свечения бариевого облака сферической формы при его фотоионизации широкополосным солнечным светом, которые позволили объяснить сложную цветовую окраску облаков, наблюдаемую в натурных ионосферных экспериментах.

Практическая значимость. Проведение натурных ионосферных экспериментов в верхней атмосфере Земли и открытом Космосе путем инжекции искусственных светящихся облаков определяет необходимость развития методов их дистанционной оптической диагностики. Построенные модели и алгоритмы позволяют рассчитывать спектрально-яркостных характеристики свечения БО, на основании которых можно проводить интерпретацию экспериментальных данных, планировать будущие эксперименты с одновременным прогнозом ожидаемых результатов.

Расчетно-теоретическая модель позволяет по интенсивности рассеянного излучения исследовать эффективное время жизни возбужденных атомов в широком диапазоне физических условий, по сравнению с асимптотической моделью Бибермана – Холстейна. Полученная зависимость от оптической толщины фактора пленения паров натрия и лития позволила указать на то, что немонотонное поведение кривых в эксперименте5

не следует приписывать эффектам частичного перераспределения по частотам.

По динамике уменьшения во времени интенсивности свечения паров можно рассчитывать характерное время ионизации атомов под действием лазерного излучения.

Исследование роли переноса вторичной радиации вне облучаемую зону среды указало на то, что в ряде экспериментов по резонансному лазерному воздействию на пары металлов, «аномальная» флуоресценция может быть объяснена только расширением области сильного насыщения среды.

Методика расчета скорости разлета ультрахолодной плазмы по измерению интенсивности проходящего лазерного излучения в крыле линии может быть использована на практике.

Численная модель скоростных уравнений баланса населенностей и переноса излучения дополняет и развивает асимптотическую приближенную модель переноса излучения Бибермана – Холстейна. Она, при одновременном учете целого ряда сопутствующих физических условий, привносит значительный фундаментальный вклад в классическую теорию переноса излучения в атмосферах звезд и планет.

Разработанные алгоритмы и методы математического моделирования процесса распространения резонансного излучения в оптически плотных средах сложных геометрических конфигураций могут быть использованы для решения целого класса задач радиационно-ударной кинетики космической и лабораторной плазмы и газов.

Реализация результатов. Расчетные данные по свечению бариевых облаков при их фотоионизации солнечным излучением переданы Институту прикладной геофизики им. академика Е.К. Федорова в форме научно-исследовательского отчета по теме «Моделирование газодинамических, радиационных и электродинамических явлений при образовании ионных облаков в ионосфере».

Апробация работы. Основные результаты и положения диссертационной работы докладывались и обсуждались на конференциях: III рабочее совещание по моделированию космических явлений в лабораторной плазме (г. Новосибирск, 1990); II Всесоюзном симпозиуме по радиационной плазмодинамике (г. Москва, 1991); XI Всесоюзном и XII Межреспубликанском симпозиумах по распространению лазерного излучения в атмосфере и водных средах (г. Томск, 1991 и 1993); I, II, III и IV Межреспубликанских симпозиумах «Оптика атмосферы и океана» (г. Томск, 1994; 1995; 1996 и 1997); II and VI Chine-Russian Symposium on Laser Physics and Laser Technology (Harbin, China, 1995; China, 2002); III, V, VII and IX Russian-Chinese Symposium on Laser Physics and Laser Technology (Krasnoyarsk, Russia, 1996; Tomsk, 2000; Tomsk 2004; Tomsk 2008); Международных конференциях «Математические модели и методы их исследования» (Красноярск, 1997; 1999; 2001); III и IV Всероссийских семинарах «Моделирование неравновесных систем» (г. Красноярск, 2000; 2001); VII International Conferences «Atomic and Molecular Pulsed Lasers» (Tomsk, Russia, 2005).

Материалы диссертации также докладывались на научных семинарах организаций: ИВМ СО РАН (г. Красноярск), ИФ СО РАН им. Л.В. Киренского (г. Красноярск), КрасГУ (г. Красноярск), ИОА СО РАН (г. Томск), ТГУ (г. Томск).

Личный вклад автора. Первоначальные результаты по свечению бариевых облаков были получены совместно со Шкедовым И.М., а по резонансному оптическому разряду и взаимодействию лазерного излучения с УП с Шапаревым Н.Я. Определяющая часть исследований, составляющих содержание работы, получена лично автором.

Структура и объем диссертации. Работа состоит из введения, пяти глав, заключения и библиографии. Содержание диссертационной работы изложено на 244 страницах, включая 90 рисунков. Список используемых источников содержит 175 наименований.

Краткое содержание работы

Во введении проведен краткий обзор теоретического и экспериментального материала в области исследований, касающейся переноса резонансного излучения в газах и плазме. Обоснована актуальность, научная новизна и практическая значимость работы. Описан личный вклад автора и апробация результатов исследований.

Первая глава посвящена изложению численных методов и алгоритмов, которые лежат в основе всей работы. В первом параграфе выписаны скоростные уравнения баланса населенностей в условии полного перераспределения фотонов по частотам. Для двухуровневого атома совместно с начальными условиями они имеют следующий вид

       (1)

       (2)

       (3)

где и – населенности основного и возбужденного состояний, и A21 – коэффициенты Эйнштейна для вынужденного фотовозбуждения, фототушения и спонтанного распада. Начальные условия (3) предполагают, что в момент времени t=0 концентрация атомов в состоянии 1 равна , а в состоянии 2 - . Величина в уравнениях (1) - (2) представляет собой усредненную по углам и частотам интенсивность излучения в любой точке среды , в момент времени t и определяется из выражения

,        (4)

в котором - контур линии поглощения, а - интенсивность излучения в точке среды , в направлениях , , на частоте , в момент времени t.

Во втором параграфе выписано уравнение переноса излучения и выражения, определяющие контур линий поглощения для различных типов уширения линии. Для модели ППЧ вид уравнения переноса для спектральной интенсивности излучения следующий

,        (5)

где коэффициент поглощения в центре линии и функция источников S определяются выражениями

,        (6)

в которых с – скорость света, h – постоянная Планка, и статвеса состояний 1 и 2, - центральная частота линии.

Граничные условия для уравнения (5) можно записать в виде

       (7)

где - угол между осью Z и направлением внешнего излучения (рис.1), а - частотно-временная его форма. Cамосогласованная система уравнений баланса и переноса излучения (1) – (7) представляет собой систему интегродифференциальных уравнений (ИДУ), которая может быть решена только численно.

В следующем параграфе описаны численные методы расчета интегральных величин для среды, геометрия которой моделируется плоскопараллельным слоем, сферой и цилиндром. В основе численных алгоритмов лежали дискретно-разностные методы, сводящие вычисление интегралов (4) к суммированию интенсивностей по специально разработанным квадратурным формулам. Опишем кратко построенный алгоритм для сферической геометрии газовой среды. На рис.1 изображена расчетная сетка в центральном сечении шара. Полученная таким путем окружность делилась на заданное число равных дуг. Через точки касания этих дуг проводились прямые параллельные осям Х и Y. Узлы их пересечения  принимались  в  качестве  узлов  расчетной  сетки.  Для  расчета

Рис.1. Геометрия задачи и расчетной сетки центральной прицельной плоскости.

населенностей уровней, используя осевую симметрию, решение (1) – (2) достаточно получить в пространственных точках верхней половины круга центрального сечения. Интеграл (4) должен учитывать рассеянное излучение, приходящее в расчетный узел из всего объема среды и внешнее излучение с учетом его поглощения. Интегрирование по частоте и угловым переменным и проводилось по квадратурным формулам Эрмита, Гаусса и трапеций соответственно. Для этого использовались значения интенсивности рассчитанной в узлах расчетных сеток для набора прицельных плоскостей, на которые был поделен сферический объем по подобию центральной прицельной плоскости рис.1. Алгоритм был построен на основе формального решения уравнения переноса в интегральной форме. Таким образом, интеграл (4) рассчитывался через значения интенсивностей в узлах расчетных сеток, полностью охватывающих сферический объем.

В четвертом параграфе построена модель переноса излучения для двухуровневых атомов в условии частичного перераспределения по частотам (ЧПЧ). Модель ППЧ является приближенной. На самом деле частоты испущенного и поглощенного фотонов коррелируют между собой вполне определенным образом. Вид этой корреляции определяется функцией перераспределения по углам и частотам . Предположим, что атом поглощает фотон, который распространяется в направлении внутри телесного угла и имеет частоту в интервале . Излучается фотон с частотой в интервале и в направлении рассеяния внутри . Тогда функция перераспределения по углам и частотам определяется таким образом, что произведение есть совместная вероятность излучения фотона с заданными частотно-угловыми характеристиками. Наибольший практический интерес представляют усредненные по телесным углам функции перераспределения по частотам и

,

.        (8)

В выражениях (8) максимальная, а минимальная из частот и (где - безразмерная частота в единицах допплеровской ширины ). Функция описывает резонансную линию, возбужденное состояние которой уширено только радиационным механизмом. Следовательно, рассеяние фотонов в системе отсчета атома для функции остается полностью когерентным. Функция перераспределения предполагает, что возбужденный уровень уширен благодаря действию радиационного и ударного механизмов. Предполагается, что за промежуток времени между поглощением и излучением фотона все электроны, из-за соударений, полностью перераспределены по частотам возбужденного состояния. Поэтому модель рассеяния описывает процесс полного перераспределения фотонов в  системе  отсчета  атома.  Поведение  графиков

Рис. 2. Вероятность излучения фотона с частотой x для функции перераспределения при . Подписи к кривым указывают частоты поглощенных фотонов .

и показано на рис.2,3, где и получены численным интегрированием выражений (8). Отношение представляет собой вероятность излучения фотона на частоте при условии, что поглощенный фотон имел частоту . Для модели рассеяния вероятность излучения максимальна при близкой к . Эта вероятность резко уменьшается при удалении от , если , рис.2. Для фотонов поглощенных на частотах количество актов излучения больше в центре контура линии, поскольку здесь всегда больше количество атомов по сравнению с крылом линии. Такое поведение графика является следствием полной когерентности процесса рассеяния в системе отсчета атома для функции . Для графика характерно то, что при поглощении фотонов с частотами переизлученные фотоны попадают в частотный диапазон (функция симметрична относительно оси ординат, рис.3). Для фотонов поглощенных на частотах вероятность излучения близка к случаю полного перераспределения (пунктирная кривая рис.3). Количество фотовозбуждений определяется произведением , где усредненная по углам на частоте интенсивность. Это произведение помножим на совместную вероятность поглощения фотона на частоте и излучения на , а затем проинтегрируем по частотам всех поглощенных фотонов. Тогда  получится  скоростной  коэффициент,  который  отвечает  за

Рис.3. То же, что и на рис.2, только для функции при . Пунктирная кривая - вероятность излучения при естественном возбуждении .

возбуждение верхнего состояния, при распаде которого фотон будет иметь частоту . Обозначим за все такие населенности верхних подуровней, спонтанный распад которых дают фотоны частоты . Тогда для переменной , которая связана с полной населенностью возбужденного уровня соотношением , получим уравнение баланса вида

,

где коэффициенты и определяются из выражений

.

Уравнение баланса населенностей для основного уровня имеет вид (1). Показано, что для модели ЧПЧ возбужденное состояние расщепляется на систему подуровней. Поэтому при численной реализации такой задачи, двухуровневая модель атома по структуре подобна многоуровневой.

В следующем параграфе проведен анализ системы ИДУ. Учет осевой симметрии задачи позволил уменьшить ее размерность в 2n раз (n – число узлов по каждой из пространственных переменных X, Y, Z). Далее описан численный метод решения системы ИДУ в условии полного перераспределения по частотам, который сводит полученную задачу Коши к системе обыкновенных дифференциальных уравнений. Для модели ЧПЧ полученная система уравнений подобна двухуровневой модели и может быть решена тем же самым численным алгоритмом.

В шестом параграфе построена модель переноса излучения в ультрахолодной лазерной плазме при учете сферической геометрии среды, неоднородного пространственного распределения концентрации и температуры ионов, их радиального разлета со скоростью . На основании данных6

, радиальное распределение и моделировалось гаусовой формой, а скорость линейно зависящей от радиальной координаты:

, ,

,        (9)

где и - концентрация и температура ионов в центре капли в начальный момент времени, r – текущая координата, t – время, а - начальный характерный радиус гауссова распределения ионов. Характерный радиус гауссова распределения меняется со временем по закону

,

где - тепловая скорость электронов.

Для описания переноса излучения при разлете вещества использован подход, основанный на системе отсчета, связанной с наблюдателем. При этом в контуре линии поглощения оказываются взаимосвязанными угловые и частотные переменные. Так если тепловое движение частиц является единственным механизмом уширения линии, то допплеровский контур линии поглощения описывается выражением

.        (10)

Здесь определяет ширину линии в единицах допплеровской ширины в центре шара, т.е. . - есть ширина допплеровского контура в произвольной точке среды r, которая зависит от температуры и скорости движения частиц. В выражении (10) переменная обозначает смещение частоты от центральной частоты линии в единицах допплеровской ширины ; , где - тепловая скорость ионов в центре шара; - есть косинус угла между направлением движения атома и направлением рассеяния фотона. Численный алгоритм решения полученной задачи Коши учитывал описанные выше особенности, вызванные макроскопическим движением вещества и был построен по подобию стационарной среды.

В следующем параграфе изложена методика расчета скоростных радиационно-столкновительных коэффициентов ионизации, рекомбинации, фоторекомбинации, возбуждения и де-возбуждения атомов электронами. Эти скоростные коэффициенты рассчитывались в борновском приближении с учетом реальной структуры атомов (ионов) и основных их эмпирических параметров, таких, как энергия состояния, тип связи и уровня, строение внешней оболочки и др.

В заключительном параграфе приводятся результаты тестирования численных алгоритмов на аналитических решениях.

Во второй главе описаны результаты расчетов поглощения лазерного излучения, резонансного рассеяния и испускания света плотными парами натрия и лития. Постановки задач соответствуют моделям полного и частичного перераспределения по частотам и двухуровневому атому. Моделировалось воздействие на газовый объем, имеющий форму сферы и цилиндра, лазерного пучка при полном облучении среды и локальном облучении некоторого малого ее объема. Частотно–временная форма интенсивности входного лазерного излучения задавалась выражением

,        (11)

в котором - определяет длительность импульса, - частотная ширина лазерного излучения, имеющего лоренцовскую зависимость от частоты. Зависимость интенсивности от поперечной координаты задавалась прямоугольным и гауссовским распределением, где характерный размер луча.

В первом параграфе обсуждаются результаты по поглощению лазерного импульса парами натрия на длине волны . Изучено влияние оптической толщины среды на деформацию частотно-временной формы импульса, динамику фотовозбуждения и свечения паров Na.

Во втором параграфе описана картина фотовозбуждения и свечения паров натрия при воздействии на них лазерным пучком, размер которого меньше размера среды. Показано, что перенос излучения из области возбуждения на периферию приводит к фотовозбуждению и флуоресценции этой области среды.

Следующий параграф посвящен интерпретации экспериментальных результатов по лазерно-индуцированной флуоресценции паров атомов натрия7

(), в котором поведение интенсивности флуоресценции охарактеризовано, как «аномальное». Аномальность по мнению авторов7 следует понимать в смысле отсутствия насыщения интенсивности флуоресценции с ростом лазерной интенсивности (рис.4, кривая 2). Анализ условий эксперимента7 показал, что для описанных в нем параметров лазерного излучения и паров натрия когерентные эффекты проявляться не должны. Для объяснения результатов работы7 проводилось численное моделирование лазерно-индуцированной флуоресценции паров натрия для условий данного эксперимента. Исследовалось влияние переноса радиации из зоны возбуждения на периферию. Расчетные данные указали на то, что расширение зоны насыщения являться главным фактором «аномальной» флуоресценции в работе7 (рис.4, кривая 1).

Рис.4. Зависимость интенсивности флуоресценции от лазерной интенсивности, полученная в эксперименте7 (кривая 2) и численно (кривая 1) при =.

В четвертом параграфе решены задачи о резонансной флуоресценции паров натрия () и лития () в условии частичного перераспределения по частотам. Проведен анализ влияния различных функций перераспределения по частотам и геометрий газовой среды на форму спектрального контура линии и на эффективное время высвечивания паров (рис.5). Получено, что цилиндры, для которых отношение высоты к диаметру основания (кривая 2, рис.5) с точки зрения процессов переноса резонансного излучения эквивалентны бесконечно-длинным цилиндрам: выход излучения для них осуществляется через боковую поверхность. Кривая 3 на рис.5 рассчитана при использовании в численной модели полной функции перераспределения в виде линейной комбинации

,        (12)

где параметр дает долю возбужденных атомов, которые излучают когерентно , - радиационная, а - столкновительная ширины линии поглощения. Остальная доля () атомов благодаря столкновениям испытывает полное перераспределение по частотам в атомной системе отсчета. Использование линейной комбинации (12) в численной модели позволяет осуществить плавный переход от модели к путем увеличения начальной концентрации атомов. Следует отметить, что кривая 3 на рис.5 близка к кривой 4, полученной по асимптотической модели Бибермана-Холстейна для бесконечного цилиндра и допплеровского контура линии.

Рис. 5. Зависимости от оптической толщины среды фактора пленения Бибермана –Холстейна. - кривая 1; - 2 (модели ППЧ); - кривая 3 (полная функция перераспределения, выр.12); модель Бибермана – Холстейна – 4.

Рис.6. Зависимости от оптической толщины фактора пленения Бибермана - Холстейна паров лития для цилиндра с и различных функций перераспределения по частотам.

На рис.6 показано поведение от оптической толщины фактора пленения Бибермана-Холстейна паров лития для различных функций перераспределения по частотам. Кривая 1 соответствует функции , 2 - , 3 - . Поскольку функция представляется в виде линейной комбинации функций и , то для нее соответствующая кривая 1, с ростом оптической толщины плавно отделяется от кривой 2 и в дальнейшем приближается к кривой 3. Например, при параметр и . Для оптической толщины получим, и . Таким образом, поведение зависимости от рис.6 объясняется видом функций и , которые в численной модели через населенности уровней формируют спектральный контур линии рассеянного излучения. Динамика его изменения во времени определяет величину на достаточно больших временах, когда декремент интегральной по частоте интенсивности послесвечения паров выходит на стационарное значение.

Последний параграф посвящен анализу результатов эксперимента5. В этой работе эффективное время жизни паров лития () и натрия () измерено, используя декремент затухания во времени интегральной по частоте интенсивности послесвечения. Лазерным импульсом длительностью осуществлялась накачка паров на соответствующих резонансных переходах. Флуоресцентный сигнал фиксировался перпендикулярно направлению лазерного излучения в кювете цилиндрической формы. Зависимость фактора пленения q от оптической толщины среды для натриевых паров построена в диапазоне , рис.7. Монотонное уменьшение фактора q с ростом оптической толщины, при в эксперименте сменяется резким увеличением его значения, на что указывают экспериментальные маркеры рис.7. Такое поведение фактора q авторы5 объясняют проявлением эффектов частичного перераспределения по частотам.

Рис.7. Зависимости фактора пленения q от оптической толщины паров натрия () и лития (), полученные в эксперименте5.

Результаты численного моделирования зависимости подобной рис.7 для атомов натрия и лития показаны на рис.8 для цилиндра с . Кривая 1 соответствует полной функции перераспределения , 2 - , 3 - . В виде точек изображены данные, рассчитанные по модели ППЧ, которой соответствует функция . Для того, чтобы сопоставить расчетные данные с экспериментальными, на рис.8 так же как и на рис.7 построена асимптотическая кривая (пунктирная линия) для бесконечного цилиндра и допплеровского контура линии поглощения в соответствии с выражением q. Поведение точечной кривой и кривых 2, 3  на  рис.8  объясняется  видом  функций  перераспределения    и  .

Рис. 8. Расчетные зависимости от оптической толщины фактора пленения q паров натрия для различных функций перераспределения по частотам. Цилиндр с .

Зависимость q от для функции (кривая 1) примерно при отделяется от кривой 2, а в дальнейшем с ростом еще больше отходит от последней, плавно приближаясь к кривой 3. Подобные рис.9 зависимости были рассчитаны и для паров атома лития (). Соответствующая полной функции перераспределения (12) зависимость q от для лития помечена на рис.8 круглыми маркерами. Обращает на себя внимание тот факт, что последняя ближе к асимптотической кривой, на которой сосредоточены экспериментальные маркеры. Отличие расчетных кривых для лития и натрия объясняется следующим образом. Параметры атома лития и натрия таковы, что при одинаковой оптической толщине среды, температуры их паров и концентрации атомов отличаются (см. таблицу).

Параметры атомов натрия и лития

при одинаковой оптической толщине вдоль радиуса цилиндра .

Атом Li

Атом Na

Температура паров ,

694

Концентрация атомов ,

Допплеровская ширина ,

Параметр

0.84

0.94

Параметр

Данные таблицы указывают на то, что эффекты пленения для паров лития должны проявляться сильней, поскольку параметр у лития меньше, чем у натрия. Это следует из вида функции Фойгта : с ростом усиливается роль крыльев спектральной линии. В условии самопоглощения излучения характерное время выхода фотонов из среды на частотах далеких от центра линии, куда они перераспределяются при частичном перераспределении по частотам, уменьшается. Поэтому кривая зависимости q от для паров натрия лежит выше, чем для лития, а последняя, в свою очередь, ближе к асимптотической кривой, соответствующей допплеровскому контуру линии. Таким образом, в работе показано, что наблюдаемое в эксперименте5 немонотонное поведение зависимости фактора пленения Бибермана-Холстейна от оптической толщины паров натрия не следует приписывать эффектам частичного перераспределения по частотам. Следовательно, нужно искать другие причины такого поведения зависимости фактора q от в эксперименте5.

Третья глава содержит численные результаты моделирования процесса лазерной резонансной ионизации атомов натрия. Геометрия среды моделировалась цилиндром, а радиальное распределение интенсивности лазерного пучка имело однородное и гауссово распределение.

В первом параграфе построена модель процесса лазерной резонансной ионизации двухуровневых атомов натрия на переходе с длиной волны 589нм. Сначала в среде появляются затравочные электроны благодаря механизму ассоциативной ионизации. Затем они набирают энергию в сверхупругих процессах, дезактивируя возбужденные атомы, что и приводит к возникновению лавинной электронной ионизации газа. Такой механизм ионизации газа, получивший название фотоионизационного просветления, был впервые описан в работе8 для объяснения экспериментальных данных9

.

Для двухуровневого атома натрия система скоростных уравнений баланса и начальные условия при учете радиационно-столкновительных процессов имеет вид

         (13)

,.        (14)

В уравнениях (13) и (14) использованы следующие обозначения: – концентрация электронов, и – частоты радиационного фотовозбуждения основного и фототушения возбужденного уровней

.

Для концентрации электронов и температуры имеем уравнения

,        (15)

,        (16)

где и - потенциалы ионизации с основного и возбужденного уровней атома, - разность энергий между возбужденным и основным уровнями, и – коэффициенты электронного возбуждения и девозбуждения уровней, и - коэффициенты ионизации атомов электронным ударом из состояний 1 и 2, – константа скорости ассоциативной ионизации, и - константы скоростей трехчастичной рекомбинации на основной и возбужденный уровни, и – коэффициенты фоторекомбинации на соответствующие состояния 1 и 2. и - скорости передачи энергии при упругих столкновениях электронов с атомами и ионами, соответственно, которые при максвелловском распределении частиц по скоростям описываются выражениями

,

,        (17)

где - сечение электронно-атомных столкновений, - заряд электрона, и - температура ионов и атомов, , и - массы ионов, электронов и атомов, соответственно, - постоянная Больцмана.

В следующем параграфе проведен анализ численных результатов, полученных на основе описанной выше двухуровневой модели атома натрия. Частотно-временная форма лазерного излучения моделировалась выражением (11). При облучении среды импульсом с однородным в поперечном сечении распределением интенсивности процессы поглощения излучения приводят к неоднородной ионизации атомов по объему. У теневого торца цилиндра радиальная неоднородность проявляется сильнее, что объясняется поглощением лазерного излучения плотной средой. При облучении малой части паров натрия лазерным пучком исследовано влияние вторичной радиации на расширение ионизованного объема среды, вызванного переносом фотонов из области возбуждения к боковой поверхности цилиндра.

Рис.9. Поведение от времени электронной плотности и температуры, плотности атомов в основном и возбужденном состояниях. 2880, мкс, =4, Вт/см2, , см.

В третьем параграфе более детально изучен вопрос о влиянии переноса излучения на ионизационное просветление газа для многоуровневой модели атома натрия, которая включала 8 уровней: - уровни и ионизационное состояние. Рис.9 показывает поведение от времени электронной температуры и плотности населенности основного – 1 и резонансного - 2 уровней. Сплошные кривые получены для облученного торца цилиндра в точке соответствующей оси Z=, пунктирные – для осевой точки на теневом торце цилиндра Z=. В целом динамика развития резонансного лазерного разряда в парах натрия развивается следующим образом. Под действием лазерного излучения населенности уровней за короткое время (мкс) достигают значений близких к насыщающим (). Электроны, появившиеся благодаря ассоциативному механизму, быстро набирают энергию в сверхупругих процессах до величины эВ. Такое состояние с установившимися заселенностями уровней и энергией электронов длится до тех пор, пока не инициируется лавинная электронная ионизация атомов. При этом плотность электронов резко возрастает при мкс. С развитием электронной лавины заселенности резонансных, а так же и всех остальных возбужденных уровней падают из-за общего уменьшения атомов. В таком переходном режиме температура электронов растет до величины эВ.

Рис.10. Зависимость от времени коэффициента пропускания Т, представляющего собой отношение проходящей к падающей лазерной интенсивности интегральной по частоте, для различных значений : кривая 1 – ; 2 – ; 3 – ; 4 – , см-3. Параметры импульса соответствуют рис.10.

В начале переходного процесса, когда инициируется электронная лавина, среда «мутнеет», так как появляется дополнительный канал потерь, обусловленный нагревом электронов. На рис.10 изображены зависимости коэффициента пропускания парами лазерного излучения от времени для различных значений начальной концентрации атомов. Видно, что чем больше величина , тем раньше начинается ионизация и больше эффект «помутнения» среды, кривая 1. При маленьких концентрациях атомов процесс ионизации за время лазерного воздействия не успевает начаться вообще, кривая 4. Уменьшение коэффициента пропускания здесь вызвано ослаблением интенсивности входного импульса после прохождения им максимального значения при . Эффект потемнения вызван тушением электронами резонансных уровней. Данное явление было предсказано в работе8. Результаты моделирования показали, что при воздействии на газ резонансным излучением могут реализовываться два режима. Сначала, после выравнивания населенностей резонансных уровней, рождаются затравочные электроны. Их энергетический спектр устанавливается, а дальнейшее увеличение электронной плотности обязано ионизации атомов в электронных столкновениях. В условиях нарастания концентрации электронов, если скорость тушащих электронных процессов больше скорости ионизующих столкновений, лазерное излучение поглощается сильнее. Второй режим реализуется при условии, когда скорость ионизации электронным ударом превосходит скорость тушащих столкновений, в результате среда просветляется. Волна ионизации, в условии усиления поглощения лазерного излучения, формируется на облученной поверхности паров. Затем она распространяется вглубь среды, вызывая ее дополнительное просветление. Газ при этом может быть ионизован полностью. Обратим внимание также на следующий факт. Поведение кривых 1-3 рис.10 показывает, что на выходе из паров лазерный импульс разбивается на два подимпульса. Это связано с тем, что за счет ионизационного просветления среды количество поглощающих частиц становится меньше и уменьшается коэффициент поглощения. Следовательно, уменьшаются потери лазерной энергии на спонтанный распад и ее рассеяние в окружающее пространство. Для передачи энергии лазерного излучения через плотные газы этот факт может оказаться чрезвычайно важным.

Рис. 11. Пространственное распределение плотности возбужденных атомов в момент времени t=. Параметры модели следующие: цилиндр с отношением , , ; ; мкс; .

При воздействии на газ лазерным пучком размер зоны с высокой плотностью возбужденных атомов значительно больше характерного радиуса луча , рис.11. Расширение возбужденного объема обязано переносу фотонов из накачиваемой области за ее границу в радиальном направлении цилиндра. По мере приближения радиальной координаты к боковой поверхности цилиндра плотность возбужденных атомов уменьшается, что обязано выходу фотонов из среды. Просветление газа, вызванное ионизацией атомов электронным ударом, начинается в канале лазерного луча. Электронная лавина формируется в слоях близких к облученному торцу цилиндра. Далее она распространяется вглубь среды, обеспечивая дополнительное ионизационное просветление паров натрия, рис.12 а. Происходит постепенная ионизация области лазерного канала, рис.12 б. Кроме ионизационного просветления паров натрия в канале лазерного луча, ионизованный объем значительно расширяется. Характерный размер плазменного канала значительно больше размера лазерного луча . Сильное фотовозбуждение приканальной зоны, обязанное переносу фотонов в радиальном направлении цилиндра, приводит здесь к появлению свободных электронов, которые затем набирают энергию в сверхупругих столкновениях. Далее, когда инициируется электронная лавина в канале луча, процесс ионизации атомов электронным ударом развивается и вне его зоны. Таким образом, процессы переноса излучения приводят к сильному фотовозбуждению среды за пределами лазерного луча, в результате чего происходит значительное расширение ионизованного объема.

Рис. 12. Пространственное распределение плотности электронов – а) и возбужденных атомов – б), мкс. Параметры модели соответствуют рис.11.

В последнем параграфе проведены расчеты характерного времени ионизации паров натрия , используя интенсивность рассеянного излучения. При развитии лавинной ионизации атомов их концентрация падает, поэтому интенсивность линии со временем будет уменьшаться. Тогда может быть определено численно по декременту затухания интегральной по частоте интенсивности из соотношения предполагая, что на достаточно больших временах уменьшение интенсивности излучения во времени носит одноэкспоненциальный характер. Расчеты показали, что при отсутствии полной ионизации объема за время воздействия импульсом, рассчитанные по различным направлениям имеют большой разброс по значению. При полной ионизации паров величины , соответствующие различным оптическим путям, сходятся к одному и тому же значению, которое и представляет собой характерное время ионизации. Следовательно, по декременту затухания интенсивности рассеянного излучения можно оценивать характерное время ионизации паров.

В четвертой главе содержится описание расчетных данных о поглощении, рассеянии и испускание света в ультрахолодной лазерной плазме. В первом параграфе описаны параметры УП, соответствующие численной модели главы 1. Контур линии поглощения моделировался фойгтовской зависимостью от частоты, а возбуждающий лазерный импульс имел однородное распределение интенсивности в поперечном сечении пучка, которое больше размеров плазмы.

В следующем параграфе описаны результаты расчета поглощения лазерного излучения кальциевой плазмой на переходе с длиной волны . На рис.13 показана форма проходящего излучения при , , и радиусе УП . Кривые 1-6 здесь соответствуют точкам шара с координатами относительно центра: –0.025, –0.013, 0 (центр шара), 0.013, 0.025, 0.05см. Пунктирная кривая соответствует проходящему излучению для малой лазерной интенсивности (, из которой видно, что его частотная зависимость имеет симметричный по частоте характер. При больших интенсивностях проявляется нелинейное поглощение (кривые 1-6). В начале трассы и далее до центра шара идёт поглощение красного крыла линии движущимся навстречу слоем плазменной среды. Однако далёкие крылья поглощения в этой части среды позволяют поглощать и синее крыло спектра. Так как на входе в среду излучение является максимальным, нелинейный эффект проявляет себя в полной мере. После прохождения центра шара наиболее интенсивно идёт поглощение в синей области спектра, однако, уже ослабленные по интенсивности оба участка спектра вызывают более слабое нелинейное поглощение и синяя область спектра поглощается всё более в линейном режиме. В результате, интенсивность проходящего излучения приобретает несимметричный по частоте характер.

Рис.13. Частотная зависимость проходящего излучения. Скорость разлета =6000 см/с.

В третьем параграфе проведен анализ расчетных данных по рассеянию лазерного излучения кальциевой плазмой. Форма частотной зависимости интенсивности оказалась чувствительной к углу рассеяния. Результаты расчетов, приведенные на рис.14 подтверждают это и соответствуют: , , , 1 – рассеяние «вперед», 2 – рассеяние «назад». Кривая 3 на рис.14 соответствует рассеянию на угол 750, из вида которой следует, что спектр рассеяния почти симметризуется. Это связано с конкуренцией процессов рассеяния «назад» и «вперед». Отметим, что с ростом скорости разлета соответствующие смещения в фиолетовую и красную области усиливаются. Таким образом, предсказана частотная асимметрия в спектре рассеянного УП излучения.

Рис.14. Спектр рассеянного излучения.

В следующем параграфе приводятся результаты расчета интенсивности послесвечения кальциевой плазмы. Полученные результаты показали, что форма спектральной линии не зависит от угла рассеяния, а частотные зависимости интенсивности оказываются смещенные в красную область. Излучение выходящее из шара на своем пути в основном, будет поглощаться в синей области разлетающимися ионами. Излучение в красной области спектра, приходящее к границе среды из глубинных слоев, будет испытывать менее заметное поглощение.

В последнем параграфе решены задачи диагностики УП, используя спектральные характеристики. Получено, что в оптически тонком случае ширина контура линии поглощения и излучения линейно возрастает с ростом скорости разлета.

В оптически плотном случае для определения концентрации нормальных ионов в плазме определялась частота, на которой входная интенсивность после поглощения в плазменном шаре уменьшится в 10 раз. Назовем такое значение частотой отсечки и обозначим ее в относительных допплеровских единицах . На этой частоте показатель экспоненты в законе Бугера должен быть равен значению 2.3. В начальный момент времени плазма покоится. Тогда, используя для оценок допплеровскую форму линию (10) при , можем получить выражение для :

.        (18)

Подставляя в (18) выражение для из (9) и для из (6), выполняя интегрирование, получим для начального момента времени (t=0)

, ,        (19)

где - сечение резонансного поглощения. График зависимости от показан на рис.15 при указанных выше параметрах иона кальция.

Рис. 15. Зависимость концентрации ионов кальция от частоты отсечки для стационарной среды.

Рис. 16. Зависимость скорости разлета от частоты отсечки для различных концентраций ионов .

На начальных этапах формирования плазмы, когда скорость разлета ионов близка к нулю по частоте отсечки можно определить концентрацию ионов . В более поздние моменты времени, когда скоростью разлета УП пренебречь нельзя, соотношение (18) перепишется в виде

.        (20)

Для вычисления интеграла (20) учитывался тот факт, что по заданному лучу зрения плазма представляет собой два плоских слоя, разлетающихся в противоположных направлениях относительно центра шара. Тогда выражение (20) сводится к двум слагаемым

       (21)

где введено обозначение и учтен тот факт, что , а тепловая допплеровская ширина определяется через тепловую скорость и длину волны перехода следующим образом . В последнем уравнении концентрация может быть определена по частоте отсечки для покоящейся плазмы по графику рис.15. Тогда выражение (21) представляет собой трансцендентное уравнение, из численного решения которого можем найти зависимость скорости разлета от частоты . Для этого интегралы в (21) заменялись квадратурными суммами, соответствующими методу трапеций. На рис. 16 показаны графики зависимости от при различных значениях : 1 - ; 2 - ; 3 – ; 4 – , определяемых по рис.15. Поведение кривых на рис. 16 указывает на то, что в оптически плотной плазме с ростом скорости разлета увеличивается частота отсечки. Радиальный разлет ионов приводит к уширению контура линии поглощения и к увеличению частоты отсечки.

В пятой главе представлены результаты моделирования динамики возбуждения, ионизации и свечения искусственного бариевого облака под действием солнечного света.

В первом параграфе проводится анализ физических условий, протекающих в БО под действием солнечного света. Построены многоуровневые модели атома и иона бария на основе данных о сечениях фотовозбуждения и распределения по спектру интенсивности солнечного излучения. Получено, что для адекватного учета сложной структуры энергетических уровней бария достаточно включить в модель 15 атомных и 5 ионных состояний. При оптических плотностях облака превышающих 30 перенос радиации необходимо учитывать на атомном переходе с длиной волны 553.5нм и на трех ионных с длинами волн 455.4нм, 493.4нм и 614.2нм. Систему балансных уравнений многоуровневых атомов и ионов, появляющихся при ионизации первых, можно представить в виде

                               

       (22)

                                       

,        (23)

где , а m > n, , – частоты фотоионизации i-го уровня атома, , - число учитываемых уровней атома и иона. Система уравнений (22) – (23), описывающая кинетику многоуровневого атома и иона, совместно с уравнениями (4)-(7), определяющих коэффициенты для набора оптически плотных линий, решалась по описанному в главе 1 численному алгоритму в приближении допплеровского контура линии поглощения, модели ППЧ и стационарной среды.

Во втором параграфе описаны результаты расчетов динамики поглощения солнечного света бариевым облаком. Выявлены два режима просветления среды. Первый режим – «быстрый», который характеризуется скачкообразным увеличением проходящей солнечной интенсивности до некоторой величины. Этот режим вызван быстрым уходом атомов из основного состояния на метастабильные уровни, в результате чего скачкообразно уменьшается коэффициент поглощения. Второй режим просветления – «медленный». Для него характерно медленное уменьшение коэффициента поглощения из-за фотоионизации атомов с метастабильных состояний под действием солнечной радиации.

В третьем параграфе получена кинетика фотовозбуждения и фотоионизации бария. Процессы поглощения солнечной радиации приводят к сложной пространственно-временной картине возбуждения и ионизации бария. Для облака в форме шара, различные точки среды соответствуют различным оптическим путям, относительно проходящего солнечного света. Поэтому из-за поглощения солнечного света скорости фотовозбуждения в узловых точках будут отличаться друг от друга. В результате этого процесс фотовозбуждения протекает неоднородно по координате (рис.17).

В следующем параграфе обсуждаются результаты моделирования, касающиеся фотоионизации облака. Проводится расчет характерного времени фотоионизации БО по данным о динамике уменьшения во времени концентрации электронов. Расчетные данные указали на то, что с ростом оптической толщины среды характерное время фотоионизации бариевого облака увеличивается. Показано также, что декремент затухания (нарастания) интенсивности свечения атомной (ионной) линии хорошо описывает характерное время ионизации атомов бария солнечным светом.

Рис.17 Динамика изменения населенности возбужденного уровня атома бария в различных узлах теневой стороны сферы. Номера кривых соответствуют номерам точек рис.1.

Рис.18. Радиальное распределение интенсивности выходящего из облака излучения для атомной линии (553.5нм) в момент времени t = 1c. - диаметр облака, R - расстояние от точки выхода излучения в центральном сечении облака до линии визирования (наблюдения). Кривая 1 соответствует углу рассеяния ; 2 - ; 3 - ().

В пятом параграфе приводятся результаты расчета частотно-угловых и пространственно-временных характеристик рассеянного БО излучения. Динамика поведения зависимости интенсивности от частоты свидетельствует о наличии реабсорбции контура спектральной линии атома в начале, а иона в конце ионизации, когда облако является оптически плотным для излучения в соответствующих спектральных диапазонах. Получено, что процессы переноса радиации существенно влияют и на радиально-угловое распределение интенсивности атомной и ионной линий бария. Например, угловая зависимость рассеянной радиации сильно зависит от момента времени ионизации, что вызвано уменьшением оптической плотности БО по мере ионизации атомов. Так в начальные моменты времени, когда среда оптически плотная для излучения, рассеяние происходит преимущественно в направлении противоположном солнечному свету. По мере ионизации и уменьшения оптической плотности БО, рассеяния становится более изотропным в пространстве. Кроме того, в моменты времени, когда облако является оптически плотным, радиальное распределение интенсивности ведет себя неодинаково в зависимости от угла, рис.18. Расчетные данные позволили объяснить цветовую окраску бариевого облака и ее изменчивость от угла наблюдения, которая получена в натурном эксперименте10.

В последнем параграфе проведено сравнение расчетных и экспериментальных данных по динамике яркости БО в спектральных линиях, представляющих интерес с практической точки зрения. Динамика яркости атомной линии (553.5нм) (кривая 2 рис.19) качественно описывает ее поведение, полученное в эксперименте АМРТЕ11 (кривая 1). Такое же качественное совпадение с результатами эксперимента СПОЛОХ12 получено и для динамики яркости ионной линии бария (455.4нм).

В заключении сформулированы основные результаты работы.

Рис.19. Динамика яркости БО на атомном переходе (553.5нм) полученная численно - 2 и в эксперименте AMPTE - 1.

Основные результаты диссертационной работы

Работа посвящена численному исследованию процессов переноса излучения в многоуровневых динамических средах в задачах о радиационном воздействии на газ и плазму. В свете современных тенденций развития лазерной физики, физики плазмы, оптики и спектроскопии, разработанные математические модели и алгоритмы вносят значительный вклад в развитие научного направления «Перенос резонансного излучения в оптически плотных, динамических, газовых и плазменных средах». Используя их, получен ряд результатов и выводов, которые имеют приоритетный характер в данной области исследования и состоят в следующем.

1. Построена модель процесса взаимодействия солнечного излучения с многоуровневыми нестационарными средами, учитывающая спектральный состав внешнего излучения, реальную структуру атома и иона, перенос излучения в оптически плотных спектральных линиях.

2. Разработаны численные алгоритмы решения систем интегродифференциальных уравнений, сводящую задачу Коши к системе обыкновенных дифференциальных уравнений для газовых сред в форме плоского слоя, шара и цилиндра.

3. Численные методы решения систем ИДУ позволяют учитывать следующие физические факторы: частичное перераспределение по частотам; макроскопическое движение вещества; неоднородное пространственное распределение вещества и его температуры в газовой среде; совместное влияние лоренцевского и допплеровского механизмов уширения спектральной линии; локальное возбуждение малой части газа световым пучком.

4. Исследовано поглощение и рассеяние лазерного излучения парами натрия на резонансном переходе с длиной волны 589нм для различных геометрий газовых сред и параметров лазерного излучения. При локальном фотовозбуждении лазерным импульсом малого объема среды распространение излучения из зоны луча на периферию приводит к расширению области с высокой плотностью возбужденных атомов. Показано, что это явление является главным фактором «аномальной» флуоресценции, полученной в эксперименте7.

5. В условии ЧПЧ рассчитаны контуры линий свечения и характерные времена высвечивания паров натрия и лития для различных функций перераспределения фотонов по частотам и геометрий газовой среды. Получено, что цилиндры, для которых отношение высоты к диаметру основания с точки зрения процессов переноса резонансного излучения эквивалентны бесконечно-длинным цилиндрам, а эффекты ЧПЧ уменьшают эффективное время жизни возбужденного состояния по сравнению с моделью ППЧ.

Результаты численного исследования зависимости фактора пленения Бибермана-Холстейна от оптической толщины паров натрия показали, что наблюдаемое в эксперименте5 «аномальное» (немонотонное) его поведение не следует приписывать эффектам частичного перераспределения по частотам.

6. Проведено численное моделирование процесса лазерной резонансной ионизации паров натрия, в основе которой лежит сверхупругий нагрев свободных электронов, образованных благодаря ассоциативному механизму. Исследованы эффекты потемнения среды, вызванного тушением электронами возбужденного состояния атомов и колебание фронта ионизационно просветленной зоны.

Получено расширение плазменного канала в форме ореола поверх лазерного пучка, обусловленное процессами переноса рассеянного резонансного излучения. Предложена методика определения характерного времени ионизации среды под действием лазерного излучения, которая основана на измерении интенсивности рассеянного излучения.

7. Построена модель и проведено численное исследование испускания света, поглощение и рассеяние лазерного излучения ультрахолодной расширяющейся плазмы на резонансном переходе иона кальция (). В проходящем лазерном излучении предсказано смещение контура линии в красную частотную область. В контуре линии рассеянного излучения наблюдается смещение интенсивности в красную область спектра для излучения выходящего с теневой стороны плазмы, и в фиолетовую область - при рассеянии в противоположном лазерному излучению направлении. Предложен метод расчета скорости разлета оптически плотной плазмы, основанный на измерении интенсивности проходящего излучения в крыле линии.

8. Решена задача об ионизации и свечении искусственного бариевого облака под действием солнечного света. Многоуровневые модели атома и иона бария построены на основе данных для сечений поглощения и распределения по спектру интенсивности солнечного излучения. Получено, что с ростом оптической толщины облака характерное время его ионизации увеличивается. Спектральные характеристики рассеянного излучения объясняют цветовую окраску облаков и ее изменчивость в зависимости от направления, что наблюдается в натурных экспериментах. Расчетная зависимость динамики яркости бариевого облака на атомном (553.5нм) и ионном (455.4нм) резонансных переходах качественно согласуется с экспериментальными данными (эксперименты11,12).

Основные результаты диссертации опубликованы в следующих работах:

  1. Гольбрайх Е.И., Косарев Н.И., Николайшвили С.Ш., Шкедов И.М. Ионизация оптически прозрачного бариевого облака // Геомагнетизм и аэрономия. – 1990. – Т. 30. – № 4. – С.688 – 690.
  2. Косарев Н.И., Шкедов И.М. Распространение широкополосного излучения в бариевом слое // Оптика атмосферы. – 1991. – Т. 4. – №11. – С.1172 – 1178.
  3. Косарев Н.И., Шкедов И.М. Распространение солнечного излучения в искусственном бариевом облаке // Оптика атмосферы и океана. – 1993. – Т.6. – № 10. – С.1298 - 1306.
  4. Kosarev N.I., Shkedov I.M. Laser-induced fluorescence of the optical thick sodium vapours // Proceedings of the second Chine-Russian symposium on laser physics and laser technology. July 14 - 17, 1995. - Harbin, China, 1995. - P. 40 - 41.
  5. Косарев Н.И., Шкедов И.М Поглощение лазерного излучения плотными парами натрия // Оптика атмосферы и океана. – 1995. – Т.8. – № 12. – С.1752 - 1756.
  6. Косарев Н.И., Шапарев Н.Я., Шкедов И.М. Компьютерное моделирование радиационных эффектов в бариевых облаках // Актуальные проблемы информатики, прикладной математики и механики / Под ред. В.В. Шайдурова. – Новосибирск – Красноярск: Изд-во СО РАН, 1996. – Ч. 2. – С. 82 - 89.
  7. Kosarev N.I., Shkedov I.M. Excitation and emission of sodium vapours under action of the laser beam // Proceedings of the third Russian-Chinese symposium on laser physics and laser technology. October 8 – 10, 1996. – Krasnoyarsk, Russia, 1996. - P. 75 - 77.
  8. Kosarev N.I., Shkedov I.M. Radiation imprisonment in vapours of the sodium under photo-excitation of its atoms by the laser impulse // Proceedings of the third Russian-Chinese symposium on laser physics and laser technology. October 8 – 10, 1996. – Krasnoyarsk, Russia, 1996. – P. 78 - 81.
  9. Косарев Н.И., Шкедов И.М. Рассеяние солнечного света ионным бариевым облаком // Оптика атмосферы и океана. – 1999. – Т. 12. – № 1. – С. 30 - 35.
  10. Kosarev N.I., Shkedov I.M. Laser-induced ionization of sodium vapours // Proceedings of the 5-th Russian-Chinese Symposium on Laser Physics and Laser Technology. October 23 – 28, 2000. - Tomsk, Russia, 2000. - P. 31 - 34.
  11. Kosarev N.I., Shkedov I.M. Transfer of laser radiation in gases // Proceedings of the 5-th Russian-Chinese Symposium on Laser Physics and Laser Technology. October 23 – 28, 2000. - Tomsk, Russia, 2000. - P. 49 - 52.
  12. Kosarev N.I., Shkedov I.M. Scattering of laser radiation by sodium atoms at partial redistribution on frequencies // Proceedings of the 6-th International symposium on laser physics and technologies. August 18 – 24, 2002. – Harbin, China, 2002. – P. 121 - 126.
  13. Kosarev N.I. Laser-induced gas transparency in conditions of radiating transfer // Proceedings of the 7-th Russian-Chinese symposium on laser physics and laser technologies. December 20 - 24, 2004. – Tomsk, Russia, 2004. – P. 296 - 300.
  14. Косарев Н.И. Лазерная резонансная ионизация атомов натрия в условиях переноса излучения // Математическое моделирование. – 2005. – Т. 17. – № 5. – C. 105 - 122.
  15. Косарев Н.И., Шапарев Н.Я. Ионизационная прозрачность газа, индуцированная резонансным лазерным воздействием // Оптика атмосферы и океана. – 2006. – Т. 19. – № 2-3. – С. 216 - 220.
  16. Косарев Н.И. Комментарии к вопросу об аномальной флуоресценции паров атомов натрия на длине волны 589.6нм // Оптика и спектроскопия. – 2006. – Т. 101. – № 1. – С. 64 - 70.
  17. Косарев Н.И. Перенос излучения в искусственном бариевом облаке при его фотоионизации солнечным светом // Математическое моделирование. – 2006. – Т. 18. – № 12. – С. 67 – 87.
  18. Косарев Н.И., Шапарев Н.Я. Резонансная лазерная ионизация паров натрия при учете радиационного переноса // Квантовая электроника. – 2006. – Т. 36. – № 4. – С. 369 - 375.
  19. Kosarev N.I., Shaparev N.Ya. Resonance fluorescence of the calcium plasma // Proceedings of the 8-th Sino-Russian symposium on laser physics and laser technologies. August 10 – 15, 2006. – Harbin, China, 2006. – P. 11 - 15.
  20. Kosarev N.I., Shaparev N.Ya. Ionisation bleaching of the sodium vapours // Proceedings of the 8-th Sino-Russian symposium on laser physics and laser technologies. August 10 – 15, 2006. – Harbin, China, 2006. – P. 20 - 25.
  21. Kosarev N.I., Shaparev N.Ya. Ionization transparency of the gas induced by resonant laser influence // Proc. SPIE. – 2006. - V.6263. – Р.124-131.
  22. Косарев Н.И. Формирование контура спектральной линии при частичном перераспределении по частотам // Оптика и спектроскопия. – 2007. – Т. 102. – № 1. – С. 13 - 19.
  23. Косарев Н.И. Эффективное время высвечивания паров лития // Оптика и спектроскопия. – 2007. – Т. 102. – № 5. – С. 718 – 724.
  24. Косарев Н.И. Моделирование переноса излучения при неполном перераспределении по частотам // Математическое моделирование. – 2008. – Т. 20. – № 3. – С. 87 - 97.
  25. Косарев Н.И. Распад возбужденного состояния атомов натрия с учетом пленения излучения // Оптика и спектроскопия. – 2008. – Т. 104. – №1. – C. 5 - 8.
  26. Косарев Н.И., Шапарев Н.Я. Поглощение резонансного излучения в ультрахолодной лазерной плазме // ДАН. – 2008. – Т. 421. – № 6. – С. 1 – 3.
  27. Косарев Н.И., Шапарев Н.Я. Resonance optical characteristics of ultracold plasma // Proceedings of the 9-th Russian-Chinese symposium on laser physics and laser technologies. October 26 – 31, 2008. – Tomsk, Russia, 2008. – P. 80 - 84.
  28. Kosarev N.I. Numerical investigation of the escape factor of lithium and sodium vapour at partial frequency redistribution // J. Phys. B: At. Mol. Opt. Phys. – 2008. – V. 41. – Р. 225401-1 – 225401-8.
  29. Kosarev N.I., Shaparev N.Ya. Absorption and scattering of resonance laser radiation in ultracold optical dense plasma // J. Phys. B: At. Mol. Opt. Phys. – 2008. – V. 41. – Р. 235701 -1 – 235701 – 3.
  30. Косарев Н.И., Шапарев Н.Я. Резонансные оптические характеристики ультрахолодной лазерной плазмы // Квантовая электроника. – 2009. – Т. 39. – № 12. – С. 1112 – 1116.

Подписано в печать __________ 2010 г. Формат 60 80/16.

Усл. печ. л.2. Тираж 100 экз. Заказ № __.

Отпечатано на участке оперативной полиграфии

Сибирского юридического института МВД России.

660050, г. Красноярск, ул. Кутузова, 6.


1 Holstein T. Imprisonment of resonance radiation in gases. I, II // Phis. Rev. – 1947. – V.72. – P.1212 – 1233. – 1951. – V.83. – P.1159 - 1168.

2 Биберман Л.М. К теории диффузии резонансного излучения // ЖЭТФ. – 1947. – Т.17. – C.416-426.

3 Гаврилюк А.П., Краснов И.В., Шапарев Н.Я. Оптическое удержание низкотемпературной плазмы с резонансными ионами // Письма в ЖЭТФ. – 1996. – Т. 63. – С.316 - 321.

4 БулышевА.Е., ПреображенскийН.Г., Суворов А.Е. Перенос излучения в спектральных линиях // УФН. – 1988. – Т.156. – Вып.1. – С.153-175.

5 A. Romberg, H.-J. Kunze. Experimental investigation of the radiative transport of the resonance lines of sodium and lithium // J. Quant. Spectrosc. Radiat. Trans. – 1988. –  V.39. – No.2. – Р. 99 - 107.

6 T.C. Killian, T. Pattard, T. Pohl, et al., Ultra-cold neutral plasmas // Phys. Rep. – 2007. – V. 449. – P. 77-130.

7 J.M. Salter, D.D. Burgess, N.A. Ebrahim. Anomalous behaviour in the saturation of the sodium D lines under high power laser illumination // J. Phys. B: Atom. Molec. Phys. – 1979. – V.12. – No.24. – P.759-762.

8 Шапарев Н.Я. Ионизационное просветление газа // ЖЭТФ. – 1981. – Т. 80. – C. 957 – 963.

9 Lucatorto T.B., McIlrath T.J. Efficient laser production of a Na+ ground-state plasma column: Absorption spectroscopy and photoionization measurement of Na+ // Phys. Rev. Lett. – 1976. – V.37. – No.7. – P.428 - 432.

10 Беликов Ю,Н., Гуревич А.В., Николайшвили С.Ш. Цветовая диагностика искусственных облаков в околоземном космическом пространстве // Космические исследования. – 1993. – Т.31. – Вып.1. – C. 108 – 114.

11 Bernhardt P.A., Roussel – Dupre R.A., Pongratz M.B. et al. Observation and theory of the AMPTE magnetotial barium releases // J. Geophys. Res. – 1987. – V. 92. – P. 5777-5794.

12 Адейшвили Т.Г., Бочаров А.А., Дорофеев В.С., и др. Фотометрические и телевизионные измерения в ракетном эксперименте «СПОЛОХ» с инжекцией бария в ионосфере Земли. – М., 1977. – 18 с. (Препринт института космических исследований АН СССР: №342).

 



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.