Одномерные вариационные и полностью консервативные разностные схемы МГД в смешанных эйлерово-лагранжевых переменных тема автореферата и диссертации по математике, 01.01.07 ВАК РФ

ВВЕДЕНИЕ

ГЛАВА. I. ВАРИАЦИОННЫЙ ПОДХОД К ПОСТРОЕНИЮ ДИСКРЕТНЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ГАЗОВОЙ ДИНАМИКИ.И. МАГНИТНОЙ ГИДРОДИНАМИКИ В СМЕШАННЫХ ЭШШРОВО-ЛАГРАНЖЕВЫХ ПЕРЕМЕННЫХ

§ I. Вариационный принцип. Дифференциальные уравнения.

§ 2. Пространственная дискретизация. Дифференциально-разностные уравнения.

§ 3.' Законы сохранения.

§ 4, Простейшая дискретная по пространству модель сплошной среды. Искусственная вязкость

§ 5. Каскадная форма аппроксимации временных производных

Широкий круг вопросов современной науки и техники связан с решением задач газовой динамики (ГД) и магнитной гидродинамики (МГД).

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

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

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

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

При разработке численных алгоритмов первым шагом является выбор способа аппроксимации сплошной среды некоторым дискретным ее аналогом [l-з].

Пространственная сетка может быть эйлеровой, лагранжевой или смешанной эйлерово-лагранжевой (СЭЛ).

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

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

