Турбулентность и разрывы в сложных гидродинамических течениях жидкости и плазмы тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

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

ПЕТРОСЯН Аракел Саркисович

ТУРБУЛЕНТНОСТЬ И РАЗРЫВЫ В СЛОЖНЫХ ГИДРОДИНАМИЧЕСКИХ ТЕЧЕНИЯХ ЖИДКОСТИ И ПЛАЗМЫ

01 04 02 - теоретическая физика

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

Москва 2009

003476611

Работа выполнена в Институте космических исследований Российской академии наук (ИКИ РАН)

Официальные оппоненты.

д ф -м н, профессор, Жмур Владимир Владимирович (Институт океанологии им П П Ширшова РАН),

чл -корр РАН, д ф -м н, профессор, Петров Олег Федорович (Институт теплофизики экстремальных состояний Объединенный институт высоких температур РАН)

д ф -м н, профессор, Соколов Дмитрий Дмитриевич (Московский государственный университет им М В Ломоносова)

Ведущая организация: Институт физики атмосферы им А М Обухова РАН (г Москва)

Защита состоится 11 ноября 2009 г в 11 ч 00 мин на заседании Диссертационного Совета Д 002 113 03 ИКИ РАН по адресу Москва, ул Профсоюзная 84/32, 2-й подъезд, конференц-зал

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

Ученый секретарь Диссертационного совета,

кандидат физ-мат наук с/^/^?"-- БуринскаяТМ

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

Актуальность темы

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

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

Сжимаемая магнитогидродинамическая (МГД) турбулентность является широко распространенным состоянием космической плазмы во многих астрофизических, гелиофизических, геофизических процессах Например, в аккреционных дисках МГД-турбулентность вызвана магнито-ротационной неустойчивостью Большинство турбулентных явлений в физике Солнца описываются в рамках МГД солнечный ветер,

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Наиболее подробную информацию о турбулентном течении жидкости можно получить с помощью прямого численного моделирования - DNS (Direct numerical simulation), которое заключается в численном решении полной нестационарной системы магнитогидродинамических уравнений При таком подходе разрешаются все масштабы движения заряженной жидкости Метод DNS не требует специальных замыканий для уравнений магнитной гидродинамики Прямой численный расчет МГД-турбулентности сталкивается с

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

В методе крупных вихрей - LES (Large Eddy Simulation) используется операция фильтрации для разложения характеристик турбулентного движения на крупномасштабную и мелкомасштабную части, что связано с изотропностью, однородностью и универсальностью мелких масштабов турбулентного движения Мелкомасштабное движение исключается из исходной системы уравнений движения с применением процедуры фильтрации и дальше их влияние на движение моделируется с использованием подсеточных моделей SGS (subgrid scale) (или, другое название, SFS (subfilter scale)), выраженных через отфильтрованные параметры турбулентных течений Крупномасштабное движение рассчитывается из решения отфильтрованных нестационарных уравнений магнитной гидродинамики LES является методом для моделирования течений с большими числами Рейнольдса, так как в методе крупных вихрей предполагается, что энергия переносится от больших масштабов к малым только внутри инерционного интервала, поэтому число степеней свободы будет меньше, чем в методе DNS, следовательно, LES требует значительно меньших вычислительных затрат по сравнению с DNS

Изначально метод крупных вихрей развивался для моделирования гидродинамической турбулентности нейтральной жидкости для изучения задач метеорологии и океанологии Большая часть работ была выполнена для несжимаемых течений Метод LES к течениям электропроводящей жидкости применялся и исследовался крайне мало Все предыдущие работы в этом направлении ограничивались только рассмотрением несжимаемой жидкости для решения инженерных задач Тем не менее, многие течения электропроводящей жидкости не могут быть описаны в рамках несжимаемой среды, а необходимо рассматривать уравнения сжимаемой жидкости в приближении политропного соотношения между плотностью и давлением или теплопроводящую жидкость с использованием уравнения энергии Применение LES для сжимаемой теплопроводящей МГД-жидкости значительно усложняется из-за того, что нужно решать уравнение энергии, в котором появляются дополнительные слагаемые из-за наличия магнитного поля К тому же, после фильтрации появляются добавочные подсеточные члены, требующие разработки теории для их параметризации

Исследование сжимаемой турбулентности, как в гидродинамике нейтральной жидкости, так и в магнитогидродинамике, - трудная задача, поскольку не существует аналитической или приближенной теории таких явлений Однако, несмотря на существенную роль сжимаемости в космической плазме, целый ряд наблюдений показывает воспроизведение колмогоровского спектра флуктуаций плотности Для интерпретации таких наблюдений была использована асимптотическая теория "почти несжимаемой" (nearly incompressible) среды, которая описывает флуктуации плотности в гидродинамике нейтрального и магнитного газа в режиме переноса пассивного скаляра Однако в силу ограничений метода прямого численного моделирования ранее не удавалось получить спектры плотности и кинетической энергии и показать их совпадения и реализацию пассивного режима для флуктуаций плотности в сжимаемой МГД-турбулентности на основе полной трехмерной системы уравнений магнитной гидродинамики

Цель работы

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

В соответствии с поставленной целью в работе решаются следующие задачи

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

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

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

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

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

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

6 Осуществить численное моделирование методом крупных вихрей затухающей сжимаемой МГД-турбулентности в политропном газе и теплопроводящей плазме Сравнить результаты моделирования метода крупных вихрей с результатами, полученными методом прямого численного моделирования Исследовать спектры флуктуаций плотности и энергии в сжимаемой магнитогидродинамической турбулентной плазме для политропного уравнения состояния Показать возможность существования режима слабосжимаемых турбулентных пульсаций, когда флуктуации плотности являются пассивным скаляром Исследовать свойства анизотропии сжимаемой турбулентности и намагниченности плазмы в слабосжимаемом режиме

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

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

потоков Найдены аналитические решения для уравнений мелкой воды над наклонной плоскостью и неоднородных потоков

Разработан новый численный алгоритм, позволяющий вести сквозной расчет гидродинамических течений над сложной подстилающей поверхностью без специального выделения зон обмеления для нестационарной многосвязной расчетной области Для расчета течений мелкой воды над произвольной поверхностью в присутствии силы Кориолиса разработан модернизированный метод Годунова, основанный на квазидвухслойном представлении В отличие от известных методов расчета течений вращающейся мелкой воды предложенный метод адаптируется к значению текущих параметров ^течения

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

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

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

Практическая и научная ценность работы

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

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

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

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

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

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

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

Результаты могут найти применение в исследованиях атмосферы Земли и других планет, ведущихся в ИКИ РАН, ИПМ РАН, ФИ РАН, ИО РАН, ИФЗ РАН, ИФА РАН, ИВМ РАН, ИНАС РАН, ИЗМИР РАН, ИАЭ им Курчатова, ИГ СО РАН, НИРФИ, МИФИ, МФТИ, МГУ, ИТЭС ОИВТ РАН

Обоснованность и достоверность полученных результатов

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

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

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

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

исследования слабосжимаемого режима МГД-турбулентности обеспечивается сравнением с имеющимися данными наблюдений и приближенными теориями

Положения, выносимые на защиту

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

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

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

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

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

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

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

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

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

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

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

конференции МФТИ «Современные проблемы фундаментальных и прикладных наук» - 1999, Международной конференции МСС-04 «Трансформация волн, когерентные структуры и турбулентность» Москва, 2004, XLVH научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук Часть VIII Физика и энергетика» Москва, 2004, XLVIII научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук Часть VIII Физика и энергетика» Москва, 2005, XVI Научная сессия Совета РАН по нелинейной динамике Москва, 2007, XVII Научная сессия Совета РАН по нелинейной динамике Москва, 2008, Международная конференция «Современные проблемы газовой и волновой динамики», Москва, 2009, Международная конференция «Потоки и структуры в жидкостях физика геосфер», Москва 2009, Kiev International Workshop on Non-linear and Turbulent Processes in Physics 1989, International Symposium "Generation of large-scale structures in continuous media", Perm-Moscow, 1990, 1-st Liquid Matter Conference, Lyon , 1990, 4-st European Turbulence Conference, Delft , 1992, IUTAM Symposium " Eddy Structures Identification in Free Turbulent Shear flows, Poitiers, 1992, COSPAR assembly, Hamburg, 1994, EGS Assembly, Hamburg, 1995, ESA workshop on MARS, Capri, 1995, ERCOFTAC workshop on "Compressible turbulence" Chatillon sous Bagneux, France, 1995, EGS Assembly, the Hague, 1996, EUROMECH colloquium, Flow Control, Berlin, 1997, EGS Assembly, Vienna, 1997, 2nd Int Symposium on Turbulence, Heat and Mass transfer, Delft, 1997, EUROMEH colloquium on Statistics and Dynamics of Vorticity, Marseille, 1997, European Fluid Mechanics Conference, Gettingen, 1997, EGS Assembly, Nice, 1998, International Conference on Air Pollution Modelling, Champs-sur-Marne, 1998, International Conference on Control in Fluid Dynamics, Cachan, France, 1998, EGS Assembly, Den Haag, 1999, International Conference " Systems Sciences 2000 ", Osnabruek, 2000, EGS Assembly, Nice,2000, 8-th European Turbulence Conference, Barcelona, 2000, EGS Assembly, Nice, 2001, Second International Conference on Air Pollution Modelling, Champs-sur-Marne, 2001, Climate Conference 2001, Utrecht, 2001, EUROMRCH colloquium, Strongly Coupled Multiphase Flows, Grenoble, 2001, EGS Assembly, Nice,2002, EGS-AGU-EGU Assembly, Nice, 2003, EGU Assembly, Nice, 2004, EGU Assembly, Vienna, 2005, EMS meeting, Utrecht, 2005, EGU Assembly, Vienna, 2006, NATO Advanced Research Workshop on "Atmospheric Boundary Layers Modelling and Applications for Environmental Security", Dubravnik, 2006, Workshop on Environmental fluid mechanics as elements in agro meteorological modeling, As, Norway, 2006, International Conference " Turbulence and Interactions", Porqueroles, France, 2006, International Conference "Mathematical Hydrodynamics", Moscow, 2006, IUTAM Symposium " Hamiltonian Dynamics, Vortex Structures, Turbulrnce", Moscow, 2006, EGU Assembly, Vienna, 2007, 8th International School/Simposium for Space Simulations

