автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование теплового баланса, нейтрального и ионного состава мезосферы и нижней термосферы

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

Жаркова Юлия Сергеевна

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛОВОГО БАЛАНСА, НЕЙТРАЛЬНОГО И ИОННОГО СОСТАВА МЕЗОСФЕРЫ И НИЖНЕЙ ТЕРМОСФЕРЫ

05.13.18 — математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Калининград, 2008

Работа выполнена в Российском государственном университете имени Иммануила Канта

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

Медведев Владимир Васильевич

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

профессор Лобанов А.И.

доктор физико-математических наук, профессор Кореньков Ю.Н.

Ведущая организация: Институт Математического моделирования РАН,

г. Москва

Защита состоится «/ ? » —СШ^ХАЛ/Л . 2008 г. в № час мин

на заседании диссертационного совета К 212.084.10 при Российском государственном университете имени Иммануила Канта, г. Калининград, ул. А.Невского, 14, ауд. /С%

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

С диссертацией можно ознакомиться в библиотеке РГУ им. И. Канта (г.Калининград, ул. А.Невского, 14)

Автореферат разослан «/?» ЫлЯ^МА 2008 г.

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

диссертационного совета Токарь В.Г.

А.

к

Актуальность работы

Как известно, область высот 50-250 км труднодоступна для экспериментальных исследований В настоящее время разработаны различные методики экспериментов, но проведение этих экспериментов несет эпизодический характер и дает различные ошибки измерений Несмотря на значительные экспериментальные успехи, достигнутые в последнее время в области исследования верхней атмосферы, математическое моделирование остается основным (достаточно дешевым) методом исследования этой области высот

В теоретическом плане изучение этой области высот затруднено необходимостью учета сложных динамических и фотохимических процессов, таких, как турбулентное перемешивание, переходящее в молекулярную диффузию, поглощение нейтральным составом солнечного излучения и его эмиссия, большая плотность и многокомпонентность состава, малые компоненты О, Оэ, С02, Ог(1Д8), Н20, N0, концентрация которых существенно меньше основных N2, 02, но которые могут играть существенную роль как в тепловом балансе, так и в образовании ионосферы

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

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

Цель работы

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

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

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

4 Исследовать возможности применения результатов модели в задачах распространения электромагнитных волн

Направления исследования

Уравнение теплопроводности, источники и стоки тепла, совместное решение уравнений непрерывности для ионов и нейтралов

Методы исследования, достоверность и обоснованность результатов

Метод исследования основан на численном решении систем нелинейных, связанных уравнений гидродинамики (уравнений непрерывности, движения, теплопроводности)

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

Правильность выбранных методов проверяется на контрольных примерах Полученные модельные расчеты проверяются путем сравнения их с имеющимися экспериментальными данными

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

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

1 Расчет источников и стоков тепла, обусловленных поглощенным солнечным излучением, фотохимическими и динамическими процессами мезосферы и нижней термосферы

3 Результаты самосогласованных расчетов высотно-временного распределения температуры нейтрального газа с нейтральным и ионным составом

4 Возможное теоретическое обоснование инверсии температуры в области турбопаузы

5 Аналитическое выражение высотного распределения окиси азота в области мезосферы

Научная новизна

1 Разработана модель теплового режима мезосферы и нижней термосферы, позволяющая рассчитать высотно-временное распределение состава и температуры

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

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

4 Методом вычислительного эксперимента получена инверсия температуры в области турбопаузы и дано ее возможное теоретическое объяснение

5 Получено новое аналитическое распределение окиси азота в области мезосферы

Практическая полезность работы

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

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

Апробация работы

Результаты работы обсуждались и докладывались на следующих конференциях и семинарах

Международной научной конференции, приуроченной к 200-летию со дня рождения великого немецкого математика Карла Густава Якоби и 750-летию со дня основания г Калининграда(Кенигсберга) «Избранные вопросы современной математики», г Калининград, 2005г, 4-8 апреля;

Physics of Auroral Phenomena 29th Annual Seminar, Polar Geophysical Institute, Apatity, 2006r,