К недостаткам эйлеровых методов следует также отнести обычно сопутствующую им чрезмерно высокую аппроксимационную вязкость, затрудняющую расчеты при наличии контактных разрывов. Для преодоления этих недостатков используется метод коррекции потоков [IIJ и различные его модификации [12-14J . Так в [13 ] идея коррекции потоков реализована на основе сравнения схем первого и более высокого порядков аппроксимации.

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

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

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

Лагранжевы сетки при расчете двумерных задач состоят обычно из четырехугольных ячеек [15-17], поскольку такие ячейки связаны между собой простыми логическими связями. Однако лагранжевы сетки могут быть образованы и из ячеек произвольной формы. Использование треугольных ячеек [16,18,19] , по сравнению с четырехугольными, имеет некоторое преимущество, потому что они всегда выпуклы. Построение разностных схем на лагран-жевой сетке не так просто, как для регулярной эйлеровой сетки [20 ] . Однако несомненным достоинством лагранжевых методов является их большая, по сравнению с эйлеровыми методами, точность расчета контактных границ, обусловленная отсутствием членов, соответствующих конвективному переносу.

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

СЭЛ методы обладают большей гибкостью по сравнению с лаг-ранжевыми и эйлеровыми методами, что позволяет приспособить вычисления к конкретной задаче [23-30] . Использование СЭЛ подхода может давать большие преимущества. Во многих задачах можно исключить за счет выбора СЭЛ разностной сетки движение среды как целого и тем самым уменьшить конвективный перенос. При использовании подвижной разностной сетки могут быть легко выделены границы раздела, а в области, где требуется хорошее разрешение, узловые точки можно "сгустить". Благодаря своей гибкости, СЭЛ методы находят широкое применение.

Помимо методов, основанных на использовании эйлеровых, лагранжевых или СЭЛ сеток, существует довольно обширный класс алгоритмов, которые нельзя отнести ни к одному из рассмотренных выше типов. Это методы, основанные на представлении сплошной среды в виде "частиц" [31-34] или "свободных точек" [35-37] .

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

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

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

Существенный шаг в разработке эффективных разностных схем

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

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

Для упрощения и формализации этой процедуры разработан достаточно гибкий вариационный подход, на основе которого были получены двумерные и трехмерные полностью консервативные разностные схемы, среди которых есть схемы в произвольных криволинейных системах координат [51-58] , а также полностью консервативные разностные схемы на треугольных лагранжевых сетках [59-60] . Примеры численных расчетов сложных двумерных прикладных задач полностью подтвердили высокую эффективность полностью консервативных вариационно-разностных схем [б1-7з] .

Один из наиболее общих и фундаментальных подходов к изучению задач теоретической и математической физики основан на использовании вариационных принципов, среди которых следует упомянуть закон минимальности энергии для стационарных механических систем [74] , закон наименьшего действия для динамических систем с конечным числом степеней свободы [75]. Известны вариационные формулировки уравнений электродинамики [7б], диссипа-тивных явлений [77] и уравнений ГД и МГД [78,79].

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

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

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

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

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

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

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

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

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

Решению первой из поставленных задач посвящена глава I.

Глава состоит из пяти параграфов и приложения. Материал, содержащийся в этой части работы, опубликован в [83-85^). В главе I подробно изложен вариационный подход к построению разностных схем для задач ГД в СЭЛ переменных. Вариационные разностные модели для задач МГД в СЭЛ переменных описаны в приложении.

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

Требование равенства нулю первой вариации функционала действия, аналогичного функционалу действия по Гамильтону для механических систем с конечным числом степеней свободы [75], приводит в отличии от лагранжева случая при использовании СЭЛ переменных к двум уравнениям. Варьирование проводится с учетом дифференциальных соотношений, представляющих собой обобщение на СЭЛ переменные закона сохранения массы, первого начала термодинамики и кинематического соотношения, описывающего процесс переноса лагранжевой координаты. Одно из полученных уравнений -известное дифференциальное уравнение, описывающее изменение импульса. Второе уравнение определяет величину члена7"Jjf » обусловленного наличием конвективного переноса энтропии в уравнении изменения внутренней энергии: через другие параметры течения: скорость, энтальпию и у - PJL А ~ Vt

Здесь % - функция, описывающая характер перемещения опорных л координат (СЭЛ координат) относительно эйлеровых;- оператор дифференцирования по времени при постоянной | - опорной координате; £ - удельная внутренняя энергия; Р - давление; 7" - температура; 5 - энтропия.

Таким образом применение вариационного формализма приводит к уравнению изменения внутренней энергии в некотором специальном виде.

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

В § 3 показана консервативность построенных вариационных дифференциально-разностных моделей, а в § 4 подробно рассмотрена одна из простейших схем рассматриваемого класса.

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

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

Глава П настоящей работы содержит восемь параграфов и приложение.

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

Случай аксиальной симметрии течения, когда H—(0,Hfi0) описан в приложении.

Предложенный метод основывается на специальном представлении, названном "каскадным", производных по времени и пространству соответствующей дифференциальной модели (§ I, § 3).

Опираясь на каскадную форму записи пространственных производных, в § 2 предложена согласованная аппроксимация конвективных потоков магнитного поля, массы, импульса и энергии, не нарушающая баланса по отдельным видам энергии. А в § 3, исходя из каскадного представления временных производных , - согласованная временная дискретизация.

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

Более подробно рассмотрен случай, когда класс полностью консервативных схем ГД и МГД в СЭЛ переменных, имеющих произвольный наперед заданный порядок аппроксимации по времени и пространству на достаточно гладких решениях (§ 4), причем показано, что дискретные аналоги временных и пространственных производных могут выбираться совершенно независимо.

В § 5 рассмотрены два варианта полностью консервативной аппроксимации временных производных. Оба эти варианта приводят к трехслойным разностным схемам, что находится в полном соответствии с результатами, полученными в [81, 82].

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

В § 6 рассмотрено два варианта полностью консервативной аппроксимации конвективных потоков: симметричная аппроксимация и аппроксимация потоков массы и внутренней энергии "вверх по потоку".

На основе этих двух вариантов аппроксимации конвективных членов предложен полностью консервативный вариант метода коррекции потоков (§ 7).

Основная идея коррекции потоков реализована аналогично [13] на основе сравнения схем первого порядка аппроксимации со схемами более высокого порядка.

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

Материал, изложенный в главе П, опубликован в работах [88-86].

Глава Ш содержит четыре параграфа.

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

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

В § 3 обсувдаются вопросы устойчивости.

§ 4 содержит описание расчетов тестовых задач, на которых были апробированы: простейшая вариационная схема, описанная в главе I, полностью консервативная схема с аппроксимацией конвективных членов в уравнениях изменения плотности и внутренней энергии "вверх по потоку", а также полностью консервативный вариант метода коррекции потоков (временная аппроксимация у всех схем двухслойна по всем переменным кроме скорости).

Вычислительные особенности этих схем рассматривались в сравнении с результатами расчетов по ряду других численных методик. Сравнение проводилось со схемой Годунова [92], схемой SHASTX Бориса-Бука с коррекцией потоков [12,II], схемой Мак--Кормака со сглаживанием (схема Жмакина-Фурсенко) [14] и некоторой консервативной схемой, построенной на основе интегро-интерполяционного подхода [2].

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

- Использование полностью консервативной схемы дает существенные преимущества при решении задач о распаде разрывов по сравнению с известными схемами: Годунова и Бориса-Бука с коррекцией потоков (SHASTA);

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

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

- Полностью консервативная схема не приводит к возникновению нефизических осцилляций вычисляемых величин;

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

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

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

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

Результаты, вошедшие в главу Ш, опубликованы в работах [85-89].

ЗЕ £ эе

Таким образом в работе:

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

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

- На основе сравнения полностью консервативных схем первого порядка аппроксимации со схемами более высокого порядка, предложен полностью консервативный вариант метода коррекции потоков в СЭЛ переменных.

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

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

Материалы диссертационной работы докладывались:

- на семинаре кафедры вычислительных методов факультета ВМК МГУ под руководством академика А.А. Самарского;

- на научно-исследовательском семинаре отдела № 3 ИПМ им. М.В. Келдыша АН СССР под руководством академика А.А. Самарского и д.ф.-м.н. А.П. Фаворского;

- на семинаре ВЦ АН СССР под руководством д.ф.-м.н. Ю.Д. Шмыглевского;

- на Всесоюзной школе "Теоретические и прикладные проблемы вычислительной математики и математической физики" (г.Рига, 1982 г.);

- на Всесоюзной школе "Численные методы решения задач математической физики" (г. Львов, 1983 г.). X

Основные результаты диссертации опубликованы в работах [83-88].

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

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

Заключение

Анализ полученных результатов позволяет сделать ряд выводов:

- Использование полностью консервативной разностной схемы дает существенные преимущества при решении поставленных задач по сравнению с известными схемами: Годунова и Бориса-Бука с коррекцией потоков (SHASTX),

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

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

- Полностью консервативная схема не приводит к возникновению нефизических осцилляций вычисляемых величин.

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

- Хорошие результаты дает полностью консервативная схема с коррекцией потоков. Она хорошо воспроизводит контактный раз

-v

5) i V v

Рис. 11 рыв, но для нее характерны недостатки, присуще лагранжевым схемам, а именно, наличие "энтропийного следа" [3] , правда, более слабо выраженного.

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

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

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

1. Самарский А.А. Введение в теорию разностных схем.- М., "Наука", 1971.

2. Самарский А.А. Теория разностных схем.- М., "Наука", 1977.

3. Самарский А.А., Попов Ю.П. Разностные методы решения задач газовой динамики.- М., "Наука", 1980.

4. Рахтмайер Р.Д., Мортон К. Разностные методы решения краевых задач.- М., "Мир", 1972.

5. Covzvni Я., IsaakSon On the solutionsof nonCineax hypexSoCic differential equationsdifferences. Comm.Pure /!рр£- Wotkf , tSb^l6. ой?х U/zndwtf 6. Sisttms of consezircitiOM

6. Cau/s . Comm. Pu?e ///>pt fflatt., П ,№o

7. Яненко H.H. Методы дробных шагов решения многомерных задач математической физики,- Новосибирск, "Наука", 1967.

8. Брушлинский К.В., Морозов А.И., Палейчик В.В. Расчет двумерных нестационарных течений плазмы в каналах,- В сб. "Плазменные ускорители".- М., Машиностроение, 1973, с.251-254.

9. Жмакин А.И., Фурсенко А.А. Об одной монотонной разностной схеме сквозного счета. ЖВМ и МФ, т.20, 4, 1980, с. 10211031.

10. Шулъц У.Д. Двумерные конечно-разностные гидродинамические уравнения в переменных Лагранжа.- В кн. "Вычислительные методы в гидродинамике".- М., "Мир", 1967.

11. Софронов И.Д., Дмитриев Л.В., Малиновская Е.В. Методика расчета нестационарных двумерных задач газодинамики в переменных Лагранжа.- М., Препринт ИПМ им. М.В.Келдыша1. АН COOP, & 59, 1976.

12. Тишкин В.Ф., Тюрина Н.Н., Фаворский А.П. Разностные схемы для расчета гидродинамических течений в цилиндрических координатах.- М., Препринт ИПМ им. М.В.Келдыша АН СССР,23, 1978.

13. Рождественский Б.Л., Левитан Ю.Л., Моисеенко Б.Д. Численное решение уравнений гидродинамики в переменных Лагранжа.-М., Препринт ИПМ им. М.В.Келдыша АН СССР, В 25, 1971.

14. Гольдин В.Я., Калиткин Н.Н., Левитан Ю.Л., Рождественский Б.Л. Расчет двумерных течений с детонацией. ЖВМ и МФ,т. 12, & 6, 1972, с. I606-I6II.

15. J. ЛррС fitys., tSSOsZi.A/z , гъг-гг* .

16. Самарский A.A., Арсенин В.Я. О численном решении уравнения газодинамики с различными типами вязкости,- ЖЕМ и МФ, т.1, 1961.

17. Франк P.M., Лазарус Р.Б. Смешанный метод, использующий переменные Эйлера и Лагранжа.- Сб. "Вычислительные методы в гидродинамике",- М., "Мир", 1967, с. 55-75.

18. Нох В.Д. СЭЛ совместный эйлерово-лагранжевый метод для расчета нестационарных двумерных задач.- Сб. "Вычислительные методы в гидродинамике".- М., "Мир", 1967, с.128-164.

19. Херт С. Произвольный лагранжево-эйлеров численный метод.-Сб. "Численные методы в механике жидкостей".- М., "Мир", 1973.

20. Годунов С.К., Прокопов Г.П. Об использовании подвижных сеток в газодинамических расчетах.- ЖВМ и МФ, 1972, 12, № 2, с 429-440.

21. Ковеня В.М., Яненко Н.Н. Разностная схема на подвижных сетках для решения уравнения вязкого газа.- ЖВМ и МФ, т. 19, № I, с.173-178.

22. Данаев Н.Т. Об одном способе построения криволинейных сеток, сгущающихся в области больших градиентов.- В сб.

23. Численные методы механики сплошной среды",- Новосибирск, 1979, т.10, & 4, с.60-74.

24. Харлоу Ф.Х, Численный метод частиц в ячейках для задач гидродинамики.- В сб. "Вычислительные методы в гидродинамике",- М., "Мир", 1967.

25. Петренко В.Е., Воронцов Е.В. Применение частиц-слоев при расчетах по методу частиц в ячейках.- Численные методы механики сплошной среды. Т.4, Jfe 2, 1973, с. I32-I4I.

26. Анучина Н.Н. Решение нестационарных задач газовой динамики методом "частиц в ячейках".- М., Препринт ИПМ им. М.В.Келдыша АН СССР, £ 101, 1976.

27. Таран М.Д. Метод частиц конечного размера для моделирования многообластных двумерных задач газовой динамики.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, й 177, 1979.

28. Дьяченко В.Ф. Об одном новом методе численного решения нестационарных задач газовой динамики,- ЖВМ и МФ, т.5, & 3, 1965, с. 680-690,

29. Подливаев И.Ф. Некоторые вопросы метода расчета двумерных задач газодинамики на переменных сетках ("Медуза").- М., Препринт ИПМ им. М.В.Келдыша АН СССР, № 52, 1976,

30. Pasta у. я., tidam s. Heuiistic numeztco£ шогк.in some pzog£ems of hydzodt/ncitnccz. Math. TaS&s Aids. Com put, /$?Ъ,/3 .

31. Попов Ю.П., Самарский А.А. Полностью консервативные разностные схемы газодинамики и магнитной гидродинамики.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, Jfc 16, 1969.

32. Попов Ю.П., Самарский А.А. Полностью консервативные разностные схемы.- ЖЕМ и МФ, т.9, № 4, 1969, с.953-958.

33. Попов Ю.П., Самарский А.А. Полностью консервативные разностные схемы для уравнений магнитной гидродинамики.-IBM и МФ, т.10, Jfc 4, 1970, с.990-998.

34. Гольдин В.Я., Ионкин Н.И., Калиткин Н.Н. Об энтропийной схеме расчета газодинамики.- ЖВМ и МФ, т.9, I 6, 1969, C.I4II-I4I3.

35. Майорова О.С., Попов Ю.П. К расчету двумерных уравнений газовой динамики.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, В 55, 1977.

36. Мажорова О.С., Попов Ю.П. О полностью консервативных схемах для двумерных уравнений газовой динамики,- Дифференциальные уравнения, 1979, 15, Jfe 7, с.1318-1331.

37. Маламуд В.Г., Попов Ю.П., Шапиро В.М. Применение полностью консервативных схем к расчету задач гидродинамики и магнитной гидродинамики.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, 18, 1970.

38. Арделян Н.В., Бисноватый-Коган Г.С., Попов Ю.П. Исследование магниторотационного взрыва сверхновой в цилиндрической модели.- Астрономический журнал, 1979, 56, №6, с.1244--1256.

39. Арделян Н.В., Попов Ю.П. Особенности применения численных методов в задаче о магниторотационном взрыве сверхновой.-М., Препринт ИПМ им. М.В.Келдыша АН СССР, № 107, 1979.

40. Косовичев А.Г., Попов Ю.П. Некоторые особенности распространения ударных волн в атмосфере Солнца.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, Я 73, 1978.

41. Гущин И.С., Попов Ю.П. К расчету задач магнитной гидродинамики с учетом фазового перехода.- ЖВМ и МФ, 1977, т.17, № 5, с. 1248-1255.

42. Данилова Г.В., Курдюмов С.П., Попов Ю.П. и др. Расчет сильноточечных разрядов с учетом эффекта вторичного пробоя.-М., Препринт ИПМ им. М.В.Келдыша АН СССР, 6, 1974.

43. Самарский А.А., Дородницын В.А., Курдюмов С.П., Попов Ю.П. Образование Г-слоев при торможении плазмы магнитным полем. -Докл. АН СССР, 216, В 6, 1974, с. 1254-1257.

44. Головизнин В.М., Самарский А.А., Фаворский А.П. Вариационный подход к построению конечно-разностных моделей в гидродинамике.- Докл. АН СССР, т.235, Jfc 6, 1977, с. 1285-1288.

45. Головизнин В.М., Самарский А.А., Фаворский А.П. Об использовании принципа наименьшего действия для построения дискретных математических моделей в магнитной гидродинамике.-Докл. АН СССР, т.246, № 5, 1979, с.1083-1088.

46. Головизнин В.М., Коршия Т.К., Самарский А.А., Фаворский А.П. Двумерные вариационно-разностные схемы магнитной гидродинамики с тремя компонентами скоростей и магнитного поля.-М., Препринт ИПМ им. М.В.Келдыша АН СССР, № 41, 1978.

47. Головизнин В.М., Самарский А.А., Фаворский А.П. Вариационный подход к построению сбалансированных дискретных моделей в многомерной гидродинамике.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, Ш 141, 1980.

48. Фаворский А.П. Вариационно-дискретные модели уравнений гидродинамики.- М., Препринт ИБМ им. М.В.Келдыша АН СССР, П> 159, 1979.

49. Головизнин В.М., Коршия Т.К., Самарский А.А., Фаворский А.П. Вариационные схемы магнитной гидродинамики в произвольной системе координат.- ЖВМ и МФ, т.21, № I, 1981, с.54-68.

50. Головизнин В.М. Трехмерные дифференциально-разностные схемы МГД.- Дифференциальные уравнения, т.17, № 7, 1981,с. 1222-1227.

51. Головизнин В.М., Коршунов В.К. 0 вариационно-разностных схемах газовой динамики на треугольных сетках.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, Ш 17, 1980.

52. Головизнин В.М., Коршунов В.К. 0 двумерных вариационно-разностных схемах МГД на треугольных сетках.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, В 99, 1980.

53. Волкова Р.А., Головизнин В.М. и др. Численное моделированиеобжатия магнитного поля кумулирующим лайнером.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, В III, 1976.

54. Гасилов В.А., Головизнин В.М. и др. Численное решение одной модельной задачи о релей-тейлоровской неустойчивости.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, Jfc 119, 1977.

55. Головизнин В.М., Коршия Т.К. и др. Численное исследование разлета плазмы в магнитном поле.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, В 61, 1978.

56. Гасилов В.А., Головизнин В.М. и др. О численном моделировании релей-тейлоровской неустойчивости в несжимаемой жидкости.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, № 70, 1979.

57. Гасилов В.А., Головизнин В.М. и др. Численное моделирование компрессии тороидальной плазмы квазисферических лайнеров.-М., Препринт ИПМ им. М.В.Келдыша АН СССР, Jfe 71, 1979.

58. Гасилов В.А., Головизнин В.М. и др. Численное моделирование динамики квазисферического жидкометаллического лайнера.-Доклад на П Международной конференции по мегагауссным магнитным полям, Вашингтон, 1979.

59. Гасилов В.А., Головизнин В.М. и др. Модельное исследование квазисферической лайнерной компрессии тороидальной плазмы.-Доклад на IX Международной конференции по УТС, 1980.

60. Гасилов В.А., Коршия Т.К. и др. К расчету уравнений магнитной гидродинамики для среды с конечной электропроводноетью.-М., Препринт ИШ им. М.В.Келдыша АН СССР, Я 6, 1980.

61. Гасилов В.А., Головизнин В.М. и др. Применение метода Ньютона для решения разностных уравнений гидродинамики.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, & 100, 1978.

62. Волосевич П.П. и др. Процесс сверхвысокого сжатия веществаи инициирование термоядерной реакции мощным импульсом лазерного излучения.- Физика плазмы, т.2, J& 6, 1976, с. 883-897.

63. Волосевич П.П. и др. Двумерные эффекты при лазерном сжатии стеклянных оболочек.- Письма в ЖЭТФ, т.24, вып.5, 1976.

64. Гамалий Е.Г. и др. Гидродинамическая устойчивость сжатия сферических лазерных мишеней,- М., Препринт ИПМ им. М.В.Келдыша АН СССР, 117, 1978.

65. Афанасьев Ю.В. и др. Нагрев и сжатие лазерных мишеней.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, # 167, 1977.

66. Морс Ф.М., Фешбах Г* Методы теоретической физики,- М., ГТТИ, 1951.

67. Парс Л. Аналитическая динамика,- М., "Наука", 1971.

68. Ландау Л.Д., Лифшиц Л.М. Теория поля.-М,, "Физматгиз", 1962.

69. Вио М. Вариационные принципы в теории теплообмена.- М., "Энергия", 1975.

70. Se€en(jerz R.X., Whit ham С. В. Vaiiaiionot pit п tipбе$ cn continuum mechanics. Рюс Ray Sos /floS^-iSM

71. Головизнин В.М., Самарский А.А., Фаворский А.П. Вариационный принцип получения уравнений магнитной гидродинамики в смешанных эйлерово-лагранжевых переменных.- ЖВМ и МФ, т.21, Jfe 2, с. 409-422, 1981.

72. Самарский А.А,, Тишкин В.Ф., Фаворский А.П., Шашков М.Ю. Построение полностью консервативных разностных схем для уравнений газовой динамики в эйлеровых координатах на основе операторного подхода.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, J& 63, 1981.

73. Кузьмин А.В., Макаров В.Л., Меладзе Г.В. Об одной полностью консервативной разностной схеме для уравнений газовой динамики в переменных Эйлера.- ЖВМ и МФ, 1980, т.20, Jfe I,с. I7I-I8I.

74. Кузьмин А.В., Макаров Б.Л. Об одном алгоритме построения полностью консервативных разностных схем.- ЖВМ и МФ, т.22,1. 1982, с. 123-133.

75. Коптелова Н.А., Сороковикова O.G. Вариационный подход к получению уравнений акустики в смешанных эйлерово-лагранже-вых переменных для адиабатических течений магнитной гидродинамики.- М., Препринт ИАЭ им. И.В.Курчатова АН СССР, ИАЭ -3568/6, 1982.

76. Головизнин В.М., Сороковикова О.С. О гамильтоновой форме уравнений МГД в смешанных эйлерово-лагранжевых переменных.-М., Препринт ИАЭ им. И.В.Курчатова АН СССР, ИАЭ-3567/16, 1982.

77. Гасилов В.А., Головизнин В.М., Сороковикова О.С. Вариационный подход к построению дискретных математических моделей газовой динамики в смешанных эйлерово-лагранжевых переменных.- М., Препринт ИПМ им. М.В.Келдыша АН СССР, № 35, 1983.

78. Головизнин В.М., Рязанов М.А., Сороковикова О.С. Одномерные полностью консервативные разностные схемы МГД в смешанных эйлерово-лагранжевых переменных.- М., Препринт ИАЭ им. И.В.Курчатова, ИАЭ-3747/16, 1983.

79. Головизнин В.М., Рязанов М.А., Самарский А.А., Сороковикова О.С. Полностью консервативная коррекция потоков в задачах газовой динамики. Докл. АН СССР, т. 274, А^ 3, 1984,с. 524-528.

80. Головизнин В.М., Рязанов М.А., Сороковикова О.С. Об одном классе полностью консервативных схем МГД в смешанных эйлерово-лагранжевых переменных. ЖВМ и МФ, т.24, J 4, 1984,с. 520-533.

81. Головизнин В.М., Рязанов М.А., Сороковикова О.С., Самарский А.А. Разностные схемы газовой динамики со сбалансированными аппроксимациями конвективных потоков.- Препринт ИПМ им. М.В.Келдыша АН СССР, J& 56, 1984.

82. Головизнин В.М., Коршия Т.К. и др. 0 некоторых свойствах вариационно-разностных уравнениях МГД и искусственных дис-сипативных процессах типа "псевдовязкости".- Препринт ИПМ им. М.В.Келдыша АН СССР, J£ 25, 1978.

83. Головизнин В.М. Об одном способе введения искусственных диссипаторов в вариационно-разностные схемы магнитной гидродинамики.- Препринт ИПМ им. М.В.Кедцыша АН СССР, 50, 1980.

84. Годунов С.К. Разностный метод численного расчета разрывных решений гидродинамики.- Матем. сб., 1959, 47(89), с.271--306.

85. Головизнин В.М., Коршунов В.К., Таран М.Д. Об аппроксимации дифференциальных операторов на нерегулярных косоугольных сетках.- Препринт ИПМ им. М.В.Келдыша АН СССР, Jfe 157, 1980.

86. Ланцош К. Вариационные принципы механики. М., "Мир", 1965.

87. Головизнин В.М., Рязанов М.А., Сороковикова О.С. Полностью консервативные дифференциально-разностные схемы газовой динамики в смешанных эйлерово-лагранжевых переменных.-Препринт ИПМ им. М.В.Келдыша АН СССР, й 19, 1982.

88. Гантмахер Ф.Р. Теория матриц. М., "Наука", 1967.