(ISSS-8), Kauai, USA, 2007, EMS meeting, SanLorenzo, Spain, 2007, International Sympoposium on " Environmental Hydraulics, Tempe, Arizona,

2007, EGU Assembly, Vienna, 2008, EUROMECH Colloquium on " Mixing of Coastal, Estuarine and Riverine Shallow Flows", Ancona, Italy, 2008, Lecture Series 2008-06 Atmospheric boundary layer flows in air pollution modelling Von Karman Instituten for fluid dynamics, Rhode Saint Genese - Belgium, ERCOFTAC workshop "Direct and Large Eddy Simulations 7", Trieste, Italy,

2008, International conference "Meso-scale Meteorology and Air Pollution" Odessa, Ukraine, 2008, EUROMEH Fluid Mechanics Conference, Manchester, UK, 2008, EMS meeting, Amsterdam 2008, EGU Assembly, Vienna 2009

Публикации по теме диссертации

Основные результаты работы опубликованы в 21 статьях в реферируемых российских и международных журналах, в 5 статьях, опубликованных в рецензируемых трудах российских и международных конференций Результаты работы представлены в 45 тезисах докладов российских и международных конференций

Личный вклад автора

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

Структура и объем диссертации

Диссертационная работа состоит из введения, шести глав, заключения и списка цитированной литературы Объем диссертации -371 страница Библиография включает 297 наименований Диссертация содержит 157 рисунков, 2 таблицы

Содержание работы

Во введении обсуждается современное состояние теоретических и численных исследований разрывных и турбулентных течений в

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

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

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

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

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

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

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

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

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

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

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

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

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

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

В разделе 2 4 предложена теория для течений мелкой воды на ступенчатой границе с учетом вертикальной неоднородности поля скорости, основанная на двухслойной аппроксимации уравнений Эйлера вблизи ступенчатой границы Идея квазидвухслойной модели состоит в разбиении потока жидкости или газа на два слоя, так что для верхнего, с параметрами потока /г2 и «2, влияние ступенчатой границы, высоты а косвенно, а для нижнего, с параметрами потока Л, и ип ступенчатая граница влияет непосредственно, останавливая перенос массы в этом слое Система уравнений для двухслойной мелкой воды выглядит следующим образом

Э к д —1 +—

Э/ дх

+—(/щ) = 0,

(Л2„2)=О,

от ох ^

±{!щ)+1 ^ )+8А=0г

д1 дх 2 дх

+^=О

Ступенчатая граница полностью тормозит нижний слой », = 0, при д = 0, />0 Условия Коши для системы (1) имеют вид /г, = Л*, /г, = Л, - Л", и] = и,, и2 = и, при х < 0, /г2 = Ьг, и2 - иг при х > 0

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

членом ^Л,— в третьем уравнении системы (1) и, выбирая временной

дх

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

нижний, пренебречь членом Следовательно, система (1)

дх

к И II, вычисляются

независимо от переменных Л2 и иг Таким образом, влияние нижнего слоя на верхний будет осуществляться за счет изменения начальных условий, в частности условия Л2 = /!,-/»* при х<0, определяемого изменением глубины нижнего слоя при торможении потока

В результате пара параметров А, и их связывается с параметрами Л, и г/, через начальный параметр /г' Следовательно, проблема нахождения значения Л* сводится к решению обратной задачи Дирихле

СИ ОЛ

с)/ ах 2

Л, = /г* при х<0, Г = 0, и, =и, при х < 0, ¿ = 0, г/, = 0 при л = 0, г>0, /г, = я при х = 0, ?>0

Решение системы (2) зависит от направления течения жидкости, то есть скорости и1

а) и, > 0, картина течения имеет вид отраженной влево ударной волны, слева от которой параметры потока Л*и н, >0, а справа жидкость покоится, то есть 1г = а и » = 0, и глубина запирания, определяемая влиянием ступенчатой границы, вычисляется следующим образом

Решение определяется непосредственно из алгебраической формулы (3), вытекающей из соотношений Гюгонио на гидравлическом прыжке, б) и,< 0, картина течения имеет вид уходящей влево волны

разрежения, слева от которой параметры потока /¡* и ¡^<0, и справа так же /г = а и и-0, следовательно, решение определяется из постоянства соответствующего инварианта Римана, описывающего волну разрежения

Выражение (4) в явном виде определяет глубину нижнего слоя, полностью останавливаемого влиянием ступенчатой границы

Принципиально отличная ситуация будет складываться в случае, если /г* лежит вне интервала (0,/;,] Физически это означает либо

равенство нулю глубины /<*= 0 и, следовательно, отсутствие нижнего слоя, а вместе с ним и жидкости и/или изменения высоты подстилающей поверхности, либо в случае превышения Л* > /;, соответствует полному торможению потока слева, то есть в такой ситуации фиктивная поверхность является непротекаемой границей для всей жидкости слева, и следовательно, необходимо принять Лг = 0 при г<0 Тогда значение Л, при г = 0 может быть найдено как решение задачи о взаимодействии с непротекаемой границей всей жидкости, находящейся слева

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

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

НТ = Н'+тх

X У

2

(5)

