close

Вход

Забыли?

вход по аккаунту

?

Малые возмущения идеального газа волны инварианты и задача Коши.

код для вставкиСкачать
Наука и Образование. МГТУ им. Н.Э. Баумана.
Электрон. журн. 2014. № 12. С. 372{391.
DOI: 10.7463/1214.0744956
Представлена в редакцию:
Исправлена:
28.06.2014
02.12.2014
c МГТУ им. Н.Э. Баумана
УДК 517.956.32+533.72
Малые возмущения идеального газа:
волны, инварианты и задача Коши
Феоктистов В. В.1,* , Мякинник О. О.1
*
1
МГТУ им. Н.Э. Баумана, Москва, Россия
Модель волнового взаимодействия, разработанная авторами для нормальной системы линейных
однородных уравнений в частных производных 1-го порядка, которая задана числовыми матрицами порядка n при производных, рассмотрена в приложении к задаче о распространении малых
возмущений в идеальном газе при одномерном течении. Модель основана на понятии n-мерной
бегущей волны, коэффициент которой является матричным коэффициентом системы. Найдено
решение задачи с использованием традиционного подхода к решению гиперболической системы
и средствами модели волнового взаимодействия. Установлена связь между волнами, распространяющимися с характеристическими скоростями, и n-мерными бегущими волнами. Решение
задачи Коши представлено как разложение в ряд по n-мерной бегущей волне с коэффициентами,
равными коэффициентам разложения аналитической начальной вектор-функции в ряд Тейлора.
Ключевые слова: гиперболическая система; бегущие волны; инварианты Римана; n-мерные
бегущие волны; разложение по бегущим волнам; решение задачи Коши; идеальный газ; малые
возмущения
1. Введение
Привлечение средств функционального анализа для аналитического решения уравнений
математической физики является актуальной задачей. Предложенный авторами аппарат
аналитических функционалов специального вида имеет применение к исследованию систем
уравнений в частных производных 1-го порядка, которые заданы произвольными некоммутативными, вообще говоря, числовыми матрицами. В работе рассматривается приложение
этого аппарата к уравнениям газовой динамики.
Нестационарные уравнения газовой динамики, выражающие законы сохранения, являются уравнениями гиперболического типа. Свойство гиперболичности тесно связано с существованием решений специального вида, которые обладают свойством инвариантности.
Такие решения называют бегущими волнами [1, 2, 3, 4].
Напомним, что гиперболическая система n, n ? 1, однородных дифференциальных уравнений в частных производных 1-го порядка относительно независимых переменных (t, x),
Наука и Образование. МГТУ им. Н.Э. Баумана
387
x = xs , и однородное уравнение в частных производных порядка m, m ? 1, относительно
переменных (t, x), x = (x1 , . . . , xs ), s ? 1, имеют решения, которые сохраняют свои значения на характеристических прямых и гиперплоскостях соответственно [2, 3]. Случай n = 1,
s = 1 и m = 1, s = 1 соответствует простейшему уравнению в частных производных,
решением которого является элементарная волна [4, 5].
Для построения волны следует определить коэффициенты при независимых переменных
в выражении (?0 t+?·x), ?·x = ?1 x1 +. . .+?s xs , решив для этого алгебраическое уравнение,
построенное по определенному правилу из числовых коэффициентов дифференциального
уравнения или системы уравнений.
Хорошо известными примерами являются: система уравнений акустики (n = 2) и ее
решение в виде элементарных волн (xs ? at) и (xs + at), распространяющихся со скоростью
звука a вдоль координатной оси OXs [1, 2, 6]; волновое уравнение (m = 2) в пространстве s
переменных, которое порождает волну f (?0 t + ? · x) (здесь f | произвольная гладкая
функция), распространяющаяся в направлении вектора (?0 , ?) [1, 3].
В работах [7, 8, 9] для системы
s
~
~
X
?B
?B
E
=
Ak
,
?t
?xk
k=1
~ = (B1 (t, x), . . . Bn (t, x))т ,
B
x = (x1 , . . . , xs ),
(1)
где Ak , k = 1, s, | числовые квадратные матрицы порядка n, E | единичная матрица,
безотносительно к определяемому матрицами-коэффициентами типу рассмотрены решения
(Exk +Ak t), k = 1, s, и функции от них. По аналогии с элементарными волнами такие решения названы бегущими волнами с матричными коэффициентами или n-мерными бегущими
волнами. При s = 1 в случае гиперболичности системы (матрица As подобна диагональной жордановой матрице) введенная n-мерная бегущая волна имеет полную интерпретацию
через n элементарных волн, которые рассмотрены как скалярные [10].
Несмотря на существование численных методов решения системы (1) в приложении к
отдельным физическим процессам (например, [6]), актуальным является построение аналитического метода решения этой системы и связанной с ней задачи Коши при произвольных n, s и произвольных матричных коэффициентах Ak , k = 1, s. Для системы (1), которая
характерна для различных областей математической физики [11], в работах [7, 8, 9] авторами построена модель волнового взаимодействия. Элементами модели являются: n-мерные бегущие волны, содержащие коэффициенты системы Ak , k = 1, s, которые распространяются по координатным осям и находятся во взаимодействии; оператор волнового
взаимодействия, характеризующийся размерностью s; отвечающая системе операторная
экспонента.
В рамках модели общее решение системы (1) представлено как разложение в ряд по взаимодействующим бегущим волнам с неопределенными числовыми вектор-коэффициентами.
Найденное решение в явном виде содержит матрицы-коэффициенты Ak , k = 1, s, задаюНаука и Образование. МГТУ им. Н.Э. Баумана
388
щие систему. Для получения решения задачи Коши достаточно в общем решении заменить
неопределенные коэффициенты на векторные коэффициенты разложения аналитической начальной функции в ряд Тейлора [8].
В качестве физического процесса для интерпретации введенных в [7, 8, 9] понятий и выполненных построений выберем течение идеального газа, линейные и нелинейные волновые
процессы в котором изучены [1, 5, 12]. Выбранному процессу соответствует квазилинейная система дифференциальных уравнений в частных производных 1-го порядка [2, 5, 12],
которые представляют собой законы сохранения. В качестве замыкающих уравнений [12]
для уранений сохранения скорости и импульса рассмотрим уравнение сохранения энергии,
выразив ее через энтальпию, и уравнение состояния, связывающее энтальпию с плотностью и давлением газа. Таким образом, результирующая квазилинейная система записана
относительно n-мерного вектора, 3 ? n ? 5, единообразно содержащего в качестве своих
компонент компоненты скорости газа u(t, x), v(t, x), w(t, x), его плотности ?(t, x) и энтальпии I(t, x), где x = (x1 , . . . , xs ), s ? 3. Поскольку энтальпия посредством коэффициента
~ x)
теплоемкости связана с температурой газа, каждая из компонент искомого вектора B(t,
является величиной, которая может быть измерена в ходе эксперимента.
Чтобы из квазилинейного процесса выделить линейные волны, порождаемые системой (1), применим метод малых возмущений (например, [2]). Возникающую при этом
задачу о распространении малых возмущений, которая представляет самостоятельный интерес [2, 5, 13], рассмотрим как тестовую для модели волнового взаимодействия на примере
движения идеального газа при заданных n и s.
По структуре построенной модели общее решение системы
E
~
~
?B
?B
= Ak
,
?t
?xk
~ = (B1 (t, xk ), . . . Bn (t, xk ))т ,
B
k = 1, s,
(2)
представленное как разложение по волне (Exk + Ak t), является частным решением системы (1) и составной частью ее общего решения, представленного в виде разложения по s
взаимодействующим волнам [7, 9].
Пусть (2) отвечает распространению малых возмущений в идеальном газе при одномерном течении. Найдем решение этой системы с использованием традиционного подхода к
одномерной гиперболической системе [2, 6] и средствами модели волнового взаимодействия
без использования условия гиперболичности [7]. Покажем тождественность полученных решений и преимущества волнового подхода при исследовании системы (2).
Целью работы является рассмотрение предлагаемого авторами нового аналитического
метода решения нормальной системы линейных однородных уравнений в частных производных 1-го порядка на примере задачи о распространении малых возмущений в идеальном
газе при одномерном течении с последующим использованием полученного решения для
построения решения задачи при течении газа в пространстве.
Наука и Образование. МГТУ им. Н.Э. Баумана
389
2. Уравнения распространения малых возмущений в идеальном газе
Квазилинейные уравнения одномерного движения идеального газа при отсутствии массовых сил в декартовой системе координат имеют вид [2, 5, 13]:
?
? ?u
?
?u?
1 ? P?
+ u?
=?
,
?
? x?
?? ? x?
?
t?
?
?
? ??
?(?? u?)
? t?
? x?
? ? t? + ? x? = 0,
?
?
?
? ?I + u? ?I = 1 ? P? + u? ? P? ,
?
??
x? ? R1 ,
(3)
? x?
? t?
? t?x?) | скорость движения газа, его плотность, давление и
где u?(t?, x?), ??(t?, x?), P? (t?, x?), I(
энтальпия соответственно. Чтобы сделать систему трех дифференциальных уравнений относительно четырех неизвестнах функций определенной, дополним ее уравнением состояния [12]
??1 ? ?
cP
??I, I = cP T? , ? =
,
(4)
?
cV
которое соответствует баротропному движению идеального газа и связывает давление газа с
P? =
его плотностью и энтальпией (T? (t?, x?) | температура газа, cP | удельная теплоемкость газа
при постоянном давлении, cV | удельная теплоемкость газа при постоянном объеме). Вос?
пользуемся уравнением состояния (4) и запишем систему относительно функций (u,
? ??, I):
?u?
?u? ? ? 1 ? ? ?? ? ? 1 ?I?
?
?
= ???u?
I
?? ,
? x?
?
? x?
?
? x?
? t?
? ??
?u?
? ??
? u? ,
= ???
? x?
? x?
? t?
?I?
?u?
?I?
? u? .
= ?(? ? 1)I?
? x?
? x?
? t?
??
Выразим зависимые и независимые переменные через одноименные безразмерные переменные:
u? = U0 u,
?? = ?0 ?,
I? = I0 I,
t? = T0 t,
x? = L0 x,
где U0 , ?0 , I0 , T0 , L0 | размерные постоянные, которые являются характерными для задачи
скоростью, плотностью, энтальпией, временем и длиной соответственно.
Запишем систему одномерного движения идеального газа относительно безразмерных
функций (u, ?, I):
?
?u
1
??
1
?I
?u
?
?
?
St ?
= ??u
I
? ,
2
2
?
?t
?x ?M
?x ?M
?x
?
?
??
?u
??
? St ?t = ?? ?x ? u ?x ,
?
?I
?u
?I
?
? St
?u .
= ?(? ? 1) I
?t
?x
Наука и Образование. МГТУ им. Н.Э. Баумана
(5)
?x
390
В (5) St | число Струхала; St = L0 /U0 T0 ; M | число Маха; M = U0 /a, a | скорость звука
в газовой среде, при этом a2 = (? ? 1)I0 . Поскольку среди размерных величин, задающих
число Струхала, две являются независимыми, можно выбрать St = 1.
Рассмотрим задачу о распространении малых возмущений в идеальном газе при одномерном параллельном оси OX движении [2, 5]. Пусть в движущемся со скоростью U0
идеальном газе, который имеет плотность ?0 и энтальпию I0 , создаются малые возмущения
скорости u? (t, x), плотности ?? (t, x) и энтальпии I ? (t, x). Тогда движению возмущенного
газа отвечают функции
u(t, x) = 1 + u? (t, x),
?(t, x) = 1 + ?? (t, x),
I(t, x) = 1 + I ? (t, x).
(6)
Чтобы нелинейная система (5) имела решение вида (6), функции, задающие малые возмущения, должны удовлетворять системе линейных уравнений [2]:
?
?u
1 ??
1 ?I
? ?u
?
?
=?
,
2
? ?t
?x ?M ?x ?M 2 ?x
?
??
?u
??
(7)
? ?t = ? ?x ? ?x ,
?
?I
?u ?I
?
?
= ?(? ? 1)
?t
?x ?x
(для упрощения обозначений звездочки при искомых функциях опущены). Выполненная
процедура называется линеаризацией системы (5) около состояния (U0 , ?0 , I0 ) [2].
3. Приведение к независимым волновым уравнениям и нахождение инвариантов
Запишем полученную линейную систему (7) в матричном виде относительно искомой
~ x) = (u(t, x), ?(t, x), I(t, x))т , x ? R1 :
вектор-функции B(t,
?
~
~
?B
?B
E
=A
,
?t
?x
?
A=?
?
?1
?
1
?M 2
?
1 ?
?M 2 ?
?1
?1
0
?(? ? 1)
0
?1
(8)
?,
?
где E | единичная матрица порядка 3. Система (8) представляет собой систему линейных
однородных дифференциальных уравнений в частных производных 1-го порядка с постоянными коэффициентами, которая записана в нормальной форме. Свойства системы и ее
решений определяются алгебраическими свойствами матрицы A.
Матрица A является невырожденной, имеет три различные действительные собственные
значения
1?M
1+M
, ?3 = ?
M
M
и подобна диагональной жордановой матрице J c матрицей подобия S (J = S ?1 AS), где
?1 = ?1,
?
0
?
S=? 1
?
?2 =
1
1
?
?M
M
?
?,
?
?1 ?(? ? 1)M (? ? 1)M
Наука и Образование. МГТУ им. Н.Э. Баумана
?
S ?1
?
0 (? ? 1)/?
?1/? ?
?
= ? 1/2 ?1/2?M ?1/2?M ? .
?
?
1/2 1/2?M
1/2?M
391
Матрица S составлена из собственных векторов матрицы A, матрица S ?1 найдена как обратная к S. Согласно определению (например, [2]) система (8) (система (7)) является гиперболической (в узком смысле).
Выполним замену переменных
~ = S V~ ,
B
т
V~ (t, x) = (V1 (t, x), V2 (t, x), V3 (t, x)) ,
(9)
и рассмотрим систему уравнений
E
? V~
? V~
= S ?1 AS
.
?t
?x
(10)
Система (10) состоит из независимых уравнений
?Vi
?Vi
= ?i
,
?t
?x
i = 1, 3,
(11)
каждое из которых содержит одну неизвестную функцию и имеет решение
Vi = fi (x + ?i t),
i = 1, 3,
(12)
где fi | произвольные непрерывно-дифференцируемые функции. Прямые
x + ?i t = c i ,
ci = const,
i = 1, 3,
(13)
являются характеристиками системы и при пересечении образуют в плоскости (t, x) треугольник, называемый характеристическим [1]. Основание треугольника, отвечающее прямой (13) при i = 1, не зависит от параметров газа. Данный факт согласуется с законом
распространения малых возмущений в покоящемся газе [1] и при относительном движении
газа, которое рассматривается в разделе 4. Напомним, что для уравнений акустики в качестве
основания характеристического треугольника выбирается интервал на оси OX [1].
Функции fi , i = 1, 3, сохраняют значения на соответствующих параллельных прямых (13) и вместе со своими аргументами называются бегущими волнами [1, 3]. Уравнение в
частных производных 1-го порядка с двумя независимыми переменными (10) моделирует одномерный волновой процесс (волновой процесс в одномерной среде) и является самым простым из уравнений, к которым применяется определение <волновое>. Поэтому аргументы
функций fi называют также элементарными волнами [4] (в [5] | простыми волнами). В [2]
для волн (12) применен термин "бегущие волны малых возмущений", который, отражая
связь этих волн с исходной квазилинейной системой и методом ее линеаризации, указывает на важность задачи о распространении малых возмущений для исследования волновых
процессов. Бегущие волны малых возмущений распространяются с характеристическими
скоростями.
Выполнение преобразования подобия (9) по отношению к гиперболической системе (8)
привело к независимым уравнениям и появлению элементарных бегущих волн в качестве
решения.
Наука и Образование. МГТУ им. Н.Э. Баумана
392
Найденные выражения
~ ,
Vi ? S ?1 B
i
i = 1, 3,
записанные через физически значимые функции, сохраняют свои значения на характеристиках и, согласно определению [1, 2], представляют собой инварианты Римана системы (8):
(? ? 1)
1
?(t, x) ? I(t, x) = f1 (x ? t),
?
?
1
1
1
(1 ? M )
u(t, x) ?
?(t, x) ?
I(t, x) = f2 x +
t ,
2
2?M
2?M
M
1
1
1
(1 + M )
u(t, x) +
?(t, x) +
I(t, x) = f3 x ?
t .
2
2?M
2?M
M
(14)
(15)
(16)
Каждый из трех инвариантов (14) { (16) распространяется без искажения формы с характеристической скоростью ?i , i = 1, 3.
В связи с постановкой задачи (3){(4) и нахождением инвариантов (14){(16) отметим
работу [14], в которой доказано, что система одномерных уравнений газовой динамики
имеет три иварианта Римана, если в качестве замыкающих выбрано уравнение сохранения
энергии, выраженное через энтропию s, и уравнение состояния в виде
V = G(P ) + A(s),
1
V = ,
?
G0 (P ) = ?
1
g 2 (P )
,
c=
g(P )
,
?
(17)
где a | скорость звука; g и A | произвольные функции давления и энтропии.
Из (9) и (14) { (16) следует, что решение нормальной линейной системы уравнений в частных производных 1-го порядка (8), которая отвечает распространению малых возмущений в
идеальном газе при одномерном движении, имеет вид
?
(1 ? M )
(1 + M )
?
u(t,
x)
=
f
x
+
t
+
f
x
?
t
,
2
3
?
M
M
?
(1 ? M )
(1 + M )
x?
t ,
M
t + M f3
? ?(t, x) = f1 (x ? t) ? M f2 x +
M
?
(1 ? M )
(1 + M )
? I(t, x) = ?f1 (x ? t) ? (? ? 1)M f2 x +
t + (? ? 1)M f3 x ?
t .
M
(18)
M
Отметим, что скорость распространения малых возмущений u(t, x) выражена через волны,
отвечающие двум из трех характеристик (13). Данный факт согласуется с представлением
скорости распространения звука, найденной из уравнений акустики.
Каждая из трех искомых функций представляет собой линейную комбинацию бегущих волн, заданных произвольными непрерывно-дифференцируемыми функциями, которые
распространяются без изменения формы со скоростями, равными собственным значениям
матрицы-коэффициента A, задающего систему (8).
Наука и Образование. МГТУ им. Н.Э. Баумана
393
4. Распространение малых возмущений при относительном движении
В разделе 3 получено решение задачи для абсолютного движения газа. Поскольку абсолютное по отношению к неподвижной системе координат движение газа состоит из переносного и относительного движения, следует рассмотреть бегущие волны, вызванные
малыми возмущениями при относительном движении газа.
Пусть возмущения привносятся в газ с элементами (u ? U0 , ?0 , I0 ) [13]. Этот случай
соответствует распространению возмущений в системе координат, которая связана с частицей газа, движущейся со скоростью U0 . Тогда в условии (6) скорость возмущенного газа
совпадает с возмущением его скорости:
u(t, x) = u? (t, x).
(19)
Выполнив линеаризацию системы (5) около состояния (0, ?0 , I0 ), по аналогии с системой (8)
получим систему
0
1
?
?M 2
?1
0
0
?(? ? 1)
0
0
?
~
~
?B
?B
E
= A?
,
?t
?x
?
A =?
?
?
1
?
?M 2 ?
?
(20)
?,
?
где в силу близости физических процессов матрицы A? и A имеют совпадающие элементы.
Матрица A? является невырожденной, имеет три различных действительных собственных значения
?1 = 0,
?2 = 1/M,
?3 = ?1/M,
и подобна диагональной жордановой матрице J ? c матрицей подобия S ? , J ? = S ? ?1 A? S ? ,
где
?
0
?
S? = ? 1
?
?1/M
1
?
1/M ?
,
1 ?
?
?
S ? ?1
?1 (? ? 1) (? ? 1)
?
(? ? 1)/? ?1/? ?
?
= ? ?M/2
.
1/2?
1/2? ?
?
?
M/2
1/2?
1/2?
0
Если возмущения рассматриваются при относительном движении, то инварианты принимают вид
(? ? 1)
1
?(t, x) ? I(t, x) = C,
?
?
M
1
1
1
? u(t, x) + ?(t, x) + I(t, x) = f2 x +
t ,
2
2?
2?
M
M
1
1
1
u(t, x) + ?(t, x) + I(t, x) = f3 x ?
t ,
2
2?
2?
M
где C | связанная с интегрированием произвольная постоянная. Скорость, плотность и
энтальпия представляют собой линейную комбинацию бегущих волн, распространяющихся
Наука и Образование. МГТУ им. Н.Э. Баумана
394
с характеристическими скоростями 1/M и ?1/M :
?
1
1
1
1 ?
u(t, x) = ? f2 x +
t + f3 x ?
t ,
?
M
M
M
M
?
1
1
?(t, x) = C + f
x+
x?
t +f
t ,
2
3
?
M
M
?
1
1
? I(t, x) = ?C + (? ? 1)f2 x +
t + (? ? 1)f3 x ?
t .
M
(21)
M
Изменение (19) в условиях на компоненты возмущенного газа (6) привело к следующему
соотношению между матрицами, задающими результирующую линейную систему:
A? = A + E,
(22)
где E | единичная матрица. Схожие структура и свойства матриц A и A? , отвечающих
системам (8) и (20), приводят к схожести решений этих систем (18) и (21), построенных при
помощи элементарных бегущих волн. Возникает вопрос: можно ли отразить данный факт в
решении, не обращаясь к элементарным волнам и, следовательно, к отысканию собственных
значений и собственных векторов каждой из матриц A и A? в отдельности (см. раздел 6).
5. n-мерные бегущие волны
При изложении метода решения системы (8), отвечающей распространению малых возмущений в идеальном газе при одномерном течении, нигде не использовались размерность
искомой функции и порядок матрицы. Следовательно, решение произвольной гиперболи~ x) | n-мерная вектор-функция, A | числовая матрица
ческой системы вида (8), где B(t,
порядка n, может быть представлено как линейная комбинация n элементарных бегущих
волн.
Наряду с элементарной волной
x + ?i t,
(23)
i = 1, n,
(здесь ?i | собственные значения матрицы A), которая отвечает дифференциальному оператору
?
?
? ?i
(24)
?t
?x
(возникающему в системе независимых уравнений, содержащих одну функцию), рассмо??i ?
трим выражения
Ex + J t
и
Ex + A t,
(25)
которые отвечают дифференциальным операторам
?J ? E
?
?
?J
?t
?x
и
Наука и Образование. МГТУ им. Н.Э. Баумана
?A ? E
?
?
?A .
?t
?x
(26)
395
Размерность функций, на которые действуют операторы, задана размерностью коэффициента при производных: оператор (24) действует на скалярную функцию Vi , операторы (26)
~ соответственно.
действуют на n-мерные вектор-функции V~ и B
По аналогии с бегущими волнами (23), выражения (25) рассмотрены в работах [7, 8, 9] как
n-мерные бегущие волны (или бегущие волны с матричными коэффициентами), отвечающие
системам дифференциальных уравнений в частных производных 1-го порядка (8) и (10). К
волнам (23) в некоторых случаях будем добавлять определение "скалярные". В этом смысле
~ i будем применять определение <скалярный>, а выражение V~ ,
к инвариантам Римана (S ?1 B)
~ рассматривать как n-мерный инвариант.
V~ = S ?1 B,
По построению волны (25) связаны посредством матрицы S. Имеет место следующее
утверждение.
Утверждение 1. Если произвольная числовая матрица A порядка n подобна матрице J
с матрицей подобия S, то n-мерная бегущая волна (Ex + At), где E | единичная матрица
порядка n, подобна бегущей волне (Ex + Jt) с матрицей подобия S ?1 :
J = S ?1 AS
=?
(Ex + At) = S(Ex + Jt)S ?1 .
(27)
Если матрица J является диагональной жордановой формой матрицы A, то из соотношения (27) следует представление n-мерной бегущей волны через скалярные волны (x + ?i t),
построеннные при помощи собственных значений матрицы A, i = 1, n, [10].
Отметим, что модель одномерного движения идеального газа является примером модели,
в которой имеет место полная интерпретация бегущей волны с матричным коэффициентом
посредством скалярных волн. Данный факт отражает гиперболичность системы уравнений (8) (системы (7)).
6. Решение нормальной системы линейных однородных уравнений
в частных производных 1-го порядка, заданной числовой матрицей
Пусть в (12) функции fi , fi ? C 1 , заданы разложением в ряд
fi (x + ?i t) =
?
X
(x + ?i t)? ??i ,
i = 1, 3,
(28)
?=0
по своим волнам-аргументам, где ??i | числовые коэффициенты. Это означает, что компоненты искомой вектор-функции
т
V~ = (f1 (x + ?1 t), f2 (x + ?2 t), f3 (x + ?3 t))
представлены как разложение по соответствующим скалярным бегущим волнам. Запишем
данный факт посредством n-мерной волны, отвечающей диагональной матрице J:
V~ ? V~ (t, x) =
?
X
(Ex + J t)? ?~? ,
т
?~? = (??1 , ??2 , ??3 ) .
(29)
?=0
Наука и Образование. МГТУ им. Н.Э. Баумана
396
Если выражения (28) и (29) являются решениями уравнений (11) и (10), то в качестве решения
системы (8) естественно рассмотреть выражение
~ ? B(t,
~ x) =
B
?
X
(Ex + A t)?~?? ,
т
(30)
~?? = (??1 , ??2 , ??3 ) ,
?=0
где ~?? | числовой вектор-коэффициент. По построению каждый член ряда (28), (29) и (30),
а также любая частичная сумма этих рядов являются решением соответствующих уравнений [7].
Особенностью построенных решений (28){(30) является то, что однозначно задающий
уравнение или систему коэффициент соответствующей размерности при производной по
пространственной переменной x входит в решение в качестве коэффициента при независимой переменной t.
По отношению к системам (8) и (20) это, в частности, означает, что что решение одномерной задачи о распространении малых возмущений в идеальном газе при относительном
движении может быть получено (с точностью до ~?? ) из решения задачи при абсолютном
движении при замене в (30) матрицы A на матрицу A? , A? = A + E.
Построение решения (30) не требует нахождения собственных чисел и собственных векторов матрицы A. Вместе с тем, их можно ввести в (30) посредством отношения подобия (27)
между трехмерными волнами (Ex + At) и (Ex + Jt):
~ x) =
B(t,
?
X
S(Ex + J t)? S ?1~?? ,
J = S ?1 AS.
(31)
?=0
Поскольку выражения (30) и (31) могут быть поставлены в соответствие системе (8)
(системе (2)) для произвольной числовой матрицы A, условие гиперболичности, которое для
рассматриваемой физической модели является следствием законов сохранения, не является
существенным при любом его определении.
Поскольку при построении выражений (30) и (31) не используется размерность n функ~ x), каждое из них является общим решением системы (2) и частным решением
ции B(t,
системы (1) для произвольного n.
Будем говорить о волновом подходе к линейной системе уравнений в частных производных 1-го порядка (8), которая рассматривается в качестве частного случая системы (1),
понимая под этим разложение решения по n-мерным бегущим волнам.
7. Представление решения задачи Коши
Напомним [3], что система
E
~ x)
~ x)
? B(t,
? B(t,
=A
,
?t
?x
~ x) = (B1 (t, x), . . . , Bn (t, x)),
B(t,
x ? R1 ,
(32)
~ x) | n-мерная вектор-функция, A | произвольная числовая матрица с элеменгде B(t,
тами {aij }ni,j=1 , является системой Ковалевской. Для такой системы естественным образом
возникает задача Коши.
Наука и Образование. МГТУ им. Н.Э. Баумана
397
Пусть n-мерная вектор-функция ?
~ (x) является аналитической функцией x в некоторой
области G0 плоскости t = 0, содержащей точку x0 , x0 = 0, и представлена в G0 рядом
Тейлора
?
X
?
~ (x) =
?
~c? · x ,
?=0
1 ??
~c? =
?
~ (x)
.
?! ?x?
x=0
(33)
Рассмотрим начальное условие по переменной t:
~ x)
B(t,
x=0
т
=?
~ (x),
(34)
?
~ (x) = (?1 (x), . . . , ?n (x)) .
Согласно теореме Коши | Ковалевской задача Коши (32){(34) имеет в области G, G0 ?
G, единственное аналитическое решение, которое выражается степенным рядом (например, [3]):
?
X
Bi (t, x) =
(35)
??i 1 , ?2 x?1 t?2 .
?1 ,?2 =0
Числовые коэффициенты ??i 1 ,?2 однозначно определяются по начальным данным из системы
уравнений (32) в результате реккурентной процедуры по индексу ?2 и являются полиномами
элементов матрицы A и коэффициентов ??i 1 ,l , l ? ?2 ? 1.
Рассмотрим понятие бегущей волны в применении к решению задачи Коши. Начнем со
скалярного случая. Пусть n = 1. Тогда, поскольку уравнение
?B(t, x)
?B(t, x)
=a
,
?t
?x
a11 = a,
B1 (t, x) = B(t, x),
(36)
имеет решение B = f (x + at), где f ? C 1 , для решения задачи Коши в качестве функции f
следует выбрать функцию, которая задает начальное условие [4]:
B(t, x) = ?(x + at),
(37)
?(x) = ?1 (x).
Очевидно (например, [2, 6]), что для гиперболической системы (32) вид функций fi , i = 1, n,
которые определяют сохраняющуюся форму i-й бегущей волны, устанавливается при помощи матрицы подобия:
fi (0, x) = (S ?1 ?
~ (x))i ,
S ?1 AS = J,
J = diag(?1 , . . . , ?n ).
(38)
Тогда при замене в i-й компоненте вектора ?
~ (x) аргумента x на волну (x + ?i t) из (9)
получим следующее решение решение задачи Коши для системы (8), соответствующей
распространению малых возмущений в идеальном газе:
u(t, x) =
3
X
i=1
s1i ?i (x + ?i t),
?(t, x) =
3
X
s2i ?i (x + ?i t),
i=1
I(t, x) =
3
X
s3i ?i (x + ?i t),
(39)
i=1
где {sij }3i,j=1 | элементы матрицы S, столбцами которой являются собственные векторы
матричного коэффициента этой системы. Решение (39) указывает на то, каким образом
в точку (t, x) вдоль проходящих через нее трех характеристик (13) приносятся значения
соответствующих функций ?i (x), i = 1, 3.
Наука и Образование. МГТУ им. Н.Э. Баумана
398
Например, выберем начальные условия в виде
т
?
~ (x) = (? sin x, 0, µ cos x) ,
?, µ 1,
(40)
?? < x < +?,
где ? и µ | малые параметры, задав начальные возмущения скорости и температуры (посредством энтальпии) при одномерном движении идеального газа. Из (39) получим следующее
представление решения задачи Коши:
?
(1 ? M )
µ
(1 ? M )
?
?
t ?
cos x +
t +
u(t, x) = sin x +
?
2
M
2?M
M
?
?
(1 + M )
µ
(1 + M )
?
+ sin x ?
t +
cos x ?
t ,
2
M
2?M
M
?
?
? ?(t, x) = ? µ cos(x ? t) ? ?M sin x + (1 ? M ) t +
?
2
M
?
?
µ
(1 ? M )
?M
(1 + M )
µ
(1 + M )
+
cos
x
+
t
+
sin
x
?
t
+
cos
x
?
t
,
?
2?
M
2
M
2?
M
?
µ
?(? ? 1)M
(1 ? M )
?
I(t, x) = cos(x ? t) ?
sin x +
t +
?
?
2
M
?
µ(? ? 1)
(1 ? M )
?(? ? 1)M
(1 + M )
?
+
cos x +
t +
sin x ?
t +
2?
M
2
M
?
?
µ(? ? 1)
(1 + M )
?
+
cos x ?
t .
2?
(41)
M
В (41), в частности, представлено, каким образом (с какой амплитудой и с какой частотой) начальное возмущение скорости газа ? sin x скажется на возмущении его плотности и
температуры (посредством энтальпии).
Воспользуемся теперь разложением начальной аналитической функции в ряд Тейлора (33). Пусть функции fi , сохраняющие свои значения на волнах-аргументах (x + ?i t),
представляют собой разложение (28). Тогда при переходе к матричной записи получим
разложение n-мерного инварианта V~ в ряд (29) по n-мерным бегущим волнам. Из (9)
следует начальное условие на функцию V~ :
V~ (t, x)
t=0
=
?
X
(Ex +
?=0
Jt)? ?~?
t=0
=S
?1
?
X
~c? · x? .
?=0
Следовательно, имеет место соотношение
~т
?1
?? = S ~c? ,
1 ??
?
~ (x)
,
~c? =
?
?! ?x
x=0
(42)
между коэффициентами разложения по n-мерным бегущим волнам векторного инварианта V~
и коэффициентами разложения в ряд Тейлора начальной вектор-функции ?
~ (x).
В то же время, рассмотрев при t = 0 общее решение системы (32) в виде (30), приходим
к выводу, что для получения решения задачи Коши (32){(33) в качестве коэффициентов ~??
Наука и Образование. МГТУ им. Н.Э. Баумана
399
при (Ex + At)? следует выбрать коэффициенты ряда (33) при x? . Тогда решение задачи
Коши принимает вид [8]
~ x) =
B(t,
?
X
(Ex + At)?
?=0
1 ??
?
~ (x) .
?! ?x?
x=0
(43)
Решение (43), которое записано без дополнительных вычислений, содержит матричный коэффициент системы в явном виде, а начальные данные | в виде коэффициентов разложения
аналитической начальной функции в ряд Тейлора. Ряд (43) организован по волновому принципу, отличному от принципа организации ряда (35).
Например, чтобы получить решение рассматриваемой задачи Коши для уравнений распространения малых возмущений в идеальном газе при одномерном движении, следует
подставить в (43) матричный коэффициент A системы (8) и векторные коэффициенты
т
~?0 = (0, 0, µ) ,
µ т
~?2 = 0, 0, ?
,
2!
т
~?1 = (?, 0, 0) ,
т
?
~?3 = ? , 0, 0 ,
3!
...,
которые отвечают разложению функции (40) в ряд Тейлора в окрестности нуля. Сравнив ком~ x), B(t,
~ x) = (u(t, x), ?(t, x), I(t, x))т , с решением (41) при
поненты полученного вектора B(t,
разложении тригонометрических функций в ряд Тейлора, убеждаемся в их тождественности.
Напомним при этом, что любая частичная сумма ряда (43) с указанными коэффициентами A
и ~?? является решением системы (8).
Запишем систему уравнений (32) в операторном виде:
E
~ x)
? B(t,
~ x),
= A0 B(t,
?t
A0 = A
?
.
?x
(44)
По аналогии с обыкновенным дифференциальным уравнением представим решение задачи
Коши для системы (44) через операторную экспоненту:
~ x) = exp(A0 t) ?
B(t,
~ (x).
(45)
Под результатом действия оператора exp(A0 t) на аналитическую функцию ?
~ (x) понимают
функцию, к которой сходится степенной ряд
exp(A0 t) ?
~ (x) =
? i
X
t
i=0
i!
(A0 )i ?
~ (x),
exp(A0 t)
t=0
(46)
= E.
Приравняв формальное решение задачи Коши в виде операторной экспоненты (45) и найденное решение (43), получим формулу
?
?
X
1 ??
? ?1 1 ?
exp(A t)~
?(x) ?
(Ex + At)
?
~
(x)
=
S(Ex
+
Jt)
S
?
~ (x) , (47)
?! ?x?
?! ?x?
x=0
?=0
?=0
x=0
0
?
X
?
где
A0 = A
?
,
?x
Наука и Образование. МГТУ им. Н.Э. Баумана
J = S ?1 AS.
400
Таким образом, в рамках модели волнового взаимодействия решение задачи Коши представлено как разложение в ряд (сходящийся в некоторой окрестности начальной точки) по
соответствующей задаче n-мерной бегущей волне. Векторные коэффициенты ряда и матричный коэффициент волны известны. Общий член ряда может быть представлен через
собственные значения и собственные векторы матричного коэффициента волны.
Отметим также, что задача Коши для уравнения (44) позволяет трактовать n-мерные
бегущие волны как результат действия соответствующей этому уравнению операторной
экспоненты exp(A0 t) на n-мерную аналитическую функцию ?
~ (x).
Разложение решения задачи Коши (32){(34) в ряд (47) имеет место как для гиперболической системы (матрица A подобна диагональной жордановой матрице J), так и для
системы (32), заданной произвольной матрицей (A подобна жордановой матрице J).
8. Заключение
В рамках предложенной авторами модели волнового взаимодействия в качестве тестовой
задачи сформулирована и решена задача о распространении малых возмущений в идеальном
газе при одномерном движении.
С использованием традиционного подхода к решению гиперболической системы малые
возмущения скорости, плотности и энтальпии газа единообразно представлены как линейная
комбинация произвольных непрерывно-дифференцируемых функций, аргументами которых
являются элементарные бегущие волны.
Число элементарных волн и число инвариантов Римана, отвечающих системе, равно
трем | порядку задающей систему матрицы, что рассматривается как одно из определений гиперболичности. Три пересекающиеся характеристики образуют на плоскости (t, x)
характеристический треугольник. Это означает, что в идеальном газе при малых возмущениях скорости, плотности и энтальпии (температуры) возникает движущаяся замкнутая
структура, которая однозначно задана параметрами газа (числом Маха).
Представление о бегущей волне как решении гиперболической системы, которое обладает свойством инвариантности, перенесено на n-мерный (матричный) аналог и связано с
решением нормальной системы уравнений в частных производных 1-го порядка, которая
задана произвольной числовой матрицей. Указано представление n-мерной волны через
элементарные, которые рассмотрены как скалярные волны (n = 1).
Решение задачи Коши в окрестности точки x0 , x0 = 0, представлено как разложение
в ряд по n-мерной бегущей волне с коэффициентами, равными n-мерным векторным коэффициентам разложения аналитической начальной функции в ряд Тейлора. Структура
построенного ряда отлична от структуры ряда, фигурирующего в теореме КошиЁ| Ковалевской. Установлено, что n-мерные бегущие волны, отвечающие системе (1) при s = 1
(системе (2)), порождаются действием отвечающей этой системе операторной экспоненты
на аналитическую функцию.
Наука и Образование. МГТУ им. Н.Э. Баумана
401
В продолжение исследований в качестве физической модели для системы (1) при s > 1
предлагается рассмотреть задачу о распространении малых возмущений в идеальном газе
при двумерном (плоскопараллельном) движении и движении в пространстве. Это потребует
применения соответствующего оператора волнового взаимодействия [7, 8]. Взаимодействие
осуществляется между волнами (скалярными или n-мерными), 3 < n ? 5, построенными
для каждой матрицы-коэффициента Ak , s = 1, 3, системы (1). Полученное решение (30),
связанное с волнами, которые распространяются вдоль выбранной оси OXk , k = 1, 3, является частным решением задачи о распространении малых возмущений в идеальном газе при
движении в пространстве и составной частью ее общего решения.
Список литературы
1. Годунов C.К. Уравнения математической физики. М.: Наука, 1979. 392 c.
2. Куликовский А.Г., Свешникова Е.И., Чугайнова А.П. Лекционные курсы НОЦ. Вып. 16.
Математические методы изучения разрывных решений нелинейных гиперболических
систем уравнений / Математический институт им. В.А. Стеклова РАН (МИАН). М.:
МИАН, 2010. 122 с.
3. Курант Р. Уравнения с частными производными: пер. c англ. / под ред. О.А. Олейник. М.:
Мир, 1964. 830 с. [Courant R. Methods of mathematical physics. Vol. 2. Partial Differential
Equations. New York-London: Interscience Publ., 1962].
4. Бхатнагар П. Нелинейные волны в одномерных дисперсных системах: пер. с англ. /
под ред. П.Е. Краснушкина, Н.Р. Сибгатуллина. М.: Мир, 1983. 136 c. [Bhatnagar P. L.
Nonlinear Waves in One-dimensional Dispersive Systems. Oxford: Clarendon Press, 1979].
5. Лойцянский Л. Г. Механика жидкости и газа. М.: Наука, 1987. 840 с.
6. Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer Berlin
Heidelberg, 2009. 724 p. DOI: 10.1007/b79761
7. Феоктистов В.В., Мякинник О.О. Структура ряда для решения системы уравнений с
частными производными 1-го порядка // Вестник МГТУ им. H.Э. Баумана. Сер. Естественные науки. 2009. № 4. С. 3{22.
8. Феоктистов В.В., Мякинник О.О. Оператор волнового взаимодействия и нормальная
форма системы линейных уравнений в частных производных 1-го порядка // Воронежская весенняя математическая школа <Современные методы теории краевых задач. Понтрягинские чтения | XXI>: матер. Воронеж: Изд-во ВГУ, 2010. С. 232{233.
9. Феоктистов В.В., Мякинник О.О. Бегущие волны с матричными коэффициентами //
Воронежская зимняя математическая школа <Современные методы теории функций и
смежные проблемы>: матер. Воронеж: Изд-во ВГУ, 2013. С. 255{256.
Наука и Образование. МГТУ им. Н.Э. Баумана
402
10. Феоктистов В.В., Мякинник О.О. Представление n-мерных бегущих волн через скалярные // Воронежская весенняя математическая школа <Современные методы теории
краевых задач. Понтрягинские чтения | XXV>: матер. Воронеж: Изд-во ВГУ, 2014.
С. 178{180.
11. Geroch R. Partial Differential Equations of Physics // In: General Relativity: Proceedings of
the 46th Scottish Universities Summer School in Physics (Aberdeen, July 1995). Edinburgh:
SUSSP Publ., 1996. P. 19{60.
12. Ковеня В.М., Яненко Н.Н. Метод расщепления в задачах газовой динамики. Новосибирск: Наука, 1981. 304 c.
13. Соболев С. Л. Уравнения математической физики. М.: Наука, 1992. 432 c.
14. Капцов О.В. Инварианты характеристик систем уравнений с частными производными //
Сибирский математический журнал. 2004. Т. 45, № 3. C. 577{591.
Наука и Образование. МГТУ им. Н.Э. Баумана
403
Science and Education of the Bauman MSTU,
2014, no. 12, pp. 372{391.
DOI: 10.7463/1214.0744956
Received:
Revised:
28.06.2014
02.12.2014
c Bauman Moscow State Technical University
Ideal Gas Small Perturbations:
Waves, Invariants and Cauchy Problem
Feoktistov V. V., Miakinnik O. O.
*
Bauman Moscow State Technical University
Keywords: hyperbolic system, traveling waves, Riemann invariants, n-dimensional traveling waves,
expansion in terms of traveling waves, Cauchy problem solution, ideal gas, small perturbations
The wave interaction model the authors have developed for the first-order linear homogenous
Cauchy system, which is specified by numerical matrices at the spatial derivatives, is considered
as applied to the one-dimensional ideal gas flow. The linear waves, meeting model requirements,
are separated from quasi-linear process by the method of small perturbations. The resulting
linear hyperbolic system is specified by the matrix containing gas parameters and written for
the unknown vector, components of which are small perturbations of gas velocity, density, and
enthalpy.
Following a common practice for a hyperbolic system, the method for solving the system via
traveling waves, propagated axially with characteristic speeds, is analyzed for the physical process
involved. Three Riemann invariants are found. The general solution is represented as a linear
combination of three arbitrary functions of wave-arguments, specified by characteristic speeds.
The initial conditions are imposed on the three functions, which can be measured experimentally,
to determine the unchangeable form of traveling waves. The formation of characteristic triangle,
uniquely determined by gas parameters, is connected with the existence of some closed structure
moving in the ideal gas.
To hold the matrix form, a 3-by-3 diagonal matrix is formatted from the traveling waves as its
entries. This matrix is considered as a three-dimensional wave, specified by coefficient matrix with
3 characteristic speeds along its diagonal. Using a similarity transformation, the diagonal matrix
coefficient is replaced by the matrix specifying the system to be solved.
By analogy with hyperbolic case, an n-dimensional wave, specified by an arbitrary numerical nby-n matrix, is considered as a solution of a first-order Cauchy system specified by this matrix. The
Cauchy problem solution is represented as a series expansion in terms of n-dimensional wave with
coefficients equal to the coefficients of n-dimensional initial analytic vector Taylor's expansion.
Science and Education of the Bauman MSTU
404
It was found that an n-dimensional wave is generated by the action of differential operator
exponential on an n-dimensional analytic initial vector via Cauchy problem for the related system.
It is proposed to apply the model to the two-dimensional gas flow using the notion of interaction
between the waves propagated along each of two coordinate axes.
References
1. Godunov S.K. Uravneniya matematicheskoy fiziki [Equations of Mathematical Physics].
Moscow, Nauka Publ., 1979. 392 p. (in Russian).
2. Kulikovskiy A.G., Sveshnikova E.I., Chugaynova A.P. Lektsionnye kursy NOTs. Vyp. 16.
Matematicheskie metody izucheniya razryvnykh resheniy nelineynykh giperbolicheskikh sistem
uravneniy [Lecture courses of Scientific and Educational Center. Is. 16. Mathematical methods
for studying discontinuous solutions of nonlinear hyperbolic systems of equations]. Moscow,
Steklov Mathematical Institute of RAS Publ., 2010. 122 p. (in Russian).
3. Courant R. Methods of mathematical physics. Vol. 2. Partial Differential Equations. New
York-London: Interscience Publ., 1962. (Russ. ed.: Courant R. Uravneniya s chastnymi
proizvodnymi. Moscow, Mir Publ., 1964. 830 p.).
4. Bhatnagar P.L. Nonlinear Waves in One-Dimensional Dispersive Systems. Oxford, Clarendon Press, 1979. (Russ. ed.: Bhatnagar P.L. Nelineynye volny v odnomernykh dispersnykh
sistemakh. Moscow, Mir Publ., 1983. 136 p.).
5. Loytsyanskiy L.G. Mekhanika zhidkosti i gaza [Fluid Mechanics]. Moscow, Nauka Publ.,
1987. 840 p. (in Russian).
6. Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer Berlin
Heidelberg, 2009. 724 p. DOI: 10.1007/b79761
7. Feoktistov V.V., Myakinnik O.O. Structure of Series for Solving System of 1st Order Equations
with Partial Derivatives. Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki = Herald
of the Bauman MSTU. Ser. Natural science, 2009, no. 4, pp. 3{22 (in Russian).
8. Feoktistov V.V., Myakinnik O.O. Operator of wave interaction and normal form of system of
linear equations in partial derivatives of the 1st order. Voronezhskaya vesennyaya matematicheskaya shkola \Sovremennye metody teorii kraevykh zadach. Pontryaginskie chteniya | 21":
mater. [Voronezh Spring School in Mathematics \Modern methods of the theory of boundary
value problems. Pontryagin reading | 21": proc.]. Voronezh, VSU Publ., 2010, pp. 232{233
(in Russian).
9. Feoktistov V.V., Myakinnik O.O. Traveling waves with matrix coefficients. Voronezhskaya
zimnyaya matematicheskaya shkola \Sovremennye metody teorii funktsiy i smezhnye problem": mater. [Voronezh Winter School in Mathematics \Modern methods of function theory
and related problems"]. Voronezh, VSU Publ., 2013, pp. 255{256 (in Russian).
Science and Education of the Bauman MSTU
405
10. Feoktistov V.V., Myakinnik O.O. N -dimensional traveling waves representation in terms
scalar waves. Voronezhskaya vesennyaya matematicheskaya shkola \Sovremennye metody
teorii kraevykh zadach. Pontryaginskie chteniya | 25": mater. [Voronezh Spring School
in Mathematics \Modern methods of the theory of boundary value problems. Pontryagin
reading | 25": proc.]. Voronezh, VSU Publ., 2014, pp. 178{180 (in Russian).
11. Geroch R. Partial Differential Equations of Physics. In: General Relativity: Proceedings of
the 46th Scottish Universities Summer School in Physics, Aberdeen, July 1995. Edinburgh,
SUSSP Publ., 1996, pp. 19{60.
12. Kovenya V.M., Yanenko N.N. Metod rasshchepleniya v zadachakh gazovoy dinamiki [Splitting
Methods for Gas Dynamics Problems]. Novosibirsk, Nauka Publ., 1981. 304 p. (in Russian).
13. Sobolev S.L. Uravneniya matematicheskoy fiziki [Equations of Mathematical Physics].
Moscow, Nauka Publ., 1992. 432 p. (in Russian).
14. Kaptsov O.V. Invariants of characteristics of some systems of partial differential equations. Sibirskiy matematicheskiy zhurnal, 2004, vol. 45, no. 3, pp. 577{591. (English
translation:
Siberian Mathematical Journal, 2004, vol. 45, is. 3, pp. 475{487. DOI:
10.1023/B:SIMJ.0000028612.51060.32).
Science and Education of the Bauman MSTU
406
ия, являются уравнениями гиперболического типа. Свойство гиперболичности тесно связано с существованием решений специального вида, которые обладают свойством инвариантности.
Такие решения называют бегущими волнами [1, 2, 3, 4].
Напомним, что гиперболическая система n, n ? 1, однородных дифференциальных уравнений в частных производных 1-го порядка относительно независимых переменных (t, x),
Наука и Образование. МГТУ им. Н.Э. Баумана
387
x = xs , и однородное уравнение в частных производных порядка m, m ? 1, относительно
переменных (t, x), x = (x1 , . . . , xs ), s ? 1, имеют решения, которые сохраняют свои значения на характеристических прямых и гиперплоскостях соответственно [2, 3]. Случай n = 1,
s = 1 и m = 1, s = 1 соответствует простейшему уравнению в частных производных,
решением которого является элементарная волна [4, 5].
Для построения волны следует определить коэффициенты при независимых переменных
в выражении (?0 t+?·x), ?·x = ?1 x1 +. . .+?s xs , решив для этого алгебраическое уравнение,
построенное по определенному правилу из числовых коэффициентов дифференциального
уравнения или системы уравнений.
Хорошо известными примерами являются: система уравнений акустики (n = 2) и ее
решение в виде элементарных волн (xs ? at) и (xs + at), распространяющихся со скоростью
звука a вдоль координатной оси OXs [1, 2, 6]; волновое уравнение (m = 2) в пространстве s
переменных, которое порождает волну f (?0 t + ? · x) (здесь f | произвольная гладкая
функция), распространяющаяся в направлении вектора (?0 , ?) [1, 3].
В работах [7, 8, 9] для системы
s
~
~
X
?B
?B
E
=
Ak
,
?t
?xk
k=1
~ = (B1 (t, x), . . . Bn (t, x))т ,
B
x = (x1 , . . . , xs ),
(1)
где Ak , k = 1, s, | числовые квадратные матрицы порядка n, E | единичная матрица,
безотносительно к определяемому матрицами-коэффициентами типу рассмотрены решения
(Exk +Ak t), k = 1, s, и функции от них. По аналогии с элементарными волнами такие решения названы бегущими волнами с матричными коэффициентами или n-мерными бегущими
волнами. При s = 1 в случае гиперболичности системы (матрица As подобна диагональной жордановой матрице) введенная n-мерная бегущая волна имеет полную интерпретацию
через n элементарных волн, которые рассмотрены как скалярные [10].
Несмотря на существование численных методов решения системы (1) в приложении к
отдельным физическим процессам (например, [6]), актуальным является построение аналитического метода решения этой системы и связанной с ней задачи Коши при произвольных n, s и произвольных матричных коэффициентах Ak , k = 1, s. Для системы (1), которая
характерна для различных областей математической физики [11], в работах [7, 8, 9] авторами построена модель волнового взаимодействия. Элементами модели являются: n-мерные бегущие волны, содержащие коэффициенты системы Ak , k = 1, s, которые распространяются по координатным осям и находятся во взаимодействии; оператор волнового
взаимодействия, характеризующийся размерностью s; отвечающая системе операторная
экспонента.
В рамках модели общее решение системы (1) представлено как разложение в ряд по взаимодействующим бегущим волнам с неопределенными числовыми вектор-коэффициентами.
Найденное решение в явном виде содержит матрицы-коэффициенты Ak , k = 1, s, задаюНаука и Образование. МГТУ им. Н.Э. Баумана
388
щие систему. Для получения решения задачи Коши достаточно в общем решении заменить
неопределенные коэффициенты на векторные коэффициенты разложения аналитической начальной функции в ряд Тейлора [8].
В качестве физического процесса для интерпретации введенных в [7, 8, 9] понятий и выполненных построений выберем течение идеального газа, линейные и нелинейные волновые
процессы в котором изучены [1, 5, 12]. Выбранному процессу соответствует квазилинейная система дифференциальных уравнений в частных производных 1-го порядка [2, 5, 12],
которые представляют собой законы сохранения. В качестве замыкающих уравнений [12]
для уранений сохранения скорости и импульса рассмотрим уравнение сохранения энергии,
выразив ее через энтальпию, и уравнение состояния, связывающее энтальпию с плотностью и давлением газа. Таким образом, результирующая квазилинейная система записана
относительно n-мерного вектора, 3 ? n ? 5, единообразно содержащего в качестве своих
компонент компоненты скорости газа u(t, x), v(t, x), w(t, x), его плотности ?(t, x) и энтальпии I(t, x), где x = (x1 , . . . , xs ), s ? 3. Поскольку энтальпия посредством коэффициента
~ x)
теплоемкости связана с температурой газа, каждая из компонент искомого вектора B(t,
является величиной, которая может быть измерена в ходе эксперимента.
Чтобы из квазилинейного процесса выделить линейные волны, порождаемые системой (1), применим метод малых возмущений (например, [2]). Возникающую при этом
задачу о распространении малых возмущений, которая представляет самостоятельный интерес [2, 5, 13], рассмотрим как тестовую для модели волнового взаимодействия на примере
движения идеального газа при заданных n и s.
По структуре построенной модели общее решение системы
E
~
~
?B
?B
= Ak
,
?t
?xk
~ = (B1 (t, xk ), . . . Bn (t, xk ))т ,
B
k = 1, s,
(2)
представленное как разложение по волне (Exk + Ak t), является частным решением системы (1) и составной частью ее общего решения, представленного в виде разложения по s
взаимодействующим волнам [7, 9].
Пусть (2) отвечает распространению малых возмущений в идеальном газе при одномерном течении. Найдем решение этой системы с использованием традиционного подхода к
одномерной гиперболической системе [2, 6] и средствами модели волнового взаимодействия
без использования условия гиперболичности [7]. Покажем тождественность полученных решений и преимущества волнового подхода при исследовании системы (2).
Целью работы является рассмотрение предлагаемого авторами нового аналитического
метода решения нормальной системы линейных однородных уравнений в частных производных 1-го порядка на примере задачи о распространении малых возмущений в идеальном
газе при одномерном течении с последующим использованием полученного решения для
построения решения задачи при течении газа в пространстве.
Наука и Образование. МГТУ им. Н.Э. Баумана
389
2. Уравнения распространения малых возмущений в идеальном газе
Квазилинейные уравнения одномерного движения идеального газа при отсутствии массовых сил в декартовой системе координат имеют вид [2, 5, 13]:
?
? ?u
?
?u?
1 ? P?
+ u?
=?
,
?
? x?
?? ? x?
?
t?
?
?
? ??
?(?? u?)
? t?
? x?
? ? t? + ? x? = 0,
?
?
?
? ?I + u? ?I = 1 ? P? + u? ? P? ,
?
??
x? ? R1 ,
(3)
? x?
? t?
? t?x?) | скорость движения газа, его плотность, давление и
где u?(t?, x?), ??(t?, x?), P? (t?, x?), I(
энтальпия соответственно. Чтобы сделать систему трех дифференциальных уравнений относительно четырех неизвестнах функций определенной, дополним ее уравнением состояния [12]
??1 ? ?
cP
??I, I = cP T? , ? =
,
(4)
?
cV
которое соответствует баротропному движению идеального газа и связывает давление газа с
P? =
его плотностью и энтальпией (T? (t?, x?) | температура газа, cP | удельная теплоемкость газа
при постоянном давлении, cV | удельная теплоемкость газа при постоянном объеме). Вос?
пользуемся уравнением состояния (4) и запишем систему относительно функций (u,
? ??, I):
?u?
?u? ? ? 1 ? ? ?? ? ? 1 ?I?
?
?
= ???u?
I
?? ,
? x?
?
? x?
?
? x?
? t?
? ??
?u?
? ??
? u? ,
= ???
? x?
? x?
? t?
?I?
?u?
?I?
? u? .
= ?(? ? 1)I?
? x?
? x?
? t?
??
Выразим зависимые и независимые переменные через одноименные безразмерные переменные:
u? = U0 u,
?? = ?0 ?,
I? = I0 I,
t? = T0 t,
x? = L0 x,
где U0 , ?0 , I0 , T0 , L0 | размерные постоянные, которые являются характерными для задачи
скоростью, плотностью, энтальпией, временем и длиной соответственно.
Запишем систему одномерного движения идеального газа относительно безразмерных
функций (u, ?, I):
?
?u
1
??
1
?I
?u
?
?
?
St ?
= ??u
I
? ,
2
2
?
?t
?x ?M
?x ?M
?x
?
?
??
?u
??
? St ?t = ?? ?x ? u ?x ,
?
?I
?u
?I
?
? St
?u .
= ?(? ? 1) I
?t
?x
Наука и Образование. МГТУ им. Н.Э. Баумана
(5)
?x
390
В (5) St | число Струхала; St = L0 /U0 T0 ; M | число Маха; M = U0 /a, a | скорость звука
в газовой среде, при этом a2 = (? ? 1)I0 . Поскольку среди размерных величин, задающих
число Струхала, две являются независимыми, можно выбрать St = 1.
Рассмотрим задачу о распространении малых возмущений в идеальном газе при одномерном параллельном оси OX движении [2, 5]. Пусть в движущемся со скоростью U0
идеальном газе, который имеет плотность ?0 и энтальпию I0 , создаются малые возмущения
скорости u? (t, x), плотности ?? (t, x) и энтальпии I ? (t, x). Тогда движению возмущенного
газа отвечают функции
u(t, x) = 1 + u? (t, x),
?(t, x) = 1 + ?? (t, x),
I(t, x) = 1 + I ? (t, x).
(6)
Чтобы нелинейная система (5) имела решение вида (6), функции, задающие малые возмущения, должны удовлетворять системе линейных уравнений [2]:
?
?u
1 ??
1 ?I
? ?u
?
?
=?
,
2
? ?t
?x ?M ?x ?M 2 ?x
?
??
?u
??
(7)
? ?t = ? ?x ? ?x ,
?
?I
?u ?I
?
?
= ?(? ? 1)
?t
?x ?x
(для упрощения обозначений звездочки при искомых функциях опущены). Выполненная
процедура называется линеаризацией системы (5) около состояния (U0 , ?0 , I0 ) [2].
3. Приведение к независимым волновым уравнениям и нахождение инвариантов
Запишем полученную линейную систему (7) в матричном виде относительно искомой
~ x) = (u(t, x), ?(t, x), I(t, x))т , x ? R1 :
вектор-функции B(t,
?
~
~
?B
?B
E
=A
,
?t
?x
?
A=?
?
?1
?
1
?M 2
?
1 ?
?M 2 ?
?1
?1
0
?(? ? 1)
0
?1
(8)
?,
?
где E | единичная матрица порядка 3. Система (8) представляет собой систему линейных
однородных дифференциальных уравнений в частных производных 1-го порядка с постоянными коэффициентами, которая записана в нормальной форме. Свойства системы и ее
решений определяются алгебраическими свойствами матрицы A.
Матрица A является невырожденной, имеет три различные действительные собственные
значения
1?M
1+M
, ?3 = ?
M
M
и подобна диагональной жордановой матрице J c матрицей подобия S (J = S ?1 AS), где
?1 = ?1,
?
0
?
S=? 1
?
?2 =
1
1
?
?M
M
?
?,
?
?1 ?(? ? 1)M (? ? 1)M
Наука и Образование. МГТУ им. Н.Э. Баумана
?
S ?1
?
0 (? ? 1)/?
?1/? ?
?
= ? 1/2 ?1/2?M ?1/2?M ? .
?
?
1/2 1/2?M
1/2?M
391
Матрица S составлена из собственных векторов матрицы A, матрица S ?1 найдена как обратная к S. Согласно определению (например, [2]) система (8) (система (7)) является гиперболической (в узком смысле).
Выполним замену переменных
~ = S V~ ,
B
т
V~ (t, x) = (V1 (t, x), V2 (t, x), V3 (t, x)) ,
(9)
и рассмотрим систему уравнений
E
? V~
? V~
= S ?1 AS
.
?t
?x
(10)
Система (10) состоит из независимых уравнений
?Vi
?Vi
= ?i
,
?t
?x
i = 1, 3,
(11)
каждое из которых содержит одну неизвестную функцию и имеет решение
Vi = fi (x + ?i t),
i = 1, 3,
(12)
где fi | произвольные непрерывно-дифференцируемые функции. Прямые
x + ?i t = c i ,
ci = const,
i = 1, 3,
(13)
являются характеристиками системы и при пересечении образуют в плоскости (t, x) треугольник, называемый характеристическим [1]. Основание треугольника, отвечающее прямой (13) при i = 1, не зависит от параметров газа. Данный факт согласуется с законом
распространения малых возмущений в покоящемся газе [1] и при относительном движении
газа, которое рассматривается в разделе 4. Напомним, что для уравнений акустики в качестве
основания характеристического треугольника выбирается интервал на оси OX [1].
Функции fi , i = 1, 3, сохраняют значения на соответствующих параллельных прямых (13) и вместе со своими аргументами называются бегущими волнами [1, 3]. Уравнение в
частных производных 1-го порядка с двумя независимыми переменными (10) моделирует одномерный волновой процесс (волновой процесс в одномерной среде) и является самым простым из уравнений, к которым применяется определение <волновое>. Поэтому аргументы
функций fi называют также элементарными волнами [4] (в [5] | простыми волнами). В [2]
для волн (12) применен термин "бегущие волны малых возмущений", который, отражая
связь этих волн с исходной квазилинейной системой и методом ее линеаризации, указывает на важность задачи о распространении малых возмущений для исследования волновых
процессов. Бегущие волны малых возмущений распространяются с характеристическими
скоростями.
Выполнение преобразования подобия (9) по отношению к гиперболической системе (8)
привело к независимым уравнениям и появлению элементарных бегущих волн в качестве
решения.
Наука и Образование. МГТУ им. Н.Э. Баумана
392
Найденные выражения
~ ,
Vi ? S ?1 B
i
i = 1, 3,
записанные через физически значимые функции, сохраняют свои значения на характеристиках и, согласно определению [1, 2], представляют собой инварианты Римана системы (8):
(? ? 1)
1
?(t, x) ? I(t, x) = f1 (x ? t),
?
?
1
1
1
(1 ? M )
u(t, x) ?
?(t, x) ?
I(t, x) = f2 x +
t ,
2
2?M
2?M
M
1
1
1
(1 + M )
u(t, x) +
?(t, x) +
I(t, x) = f3 x ?
t .
2
2?M
2?M
M
(14)
(15)
(16)
Каждый из трех инвариантов (14) { (16) распространяется без искажения формы с характеристической скоростью ?i , i = 1, 3.
В связи с постановкой задачи (3){(4) и нахождением инвариантов (14){(16) отметим
работу [14], в которой доказано, что система одномерных уравнений газовой динамики
имеет три иварианта Римана, если в качестве замыкающих выбрано уравнение сохранения
энергии, выраженное через энтропию s, и уравнение состояния в виде
V = G(P ) + A(s),
1
V = ,
?
G0 (P ) = ?
1
g 2 (P )
,
c=
g(P )
,
?
(17)
где a | скорость звука; g и A | произвольные функции давления и энтропии.
Из (9) и (14) { (16) следует, что решение нормальной линейной системы уравнений в частных производных 1-го порядка (8), которая отвечает распространению малых возмущений в
идеальном газе при одномерном движении, имеет вид
?
(1 ? M )
(1 + M )
?
u(t,
x)
=
f
x
+
t
+
f
x
?
t
,
2
3
?
M
M
?
(1 ? M )
(1 + M )
x?
t ,
M
t + M f3
? ?(t, x) = f1 (x ? t) ? M f2 x +
M
?
(1 ? M )
(1 + M )
? I(t, x) = ?f1 (x ? t) ? (? ? 1)M f2 x +
t + (? ? 1)M f3 x ?
t .
M
(18)
M
Отметим, что скорость распространения малых возмущений u(t, x) выражена через волны,
отвечающие двум из трех характеристик (13). Данный факт согласуется с представлением
скорости распространения звука, найденной из уравнений акустики.
Каждая из трех искомых функций представляет собой линейную комбинацию бегущих волн, заданных произвольными непрерывно-дифференцируемыми функциями, которые
распространяются без изменения формы со скоростями, равными собственным значениям
матрицы-коэффициента A, задающего систему (8).
Наука и Образование. МГТУ им. Н.Э. Баумана
393
4. Распространение малых возмущений при относительном движении
В разделе 3 получено решение задачи для абсолютного движения газа. Поскольку абсолютное по отношению к неподвижной системе координат движение газа состоит из переносного и относительного движения, следует рассмотреть бегущие волны, вызванные
малыми возмущениями при относительном движении газа.
Пусть возмущения привносятся в газ с элементами (u ? U0 , ?0 , I0 ) [13]. Этот случай
соответствует распространению возмущений в системе координат, которая связана с частицей газа, движущейся со скоростью U0 . Тогда в условии (6) скорость возмущенного газа
совпадает с возмущением его скорости:
u(t, x) = u? (t, x).
(19)
Выполнив линеаризацию системы (5) около состояния (0, ?0 , I0 ), по аналогии с системой (8)
получим систему
0
1
?
?M 2
?1
0
0
?(? ? 1)
0
0
?
~
~
?B
?B
E
= A?
,
?t
?x
?
A =?
?
?
1
?
?M 2 ?
?
(20)
?,
?
где в силу близости физических процессов матрицы A? и A имеют совпадающие элементы.
Матрица A? является невырожденной, имеет три различных действительных собственных значения
?1 = 0,
?2 = 1/M,
?3 = ?1/M,
и подобна диагональной жордановой матрице J ? c матрицей подобия S ? , J ? = S ? ?1 A? S ? ,
где
?
0
?
S? = ? 1
?
?1/M
1
?
1/M ?
,
1 ?
?
?
S ? ?1
?1 (? ? 1) (? ? 1)
?
(? ? 1)/? ?1/? ?
?
= ? ?M/2
.
1/2?
1/2? ?
?
?
M/2
1/2?
1/2?
0
Если возмущения рассматриваются при относительном движении, то инварианты принимают вид
(? ? 1)
1
?(t, x) ? I(t, x) = C,
?
?
M
1
1
1
? u(t, x) + ?(t, x) + I(t, x) = f2 x +
t ,
2
2?
2?
M
M
1
1
1
u(t, x) + ?(t, x) + I(t, x) = f3 x ?
t ,
2
2?
2?
M
где C | связанная с интегрированием произвольная постоянная. Скорость, плотность и
энтальпия представляют собой линейную комбинацию бегущих волн, распространяющихся
Наука и Образование. МГТУ им. Н.Э. Баумана
394
с характеристическими скоростями 1/M и ?1/M :
?
1
1
1
1 ?
u(t, x) = ? f2 x +
t + f3 x ?
t ,
?
M
M
M
M
?
1
1
?(t, x) = C + f
x+
x?
t +f
t ,
2
3
?
M
M
?
1
1
? I(t, x) = ?C + (? ? 1)f2 x +
t + (? ? 1)f3 x ?
t .
M
(21)
M
Изменение (19) в условиях на компоненты возмущенного газа (6) привело к следующему
соотношению между матрицами, задающими результирующую линейную систему:
A? = A + E,
(22)
где E | единичная матрица. Схожие структура и свойства матриц A и A? , отвечающих
системам (8) и (20), приводят к схожести решений этих систем (18) и (21), построенных при
помощи элементарных бегущих волн. Возникает вопрос: можно ли отразить данный факт в
решении, не обращаясь к элементарным волнам и, следовательно, к отысканию собственных
значений и собственных векторов каждой из матриц A и A? в отдельности (см. раздел 6).
5. n-мерные бегущие волны
При изложении метода решения системы (8), отвечающей распространению малых возмущений в идеальном газе при одномерном течении, нигде не использовались размерность
искомой функции и порядок матрицы. Следовательно, решение произвольной гиперболи~ x) | n-мерная вектор-функция, A | числовая матрица
ческой системы вида (8), где B(t,
порядка n, может быть представлено как линейная комбинация n элементарных бегущих
волн.
Наряду с элементарной волной
x + ?i t,
(23)
i = 1, n,
(здесь ?i | собственные значения матрицы A), которая отвечает дифференциальному оператору
?
?
? ?i
(24)
?t
?x
(возникающему в системе независимых уравнений, содержащих одну функцию), рассмо??i ?
трим выражения
Ex + J t
и
Ex + A t,
(25)
которые отвечают дифференциальным операторам
?J ? E
?
?
?J
?t
?x
и
Наука и Образование. МГТУ им. Н.Э. Баумана
?A ? E
?
?
?A .
?t
?x
(26)
395
Размерность функций, на которые действуют операторы, задана размерностью коэффициента при производных: оператор (24) действует на скалярную функцию Vi , операторы (26)
~ соответственно.
действуют на n-мерные вектор-функции V~ и B
По аналогии с бегущими волнами (23), выражения (25) рассмотрены в работах [7, 8, 9] как
n-мерные бегущие волны (или бегущие волны с матричными коэффициентами), отвечающие
системам дифференциальных уравнений в частных производных 1-го порядка (8) и (10). К
волнам (23) в некоторых случаях будем добавлять определение "скалярные". В этом смысле
~ i будем применять определение <скалярный>, а выражение V~ ,
к инвариантам Римана (S ?1 B)
~ рассматривать как n-мерный инвариант.
V~ = S ?1 B,
По построению волны (25) связаны посредством матрицы S. Имеет место следующее
утверждение.
Утверждение 1. Если произвольная числовая матрица A порядка n подобна матрице J
с матрицей подобия S, то n-мерная бегущая волна (Ex + At), где E | единичная матрица
порядка n, подобна бегущей волне (Ex + Jt) с матрицей подобия S ?1 :
J = S ?1 AS
=?
(Ex + At) = S(Ex + Jt)S ?1 .
(27)
Если матрица J является диагональной жордановой формой матрицы A, то из соотношения (27) следует представление n-мерной бегущей волны через скалярные волны (x + ?i t),
построеннные при помощи собственных значений матрицы A, i = 1, n, [10].
Отметим, что модель одномерного движения идеального газа является примером модели,
в которой имеет место полная интерпретация бегущей волны с матричным коэффициентом
посредством скалярных волн. Данный факт отражает гиперболичность системы уравнений (8) (системы (7)).
6. Решение нормальной системы линейных однородных уравнений
в частных производных 1-го порядка, заданной числовой матрицей
Пусть в (12) функции fi , fi ? C 1 , заданы разложением в ряд
fi (x + ?i t) =
?
X
(x + ?i t)? ??i ,
i = 1, 3,
(28)
?=0
по своим волнам-аргументам, где ??i | числовые коэффициенты. Это означает, что компоненты искомой вектор-функции
т
V~ = (f1 (x + ?1 t), f2 (x + ?2 t), f3 (x + ?3 t))
представлены как разложение по соответствующим скалярным бегущим волнам. Запишем
данный факт посредством n-мерной волны, отвечающей диагональной матрице J:
V~ ? V~ (t, x) =
?
X
(Ex + J t)? ?~? ,
т
?~? = (??1 , ??2 , ??3 ) .
(29)
?=0
Наука и Образование. МГТУ им. Н.Э. Баумана
396
Если выражения (28) и (29) являются решениями уравнений (11) и (10), то в качестве решения
системы (8) естественно рассмотреть выражение
~ ? B(t,
~ x) =
B
?
X
(Ex + A t)?~?? ,
т
(30)
~?? = (??1 , ??2 , ??3 ) ,
?=0
где ~?? | числовой вектор-коэффициент. По построению каждый член ряда (28), (29) и (30),
а также любая частичная сумма этих рядов являются решением соответствующих уравнений [7].
Особенностью построенных решений (28){(30) является то, что однозначно задающий
уравнение или систему коэффициент соответствующей размерности при производной по
пространственной переменной x входит в решение в качестве коэффициента при независимой переменной t.
По отношению к системам (8) и (20) это, в частности, означает, что что решение одномерной задачи о распространении малых возмущений в идеальном газе при относительном
движении может быть получено (с точностью до ~?? ) из решения задачи при абсолютном
движении при замене в (30) матрицы A на матрицу A? , A? = A + E.
Построение решения (30) не требует нахождения собственных чисел и собственных векторов матрицы A. Вместе с тем, их можно ввести в (30) посредством отношения подобия (27)
между трехмерными волнами (Ex + At) и (Ex + Jt):
~ x) =
B(t,
?
X
S(Ex + J t)? S ?1~?? ,
J = S ?1 AS.
(31)
?=0
Поскольку выражения (30) и (31) могут быть поставлены в соответствие системе (8)
(системе (2)) для произвольной числовой матрицы A, условие гиперболичности, которое для
рассматриваемой физической модели является следствием законов сохранения, не является
существенным при любом его определении.
Поскольку при построении выражений (30) и (31) не используется размерность n функ~ x), каждое из них является общим решением системы (2) и частным решением
ции B(t,
системы (1) для произвольного n.
Будем говорить о волновом подходе к линейной системе уравнений в частных производных 1-го порядка (8), которая рассматривается в качестве частного случая системы (1),
понимая под этим разложение решения по n-мерным бегущим волнам.
7. Представление решения задачи Коши
Напомним [3], что система
E
~ x)
~ x)
? B(t,
? B(t,
=A
,
?t
?x
~ x) = (B1 (t, x), . . . , Bn (t, x)),
B(t,
x ? R1 ,
(32)
~ x) | n-мерная вектор-функция, A | произвольная числовая матрица с элеменгде B(t,
тами {aij }ni,j=1 , является системой Ковалевской. Для такой системы естественным образом
возникает задача Коши.
Наука и Образование. МГТУ им. Н.Э. Баумана
397
Пусть n-мерная вектор-функция ?
~ (x) является аналитической функцией x в некоторой
области G0 плоскости t = 0, содержащей точку x0 , x0 = 0, и представлена в G0 рядом
Тейлора
?
X
?
~ (x) =
?
~c? · x ,
?=0
1 ??
~c? =
?
~ (x)
.
?! ?x?
x=0
(33)
Рассмотрим начальное условие по переменной t:
~ x)
B(t,
x=0
т
=?
~ (x),
(34)
?
~ (x) = (?1 (x), . . . , ?n (x)) .
Согласно теореме Коши | Ковалевской задача Коши (32){(34) имеет в области G, G0 ?
G, единственное аналитическое решение, которое выражается степенным рядом (например, [3]):
?
X
Bi (t, x) =
(35)
??i 1 , ?2 x?1 t?2 .
?1 ,?2 =0
Числовые коэффициенты ??i 1 ,?2 однозначно определяются по начальным данным из системы
уравнений (32) в результате реккурентной процедуры по индексу ?2 и являются полиномами
элементов матрицы A и коэффициентов ??i 1 ,l , l ? ?2 ? 1.
Рассмотрим понятие бегущей волны в применении к решению задачи Коши. Начнем со
скалярного случая. Пусть n = 1. Тогда, поскольку уравнение
?B(t, x)
?B(t, x)
=a
,
?t
?x
a11 = a,
B1 (t, x) = B(t, x),
(36)
имеет решение B = f (x + at), где f ? C 1 , для решения задачи Коши в качестве функции f
следует выбрать функцию, которая задает начальное условие [4]:
B(t, x) = ?(x + at),
(37)
?(x) = ?1 (x).
Очевидно (например, [2, 6]), что для гиперболической системы (32) вид функций fi , i = 1, n,
которые определяют сохраняющуюся форму i-й бегущей волны, устанавливается при помощи матрицы подобия:
fi (0, x) = (S ?1 ?
~ (x))i ,
S ?1 AS = J,
J = diag(?1 , . . . , ?n ).
(38)
Тогда при замене в i-й компоненте вектора ?
~ (x) аргумента x на волну (x + ?i t) из (9)
получим следующее решение решение задачи Коши для системы (8), соответствующей
распространению малых возмущений в идеальном газе:
u(t, x) =
3
X
i=1
s1i ?i (x + ?i t),
?(t, x) =
3
X
s2i ?i (x + ?i t),
i=1
I(t, x) =
3
X
s3i ?i (x + ?i t),
(39)
i=1
где {sij }3i,j=1 | элементы матрицы S, столбцами которой являются собственные векторы
матричного коэффициента этой системы. Решение (39) указывает на то, каким образом
в точку (t, x) вдоль проходящих через нее трех характеристик (13) приносятся значения
соответствующих функций ?i (x), i = 1, 3.
Наука и Образование. МГТУ им. Н.Э. Баумана
398
Например, выберем начальные условия в виде
т
?
~ (x) = (? sin x, 0, µ cos x) ,
?, µ 1,
(40)
?? < x < +?,
где ? и µ | малые параметры, задав начальные возмущения скорости и температуры (посредством энтальпии) при одномерном движении идеального газа. Из (39) получим следующее
представление решения задачи Коши:
?
(1 ? M )
µ
(1 ? M )
?
?
t ?
cos x +
t +
u(t, x) = sin x +
?
2
M
2?M
M
?
?
(1 + M )
µ
(1 + M )
?
+ sin x ?
t +
cos x ?
t ,
2
M
2?M
M
?
?
? ?(t, x) = ? µ cos(x ? t) ? ?M sin x + (1 ? M ) t +
?
2
M
?
?
µ
(1 ? M )
?M
(1 + M )
µ
(1 + M )
+
cos
x
+
t
+
sin
x
?
t
+
cos
x
?
t
,
?
2?
M
2
M
2?
M
?
µ
?(? ? 1)M
(1 ? M )
?
I(t, x) = cos(x ? t) ?
sin x +
t +
?
?
2
M
?
µ(? ? 1)
(1 ? M )
?(? ? 1)M
(1 + M )
?
+
cos x +
t +
sin x ?
t +
2?
M
2
M
?
?
µ(? ? 1)
(1 + M )
?
+
cos x ?
t .
2?
(41)
M
В (41), в частности, представлено, каким образом (с какой амплитудой и с какой частотой) начальное возмущение скорости газа ? sin x скажется на возмущении его плотности и
температуры (посредством энтальпии).
Воспользуемся теперь разложением начальной аналитической функции в ряд Тейлора (33). Пусть функции fi , сохраняющие свои значения на волнах-аргументах (x + ?i t),
представляют собой разложение (28). Тогда при переходе к матричной записи получим
разложение n-мерного инварианта V~ в ряд (29) по n-мерным бегущим волнам. Из (9)
следует начальное условие на функцию V~ :
V~ (t, x)
t=0
=
?
X
(Ex +
?=0
Jt)? ?~?
t=0
=S
?1
?
X
~c? · x? .
?=0
Следовательно, имеет место соотношение
~т
?1
?? = S ~c? ,
1 ??
?
~ (x)
,
~c? =
?
?! ?x
x=0
(42)
между коэффициентами разложения по n-мерным бегущим волнам векторного инварианта V~
и коэффициентами разложения в ряд Тейлора начальной вектор-функции ?
~ (x).
В то же время, рассмотрев при t = 0 общее решение системы (32) в виде (30), приходим
к выводу, что для получения решения задачи Коши (32){(33) в качестве коэффициентов ~??
Наука и Образование. МГТУ им. Н.Э. Баумана
399
при (Ex + At)? следует выбрать коэффициенты ряда (33) при x? . Тогда решение задачи
Коши принимает вид [8]
~ x) =
B(t,
?
X
(Ex + At)?
?=0
1 ??
?
~ (x) .
?! ?x?
x=0
(43)
Решение (43), которое записано без дополнительных вычислений, содержит матричный коэффициент системы в явном виде, а начальные данные | в виде коэффициентов разложения
аналитической начальной функции в ряд Тейлора. Ряд (43) организован по волновому принципу, отличному от принципа организации ряда (35).
Например, чтобы получить решение рассматриваемой задачи Коши для уравнений распространения малых возмущений в идеальном газе при одномерном движении, следует
подставить в (43) матричный коэффициент A системы (8) и векторные коэффициенты
т
~?0 = (0, 0, µ) ,
µ т
~?2 = 0, 0, ?
,
2!
т
~?1 = (?, 0, 0) ,
т
?
~?3 = ? , 0, 0 ,
3!
...,
которые отвечают разложению функции (40) в ряд Тейлора в окрестности нуля. Сравнив ком~ x), B(t,
~ x) = (u(t, x), ?(t, x), I(t, x))т , с решением (41) при
поненты полученного вектора B(t,
разложении тригонометрических функций в ряд Тейлора, убеждаемся в их тождественности.
Напомним при этом, что любая частичная сумма ряда (43) с указанными коэффициентами A
и ~?? является решением системы (8).
Запишем систему уравнений (32) в операторном виде:
E
~ x)
? B(t,
~ x),
= A0 B(t,
?t
A0 = A
?
.
?x
(44)
По аналогии с обыкновенным дифференциальным уравнением представим решение задачи
Коши для системы (44) через операторную экспоненту:
~ x) = exp(A0 t) ?
B(t,
~ (x).
(45)
Под результатом действия оператора exp(A0 t) на аналитическую функцию ?
~ (x) понимают
функцию, к которой сходится степенной ряд
exp(A0 t) ?
~ (x) =
? i
X
t
i=0
i!
(A0 )i ?
~ (x),
exp(A0 t)
t=0
(46)
= E.
Приравняв формальное решение задачи Коши в виде операторной экспоненты (45) и найденное решение (43), получим формулу
?
?
X
1 ??
? ?1 1 ?
exp(A t)~
?(x) ?
(Ex + At)
?
~
(x)
=
S(Ex
+
Jt)
S
?
~ (x) , (47)
?! ?x?
?! ?x?
x=0
?=0
?=0
x=0
0
?
X
?
где
A0 = A
?
,
?x
Наука и Образование. МГТУ им. Н.Э. Баумана
J = S ?1 AS.
400
Таким образом, в рамках модели волнового взаимодействия решение задачи Коши представлено как разложение в ряд (сходящийся в некоторой окрестности начальной точки) по
соответствующей задаче n-мерной бегущей волне. Векторные коэффициенты ряда и матричный коэффициент волны известны. Общий член ряда может быть представлен через
собственные значения и собственные векторы матричного коэффициента волны.
Отметим также, что задача Коши для уравнения (44) позволяет трактовать n-мерные
бегущие волны как результат действия соответствующей этому уравнению операторной
экспоненты exp(A0 t) на n-мерную аналитическую функцию ?
~ (x).
Разложение решения задачи Коши (32){(34) в ряд (47) имеет место как для гиперболической системы (матрица A подобна диагональной жордановой матрице J), так и для
системы (32), заданной произвольной матрицей (A подобна жордановой матрице J).
8. Заключение
В рамках предложенной авторами модели волнового взаимодействия в качестве тестовой
задачи сформулирована и решена задача о распространении малых возмущений в идеальном
газе при одномерном движении.
С использованием традиционного подхода к решению гиперболической системы малые
возмущения скорости, плотности и энтальпии газа единообразно представлены как линейная
комбинация произвольных непрерывно-дифференцируемых функций, аргументами которых
являются элементарные бегущие волны.
Число элементарных волн и число инвариантов Римана, отвечающих системе, равно
трем | порядку задающей систему матрицы, что рассматривается как одно из определений гиперболичности. Три пересекающиеся характеристики образуют на плоскости (t, x)
характеристический треугольник. Это означает, что в идеальном газе при малых возмущениях скорости, плотности и энтальпии (температуры) возникает движущаяся замкнутая
структура, которая однозначно задана параметрами газа (числом Маха).
Представление о бегущей волне как решении гиперболической системы, которое обладает свойством инвариантности, перенесено на n-мерный (матричный) аналог и связано с
решением нормальной системы уравнений в частных производных 1-го порядка, которая
задана произвольной числовой матрицей. Указано представление n-мерной волны через
элементарные, которые рассмотрены как скалярные волны (n = 1).
Решение задачи Коши в окрестности точки x0 , x0 = 0, представлено как разложение
в ряд по n-мерной бегущей волне с коэффициентами, равными n-мерным векторным коэффициентам разложения аналитической начальной функции в ряд Тейлора. Структура
построенного ряда отлична от структуры ряда, фигурирующего в теореме КошиЁ| Ковалевской. Установлено, что n-мерные бегущие волны, отвечающие системе (1) при s = 1
(системе (2)), порождаются действием отвечающей этой системе операторной экспоненты
на аналитическую функцию.
Наука и Образование. МГТУ им. Н.Э. Баумана
401
В продолжение исследований в качестве физической модели для системы (1) при s > 1
предлагается рассмотреть задачу о распространении малых возмущений в идеальном газе
при двумерном (плоскопараллельном) движении и движении в пространстве. Это потребует
применения соответствующего оператора волнового взаимодействия [7, 8]. Взаимодействие
осуществляется между волнами (скалярными или n-мерными), 3 < n ? 5, построенными
для каждой матрицы-коэффициента Ak , s = 1, 3, системы (1). Полученное решение (30),
связанное с волнами, которые распространяются вдоль выбранной оси OXk , k = 1, 3, является частным решением задачи о распространении малых возмущений в идеальном газе при
движении в пространстве и составной частью ее общего решения.
Список литературы
1. Годунов C.К. Уравнения математической физики. М.: Наука, 1979. 392 c.
2. Куликовский А.Г., Свешникова Е.И., Чугайнова А.П. Лекционные курсы НОЦ. Вып. 16.
Математические методы изучения разрывных решений нелинейных гиперболических
систем уравнений / Математический институт им. В.А. Стеклова РАН (МИАН). М.:
МИАН, 2010. 122 с.
3. Курант Р. Уравнения с частными производными: пер. c англ. / под ред. О.А. Олейник. М.:
Мир, 1964. 830 с. [Courant R. Methods of mathematical physics. Vol. 2. Partial Differential
Equations. New York-London: Interscience Publ., 1962].
4. Бхатнагар П. Нелинейные волны в одномерных дисперсных системах: пер. с англ. /
под ред. П.Е. Краснушкина, Н.Р. Сибгатуллина. М.: Мир, 1983. 136 c. [Bhatnagar P. L.
Nonlinear Waves in One-dimensional Dispersive Systems. Oxford: Clarendon Press, 1979].
5. Лойцянский Л. Г. Механика жидкости и газа. М.: Наука, 1987. 840 с.
6. Toro E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer Berlin
Heidelberg, 2009. 724 p. DOI: 10.1007/b79761
7. Феоктистов В.В., Мякинник О.О. Структура ряда для решения системы уравнений с
частными производными 1-го порядка // Вестник МГТУ им. H.Э. Баумана. Сер. Естественные науки. 2009. № 4. С. 3{22.
8. Феоктистов В.В., Мякинник О.О. Оператор волнового взаимодействия и нормальная
форма системы линейных уравнений в частных производных 1-го порядка // Воронежская весенняя математическая школа <Современные методы теории краевых задач. Понтрягинские чтения | XXI>: матер. Воронеж: Изд-во ВГУ, 2010. С. 232{233.
9. Феоктистов В.В., Мякинник О.О. Бегущие волны с матричными коэффициентами //
Воронежская зимняя математическая школа <Современные методы теории функций и
смежные проблемы>: матер. Воронеж: Изд-во ВГУ, 2013. С. 255{256.
Наука и Образование. МГТУ им. Н.Э. Баумана
402
10. Феоктистов В.В., Мякинник О.О. Представление n-мерных бегущих волн через скалярные // Воронежская весенняя математическая школа <Современные методы теории
краевых задач. Понтрягинские чтения | XXV>: матер. Воронеж: Изд-во ВГУ, 2014.
С. 178{180.
11. Geroch R. Partial Differential Equations of Physics // In: General Relativity: Proceedings of
the 46th Scottish Universities Summer School in Physics (Aberdeen, July 1995). Edinburgh:
SUSSP Publ., 1996. P. 19{60.
12. Ковеня В.М., Яненко Н.Н. Метод расщепления в задачах газовой динамики. Новосибирск: Наука, 1981. 304 c.
13. Соболев С. Л. Уравнения математической физики. М.: Наука, 1992. 432 c.
14. Капцов О.В. Инварианты характеристик систем уравнений с частными производными //
Сибирский математический журнал. 2004. Т. 45, № 3. C. 577{591.
Наука и Образование. МГТУ им. Н.Э. Баумана
403
Science and Education of the Bauman MSTU,
2014, no. 12, pp. 372{391.
DOI: 10.7463/1214.0744956
Received:
Revised:
28.06.2014
02.12.2014
c Bauman Moscow State Technical University
Ideal Gas Small Perturbations:
Waves, Invariants and Cauchy Problem
Feoktistov V. V., Miakinnik O. O.
*
Bauman Moscow State Technical University
Keywords: hyperbolic system, traveling waves, Riemann invariants, n-dimensional traveling waves,
expansion in terms of traveling waves, Cauchy problem solution, ideal gas, small perturbations
The wave interaction model the authors have developed for the first-order linear homogenous
Cauchy system, which is specified by numerical matrices at the spatial derivatives, is considered
as applied to the one-dimensional ideal gas flow. The linear waves, meeting model requirements,
are separated from quasi-linear process by the method of small perturbations. The resulting
linear hyperbolic system is specified by the matrix containing gas parameters and written for
the unknown vector, components of which are small perturbations of gas velocity, density, and
enthalpy.
Following a common practice for a hyperbolic system, the method for solving the system via
traveling waves, propagated axially with characteristic speeds, is analyzed for the physical process
involved. Three Riemann invariants are found. The general solution is represented as a linear
combination of three arbitrary functions of wave-arguments, specified by characteristic speeds.
The initial conditions are imposed on the three functions, which can be measured experimentally,
to determine the unchangeable form of traveling waves. The formation of characteristic triangle,
uniquely determined by gas parameters, is connected with the existence of some closed structure
moving in the ideal gas.
To hold the matrix form, a 3-by-3 diagonal matrix is formatted from the traveling waves as its
entries. This matrix is considered as a three-dimensional wave, specified by coefficient matrix with
3 characteristic speeds along its diagonal. Using a similarity transformation, the diagonal matrix
coefficient is replaced by the matrix specifying the system to be solved.
By analogy with hyperbolic case, an n-dimensional wave, specified by an arbitrary numerical nby-n matrix, is considered as a solution of a first-order Cauchy system specified by this matrix. The
Cauchy problem solution is represented as a series expansion in terms of n-dimensional wave with
coefficients equal to the coefficients of n-dimensional initial analytic vector Taylor's expansion.
Science and Education of the Bauman MSTU
404
It was found that an n-dimensional wave is generated by the action of differential operator
exponential on an n-dimensional analytic initial vector via Cauchy problem for the related system.
It is proposed to apply the model to the two-dimensional gas flow using the notion of interaction
between the waves propagated along each of two coordinate axes.
References
1. Godunov S.K. Uravneniya matematicheskoy fiziki [Equations of Mathematical Physics].
Moscow, Nauka Publ., 1979. 392 p. (in Russian).
2. Kulikovskiy
Документ
Категория
Без категории
Просмотров
8
Размер файла
616 Кб
Теги
кошик, волна, возмущений, малыш, инвариантов, идеального, задачи, газа
1/--страниц
Пожаловаться на содержимое документа