CC BY
38
Вычислительные технологии
Том 10, № 1, 2005
ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ А. Н. КРАЙКО
Ю. В. Николаев Уpальский государственный университет путей сообщения,
Екатеринбург, Россия e-mail: YNikolaev@math.usart.ru
The problem, suggested by A.N. Kraiko in the series of articles [1-5], is investigated. In these articles the process of a strong compression of a one-dimensional gas layer "from quiescent to the quiescent state" without formation of shock waves is considered. Other authors investigate similar problems on a compression of a quiescent gas into a moving uniform flow. In this case the energy of external compressing power is spent both on the compression and acceleration of gas. The way of compression, suggested by A.N. Kraiko, is optimal from the viewpoint of energy losses: quiescent gas is compressed to the quiescent one, therefore energy is spent on compression only.
1. Постановка задачи
В работе численно исследуется задача, предложенная А.Н. Крайко в серии статей [1-5].
В пространстве физических переменных t, r будем рассматривать одномерные изоэнтро-
пические течения, возникающие в политропном газе с уравнением состояния p = pY/7.
Рассматриваемые течения являются решениями системы уравнений
пространственная переменная; 7 — показатель адиабаты; р — давление; V — параметр геометрии: V = 0 в случае плоских течений, V =1 — цилиндрических, V = 2 — сферических.
Пусть в начальный момент времени Ь = Ь0 газ является покоящейся (и = 0) однородной средой с плотностью р =1. Это состояние газа назовем состоянием 1. Состояние газа в момент времени Ь = Ь* > Ь0 такое: газ однороден, его плотность р = р* > 1, при этом газ покоится. Это состояние газа назовем состоянием 2.
Тогда задача Крайко формулируется следующим образом: требуется найти течения газа, возникающие при безударном переходе одномерного газового слоя из состояния 1 в состояние 2. Другими словами, требуется безударным способом сжать покоящийся однородный газовый слой с плотностью р = 1 в покоящийся однородный газовый слой с плотностью р = р* > 1. В таком виде задача Крайко сформулирована С.П. Баутиным [6].
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2005.
(1.1)
для скорости газа и и скорости звука с = р(т 1)/2. Здесь Ь — время; г
vs
S X > 0
В [7] для плоскосимметричного случая приведен пример составного течения газа, в котором с помощью двух центрированных волн и постоянного потока осуществляется переход из однородного состояния покоя газа (состояние 1) в другое состояние покоя (состояние 2). Указанное течение вызвано движением двух непроницаемых поршней. В [7. С. 216] доказано, что решений такого типа не существует в том случае, когда один из поршней остается в покое. В [8] для случаев V =1, 2 доказано существование ненулевой массы газа, которую можно безударно сжать, используя начальную стадию конфигурации из [7]. А в [9, 10] приведено численное решение такой задачи при V = 0.
В [11] рассмотрена задача об оптимальном движении непроницаемого поршня, совершающего максимальную работу при заданных ограничениях на перемещение и время перемещения. Для плоскосимметричного случая приведены результаты численных расчетов и построена аналогия между этим оптимальным движением поршня и известным движением в двумерном сверхзвуковом сопле максимальной тяги [12].
В [1-5] построены и используются, в том числе для описания неограниченной кумуляции газа, составные плоско-, цилиндрически- и сферически-симметричные нестационарные течения газа, описывающие изоэнтропический переход идеального газа из состояния 1 в состояние 2 с большим или меньшим значением плотности р (волна сжатия или разрежения соответственно). Эти течения отличаются от конфигурации, представленной в [7]. В [2, 4] рассмотрены течения, когда сжимающий поршень приходит в точку, из которой распространяется центрированная волна сжатия. Тогда в состоянии 2 плотность постоянна и скорость газа равна нулю. В [1-5] для раскрытия особенностей в решении приведены приближенные формулы, а построение течения "в целом" сведено к численному решению системы уравнений газовой динамики методом характеристик.
В [6, 13, 14] доказано существование решения задачи Крайко в случае, когда газ сжимается двумя поршнями, один из которых покоится в точке с г = г0 > 0, а второй сжимает газ снаружи, т.е. координата второго поршня все время больше г0. Там же показано, что поскольку г = го (координата неподвижного поршня) строго положительна, то факт существования решения задачи Крайко имеет место и при сжатии газового слоя изнутри.
Однако доказанные теоремы имеют локальный характер и не указывают конкретных размеров области существования решения, что не позволяет установить количественные значения массы газа, для которых возможно существование решения задачи Крайко. Цель данной работы — решив задачу Крайко численно, указать значения массы газа, сжимаемого безударным способом.
2. Начально-краевые задачи
Система уравнений газовой динамики инвариантна относительно сдвига по времени а также инвариантна при одновременной смене знака у скорости газа и у времени. Это свойство позволяет вместо задачи о сжатии газа исследовать задачу о разрежении газа. Поэтому без нарушения общности далее рассматривается случай волны разрежения. Полученная волна разрежения при рассмотрении в обратном направлении изменения времени является искомой волной сжатия газа.
Рассмотрим плоский, цилиндрический или сферический слой однородного газа (р = р*), который в момент £ = 0 покоится между двумя непроницаемыми стенками, расположенными в точках 01(г = г0) и А(г = г*) (рис. 1). Точка 01 находится левее точки А: г0 < г*. При V =1 или V = 2 точка О1 — строго правее оси или центра симметрии соответственно: г0 > 0.
t
D
С
F
О
о]
А
г
Рис. 1.
Для системы (1.1) сформулируем три начально-краевые задачи, последовательное решение которых и даст решение задачи Крайко.
Задача 1 (о выдвижении поршня). Пусть при I > 0 из точки А выдвигается поршень (см. рис. 1). Его траектория (линия А¥) задается уравнением г = гр(£). В зависимости от способа выдвижения поршня возможны три случая: первый случай — поршень выдвигается с нулевой скоростью; второй — поршень выдвигается резко, но не быстрее скорости истечения газа в вакуум (именно эта ситуация приведена на рис. 1), третий — поршень выдвигается резко, со скоростью, большей скорости истечения газа в вакуум. В [6] показано существование единственного решения задачи 1 для указанных трех случаев. В данной работе рассмотрен случай, когда поршень выдвигается резко, но скорость выдвижения меньше, чем скорость истечения газа в вакуум:
жения поршня решение такой задачи существует в некоторой окрестности точки A, оно единственно в классе кусочно-аналитических функций и состоит из трех течений, разделенных звуковыми характеристиками. Первое из течений, расположенное в области AOiB, соответствует однородному состоянию покоя. Второе при v = 0 является центрированной волной Римана, а при v = 1,2 — ее обобщением [8] и имеет конкретную особенность в начальный момент времени (область ABC). Третье течение фактически представляет собой решение задачи о плавном выдвижении поршня из заданного течения (в рассматриваемом случае — из обобщения центрированной волны Римана) и не имеет особенностей в некоторой окрестности точки A (область правее линии AF).
Итак, в физическом пространстве решение, соответствующее волне разрежения, определяется однозначно, через звуковую характеристику AB непрерывно склеивается с решением, соответствующим состоянию исходного однородного покоя, и в некоторой окрестности точки B(tB,r0) задается аналитическими функциями, за исключением точки A. В точке A особенность течения известна (см. [2, 6]) и подобна особенности в центрированной волне.
Задача 2 (об отражении волны разрежения от жесткой стенки). В решении задачи 1 скорость газа при r = r0, t > tB строго больше нуля, поэтому для него условие
непротекания на стенке 0102 выполняться не будет. Следовательно, в результате взаимодействия течения из задачи 1 и стенки 0102 возникнет новое течение, отделенное от первого звуковой характеристикой ВС из семейства С+-характеристик (см. рис. 1). Новое течение расположено в области СВ02. Функция г = задающая характеристику ВС, и значения параметров газа на ней
с|г = <^) = с1(Ь)> и|г = <^) = и1(Ь) (2-1)
однозначно определяются решением задачи 1 и в некоторой окрестности точки Ь = Ьв являются аналитическими функциями. При этом = и1(Ь) + с1(Ь), р(Ьв) = г0. Новое течение из области СВ02 на характеристике ВС должно удовлетворять условиям (2.1) непрерывного примыкания к решению задачи 1. На стенке 0102 для этого течения должно выполняться условие непротекания
и|г = го = 0. (2.2)
Задача (1.1), (2.1), (2.2) является задачей 2. Справедлива следующая теорема (доказательство приведено в [6]):
Теорема 2. Задача 2 в некоторой окрестности точки (Ь = Ьв, г = г0) имеет единственное аналитическое решение.
Задача 3. На стенке 0102 выбирается какая-либо точка Е с координатами (Ье, го), Ье > Ьв (см. рис. 1), лежащая в области определения решения задачи 2. В выбранной точке Е однозначно определяется значение скорости звука газа с0 из решения задачи 2: с(Ье, г0) = С0, а следовательно, и значение плотности газа р0 в этой точке Е: р0 = с0/(т 1). Значение р = р0 является вторым произвольным элементом, который в совокупности со значением г* — г0 (первый произвольный элемент) однозначно определяет решение задачи Крайко в случае волны разрежения. Таким образом, расчетами будет искаться ответ на вопрос: какой ширины слой газа можно сжать в задаче Крайко и во сколько раз.
Выбор точки Е также единственным образом определяет аналитические функции г = и г = <^+(Ь) = г0 + с0(Ь — Ье). Эти функции задают траектории звуковых характеристик. Линия ЕД — характеристика г = семейства С--характеристик решения задачи 2, определенного в области ЕВС. Прямая ЕЕ есть характеристика семейства С+, отделяющая однородный покоящийся газ с плотностью р0 и расположенный в области 02ЕЕ (см. рис. 1) от течения в области ДЕЕ, являющегося решением задачи 3. Далее единственным образом определяются аналитические функции
с|г = р-(Ь) = с— (Ь)'с— (ье) = с0, и|г = = и-(ь),и-(ье) = 0; (2-3)
с|г = = с+ (Ь) = С0, и|г = = и+ (Ь) = 0, (2.4)
которые задают соответственно значения параметров состояния газа на характеристике ЕД уравнения (2.3) и на характеристике ЕЕ уравнения (2.4). При этом ф±(Ь) = и±±с±(Ь), ) = г0.
Итак, задача 3 (1.1), (2.3), (2.4) сводится к решению задачи Гурса (задача с данными на двух пересекающихся характеристиках ЕД и ЕЕ) в области ДЕЕ (см. рис. 1). Для доказательства существования решения задачи 3 последняя сводится к задаче о распаде слабого разрыва (см. [6]), для которой в классе кусочно-аналитических функций доказаны соответствующие теоремы.
Для расчета следующих слоев берутся п значений скорости звука с, распределенных равномерно от с* = р!7 1)/2 до 0, по формуле
. С*
С? — С* 7 ,
п
для 7 = 1,..п. Для каждого значения с? рассчитывается значение скорости газа и из формулы (3.1), верной для всех трех случаев V = 0,1, 2. Таким образом, получается п пар значений (с?, и?).
Для каждой пары (с? , и?) вычисляется угол наклона характеристики С+, выходящей из точки (¿*, г*), по формуле
? = и? + с?, 7 — 1,...,п-
Таким образом раскрывается особенность течения в точке (£ = ¿*,г = г*) и рассчитывается пучок характеристик, выходящих из точки А, необходимых для дальнейших расчетов.
Построение первого слоя происходит следующим образом (рис. 2, а): из точки (¿1,г0) на характеристике С+ выпускается характеристика С-. Затем по формулам метода характеристик находится точка пересечения С+-характеристики с С--характеристикой и значения скорости и плотности газа в найденной точке. Из точки (¿0,г0) на характеристике С+ выпускается характеристика С—. Характеристика С+ продолжается до ее пересечения с характеристикой С—. В найденной точке пересечения этих двух характеристик по формулам метода характеристик рассчитываются значения параметров газа. И так далее: продолжаем С+-характеристику до пересечения с С^-характеристикой, выходящей из т-й точки на характеристике С+, и в этой точке определяем параметры газа. Завершает построение первого слоя нахождение точки пересечения С+-характеристики с вертикальной прямой г = 1. Учитывая, что скорость газа и на неподвижной стенке 0102 равна нулю, то легко получить
Ь = 2и - Я = -Я. (3.2)
Таким образом, для расчета параметров газа в последней точке слоя достаточно найти изменение инварианта Я вдоль С+-характеристики.
Найденные точки пересечения вместе со значениями инвариантов Римана в них образуют слой с номером 1.
Слой с номером 2 строится по точкам слоя 1 аналогично. Для его построения из точки А выпускается характеристика С+ (с угловым коэффициентом, вычисленным по известным значениям (с2,и2)) до пересечения с характеристикой С-. Затем характеристика С+ продолжается до пересечения с характеристикой С2- и т. д. Заканчивается построение слоя с номером 2 нахождением точки пересечения характеристики С+ и вертикальной прямой г =1.
Подобным образом рассчитываются следующие слои характеристической сетки. Значение плотности в последней точке слоя уменьшается с увеличением номера слоя. Таким образом, построение каждого нового слоя заканчивается достижением меньшего значения плотности на стенке 0102. Слои строятся до тех пор, пока не будет достигнуто заданное значение плотности р0 = 1.
Пусть точка (¿я, го = 1) — последняя рассчитанная точка на линии г =1 ив ней достигнуто требуемое значение плотности р0, а N — номер последнего построенного слоя. Из точки (¿я, г0 = 1) (точка Е на рис. 2, б) выпускаем С—-характеристику (индекс "н"
Рис. 2.
означает, что характеристика отделяет рассчитываемое течение от несжатого покоящегося газа), заданную следующим уравнением:
Г = Го + (¿е - 1).
На этом первый этап расчета характеристической сетки завершен.
На втором этапе построение характеристической сетки сопровождается одновременным расчетом траектории движения сжимающего поршня. Построение этой траектории происходит в обратном направлении изменения времени, при 1 < 1*. Из точки (1 = 1*, г = г *) выпускаем траекторию поршня (см. рис. 2, б), определяемую из следующей задачи Коши:
^ = ,г00) + Ь(1$ ,г°0 ) 2 '
ги* =Г*.
Затем на С— -характеристике с шагом по времени Д12 расставляются точки, координаты которых вычисляются по следующим формулам:
1к = 1е - кД12, Гк = Го - Со кД12, к = 0,1,...
Из найденных точек выпускаются С+-характеристики. При этом значения плотности газа и скорости на С— -характеристике известны, поэтому значения инвариантов Римана на ней тоже известны.
Используя точки слоя с номером N, рассчитывается новый слой — характеристика С++1 (см. рис. 2, б). Новый слой (характеристика семейства С +) строится до пересечения с С--характеристикой, выходящей из точки А, значения параметров газа на указанной С--характеристике в точке А равны значениям параметров газа на характеристике С+.
Затем из точки А выпускается траектория сжимающего поршня (скорость поршня равна скорости газа на характеристике С+ в точке А), рассчитываются точки пересечения траектории поршня с линиями характеристической сетки, пока не будет рассчитана точка пересечения траектории поршня с характеристикой С++1. Затем строится новый слой — характеристика С++2. Новый слой строится до ячейки, в которой лежит траектория движения поршня, траектория поршня продолжается до пересечения с характеристикой С++2. Аналогично строится следующий слой — характеристика С++3 и т. д.
Скорость сжимающего поршня рассчитывается следующим образом. В узлах характеристической сетки скорость известна, поэтому в найденной точке пересечения характеристики с линией траектории поршня скорость газа и можно найти с помощью линейной интерполяции значений скорости. Построенная таким образом траектория движения поршня должна достичь характеристики фонового течения — С—-характеристики. Тогда на этом закончатся построение траектории сжимающего поршня и решение всей задачи Крайко "в целом", поскольку определится точка (£ = ¿8, г = г8). Эта точка лежит на С— -характеристике, и из нее в момент времени £ = ¿3 стартует непроницаемый поршень, безударно сжимающий газ в задаче Крайко.
Заметим, что найденное значение ¿а в общем случае отлично от нуля, но это легко поправить: достаточно перенести систему координат ¿Ог параллельно оси О£.
4. Результаты расчетов и выводы
Описанный выше алгоритм был реализован в виде программы для персонального компьютера класса 1ВМ РС. Написанная программа хранит в памяти два массива: первый массив — один полностью рассчитанный слой характеристической сетки, другой массив — второй слой, который рассчитывается по данным слоя из первого массива. После заполнения второго массива данные из него копируются в первый массив, а затем во втором массиве рассчитывается следующий слой.
Для измерения числовых значений в физической модели поставленной задачи были введены следующие единицы измерения. За единицу расстояния берется расстояние от оси симметрии (в случае цилиндрической симметрии) или от центра симметрии (в случае сферической симметрии) до стенки О102. Плотность несжатого газа считается равной 1, тогда значение р* показывает, во сколько раз сжат газ. За единицу измерения скорости берется скорость звука в несжатом газе с0 = 1. После этого масштаб измерения времени определяется однозначно.
Проверка точности расчетов. Точность выполняемых программой расчетов проверяется по значениям двух величин: относительной погрешности масс ¿т и по величине работы Ар, совершенной сжимающим поршнем.
Относительная погрешность массы газа: ¿т = -0' .100% , где т* — масса сжатого
т*
газа, т0 — масса несжатого газа. Так как сжатие производится непроницаемыми поршнями, то массы газа до и после сжатия должны совпадать. Но из-за ошибок численного расчета (округление чисел компьютером, приближенные формулы метода) наблюдается различие масс газа до и после сжатия. Поэтому, чем меньше значение ¿т, тем точнее выполнены расчеты траектории сжимающего поршня и всего течения в целом.
Вторая величина, также характеризующая точность вычислений, — это значение работы, выполненной сжимающим поршнем. Работа по сжатию газа, при фиксированных
т*, 7, р*, для случаев, когда течения газа плоско, цилиндрически или сферически симметричные, должна быть одинаковой, так как работа тратится только на сжатие газа — в начальный и конечный момент времени рассматриваемые объемы заполнены неподвижным газом. Таким образом, количество работы по сжатию газа не зависит от вида симметрии, и можно сравнивать точность расчетов по степени совпадения работы по сжатию газа (для фиксированных т*,7,р*) при различных видах симметрии V.
Работа рассчитывается в соответствии с общей формулой
А = J р^у.
VI
Также точность расчетов проверяется по степени совпадения работы, выполненной поршнем Ap, и изменения внутренней энергии газа — Е3 (так как вся работа сжимающего поршня тратится на сжатие газа, т. е. на изменение внутренней энергии), здесь — внутренняя энергия газа в финальный момент времени, а Е3 — внутренняя энергия газа в начальный момент времени. Изменение внутренней энергии газа рассчитывается по формуле
т(р*-1 - 1)
Ef - Ез = -1-—.
7 (7 - 1)
На точность результатов влияют значения параметров счета Д11, Д12, п. Эти параметры подбирались для каждого варианта расчетов отдельно, опытным путем, так, чтобы относительная погрешность масс была менее 1 % и работы сжимающего поршня Ар для случаев разных симметрий также отличались не более чем на 1 %.
Данные расчетов траектории сжимающего поршня. Траектория сжимающего поршня получена расчетами в табличном виде: набор (£, г) — координат точек пространства физических переменных, в этих точках также сосчитаны плотность и скорость газа.
Полученные данные о траектории движения сжимающего поршня приведены в таблицах 1,2. Все данные получены численными расчетами.
В табл. 1 указаны значения параметров газа, для которых проведены расчеты, и значения параметров счета Д11, Д12, п, при которых были сделаны расчеты. Также указаны значения величин, характеризующих точность выполненных расчетов: относительной погрешности масс работы Ар сжимающего поршня по сжатию газа и изменения внутренней энергии газа Ef - Е3.
В табл. 2 для различных значений параметров расчетов приведены данные о модуле скорости сжимающего поршня и плотности газа в некоторых точках траектории поршня: модуль скорости Uf и плотность газа рf в финальной точке траектории поршня (1f ,Гf); максимальное по модулю значение скорости поршня и, достигаемое в точке (1ад, гм). Также указано значение пространственной координаты г3 точки, из которой в момент времени 13 = 0 стартует сжимающий поршень. Во всех расчетах в момент старта модуль скорости поршня и3 равен нулю, плотность газа на поверхности поршня р3 равна единице.
На рис. 3 для трех различных вариантов приведены графики изменения скорости поршня, плотности газа на поршне и плотности газа на стенке 0102 от момента старта поршня 1 = 0 до момента финиша 1 = tf. Первый случай (см. рис. 3, а — в) соответствует варианту № 1; второй случай (см. рис. 3, а, г, д) — варианту № 2; третий случай (см. рис. 3, а, е, ж) — варианту № 5.
В варианте 1 сжимающий поршень монотонно увеличивает скорость до некоторого значения, а затем скорость поршня постоянна до получения финального распределения
Таблица!
Номер т* р* 7 V Д^2 п 5т Ар Ef — Es
1 1000 10000 1.4 0 0.00100 0.010 1000 0.11 69382 69305
2 1000 10000 1.4 1 0.00010 0.001 1000 0.73 69821 69305
3 1000 10000 1.4 2 0.00001 0.001 1000 0.30 69302 69305
4 10000 10000 1.4 0 0.01000 0.010 1000 0.04 693189 693048
5 10000 10000 1.4 1 0.00010 0.010 1000 0.55 696129 693048
6 10000 10000 1.4 2 0.00010 0.001 1000 0.15 689564 693048
7 1000 100000 1.4 0 0.01000 0.010 1000 0.34 177435 176786
8 1000 100000 1.4 1 0.01000 0.001 1000 0.28 176927 176786
9 1000 100000 1.4 2 0.01000 0.005 1000 0.19 171697 176786
10 10000 100000 1.4 0 0.10000 0.100 10000 0.14 1770384 1767857
11 10000 100000 1.4 1 0.01000 0.010 1000 0.66 1764950 1767857
12 10000 100000 1.4 2 0.01000 0.010 1000 0.61 1726041 1767857
13 1000 10000 5/3 0 0.01000 0.0100 1000 0.27 417700 416843
14 1000 10000 5/3 1 0.01000 0.0010 1000 0.03 414943 416843
15 1000 10000 5/3 2 0.01000 0.0001 1000 0.15 417279 416843
16 1000 100000 5/3 0 0.01000 0.0100 1000 0.85 1949063 1938091
17 1000 100000 5/3 1 0.01000 0.0100 1000 0.65 1877408 1938091
18 1000 100000 5/3 2 0.01000 0.0001 1000 0.30 1943278 1938091
19 1000 10000 2.0 0 0.10000 0.1000 10000 0.52 5004977 4999500
20 1000 10000 2.0 1 0.01000 0.0001 10000 0.60 5038625 4999500
21 1000 10000 2.0 2 0.00100 0.0001 1000 0.01 4984690 4999500
Таблица2
Номер Тз Тп и Ь rf Uf Pf
1 1002.11 1062.36 3.2853 13.28 1062.53 1.1000 13.28 652
2 17.93 19.29 1.0158 11.29 19.29 1.0158 11.29 1090
3 6.22 6.20 1.0079 11.10 6.20 1.0079 11.10 1143
4 10005.15 10616.66 23.8533 13.27 10618.30 2.0000 13.27 652
5 56.58 61.11 2.3829 9.35 61.24 1.1482 8.83 1935
6 13.37 14.41 1.5807 8.47 14.47 1.0740 8.23 2204
7 1004.39 1041.15 1.2832 22.50 1041.16 1.0100 22.50 5033
8 17.89 18.71 1.0016 20.85 18.71 1.0016 20.85 6772
9 6.21 6.11 1.0008 20.75 6.11 1.0008 20.75 7013
10 10015.01 10391.43 3.8318 22.50 10391.55 1.1000 22.50 5033
11 56.61 60.55 1.0158 18.35 60.55 1.0158 18.35 10163
12 13.39 14.25 1.0079 17.90 14.25 1.0079 17.90 10906
13 1003.72 1026.61 1.1023 30.83 1026.62 1.1000 30.83 1432
14 17.87 18.16 1.0158 28.05 18.16 1.0158 28.05 1813
15 6.22 5.75 1.0079 27.86 5.75 1.0079 27.86 1842
16 1009.46 1021.95 1.0153 68.23 1021.96 1.0100 68.23 13325
17 17.81 17.85 1.0016 66.28 17.85 1.0016 66.28 14680
18 6.22 5.63 1.0008 66.14 5.63 1.0008 66.14 14569
19 1006.21 1011.56 1.6517 99.00 1011.57 1.1000 99.00 2550
20 17.92 17.40 1.0158 94.32 17.40 1.0158 94.32 2792
21 6.21 5.47 1.0079 94.20 5.47 1.0079 94.20 2798
(см. рис. 3, б); похожим образом изменяется плотность газа на поверхности поршня (см. рис. 3, в). В варианте 2 скорость сжимающего поршня монотонно увеличивается до момента получения финального распределения (см. рис. 3, г); похожим образом изменяется плотность газа на поверхности поршня (см. рис. 3, д). В варианте 3 скорость сжимающего поршня монотонно увеличивается до некоторого значения, а затем монотонно убывает — поршень притормаживает (см. рис. 3, е); плотность газа на поршне изменяется, как в варианте 2: плотность монотонно увеличивается от момента старта поршня до момента финиша (см. рис. 3, ж).
Во всех случаях плотность газа на неподвижной стенке 0102 сначала постоянная (отрезок 01В на рис. 1), затем монотонно растет (отрезок ЕВ на рис. 1) и в конце плотность принимает постоянное значение (отрезок 02Е на рис. 1). Характер изменения плотности газа на неподвижной стенке для трех рассмотренных случаев изображен на рис. 3, а.
По результатам расчетов можно сделать следующий вывод. В случае плоской симметрии (V = 0) модуль скорости сжимающего поршня будет изменяться, как показано на рис. 3, б, монотонно возрастая от нуля до некоторого значения, а затем до момента получения требуемого финального распределения поршень будет двигаться с постоянной скоростью. Расчетами установлено, что на последней стадии сжатия в случае плоской симметрии между центрированной волной Римана и поршнем находится однородный движущийся с постоянной скоростью газ.
Для случев цилиндрической (V = 1) и сферической (V = 2) симметрии можно подобрать такое значение массы газа т = т0 (для каждого вида симметрии свое значение т0), что если сжимать газ меньшей массы, то характер изменения модуля скорости сжимающего поршня будет, как показано на рис. 3, г, т. е. монотонно увеличиваться. Если сжимать газ большей массы, то характер изменения модуля скорости сжимающего поршня будет, как показано на рис. 3, е, т. е. на финальной стадии — не монотонный. Характер изменения модуля скорости сжимающего поршня исчерпывается указанными тремя вариантами.
Расчетами установлено, что в случае, когда сжимающий поршень притормажимает перед получением финального распределения, он никогда не останавливается до скорости
Рис. 3.
и = 0. То есть в момент времени 1 = tf поршень входит в однородный покоящийся сжатый газ с ненулевой скоростью, это получено во всех рассчитанных вариантах.
Также расчетами установлено, что плотность на сжимающем поршне растет монотонно от 1 до рf при увеличении 1, при этом рf всегда меньше р*. Таким образом плотность газа на поршне всегда меньше плотности сжатого газа (аналогичное утверждение справедливо для давления).
Параметры газа на неподвижной стенке. По условиям поставленной задачи скорость газа на неподвижной стенке и = 0. На рис. 3, а показано поведение плотности газа на стенке г0 = 1 в зависимости от времени для вариантов 1, 2, 5. По данным, полученным численными расчетами, можно утверждать, что характер распределения плотности на неподвижной стенке для других вариантов аналогичен: плотность на стенке увеличивается монотонно от 1 до р*, значение плотности максимально в области покоящегося однородного с плотностью р = р* газа (состояние 2). Вместе с плотностью монотонно увеличивается давление на стенку.
Обсуждение результатов
В основу метода расчета положен проверенный метод характеристик, позволяющий с достаточной степенью точности восстанавливать течения в газе. В частности, для проведенных расчетов относительная погрешность масс газа до сжатия и после сжатия получилась меньше 1 %.
Для раскрытия особенности в течении использовано свойство решения, полученное строго математически: особенность описывается сходящимся рядом, первые слагаемые которого известны. Расчеты ведутся в обратном направлении изменения времени, течения рассчитываются в следующей последовательности: сначала рассчитываются течения в точке с особенностью, затем методом характеристик делается расчет в области, содержащей траекторию сжимающего поршня, где особенностей не возникает. Предложенная методика расчета была использована в [9, 10], полученные в этих работах результаты подтверждают надежность используемой методики расчетов. Все сказанное указывает на достоверность полученных результатов.
Численно построены течения, возникающие в газе при переходе из состояния 1 в состояние 2. В качестве одного из искомых элементов исследуемой задачи рассчитана траектория сжимающего поршня. Анализ рассчитанных течений показал, что полученное численное решение описывает процесс безударного сжатия одномерного слоя газа из покоя в покой, т. е. в случае конфигурации Крайко.
Самое главное: удалось получить конкретное значение массы газа, которую можно сжать безударно до плотности р*, в 104-105 раз превышающей исходную. Тем самым дан ответ на вопрос о возможности безударного сжатия газа в случае конфигурации Крайко достаточно больших масс газа до больших плотностей.
Выражаю благодарность своему руководителю, профессору С.П. Баутину за ценные замечания по данной статье, а также за его научное руководство в целом.
Список литературы
[1] Крайко А.Н. О свободном нестационарном расширении идеального газа // Изв. РАН. Механика жидкости и газа. 1993. Т. 4. С. 155-163.
[2] КрАйко А.Н. Вариационная задача об одномерном изоэнтропическом сжатии идеального газа// Прикл. математика и механика. 1993. Т. 57, вып. 5. C. 35-51.
[3] Крайко А.Н. Асимптотические закономерности нестационарного расширения идеального газа в пустоту // Прикл. математика и механика. 1994. Т. 58, вып. 4. C. 70-80.
[4] КрАйко А.Н. О неограниченной кумуляции при одномерном нестационарном сжатии идеального газа // Там же. 1996. Т. 60, вып. 6. C. 1000-1007.
[5] КрАйко А.Н., ТилляЕвА Н.И. Автомодельное сжатие идеального газа плоским, цилиндрическим или сферическим поршнем // Теплофизика высоких температур. 1998. Т. 36, № 1. C. 120-128.
[6] Баутин С.П. О существовании решений задачи А.Н. Крайко // Прикл. механика и техническая физика. 2000. Т. 41, № 3. С. 48-55.
[7] Мизес Р. Математическая теория течений сжимаемой жидкости. М.: Изд-во иностр. лит., 1961.
[8] Баутин С.П. Математическая теория безударного сильного сжатия идеального газа. Новосибирск: Наука, 1997.
[9] Баутин С.П., Николаев Ю.В. Об одном методе расчета безударного сильного сжатия одномерных слоев газа // Вычисл. технологии. 2000. Т. 5, № 4. С. 3-12.
[10] Николаев Ю.В. О численном решении задачи безударного сильного сжатия одномерных слоев газа // Там же. 2001. Т. 6, № 2. С. 104-109.
[11] Рылов А.И. Вариационная задача одномерной нестационарной газовой динамики // Изв. АН СССР. Механика жидкости и газа. 1984. Т. 4. C. 171-175.
[12] КрАйко А.Н. Вариационные задачи газовой динамики. М.: Наука, 1979.
[13] Баутин С.П. О возможности изоэнтропического перехода от однородного покоя в другое однородное покоящееся состояние идеального газа // Докл. РАН. 1998. Т. 362, № 5. С. 621-624.
[14] Баутин С.П. Математическое исследование безударного сжатия газа // Успехи механики. 2002. Т. 1, № 2. С. 3-36.
[15] Жуков А.И. Применение метода характеристик к численному решению одномерных задач газовой динамики // М.: Изд-во АН СССР, 1960. (Тр. Матем. ин-та им. В.А. Стеклова; Т. 58).
[16] Рождественский Б.Л., Яненко Н.Н. Системы квазилинейных уравнений и их приложения к газовой динамике. М.: Наука, 1968.
Поступила в редакцию 25 ноября 2003 г., в переработанном виде — 20 января 2004 г.
