УДК 532.5:533.6.011.5
ГОМЭНТРОПИЧЕСКАЯ МОДЕЛЬ ОТРАЖЕНИЯ СФЕРИЧЕСКОЙ УДАРНОЙ ВОЛНЫ ОТ ЦЕНТРА СХОЖДЕНИЯ
И.А. Чернов
Саратовский государственный университет, кафедра вычислительного эксперимента в механике E-mail: chernov-ia@yandex.ru
Обсуждается частный случай движения имплозивной ударной волны по покоящемуся газу с нулевым давлением, но с переменной плотностью. Плотность описывается степенной зависимостью от расстояния до точки фокусировки ударной волны. Предлагается такой выбор показателя степени в этой зависимости, чтобы энтропия во всей области течения после прохождения ударной волны была постоянной (гомэнтропичность). При этом получается качественно другое по сравнению с классическим случаем Гудерлея - Ландау - Станюковича поведение температуры.
Ключевые слова: одномерные течения, автомодельные течения, сходящаяся ударная волна, гомэнтропическая модель.
Homentropic Model of Spherical Shock Wave Reflection from the Center of Convergence
I.A. Chernov
Saratov State University,
Chair of Mechanics Computational Experiment
E-mail: chernov-ia@yandex.ru
An implosive shock wave on a based gas the particular case of motion with zero pressure, but with variable density is discussed. The density is described by degree relation to distance up to a point of focusing of a shock wave. Such selection of an exponent in this relation that the entropy in all area of flow after passage of a shock wave was a constant (homentropic case) is offered. Thus qualitatively different behaviour of temperature in comparison with classical case Guderley - Landau - Stanjukovich is obtained.
Key words: one-dimensional flows, self-similar flows, converging shock wave, homentropic model.
ВВЕДЕНИЕ
Решение задачи о сходящейся к центру ударной волне (УВ) было получено Гудерлеем [1] и независимо от него Л.Д. Ландау и К.П. Станюковичем [2, 3]. Эта задача имеет большое прикладное значение, поскольку описывает один из способов создания экстремальных условий в малой области пространства. Существование автомодельных решений степенного вида предопределено наличием двухпараметрической группы подобных преобразований соответствующей системы уравнений газовой динамики идеального газа [4]. Эти решения описываются нелинейной системой обыкновенных дифференциальных уравнений (ОДУ) первого порядка, в которую входят три независимых параметра: Y = Cp/Cv — показатель адиабаты Пуассона (p/p7 = const); ш — показатель степени в автомодельном представлении плотности p = p0twR (£) (p0 = const); n — показатель автомодельности, где £ = r/ (Ktn) — независимая автомодельная переменная (t — время, r — пространственная координата, K — масштабная константа). Выбор показателя 7 определяет свойства изучаемого газа. Нахождение n осуществляется решением специальной переопределенной краевой задачи названной системы ОДУ и находится на первом временном этапе (до момента фокусировки УВ). При выборе ш используют [1-3] модель сильной УВ (СУВ) и поскольку УВ распространяется по однородному газу, то полагают dp (t,£) /dt = 0, как до, так и после нее, тогда как p (t,£) — функция двух независимых переменных t, £. В результате энтропийная функция s = p/p7 для течения после ударной волны является функцией лагранжевой координаты, то есть сохраняется вдоль траектории частицы — это свойство изэнтропии, но не гомэнтропии (постоянство энтропии в двумерной области (r, t) — англояз. термин).
Ниже обсуждается другая теоретическая возможность: показатель ш выбирается из условия ds (t, £) /dt = 0 — это оказывается (вместе с условием адиабатичности) достаточным условием гомэнтропии (соответственно гомэнтропическая УВ (ГУВ)). Это означает, что энтропия меняется скачком на ГУВ одинаково для всех жидких частиц.
В работе Хантера [5], посвященной проблеме схлопывания пустой сферической полости (каверны) в воде, для изучения эффекта сжимаемости среды использовалась газодинамическая модель течения с показателем 7 = 7. Так как это течение начинается как гомэнтропическое, то автор использовал для определения давления уравнение состояния среды в конечной форме и рассмотрел систему двух дифференциальных законов: неразрывности и количества движения вместе с интегралом адиабаты.
Этот подход распространялся и на вторую стадию течения, когда возникала отраженная УВ (от центра после схлопывании каверны). Это пример течения, в котором энтропия постоянна всюду (не меняется на УВ), то есть выполняется глобальная гомэнропичность. Работа интересна, в частности, тем, что в ней описана техника построения ударной волны, которая возникает после схлопывания полости.
В дальнейшем это изучение для других значений 7, а также при использовании кусочно-постоянных значений энтропии до и после отраженной УВ (гомэнтропическая модель) продолжил Лазарус [6]. Задача о каверне изучалась параллельно с задачей о сходящейся к центру УВ с последующим отражением от центра, но здесь автор использовал классическую модель СУВ.
Ниже повторяется соответствующий результат Лазаруса для частного случая 7 = 7/5 и результаты сравниваются с новым расчетом при том же 7, но по гомэнтропической модели.
Если оставаться в рамках классического подхода, то представленный результат можно рассматривать как относящийся к частному случаю движения сильной УВ по покоящемуся газу с нулевым давлением, но с переменной по пространству начальной плотностью. Подобная трактовка используется в астрофизике [7].
Таким образом, цель работы — сравнение результатов расчета процесса схождения - отражения сферической УВ по двум моделям: классической [1-3] и по модели ГУВ для одного и того же показателя 7 = 7/5.
1. ОСНОВНЫЕ УРАВНЕНИЯ
Изучаются одномерные неустановившиеся адиабатические течения идеального совершенного газа для случая сферической симметрии течения. Оно характеризуется скоростью и = и(г,£), плотностью р = р(г, £) и квадратом локальной скорости звука с2 = с2(г, £), который с точностью до коэффициента совпадает с абсолютной температурой. Основные уравнения таковы:
d д (V — 1) ри „ d д ^
—р + р—и + ^-— = 0, 7-г: и + с2 + с2—1п(р)=0,
dt дг г dt дг дг
d / с2 \ d д д ( )
= 0, — = — + и-
dí VP7-V ' dt д£ дг'
Предполагая автомодельность, зададим искомые функции в форме (т = 1/п, п = 1/Сп):
Ктг1-т и 1(п) К2тг2-2т^1(п) I = К-тпгт,и = —-и1п), р = рог™Я1(п), с2 = К---2 ^1(п). (2)
т т2
Функции {и 1, Ш, Z1} называются автомодельными представителями (АП) (по г) функций {и, р, с2} соответственно. Они описывают поведение {и, р, с2} при фиксированном г в зависимости от £.
После подстановки заданного вида (2) в основную систему (1) получается три ОДУ для {и 1(п), Я1(п), Z 1(п)}.
Эта система приводится к нормальному виду (разрешенному относительно производных), затем она записывается в виде автономной системы четырех ОДУ (без знаменателей в правых частях с независимой переменной х — искусственным временем)
^ 1 (х) = Z 1(х) [(кти 1(х) + кцо) Z 1(х) + кюзи 1(х)3 + кю2и 1(х)2 + кти 1(х)] , -^Хи 1(х) = [—1 + и 1(х)п(х)] [(кз11 и 1(х) + кз1о) Z 1(х) + кзози 1(х)3 + кзо2и 1(х)2] ,
-х (3)
-х^1(х) = Я1(х) [(к211и 1(х) + к21о) Z1 (х) + к2о2и 1 (х)2 + к2о1и 1 (х)] , еdХn(x) = т7 [—1 — 2и 1(х)п(х) + и 1(х)2п(х)2 — Z 1(х)п(х)2] [—1 + и 1(х)п(х)] .
Коэффициенты к зависят от основных параметров задачи {т, т = о>/п, 7, V =1, 2,3 для плоского, цилиндрического, сферического случаев симметрии течения} и автомодельной независимой переменной п: к11о = —п (2т — 2 — т) (7 — 1), к111 = 2п27 (т — 1), к1оз = 7п2 (—7 + 3 — 2т + 7V — V), к1о2 = 7V (7 — 5 + 3т — 27V + 2v + 7т), к1о1 = —7 (7т — 7V — 2 + т + V), к21о = —п (2т — 2 — т),
&211 = —П2YW, ^201 + Y (-w + m - v), ^202 = Yn (1 - 2w + m - 2v), &203 = П2Y (-1 + w + v), ksio = 2m - 2 - w, k3ii = Yn (-v + m), ^302 = Y (m - 1) ^303 = -Yn (m - 1).
Эта система удобна для расчетов и в том частном случае, когда приходится проходить точку с П = 0, что соответствует t = 0 (момент фокусировки УВ), которая не является особой точкой для системы (3).
Чтобы иметь возможность сравнения результатов с результатами классического подхода, рассмотрим наряду с (2) традиционное представление основных величин [8]
r = Ktn £, u = nKtn-1£U (£), p = p0tw R(£), c2 = n2K2 t2n-2 £2 Z (£). (4)
Отметим различие: здесь U(£), Z(£) — это фазовые переменные, тогда как {U2 = £U(£), Z2 = £2Z(£), R2 = R(£)} — это АП (по t) искомых функций {u, p, c2}. Две функции U(£), Z(£) параметрически представляют решение ОДУ первого порядка для функции Z = Z(U) (фазовая плоскость {U,Z })
dUz (u )=nem (5)
num = -Z (5Uy + 2YnU2v - U3Ynv - 2Y2nU2v + U3Y2nv + y2nUv+ + Z(UЬш - UY2 + U2Y2 + 2YnZ(U) + 2YnU - 5YnU2 + 3U3Yn + Y2nU2--U3y2n - 2YnUZ(U) - YnUv + 2Z(U) - Z(U)ш - 2nZ(U) - 3U2y - 2y) ,
den = (-U3nY + U2y + nU2y + Z(U)YnUv - Uy - 2Z(U) + Z(U)ш + 2nZ(U)) (U - 1).
Качественный анализ этого ОДУ содержится в [8]. После того как решение Z = Z(U) найдено, следует осуществить две квадратуры для определения £ = £(U) и R = R(U).
Сравнивая два представления (2) и (4), легко получить соотношения:
£ = n-n, и (£ ) = nU 1(n), Z (£ ) = n2 Z 1(n), R(£) = ■
Если в течении есть автомодельная ударная волна nS = n1 = n2 = const, то при ее переходе выполняются три закона сохранения — массы, импульса и энергии. Считая параметры с индексом 1 известными, находят параметры с индексом 2 (U1 = U11, Z1 = Z11, R1 = R11, ...; U2 = U12, . ..) — для значений представителей сразу за скачком
U2 = YU12 - U12 + 3U 1n1 - YU 1n1 - 2n12 + 2Z1
= YU 1 + U1 - Yn1 - n1 '
= R1 (Yn12 + YU 12 - 2YU 1n1 + n12 + U12 - 2U 1n1) ()
R2 = Yn12 - n12 + 2Z1 + YU 12 - 2YU 1n1 - U12 + 2U 1n1' ( )
Z2 = (YU 1 + U1 - YnS - nS)2 (12Y2U12n12 - 12YU12n12 - Y2U12Z1 + 6Yn12Z1 - 2YZ12 - 2Yn14+
+6YU 12Z1 + 8YU 1n13 + 2U 1n1Z 1 - 12YU 1n1Z 1 + 2Y2U 1n1Z 1 - 8Y2U 13n1 - 8Y2U 1n13--Y2n12Z1 + 2Y2U14 - 2YU 14 - U12Z1 - n12Z 1 + 2Z12 + 2Y2n14 + 8YU 13n1) ■
Если УВ распространяется по неподвижному газу с {U10 = 0, Z10 = 0, R10 = p0}, то позади нее должно быть
U20 = 2L, Z20 = 2Yn12 (Y - 1), R20 = p0 (y + 1).
Y + 1 (Y + 1)2 Y - 1
p
Введем в рассмотрение энтропийную функцию s = —. Вдоль траектории частицы s имеет посто-
pY
янное значение в силу третьего уравнения из (1). Для гомэнтропического потока эта постоянная одна для всех траекторий. Выражая s в автомодельном случае через (t,£) с использованием формул (4), получим s = s0tw(1-Y)+2n-2S(£), где s0 = p0-Yn2K2, S = £2R(£)1-YZ(£).
Условием того, что 5 явно не зависит от времени, является равенство
и =
2 (n - 1)
Y - 1
Если предполагается гомэнтропичность течения в какой-то области, то это означает использование уравнения адиабаты Пуассона, которая определяет квадрат скорости звука через плотность. Вместо основных уравнений (1) можно использовать первые два из них вместе с интегралом адиабаты, при этом значение двух разных постоянных в интеграле вычисляется по (6).
2. ПОСТАНОВКА ЗАДАЧИ
Чтобы описать изучаемое явление, полезно рассмотреть плоскость (r,t) (рис. 1). Отрицательному времени соответствует фаза схождения УВ к центру. УВ1 изображает ее траекторию, ПХ — траекторию предельной характеристики, УВ2 — отраженной УВ. Каждая из них описывается обобщенной параболой г] = const. Полуплоскость г = 0 разделена на 4 подобласти, обозначенные цифрами 1-4. Подобласть 1 изображает покоящийся газ, траектория частицы — пунктирная вертикальная прямая, идущая вверх.
После пересечения ее с УВ1 возникает движение частиц к центру фокусировки r = 0 (ЦФ). Подобласть 2 является зоной влияния на УВ1: ПХ — последняя из центростремящихся характеристик, приходящих на УВ1. УВ2 изображает отраженную от ЦФ ударную волну, которая идет по газу, движущемуся навстречу. Рис. 1 может рассматриваться как гипотетическая картина явления в малой окрестности ЦФ, но в данном случае это также результат расчета по ГУВ модели.
Математическое моделирование включает:
1) обоснование факта использования автомодельных решений (2);
2) описание краевой задачи для соответствующей системы ОДУ как автомодельной задачи второго рода [8, 9], в которой показатель автомодельности определяется методом пристрелки из условия аналитического прохода известной интегральной кривой уравнения (5) через особую точку на фазовой плоскости, являющейся образом ПХ. На рис. 2 показаны две кривые (1 — по классической модели и 2 — по ГУВ) для y = 1.4. Особая точка ПХ лежит на звуковой линии (ЗЛ);
3) построение ударного перехода (см. рис. 2) от точки УВ2-к точке УВ2+ осуществляется таким образом, чтобы фазовая траектория могла пройти из точки УВ+ вверх к особой точке (U = (2 — 2n — и) / (nY), Z ^ которая является седлом и приход в которую обеспечивает продолжимость течения до точки r = 0 при t ^ для подобласти 4 рис. 1.
Дополнительные детали по методике расчетов есть в [7-9]. рис 2
ПХ УВ1
Рис. 1
-0,8
3. РЕЗУЛЬТАТЫ РАСЧЕТОВ
Первым этапом построения решения является определение показателя автомодельности при заданных 7, ш и нахождение АП для трех параметров течения в подобласти 2 рис. 1. Поскольку автомодельное решение определено с точностью до масштабной постоянной К в (2), это позволяет положить П = —1 на ПХ (такова предварительная нормировка). Далее можно построить решение системы (3) в виде степенных рядов по (п + 1) с учетом того, что нас интересует окрестность одной (с меньшим значением ^) из двух особых точек, которые соответствуют ПХ. Сделав небольшое отступление от п = —1 в сторону уменьшения п, можно получить начальные условия для интегрирования системы
(3) при выбранном 7 = 7/5, ш = 0 и некотором пробном значении п. Решая задачу Коши для системы в сторону уменьшения п, следят за поведением соответствующей фазовой траектории на плоскости (и, ^) с тем, чтобы попасть в известную неособую точку — образ УВ1+. Такое попадание возможно лишь при определенном подборе п и определенном размере интервала интегрирования по х. Так были получены первая и вторая строки в таблице и известное значение п = 0.7171745.
Модель Характерная точка п и 1 Я1 ^ 1
СУВ ш = 0 п = 0.7171745 7 = 7/5 УВ1+ -1.0 -5/6 6.0 7/36
ПХ- -0.845089 -0.773168 8.49174 0.1682177
ПХ+ -0.845065 -0.773182 8.49062 0.168228
УВ2- 2.68841 -0.287893 64.3086 0.144550
УВ2+ 2.68841 0.0794383 145.062 0.214518
16.7345 0.0112545 54.4395 0.185998
п ^ 50.03389 0.00373867 31.5279 0.172507
ГУВ ш = 0.7900640 п = 0.8419872 7 = 7/5 УВ1+ -1.0 -5/6 6.0 7/36
ПХ- -0.821996 -0.762431 6.95182 0.206241
ПХ+ -0.821976 -0.762424 6.95178 0.206242
УВ2- 1.44198 -0.158958 18.3459 0.304056
УВ2+ 1.44198 0.254178 35.5986 0.410791
16.0763 0.0195063 4.36336 0.177409
п ^ 50.4683 0.00619525 1.76270 0.123458
Аналогичное вычисление зависимости п = п(7) проведено для модели ГУВ, сравнение двух графиков (кривая 1 по модели СУВ, 2 — по ГУВ) дано на рис. 3.
Перейдем к описанию второго этапа построения решения. Используя найденные разложения для АП в окрестности п = —1 (их коэффициенты стали числами), находят начальные условия в задаче Коши для построения решения в подобласти 3 рис. 1. Интегрируют систему (3) в сторону увеличения п, так чтобы пройти точку с п = 0. Далее интегральная кри-
0,9-
0,8-
0,7-
0,6-
вая на плоскости (и, Z) выходит на ЗЛ, где Z = (и — 1)2
1
3 5
Рис. 3
7 У
точке, отличной от ПХ. Это означает предельную линию на (г, £)-плоскости, что недопустимо. Возникает задача определения местоположения отраженной УВ (см. п. 3 вышеприведенного математического моделирования) с тем, чтобы с помощью
скачка УВ2--> УВ2+ перепрыгнуть ЗЛ.
В таблице приведены характерные точки, полученные при интегрировании системы ОДУ для УВ1+, ПХ, УВ2-, УВ2+ и больших значений п- Расчеты были проведены для 2-х моделей (СУВ и ГУВ). Их использование позволяет легко воспроизвести АП и построить параметры соответствующих течений. В таблице была проведена вторая нормировка, так что значению УВ1 соответствует п = —1.
4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
На рис. 4-7 приведены АП (по г) для скорости частиц, плотности, квадрата локальной скорости звука и энтропийной функции параллельно для 2-х обсуждаемых моделей. Всюду кривые 1 относятся к классической СУВ-модели, 2 — к ГУВ. Напомним, что эти функции при изменении п от до описывают поведение соответствующего параметра при фиксированном г и изменении £ от до
в
На рис. 4 вертикальным пунктиром изображены УВ1 и УВ2. Область между УВ+ и ПХ- изображена утолщенной линией. Отметим, что сила (по и 1) отраженной УВ почти в 2 раза меньше, чем падающей. Отраженная УВ по модели ГУВ возникает раньше, чем по СУВ. Качественное поведение АП(и 1) одинаковое.
Рис. 4
На рис. 5 показан АП(Л1) (по г) для плотности. Сила (по Л1) отраженной УВ по СУВ-модели более чем в 3 раза превосходит по ГУВ-модели. Качественное поведение АП(Л1) одинаковое.
150 т N. 1 \
100- 1 |УВ2 1
50- ^___2
УВ1|-* 1 1 1 1 1
Рис. 5
На рис. 6 показан АП(^ 1) (по г) для квадрата локальной скорости звука (температуры с точно-чтью до постоянного множителя). Кривая 1 показывает, что температура в фиксированной точке г подскакивает на приходящей УВ1, затем начинает падать до минимума, затем слегка возрастает, подскакивает в отраженной УВ2 с тем, чтобы уменьшаться со временем. Таким образом, АП температуры — немонотонная функция в интервале между приходящей и отраженной УВ. Это означает наличие характерного размера. Кривая 2 имеет другое качественное поведение: температура продолжает расти в интервале между двумя УВ. Это можно объяснить тем, что газ адиабатически сжимается, так как частицы двигаются к ЦФ после сходящейся УВ. После отраженной УВ газ адиабатически расширяется, что приводит к его охлаждению.
На рис. 7 показаны АП(Б 1) (по г) для энтропийной функции. Кривая 1 после УВ1 монотонно и сильно убывает до появления УВ2, что представляется качественно правдоподобным, поскольку частицы, которые в начальный момент располагались ближе к ЦФ, проходят через более сильную УВ и получают большее значение энтропии позади нее. За отраженной УВ2 наблюдается монотонное возрастание энтропии, что не противоречит высказанному объяснению, так как скорость частиц в зоне 4 рис. 1 меняет знак. Кривая 2 демонстрирует кусочно постоянные значения Б1, что соответствует гомэнтропической модели.
Рис. 7
ЗАКЛЮЧЕНИЕ
Классическое решение [1-3] — это решение хорошо поставленной задачи с краевыми и начальными условиями. Такая задача является полностью автомодельной [10]: основные уравнения, начальные и граничные условия инвариантны относительно группы подобных преобразований. Решение такой задачи дает локальную асимптотику (фазу схождения УВ к ЦФ) и от него нельзя требовать правильного описания второй фазы (отражения УВ от ЦФ). Физически интересными могут оказаться частично автомодельные решения (уравнения инвариантны, но часть условий — нет). Иногда такие решения выступают в роли промежуточных асимптотик изучаемого явления, они, в частности, способны «забывать» начальные условия. Приведем простой пример: если в лужу бросить прямоугольный кирпич, то от него почти сразу пойдет волна в виде расходящегося круга, которая описывается частично автомодельным решением, «забывшим» о «прямоугольных» начальных условиях. Возможно, что представленное новое решение является промежуточной асимптотикой для схождения - отражения УВ по газу постоянной начальной плотности. Для окончательного вывода требуется сравнение с результатами анализа асимптотик в конечноразностных решениях соответствующих задач.
Автор благодарит В.С. Кожанова за помощь в оформлении статьи.
Библиографический список
1. Guderley, G. Starke kugelige und zylindrische Ver-dichtungsstoe in der Nahe des Kugelmittelpunktes bzw. der Zylinderachse / G. Guderley // Luftfahrtforschung. -1942. - B. 19, lfg. 9. - S. 302-312.
2. Ландау, Л.Д. Теоретическая физика. Т. VI. Гидродинамика / Л.Д. Ландау, Е.М. Лифшиц. - М.: Физмат-лит, 2003. - 736 с.
3. Станюкович, К.П. Неустановившиеся движения сплошной среды / К.П. Станюкович. - М.: Наука, 1971.
- 856 с.
4. Овсянников, Л.В. Групповой анализ дифференциальных уравнений / Л.В. Овсянников. - Новосибирск: Изд-во СО АН СССР, 1962. - 236 с.
5. Хантер, К. О захлопывании пустой полости в воде / К. Хантер // Механика: период. сб. пер. иностр. ст.
- 1961. - № 3 (67). - С. 77-100.
6. Lazarus, R.B. Self-Similar Solutions for Converging Shocks and Collapsing Cavities / R.B. Lazarus // SIAM J. Numer. Anal. - 1981. - V. 18, iss. 2. - P. 316-371.
7. Зельдович, Я.Б. Физика ударных волн и высокотемпературных гидродинамических явлений / Я.Б. Зельдович, Ю.П. Райзер. - М.: Наука, 1966. - 688 с.
8. Седов, Л.И. Методы подобия и размерностей в механике / Л.И. Седов. - М.: Наука, 1967. - 428 с.
9. Брушлинский, К.В. Об автомодельных решениях некоторых задач газовой динамики / К.В. Брушлинский, Я.М. Каждан // УМН. - 1963. - Т. 18, вып. 2 (110). - С. 3-23.
10. Баренблатт, Г.И. Подобие, автомодельность, промежуточная асимптотика / Г.И. Баренблатт. - Л.: Гид-рометеоиздат, 1978. - 207 с.