/ ХН[ +

и\и[

Н'+|

где г - шаг по времени, Х - шаг по пространству, Н[ - глубина жидкости, и[ - скорость жидкости Нижний индекс х - значение функции, отнесенное к центру масс ячейки с номером х Полунижние индексы л ±1/2 - значения величин на границе между ячейками с номерами л и а±1 соответственно Верхний индекс t обозначает номер шага по времени

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

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

имеют вид ——, —---—, где слагаемое гха описывает

2 2

работу перепада высот подстилающей поверхности над потоком жидкости Здесь г принимает значение либо 0 в случае отрицательного перепада вь.сот подстилающей поверхности, либо 0</<1 в случае положительного перепада Переменная г принимает значение равное 1, в случае, если //", определяемое соотношениями (3),(4) для соответствующей грани ячейки х, не превышает значения глубины внутри нее В противном случае значение I на соответствующей грани является отношением глубины, формирующейся при полном торможении потока на указанной грани, к соответствующему перепаду высот подстилающей поверхности Значения //^±|/2, и'ла/2 на гранях вычисляются путем решения соответствующих задач Римана, на основе квазидвухслойной модели

В разделе 2 5 даны определения решений типа простых волн и ударных волн, приведена идея построения аналитического решения задачи Римана как для классической мелкой воды, так и для случая ступенчатой границы, и показаны возможные конфигурации течения при наличии ступеньки Аналитическое решение задачи принципиально важно для анализа получаемых в разделе 2 6 численных решений

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

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

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

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

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

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

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

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

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

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

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

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

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

линия, численно на основе \Л/АР - метода - прерывная линия, сплошная черная линия - на основе предложенного квазидвухслойного метода

0 6^

(¡"^(з) 1, ЧЦз)

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

точках

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

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

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

Аналогичный механизм генерации вихревых структур исследуется в разделе 4 3 для жидкости с распределенными газовыми пузырьками В такой среде асимметрия напряжения Рейнольдса обеспечивается

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

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

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

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

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

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

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

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

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

В разделе 5 2 приведены исходные уравнения и сформулирован метод крупных вихрей для политропной плазмы Для упрощения уравнений, описывающих турбулентное МГД-движение с переменной плотностью, применялась фильтрация по Фавру (другое название -средневзвешенная фильтрация), для того чтобы избежать появления дополнительных подсеточных слагаемых / = р f / р В работе получены следующие отфильтрованные МГД-уравнения в безразмерном виде до Э рй

+ о. (6)

Э/ dXj v '

дрй д pr _ 1 . В2 1 -- дт"

dt дх^'' уМ] ,J Re " 2МЛ М^ ' '' Эх' К '

Э - 1 Э2Д dfn

dt

Э В,

(1ЪВ,-В/0-—-= (8)

dt dXj J ' J ' Re„, Эх; d\j

Эх.

= 0, (9)

где р - плотность, скорость в направлении _

2 1 Эи Эм

ач = - вязкий тензор напряжений, = - тензор

скорости деформации, В1 - магнитное поле в направлении М-динамический (молекулярный) коэффициент вязкости, 8а - символ

Кронекера В системе МГД уравнений (6) - (9) использовались следующие безразмерные параметры подобия Лс- гидродинамическое число Рейнольдса, Яеи- магнитное число Рейнольдса, М,- число Маха, альфвеновское (магнитное) число Маха

Вследствие фильтрации по Фавру системы уравнений МГД политропной плазмы появляются подсеточные слагаемые в правой

части уравнений (7) и (8) г," = /?((»_,/',)' - нд)--^{В.В; - В)В,) -

"'л

подсеточный тензор напряжений и ^ ~-

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

В разделе 5 3 обобщены различные модели замыкания подсеточных слагаемых на случай сжимаемой МГД-турбулентности, а именно модель Смагоринского для МГД-турбулентности, модель Колмогорова для МГД-турбулентности, модель, основанная на перекрестной спиральности скорости и магнитного поля, для МГД-турбулентности, модель подобия масштабов для МГД-турбулентности и смешанная модель для МГД-турбулентности Все модели (кроме модели подобия масштабов) содержат константы, которые необходимо определить Для того чтобы решить проблему, связанную с выбором значений констант, была использована динамическая процедура, обобщенная в данной работе на случай сжимаемых МГД-уравнений, которая определяет значение модельной константы на каждом временном шаге

В разделе 5 4 представлены численные методы, которые использовались в работе при моделировании сжимаемой МГД-турбулентности Для интеграции по времени использовался модифицированный явный метод Рунге-Кутта третьего порядка точности, который требует меньше ресурсов оперативной памяти по сравнению со стандартным методом Рунге-Кутта Был разработан и создан численный код с конечно-разностными схемами четвертого порядка точности для системы уравнений магнитной гидродинамики, записанной в консервативной форме Однако для нелинейных слагаемых применяется кососимметричная форма, 'которая обеспечивает более точные результаты, так как уменьшаются ошибки, связанные с дискретизацией при использовании конечно-разностного подхода для моделировании турбулентных течений Для разделения

турбулентного течения на крупномасштабные и мелкомасштабные вихри применялся фильтр Гаусса четвертого порядка точности Так как в работе рассматривается трехмерное сжимаемое турбулентное МГД-течение, то использовалось последовательное применение одномерных фильтров

Для оценки эффективности различных предложенных подсеточных замыканий LES в разделе 5 5 использовались результаты, полученные прямым численным моделированием МГД-турбулентности в различных диапазонах магнитного числа Рейнольдса, гидродинамического числа Рейнольдса и числа Маха Всего было рассмотрено семь различных случаев с варьированием начальных параметров вычислений Приведены результаты расчетов сжимаемых МГД-течений методом LES с использованием пяти подсеточных параметризаций

Показано, что при уменьшении магнитного числа Рейнольдса разница между подсеточными моделями уменьшается для магнитной энергии и все предложенные подсеточные модели демонстрируют хорошее согласование с результатами DNS при малых значениях магнитного числа Рейнольдса При увеличении магнитного числа Рейнольдса растет роль подсеточных замыканий в моделировании сжимаемой МГД-турбулентности и уменьшается скорость диссипации магнитной энергии Наилучшие результаты показали модели Смагоринского, Колмогорова и модель, основанная на перекрестной спиральности для временной эволюции магнитной энергии Такое же поведение наблюдается и для перекрестной спиральности роль подсеточных параметризаций растет при увеличении магнитного числа Рейнольдса Для кинетической энергии при уменьшении магнитного числа Рейнольдса наблюдалось большее расхождение в результатах LES при применении различных подсеточных параметризаций Модель подобия масштабов продемонстрировала наихудшие результаты, остальные подсеточные замыкания увеличили точность расчетов Для временной динамики турбулентных напряжений, как магнитных, так и кинетических, характерно увеличение влияния на результаты моделирования МГД-турбулентности подсеточных параметризаций при увеличении магнитного числа Рейнольдса Роль анизотропии в расчетах и расхождение LES- и DNS-результатов для анизотропии турбулентности увеличивались при уменьшении магнитного числа Рейнольдса

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

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

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

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

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

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

В разделе 6 2 представлена полная система уравнений сжимаемой магнитной гидродинамики электро- и теплопроводящей жидкости, исходные уравнения МГД приводятся к безразмерному виду Получены отфильтрованные по Фавру уравнения для полной системы МГД-уравнений теплопроводящей жидкости, которая используется в методе

крупных вихрей для моделирования трехмерной сжимаемой турбулентности

' (10)

дрй, э рг _ 1 _ л2 1 — . Эг;

Э/ Эх, ; ' 11е„, Эх; Эх, ' и ;

+-[(Е + Р)и,--ТВВ и ]+-Г-^---а „и, ]-

Э/ Эх/' М2 ' Эх/РгЯеМ^г-!) Яе

Э т] - дВ ЬВ) _ Э 1 1 1 1 (13)

ът

Здесь £ - полная энергия, усредненная по Фавру, г =-А(—) -

Эх,

диссипация, связанная с теплопроводностью (закон Фурье) Для замыкания системы используется уравнение состояния в виде

-= тр Р /М '

Наличие уравнения для полной энергии в системе МГД-уравнений значительно усложняет решение задачи методом крупных вихрей Влияние мелкомасштабной турбулентности на отфильтрованную часть уравнений магнитной гидродинамики определяется через подсеточные слагаемые в правой части уравнений (11) - (13) Для определения тензора подсеточных напряжений г," и магнитного подсеточного тензора

напряжений т* в разделе 6 3 применялась модель Смагоринского для МГД-турбулентности, подробно рассмотренная в пятой главе Для определения подсеточного потока тепла QJ = р((и,Г)~ -й^Г) и

турбулентной подсеточной диффузии JJ=p((uJukukУ-йJ(uilllУ) в уравнении полной энергии используются параметризации в виде

(14)

. (15)

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

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

итоге для потока подсеточной магнитной энергии — {ВкВки]

и подсеточной энергии взаимодействия магнитного натяжения и

скорости - {В]Вкик -йкВкВк) была получена следующая

подсеточная модель

х\ (16)

J_ 2

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

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

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

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

В данной главе показано, что 1_Е5-метод с использованием явной средневзвешенной фильтрации демонстрирует хорошие результаты при моделировании теплопроводящей плазмы в сжимаемой МГД-турбулентности при различных числах Маха

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

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

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

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

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

Рис 2 Спектр кинетической энергии (слева) Нормализованный (умноженный на к5") сглаженный спектр кинетической энергии (справа) Видно, что степенной показатель спектра близок к к'1 для большей части турбулентного каскада Однако существует четко выраженный инерционный интервал колмогоровского типа уГ5"

Рис 3 Спектр плотности - сплошная линия, спектр флуктуаций плотности -пунктирная линия (слева) Нормализованный (умноженный на к5П сглаженный спектр флуктуаций плотности (справа) На рисунке (слева) оба графика имеют показатель спектра, близкий к Также существует четко выраженный инерционный интервал колмогоровского типа к~5'ъдля флуктуаций плотности (справа), что подтверждает наблюдательные данные

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

интервал, что соответствует уменьшению со временем гидродинамического числа Рейнольдса в затухающей турбулентности

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

Продемонстрировано, что МГД-турбулентность в условиях локальной межзвездной среды является анизотропной, что подтвердили данные наблюдений Получено различное поведение компонент скорости в спектральном каскаде при малых волновых числах для х-, у-и z - компонент поля скорости и почти отсутствие различий при больших фурье-модах, что означает возникновение турбулентных анизотропных каскадов, главным образом, на больших масштабах Также показано, что при больших значениях волнового числа, что соответствует мелким масштабам, спектр является изотропным

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

Публикации по теме диссертации

1 Belyan А V , Moiseev S S , Petrosyan A S Large-scale structures in turbulent multiphase flows// J Phys Condensed Matter 1990 SA P 469-475

2 Belyan A V , Moiseev S S , Petrosyan A S Sound transformation in a turbulent medium//Physics letters A 1991 V 155 P 181-184

3 Chkhetiani О G , Moiseev S S , Petrosyan A S , Sagdeev R Z The „ inverse energy cascade and self-organization in homogeneous turbulent

shear flow//Applied Scientific Research 1993 V 51 P 67-72

4 Chkhetiani О G , Moiseev S S , Petrosyan A S , Sagdeev R Z The large scale stability and self-organization in homogeneous turbulent shear flow// Physica Scripta 1994 V 49 P 214-220

5 Karelsky К V , Petrosyan A S Numerical simulations of near-surface phenomena on Mars//Adv Space Res 1995 V 16 P 645-648

6 Petrosyan A S , Polunine D Yu Reynolds stresses in rotating fluids// In Turbulence, Heat and Mass Transfer/ Ed К Hanjalic, T W J Peters Delft University Press 1997 ISBN 90-407-1465-7 P 549-556

7 Karelsky К V , Papkov V V , Petrosyan A S Particular solutions of shallow water equations on slopes// Physics Letters A 2000 V 271 P 341-348

8 Karelsky К V , Papkov V V , Petrosyan A S The initial discontinuity decay problem for shallow water equqtions on slopes// Physics Letters A 2000 V 271 P 349-357

9 Chernyshov A A , Karelsky К V , Petrosyan A S Subgrid-scale modeling in Large-Eddy Simulations of Compressible Magnetohydrodynamic Turbulence// Russian Journal of Numerical Analysis and Mathematical Modeling 2006 V 21 №1 P 1-20

10 Chernyshov A A , Karelsky К V , Petrosyan A S Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Compressible Fluid// Physics of Plasmas 2006 V 13 Issue 3 P 032304-032304-9

11 Karelsky К V , Petrosyan A S Particular solutions and Riemann problem for modified shallow water equations// Fluid Dynamics Research 2006 V 38 P 339-358

12 Карельский К В , Петросян А С Задача о стационарном обтекании ступеньки в приближении мелкой воды Известия РАН, Механика Жидкости и Газа, 2006, №1, С 15-24

13 Karelsky KV, Petrosyan AS, Slavin AG Quasi-Two-layer model for numerical analysis shallow-water flows on step// Russian Journal of Numerical Analysis and Mathematical Modeling 2006 V 21 №6 P 539-559

14 Chernyshov A A , Karelsky К V , Petrosyan A S Subgrid-scale modeling of compressible MHD turbulence in heat-conducting plasma// Physics of Plasmas 2006 V 13, Issue 10 P 104501-104501-4

15 Chernyshov A A , Karelsky К V , Petrosyan A S Large Eddy Simulation of Compressible Magnetohydrodynamic Turbulence in Heat-conducting Plasma// In Advances in Turbulence XI Proceedings of the EUROMECH European Turbulence Conference 2007 ISSN 0930-8989 ISBN 978-3540-72603-6 Springer, Berlin, Heidelberg, New York P 20-22

16 Chernyshov A A , Karelsky К V , Petrosyan A S Development of Large-Eddy Simulation for Modeling of Decaying Compressible MHD turbulence//Physics of Fluids 2007 V 19 Issue 5 P 055106-05510614

17 Karelsky К V , Petrosyan A S , Slavin A G Numerical Simulations of free-surface hydrodynamic flows on arbitrary bed// Russian journal of Numerical Analysis and Mathematical modeling 2007 V 22 №6 P 543-565

18 Karelsky К V , Petrosyan A S , Smirnov I V A new model for boundary layer flows interacting with particulates in land surface on complex terrain// Quarterly Journal of the Hungarian Meteorological Service 2007 V 11 №2-3 P 149-159

19 Chemyshov A A , Karelsky К V , Petrosyan A S Large Eddy Simulation of Compressible MHD Turbulence// Flow, Turbulence and Combustion 2008 V 20 №1 P 21-35

20 Chemyshov A A, Karelsky К V, Petrosyan A S Modeling of compressible magnetohydrodynamic turbulence in electrically and heat conducting fluid using large eddy simulation// Physics of Fluids 2008 V 20 Issue 8 P 085106-085106-13

21 Chemyshov A A , Karelsky К V , Petrosyan A S Three-dimensional modeling of compressible magnetohydrodynamic turbulence in the local interstellar medium//Astrophysical Journal 2008 V 686 P 1137

22 Karelsky К V and Petrosyan A S Modified Shallow water equations Simple waves and Riemann problem// Proceedings IUTAM Symposium on Hamoltonian Dynamics, Vortex Structures, Turbulence 2008 Springer ISBN 978-1-4020-6743-3 P 383-392

23 Petrosyan ASA new model for aerosols transportation in boundary layer on complex terrain// In Atmospheric boundary layer flows in air pollution modelling Von Karman Instituten for fluid dynamics/Ed J van Beeck, J-M Buchlin, A Petrosyan 2008 ISBN-13 978-2-930389-85-0

24 Petrosyan A S Numerical simulation of flows of a heavy aerosols transportation over a bed surface with an arbitrary profile// In Atmospheric boundary layer flows in air pollution modelling Von Karman Instituten for fluid dynamics/ Ed J van Beeck, J -M Buchlin, A Petrosyan 2008 ISBN-13 978-2-930389-85-0

25 Karelsky К V, Petrosyan A S, Slavin A G Finite-difference presentation of the Coriolis force in numerical Godunov-type models for flows of rotating shallow water// Russian journal of Numerical Analysis and Mathematical modeling 2009 V 24(3) P 229-259

26 Карельский К В , Петросян А С , Славин А Г Численный метод для исследования течений мелкой воды над сложным профилем дна в присутствии внешней силы// Математическое моделирование 2009 Т 21 №6 С 41-58

055(02)2

Ротапринт ИКИ РАН Москва, 117997, Профсоюзная, 84/32

Подписано к печати 14 07 2009

Заказ 2184 Формат 70x108/32 Тираж 100 2,0 уч -изд л

ВВЕДЕНИЕ.

ГЛАВА 1. ЧАСТНЫЕ РЕШЕНИЯ И ЗАДАЧА РАСПАДА ПРОИЗВОЛЬНОГО РАЗРЫВА ДЛЯ УРАВНЕНИЙ МЕЛКОЙ ВОДЫ НАД НЕОДНОРОДНОЙ ПОВЕРХНОСТЬЮ

1.1. Введение.

1.2. Инварианты Римана для уравнений мелкой воды над неоднородной поверхностью.

1.3. Частные решения УМВ на наклонной плоскости. Преобразование уравнений Сен-Венана к классическим уравнениям мелкой воды.

1.4. Задачи Римана для уравнений «мелкой воды» на наклонной плоскости.

1.5. Модернизированные уравнения мелкой воды. Инварианты Римана.

1.6. Непрерывные и разрывные решения модернизированных УМВ. Анализ полученных решений.!.

1.7. Задачи Римана для модернизированных уравнений мелкой воды.

1.8. Выводы.

ГЛАВА 2. ЗАДАЧА РИМАНА ДЛЯ ТЕЧЕНИЙ МЕЛКОЙ ВОДЫ НАД СТУПЕНЧАТЫМ ПРОФИЛЕМ ДНА

2.1. Введение.

2.2. Стационарные уравнения мелкой воды. Области допустимых решений.

2.3. Анализ возможных стационарных режимов обтекания ступеньки.

2.4. Формулировка квазидвухслойной модели течений на уступе дна.

2.5. Нестационарные волновые картины над уступом дна.

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

2.7. Обсуждение аналитического подхода и квазидвухслойной модели. Анализ численных решений.

2.8. Выводы.

ГЛАВА 3. ИСПОЛЬЗОВАНИЕ КВАЗИ-ДВУХСЛОЙНОГО ПОДХОДА ДЛЯ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ТЕЧЕНИЙ МЕЛКОЙ ВОДЫ НАД СЛОЖНОЙ ПОСТИЛАЮЩЕЙ ПОВЕРХНОСТЬЮ В ПРИСУТСТВИИ ВНЕШНИХ СИЛ

3.1. Введение.

3.2. Конечно-разностная схема для уравнений мелкой воды над подстилающей поверхностью произвольного профиля.

3.3. Численное моделирование течений жидкости над подстилающей поверхностью сложного профиля.

3.4. Представление силы Кориолиса в методах годуновского типа на основе квазидвухслойной модели.

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

3.6. Результаты численного моделирования вращающейся мелкой воды.

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

3.8. Исследование течений мелкой воды над произвольным профилем дна в присутствии внешней силы.

3.9. Выводы.

ГЛАВА 4. ТЕОРИЯ И МОДЕЛИ МНОГОФАЗНЫХ И СДВИГОВЫХ ТУРБУЛЕНТНЫХ ТЕЧЕНИЙ НЕЙТРАЛЬНОЙ ЖИДКОСТИ

4.1. Введение.

4.2. Крупномасштабные структуры в жидкости с твердыми частицами.

4.3. Гидродинамический а-эффект в жидкости с осциллирующими пузырьками газа.

4.4. Односкоростная трехмерная модель переноса атмосферных примесей течением ветра в областях со сложной границей.

4.5. Рассеяние звука спиральной турбулентностью.

4.6. Генерация крупномасштабных структур на сдвиговых течениях.

4.7. Выводы.

ГЛАВА 5. ПОДСЕТОЧНЫЕ МОДЕЛИ СЖИМАЕМОЙ МАГНИТОГИДРОДИНАМИЧЕСКОЙ ТУРБУЛЕНТНОСТИ В КОСМИЧЕСКОЙ ПЛАЗМЕ

5.1. Введение.

5.2. Уравнения магнитной гидродинамики политропного газа. Формулировка метода крупных вихрей.

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

5.4. Численные методы, используемые при моделировании сжимаемой МГД-турбулентности.

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

5.6. Выводы.

ГЛАВА 6. ИССЛЕДОВАНИЕ ТУРБУЛЕНТНОСТИ СЖИМАЕМОГО ТЕПЛОПРОВОДНОГО ГАЗА И КОСМИЧЕСКОЙ ПЛАЗМЫ МЕТОДОМ КРУПНЫХ ВИХРЕЙ

6.1. Введение.

6.2. Отфильтрованные уравнения магнитной гидродинамики теплопроводящей жидкости.

6.3. Параметризации подсеточных слагаемых для сжимаемой МГД-турбулентности теплопроводящего газа.

6.4. Численное моделирование и анализ результатов расчетов затухающей МГД-турбулентности теплопроводящего газа.

6.5. МГД-модель и численное исследование локальной межзвездной среды.

6.6. Анализ результатов моделирования и теоретическая интерпретация.

6.7. Выводы.

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

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

Сжимаемая магнитогидродинамическая (МГД) турбулентность является широко распространенным состоянием космической плазмы во многих астрофизических, гелиофизических, геофизических процессах. Например, в аккреционных дисках МГД турбулентность вызвана магниторотационной неустойчивостью. Большинство турбулентных явлений в физике Солнца описываются в рамках МГД: солнечный ветер, расширение солнечной короны, конвективная зона, фотосфера, солнечный тахоклин. Явления турбулентности наблюдаются в околоземном пространстве как в солнечном ветре, так и в в различных областях магнитосферы Земли, в частности, в дальней области геомагнитного хвоста наблюдаемые спутниками свойства космической плазмы адекватно можно объяснить только в рамках теории и моделей сжимаемой турбулентности. Сжимаемые турбулентные течения в магнитном поле также широко распространены и в прикладных областях. Среди инженерных применений можно указать возможность управления пограничным слоем и снижение сопротивления потоку.

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

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

Необходимость редукции исходных уравнений в классе задач со свободной поверхностью привела после предположения гидростатичности распределения давления, усреднению поля скорости по глубине потока и пренебрежения изменением горизонтальных скоростей вдоль линий коллинеарных вектору силы тяжести, к построению Стокером математической модели более низкого порядка. Впоследствии данная модель была названа моделью- «мелкой воды» [1,2], поскольку редукция осуществлялась асимптотическим разложением по малому параметру, определяемому отношением глубины жидкости к характерному линейному размеру. При наличии внешнего источника, например, силы Кориолиса, область применения модели не ограничена условием малости глубины жидкости по сравнению с характерными линейными размерами, поскольку в этом случае возможность редукции трехмерных уравнений к двумерным есть следствие вращения, а не тонкости слоя. Уравнения мелкой воды, являясь системой нелинейных гиперболических уравнений, аппроксимируют полную систему уравнений Эйлера, описывающую течения несжимаемой жидкости со свободной поверхностью в поле силы тяжести при пренебрежении эффектами вертикальной неоднородности горизонтального поля скорости.

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

Традиционно, океанические течения являются важным направлением применения уравнений мелкой воды. Например, в работе [4] приведены результаты моделирования крупномасштабных течений на всем земном шаре, учитывающие влияние сил, вызываемых воздействием Солнца'и Луны. Интересный,пример применения приближения мелкой воды приведен в работе [5], в которой решены усредненные уравнения > мелкой воды для остаточных циркуляций. Уравнения работы [5] усреднялись по периоду течения, нелинейные члены при таком усреднении описывали суммарный- поток, импульса (радиационное напряжение), который определяется волновой компонентой течения. Важным приложением; уравнений мелкой воды является предсказание штормовой картины, а именно, генерация течений1 и вариаций уровня воды, вызванных отношением атмосферного давления, и напряжением ветра на водной поверхности [6,7].

Если жидкость расслаивается по причине разной солености, то полученный в результате слоистый поток очень похож на.течение мелкой воды. Пример использования многослойной модели мелкой воды приведен в работе [8]. Схожая многослойная модель использовалась в работе [9] для моделирования и описания Большого Красного Пятна в атмосфере Юпитера.

Также следует отметить актуальные применения приближений мелкой воды для описания гидравлических течений [10,11], береговых течений [12], течений в реках [13] и озерах [14], течений в водозаборниках, технических сужениях и лотках [1], моделирования, цунами [15,16], распространения волн прорыва и приливных бор в реках [17], распространения тяжелых газов и примесей в атмосферах планет [9], атмосферных движений крупных масштабов, используемых при предсказании погоды [3,18,19].

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

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

Простые автомодельные решения гиперболических систем уравнений [16,22,23,24,25,26,27,28] являются основополагающими в исследовании нелинейных волновых явлений, поскольку позволяют найти точное решение задачи распада произвольного разрыва. Задача Коши о распаде произвольного разрыва кусочно-постоянных начальных условий, впервые возникшая в газовой динамике (задача Римана) [29,30], имеет фундаментальное, значение. Ее решение облегчает понимание множества нелинейных явлений в течениях несжимаемой жидкости со свободной границей в рамках приближения мелкой воды. Решение такой задачи на ступенчатой границе особенно важно при изучении течений над сложной подстилающей поверхностью, аппроксимируемой системой таких ступенчатых границ. Существующий аналитический подход для решения задачи Римана, основанный на предположении о наличии стационарной зоны в окрестности ступенчатой границы, с одной стороны не учитывает потери кинетической энергии на турбулентное перемешивание вблизи ступеньки, с другой - накладывает ограничения на возможные глубины натекающего потока для преодоления ступенчатой границы. Решение этой задачи является важнейшим компонентом в разработке численных методов сквозного счета, основанных, как на аналитических, так и на приближенных решениях задачи Римана. Основная идея методов [21,31,32,33,34,35,36] годуновского типа [37], основанных на решении задачи Римана, состоит в расщеплении решения многомерной задачи на набор одномерных подзадач, возникающих после разбиения расчетной области на ячейки и записи соответствующих интегральных соотношений для всех элементов (ячеек), посредством которых осуществлялось разбиение. Данные методы особенно часто находят применение в численном моделировании, поскольку позволяют получать решение, не только в области непрерывного течения, но и в областях разрыва решений, без специального выделения и отслеживания поверхностей разрыва [10,38,39,40,41,42,43,44,45,46]. Кроме того, эти методы хорошо адаптируются к сложным граничным условиям, характерным для большинства постановок задач, описывающих реальные природные течения [45].

При численном интегрировании гиперболической системы балансовых уравнений ключевой задачей является получение конечно-разностных схем, удовлетворяющих сохранению устойчивых состояний, таких как, например, равновесие покоящейся воды. Схемы, удовлетворяющие таким свойствам, называют well-balanced схемами. В'работе [47] для выполнения условий сохранения устойчивых состояний было предложено использовать модифицированную схему Роу [48]. Дальнейшее развитие идеи этой работы получило в статьях [49,50], в работах [51,52] рассмотрены явления турбулентности и массовые источники в задачах, включающих неоднородную топографию, в работах [53,54] применена улучшающая реконструкция. В-работах [34,35] было предложено использовать либо точные, либо модернизированные схемы годуновского типа (Riemann-solvers). Такой подход хорошо себя зарекомендовал, и данная идея получила широкое развитие для- решения задач с источниковыми членами различного происхождения. В частности, в работах [21,55] предложен метод, компенсирующий влияние подстилающей поверхности, изменением значений глубин жидкости. В работах [31,56,57] предложен метод, использующий схему гидростатической реконструкции подсчета потоков на гранях ячеек. С помощью этой схемы промоделированы явления, обусловленные наличием силы Кориолиса, в частности явление геострафической адаптации. Изучению этих же явлений посвящена работа [33], в которой применен метод Левека [34,35]. Распространение волн цунами и процессы затопления изучены в работах [58,59]. В работе [60] был предложен метод поверхностных градиентов -способ интерпретации источниковых членов в уравнениях мелкой воды, основанный на точной реконструкции консервативных переменных на гранях ячеек. Метод разделения потоковых величин, который, когда это необходимо, поддерживает точный баланс между градиентами потока и источниковыми членами, предложен в работе [61]. В работе [62] предложена нелинейная реконструкция параметров, обеспечивающая высокий порядок точности разностной схемы для уравнений мелкой воды. Точные частные решения задачи, включающую ступенчатую границу на дне, представлены в работах [63]. Точные частные решения для задачи, включающую неоднородную топографию, представлены в работе [64].

Также следует отметить работы [10, 32,36,41,56,65,66,67,68,69,70], внесшие существенный вклад в решение задач моделирования потоков жидкости над неоднородной поверхностью с источниковыми членами различного происхождения.

Для решения многих реальных задач недостаточно учитывать только неоднородность подстилающей поверхности. Необходимость рассмотрения дополнительных воздействий, определяемых конкретными условиями течения, приводит к появлению в системе уравнений мелкой воды дополнительных членов. В частности, для расчета крупномасштабных атмосферных и океанических задач следует принимать во внимание эффекты, определяемые планетарным вращением. Наличие хорошо разработанного и апробированного численного аппарата вкупе с многократно протестированной программной реализацией сделало особенно привлекательным сведение решения задачи о вращающейся «мелкой воде» над ровной подстилающей поверхностью к решению- задачи о течениях мелкой воды над комплексной нестационарной границей. Представление силы Кориолиса фиктивной нестационарной границей* создает важные преимущества при моделировании течений на неровной границе, сводя задачу к моделированию течений мелкой воды над нестационарной эффективной поверхностью [31,33,34,35,57]. Однако, применение такого представления в расщепляющихся численных методах затрудняется отсутствием одномерной постановки-задач для вращающейся жидкости. Формальное постановка одномерной задачи, определяемой' отказом от частных производных по одному из пространственных направлений, делает особенно актуальным нахождение горизонтальной неоднородности трансверсальной; составляющей вектора скорости, определяющей консервативность силы Кориолиса, в зависимости от вертикальной структуры течения;

В последнее время, большое внимание уделяется проблемам возникновения, структур в неравновесных средах. Особый интерес вызывают крупномасштабные вихри, возникающие в турбулентной Жидкости, - когерентные структуры. Такие структуры наблюдаются экспериментально [71, 72]. Когерентные структуры наблюдаются также и в природных овиях: вихри Россби в атмосфере Земли и Юпитера, тропические циклоны, [73-75]. При изучении явления, турбулентности важную роль играет модель однородной изотропной и стационарной турбулентности. Возникает вопрос, может ли такая турбулентность усиливать крупномасштабные возмущения: В магнитной гидродинамике развита теория усиления, и поддержания магнитных полей мелкомасштабной турбулентностью [76-80]. Краузе и Рэдлер [78] показали, что исходно однородная, изотропная и зеркально-симметричная турбулентность не может усиливать крупномасштабные магнитные поля. Однако, если зеркальная симметрия турбулентности нарушена, то такая среда может усиливать крупномасштабное магнитное поле. Это явление получило название магнитного динамо. Простейший пример нарушения зеркальной симметрии представляет собой поле скоростей, в котором величина средней спиральности отлична от куля. Такое турбулентное поле скоростей характерно тем, что правовинтовые и левовинтовые вихри наблюдаются с разной вероятностью, то есть вихрей одного знака больше, чем другого.

Хорошо известно, что уравнение для завихренности подобно уравнению индукции магнитного поля [81,82]. Поэтому вполне естественно желание изучать эволюцию крупномасштабных вихревых возмущений, тем более, что существование инварианта спиральности в невязкой гидродинамике обеспечивает ограничение потока энергии от крупных масштабов к меньшим [83,84]. Действительно, дробление масштабов определяется нелинейным слагаемым уравнения гидродинамики. В работе [85] показано, что увеличение спиральности уменьшает слагаемое, определяющее взаимодействие мод в гидродинамических уравнениях. Важную роль спиральности в возникновении крупномасштабных вихревых структур в атмосфере отмечают авторы [86]. Однако, несмотря на отмеченную выше аналогию между уравнением для-завихренности и индукции магнитного поля в работе [87] Краузе и Рюдигер в корреляционном приближении второго порядка исследовали напряжения Рейнольдса для- случая однородной, изотропной и спиральной турбулентности и показали, что в усредненных уравнениях мелкомасштабная турбулентность модифицирует вязкость. Этот результат является следствием симметрии напряжений Рейнольдса в несжимаемой жидкости. Тем* не менее обратный энергетический каскад в спиральной, турбулентности возможен. Для» этого требуется дополнительное нарушение симметрии напряжений Рейнольдса.-В работе [88]?(см. также [89, 90]) впервые получена генерация, крупномасштабных вихрей в сжимаемой жидкости и получено осредненное уравнение, описывающее эволюцию завихренности. В- [92] подробно исследуется роль сжимаемости среды на генерационные свойства турбулентности, а в [92] изучается усиление вихревых возмущений в сжимаемой самогравитирующей среде. Для описания крупномасштабных процессов в атмосфере представляет интерес эволюция вихрей в несжимаемой жидкости. В работе [93] исследована, возможность генерации вихрей в спиральной турбулентности с неоднородным потоком и показано, что добавления потенциального потока уже достаточно для появления нетривиальных слагаемых в усредненных уравнениях. Для бессдвигового и линейного по координатам потока найдено решение, описывающее растущие вихревые возмущения поля скорости. В* работах [94-96] исследовано вихревое динамо конвективной среде со спиральной турбулентностью и показано, что учет спиральной турбулентности приводит к существенному изменению характера конвекции и к положительной обратной связи тороидальных и полоидальных полей. В результате, при значениях спиральности выше критических возникает неустойчивость. Теория, изложенная в [94-96] имеет прямое отношение к зарождению тайфунов в атмосфере Земли. Для применения этой теории к описанию возникновения тропических циклонов в работе [97] изучается процесс зарождения крупномасштабной неустойчивости в конвективной турбулентности с учетом фазовых переходов влаги. Показано, что процесс развития при влажной конвекции существенно отличается от случая сухой конвекции. Учет фазовых переходов влаги уменьшает порог неустойчивости, то есть это эквивалентно появлению дополнительного источника эиергии.

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

Система уравнений двухскоростного движения бесстолкновительной дисперсной смеси негиперболична вследствие недостаточно полного описания межфазного взаимодействия ц взаимодействия между дисперсными частицами. Но, несмотря на это, существует достаточно широкий класс задач, когда эта система уравнений правильно отражает физику процесса и физическую неустойчивость некоторых течений. Для таких задач постановка задачи Коши «условно корректна» в определенном классе функций, т. е. в этом классе функций решение задачи Коши существует, единственно и непрерывно зависит от начальных условий. Зависимость скорости развития неустойчивости от длины волны может накладывать ограничение снизу на выбор шага расчетной сетки при численном решении задачи на основе уравнений двухскоростного движения. Структура стационарных ударных волн с передним скачком в газе с твердыми частицами, когда отсутствует фазовый переход, была исследована в работе [98]. В работе [99] разработан алгоритм сквозного счета дифференциальных уравнений одномерного нестационарного движения двухскоростной среды в эйлеровых переменных с использованием разностных схем метода «крупных частиц» [100] и метода [101]. В литературе для расчета одномерных нестационарных течений газовзвесей рассмотрены различные методы [102]. В работе [103] предложен конечно-разностный метод интегрирования уравнений в лагранжевых переменных г, s газовой фазы, где г —• координата материальных точек газовой фазы в начальный момент времени. Естественно, что дисперсные частицы, имея отличную от газа скорость, перемещаются относительно лагранжевых координат газа. Некоторые результаты расчетов представлены также в статье [104].

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

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

Наиболее подробную информацию о турбулентном течении жидкости можно получить с помощью прямого численного моделирования - DNS (Direct numerical simulation) [89], которое заключается в численном решении полной нестационарной системы магнитогидродинамических уравнений. При таком подходе разрешаются все масштабы движения заряженной жидкости. Метод DNS не требует специальных замыканий для уравнений магнитной гидродинамики. Прямой численный расчет МГД турбулентности сталкивается с принципиальными трудностями, связанными с большими гидродинамическими и магнитными числами Рейнольдса, которые характерны для исследуемых процессов, так как в этом случае число степеней свободы турбулентного движения велико и минимальное количество узлов на численной сетке должно быть столь большим, что ограничивает применение прямого численного моделирования для изучения турбулентных течений с реальными характерными числами Рейнольдса.

Осборн Рейнольде предложил статистический подход для исследования турбулентных течений [105], который заключается в осреднении уравнений движения - RANS (Reynolds averaged Navier-Stokes). В методе RANS все параметры движения разлагаются на среднюю и турбулентную составляющие. В уравнении Навье-Стокса появляются рейнольдсовские напряжения, которые необходимо замкнуть. Следовательно, вся турбулентность моделируется (например, к-е модель [106,107]), а не высчитывается, как в методе DNS. Метод RANS обычно используется для теоретических исследований средних течений. Этот подход не содержит информации о динамике турбулентности.

Метод крупных вихрей LES (Large eddy simulation) - это метод, который описывает приближенную динамику турбулентности, где крупномасштабная часть турбулентного потока высчитывается непосредственно, а мелкомасштабная - моделируется, то есть LES является промежуточным подходом к изучению турбулентности между DNS и RANS. Это видно на рисунке 1, который показывает различие между тремя вычислительными методами, применяемыми для исследования турбулентности.

В методе LES используется операция фильтрации для разложения характеристик турбулентного движения на крупномасштабную и мелкомасштабную части, что связано с достаточной изотропностью, однородностью и универсальностью мелких масштабов турбулентного движения. Мелкомасштабное движение исключается из исходной системы уравнений движения с применением процедуры фильтрации и дальше их влияние на движение моделируется с использованием подсеточных моделей SGS (subgrid scale) (или, другое название, SFS (subfilter scale)), выраженных через отфильтрованные параметры турбулентных течений. Крупномасштабное движение рассчитывается из решения отфильтрованных нестационарных уравнений магнитной гидродинамики. LES является методом для моделирования течений с большими числами Рейнольдса, так как в методе крупных вихрей предполагается, что энергия переносится от больших масштабов к малым только внутри инерционного интервала, поэтому число степеней свободы будет меньше, чем в методе DNS, следовательно, LES требует значительно меньших вычислительных затрат по сравнению с DNS.

Изначально метод крупных вихрей развивался для моделирования гидродинамической турбулентности нейтральной жидкости [108, 109; 110, 111, 112, 113] для изучении задач метеорологии и океанологии. Большая часть работ была выполнена для несжимаемых течений. Использования LES к сжимаемым средам встречается значительно реже, вследствие увеличения сложности задачи из-за необходимости решения уравнения энергии. В отфильтрованном уравнении энергии появляются сразу несколько дополнительных подсеточных (SGS) слагаемых, которые необходимо параметризовать. Впервые LES рассматривался и использовался к сжимаемой жидкости в работе [114]. Первые применения метода LES к сжимаемым течениям рассматривали транспортное уравнения для внутренней энергии на единицу массы [115, 116], для давления [117] или для удельной энтальпии [118, 119]. В работах [120, 121] было предложено использовать уравнение полной энергии для замыкания системы гидродинамических уравнений нейтральной жидкости, причем некоторые подсеточные слагаемые имеют точно такой же вид, как и в уравнениях внутренней энергии или энтальпии. Подробная информация о различных подсеточных моделях метода крупных вихрей для случая сжимаемой жидкости содержится в работе [122]. В этой статье авторы рассматривают и тестируют параметризации для различных видов уравнений энергии: внутренней энергии, энтальпии и полной энергии.

Метод LES к течениям электропроводящей жидкости применялся и исследовался крайне мало. Все предыдущие работы в этом направлении ограничивались только рассмотрением несжимаемой жидкости для решения индустриальных задач. В работах [123, 124, 125, 126, 127, 128] авторы использовали LES для изучения несжимаемой МГД турбулентности, в статье [129] рассматривалось влияние магнитного поля в методе LES на течение несжимаемой проводящей жидкости при низких значениях магнитного числа Рейнольдса без использования уравнения для магнитной индукции. Во всех вышеупомянутых статьях несжимаемая система МГД уравнений рассматривалась без использования уравнения энергии. Для замыкания системы уравнений магнитной гидродинамики предполагалась политропность (или адиабатичность) процесса, или давление рассматривалось лишь как пассивная величина, которая обеспечивала несжимаемость МГД турбулентности. Тем не менее, многие течения электропроводящей жидкости не могут быть описаны в рамках несжимаемой среды, или сжимаемыми уравнениями в приближении политропии, а необходимо рассматривать теплопроводящую жидкость с использованием уравнения энергии. Применение LES для сжимаемой теплопроводящей МГД жидкости значительно усложняется из-за того, что нужно решать уравнение энергии, в котором появляются дополнительные слагаемые из-за наличия магнитного поля. К тому же, после фильтрации появляются добавочные подсеточные члены, требующие разработки теории для их параметризации.

Вообще, полные нелинейные трехмерные уравнения магнитной гидродинамики (включая диссипативные, тепловые, диффузионные и сжимаемые эффекты) столь сложны, что поддаются лишь приближенному численному решению. Однако, из-за того, что для космических МГД течений характерны большие числа Рейнольдса и числа Маха отличны от нуля, моделирование сжимаемой МГД турбулентности ограничено вычислительными ресурсами и встречается намного реже, чем для несжимаемых сред. Поэтому зачастую пользуются упрощенными моделями, пренебрегая некоторыми эффектами. Например, моделирования идеальных МГД течений, когда пренебрегают диссипацией и теплопроводностью и считают, что проводимость плазмы бесконечная [130]. В этом случае система МГД уравнений становится гиперболической, а не параболической, как для диссипативной системы уравнений, что упрощает численные решения, так как можно использовать хорошо развитые годуновские схемы различного порядка точности (в результате решение трехмерной задачи сводится к решению серии одномерных задач, численные потоки в каждом пространственном направлении вычисляются на основе соответствующей одномерной задачи Римана о распаде произвольного разрыва [55]). Существуют несколько работ, где используют TVD схемы, ENO/WENO схемы, схем на основе принципа минимального значения производной и т.д. для решения уравнений вязкого сжимаемого газа для МГД случая [131] путем добавления численных вязких потоков к соответствующим невязким потоках!, однако это часто нарушает монотонность разностной схемы и может привести к некорректным результатом, обычно подобного рода реконструкции используются для невязкой жидкости. Иногда используют квазиупругое приближение (anelastic approximation) [132, 133, 134] для МГД моделирования, когда предполагается, что звуковые моды отсутствуют или стационарны, однако такого рода приближение используются в основном только при моделировании конвективных зон Солнца и звезд. Часто межзвездную и межпланетную среду, а также солнечный и звездный ветер, рассматривают предполагая политропное (или адиабатическое) соотношение между плотностью и давлением для замыкания системы уравнений, в этом случае соображения о температуре процесса не являются основными и система сжимаемых МГД уравнений решается без уравнения энергии [135, 136, 137, 138, 139, 140]. Еще одним упрощением является рассмотрение сжимаемой двухмерной МГД турбулентности [141, 142, 143, 144], причем в работах [141, 142] в качестве начальных условий для скорости и магнитного поля использовалось детерминированное, неслучайное распределение (так называемый вихрь Орсзага-Танга [145, 146]), однако случайное распределение начальных значений скорости и магнитного поля является более подходящим условием для космических применений в МГД моделировании. Подробное влияние магнитного числа Рейнольдса на двумерное магнитогидродинамическое течение при различных начальных условиях описывается в работе [144]. Двухмерная МГД турбулентность существенно отличается от трехмерной турбулентности, так как в двухмерном потоке (если пренебречь вязкостью) сохраняется средняя завихренность, в то время как в трехмерном вихревые трубки деформируются и завихренность не является инвариантом движения, также турбулентные динамо-процессы и генерация крупномасштабного магнитного поля возможны только в трехмерном случае МГД турбулентности [147]. Иногда при вычислении космических течений предполагают, что плазменная бета (отношение давления плазмы к энергии магнитного поля) столь велико, что пренебрегают магнитным полем и задача сводится к гидродинамической и решается система уравнений для движения нейтральной жидкости [148, 149, 150].

Исследование сжимаемой турбулентности как в гидродинамике нейтральной жидкости, так и в магнитогидродинамике является трудной задачей, поскольку не существует аналитической или приближенной теории таких явлений. Однако, несмотря на существенную роль сжимаемости в космической плазме, целый ряд наблюдений показывает воспроизведение колмогоровского спектра флуктуаций плотности [151, 152, 153, 154] (этот факт продемонстрирован на рисунке 2). На рисунке 2 приведен трехмерный спектр плотности электронов в межзвездной среде, полученный с помощью различных прямых и косвенных наблюдений [151], в работе [151] поясняются различные наблюдательные методы, изображенные на этом рисунке. Для интерпретации таких наблюдений была предложена теория "почти несжимаемой" (nearly incompressible) среды, которая описывает флуктуации плотности в гидродинамике нейтрального [155] и магнитного [156] газа в режиме переноса пассивного скаляра. В работе [150] аналитическая теория для сжимаемого нейтрального газа была подтверждена прямым численным моделированием только для двухмерного случая. В работе [140] прямым численным моделированием было показано, что в сжимаемой магнитогидродинамике существует аналогичный эффект уменьшения локального турбулентного числа Маха со сверхзвукового режима в дозвуковой, что соответствует режиму преобразования сверхзвуковых турбулентных флуктуаций в дозвуковые. Однако в этой работе в силу ограничений метода прямого численного моделирования не удалось получить спектры плотности и кинетической энергий и показать их совпадения и реализацию пассивного режима для плотности в сжимаемой МГД турбулентности. Несмотря на это авторы этой работы [140] используют эти результаты для интерпретации спутниковых данных о солнечном ветре и локальной межзвездной среде.

Цель работы

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

В соответствии с поставленной целью в работе решаются следующие задачи:

1. Изучить точные нелинейные решения уравнений мелкой воды над неровной поверхностью pi над наклонным дном. Обобщить классическую теорию мелкой воды на случай неоднородных потоков.

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

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

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

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

6. Осуществить численное моделирование методом крупных вихрей затухающей сжимаемой МГД турбулентности в политропном газе и теплопроводящей плазме. Сравнить результаты моделирования метода крупных вихрей с результатами, полученные методом прямого численного моделирования. Исследовать спектры флуктуаций плотности и энергии в сжимаемой магнитогидродинамической турбулентной плазме для политропного уравнения состояния. Показать возможность существования режима слабосжимаемых турбулентных пульсаций, когда флуктуации плотности являются пассивным скаляром. Исследовать свойства анизотропии сжимаемой турбулентности и намагниченности плазмы в слабосжимаемом режиме.

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

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

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

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

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

Практическая и научная ценность работы.

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

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

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

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

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

Результаты могут найти применение в исследованиях атмосферы Земли и других планет, ведущихся в ИКИ РАН, ИПМ РАН, ФИ РАН, ИО РАН, ИФЗ РАН, ИФА РАН, ИВМ РАН, ИНАС РАН, ИЗМИР РАН, ИАЭ им. Курчатова, ИГ СО РАН, НИРФИ, МИФИ, МФТИ, МГУ, ИВТАН.

Обоснованность и достоверность полученных результатов.

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

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

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

Положения выносимые на защиту

1. Показано, что обобщение частных автомодельных решений уравнений мелкой воды над ровным- дном на случай неоднородной границы имеет место только для тех типов поверхностей, для которых существует решение типа «простой волны», а именно, для наклонной плоскости. Найдены все частные автомодельные решения уравнений Сен-Венана и показано, что характеристики уравнений мелкой воды на наклонной плоскости являются семейством парабол, имеющих точку касания 2-го порядка с соответствующими характеристиками классических уравнений мелкой воды над ровным дном. Решена задача о распаде произвольного разрыва для уравнений Сен-Венана на наклонной плоскости. Найденное решение описывает новые физические явления, определяемые дополнительной скатывающей силой. Для течений мелкой воды на ровной поверхности при наличии слабых вертикальных неоднородностей начальных условий предложены модифицированные уравнения "мелкой воды", найден новый безразмерный параметр, ограничивающий пределы применимости классических уравнений мелкой воды. Найдено решение задачи Римана для модифицированных уравнений "мелкой воды". Показано, что учет вертикальных неоднородностей исключает зону вакуума из конфигураций, реализуемых для классических уравнений "мелкой воды".

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

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

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

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

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

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

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

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

Основные положения диссертационной работы и полученные результаты докладывались на российских и международных конференциях и симпозиумах: Всесоюзной научной конференции по нелинейным явлениям, Звенигород 1989; Всесоюзной школе-семинаре по нелинейной математической физике, Светлогорск, 1989; XLII научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук» -1999; Международной конференции МСС-04 «Трансформация волн, когерентные структуры и турбулентность». Москва, 2004; XLVII научной конференции МФТИ. Современные проблемы фундаментальных и прикладных наук. Часть VIII. Физика и энергетика. Москва, 2004; XLVIII научной конференции МФТИ. Современные проблемы фундаментальных и прикладных наук. Часть VIII. Физика и энергетика. Москва, 2005; XVI Научная сессия Совета РАН по нелинейной динамике. Москва, 2007; XVII Научная сессия Совета РАН по нелинейной динамике. Москва, 2008; Международная конференция «Современные проблемы газовой и волновой динамики», Москва 2009; Международная конференция «Потоки и структуры в жидкостях: физика геосфер», Москва 2009; Kiev International Workshop on Non-linear and Turbulent Processes in Physics. 1989; International Symposium "Generation of large-scale structures in continuous media", Perm-Moscow, 1990; 1-st Liquid Matter Conference, Lyon , 1990; 4-st European Turbulence Conference, Delft , 1992; IUTAM Symposium " Eddy Structures Identification in Free Turbulent Shear flows, Poitiers, 1992; COSPAR assembly, Hamburg, 1994; EGS Assembly, Hamburg, 1995; ESA workshop on MARS, Capri, 1995; ERCOFTAC workshop on "Compressible turbulence". Chatillon sous Bagneux, France, 1995; EGS Assembly, the Hague, 1996; EUROMECH colloquium, Flow Control, Berlin 1997; EGS Assembly, Vienna, 1997; 2nd Int. Symposium on Turbulence, Heat and Mass transfer, Delft, 1997; EUROMEH colloquium on Statistics and Dynamics of Vorticity, Marseille, 1997; European Fluid Mechanics Conference, Gettingen, 1997; EGS Assembly, Nice, 1998; International Conference on Air Pollution Modelling, Champs-sur-Marne, 1998; International Conference on

Control in Fluid Dynamics, Cachan, France, 1998; EGS Assembly, Den Haag, 1999; International Conference " Systems Sciences 2000 " Osnabruek, 2000; EGS Assembly, Nice,2000; 8-th European Turbulence Conference, Barcelona, 2000; EGS Assembly, Nice, 2001; Second International Conference on Air Pollution Modelling, Champs-sur-Mame, 2001; Climate Conference 2001, Utrecht, 2001; EUROMRCH colloquium, Strongly Coupled Multiphase Flows, Grenoble, 2001; EGS Assembly, Nice,2002; EGS-AGU-EGU Assembly, Nice, 2003; EGU Assembly, Nice, 2004; EGU Assembly, Vienna, 2005; EMS meeting, Utrecht, 2005; EGU Assembly, Vienna, 2006; NATO Advanced Research Workshop on "Atmospheric Boundary Layers: Modelling and Applications for Environmental Security", Dubravnik, 2006; Workshop on Environmental fluid mechanics-as elements in agro meteorological modeling, As, Norway, 2006; International Conference " Turbulence and Interactions", Porqueroles, France, 2006; International Conference "Mathematical Hydrodynamics", Moscow, 2006; IUTAM Symposium " Hamiltonian Dynamics, Vortex Structures, Turbulrnce" Moscow 2006; EGU Assembly, Vienna 2007; 8th International.School/Simposium for Space Simulations (ISSS-8), Kauai, USA, 2007; EMS meeting, SanLorenzo, Spain 2007; International Sympoposium on " Environmental Hydraulics, Tempe, Arizona, 2007; EGU Assembly, Vienna 2008; EUROMECH Colloquium on " Mixing of Coastal , Estuarine andRiverine Shallow Flows", Ancona, Italy, 2008; Lecture Series 2008-06 Atmospheric boundary layer flows in air pollution modelling. Von Karman Instituten for fluid dynamics, Rhode Saint Genese - Belgium; ERCOFTAC workshop "Direct and*Large Eddy Simulations 7", Trieste, Italy, 2008; International conference "Meso-scale Meteorology and Air Pollution" Odessa, Ukraine, 2008; EUROMEH Fluid Mechanics Conference, Manchester, UK, 2008; EMS meeting, Amsterdam 2008!

Публикации по теме диссертации.

Основные результаты работы опубликованы в 21 статьях в реферируемых российских и международных журналах, рекомендованных ВАК, в 5 статьях опубликованных в рецензируемых трудах российских и международных конференций. Результаты работы представлены в 45 тезисах докладов российских и международных конференций.

Личный вклад автора

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

Структура и объем диссертации.

Диссертационная работа состоит из введения, шести глав, заключения и списка цитированной литературы. Объем диссертации - 371 страница. Библиография включает 297 наименований. Диссертация содержит 157 рисунков, 2 таблицы.

6.7. Выводы

В данной главе разработан и изучен LES метод для магнитогидродинамической сжимаемой турбулентности электро- и теплопроводящей жидкости. Получена система отфильтрованных уравнений магнитной, гидродинамики при наличии уравнения полной энергии. В отличии от политропного случая (детально рассмотренного в главе ) вследствие присутствия уравнения энергии появляются новые подсеточные слагаемые в уравнении полной энергии. При рассмотрении полной системы уравнений сжимаемой МГД без предположении о политропности или адиабатичности.процесса, увеличивается количество подсеточных слагаемых, которые необходимо параметризовать. По сравнению с гидродинамикой нейтральной жидкости, в магнитогидродинамике появляются сразу несколько новых подсеточных членов. Предложены новые подсеточные модели для SGS слагаемых, появляющиеся после операции фильтрации, в уравнении полной энергии при наличии магнитного поля. Эти модели были применены для LES-расчетов со следующими числами Маха: Ms= 0.38, Ms= 0.65, Ms = 1.45 и результаты сравнивались с DNS результатами, а также с LES, но без использования подсеточных параметризаций в уравнении полной энергии, а только в уравнениях сохранения магнитной индукции и количества движения. На кинетическую и магнитную энергию учет подсеточных слагаемых в уравнении полной энергии почти не оказывает никакого эффекта, даже при высоких числах Маха, в то же время для температуры (соответственно и для внутренней энергии) наличие SGS моделей в уравнении полной энергии является важным условием для повышения точности вычислений термодинамических величин. При увеличении значения числа Маха увеличиваются осцилляции кинетической энергии и температуры.

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

В главе исследовалась сжимаемая МГД турбулентность в локальной межзвездной среде методом LES для решения системы магнитогидродинамических уравнений. Несмотря на то что для локальной межзвездной среды характерны сверхзвуковые течения с высокими крупномасштабными числами Маха, также существуют дозвуковые флуктуации слабосжимаемых компонент межзвездной среды. Именно эти слабосжимаемые дозвуковые флуктуации отвечают за появления спектра колмогоровского типа в локальной межзвездной турбулентности, который наблюдается из экспериментальных данных. В этой главе показано, что флуктуации плотности являются пассивным скаляром в поле скорости в умеренно сжимаемой магнитогидродинамической турбулентности и демонстрируют колмогоровский спектр в диссипативном интервале энергетического каскада. Показатель степенного спектра для плотности и кинетической энергии почти совпадают и близок к к~ъ для сжимаемой затухающей МГД турбулентности. Также продемонстрировано изменение со временем спектра кинетической энергии, показывающей уменьшение энергосодержащих крупных вихрей и инерционного интервала, и увеличение диссипативного масштаба. Показано, что турбулентное число Маха уменьшается значительно со сверхзвукового режима турбулентности (м5> 1), где рассматриваемая среда сильно сжимаемая до дозвукового значения (м5<1), характеризующее слабосжимаемое течение. Этот вывод, об уменьшении роли сжимаемости в турбулентных флуктуациях, подтверждается временной эволюцией дивергенции скорости, которая уменьшается и стремится к нулю (но не нуль). В локальной межзвездной среде переход плазмы от существенно сжимаемого МГД турбулентного течения к умеренно сжимаемому течению не только преобразовывает сверхзвуковое движение в дозвуковое, но также приводит к ослаблению намагниченности плазмы, что было показано в данной главе, так как плазменная бета р увеличивается со временем, таким образом роль магнитной энергии падает по сравнению с давлением плазмы.

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

ЗАКЛЮЧЕНИЕ

1. Показано, что обобщение частных автомодельных решений уравнений мелкой воды над ровным дном на случай неоднородной границы имеет место только для тех типов поверхностей, для которых существует решение типа «простой волны», а именно, для наклонной плоскости. Найдены все частные автомодельные решения уравнений Сен-Венана и показано, что характеристики уравнений мелкой воды на наклонной плоскости являются семейством парабол, имеющих точку касания 2-го порядка с соответствующими характеристиками классических уравнений мелкой воды над ровным дном. Как следствие, в отличие от классического случая, сильные разрывы распространяются по параболической траектории с постоянным ускорением. Показано, что конечно-разностные соотношения на поверхности разрыва не зависят от угла наклона плоскости. Найдено невырожденное преобразование зависимых и независимых переменных, сводящее уравнения Сен-Венана для наклонной плоскости к классическим уравнениям мелкой воды над ровным дном. Решена задача о распаде произвольного разрыва для уравнений Сен-Венана на наклонной плоскости. Найденное решение описывает новые физические явления, определяемые дополнительной скатывающей силой.

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

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

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

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

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

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

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

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

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

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

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

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

1. Stoker J.J. Water Waves: The Mathematical Theory with Applications. New York: Interscience, 1957.

2. Овсянников Л. В. К обоснованию теории мелкой воды, Сборник «Динамика сплошной среды» вып. 15, СО АН СССР, 1973.

3. Charney J.G, Fjortoft R, J. Von Neumann. Numerical integration of the barotropic vorticity equation. Tellus. 1950. 2, 237-254.

4. Hendershott M.C. Long waves and ocean tides. Ch. 10 in: B.A. Warren & C. Wunsch (eds.) Evolution of physical oceanography, MIT Press, 1981. 292-341.

5. Nihoul J.C.J, Ronday F.G. The influence of tidal stresses on the residual circulation, Tellus 29 (1975) 484-490.

6. Dube S.K, Sinha P.C, Roy G.D. The numerical simulation of storm surgers along the Bangla Desh coast, Dyn. Atm. Oceans 9 (1985), 121-133.

7. Johns B, Dube S.K, Sinha P.C, Mohanty U.C, Rao A.D. The simulation of continuously deforming lateral boundaries in problems involving the shallow-water equations, Сотр. And' Fluids 10, 2,(1982) 105-116.

8. Garvine R.W. Estuary plume and fronts in shelf waters: a layer modele, J. Phys. Oceanogr. 17 (1987) 1877-1896.

9. Dowling T.E, Ingersoll A.P. Jupiter's Great Red Spot as a shallow-water system // Atmosph. Sci. 1989. 46, 21, P. 3256-3278.

10. Alcrudo F, Garcia-Navarro P. A high resolution Godunov-type scheme in finite volumes for the 2D shallow water equation, Int. J. Numer. Meth. Fluids 16 (1993) 489-505.

11. Stelling G.S. On the construction of the computational methods for shallow-water flow problems, PhD thesis Delft Univ. of Technology (1983).

12. Wind H.G, Vreugdenhil C.B. Rip-current generation near structures // J. Fluid Mech. 1986. 171,459-476.

13. Ogink H.J.M, Grijsen J.G. Wijbenga A.J.H. Aspects of flood level computations. Int. Symp. Flood Frequency and Risk Analysis, Baton Rouge, USA, 1986 (also Delft Hydraulics Comm. 357).

14. Platzman G.W. Two-dimensional free oscillations in natural basins // J. Phys. Oceanogr. 1972. 2, 2, 117-138.

15. Shokin Y.I, Chubarov L.B. Finite-difference simuletion of tsunami propogation, in (U. Muller et al, eds) Theoretical and experimental fluid mechanics, Springer, Berlin (1980) 599-606.

16. Марчук А.Г, Чубаров Л.Б, Шокин Ю.И. Численное моделирование цунами.17