6-th International conference Problems of Geocosmos, Saint-Petersburg,

2006,

семинарах математического факультета РГУ им И Канта

Публикации

Основное содержание диссертации отражено в 11 научных работах

Структура и объем работы

Диссертация состоит из введения, четырех глав, заключения, библиографического списка, включающего 67 наименований Работа изложена на 130 листах машинописного текста, содержит 20 рисунков, 7 таблиц

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

Во введении обоснована актуальность проблемы теоретических исследований верхней атмосферы и ионосферы

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

полученные к настоящему моменту по атомарному и молекулярному кислороду, озону, электронно-возбужденному молекулярному кислороду, азотным компонентам, ионному составу и электронной концентрации, на основе современной литературы

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

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

- накоплен большой теоретический и экспериментальный материал, который позволяет более надежно, чем в области Б моделировать нейтральный и ионный состав Е — области,

- что касается О — области, то существующие в настоящее время экспериментальные данные по высотному распределению малых компонент О, 03, 02('/1я), N0 и по ионному составу, малочисленны, нерегулярны и разнородны, получены для различных широт, сезонов и моментов времени,

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

- существенным недостатком теоретических моделей является отсутствие надежной теории, способной объяснить экспериментальные высотные распределения окиси азота в области мезопаузы и Оз, 02('Л^ в области 85-100 км Данное обстоятельство приводит к определенным трудностям в выяснении относительной роли источников ионизации, что в

конечном итоге затрудняет получение высотного профиля Ые, удовлетворительно согласующегося с экспериментальными данными,

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

В главе 2 представлена постановка задачи математического моделирования нейтрального и ионного состава мезосферы и нижней термосферы (области О, Е ионосферы) Предлагается одномерная, нестационарная модель мезосферы и нижней термосферы с учетом молекулярного и турбулентного переноса для области высот 50-250 км Модель позволяет рассчитать пространственно-временное распределение следующих компонент

ы2,о,о2,ог,оСо),о{х8),огС^),о2(}^ ),н,н2,он,н2о,но2,н2о2, Ж(4 Я), Щ2П), N0, м?2, со, со2, о+, 02, , 04+, , щ, , о2, Г+, У,

где У+ - суммарная концентрация положительных ионов-связок, V -суммарная концентрация отрицательных ионов Рассматриваются уравнения непрерывности для ионов и нейтральных частиц Выбрана фотохимическая схема для их описания, основанная на современных данных о скоростях химических реакций Основными входными параметрами представляемой модели являются начальные, верхние и нижние граничные условия, сечения фотоионизации и фотопоглощения, поток ультрафиолетового излучения Солнца, коэффициенты химических реакций, частоты столкновений, коэффициент турбулентного перемешивания На примере расчета высотного распределения [О] и [N0] показана возможность использования распределенных граничных условий для определения граничных значений, и отмечено, что в случае малого количества информации об искомой функции, возможно получение результатов, не соответствующих реально существующим

Проведено сравнение результатов численного расчета высотного распределения [О] и [<Э2] различными вариантами метода прогонки потоковой, немонотонной, матричной, с заменой переменной, обыкновенной Показано, что рассчитанные [О] и [С?2] этими методами находятся в хорошем согласии между собой, что говорит о правильности построенных алгоритмов

На основе результатов математического моделирования высотно-временного распределения [Лг(45')] и [М?] показано выполнение условий, близких к фотохимическому равновесию Данное обстоятельство дало возможность составить несложную систему уравнений непрерывности для [ 7У(45')] и [N0] Решение этой системы позволило получить простое аналитическое выражение для расчета высотного распределения [Ж)], зависящее от температуры нейтрального газа и концентрации молекулярного кислорода Полученное выражение объясняет высотное поведение [М7] как для спокойных, так и возмущенных условий На основе полученного выражения объясняется увеличение [N0], уменьшение доли ионов-связок и эффективного коэффициента рекомбинации во время зимней аномалии области Б, а также увеличения концентрации ионов ЫО+

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

