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

Садов, Алексей Павлович

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

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

кандидата физико-математических наук: Садов, Алексей Павлович
город: Москва
год: 2008
специальность ВАК РФ: 05.13.18

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

Автореферат диссертации по теме "Численно-аналитическое моделирование одномерных течений теплопроводного невязкого газа"

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

САДОВ Алексей Павлсши1

ЧИСЛЕННО-АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОДНОМЕРНЫХ ТЕЧЕНИЙ ТЕПЛОПРОВОДНОГО НЕВЯЗКОГО ГАЗА

Специальность 05.13.18 - математическое моделирование, численные методы

и комплексы программ

АВТОРЕФЕРАТ

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

иис(4Ь40е 1

Москва - 2008

003454061

Работа выполнена в Уральском государственном университете путей сообщения на кафедре "Прикладная математика"

Научный руководитель:

-доктор физико-математических; наук, профессор Дерябин С.Л.

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

-доктор физико-математических наук, профессор Филимонов А.М (Московский Государственный Университет Путей Сообщсния(МИИТ)) -кандидат физико-математических наук Козырев О.М. (Российский Федеральный Ядерный Центр - Всероссийский Научно -Исследовательский Институт Технической Физики)

Ведущая организация'

Уральский государственный университет им. А.М. Горького

Защита состоится 03 ^^аг^а. 2008 года в ¿5\3о, ауд. на заседании диссертационного совета Д 218.005.10 по присуждению ученой степени кандидата физико-математических наук при Московском государственном университете путей сообщения (127994, ГСП-4, г.Москва, ул Образцова,

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

Д15)

Автореферат разослан ЗОо^^а 2008 года. Ученый секретарь диссертационного совета кандидат технических наук, профессор

ВВЕДЕНИЕ

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

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

В задачах о получении больших значений температуры и плотности для более адекватного описания возникающих течений необходимо учитывать равновесное излучение. Этот физический процесс описан в книгах Е.И. Забабахина, И.Е. Забабахина1, Я.Б. Зельдовича, Ю.П. Райзера2. В реальных физических экспериментах эффекты лучистой теплопроводности проявляются при очень больших значениях температуры нагретого газа Поэтому в физической и математической литературе при моделировании подобных течений часто полагают температуру фона равной нулю.

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

Математическое моделирование сильного сжатия газа ведется в различ-

'Забабахин Е.И., Забабахин И Б Явления неограниченной кумуляции М Наука, 1988

2Зелъдович Я Б , Райзер Ю П Физика ударных волн и высокотемпературных гидродинамических явлений М Наука, 1966.

ных направлениях.

if = ~Ü ^—Y^j , tt, х{ - const,

где t, Х\ - независимые переменные, вектор ~Ü задает параметры течения газа. При t < tt центрированная волна Римана описывает течение, возникающее при безударном сильном сжатии плоского слоя в идеальном газе, что описано в работе К.П. Станюковича. Е.И. Забабахиным, И.Е. Забабахиным1 указано, что впервые центрированная волна Римана применена Гюгонио и Релеем для описания сжатия плоского слоя газа до сколь угодно большой плотности. В случае цилиндрически- и сферически-симметричных течений особенностью, аналогичной особенности в центрированной волне Римана, обладают автомодельные решения Л.И. Седова3, с помощью которых осуществляется описание безударного сильного сжатия первоначально однородного и покоящегося в цилиндре или шаре идеального политропного газа. Интерпретации этих решений для задач о безударном сильном сжатии газа посвящены работы Я.М. Каждана, И.Е. Забабахина, В.А Симоненко, А.Н. Крайко

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

3Седов Л И. Методы подобия и размерности в механике - М.' Наука, 1981 - 448 с

^Сидоров А.Ф Избранные труды. Математика Механика - М ФИЗМАТ Л ИТ, 2001 - 576 с.

6Долголева Г В , Забродин А В Кумуляция энергии в слоистых системах и реализация безударного сжатия М • Физматлит, 2004

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

В монографии С.П. Баутина6 предложен и в последующем развит (см. монографию Баутина С.П. 7) единый подход к математическому моделированию сильного сжатия газа. Прошедшее с 1997 года десятилетие показало эффективность предложенного подхода. При этом подходе сначала ставятся начально-краевые задачи, описывающие процесс безударного сильного сжатия произвольного, локально-аналитического фонового течения на произвольной локально-аналитической поверхности. Для поставленных начально-краевых задач доказываются теоремы существования и единственности аналитических и кусочно-аналитических решений. Исследуются свойства решений, в том числе, устанавливаются асимптотические законы поведения газодинамических параметров при неограниченном росте плотности. С использованием строго полученных математических свойств исследуемых процессов далее численно и аналитически решаются конкретные задачи, моделирующие процессы сильного сжатия газа, в том числе, для специальных конфигураций сжимаемых объемов, используемых в физических экспериментах, а также при учете дополнительных физических факторов (излучение, теплопроводность) и реальных уравнений состояния. Многие полученные ранее точные решения, в том числе, центрированная волна Римана вкладываются в решения, найденные с использованием методики С.П. Баутина6'7, как частные случаи, при которых обрываются соответствующие ряды и получаются конечные формулы. Этот подход получил дальнейшее развитие при математическом моделировании течений газа в работах С.П. Баутина и его учеников: С.Л. Дерябина, A.JI. Казакова, Ю.В. Николаева, A.B. Рощупкина, Ю Ю. Чернышова, С.А. Ягупова и автора. В том числе, С.П. Баутиным совместно с C.JI. Дерябиным при t > U рассмотрено одно- и двумерное истечение газа в вакуум в различных ситуациях, в том числе и в случае произвольного уравнения состояния р — р1 }{р, S) с аналитической функцией / Обобщение центрированной волны на случай уравнения состояния в виде р = р(р, е) принадлежит

6Баутин С П Математическая теория безударного сильного сжатия идеального газа - Новосибирск Наука, 1997 -160 с.

тБаутин С П Математическое моделирование сильного сжатия газа - Новосибирск Наука, 2007-312 с

С.П. Баутину. С.А. Ягуповым построены течения, обобщающие свойства центрированной волны Римана на случай изэнтропических течений водорода с реальной изэнтропой, полученных в физических экспериментах. Ю.В. Николаевым при решении различных задач численно построены обобщения центрированной волны в цилиндрически- и сферически-симметричных случаях, описывающие безударное сжатие газа до плотностей, в десятки-сотни тысяч раз превышающих исходную плотность однородного покоящегося газа. A.B. Рощупкиным рассчитаны двумерные течения с соответствующей особенностью до плотностей, в 6-8 раз превышающих исходную плотность однородного покоящегося газа. Естественно, такие степени сжатия недостаточны для моделирования требуемого в экспериментах сильного сжатия газа. Тем не менее рассчитанные степени двумерного сжатия больше, чем в бесконечно сильной ударной волне.

Влияние лучистой теплопроводности на процесс кумуляции плоского слоя газа численно было исследовано М.Г. Анучиным 8 с использованием комплекса программ "Тигр".

Математическое моделирование тепловых волн проводилось многими авторами. Для нелинейного уравнения теплопроводности Г.И Баренблатом впервые были приведены примеры тепловых волн, распространяющихся по холодной среде с конечной скоростью, то есть сжимаемость среды и ее движение не учитывалось. Применение методологии характеристических рядов к построению тепловых волн было предложено А.Ф. Сидоровым4, в частности, в одномерном плоскосимметричном случае доказано существование и единственность аналитического решения при специальном виде фронта тепловой волны. С.П. Баутиным9 доказано существование и единственность аналитической тепловой волны у нелинейного уравнения теплопроводности как при задании аналитического фронта тепловой волны, так и при задании краевого температурного режима. С.П. Баутиным и A.A. Елисеевым доказаны существование и единственность аналитического решения у нелинейно го уравнения теплопроводности в многомерном случае при заданном произвольном аналитическом краевом режиме с отличной от нуля в момент t = 0 производной по времени. С.С. Титовым при рациональных значениях константы с, являющейся показателем степени в нелинейном коэффициенте теплопроводности, в плоско-, цилиндрически- и сферически-симметричных случаях в виде специальных рядов построены решения задачи о тепловой волне с заданным фронтом, обладающие на фронте тепловой волны конкретными особенностями.

8 Анучин М Г Влияние теплопроводности на неограниченное безударное сжатие плоского газового слоя // Прикладная механика и техническая физика - 1998 -T 39, N>4 -С 25-32.

9Баутин С П Аналитическая тепловая волна M Физматлит, 2003

С.С. Титовым приведены примеры точных решений нелинейного уравнения теплопроводности в виде специальных многочленов. С.С. Титовым также получены решения нелинейного уравнения теплопроводности для симметричного случая в виде сходящихся логарифмических рядов. С П. Кур-дюмовым в автомодельном виде построено плоскосимметричное течение теплопроводного невязкого газа, в котором тепловая волна распространяется по холодной среде с конечной скоростью по фону, имеющему строго положительную температуру. Более общее течение теплопроводного невязкого газа с аналогичным свойством построено С.П. Баутиным7. Также С.П. Баутиным построены тепловые волны, распространяющиеся в теплопроводном вязком газе при задании фронта их движения в случае, когда показатель степени температуры в коэффициенте теплопроводности равен единице. Аналог центрированной волны Римана в теплопроводном невязком газе описан в одномерном случае С.П. Баутиным и Ю.Ю. Черны-шовым, в двумерном - Ю.Ю. Чернышовым. Для идеального теплопроводного газа в случае холодного фона (Т0 = 0) Е.И. Забабахиным и В.А. Симоненко показано, что в бесконечно сильной волне реализуется течение с изотермическим скачком. В случае 7о > 0 в работах Я.Б. Зельдовича и Ю.П. Райзера2, Л.Д. Ландау и Е.М. Лифшица, Б.Л. Рождественского и H.H. Яненко10 установлено, что при малых значениях скорости D движения бегущей волны имеет место непрерывный переход, а при больших D - изотермический скачок. Бегущие волны в теплопроводном газе с уравнениями состояния, учитывающими равновесное излучение, были ранее рассмотрены для холодного фона с То = 0 в работах В.А. Белоконя и в уже упомянутой книге2. Также найдены два режима: при малых D изотермический скачок, при больших D - непрерывный переход.

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

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

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

г = х — Б • Ь,

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

Целями работы являются:

1. Доказательство существования плоско-, цилиндрически- и сферически-симметричных тепловых волн при учете равновесного излучения с использованием математической модели теплопроводного невязкого газа в окрестности фронта тепловой волны Доказательство проведено в двух случаях: 1.1) когда задан аналитический закон движения фронта тепловой волны, 1.2) когда задан аналитический краевой температурный режим на некоторой аналитической поверхности.

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

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

