Вычислительные технологии
Том 17, № 3, 2012
К вопросу сохранения типа системы квазилинейных уравнений в частных производных
В. В. Блшуров
Центр экологического и техногенного мониторинга, г. Трехгорный, Россия
e-mail: VVBashurov@mail.ru
Рассматривается система квазилинейных дифференциальных уравнений в частных производных, записанная в подвижной системе координат. Связь между подвижной и неподвижной системами координат осуществляется с помощью квазилинейного уравнения для вектора переносной скорости. Показано, что в общем случае тип системы зависит от вида дополнительного уравнения для вектора переносной скорости. Доказаны теоремы, дающие достаточные условия сохранения типа исходной системы дифференциальных уравнений, записанной в декартовой системе координат.
Ключевые слова: система квазилинейных уравнений в частных производных, тип дифференциального уравнения, инвариантность типа.
Введение
В связи с развитием конечноразностных методов решения дифференциальных уравнений большое внимание уделяется построению так называемых адаптирующихся сеток. Принципы, на которых основано конструирование алгоритма построения таких сеток, достаточно подробно изложены в работах [1, 2]. Все алгоритмы условно разделим на три класса: функциональные, эллиптические и эволюционные. Для первого класса характерно задание узлов сетки при помощи формулы или интерполяции между некоторым набором линий (границы области и границы раздела сред, ударные волны, слабые разрывы и т. д.) К таким методам можно отнести методы, описанные в [1, 3]. Второй класс охватывает такие методы построения сеток, в которых строится отображение реального физического пространства на параметрическое, удовлетворяющее некоторой краевой задаче, как правило, эллиптического типа (например, реализующее конформное отображение или иной геометрический принцип) [4-7]. К третьему классу отнесем методы, использующие для нахождения узлов сетки дифференциальное уравнение для вектора переносной скорости, т. е. для скорости перемещения самого узла в неподвижной координатной системе. Как правило, такое уравнение (уравнения), добавленное (добавленные) к основной системе уравнений, образует (образуют) систему уравнений в частных производных, названную нами расширенной, для которой ставится задача Коши. В процессе решения этой системы вместе с "нужными" функциями находятся узлы сетки, или отображение неподвижной системы координат на параметрическое пространство переменных, в которых записывается исходная система уравнений. Это отображение зависит от времени и никаким способом, кроме решения расширенной системы, его получить невозможно [8-11].
Общим для всех трёх типов алгоритмов является вопрос о невырожденности преобразования реального (декартова) пространства в параметрическое, задающее, например, узлы сетки. Особый интерес представляет и вопрос: не меняется ли тип исходного уравнения (уравнений), если к нему (к ним) добавляются уравнения, задающие параметризацию. Впервые, насколько нам известно, это было рассмотрено в работе [8]: если исходная система дифференциальных уравнений в частных производных, записанная в декартовых координатах, относится к гиперболическому типу и поэтому классическая задача Коши имеет решение хотя бы в некоторой окрестности многообразия, на котором заданы начальные условия. то та же система, записанная в системе координат, задаваемой решением расширенной системы, может изменить свой первоначальный тип.
На первый взгляд, в постановке такого вопроса есть что-то алогичное. Ведь до сих пор утверждалось, что невырожденное преобразование как зависимых, так и независимых переменных не меняет тип исходной системы. Вместе с тем доказывалось это утверждение исключительно для невырожденных преобразований, задаваемых функционально. Однако, как следует из [12], в общем случае при преобразованиях системы координат, задаваемых дифференциальным уравнением (уравнениями) для вектора переносной скорости, тип системы не инвариантен.
В настоящей работе получены достаточные условия, накладываемые на преобразования координат (на коэффициенты уравнений, связывающих координаты, в которых изначально были написаны уравнения, и на координаты, в которых будут решаться эти уравнения), для того чтобы исходный тип уравнений не менялся. В качестве исходной системы выбраны уравнение переноса, уравнения газовой динамики и уравнения ги-поупругости [8]. Исследован случай, когда искомые функции зависят только от двух переменных — времени t и пространственной переменной £.
1. Некоторые сведения из тензорного анализа
Введем в рассмотрение ковариантный базис {Эк} и метрический тензор дц~ в простран-
d r
стве {£к}. Базисные векторы {Эк} определены равенством Эк = г.
д£к
Чтобы сократить размер формул, воспользуемся предложенным Эйнштейном способом изображать суммирование без греческого символа Е, а именно, если в формуле у каких-либо сомножителей одинаковые индексы располагаются на "разных" этажах, то это подразумевает суммирование по данному индексу.
Вектор скорости u = ukЭк можно представить в виде суммы двух векторов — вектора относительной и вектора переносной скорости. Обозначим их соответственно u(r) = икг)Эк и u(e) = «ке)Эк. Вектор переносной скорости определяется уравнением
d r
dt
= U(e).
d
Оператор — определим следующим образом: dt
d д д _=__i «к _
dt dt + "(r) д£к.
Уравнение переноса может быть записано как
du
~dt = '
откуда
- u , -u
dU +uk') -u = 0 (1)
В неподвижной эйлеровой системе координат вектор переносной скорости равен нулю, т. е.
U(e) = 0. (2)
В лагранжевой системе координат вектор относительной скорости равен нулю, и поэтому уравнение для переносной скорости принимает вид
u = u(e). (3)
Получим ряд нужных в дальнейшем формул. При выводе уравнения
дЭк д дг д д
dt -Ск дг д£кu(e) д£к Me) Э
ч-е
было использовано равенство
-ue \
-k- + u]me)T<m^ Эе = vku(e) Э(
д Э,
rm э
-Ск
где Гmk — символ Кристоффеля. Используя определение метрического тензора gke (Эк' Э(), получим
дцы I дЭk \ I д Э(
-ж =( ж'Э') + Iя -иг
Здесь скобки означают скалярное произведение векторов. Воспользовавшись уравнением (4), получим
-gke
"д^ = Vkи™) • gme + Veu™) • gmk. (5)
Уравнения (5) по сути позволяют установить взаимно-однозначное соответствие между координатами {Ck} и декартовыми координатами {Xk}.
Группу этих уравнений будем называть уравнениями для метрики. Выполним в (1) дифференцирование по переменным {Ck}. С учётом формулы (4) получим
k , дЭ' -u
-Эк + uk^ + u7l)—— Э' + u7l)u
дuk ,д Э' -uk -.-Э'
э+uk ^+um) ^+
-uk r , -uk
-f + u(Veuke) + [um - u}J -cm + [um - u})] ue Г km ) Э' = 0.
С учётом линейной независимости векторов {Эк} уравнение переноса принимает вид
duk
— + u'V,ufe) + [um - u(] Vmuk = 0. (6)
Формулы (5) и (6) образуют систему, представляющую уравнение переноса, записанное в произвольной системе координат. Система не полна. Для её замыкания необходимо добавить уравнение для вектора переносной скорости.
Для первого типа методов, названных выше функциональными, уравнение для переносной скорости задаётся формулой. Так, эйлерово пространство, в котором введён неподвижный базис, описывается простым уравнением (1). В этом случае уравнение переноса принимает вид
duk , , „
— + umVmuk = 0. (7)
Если пространство параметризуется лагранжевыми координатами, т. е. в качестве координат выступают координаты, которые подвижная точка занимала в начальный момент времени, то уравнение для вектора переносной скорости записывается в виде (3) и уравнение переноса принимает тот же вид (7).
В общем случае систему координат, в которых записывается уравнение переноса, будем называть счётной, поскольку именно к записанным в таком виде уравнениям применяют те или другие методы численного решения.
2. Нелинейное уравнение переноса
Рассмотрим проблему сохранения типа дифференциального уравнения (системы уравнений) на примере простейшего нелинейного уравнения переноса. Напомним, что в неподвижной декартовой системе координат оно имеет вид
du du dt dx
Здесь u — скорость. В подвижной счётной системе координат, задаваемой вектором переносной скорости v, это уравнение переходит в систему двух нелинейных уравнений в частных производных (5), (6):
du du dv u2g11 dgn — + (u - v)— + u— + —--— = 0,
dt дС dC 2 dC (8)
dgii dgii 0 dv .
"dt" - v^i" - 2g11 ai = 0.
Здесь g11 и g11 — ковариантная и контрвариантная компоненты метрического тензора, i — координата в подвижной системе, v — переносная скорость. Отметим, что переносная скорость и скорость — суть контрвариантные вектора.
Пусть уравнение, задающее переносную скорость v, имеет вид
dv du dg11 dv dt- adi - - edi = f, (9)
где коэффициенты a, b, e — произвольные функции переменных u, v, g11.
Система уравнений (8)-(9) [12] является системой, эквивалентной одному уравнению переноса в декартовой неподвижной системе координат, и поэтому её по-прежнему будем называть уравнением переноса. Убедимся, что эта система в декартовой системе координат гиперболична.
Действительно, при v = 0 (8)-(9) запишется как
ii
0,
du du u2g11 дд
--\- u--\---
dt ^ д£ 2 д£
dgii
ж=0
и характеристическое уравнение имеет только действительные корни. Ставя для этой системы задачу Коши, т. е. задавая в начальный момент времени метрику д11, переносную скорость v и "истинную" скорость u, можно установить в окрестности начальных значений соответствие между декартовой неподвижной системой координат и "счётной" системой.
Введем вектор решения w = (u, g11, v). Систему (8) виде
(9) можно записать в матричном
dw dw
~at +
F.
:iq)
Здесь
A
V
u - v 0
—a
u2g11
u
2
-v -2gn
b e
F
/
0 0
f
Рассмотрим теорему.
Теорема 1. Если вектор w на прямой t = 0 аналитичен, gn(0,£) > £ > 0, то соответствие координатных систем {х} и {£}, определяемое решением задачи Коши для системы (10), не вырождено в некоторой окрестности прямой t = 0.
Доказательство следует непосредственно из того, что (10) является системой типа системы Коши — Ковалевской; следовательно, вектор решения w непрерывен в некоторой окрестности прямой t = 0, а так как д11(0, £) > 0 и якобиан преобразования х = х(£) удовлетворяет уравнению
дх
= л/эй,
то однозначность очевидна.
Поставим следующий вопрос: если исходное уравнение переноса гиперболично, что при определённых условиях на гладкость начальных данных зквивалентно разрешимости задачи Коши, то сохраняется ли тип уравнения при переходе к подвижной системе координат, задаваемой квазилинейным уравнением для переносной скорости v? Как показывает приводимый ниже пример, — не сохраняется.
Характеристическое уравнение для системы (10) имеет вид
[Л - (u - v)] [(Л + v)(A + e) + au - 2g11b] = 0. Положим a = u, e = v, b = f = 0. Тогда характеристическое уравнение запишется как
[Л - (u - v)] [(Л + v)2 + u2] =0.
Оно имеет пару комплексных корней. Уравнение переноса в системе координат, задаваемой уравнением
дv du dv
т — иЖ — vdc= 0'
является уравнением смешанного типа.
Что гарантирует сохранение типа уравнения? На этот вопрос дает ответ следующая теорема.
Теорема 2. Если уравнение (9) инвариантно относительно преобразования п = дп
п(0> = 0, то система (10) является системой гиперболического типа.
Доказательство. Рассмотрим уравнение (9) для переносной скорости в системе координат {£} и введем некоторую новую систему п = п(С). Тогда новая переносная
д£
скорость v будет связана со старой равенством v = —v, так как v и v — контрвариант-
дп
д£
ные составляющие вектора переносной скорости. Обозначим — через p. Неизвестные
дп
функции в новой системе координат определяются по формулам
u = pu gii = Agn, д11 = <p2gn. p2
Подставим новые значения искомых функций в уравнение (9). Отметим, что это уравнение носит тензорный характер и поэтому инвариантно относительно преобразования п = п(С). Выпишем уравнение (9) в новых координатах
dv v du b дд11 v dv
dt рдп p4 дп рдп
dp drq
au
2 bg
11
vv
p4
f.
Поскольку по условию теоремы уравнение (9) инвариантно относительно преобразования координат, то из предыдущего равенства следует:
1) au 2bW1 +
1) au--5—+ vv
p3
0;
2) коэффициенты a, e, f являются контрвариантными тензорами первой валентности, a, b — четвертой валентности. Условия 1, 2 по сути являются условиями теоремы 2. Поскольку p произвольно, положим p =1. Тогда условие 1 принимает вид
au — 2bg11 + ev = 0.
11)
Рассмотрим характеристическое уравнение для системы
Л — (u — v) —
0
a
u2 g11 2
Л + v b
—u
2gn
Л + e
0
Раскрывая определитель, получим
[Л — (u — v)][^ + v)^ + e) + au — 2g11b] = 0.
Отсюда вторая квадратная скобка дает уравнение для Л
Л2 + Л(и + e) + ve + au — 2g11b = 0.
В силу (11) оно имеет действительные корни.
Теорема доказана. □
3. Уравнения движения гипоупругого тела
Рассмотрим уравнения "одномерного" движения гипоупругого тела [13], записанные в произвольной системе координат, задаваемой переносной скоростью v с уравнением состояния в виде p = c2(p — ро):
du du dv c2 dp 1 dT11 1 dg11 0 ,,
m +(u - v) * + ^ = - 7g dp +1 "aT + 2Pdgr 11 - Pu2)g11+
(T 11g22 - T22g11) + (t 11g33 - T33g11) ,
dp du dp pu / 11 dg1i 22 dg22 33 dg33\ n
dt+pdi + (u-v)di + T[g ИТ +g ИГ+g =0,
dT11 8T11 du /4 \ n dv dgii n /2 n „ ,
it+<u - v) der = J - 2T11 di+- T"1 -
^ 11 22 dg22 ^ 11 33 dg33
:ug g g
3 * y di 3 * y di'
dgii dv dgn
~dt~ = 2g11 di + ,
dg22 dg22
-= v-,
dt di ,
dg33 dg33
=. (12)
Здесь c — объёмная скорость звука, p, p0 — текущая и начальная плотности, g11, g22,
33
g — контрвариантные составляющие метрического тензора, g11, g22, g33 — ковариант-ные составляющие того же тензора, T11 — контрвариантная составляющая девиатора тензора напряжений.
К системе уравнений (12) необходимо добавить соотношения, определяющие выбор системы координат. Рассмотрим полученную систему для случая v = au, где a — постоянная величина.
При a = 1 имеем лагранжеву систему координат, при a = 0 — неподвижную эйлерову. Тогда система (12) принимает вид
du du c2 „ dp 1 dT11 — + u— =--g11^ + -
dt di P di P di
( %i[2T11 - pu2]g11 + ^ [T 11g22 - T22g11] + д^3 [T11 g33 - T33g11]
2pV д£
д^ д£
др д-u др 1
Ж + +(1 - a)us£ - 2pu
дЭ1
\ g22 ^g^2 \ g33 дд33
(дЭ3
д£
дт11 дт11 дu
"дГ +(1 - a)u^ =
4
^^g11 - 2T 11а
д£
, дЭп 11
+ u
д£
2
a™11 - т 11
дд22 ^ 11 22 дд33 ^ 11 33
ug д - -^т--^д д ,
д£ 3
д£ 3
дЭ11 0 д« дЭ11
"Ж" = 2д11ад£ + au^T дд22 дд22
au
дt дд33
au-
д£ ' дд33
13)
дt д£ '
Система (13) квазилинейна и ее тип определяется собственными числами характери
стического уравнения
(au - Л)2(Л + [1 - a]u)х
T
11
х <-[2(Л - au)a - 2а(Л + [1 - a]u) + 2au] +
р
+cLg {au - Л - au} + [Л + (1 - a)u]^ + Лu - аЛ^ > = 0.
Здесь
3P \ c2
3р
Уравнение для определения характеристических чисел распадается на следующие:
(au - Л)2 = 0,
Отсюда следует
Л + (1 - a)u = 0, Л [Л2 + u2 - (1 - 2a + a2) - c^g11 - 2Лu(1 - a)] = 0.
Л1 = au, Л2 = (a - 1)u, Л3 = 0,
Л4,5 = -(1 - a)u ± c^x/g11.
Таким образом, при выборе переносной скорости пропорционально абсолютной система имеет действительные характеристики.
Поскольку уравнение переноса следует из системы (12) при c = 0 и ^ = 0 и относительно него справедливо утверждение о несохранении типа уравнения при переходе от одной системы координат к другой, в силу непрерывной зависимости корней алгебраического уравнения от коэффициентов неинвариантность типа при переходе к подвижной системе координат очевидна.
0
2
c
L
4. Уравнения газовой динамики
Рассмотрим уравнения газовой динамики, записанные для двух переменных С и t и вытекающие из уравнений (12) при ^ = 0. В этом случае систему удобно представить в виде
а,. а,. „2 а, „.2 я,, + (u — v)--1--q11--1--q11—
dt
du du c2 ^ dp u2 ,dq22 dv
—+(u — v)—+-q11 — + — q11
dt
дГ P дС
l \др •дС (u — v) дС
p
дс + = 0,
др дu др pu ( .. д.ц д.22 33 дqзз
^7 + + (u — v)^r+ — q 1 " ' 33
д.ц д.1^ дv
— — 2q22 дс
д.22 д.22 ■ — v-
2
0,
дС
+ q2
дС
+q
дС
0,
дt д.зз
~St"
v
дС
д.зз
0, 0,
дv дu др
— адё — — c дС
д.ц dv
f.
В матричном виде система запишется как
д^ д^ ж+= F
где матрица A
A
c2 c 11 —q p 0 u2 2р b
u — Л a —u e — Л 0 0 0 0
—p 0 v — u — Л pu 11 — Т q v — Л pu --q22 2 У 0 pu 2"' 0
0 2.11 0
0 0 0 0 v — Л 0
0 0 0 0 0 v
зз
V
Характеристическое уравнение для этой системы имеет вид
det^E — A) = 0.
Раскрывая определитель, получим
(Л — v)2[c2qn — (Л + u — v)2p2 — (v + ^)Л + vd + au — 2qne] = 0.
:14)
:15)
Чтобы гарантировать гиперболичность системы (14), наложим на коэффициенты уравнения, задающего переносную скорость, ограничения, которые формулируются в виде следующей теоремы.
Теорема 3. Если уравнение
дv дu др д.22 дv
:1б)
инвариантно относительно преобразования п = п(0, то система (16) имеет гиперболический тип.
Инвариантность первых пяти уравнений в системе (14) очевидна. Новая скорость v
связана с прежней переносной скоростью v соотношением
д
v = ^-v. дп
д
Обозначим производную — через p. Переменные в новой системе координат определя-
дп
ются по известным формулам
р = р, u = pu, gn = p2gn. Подставляя новые значения в уравнение (15), получим
p
дv aдu b др v дg11 d дv дt ^дп p2 дп p4 дп ^дп
p4
au 2g11e dv
p
p'
4
p
= f.
Поскольку по условию теоремы уравнение (16) инвариантно относительно преобразования п = п( ) , то имеют место два условия:
n— 2g^v^dv п
1) au---—|--= 0;
p3 p
2) коэффициенты a, d , f — контрвариантные составляющие тензора первой валентности, b — второй, v — четвёртой валентности. Эти условия по сути являются математической формой записи условий теоремы.
Поскольку p произвольно, то условие 1 при p =1 принимает вид
au — 2g11e + dv = 0.
:17)
Рассмотрим характеристическое уравнение (15). Действительность характеристик (а, значит, и гиперболичность рассматриваемой системы) связана со знаком дискриминанта D квадратного уравнения
Л2 — (v + d^ + vd + au — 2g11e = 0, D = v2 — 2vd + d2 — 4ua + 8g11e.
Подставляя взятое из (17) выражение для 2g11e, получим
D = (v + d)2 > 0.
Доказательство закончено. □
Таким образом, в работе доказаны теоремы о достаточных условиях, которым должны удовлетворять коэффициенты уравнения для переносной скорости, чтобы тип системы не менялся, для уравнений переноса и газовой динамики в случае, когда зависимые переменные зависят только от двух независимых переменных времени и одной пространственной координаты. Для гипоупругого тела отмечено, что введение подвижной системы координат может изменить тип исходной системы.
Список литературы
[1] Численное решение многомерных задач газовой динамики / Под ред. С.К. Годунова. М.: Наука, 1976.
[2] Сидоров А.Ф., Шабашова Т.И. Об одном методе расчёта оптимальных разностных сеток для многомерных областей // Численные методы механики сплошной среды. 1981. Т. 12, № 5. C. 106-123.
[3] Сараев В.А. Построение регулярных четырёхугольных сеток в областях с криволинейными границами по точкам вложенных друг в друга криволинейных четырёхугольников // Вопросы атомной науки и техники. Методики и программы численного решения задач матем. физики. 1985. Вып. 3. C. 25-38.
[4] Годунов С.К., Прокопов Г.Д. О расчётах конформных отображений и построении разностных сеток // Журн. вычисл. математики и матем. физики. 1967. Т. 7, № 5. C. 1031-1059.
[5] Breakbill J.U., Saltzman J.S. Adaptive zoning for singular problems in two dimensions // J. Comput. Phys. 1982. Vol. 46. P. 342-368.
[6] Крейс Р.И., Тэймз Ф.К., Хасан Х.А. Построение адаптирующихся сеток с помощью вариационного метода Брэкбилла —Занцмана // AIAA J. 1966. No 3. P. 404-410.
[7] Яненко Н.Н., Данаев Н.Т., ЛисЕйкин В.Д. O вариационном методе построения сеток // Численные методы механики сплошной среды. 1977. Т. 8, № 4. C. 157-163.
[8] Медведев А.Е., Фомин В.М. Метод подвижных координат в задачах механики сплошной среды // Прикладная математика и механика. 1981. Т. 45, вып. 4.
[9] Неуважаев В.Е., Фролов В.Д., Яненко Н.Н. Уравнение движения теплопроводного газа в смешанных эйлеро-лагранжевых координатах // Численные методы механики сплошной среды. 1972. Т. 3, № 1. C. 90-96.
[10] Дарьин Н.А., Мажукин В.И. Об одном подходе к построению адаптивных разностных сеток // Докл. АН. 1988. T. 298, № 1. C. 64-68.
[11] Гриднева В.А., Меркулова Н.Н. О построении подвижных сеток // Численные методы механики сплошной среды. 1983. T. 1, № 4. C. 34-44.
[12] Башуров В.В. Инвариантные свойства модельного уравнения переноса при преобразованиях системы координат // Вопросы атомной науки и техники. Математическое моделирование физических процессов. 1989. Вып. 3. С. 24-28.
[13] Башуров В.В., Свидинский В.А. Об уравнениях динамики упругих сред в произвольных системах координат // Численные методы механики сплошной среды. 1982. Т. 13, № 2. C. 14-29.
Поступила в 'редакцию 17 февраля 2012 г.
с доработки — 13 апреля 2012 г.
ЖВТ. 2Q12. Т. 17, № 3. С. 13-23.
On invariance of type of quasi-linear partial differential equations
Bashurov V.V.
We consider a system of linear partial differential equations written in a moving coordinate system. Relation between the moving and fixed coordinate systems is realized by a quasi-linear equation for the moving velocity vector. It is shown that the type of system generally depends on the additional equation for the velocity of the moving vector. A theorem on sufficient conditions to conserve the type of initial system of differential equations written in Cartesian coordinates is proved.
Keywords: system of quasi-linear partial differential equations, differential equation type, type invariance.