- нагревание за счет фотоионизации нейтральных компонент,

- нагревание за счет фотодиссоциации 02 в континууме Шумана-Рунге,

- нагревание за счет фотоэлектронов,

- охлаждение за счет излучения О на длине волны 630 нм,

- нагревание за счет диссоциативной рекомбинации (химическая энергия),

- нагревание за счет диссипации турбулентной энергии,

- охлаждение за счет инфракрасного излучения молекул СО2 на длине волны 150 нм ()1к,

- нагрев за счет рекомбинации кислородных составляющих

- охлаждение за счет излучения N0 на длине волны 53 нм О^о

На уровнях верхней мезосферы и нижней термосферы, где чрезвычайно велика роль процессов переноса, таких, как вертикальные упорядоченные макроскопические движения, турбулентность, формируется качественный и количественный состав вышележащих слоев нейтральной атмосферы, поскольку распределение концентраций нейтральных компонент выше примерно 140 км в большей степени подчиняется условию диффузионного равновесия Развитие современных средств экспериментального исследования параметров околоземного космического пространства, успехи лабораторного изучения элементарных процессов, проведение комплексных экспериментов позволили получить более надежные сведения о верхней атмосфере Земли Накопление

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

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

скоростей движения нейтрального состава использовалось уравнение движения нейтрального газа, в результате численного решения которого показано, что в уравнении для вертикальной составляющей скорости

дР

нейтрального ветра члены — и g являются определяющими и

дг

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

— в области нижней термосферы основным источником нагрева является солнечное излучение в континууме Шумана - Рунге,

— несовпадение высот максимума выхолаживания за счет турбулентного перемешивания и максимума нагрева за счет континуума Шумана-Рунге приводит к инверсии температуры в области турбопаузы

— в области термосферы источник нагрева и охлаждения -фотохимические реакции

В главе 4 приводятся варианты разностных аппроксимаций моделирующих уравнений и методы их решения

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

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

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

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

Для проверки качества прогноза распространения КВ-радиосигналов необходимо

- спрогнозировать состояние ионосферы динамической моделью,

- рассчитать трассы,

- результаты по расчету трасс сравнить с экспериментальными данными

В главе были рассмотрены следующие вопросы

- рассмотрены различные варианты решения уравнения непрерывности,

- построены различные варианты потокового метода прогонки для уравнений непрерывности и теплопроводности;

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

- для реализации этих методов был составлен алгоритм и написана ФОРТР АН-программа

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

В заключении приведены основные результаты, полученные в работе

Заключение

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

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

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

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

• разработаны алгоритмы потокового варианта метода прогонки для уравнения теплопроводности

3 Проведены вычислительные эксперименты по самосогласованному расчету температуры нейтрального газа, ионного и нейтрального состава для средних широт (45° К), средней солнечной активности (Г^ 7 =150), с шагом интегрирования по высоте А=2 км и времени т =5 мин , в результате которого

- выявлены основные источники нагрева и охлаждения мезосферы и нижней термосферы,

- показана инверсия температуры в области турбопаузы, объяснена природа этой инверсии,

- получено новое аналитическое распределение окиси азота в области мезосферы,

- рассчитаны высотно-временные распределения электронной концентрации в области D ионосферы, удовлетворительно совпадающие с последними экспериментальными данными

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

Список опубликованных работ по теме диссертации

1 Игианов С А , Медведев В В, Залесская В А , Жаркова Ю С Математическое моделирование ионосферных процессов в целях распространения радиоволн//Математическое моделирование, 2008 Т 20, №4, С 3-7 (опубликована в журнале, включенном в перечень ведущих рецензируемых научных журналов и изданий, в которых должны быть опубликованы основные научные результаты диссертации на соискание ученой степени кандидата наук)

2 Ишанов С А , Медведев В В, Захаров Л П, Залесская В А, Жаркова Ю С Эффекты возмущения нейтральных ветров//Вестник Калининградского государственного университета Калининград Серия информатика и телекоммуникации, 2005 № 1-2 С 54-59