4. Приближенное построение полей течений во всей рассматриваемой области: от сжимающего поршня до фронта тепловой волны; от точки краевого температурного режима до фронта тепловой волны.

Научная новизна работы заключается в следующем.

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

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

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

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

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

Теоретическая ценность работы состоит в следующем.

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

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

С помощью конечных отрезков рядов смоделированы тепловые волны с заранее заданным фронтом.

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

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

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

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

Основные результаты, выносимые на защиту

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

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

4. В виде бегущих волн численными расчетами решений системы обыкновенных дифференциальных уравнений смоделирована тепловая волна, идущая по нехолодному фону. Расчетами установлены значения констант Дь .Оь Х?2 и наличие трех режимов распространения фронта бегущей волны в зависимости от величины Б - скорости движения бегущей волны: при

Dq < D < Di - непрерывный переход; при Di < D < D<¿ - изотермический скачок; при Di< D снова непрерывный переход.

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

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

IX Всероссийский съезд по теоретической и прикладной механике, Нижний Новгород, 2006 г.;

Международная конференция "Лаврентьевские чтения по математике, механике и физике", посвященная "105-летию со дня рождения академика М.А. Лаврентьева, Новосибирск, ИгИЛ, СО РАН, 2005 г.;

Международная конференция "VII Забабахинские научные чтения", Снежинск, РФЯЦ-ВНИИТФ, 2005 г.;

Международная конференция "VIII Забабахинские научные чтения", Снежинск, РФЯЦ-ВНИИТФ, 2007 г.;

21-ая Всероссийская конференция "Аналитические методы в газовой динамике САМГАД-2006", ИгИЛ СО РАН, ИПМ РАН, Санкт-Петербург, 2006 г.;

37-ая региональная молодежная конференция "Проблемы теоретической и прикладной математики", Екатеринбург, УрО РАН, 2006 г.;

Международная научная конференция "Устойчивость, управление и моделирование динамических систем", посвященная 75-летию со дня рождения И.Я. Каца, Екатеринбург, УрГУПС, 2006 г.;

Международная научно-практическая конференция "Снежинск и наука - 2006. Трансфер технологии, инновации, современные проблемы атомной отрасли", Снежинск, СГФТА, 2006 г.;

Молодежная научно-практическая конференция "Молодые ученые -транспорту", Екатеринбург, УрГУПС, 2006 г.,

Молодежная научно-практическая конференция "Молодые ученые -транспорту", Екатеринбург, УрГУПС, 2007 г.

Публикации. По теме диссертации опубликовано 14 работ [1-14]. Из них - восемь в виде статей [X—8] (в том числе, две - в журналах из списка ВАК [1-2]), шесть в виде тезисов. Результаты докладывались на десяти научных конференциях

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

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

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

В главе 1 рассматривается система нелинейных дифференциальных уравнений с частными производными, решения которой описывают плос-KO-, цилиндрически- и сферически-симметричные течения теплопроводного невязкого газа Построены одномерною тепловые волны в невязком газе в случае задания закона движения фронта тепловой волны по холодному однородному покоящемуся газу, х — a(t) в случае v = 0 и г = a(t), где fv+1 •

г = J x\. в случае ^ = 1,2. Также построены тепловые волны, когда

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

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

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

1 Т4

р = RpT + - сгТ4 , е = cvT + с— ; а = const > 0. (1)

3 р

Здесь: р - давление; р - плотность; Т - температура; е - внутренняя энергия; R - газовая постоянная; а - константа, связанная с постоянной

Стефана-Больцмана а* соотношением а = 4ст*/с,; с, - скорость света в вакууме; с„ - удельная теплоемкость при постоянном объеме. Коэффициент лучистой теплопроводности такой

(7-1) Р

где а, - положительная константа, зависящая от выбора системы единиц, 7 — 1 = fi/cu > 0 — показатель политропы идеального газа.

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

Pt + ирх + рих = О,

Щ + иих + ± [ТРх + (р + сTiT3) Тх] = 0, (3)

