Введение к работе
Актуальность проблемы и предмет исследования. Анализ эволюции представлении о различных наблюдаемых астрофизических объектах показывает, что очень часто такие объекты, воспринимаемые изначально как сугубо вещественные, материальные образования, впоследствии оказывались результатом визуализации волновых структур. В этом плане чрезвычайно актуальным для современной астрофизики является исследование коллективных процессов, протекающих в образованной звездами галактик плазме гравитирующих частиц и в газовых подсистемах астрофизических объектов, и приводящих к возникновению таких регулярных структур. Важное место в этих исследованиях должно занимать выявление различного рода ре-зонансов, поскольку они позволяют осуществлять эффективную передачу гравитационной пли кинетической энергии неволновых движений звезд или газа в энергию волны без дополнительной подпитки энергией извне, обуславливая тем самым колебательные неустойчивости, лпбо автоколебательные процессы в исследуемых объектах.
Если важная роль динамических резонансов и их влияние на орбиты звезд в дисках плоских галактик исследованы уже достаточно полно, то применительно к резонансным явлениям в газовых и газопылевых подсистемах астрофизических объектов можно сказать, что их изучение пока только начинается. В то же время целый ряд успешных исследований убедительно показывает перспективность такого подхода. Примером прогресса, достигнутого благодаря исследованию резонансных эффектов, является предсказание (с погрешностью менее 0.5%) Н.Н. Горькавым и A.M. Фридманом орбит резонансных спутников за внешней границей колец Урана, сделанное на основе разработанной ими замкнутой теории коллективных процессов в кольцах планет [1]. В качестве второго примера можно привести работы [2-3], где показана возможность развития большого числа резонансных неустойчивых спиральных волн в газовых и плазменных дисках аккрецирующего вещества, чем может быть объяснен наблюдаемый широкий спектр нестационарностей аккреционных дисков и интенсивная турбулпзацпя вещества в нпх.
В настоящей работе рассмотрение динамики газовых подсистем астрофизических объектов самых разных масштабов — нижней хромосферы Солнца, аккреционно-струпных течений вокруг звезд на ранней стадии эволюции, аккреционных дисков в тесных двойных сп-
стемах, дисков плоских галактик, — проводится с позиций возможности согласования волновых процессов в различных областях этих подсистем резонансными эффектами.
Газовые подсистемы астрофизических объектов характеризуются значительными градиентами скоростей и термодинамических параметров газа. Наличие в них сверхзвуковых перепадов скорости делает возможным развитие неустойчпвостей, имеющих характер акустического резонанса. В то же время отсутствие твердых границ и крайне малая вязкость открывают обширное поле деятельности для теоретика, позволяя значительно упростить модель и исследовать механизмы раскачки неустойчивостей "в чистом виде", чем астрофизические задачи существенно отличаются от прикладных, в которых теплопроводность границ и вязкие эффекты сильно осложняют анализ устойчивости сверхзвуковых течений.
Подробное теоретическое исследование резонансных неустойчивостей сверхзвуковых потоков газа начато относительно недавно — в конце 70-х годов [4-6], несмотря на то, что наблюдаемые проявления этих неустойчпвостей известны уже более века. Предшествующие разработки носили сугубо прикладной, инженерный характер и, в соответствии с запросами аэродинамики, относились к течениям с характерными числами Маха М ~ 1.5 -f 3. Основным отличием инженерных задач, решаемых применительно к струям из сопел реактивных двигателей, является истечение струи из области с давлением, существенно превышающим давление окружающего газа.
К концу 70-х годов наблюдательная астрономия предоставила теоретикам великолепную модель для исследования резонансных неустойчивостей — струйные выбросы из радиогалактик, проэволюцп-онировавших двойных звездных систем и молодых звезд на стадии эволюции до главной последовательности. Струйные выбросы из молодых и проэволюционировавших двойных звездных систем являются, по современным представлениям, скорее правилом, чем исключением, и характеризуются отсутствием контраста концентраций между веществом струп и внешней средой (п{п/пех ~ 1), числом Маха М = Vin/c{„ ~ 10 -f- 40, отношением длины к диаметру L ~ 10 4- 30 и очень высокой степенью коллимации — угол раствора составляет а ~ 0 4-5 [7-8]. Анализ показывает, что эти джеты коллимпруются внешним давлением (т.е. не являются баллистическими), поэтому для исследования их устойчивости применимы баротропно равновесные модели,
чем джеты принципиально отличаются от струйных выбросов реактивных двигателей. Исследования устойчивости джетов в моделях цилиндрических и плоских струй с разрывными границами показали, что существенно сверхзвуковой характер течения приводит к появлению нового типа неустойчивости. Наряду с основными модами, имеющими винтовую и пинчевую геометрию, присутствующими и в дозвуковых струях и обусловленными механизмом неустойчивости Кельвина-Гельмгольца, в этом случае возбуждается большое число высших неустойчивых гармоник этих мод, называемых отражательными. Эти гармоники по существу являются неустойчивыми собственными модами волноводного слоя, образованного границами струи, и различаются числом узлов собственных функций между этими границами. Механизм неустойчивости отражательных гармоник подробно анализировался в [9-12]. В его основе лежит открытый в 1957 г. Майлсом и Рибнером эффект сверхотраженпя [13-14]: звуковая волна, падающая на поверхность тангенциального разрыва с перепадом скорости AV > с\ + С2, где сі и С2 — скорости звука по разные стороны от разрыва, может отражаться от него с усилением по амплитуде. При этом существуют углы падения, называемые резонансными, для которых коэффициент отражения стремится к бесконечности из-за того, что в прошедшей волне поток энергии направлен к разрыву. Последнее означает, что сверхзвуковой разрыв спонтанно излучает звуковые волны. Если испытавшая сверхотраженпе волна имеет возможность возвращаться к разрыву, т.е. если в потоке присутствует параллельная разрыву отражающая поверхность, то энергия такой волпы будет нарастать во времени в волноводном слое между разрывом и акустическим экраном, что и представляет собой неустойчивость.
Таким образом, исследование гидродинамических резонансных неустойчивостей в сверхзвуковых потоках газа начато сравнительно недавно и в последние два десятилетия находится в стадии бурного развития. В астрофизике это связано с тем, что, как оказалось, раскачка возмущений этими неустойчивостями до нелинейных амплитуд теоретически способна приводить к крупномасштабным волновым структурам, сходным с наблюдаемыми в струйных выбросах из молодых звезд и двойных звездных систем, в газовых галактических и в аккреционных дисках. Кроме того, данные неустойчивости являются многомодовыми, поддерживаемые ими высшие гармоники образуют иерархический ряд уменьшающихся пространственных масштабов и
возрастающих частот, что способно приводить к эффективной турбу-лизацип течения. Последнее может оказаться важным применительно к объяснению высокой турбулентной вязкости вещества аккреционных дисков.
Несмотря на значительное количество появившихся в последнее время публикаций по данному кругу вопросов, остается ряд нерешенных проблем, и реферируемая работа должна внести существенный вклад в их прояснение.
Цепь работы. Основной целью реферируемой диссертации являлось изучение принципиальной возможности и характерных особенностей развития резонансных эффектов и, в частности, волноводно-резононсных неустойчивостей в астрофизических струях, в аккреционных дисках и в газовых подсистемах дисков плоских галактик, всестороннее исследование их физических механизмов и параметров генерируемых ими волновых структур, выявление возможности согласования процессов, протекающих в различных газовых подсистемах, из-за резонансного обмена энергией и угловым моментом, построение последовательной теории аккреционно-струйных систем, находящихся в поле гравитационного потенциала центрального тела или распределенной массы, и объяснение наблюдаемых феноменов на основе проведенного анализа.
Научная новизна реферируемой работы состоит в том, что в ходе проведенных исследований были впервые получены следующие основные результаты:
разработана теория аккреционно-струйных систем, находящихся в гравитационном поле центрального объекта, позволящая исследовать причины формирования и эволюции крупномасштабной волновой структуры таких систем с учетом радиальных неодно-родностей их параметров;
предложено коротковолновое приближение, пригодное для описания динамики возмущений в системах, находящихся в центрально-симметричном гравитационном поле со степенной зависимостью от радиуса;
предложено объяснение синхронизации излучающих узлов в диаметрально противоположных выбросах из протозвездных систем неустойчивыми модами джетов посредством возбуждения ими конусов Маха в окружающей среде и волн в аккреционном диске;
найдена и исследована волноводно-резонансная неустойчивость
энтропийно-вихревых волн струп, находящейся в квадратичном гравитационном потенциале, отвечающем потенциалу звездного балджа;
предсказана теоретически и впервые показана численным моделированием закономерность образования широких конусов горячего вещества вокруг коллимированных выбросов (радиоджетов) из галактических ядер и формирования в таких конусах ярких волокнистых структур;
показано, что при наличии околозвездного диска одиночный сверхзвуковой выброс из протозвезды способен вызвать серию вторичных выбросов и в конечном итоге привести к формированию высококоллпмпрованной сверхзвуковой струп; определены физические причины этого эффекта;
найдена неустойчивость типа гироскопического резонанса в газовых дисках с резким изменением скорости вращения с радиусом;
предложена методика обработки результатов численного нелинейного моделирования динамики возмущений в газовых дисках, позволяющая анализировать физические причины п механизмы переноса массы и углового момента;
предсказан теоретически и впервые показан численным моделированием автоколебательный квазппериодпческий режим во внутренних областях галактических дисков, обусловленный резонансным взаимодействием спиральной волны плотностп в газовом диске с волной потенциала бара;
впервые предложено объяснение солнечных хромосферных спикул развитием коллективных процессов в плазме.
Научная и практическая значимость диссертационной работы состоит в том, что проведенный анализ возможности резонансно-волнового обмена энергией и угловым моментом и соответствующего взапмосогласованпя различных подсистем астрофизических объектов позволил объяснить ряд наблюдаемых в них феноменов, открыть принципиально новые режимы эволюции этих объектов и разработать теорию аккреционно-струйпых систем. Полученные в дпссертацпп результаты могут представлять интерес как с точки зрения фундаментальных исследований, так п с точки зрения применений для широкого круга специалистов, занимающихся проблемами астрономии, астрофизики и газодинамики в таких научных учреждениях как Институт астрономии РАН, Астрономический институт СПбГУ,
Институт космических исследований РАН, Государственный астрономический институт им. Штернберга, Специальная астрофизическая обсерватория РАН, Московский, Санкт-Петербургский, Ростовский госуниверситеты и т.д. Отдельные параграфы диссертации могут быть включены в учебные курсы по астрофизике, гидродинамике и теории гидродинамической устойчивости. Результаты первой главы могут быть использованы при подготовке программ наблюдений конусов ионизованного вещества в галактиках с активными ядрами.
В целом полученные результаты развивают важное новое направление астрофизики — теорию резонансной самоорганизации в астрофизических объектах.
На защиту выносятся следующие основные результаты:
-
Разработана теория аккреционно-струйных систем, находящихся в гравитационном поле центрального объекта. Она позволяет исследовать причины формирования и эволюции крупномасштабной волновой структуры таких систем.
-
В рамках созданной теории показано, что в молодых звездных системах, включающих аккреционные диски и биполярные струйные выбросы, глобальную волновую структуру создают неустойчивые моды струй. Предложено объяснение синхронизации излучающих узлов в диаметрально противоположных струйных выбросах такими модами посредством возбуждения ими конусов Маха в атмосфере и волн в диске.
-
Развитие резонансной неустойчивости внутренних гравитационных волн на границе коллимированных выбросов (радиоджетов) из галактических ядер приводит к образованию широких конусов ионизованного вещества вокруг радиоджета и к формированию в таких конусах наблюдаемых ярких волокнистых структур. Динамическое охлаждение высвечиванием усиливает этот эффект.
-
.При наличии околозвездного диска, одиночный сверхзвуковой выброс из протозвезды способен вызвать серию вторичных выбросов и в конечном итоге привести к формированию высококолли-мпрованной сверхзвуковой струи.
-
Развитие в аккреционных дисках магнитоакустической или акустической неустойчивости волн, резонирующих между критическими слоями в перпендикулярном к плоскости симметрии диска направлении, может создать условия, ведущие к турбулизации вещества диска, необходимой для обеспечения наблюдаемого темпа
аккреции.
-
Нелинейным численным моделированием показано, что во внутренних областях галактических дисков возможен автоколебательный квазиперподпческип режим, обусловленный резонансным взаимодействием спиральной волны плотности в газовом диске с волной потенциала бара. При этом за один квазиперпод, равный двум периодам оборота бара (3 107 4-108 лет) галактика последовательно проходит морфологические стадии спиральной, кольцевой и 0-галактпкн.
-
Впервые предложено объяснение солнечных спикул развитием коллективных процессов в плазме. Резонансное взаимодействие альфвеновскпх и внутренних гравитационных волн и развитие неустойчивости Кельвпна-Гельмгольца на границе ячейки суперкон-векцпп, либо полутени пятна в переходном слое от фотосферы к хромосфере естественным образом объясняет наблюдаемый спектр нестацпонарностей нижней хромосферы Солнца и фпла-ментарную структуру полутени солнечных пятен.
Достоверность результатов и выводов диссертации определяется физической обоснованностью используемых моделей и применением при решении поставленных задач строгих математических методов, проверкой согласования полученных в работе приближенных аналитических асимптотических решений с точными численными решениями в широких диапазонах значений параметров, сравнением результатов линейного анализа с результатами нелинейного моделирования, тщательным тестированием применяемых для этого моделирования конечно-разностных схем, а также совпадением в частных и предельных случаях полученных результатов с известными ранее и с данными наблюдений реальных объектов.
Публикации. Основное содержание диссертации изложено в 32 научных публикациях, список которых приведен в конце автореферата.
Апробация работы. Материалы настоящей диссертации докладывались на Всесоюзных научных конференциях "Классическая гравпфпзпка" (Волгоград, сентябрь 1989 г.), "Астрофизика сегодня" (Нижний Новгород, март 1991 г.), в рамках научной программы 11-го съезда Астрономического общества СССР (Москва, ноябрь 1991 г.), на Всесоюзном научном семинаре "Астрофизика - IV" рабочей
группы "Физика галактик" (Цейское ущелье, сентябрь 1990 г.), Всероссийском научном семинаре рабочей группы МЗС "Явления самоорганизации в галактиках" (Пулково, ГАО, сентябрь 1992 г.), международных научных конференциях "Physics of Gaseous and Stellar Disks of Galaxy" (Нижний Архыо, CAO, сентябрь 1993 г.), "Structure and Evolution of Stellar Systems" (г. Петрозаводск, август 1995 г.), "Современные проблемы астрофизики" (Москва, ноябрь 1996 г.), в рамках научной программы IV-ro съезда Астрономического общества (Москва, ноябрь 1997 г.), на Всесоюзных студенческих научных конференциях "Физика Космоса" (Свердловская обл., Коуровская АО, февраль 1991, 1993-2000 гг.), ВНКСФ-1 и ВНКСФ-2 (г. Заречный, Свердловская обл., апрель 1993 и 1994 гг.), научных семинарах кафедры теоретической астрофизики СПбГУ и Астрономического института РАН (Санкт-Петербург, сентябрь 1992 г.), кафедры астрономии и геодезии Уральского госуниверситета (Екатеринбург, апрель 1993 г.), Специальной астрофизической обсерватории РАН (октябрь 1996 г., апрель 1997 г. п ноябрь 1998 г.), Государственного астрономического института им. П.К. Штернберга (октябрь 1998 г.) и кафедры теоретической физики Волгоградского госуниверситета в 1985-2000 гг.
Личный вклад автора. По материалам первой главы автором предложена идея исследования, проведены все постановки задач и осуществлялось общее руководство работой. Участие в получении численных результатов и в написании программ для ЭВМ — равноправное с С.С. Храповым и К.А. Левиным. Интерпретация результатов проведена преимущественно автором.
По материалам второй главы автором выполнены общая постановка задачи, интерпретация результатов и написание статей. Математическая постановка задачи проделана совместно с С.С. Храповым. В обработке результатов моделирования участие К.А. Левина, С.С. Храпова и автора равноправное.
В описанных в третьей главе исследованиях в п. 3.1 и п. 3.2 автору принадлежит равноправная с А.В. Хоперсковым, а в остальных пунктах — лидирующая роль в постановках задач. Автором получены все аналитические результаты и примерно четвертая часть численных. Статья, составляющая содержание п. 3.1, написана А.В. Хоперсковым, в остальных пунктах интерпретация результатов п написание статей проделаны преимущественно автором.
В четвертой главе автором полностью выполнены постановки за-
дач, интерпретация результатов, написаны статьи, получены все аналитические результаты и примерно третья часть численных результатов.
В пятой главе автору принадлежат алгоритм обработки результатов численного моделирования, позволивший выявить квазиперно-дпческпй режим движения газа, интерпретация результатов и написание статей. В п.5.2 автору принадлежит постановка задачи и интерпретация результатов. В п.5.3 автор вместе с Е.А. Михайловой и В.В. Леви принимал равное участие в постановке задачи.
В главе 6 идея исследования, интерпретация результатов и написание статей выполнены совместно с А.А. Соловьевым. Автору принадлежат математическая постановка задачи и численные расчеты.
Объем и структура работы. Диссертация состоит из введения, шести глав, заключения, одного приложения п списка литературы. Общий объем диссертации составляет 352 страницы машинописного текста, включая 126 рисунков, 2 таблицы и список литературы (271 наименование).