3 Ishanov S А , Medvedev V V, Zalesskaya V А , Zharkova YS Mathematical model of the metastable species in the ionosphere and

thermosphere/Y№6paHHbie вопросы современной математики Издательство КГУ, 2005. С 139-140

4 Ishanov SA, Medvedev V V, Tokar VG, Zharkova YS Possibility of applying the mathematical models to the problems of radiowave propagation//Physics of Auroral Phenomena 29th Annual Seminar Polar Geophysical Institute Apatity, 2006 P62

5 Medvedev V V, Zharkova YS Mathematical modeling for the processes upper atmosphere and ionosphere//Physics of Auroral Phenomena 29th Annual Seminar Polar Geophysical Institute Apatity, 2006 P 65

6 Medvedev V V, Zharkova YS Numerical model of the heat budget of the Earth's upper atmosphere//Physics of Auroral Phenomena 29th Annual Seminar Polar Geophysical Institute Apatity, 2006 P 66

7. Ishanov S A , Medvedev V V, Zharkova YS Excited components m the magnetosphere and the ionosphere with the powerful ionosphere disturbances//6-th International conference Problems of Geocosmos SamtPetersburg, 2006 P 207

8 Medvedev V V, Zharkova YS Dynamical and photochemical heating and cooling algorithm used m a dynamical model of the upper atmosphere Earth//6-th International conference Problems of Geocosmos Saint-Petersburg, 2006.

P 244-245

9 Ишанов С A , Леванов E И, Медведев В В, Залесская В А , Жаркова Ю С Магнитосферно-ионосферные изменения, вызванные полетами космических аппаратов//Инженерно-физический журнал Минск, Беларусь, 2006гТ 79 №6 С 11-15

10 Жаркова Ю С, Ишанов С А , Медведев В В, Токарь В Г Использование математических моделей ионосферы для изучения распространения КВ-радиотрасс (радиосигналов)//Вестник РГУ им И Канта Калининград, 2006 №10 С 49-54

11 Ишанов С А, Медведев ВВ, Новикова ЕИ, Жаркова ЮС Влияние магнитосферно-ионосферных потоков плазмы на F-область ионосферы// Вестник РГУ им И Канта, Калининград, 2007, №10 С 15-19

Жаркова Юлия Сергеевна

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛОВОГО БАЛАНСА, НЕЙТРАЛЬНОГО И ИОННОГО СОСТАВА МЕЗОСФЕРЫ И НИЖНЕЙ ТЕРМОСФЕРЫ

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Подписано в печать 13.03.08 г. Формат 60x90

Бумага для множительных аппаратов. Ризограф. Усл.печ.л. 0,8. Уч.-изд. л. 0,7. Тираж 100 экз. Заказ 33 .

Издательство Российского государственного университета имени Иммануила Канта 236041, г. Калининград, ул. А. Невского, 14

Введение.

Глава 1. Экспериментальные данные и система гидродинамических уравнений.

1.1. Экспериментальные данные по составу средней и верхней атмосферы.

1.2. Основные физико-химические процессы и модельные уравнения.

1.3. Характерные времена химических процессов.

1.4. Времена диффузии,.

1.5. Времена среднемассового переноса.

1.6. Уравнения гидродинамики в одномерном приближении.

1.7. Краевые условия, начальные условия.

Глава 2. Уравнения непрерывности и фотохимическая схема

2.1. Уравнения непрерывности для нейтральных и заряженных частиц.

2.2. Расчеты высотного распределения [О] и [О2].