^ (р + а2Т3) (Т( + uïi) + (7 - 1)Т (р + ^Т3) = ^ (Çg) •

В этой системе: t - время; х - пространственная координата; и - скорость газа;

4 TqÔ 2сте,а„Тдо

crj = -а—-Ü-; а2 = 3(7 - 1)cti; К0 = »-

о лроо /шооРоо^оо

гДе Poo^ooi^ooi^oo ~ постоянные величины, используемые при введении

безразмерных переменных. Искомая тепловая волна непрерывно премыкает

к однородному холодному покоящемуся газу со следующими значениями

параметров

P(M)UaM = 1' u0M)lï=e«) =°» T(«.t)l,=e(i) = 0- (4)

Это состояние газа далее называется фоновым течением.

Некоторая заданная аналитическая в окрестности t ~ 0 функция a(t), удовлетворяющая неравенству

а'(0) ф 0, (5)

определяет фронт тепловой волны в виде х = a(t). Ha этом фронте тепловая волна с Т > 0 непрерывным образом стыкуется с фоновым течением.

Делается замена независимых переменных

Г 2 = X - a(t), ^

I т = i,

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

(Г)

т = т.

При заменах (6), (7) система (3) переходит в систему

" 3?рт + {и-а'(т))р1 + рщ = 0, 3¿V + (и - о'(т)) Щ + ^р [Тр( + {р + ^Т3) Те] = о, З£3 (р + а2Т3) (3етт + [и - а'(т)] Т() + зеь ~ 1)Г {р + щ =

Ко

- ^ - +

W(r,0=EWfc(r)g, IV

(8)

Условия (4) при заменах (6), (7) переходят в условия

р(г,0|£=о = 1. «(г,е)1{=о = 0, Т(т,О1«=о = 0. (9)

Решение задачи (8), (9) представляется в виде ряда по степеням £

, W = (р,и,Т). (10)

С=о

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

Сходимость ряда (10) доказана методом мажорант.

В параграфе 1.2 рассматривается задача о моделировании цилиндрически- [и = 1) и сферически-симметричных (г/ = 2) тепловых волн, распространяющихся по холодному однородному покоящемуся газу при условии задания закона движения фронта в окрестности точки, имеющей

10+1

ненулевой радиус г — а(1), а(0) = г0 > 0, г = Л ^ х2, и = 1,2.

Система уравнений газовой динамики в этом случае имеет вид

' р{ + ирт+р (иг + иЩ = О, щ + иит + ± [ТРг + (р + агТ3) Тг} = О, (р + а2Т3) (Г, + иТг) + (7 - 1 )Т (р + о<Г*) (иг + V*) =

(И)

к0

Р тт "г р лг ' г р т

Искомая тепловая волна также непрерывно стыкуется с фоновым течением (см рис. !)•

Р(Х,«) |г=вМ = 1, <г, I) |г=в(() - 0, Г (г, 0 |г=а{4) = 0.

(12)

0 'о

г=аф

л- <•*

А\

,-Л

г=а(1)

и=0

<А \

0 'о

0 го

Рис. 1

Для решения задачи (11)—(12) делаются замены (6)-(7), где вместо переменной х берется переменная г. В результате этих замен задача (11)-(12) переходит в задачу

\ 3еРг + а(т)] + [£3 + а(г)] [и - а'(г)} р(+ +р [ис (£3 + а(т)) + 3еМ = 0,

3еит + (и - а'(т)) щ + ± [Гр€ + (р + ауТ3) Ц] = 0,

3£3Р2 (£3 + а(г)) (р + а2Т3) (3(2Тт + [и - а'(т)] Т() + +3£У (7 - 1 )Т (р + ^Г3) (е3 + а(т)) + 3£2™] = = КО [(£3 + 0(т)) £рТ% - (£3 + а(г)) 2рТ3Т(+ I + а(т)) 3рЦрТ* - (е + а{т)) + 3?>Т3Т?] .

(13)

р(г,Ок=о = 1, «(г,01€=о = 0, Г(т,01с=о = 0.

(14)

Решение задачи (13), (14) представляется в виде ряда по степеням £

оо ¿к

\У(т,0= — (р,и,Т), (15)

к=0

значения производных, выводящих с фронта

где ^¿(т) = &

иС- £=0 тепловой волны.

В диссертации доказано, что формальный ряд (15), решающий задачу (13)—(14), строится однозначно и его коэффициенты являются аналитическими функциями от т. Сходимость ряда (15) доказана методом мажорант.

В параграфе 1.3 рассматривается задача о построении плоскосимметричной тепловой волны при специально заданном краевом температурном режиме. Тепловая волна, распространяющаяся по холодному однородному покоящемуся газу, строится в случае, когда прогрев среды осуществляется на движущейся за фронтом волны поверхности х — Ь(1) по заданному закону

т&х)\я=т = №, т = о, /'(о) = Л > 0- (16)

При этом, по заданной функции /(£) законы движения нагреваемой поверхности х = &(£) и фронта х = а(Ь) определяются в процессе построения решения задачи Таким образом, решается задача (16), (17), (18)

Л + ирх + рих = О,

Щ +иих + ± [Трх + (р + <тхТ3) Тх\ = 0, (17)

к [р + <72Т3) (Я + иТх) + (7 - 1 )Т (р + а{Г3) их = ко^

р(М)|*=а(0 = 1. ЦМ)|х=а(*) = °> Т(х,{)\х=ащ= 0, (18)

с неизвестными пока функциями а(£), Ь(Ь).

Далее делаются замены переменных вида (6)-(7) и система (17) переходит в систему

3?рт + {и-а\т))р1 + рЧ = 0, (19)

3?ит + (и- а'(т)) Ч + ^р [ТР(. + (р + ^Т3) Щ =

З£3 (р + <т2Т3) (зетт + [и - а'(т)} Т() + 3^(7 - 1)Т (р + ахТ3) щ = = Ко

гр2 оггЗ ггЗ о'/'2

Условия (18), (16) при заменах (6)~(7) переходят в условия

P(r,e)|f=o = l, «(т,0|£=о = О, Г(г,01{=о = 0. (20)

= (21)

В силу того, что а(0) = 6(0) = 0 задача (19)—(21) в точке t — 0 имеет особенность, поэтому в дальнейшем полагается, что

b(t) = a(t) - a3(i).

x=b(t)sa(t)-a3(t) x=a(t)

Рис. 2

Тогда условия (21) перепишутся в существенно более простом виде, повторяющем вид краевого условия из книги9:

тОЫ-а«) = /(*)• (22)

Для построения такой тепловой волны, как и в книге9, меняются ролями искомая температура Т и независимая переменная В результате задача (16), (17), (18) переходит в задачу

3£2 (Mr - &Яг) + [u - o'(t)] рт + рит = 0, 3£2 Ыт ~ 6«г) + [« - e'(*)] UT + Yp \-Трт + (р + СТ1Т3)] = °>

(23)

З£3 (р + а2Т3) & (3e2ft - u + a'(t)) - 3Ç3$(7 - 1 )Т (р + <пТ3) ит = / тз 2Тз тз ЗГ2 \ = «О Ç—ÇTT +-£г + Ç-JÇTPT - £-6г •

VP Р Р Р )

e(i,T)fr=0 = 0, u(i,T)|r=0 = 0, p(i,T)|r=0-l. (25)

Теорема. Задача (23)-(25) имеет единственное аналитическое реше-

ние

Для доказательства теоремы делается еще одна замена переменных

v = T; у — Т — /(f).

Решение строится в виде формальных рядов (26) с аналитическими коэффициентами в виде

Z(y,«0=EZ*(y)F. z=(p,u,a Zk(y)=0 (26)

k=o av

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

Поскольку соответствующие якобианы отличны от нуля, решение (26) восстанавливается в физическом пространстве с помощью обратных замен при малых t в области между х = b(t) и х — a(t) (см. рис. 2).

В главе 2 построены бегущие волны, распространяющиеся по нехолодному фону с То > О

Рассматривается система уравнений газовой динамики и условия Гюго-нио для нее в случае теплопроводного невязкого газа, когда коэффициент теплопроводности имеет вид (2) и, следовательно, положителен при Т > 0. Рассмотрение ведется в случае нестационарных плоскосимметричных течений. Для полноты изложения в диссертации также приведены подробные выкладки, проделанные С.П. Баутиным и показывающие эквивалентность различных видов системы уравнений газовой динамики и условий Гюгонио для рассматриваемого случая теплопроводного невязкого газа. В частности, условия Гюгонио приведены как в традиционно используемом виде10, так и виде, когда все значения параметров газа за волной выражены как функции от температуры фона и скорости движения волны. Именно этот вид условий Гюгонио используется в диссертации.

Численным построением решений соответствующей СОДУ описываются бегущие волны, распространяющиеся по нехолодному фону с Tq > 0. В рамках одной термодинамической модели в зависимости от D - скорости распространения фронта бегущей волны - выявлены три режима. В первом случае, когда скорость D невелика, переход от постоянных значений параметров газа по разные стороны от фронта бегущей волны осуществляется непрерывно. При увеличении D возникает изотермический скачок: плотность и скорость газа претерпевают сильный разрыв, а температура меняется непрерывно. При дальнейшем увеличении скорости D бегущей волны изотермический скачок исчезает и течение газа снова становится непрерывным. Причем, при неограниченном увеличении D температура и

10Рождественский Б Л , Янекко Н Н Системы квазилинейных уравнений и их приложения к газовой динамике М Наука, 1968

скорость газа также растут неограниченно, а значение плотности газа за бегущей волной стремится снизу к значению 7ро, где ро - плотность газа перед бегущей волной. Этот результат совпадает с полученным ранее в работе2. Исследован случай, когда температура фона стремится к нулю Т0 —> 0. Показано, что в пределе первый режим исчезает и получаются два режима, ранее описанных в работе2.

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

В параграфе 2.1 строится математическая модель бегущих волн в теплопроводном невязком газе. Рассматривается система уравнений с частными производными10

д{ри) д{р + ри2)

иг+ " м дх " ' -°>

(27)

д т

(ре + р^р)

&Н

е + 1 +

= о,

решения которой описывают плоскосимметричные течения теплопроводного невязкого газа.

Далее приводятся система уравнений с частными производными и условия Гюгонио в случае конкретных уравнений состояния (1) и коэффициентом теплопроводности (2). Традиционным образом вводятся безразмерные переменные и получается система уравнений (3), эквивалентная (27).

Для получения СОДУ, описывающей бегущие волны, в системе (3) делается замена переменных

т = г = х-И ^ (28)

и полагается д/дт — 0 (штрихом далее обозначена производная по г): {и - £>)р' + ри' = 0,

(« - Я)и' + ^ [Тр' +(р + агТ^Т] = 0,

(р + а2Т5)(и - И)Т + (7 - 1 ){рТ + <т{ГА)и' =

(тгт" _ ТгТ'р' 3Г2(Г)-

(29)

= «о - + р

Система (29) имеет три первых интеграла, которые с учетом условий при г -4- +оо

р|+оо = 1, "1+00 = О, Т|+00 = 1, Г |+оо=0 (30)

записываются в следующем виде р{Э - и) = Д

_ „ Т3<1Т -

Если в системе (31) при 2 ±оо положить йТ/йг = 0 или взять «о — = 0, то получаются традиционные условия Гюгонио для нетеплопроводного невязкого газа10, которые также можно получить из интегральных законов сохранения.

Далее показывается, что условия Гюгонио при учете уравнений состояния (1), коэффициента теплопроводности (2) эквивалентны условиям Гюгонио в традиционном виде, полученным10 для системы (27) в случае теплопроводного газа.

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

у У = 9±(П (32)

где функции 5±(Т), р±(Т) определены на полуинтервале (0, Г*] и являются двузначными

9±(Г) = 7(721)Д4СП + (7 - 1) (1 + и±(Т)--ПГ - а[£> - и±(Т)]Т4 + £> (1 + ¡а);

Р±(Т) =

Г>-и±(Г)'

(33)

(34)

а функция и±(Т) имеет вид:

и±(Т)

= чР2 + ^Г4 - (1 + ) ± у^(г) 27Я

где

Ж(Т) = [^-Т4 - (7Р2 + 1 + - 47£>2Г.

Функция \У(Т) положительна на интервале (0;Т,], где Т„ > То = 1. Поэтому функции р±(Т), и±(Т), Ч±{Т) определены от 0 до Т».

Для решения дифференциального уравнения (32) принципиальное значение имеет расположение корней функции д±(Т). При этом, значение То = 1 всегда является корнем этой кривой, лежащим на нижней ветви. Причем, <7_(1) = 0, д_(Т.) = 9+(Г.) и 9_(Г) < 9+(Т) при 0 < Т < Т*. Предельные при 2 —V +оо значения плотности р = 1 и скорости газа и — 0 принадлежат нижним ветвям соответствующих кривых р±(Т), и±(Т) при Т = 1, поэтому нулевое значение теплового потока при Т = 1 также принадлежит нижней ветви кривой <7±(Т). На рис. 3 приведена зависимость ?±(Т) при £> — 5 и при следующих значениях используемых констант: 7 = «о — 1,4842; <71 = 0, 2366 и при этом корень Т = 1 лежит на нижней ветви, а корень Т = Т\- на верхней.

100

50 0

-50

-100 __

0 12 3 ТТ 1*

Рис. 3

Поскольку бегущая волна моделирует ударное воздействие на теплопроводный невязкий газ, то скорость движения О должна быть больше со -скорости звука в однородном покоящемся при 2 = +оо нехолодном газе При указанных выше значениях 7, «о, ^х значение со & 1,00757, поэтому далее полагается О > Ио — со- При Б > Бо на кривой д±(Т) появляется второй корень Т = Х^ > 1 Сначала он располагается на нижней ветви (?-№) = 0), но с ростом I) он переходит на верхнюю ветвь {д+(Т\) — 0), а затем вновь возвращается на нижнюю. В следующем параграфе показывается, что расположение корня Т — Т\ на разных ветвях кривой <?±(Т) приводит к принципиально различным течениям газа: без изотермического скачка или с ним.

В параграфе 2.2 численно строятся и анализируются бегущие волны, распространяющиеся по нехолодному фону с То > 0. Они определяются при численном решении дифференциального уравнения (32) с начальным условием

ги = т°, (36)

2Х

где значение константы Т° должно быть выбрано в интервале (1; Tj). Например, Т° = (1 + Tj)/2. При этом в уравнении (32) и в зависимостях (33)-(35) в качестве нижнего индекса всегда берется "минус", а расчет ведется от точки z — О как в сторону возрастания z, так и в сторону убывания г. Температура газа в бегущей волне определяется при численном решении задачи Коши (32), (36), после чего с помощью формул (34)-(35) определяются скорость и плотность газа в бегущей волне во всей области течения ОТ Z — —00 ДО 2 = +00.

Численными расчетами установлено, что в случае а ф 0, «о ф 0, То > О имеется три режима бегущей волны: непрерывный переход при Dq< D < Di « 1,366; изотермический скачок при Di < D < £>2 и 13,675; вновь непрерывный переход при D > D^. Конкретные численные значения D\ и Г>2 приведены для указанных выше значений 7, ко и о\.

Пусть D0 < D < Dv Тогда корень Т = Т\ > 1, как и корень То = 1, лежит на нижней ветви: g-(Ti) = 0. Решение задачи Коши (32), (36) в этом случае непрерывно и все зависимости T(z),u(z), p(z) монотонно убывают от значений Т = Т\, и = u_(Ti), р = р-{Т\) при г -> —со до значений Т = Т0, и — 0, р = 1 при г +оо.

При D\ < D < D4 корень Т = Т\ переходит на верхнюю ветвь кривой q±(T): q+(Ti) = 0 и возникает второй режим. В этом случае при расчете от г = 0 в сторону уменьшения г уже при некотором конечном значении z — z\ функция T(z) принимает значение Т = Т\. В точке z = z\ для величин q, рай необходимо ставить условие сильного разрыва- перехода от значений q = qi — q-{T\), р = рг — Р-{Т\) и и = щ = w_(Ti) к значениям q = <71 = q+(Ti), p = pi = p+{T\) и и = u\ = u+(Ti) соответственно. При этом плотность и скорость газа в точке z\ меняются скачком, а температура в этом течении меняется непрерывно. Такие скачки в газовой динамике называются изотермическими. За изотермическим скачком все газодинамические параметры постоянны: р — pi, и = ui, Т = Т\.

При дальнейшем увеличении D, то есть при D > £>2, корень Ti лежит на нижней ветви; бегущая волна, как и в случае Do < D < D\, передает размазанный ударный переход со своими значениями параметров газа при 2 = -оо и значениями (30) при г = +оо. В третьем режиме при дальнейшем увеличении D температура и скорость газа растут неограниченно, а значение плотности за бегущей волной стремится снизу к значению 7ро, где ро - плотность газа перед бегущей волной.

Далее в параграфе проведены расчеты со значениями То, уменьшающимися вплоть до То = 10~6. Установлено, что при То —> 0 значения Dq и D\

также стремятся к нулю. Но значение £>2 стремится к величине

^ №+ 4.1(4+ .-Ш*.

а 1713 ] V 7— 1

Следовательно, в пределе при То -> 0 три режима за счет исчезновения первого режима непрерывно перейдут в два режима для бегущих волн, распространяющихся по холодному фону10.

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

Численное моделирование осуществляется двумя способами:

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

2. Численными расчетами полей течений с помощью одной разностной схемы.

Вначале приближенно определяется область применимости рядов, построенных в параграфе 1.1. Доказанные теоремы обеспечивают сходимость этих рядов в окрестности точки (т = 0,2 = 0). Метод мажорант, с помощью которого доказаны эти теоремы, позволяет получить оценки радиуса сходимости рядов. Однако, как показывает практика, эти оценки являются заниженными. Для выявления реальной области сходимости а, следовательно, и области применимости построенных рядов проводится сравнение начальных отрезков рядов с точным решением. Точное решение получается в плоскосимметричном случае, когда тепловая волна распространяется с постоянной скоростью О, то есть х — В1. В этом случае тепловая волна зависит только от одной переменной г = х — БЬ и является бегущей волной. Она строится не только с помощью бесконечных рядов из параграфа 1.1, но и численно при решении обыкновенного дифференциального уравнения. Но в отличие от параграфа 2.2, у уравнения (32) в данном случае имеется особенность на фронте тепловой волны при Т — 0. С помощью рядов эта особенность раскрыта явно в аналитическом виде, что позволяет поставить начальные данные не в особой точке. Численное решение задачи Коши, поставленной не в особой точке, для обыкновенного дифференциального уравнения в современной научной практике считается построением точного решения Сравнение численно построенного решения с начальными отрезками рядов дало возможность оценить область применимости используемых отрезков рядов.

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

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

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

Во-первых, начальные данные поставлены на фронте х = а(£), то есть на прямой 2 = 0. Но поскольку на прямой г = 0 имеют место бесконечные градиенты газодинамических параметров, то от нее отступают на некоторую прямую г = го < 0 (|го| мал). Именно на этой прямой г — го и ставятся граничные условия для всех параметров газа, которые, в свою очередь, вычисляются с помощью конечных отрезков рядов.

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

Расчеты проводятся для тепловых волн с заданным фронтом а(£) = +М2 (Л = 0, |,1,10).

Для задачи о моделировании тепловой волны, когда фронт движется с постоянной скоростью а'(£) = 1 (А = 0), имеется точное решение - бегущая по холодному фону тепловая волна. Бегущая по холодному фону тепловая волна с а(Ь) — I. позволяет протестировать как используемую разностную схему, так и приближенное решение, построенное с помощью конечных отрезков рядов.

В случае ускоренного движения фронта: а(£) = t + М2, а'(£) = 1 + 2А£ (А = 1,10) - в главе производится сравнение приближенного решения, полученного с помощью конечных отрезков рядов, с численным решением, построенным с помощью разностной схемы.

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

В параграфе 3.1 используются бегущие волны для оценки области при-

менимости построенных аналитических решений в случае закона движения фронта тепловой волны вида а(Ь) = ОЬ.

Аналитические решения строятся в виде ряда

к=О

, W = (р,и,Т), (37)

е=о

где £ — \fz. г = х — a(t). В данном случае £ = л^ж — Dt и, поскольку a'(i) = D = const, коэффициенты Wjt (т) ряда (37) являются константами, то есть не зависят от г: Wjt (г) =

Тогда начальные отрезки ряда (37) имеют такой вид

u(t,x)f*-b -^Vx-Dt, (38)

7 V «о^

FTn ,_

V Ло

если учтены слагаемые с £ и следующий вид

пм) --^ггъ-

если учтены слагаемые по £2 включительно.

Бегущие волны описываются решениями СОДУ (29), приведенной в гл 2

В параграфе 2.1 решение этой СОДУ (29) с помощью нахождения первых интегралов сведено к решению одного обыкновенного дифференциального уравнения (32). Условия непрерывного премыкания бегущей тепловой волны к холодному однородному фону имеют вид

р(т,г) |,=о = 1, и(г,2)|,=0 = 0| Т(т, г)|2=о = 0. (40)

Но тогда для уравнения (32) точка г = 0 является особой (коэффициент перед производной равен нулю). Поэтому при численном решении уравнения (32) делается отступ влево в точку го < 0, в которой по формулам

(39) вычисляются значения температуры газа. После этого расчетом от го в сторону убывания г производится численное построение решения уравнения (32) с полученным начальным данным в точке го- Это решение в последующем называется точным решением.

На рис. 4 представлены следующие графики функции Т = Т{г): 1 -точное решение, то есть численное решение задачи Коши с данными в точке г = го, 2 и 3 - приближенные решения, полученные по формулам (38) и (39) соответственно

Рис. 4

Для определения границ областей сходимости рядов (37) при о(4) = 1Н сравнивается точное решение с конечными суммами (39). Будем полагать, что отрезок ряда с тремя членами хорошо передает точное решение, если выполняется следующее условие

Т1

к=О

<6.

(41)

В диссертации берется 5%-ная величина относительной погрешности, то есть 5 = 0,05. С помощью этого критерия определяются границы применимости приближенного решения, то есть определяется значение г такое, что при |г| < г выполняется неравенство (41).

В случае произвольного закона движения тепловой волны х = а{Ь) ряд (37) зависит от двух независимых переменных т, В этом общем случае областью сходимости ряда (37) является некоторая окрестность точки (т = 0,£ = 0). В рассматриваемом случае: о(4) = БЬ - ряд (37) зависит только от одной переменной г и сходится в окрестности точки г = 0, которая задается неравенством < г. В физическом пространстве переменных х, Ь последнее неравенство задает следующую неограниченную полосу {(я,2): £>г-г<£<г + (см. рис. 5а).

Рис. 5а

Рис. 5Ь

На рис. 56 показана область восстановления составного течения в физическом пространстве. В правой полуполосе х > располагается однородный холодный покоящийся газ. В левой полуполосе х < располагается построенное приближенное течение. На фронте тепловой волны а: = ОЬ это течение непрерывно стыкуется с однородным холодным покоящимся газом.

Для сравнения численно построенного точного решения и приближенного решения (39) в диссертации приведены соответствующие таблицы и графики.

В параграфе 3.2 проводится приближенное построение закона движения сжимающего поршня и краевого режима на прямой х — 0 с помощью конечных отрезков рядов, а также определение всех газодинамических параметров на линиях х — 0 и х = хр(£) (траектория движения поршня). В качестве закона движения фронта тепловой волны берутся функции а(г) = Ь и а(г) = * + и2, где Л = 1,10.

Исследование начинается с построения траектории движения поршня х — хр{{), порождающего данную тепловую волну. За траекторию движения непроницаемого поршня принимается траектория движения частицы газа, которая в начальный момент времени при £ = 0 находится в точке х = 0, то есть в точке, из которой начинают движения сжимающий поршень и фронт движения тепловой волны. Для этого ставится задача Коши

где в качестве х) берется ряд (37) при заданном фронте в виде х =

= «(г, х), ж|(=0 = 0,

(42)

и(£,а:) = —

В частном случае при а(*) = I траектория движения поршня найдена в виде следующей неявно заданной функции

£ = 1 -75-111 ¥ Ьл

Ь2$/{х - ¿)2 + 2- 2

+1Г1п

+ 1 - л/3)(1 + л/3)

+ 1 + л/3)

Л3/--

ТУ «о

Кроме этого, задача Коши с правой частью (43) решается численно при упомянутых выше функциях а{£).

Найденная зависимость х = ¡гр(£) подставляется в формулы (43), (44) и в конечные отрезки рядов, задающих плотность и температуру:

X=Хр[1)

= 1 - ^

УУко[аЩ

(44)

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

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

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

В случае a(t) — t функция T(t,x)\x=0 — f(t), задающая краевой температурный режим на прямой х = 0 имеет вид

T(t,0) = Щс + dy/i), с, d- const,

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

u(t,0) = \/i(p + q\/t); р{г,0) = 1 + \/1{р+еЩ, (45)

p,q,c — const.

В случае a(t) = t + Xt2 функция T(t,x)|г_0 = f(t), задающая краевой температурный режим при х = 0, имеет вид

T(t, о) = П + 2Щ1 + Xt) + я^^Х)^) .

а скорость и плотность газа на прямой х — 0 такие:

1+ ? f i^W^i7'

««'°>=ffc^r^f (46)

Тем самым получены краевые режимы на прямой а; = 0, реализация которых порождает тепловую волну с заданным фронтом и определены особенности этих краевых режимов в момент времени i = 0. Эта информация является рекомендацией физикам для постановки соответствующих физических экспериментов Но при этом надо иметь в виду, что из формул (45), (46) следует, что через прямую х = 0 осуществляется приток газа.

В параграфе 3.3 проводится приближенное построение полей течений газа во всей рассматриваемой области в двух случаях: 1) от фронта тепловой волны до траектории движения сжимающего поршня х — хр(1.)\ 2) от фронта тепловой волны до прямой х — 0. В обоих случаях в качестве закона движения фронта тепловой волны берутся функции a{t) = t и a(t) — t + Xt2, где Л = 1,10.

Область применимости построенных приближенных решений определена в параграфе 3.1 и имеет вид полуполосы {(x,i): Dt — г < х < Dt} (см. рис. 5Ь), определяемой заданным законом движения фронта тепловой волны х = Dt

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

действовать с момента времени í = 0, то область применимости ограничена снизу этим моментом времени £ = 0. Сверху область применимости ограничена некоторым моментом времени £ = > 0, когда пересекаются линии краевых режимов с левой границей указанной полуполосы

В первом случае - это момент времени, когда траектория движения поршня х = хр{Ь) пересекает границу полосы сходимости (см. рис. 6), во втором случае = г/Б - момент времени, когда граница полосы сходимости пересекает прямую х = 0 (см. рис. 7).

На рис. 6 показана область, где определяется течение газа в первом случае, когда сжимающий поршень х = хр{£) порождает тепловую волну с заданным законом движения фронта х — БЬ.

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

Далее в параграфе предложена одна разностная схема для системы (1.7),

Рис. 6

х=£><-г ' х=а(1)

Рис. 7

позволяющая при заданном фронте тепловой волны рассчитать поля течений.

Во-первых, начальные данные поставлены на фронте х = а(£), то есть на прямой 2 = 0. Но поскольку на прямой г = 0 имеют место бесконечные градиенты газодинамических параметров, то от нее отступают на некоторую прямую г = го < 0 (|.г0| мал). Именно на этой прямой г = го и ставятся граничные условия для всех параметров газа, которые, в свою очередь, вычисляются с помощью конечных отрезков рядов.

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

Расчеты проводятся для тепловых волн с заданным фронтом а(Ь) — 1+ +М2 (А = 0,|,1,10).

Для задачи о моделировании тепловой волны, когда фронт движется с постоянной скоростью а'(£) = 1 (Л = 0), имеется точное решение - бегущая по холодному фону тепловая волна. Бегущая по холодному фону тепловая волна с а(£) = £ позволяет протестировать как используемую разностную схему, так и приближенное решение, построенное с помощью конечных отрезков рядов.

В случае ускоренного движения фронта: а(£) = £ + А£2, а'(£) = 1 + 2А£ (А = 1,10) - в главе производится сравнение приближенного решения, полученного с помощью конечных отрезков рядов, с численным решением, построенным с помощью разностной схемы.

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

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

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

1. Баутин С.П., Садов А.П. Бегущая волна при учете равновесного излучения // Прикладная механика и техническая физика.-2006.-Т. 47.-Щ.-С. 15-25.

2. Баутин С.П., Садов А.П. Одномерная тепловая волна в невязком газе // Вычислительные технологии. 2006- Т. 11.-К1 5.-С 11-20.

3. Баутин С.П., Садов А.П. Численное исследование бегущих волн в теплопроводном невязком газе // Проблемы прикладной математики: Сборник научных трудов.-Екатеринбург. УрГУПС,-КЦ1(124).-В 2-х т.: Т.1.-2005-2006. С. 22-43.

4■ ЬаОов А.11. Плоскосимметричная тепловая волна со специальным краевым температурным режимом // Проблемы прикладной математики и механики: Сборник научных трудов. - Екатеринбург: УрГУПС. -№58(Ц1 )/2м. - В 2-х т.: Т.1.- 2007. С. 204-221

5. Садов А.П Численно-аналитическое моделирование одномерных тепловых волн // Проблемы прикладной математики и механики: Сборник научных трудов. -Екатеринбург• УрГУПС. -М!58(Ц1)/2м.-В 2-х т.: Т.1.-2007.-С. 162-203.

6. Садов А.П. Цилиндрически- и сферически-симметричная тепловая волна в невязком газе. Деп. в ВИНИТИ.-№ 1533-В2006.-2006.-26 с.

7. Садов А.П. Тепловая волна в невязком газе // Проблемы теоретической и прикладной математики: Труды 37-й регион, молодежи, конференц. Екатеринбург: УрО РАН.-2006.- С. 241-247.

8. Садов А.П. Плоско симметричная тепловая волна со специальным краевым температурным режимом // Молодые ученые - транспорту-2007: Сборник научных трудов, посвященный 170-летию российских железных дорог,- Екатеринбург: УрГУПС.-2007.-С. 502-519.

9. Баутин С.П., Николаев Ю.В., Рощупкин A.B., Чернышов Ю.Ю., Ягупов С А. Математическое моделирование безударного сильного сжатия газа //IX Всероссийский съезд по теоретической и прикладной механике. Аннотации докладов. - Нижний Новгород -Т.2.-2006 -С. 25-26.

10. Баутин С.П, Чернышов Ю.Ю., Садов А П. Особенности течений теплопроводного невязкого газа Тезисы междунар. конф. "Лаврентьевс-кие чтения по математике, механике и физике". Новосибирск, 2005.-С. 24-25.

11 Баутин С.П., Чернышов Ю.Ю., Садов А.П. // VII Забабахинские научные чтения• Тезисы докладов междунар. конф.-Снежинск: РФЯЦ-ВНИИТФ, 2005.-С. 25-26.

12. Садов А.П. Тепловая волна при учете равновесного излучения // VIII Забабахинские научные чтения. Тезисы докладов междунар. конф,-Снежинск. -2007.-С. 229.

13. Садов А.П. Одномерная тепловая волна в невязком газе // Аналитические методы в газовой динамике "САМГАД-2006": Тезисы докладов 21-й всерос. конф -Санкт-Петербург-2006.- С. 73-74-

14■ Садов А.П. Цилиндрически- и сферически-симметричная тепловая волна в невязком газе // Устойчивость, управление и моделирование динамических систем: Материалы Международной научной конференции, посвящ. 75-летию со дня рожд. И.Я. Каца.- Екатеринбург: УрГУПС. - т4(137). - 2006. - С. 82. ■

САДОВ Алексей Павлович

ЧИСЛЕННО-АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОДНОМЕРНЫХ ТЕЧЕНИЙ ТЕПЛОПРОВОДНОГО НЕВЯЗКОГО ГАЗА

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

Типография УрГУПС МПС РФ,620034,г. Екатеринбург,ул. Колмогорова,66

Подписано в печать 23.10.2008 г. Формат 69 х 90 1/16. Объем 2 п.л. Заказ 292 Тираж 100 экз.

Оглавление автор диссертации — кандидата физико-математических наук Садов, Алексей Павлович

ВВЕДЕНИЕ.

1. АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОДНОМЕРНЫХ ТЕПЛОВЫХ ВОЛН

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

1.2 Цилиндрически- и сферически-симметричная тепловая волна с заданным фронтом

1.3 Плоскосимметричная тепловая волна со специальным краевым температурным режимом.

2. БЕГУЩИЕ ВОЛНЫ В ТЕПЛОПРОВОДНОМ НЕВЯЗКОМ ГАЗЕ.

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

2.2 Численно-аналитическое исследование бегущих волн в теплопроводном невязком газе.

3. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПЛОСКО

СИММЕТРИЧНЫХ ТЕПЛОВЫХ ВОЛН.

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

3.2 Приближенное построение закона движения сжимающего поршня и краевого режима.

3.3 Численное построение полей течений газа.

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

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

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

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

Первое направление математического моделирования сильного сжатия газа состоит в использовании точных решений систем уравнений газовой динамики для политропного газа. Для этого случая наиболее исследованными являются одномерные неустановившиеся течения. Для описания некоторых плоскосимметричных течений газа применяется центрированная волна Римана - классическое решение системы уравнений газовой динамики, обладающее особенностью jf ( El—El j 5 x* — const, \ t - U J где t, x\ - независимые переменные, вектор if задает параметры течения газа. При t < t* центрированная волна Римана описывает течение, возникающее при безударном сильном сжатии плоского слоя в идеальном газе, что описано в работе К.П. Станюковича [85]. Е.И. Забабахиным и И.Е. Забабахиным [53] указано, что впервые центрированная волна Римана применена Гюгонио и Ре'леем для описания сжатия плоского слоя газа до сколь угодно большой плотности. В случае цилиндрически- и сферически-симметричных течений особенностью, аналогичной особенности в центрированной волне Римана, обладают автомодельные решения Л.И. Седова [84], с помощью которых осуществляется описание безударного сильного сжатия первоначально однородного и покоящегося в цилиндре или шаре идеального политропного газа. Интерпретации этих решений для задач о безударном сильном сжатии газа посвящены работы Я.М. Каждана [55-57], И.Е. Забабахина и В.А. Симоненко [52], А.Н. Крайко [59-66].

Точные автомодельные решения задачи об истечении газа в вакуум для двумерных и трехмерных течений были применены А.Ф. Сидоровым [83] для описания безударного сильного сжатия до бесконечной плотности газа, который в начальный момент покоится внутри призмы или многогранника при согласованных значениях показателя политропы газа 7 и двугранных углов.

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

В монографии С.П. Баутина [10] предложен и в последующем развит

30] единый подход к математическому моделированию сильного сжатия газа. Прошедшее с 1997 года десятилетие показало эффективность предложенного подхода. При этом подходе сначала ставятся начально-краевые задачи, описывающие процесс безударного сильного сжатия произвольного, локально-аналитического фонового течения на произвольной локально-аналитической поверхности. Для поставленных начально-краевых задач доказываются теоремы существования и единственности аналитических и кусочно-аналитических решений. Исследуются свойства решений, в том числе, устанавливаются асимптотические законы поведения газодинамических параметров при неограниченном росте плотности. С использованием строго полученных математических свойств исследуемых процессов далее численно и аналитически решаются конкретные задачи, моделирующие процессы сильного сжатия газа, в том числе, для специальных конфигураций сжимаемых объемов, используемых в физических экспериментах, а также при учете дополнительных физических факторов (излучение, теплопроводность) и реальных уравнений состояния. Многие полученные ранее точные решения, в том числе, центрированная волна Римана вкладываются в решения, найденные с использованием методики С.П. Баути-на [10,30], как частные случаи, при которых обрываются соответствующие ряды и получаются конечные формулы. Этот подход получил дальнейшее развитие при математическом моделировании течений газа в работах С.П. Баутина и его учеников: C.JI. Дерябина [31-36], A.JI. Казакова, Ю.В. Николаева [39,69-71], A.B. Рощупкина [39,74], Ю.Ю. Черны-шова [39,43-45,90], С.А. Ягупова [39,48,91-92] и автора [39-42,46-47,7582]. В том числе, С.П. Баутиным совместно с C.JI. Дерябиным [31-36] при t > t* рассмотрено одно- и двумерное истечение газа в вакуум в различных ситуациях, в том числе и в случае произвольного уравнения состояния р = р1f(p, S) с аналитической функцией /. Обобщение центрированной волны на случай уравнения состояния в виде р = р(р, е) принадлежит С.П. Баутину [10]. С.А. Ягуповым [91] построены течения, обобщающие свойства центрированной волны Римана на случай изэнтропических течений водорода с реальной изэнтропой, полученных в физических экспериментах. Ю.В. Николаевым [69-71] при решении различных задач численно построены обобщения центрированной волны в цилиндрически- и сферически-симметричных случаях, описывающие безударное сжатие газа до плотностей, в десятки-сотни тысяч раз превышающих исходную плотность однородного покоящегося газа. A.B. Рощупкипым [74] рассчитаны двумерные течения с соответствующей особенностью до плотностей, в 68 раз превышающих исходную плотность однородного покоящегося газа. Естественно, такие степени сжатия недостаточны для моделирования требуемого в экспериментах сильного сжатия газа. Тем не менее рассчитанные степени двумерного сжатия больше, чем в бесконечно сильной ударной волне.

Влияние лучистой теплопроводности на процесс кумуляции плоского слоя газа численно было исследовано М.Г. Анучиным [1] с использованием комплекса программ "Тигр".

Математическое моделирование тепловых волн проводилось многими авторами. Для нелинейного уравнения теплопроводности Г.И. Баренблатом [2-4] впервые были приведены примеры тепловых волн, распространяющихся по холодной среде с конечной скоростью, то есть сжимаемость среды и ее движение не учитывалось. Применение методологии характеристических рядов к построению тепловых волн было предложено А.Ф. Сидоровым [83], в частности, в одномерном плоскосимметричном случае доказано существование и единственность аналитического решения при специальном виде фронта тепловой волны. С.П. Баутиным [26] доказано существование и единственность аналитической тепловой волны у нелинейного уравнения теплопроводности как при задании аналитического фронта тепловой волны, так и при задании краевого температурного режима. С.П. Баутиным и А.А. Елисеевым [37] доказаны существование и единственность аналитического решения у нелинейного уравнения теплопроводности в многомерном случае при заданном произвольном аналитическом краевом режиме с отличной от нуля в момент Ь — 0 производной по времени. С.С. Титовым [89] при рациональных значениях константы сг, являющейся показателем степени в нелинейном коэффициенте теплопроводности, в плоско-, цилиндрически- и сферически-симметричных случаях в виде специальных рядов построены решения задачи о тепловой волне с заданным фронтом, обладающие на фронте тепловой волны конкретными особенностями. С.С. Титовым [86-87] приведены примеры точных решений нелинейного уравнения теплопроводности в виде специальных многочленов. С.С. Титовым [88] также получены решения нелинейного уравнения теплопроводности для симметричного случая в виде сходящихся логарифмических рядов. С.П. Кур-дюмовым [67] в автомодельном виде построено плоскосимметричное течение теплопроводного невязкого газа, в котором тепловая волна распространяется по холодной среде с конечной скоростью по фону, имеющему строго положительную температуру. Более общее течение теплопроводного невязкого газа с аналогичным свойством построено С.П. Баутиным [30]. Также С.П. Баутиным построены тепловые волны, распространяющиеся в теплопроводном вязком газе при задании фронта их движения в случае, когда показатель степени температуры в коэффициенте теплопроводности равен единице [38]. Аналог центрированной волны Римана в теплопроводном невязком газе описан в одномерном случае С.П. Баутиным и Ю.Ю. Чернышевым [43-45], в двумерном - Ю.Ю. Чернышевым [90]. Для идеального теплопроводного газа в случае холодного фона (Го = 0) Е.И. Забабахи-ным и В.А. Симоненко [52] показано, что в бесконечно сильной волне реализуется течение с изотермическим скачком. В случае Го > 0 в работах Я.Б. Зельдовича и Ю.П. Райзера [54], Л.Д. Ландау и Е.М. Лифшица [68], Б.Л. Рождественского и H.H. Яненко [73] установлено, что при малых значениях скорости D движения бегущей волны имеет место непрерывный переход, а при больших D - изотермический скачок. Бегущие волны в теплопроводном газе с уравнениями состояния, учитывающими равновесное излучение, были ранее рассмотрены для холодного фона с Го = 0 в работах В.А. Белоконя [49] и в уже упомянутой книге [54]. Также найдены два режима: при малых D изотермический скачок, при больших D - непрерывный переход.

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

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

Для моделирования плоских тепловых волн используется один частный класс решений рассматриваемой системы уравнений с частными производными: бегущие как по холодному, так и по нехолодному фону волны, то есть решения, зависящие от одной независимой переменной z = х — D ■ t, у которых константа D задает скорость движения бегущей волны и, не нарушая общности, полагается D > 0, то есть бегущая волна движется слева направо. Исследование этого класса течений осуществляется с помощью построения численного решения соответствующей системы обыкновенных дифференциальных уравнений (СОДУ).

Цели работы:

1. Доказательство существования плоско-, цилиндрически- и сферически-симметричных тепловых волн при учете равновесного излучения с использованием математической модели теплопроводного невязкого газа в окрестности фронта тепловой волны. Доказательство проведено в двух случаях: 1.1) когда задан аналитический закон движения фронта тепловой волны; 1.2) когда задан аналитический краевой температурный режим на некоторой аналитической поверхности.

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

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

• IX Всероссийский съезд по теоретической и прикладной механике, Нижний Новгород, 2006 г.;

• Международная конференция "Лаврентьевские чтения по математике, механике и физике", посвященная 105-летию со дня рождения академика М.А. Лаврентьева, Новосибирск, ИгИЛ, СО РАН, 2005 г.;

• Международная конференция "VII Забабахинские научные чтения", Снежинск, РФЯЦ-ВНИИТФ, 2005 г.;

• Международная конференция "VIII Забабахинские научные чтения", Снежинск, РФЯЦ-ВНИИТФ, 2007 г.;

• 21-ая Всероссийская конференция "Аналитические методы в газовой динамике САМГАД-2006", ИгИЛ СО РАН, ИПМ РАН, Санкт-Петербург, 2006 г.;

• 37-ая Региональная молодежная конференция "Проблемы теоретической и прикладной математики", Екатеринбург, УрО РАН, 2006 г.;

• Международная научная конференция "Устойчивость, управление и моделирование динамических систем", посвященная 75-летию со дня рождения И.Я. Каца, Екатеринбург, УрГУПС, 2006 г.;

Международная научно-практическая конференция "Снежинск и наука - 2006. Трансфер технологии, ииновации, современные проблемы атомной отрасли", Снежинск, СГФТА, 2006 г.;

• Молодежная научно-практическая конференция "Молодые ученые -транспорту", Екатеринбург, УрГУПС, 2006 г.;

• Молодежная научно-практическая конференция "Молодые ученые -транспорту", Екатеринбург, УрГУПС, 2007 г.

Публикации. По теме диссертации опубликовано 14 работ [39-42, 7582]. Из них - восемь в виде статей [40-42,75,78,80-82] (в том числе, две - в журналах из списка ВАК [41-42]), шесть в виде тезисов [39,46-47,7677,79].

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

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

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

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

3. В виде бесконечных сходящихся рядов, передающих особенность на фронте, построены плоскосимметричные тепловые волны в теплопровод ном невязком газе в случае специально подобранного краевого темпера турного режима.4. В виде бегущих волн численными расчетами решений системы обык новенных дифференциальных уравнений смоделирована тепловая волна, идущая по нехолодному фону. Расчетами установлены значения констант DQ, D\, D2 и наличие трех режимов распространения фронта бегущей вол ны в зависимости от величины D - скорости движения бегущей волны: при А) < D < D\ - непрерывный переход; при D\ < D < D2 - изотермический скачок; при D<i < D снова непрерывный переход.5. С помощью всех построенных решений приближенно восстанавлива ются: закон движения сжимающего поршня, порождающего заданную теп ловую волну; краевой температурный режим в заданной точке физического пространства; параметры течения газа во всей рассматриваемой области.Все полученные в диссертации результаты являются новыми. В част ности, с использованием аналитического подхода математически смодели рованы ранее не исследовавшиеся плоско-, цилиндрически- и сферически симметричные течения газа с учетом равновесного излучения.Впервые строго математически доказано существование тепловых волн, распространяющихся по холодному теплопроводному невязкому газу, когда коэффициент теплопроводности к = щТ

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

Библиография Садов, Алексей Павлович, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Анучин М.Г. Влияние теплопроводности на неограниченное безударное сжатие плоского газового слоя // ПМТФ. - 1998. - Т. 39, № 4. - 25-32.

2. Баренблат Г.И. О некоторых неустановившихся движениях жидкости и газа в пористой среде // Прикладная математика и механика. - 1952. -Т. 16. - Вып. 1. - 67-78.

3. Баренблат Г.И. Об автомодельных движениях сжимаемой жидкости в пористой среде // Прикладная математика и механика. - 1952. - Т. 16. -Вып 6. - 679-699.

4. Баренблат Г.И., Зельдович Я.Б. Промежуточные асимптотики в математической физике // Успехи мат. наук. - 1971. - Т. 26, № 2.- 115-129.

5. Баутин С П . Исследование области сходимости специальных рядов, решающих некоторые задачи газовой динамики // Числ. методы механики сплошной среды. - 1978. - Т. 9, № 4. - 5-17.

6. Баутин С П . Охлопывание одномерной полости // Прикл. математика и механика. - 1982. - Т. 46. - Вып. 1. - 50-59.

7. Баутин С П . Одномерное истечение газа в вакуум // Числ. методы механики сплошной среды. - 1983. - Т. 14, № 4. - 3-20.

8. Баутин С П . Представление решений системы Навье-Стокса с помощью характеристических рядов // Динамика сплошной среды. - 1987. -Вып. 83. - 11-31.

9. Баутин С П . Математическая теория безударного сильного сжатия идеального газа. - Новосибирск: Наука, 1997. - 160 с.

10. Баутин С П . О возможности изэнтропического перехода от однородного покоя в другое однородное покоящееся состояние идеального газа // Докл. РАН. - 1998. - Т. 362, № 5. - 621-624.

11. Баутин С П . Безударное сверхсжатие идеального газа // Вычислит, технологии. - 1998. - Т. 3, № 6. - 3-8.

12. Баутин С П . Возникновение градиентной катастрофы при фокусировке на ось или в центр симметрии волн сжатия и разрежения; УрГАПС, Екатеринбург, 1999. - 30 с. - Деп. в ВИНИТИ 18.01.1999. № 109-В99.

13. Баутин С П . Асимптотические законы безударного сильного сжатия квазиодномерных слоев газа // Прикл. математика и механика. - 1999. -Т. 63. - Вып. 3. - 415-423.

14. Баутин С П . О задаче получения наперед заданных распределений параметров газа // Прикл. математика и механика. - 1999. - Т. 63. - Вып. 6. - 938-946.

15. Баутин С П . Некоторые применения математической теории безударного сильного сжатия идеального газа // Тр. междунар. конф. "V Забаба-хинские научные чтения". - Ч. 1. - Снежинск: РФЯЦ-ВНИИТФ, 1999. -С. 14-19.

16. Баутин С П . О существовании решений задачи А.Н. Крайко // При- кл. механика и техн. физика. - 2000. - Т. 41, № 3. - 48-55.

17. Баутин С П . Слабые разрывы в течениях теплопроводного невязкого газа // Докл. РАН. - 2001. - Т. 377, № 4. - 481-484.

18. Баутин С П . Характеристические поверхности в течениях газа // Прикладная математика и механика. - 2001. - Т. 65. - Вып. 5. - 862-875.

19. Баутин С П . Математическое исследование безударного сжатия газа // Успехи механики. - 2002. - Т.1, №2. - С 3-36.

20. Баутин С П . Существование аналитического решения у краевой задачи для нелинейного уравнения теплопроводности; УрГУПС, Екатеринбург, 2002. - 28 с. - Деп. в ВИНИТИ 24.06.2002, № 1175-В2002.

21. Баутин С П . Тепловая волна нелинейного уравнения теплопроводности. - Екатеринбург: УрГУПС, 2002. - 80 с.

22. Баутин С П . Существование тепловой волны, порожденной краевым режимом // Вычисл. технологии. - 2002. - Т. 7, совм. вып., ч. 2. - 91-96.

23. Баутин С П . Существование аналитической тепловой волны определяемой заданным краевым режимом // Сибирский журнал индустриальной математики. - 2003. - Т. VI, № 1(13). - С 3-11.

24. Баутин С П . Тепловая волна, порожденная заданным краевым режимом // Доклады РАН. - 2003. - Т. 391, № 3. - 327-330.

25. Баутин С П . Аналитическая тепловая волна. М.: Физматлит, 2003.

26. Баутин С П . Математическое описание безударного сильного сжатия одномерных объемов политропного газа // Вычисл. технологии. - 2004. -Т. 9, № 1. - 26-33.

27. Баутин С П . Бегущая волна в теплопроводном невязком газе // Докл. АН. - 2006. - Т. 407, № 3. - 335-340.

28. Баутин С П . Различные виды условий Гюгонио в теплопроводном невязком газе // Проблемы прикл. математики. - Т. 1: Сб. науч. тр. - Екатеринбург: УрГУПС, 2006. - № 41(124). - 4-21.

29. Баутин С П . Математическое моделирование сильного сжатия газа.- Новосибирск: Наука, 2007. - 312 с.

30. Баутин СП., Дерябин Л. Двумерное истечение в вакуум // Точные и приближенные методы исследования задач механики сплошной среды: Тр. ИММ УНЦ АН СССР. - Свердловск: Изд-во УНЦ АН СССР, 1983. -С 3-15. из

31. Баутин СП., Дерябин Л. Истечение идеального газа в вакуум // Докл. АН СССР. - 1983. - Т. 273, № 4. - 817-820.

32. Баутин СП., Дерябин Л. Задача об истечении в вакуум нормального газа // Динамика сплошной среды. - 1993. - Вып. 107. - 26-38.

33. Баутин СП., Дерябин С Л . Математическое моделирование истечения идеального газа в вакуум. - Новосибирск: Наука, 2005. - 390 с.

34. Баутин СП., Дерябин Л. Аналитическое моделирование истечения идеального газа в вакуум // Успехи механики. - 2006. - Т. 4, № 1. - 77-120.

35. Баутин СП., Елисеев А.А. Многомерная аналитическая тепловая волна, определяемая краевым режимом // Дифференц. уравнения. - 2006. - Т. 42, № 8. - С 1052-1062.

36. Баутин СП., Кузьминых М.Ю. Тепловая волна в невязком газе // Актуальные проблемы прикладной математики. Екатеринбург: ИММ УрО РАН, 2003. - 13-14.

37. Баутин СП., Садов А.П. Численное исследование бегущих волн в теплопроводном невязком газе // Проблемы прикладной математики: Сборник научных трудов. - Екатеринбург: УрГУПС - № 41(124). - В 2-х т.: Т. 1. - 2005-2006. - 22-43.

38. Баутин СП., Садов А.П. Бегущая волна при учете равновесного излучения // Прикладная механика и техническая физика. - 2006. - Т. 47, № 4. - С 15-25.

39. Баутин СП., Садов А.П. Одномерная тепловая волна в невязком газе // Вычислительные технологии. - 2006- Т. 11, № 5. - 11-20.

40. Баутин СП., Чернышов Ю.Ю. Аналог центрированной волны Рима- на в теплопроводном невязком газе // Докл. РАН. - 2001. - Т. 380, № 1. -С. 43-46.

41. Баутин СП., Чернышов Ю.Ю. Одно течение теплопроводного газа, аналогичное центрированной волне Римана // Прикладная математика и механика. - 2002. - Т. 66, Вып. 1. - С 87-94.

42. Баутин СП., Чернышов Ю.Ю. Аналог центрированной волны Римана в теплопроводном невязком газе // Вычисл. технологии. - 2002. - Т. 7, № 4. - 3-8.

43. Баутин СП., Чернышов Ю.Ю., Садов А.П. Особенности течений теплопроводного невязкого газа. Тезисы междунар. конф. "Лаврентьевские чтения по математике, механике и физике". - Новосибирск, 2005. - 24-25.

44. Баутин СП., Чернышов Ю.Ю., Садов А.П. // VII Забабахинские научные чтения: Тезисы докладов междунар. конф. - Снежинск: РФЯЦ-ВНИИТФ, 2005. - С 25-26.

45. Баутин СП., Ягупов А. Математическое исследование изэнтропи- ческого сильного сжатия водорода // Прикл. математика и механика. -2003. - Т. 67, вып. 1. - 42-48. - .

46. Белоконь В.А. Влияние радиации на амплитуду изотермического скачка // Докл. АН СССР. - 1972. - Т. 202, № 6. 1296-1299.

47. Васин В.В., Сидоров А.Ф. О некоторых методах приближенного решения дифференциальных и интегральных уравнений // Известия вузов. Математика. - 1983. - № 7(254). - С 13-27.

48. Долголева Г.В., Забродин А.В. Кумуляция энергии в слоистых системах и реализация безударного сжатия. М.: Физматлит, 2004.

49. Забабахин Е.И., Симоненко В.А. Сходящаяся ударная волна в теплопроводном газе // Прикладная математика и механика. - 1965. - Т. 29, вып. 2. - 334-336.

50. Забабахин Е.И., Забабахин И.Е. Явления неограниченной кумуляции. М.: Наука, 1988.

51. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М.: Наука, 1966.

52. Каждан Я.М. Сферический разлет газа к центру. - М., 1969. - (Препр. / Ин-т прикл. математики РАН; № 2).

53. Каждан Я.М. К вопросу об адиабатическом сжатии газа под действием сферического поршня // Журн. прикл. механики и техн. физики. -1977. - № 1 . - С . 23-30.

54. Каждан Я.М. Адиабатическое сжатие газа под действием цилиндрического поршня. - М., 1980. - (Препр. / Ин-т прикл. математики РАН; № 56).

55. Ковалевская С В . Научные труды. К теории дифференциальных уравнений в частных производных. - М.: Изд-во АН СССР, 1948. - 368 с.

56. Крайко А.Н. О возникновении ударных волн при неравновесных течениях // Прикл. математика и механика. - 1967. - Т. 31, вып. 5. - 794-803.

57. Крайко А.Н. Вариационные задачи газовой динамики. - М.: Наука, 1979.

58. Крайко А.Н. Вариационная задача об одномерном изэнтропическом сжатии идеального газа // Прикл. математика и механика. - 1993. - Т. 57, вып. 5, - 35-51.

59. Крайко А.Н. О свободном нестационарном расширении идеального газа // Изв. РАН. Механика жидкости и газа. - 1993. № 4. - 155-163.

60. Крайко А.Н. Асимптотические закономерности нестационарного расширения идеального газа в пустоту // Прикл. математика и механика. -1994. - Т. 58, вып. 4. - 70-80.

61. Крайко А.Н. О неограниченной кумуляции при одномерном нестационарном сжатии идеального газа // Прикл. математика и механика. -1996. - Т. 60, вып. 6. - 1000-1007.

62. Крайко А.Н. Структура течений разрежения и сжатия в окрестности точки отражения «граничной» характеристики // Тр. Мат. Ин-та им. В.А. Стеклова. - 1998. - Т. 223. - 187-195.

63. Крайко А.Н. Структура сферического течения разрежения в окрестности точки отражения «граничной» характеристики // Прикл. математика и механика. - 1999. - Т. 63, вып. 6. - 972-979.

64. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Том VI. Гидродинамика. - М.: Физматлит, 2001. - 736 с.

65. Николаев Ю.В. О численном решении задачи безударного сильного сжатия одномерных слоев газа // Вычисл. техно логии. - 2001. - Т. 6, № 2. - 104-109.

66. Николаев Ю.В. О численном решении задачи Крайко в точке г = 0 и в ее окрестности; УрГУПС, Екатеринбург, 2004. - 29 с. - Деп в ВИНИТИ 15.11.2004, № 1768-В2004.

67. Николаев Ю.В. Численное решение задачи Крайко // Вычисл. технологии. - 2005. - Т. 10, № 1. - 90-102.

68. Овсянников Л.В. Лекции по основам газовой динамики. М. - Ижевск: Институт компьютерных исследований, 2003.

69. Рождественский Б.Л., Яненко Н.Н. Системы квазилинейных уравнений и их приложения к газовой динамике. М.: Наука, 1968.

70. Рощупкин А.В. Исследование некоторых характеристических задач Коши, возникающих при решении неодномерных задач безударного сильного сжатия газа // Вычисл. технологии. - 2002. - Т. 7, совм. вып., ч. 4. -С. 96-103.

71. Садов А.П. Тепловая волна в невязком газе // Проблемы теоретической и прикладной математики: Труды 37-й регион, молодежи, конференц. Екатеринбург: УрО РАН. - 2006. - 241-247. Иб

72. Садов А.П. Одномерная тепловая волна в невязком газе // Аналитические методы в газовой динамике "САМГАД-2006": Труды 21-й всерос. конф. - Санкт-Петербург. - 2006. - 73-74.

73. Садов А.П. Цилиндрически- и сферически-симметричная тепловая волна в невязком газе; УрГУПС, Екатеринбург, 2006. - 26 с. - Деп. в ВИНИТИ 11.12.2006, № 1533-В2006.

74. Садов А.П. Тепловая волна при учете равновесного излучения // VIII Забабахинские научные чтения: Тезисы докладов междунар. конф. -Снежинск. - 2007. - 229.

75. Садов А.П. Плоскосимметричная тепловая волна со специальным краевым температурным режимом // Проблемы прикладной математики и механики: Сборн. науч. трудов. - Екатеринбург: УрГУПС. - № 58(141)/2м. - В 2-х т.: Т. 1. - 2007.- 204-224.

76. Садов А.П. Численно-аналитическое моделирование одномерных тепловых волн // Проблемы прикладной математики и механики: Сборник научных трудов. - Екатеринбург: УрГУПС. - № 58(141)/2м. - В 2-х т.: Т. 1. - 2007. - 162-203.

77. Садов А.П. Плоскосимметричная тепловая волна со специальным краевым температурным режимом // Молодые ученые - транспорту-2007: Сборник научных трудов, посвященный 170-летию российских железных дорог. - Екатеринбург: УрГУПС. - 2007. - 502-519.

78. Сидоров А.Ф. Избранные труды: Математика. Механика. - М. ФИЗМАТЛИТ, 2001. - 576 с.

79. Седов Л.И. Методы подобия и размерности в механике. - М.: Наука, 1981. - 448 с.

80. Станюкович К.П. Неустановившиеся движения сплошной среды.- М.: Гос. изд-во техн.-теор. лит-ры, 1955. - 804 с.

81. Титов С. Решение двумерного уравнения фильтрации газа в виде многочлена по пространственным переменным // Динамика многофазных сред. - Новосибирск: ИТПМ СО АН СССР, 1981. - 291-293.

82. Титов С О Точные решения многомерного симметричного уравнения фильтрации в виде обобщенного многочлена // Динамика многофазных сред. - Новосибирск: ИТПМ СО АН СССР, 1983. - 287-290.

83. Титов С О Представление решений нелинейного осесимметрического уравнения фильтрации газа в виде логарифмического ряда // Динамика сплошной среды. - 1984. - Вып. 68. - 132-144.

84. Титов С. О движении фронта нелинейной диффузии // Прикл. механика и техн. физика. - 1996. Т. 37, № 4. - 113-118.

85. Чернышов Ю.Ю. Двумерная центрированная волна в теплопроводном невязком газе // Вычисл. технологии. - 2002. - Т. 7, № 4. - 107-115.

86. Ягупов А. Безударное сильное сжатие газа с реальными уравнениями состояния // Вопр. атомной науки и техники. Сер. Математическое моделирование физических процессов. - 2001. - Вып. 2. - 49-58.

87. Ягупов А. О возможности изэнтропического сжатия водорода при больших значениях давления // Вычисл. технологии. - 2002. - Т. 7, совм. вып., ч. 4. - 342-346.

Похожие работы

Информатика, вычислительная техника и управление
05.13.00