2-3. Расчеты высотного распределения [А^(45')], [N0].

2.4. Расчеты высотного распределения [Ж)] с использованием связанных краевых условий.

Глава З.Уравнения теплопроводности и движения нейтрального газа

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

3.2. Уравнение теплового баланса нейтрального газа.

3.3. Результаты численных расчётов.

3.4. Макроскопические скорости движения нейтрального газа.

Глава 4. Варианты разностной аппроксимации моделирующих уравнений и численные методы их решения.

4.1. Численные методы решения уравнения непрерывности.

4.2. Граничные и начальные условия.

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

4.4. Потоковый вариант метода прогонки для уравнения теплового баланса.

4.5. Численные методы решения уравнения движения нейтрального газа.

4.6. Построение консервативных разностных схем для уравнений теплопроводности и диффузии.

4.7. Использование модели в целях распространения радиотрасс.

Как известно, область высот 50—250 км труднодоступна для экспериментальных исследований. В настоящее время разработаны различные методики экспериментов, проведение которых несет эпизодический характер и дает различные ошибки измерений. Несмотря на значительные экспериментальные успехи, достигнутые в последнее время в области исследования верхней атмосферы, математическое моделирование остается основным (достаточно дешевым) методом исследования этой области высот.

В теоретическом плане изучение этой области высот затруднено необходимостью учета сложных динамических и фотохимических процессов, таких, как турбулентное перемешивание, переходящее в молекулярную диффузию, поглощение нейтральным составом солнечного излучения и его переизлучение, большая плотность и многокомпонентность состава, наличие «малые компонент»: О, Оз, С02, 02('АН20, N0, концентрация которых существенно меньше основных N2, О2, но которые могут играть существенную роль как в тепловом балансе, так и в образовании ионосферы.

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

Цель работы:

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

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

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

4. Исследовать возможности применения результатов модели в задачах распространения электромагнитных волн.

Метод исследования основан на численном решении систем нелинейных связанных жестких уравнений гидродинамики (непрерывности, движения, теплопроводности).

Научная новизна работы:

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

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

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

4. Методом вычислительного эксперимента получена инверсия температуры в области турбопаузы и дано ее возможное теоретическое объяснение.

5. Получено новое аналитическое распределение окиси азота в области мезосферы.

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

Диссертация состоит из введения, четырех глав, заключения, библиографического списка. Главы делятся на пункты. Нумерация формул и рисунков в диссертации проводится по главам.

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

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

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

-разработаны алгоритмы потокового варианта метода прогонки для уравнения теплопроводности.

3. Проведены вычислительные эксперименты по самосогласованному расчету температуры нейтрального газа, ионного и нейтрального состава для средних широт (45° IV), средней солнечной активности (^0.7 =150), с шагом интегрирования по высоте к=2 км и времени т =5 мин., в результате которого: выявлены основные источники нагрева и охлаждения мезосферы и нижней термосферы; показана инверсия температуры в области турбопаузы, объяснена природа этой инверсии; получено новое аналитическое распределение окиси азота в области мезосферы; рассчитаны высотно-временные распределения электронной концентрации в области О ионосферы, удовлетворительно совпадающие с последними экспериментальными данными.

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

124

ЗАКЛЮЧЕНИЕ

1. Shimazaki Т., Layrd A.R. Seasonal effect on distributions of minor neutral constituents in the mesosphere and lower thermosphere.// Radio Sci., 1972, vol. 7,1, p. 23-43.

2. Philbrick G.R., Narcisi R.S., Good R.E., Hoffman H.S., Keneshea T.J., McLeod M.A., Zimmerman S.P., Reinisch B. W. The Alladin experiment — part 2, composition.//Space Res., 1973, vol. 13, p. 441-448.

3. Offerman D. Some resulte from the european winter anomaly campaign 1975/76. In: Dynamical and chemical coupling between the Neutral and Ionized Atmosphere, Reidel, Dordrecht, Holland, 1977, p. 235-252.

4. Good R.E., Golomb D. Atomic oxigen profiles in the lower thermosphere.//Space Res., 1973, vol. 13, p. 249-253.

5. Trinks H., Offerman D., Zahn U., Steinhauer C. Neutral composition measurements between 90 and 220 km altitude by rocketborn mass-spectrometr.//J. Geoph. Res., 1978, vol. 83, №5, p. 2169-2176.

6. Колесник А.Г. Полугодовые вариации нейтрального состава основания термосферы Земли.// Геомагнетизм и аэрономия, 1975, т.15, № 2, с. 286-291.

7. Кошелев В.В., Климов Н.Н., Сутырин НА. Аэрономия мезосферы и нижней термосферы.//М.,Наука, 1983, с. 112-124.

8. Соломонов С.В., Кропоткина Е.П., Лукин А.Н., Пономаренко Н.И., Розанов С.Б. О вариациях атмосферного озона по наземным измерениям на миллиметровых волнах.// Изв. АН. Физика атмосферы и океана, 1993, т. 29, №8, с. 525-531.

9. Гайкович К.П., Кропоткина Е.П., Соломонов С.В. Определение вертикального профиля атмосферного озона по наземным измерениям излучения в миллиметровом диапазоне.// Изв. АН. Физика атмосферы и океана, 1999, т. 35, №1, с. 86-95.

10. Брасъе Г., Соломон С. Аэрономия средней атмосферы. — JI. Гидрометеоиздат, 1987, с. 414.

11. Семенов А.И., Шефов H.H. Модель вертикального распределения концентрации атомарного кислорода в области мезопаузы и нижней термосферы.// Геомагнетизм и аэрономия, 2005, т. 45, №6, с. 844-855.

12. Barth С А. Rocket measurement of nitric oxyde in the upper atmosphere.// Planet. Space Sei., 1966; v. 14, N. 7, p. 623-630.

13. Tohmatsu Т., Iwagami N. Measurement of nitric oxide abundance in equatorial upper atmosphere.// J. Geophys. Res., 1976, v. 28, N. 5, p. 343-358.

14. Baker K.D., Nagy A.F, Olsen R.O., Oran E.S., Strobel D.F., Tohmatsu T. Measurement of the nitric oxide altitude distribution in the mid-latitude mesosphere.// J. Geophys. Res., 1977, v. 82, N. 22, p. 3281-3286.

15. Beran D., Bangert W. Trace constituents in the mesosphere and lower thermosphere during winter anomaly events.// J. Atmoph. Terr. Phys., 1979, v. 41, N. 10/11, p. 1091-1095.

16. Feldman P.D., Takats P.Z. Nitric oxide gamma and delta band emission at twilight.// J. Geophys. Res. Lett., 1974, v. 1, N. 4, p. 169-171.

17. Gerard J. С. Satellite observations of the nitric oxide.// J. Geophys. Res., 1975, v. 2, N. 5> p. 179-182.

18. Takats P.Z., Feldman P.D. Far ultraviolet atomic and molecular nitrigen emission in the daygiow// J. Geophys. Res., 1977, v. 82, N. 32, p. 5011-5023.

19. Rusch D.W., Stewart Ä.I., Hays P.В., Hoffman J.H. The N1(5500 A) Daygiow.// J. Geophys. Res., 1976, v.81, N. 1, p. 295.

20. Roble R.G., Stewart A.I., Torr M.R., Rusch D.W., Wand R.H. The calculated and observed ionospheric properties during Atmospheric Explorer-C satellite crossing jver Millstone Hill.// J. Atmosph. Terr. Phys., 1978, v. 40, N. 1,p. 21-23.

21. Narcisi R.S., Beiley A.D.,Delia Lucca L., Sherman C., Thomas D.MJI J. Atmos. Terr. Phys, 1971, v. 33, p. 1147.

22. Arnold F., Kissel J., Krankovsky D., Wieder H., Zahringer J Л J. Atmos. Terr. Phys., 1971, v. 33, p. 1169.

23. Rowe J.N., Mitra A.P., Ferraro A.J., Les H.S. An experimental and theoretical study of D-regiomll. A semiepocal model for midlatitude D-region.// J. Atmos. Terr. Phys., 1974, v. 36, № 5, p. 755-785.

24. Власов M.H., Медведев.В.В. О возможном механизме образования NO и N в нижней термосфере и мезосфере.// Геомагнетизм и аэрономия, 1981,т. 21, №5, с. 857.

25. Власов М.Н., Медведев В.В. Анализ источников ионизации в области D на основе теоретической модели.// Геомагнетизм и аэрономия, 1984, т. 24, №2, с. 187.

26. Медведев В.В., Ишанов С.А. Зенкин В.И. Самосогласованная модель нижней ионосферы.//Геомагнетизм и аэрономия, 2002, т. 42, № 6, с. 780-789.

27. Каширин А.И. Фотоионизация в ночной ионосфере//ГиА, 1986, т. 26, №4, с. 563-569

28. Фейгин A.M. Нелинейно-динамические модели атмосферных фотохимических схем: методы построения и анализа// Изв. РАН. Физика атмосферы и океана, 2002, т. 38, №5, с. 521-628

29. Каролъ И.Л., Киселев А.А. Фотохимические модели атмосферы и. их использование в исследованиях озоносферы и климата (обзор)// Изв. РАН. Физика атмосферы и океана, 2006, т. 42, №1, с. 2-34.

30. Мак-Ивен М., Филлипс Л. Химия атмосферы. М., Мир, 1978.

31. Shimazaki Т., Laird А.г. A model calculation of the diurnal variation in minor neutral constituents in the including transport effect.// J. Geophys. Res., 1970, v. 75, № 16, p. 3221-3235.

32. Strobel D.F. Diurnal variations of nitric oxide in the upper Atmosphere.// J. Geophys. Res., 1971, v. 76, № 10, p. 2441-2452.

33. Norton R.B., Barth С.A. Theory of nitric oxide in earth's atmosphere.// J. Atmos. Terr. Phys., 1970, v. 75, № 18, p. 3903-3909

34. Ogawa Т., Shimazaki T. Diurnal Variations of Odd Nitrogen and Ionic Densities in the Mesosphere and Lower Thermosphere: Simultaneous Solution of Photochemical-Diffusive Equations // J. Geophys. Res., 1975, v. 80, N. 28, p. 3945 -3958.

35. Медведев B.B., Зенкин В.И. Возможная роль NO и 02('Лё) в образовании зимней аномалии области D ионосферы.//Геомагнетизм и аэрономия, 1998, т. 38, № 2, с. 156-160.

36. Kondo J., Ogawa J. A temperature-dependent model of the thermospheric odd nitrogen//J:Geomagn. Geoelectr, 1977, v. 29, p. 65-80.

37. Данилов АД. Ледомская С.Ю. Роль фотохимии и динамики в области D ионосферы.//Ионосферные исследования, 1982, № 32, с. 78-99.

38. Solar-Terrestrial Energy program: Handbook of Ionospheric models. Edited by R.W.Schrunk, August, 1996, p. 206.

39. Семенов А.И., Шефов H.H. Вариации температуры и содержания атомарного кислорода в области мезопаузы и нижней термосферы при изменении солнечной активности //Геомагнетизм и аэрономия, 1999, т. 39, № 4, с. 87-91.

40. Перминов В.И., Кропоткина Е.П., Баканас В.В., Перцев Н.Н., Соломонов С.В, Розанов С.Б., Лукин А.Н. Определение концентраций основных и малых газовых компонентов атмосферы на высотах мезопаузы// Геомагнетизм и аэрономия, 2002, т. 42, № 6, с. 814-820.

41. Shved G. М:, Ogibalov V.P., Pogoreltsev A.I. Effect of planetary waves on cooling the upper mesosphere and lower thermosphere by the C02 15-цт emission// Annales Geophysical; 2004, № 22, p. 3383-3394.

42. Medvedev V. V., Zharkova Y.S. Numerical model of the heat budget of the Earth's upper atmosphere// Physics of Auroral Phenomena 29th Annual Seminar Polar Geophysical Institute, Apatity, 2006, p. 66.

43. Лысенко E. В., Перов С. П, Семенов А. И., Шефов Н. Н, Суходоев В. А., Гивишвши Г. В., Лещенко Л. Н. Многолетние тренды среднегодовой температуры на высотах 25-110 км// Известия Ан. Физика атмосферы и океана, 1999, т. 35, №4, с. 435-443

44. Фаделъ X. М, Семенов А. И., Шефов Н. Н., Суходоев В. А., Марцвсшадзе Н. М. Квазидвухлетние вариации температуры мезопаузы, нижней термосферы и солнечной активности// Геомагнетизм и аэрономия, 2002, т. 42, № 2, с. 203-207

45. Семенов А. О., Швед Г. М. Полуэмпирическая модель среднеглобальной температурной структуры земной термосферы для переменного содержания углекислого газа.//Известия Ран. Физика атмосферы и океана, 2004, т. 40, № 3, с. 291-305

46. Косцов В. С., Тимофеев Ю. М. Мезосферные инверсии температуры по данным эксперимента CRISTA-1.//Известия Ран. Физика атмосферы и океана, 2005, т. 41, №6, с. 814-823

47. Самарский А.А., Вабишевич П.Н. Численные методы решения задач конвекции-диффузии. М.: Едиториал УРСС. - 2003.

48. Петров И. Б., Лобанов А. И. Лекции по вычислительной математике.//М., 2006, с. 522

49. Дегтярёв А.П., Фаворский А.П. Потоковый вариант метода прогонки для разностных задач с сильно меняющимися коэффициентами.// Ж. вычисл. математики и мат. физики, 1968, т. 8, № 3, с. 211-218.

50. Латышев К. С., Медведев В.В. Варианты метода прогонки численного решения уравнений диффузии ионов в задаче моделирования ионосферы. Диагностика и моделирование ионосферных возмущений.//М. Наука, 1978,с. 108-114.

51. Медведев В.В., Латышев К.С., Коняхина Л.В. Описание модели DEF.// Материалы мирового.центра данных Б, выпуск 8, 1989, с. 8-41.

52. Медведев В.В., Зенкин В.И. Численные методы решения уравнений диффузии.// Математическое моделирование и численные методы решения интегральных уравнений. Калининград. Из-во КГТУ, 2003, с. 72-78.

53. Жаркова Ю.С., Ишанов С.А., Медведев В.В., Токарь В.Г. Использование математических моделей ионосферы в целях распространения КВ-радиотрасс (радиосигналов).// Вестник РГУ им. И. Канта, 2006, №10, с. 4954.

54. Кравцов Ю.А., Орлов Ю.И. Геометрическая оптика неоднородных сред.// М.:Наука, 1980.

55. Лукин Д. С., Спиридонов'Ю.Г. Применение метода характеристик для решения на ЭВМ задач распространения электромагнитных волн в неоднородных анизотропных средах.//сб. Лучевое приближение и вопросы распространения радиоволн. М.:Наука, 1971.

56. Токарь В.Г. Надежников Ю.И., Никитин МАРасчет параметров КВ-радиосигнолов с использованием табличных моделей ионосферы.// Геомагнетизм и аэрономия, 1980, № 4, с. 753-756.

57. Токарь В.Г., Рубинштейн Л.И., Никитин М.А. Исследование деполяризации КВ радиоволн в ионосфере в квазиизотропном приближении.//Известия Вузов-Радиофизика, 1987, т. XXX, № 1, с. 36-41.

58. Крашенинников И.В., Егоров И.Б., Коломийцев О.П., Черкашин Ю.Н. Погрешности прогнозирования ионосферного происхождения радиоволн на основе глобальной ионосферной модели.//Геомагнетизм и аэрономия, 2004,т. 44, №2, с. 221-226.

59. Жаркова Ю.С., Ишанов СЛ., Медведев В.В., Токарь В.Г. Использование математических моделей ионосферы для изучения распространения КВ-радиотрасс (радиосигналов).//Вестник РГУ им. И. Канта, Калининград, 2006, №10, с. 49-54.

60. Ишанов СЛ., Медведев В.В., Залесская В.А., Жаркова Ю.С. Математическое моделирование ионосферных процессов в целях распространения радиоволн.//Математическое моделирование, 2008, т. 20, №4, с. 3-7

61. Ишанов СЛ., Медведев В.В., Новикова Е.И., Жаркова Ю.С. Влияние магнитосферно-ионосферных потоков плазмы на Р-область ионосферы.// Вестник РГУ им. И. Канта, Калининград, 2007, №10, с. 15-19.