Редукция переопределенных систем дифференциальных уравнений математической физики Текст научной статьи по специальности «Математика»

www.volsu.ru

DOI: https://doi.org/10.15688/mpcm.jvolsu.2017.4.5

УДК 530.182, 53.02 ББК 22.19

РЕДУКЦИЯ ПЕРЕОПРЕДЕЛЕННЫХ СИСТЕМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ МАТЕМАТИЧЕСКОЙ ФИЗИКИ

Максим Леонидович Зайцев

Аспирант, соискатель,

Институт проблем безопасного развития атомной энергетики РАН mlzaytsev@gmail. com

ул. Большая Тульская, 52, 115191 г. Москва, Российская Федерация

Вячеслав Борисович Аккерман

Кандидат физико-математических наук, преподаватель, Университет Западной Вирджинии Vyacheslav. Akkerman@mail.wvu. edu WV 26506-6106 г. Моргантаун, США

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

Ключевые слова: переопределенные системы дифференциальных уравнений, урав-^ нения Эйлера, Навье - Стокса, дифференциальные уравнения на поверхности, ОДУ, раз-сЗ мерность дифференциальных уравнений, задача Коши, уравнения в частных производных.

w pq

s Введение

а

Н

•< Нелинейные дифференциальные уравнения в частных производных (нелинейные уравнения tq математической физики) часто встречаются в различных областях математики, физики, меха-S ники, химии, биологии и в многочисленных приложениях [15]. Общее решение нелинейных урав-g нений математической физики удается получить только в исключительных случаях. Поэтому Л обычно приходится ограничиваться поиском и анализом частных решений, которые принято на© зывать точными решениями [5; 13].

Точные решения дифференциальных уравнений математической физики всегда играли и продолжают играть важнейшую роль в формировании правильного понимания качественных особенностей многих явлений и процессов в различных областях науки и техники [15]. Можно сказать, что невозможность решить эти уравнения является препятствием для дальнейшего изучения многих физических явлений и применения их на практике. Точные решения нелинейных уравнений наглядно демонстрируют физику и позволяют разобраться в механизме таких сложных нелинейных эффектов, как пространственная локализация процессов переноса, множественность или отсутствие стационарных состояний при определенных условиях, существование режимов с обострением и др. [7]. Научный интерес представляет даже поиск решений дифференциальных уравнений, которые не имеют ясного физического приложения. Подобные примеры могут быть успешно использованы в качестве «тестовых» задач при проверке корректности и оценке точности различных численных, асимптотических и приближенных методов. Допускающие точные решения модельные уравнения и задачи служат основой для разработки новых численных, асимптотических и приближенных методов, которые, в свою очередь, позволяют исследовать уже более сложные задачи, не имеющие точного решения, представимого в явном виде. Точные методы и решения необходимы также для разработки и совершенствования соответствующих разделов компьютерных программ, предназначенных для аналитических (символьных) вычислений (системы MATHEMATICA, MAPLE, CONVODE, MATHCAD и др.).

Существуют различные приемы для поиска точных решений уравнений математической физики, например, метод дифференциальных связей [5; 13; 15]. В данной работе мы развиваем один из этих методов, а именно прием редукции переопределенных систем дифференциальных уравнений, предложенный ранее в работах авторов [1; 3]. Решается обратная задача. Показывается, как, решая редуцированные уравнения на поверхности, можно составлять и находить решения исходной системы дифференциальных уравнений во всем объеме. Мы приводим также аналитические, модельные примеры и преобразовываем уравнения механики и гидродинамики к переопределенным с тем, чтобы к ним можно было бы применить данный метод.

1. Теория метода. Прямая задача

Для полноты изложения представим сначала вкратце основные теоретические результаты, полученные авторами ранее в работах [1; 3]. Здесь и далее, если не оговорено особо, мы будем предполагать, что все рассматриваемые функции являются гладкими и достаточное число раз непрерывно дифференцируемыми по своим аргументам. Рассмотрим в пространстве (г, t) систему из р дифференциальных уравнений в частных производных первого порядка относительно неизвестных 5к (г,t), V = 1...p, переопределенную одним независимым уравнением (на каком-нибудь решении, например):

И,

0 8S 8S Sv, —,—-, r, t 8r 8t

_ 0 v _

, v _ 1...p, k _ 1...p,

G

8S 8S

Sv,—, —, r, t 8r 8t

_ 0

(1) (2)

Перейдем в точке М в систему координат (т1, т2, п) на некоторой неподвижной поверхности (п -нормаль к поверхности; т1 и т 2 - единичные касательные к поверхности). Тогда уравнения (1), (2) запишутся в виде

И

( 8S 8S 8S 8S

v v' 8tj 3 8т2

8n 8t

_ 0 v _ 1...p к _ 1...p

(3)

GfS ^^^^ " v' дт1 ' dr2 ' dn ' dt

= 0

(4)

Выразим из уравнений (3) нормальные производные дяк/дп в явном виде

dSi

dn

k = Fk

f dS,, dS,, dS,, ') S,

v v

дт1 дт2 дt

Подставим их выражения (5) в формулу (4). Тогда

v = l...p k = l...p

(5)

G

(i)

f dS„ dS„ dS„ ')

S

vv

dzi ' dr2 ' dt

= 0

(6)

Продифференцируем уравнение (6) в направлении п и подставим вместо дяк/дп их выражения из (5). Тогда находим, что

G

(2)

f d2 s d2 s d2 s л

S

'dr1dt'dr2dt' dt2

= 0

(7)

Проделаем аналогичную процедурур раз. Получимр уравнений на поверхности вида

G

(i)

( dS dS dS Л

vv

S

у v' dzi ' dz2' dt j

= 0

G

(p)

S..

d pS„

dpS„ dpS„

dridtp-^ dr2dtp~ dtp

= 0

(8)

Мы нашли замкнутую систему изр поверхностных дифференциальных уравнений (8) строго вдоль границы рассматриваемой поверхности и такого же количества переменных Sk, к = 1...р, эволюционирующих во времени. Таким образом, может быть высказано следующее утверждение (возможно, являющееся перенесением на случай рассматриваемой переопределенной системы известного метода «подстановки»).

Утверждение 1. Любую систему из р дифференциальных уравнений в частных производных первого порядка (1), описывающих эволюционирующие во времени их решения Sk (г,t), к = 1...р, в евклидовом пространстве и переопределенную любым дифференциальным уравнением (2) первого порядка, можно преобразовать к системе из р дифференциальных уравнений и такого же числа неизвестных Sk (г,I), к = 1...р, на любой неподвижной поверхности, если в системе координат (т1, т2, п) для любой точки на этой поверхности все нормальные производные дSk/дп из системы (1) явно выражаются, как функции от переменных, задаваемых исключительно на рассматриваемой поверхности (5).

Формально ничего не мешает подобной процедурой получать больше, чем р уравнений на поверхности (8), то есть найти переопределенную систему уравнений уже на ней, и, следовательно, сократить размерность на поверхности и т. д. вплоть до полного решения, представимого в явном виде (см. часть 5). Однако это не означает, что подобной процедурой можно найти полное решение любой переопределенной системы уравнений [3]. Кроме того, сократить размерность можно не на каждой поверхности (см. [3] и часть 3).

Введем обозначения Ау = дSv/дt, V = 1...р . Представим с учетом этого обозначения выражения (5), (6) в виде

dS

-n = Fk (Av...), v = i...p, k = i...p,

(9)

G (1)( Л ...) = 0.

(10)

Продифференцируем уравнение (10) в направлении п и в полученном разложении обозначим слагаемые, содержащие наибольшие производные по времени Н,

^ да «з2 ^ >---+... = о

7=1 дЛ дпдг

Ч-1 V!

(11)

или, учитывая (9),

^ дС(1) дF д25 (2), ч

у --^—^ +... = о(2)(...) = 0.

у 2т=1 дЛ дЛ дН2 у '

V1,V7 =1 V V-,

(12)

Проделаем аналогичную процедуру р раз. Тогда систему уравнений (8), где выделены слагаемые, содержащие наибольшие производные по времени н, можно записать в виде

^ да(1) д^ д^ д'Б ), ч

у --+... = а(1)(...) = 0, I = 1...р

^ дЛп дЛ^ дЛV дН1 К ' ' р.

(13)

Продифференцируем каждое 1-е уравнение из системы уравнений (13) по времени Н (р - I) раз. Тогда получим следующую систему:

^ дС(1) дF дF дрБ У --—... —-^ + ... = 0, I = 1...р .

1^=1 дЛ-1 дЛ-2 дЛ- дН" ' ^

(14)

Система поверхностных уравнений (14) линейна относительно старших производных по времени дРБ- /дНр , V = 1...р . Условие того, что эти производные могут быть явно выражены из (14), выглядит следующим образом [2]:

ы *0,

где

(15)

а да «д^ да(1) ,

= М2Г-Л при ' >1 и ^ = ~дЛ~ при 7 =1.

(16)

В этом случае систему (14) можно представить в виде

дрБ, „ (др-1Б др-1Б

к=а " ■

дн"

дНр-1 'дт1днр_2'"

, V = 1...р, к = 1...р .

(17)

Таким образом, может быть предложено следующее утверждение.

Утверждение 2. Если для системы из р дифференциальных уравнений в частных производных первого порядка (1), описывающих эволюционирующие во времени их решения Бк (г, н), к = 1...р, в евклидовом пространстве и переопределенной некоторым дифференциальным уравнением (2) первого порядка, на некоторой неподвижной поверхности выполняется условие (15), то ее можно преобразовать к системе из р дифференциальных уравнений и такого же числа неизвестных Бк (г,н), к = 1...р, на рассматриваемой поверхности, причем старшие производные по времени от Бк (г, н) будут явно выражены через все остальные неизвестные.

а

К системе поверхностных уравнений (17) уже можно поставить соответствующую задачу Коши, чтобы выделить интересующее заранее решение. Начальные данные у этой задачи Коши не могут быть произвольными из-за наличия уравнения связи (2) и условия (15). Согласно общей теореме Коши - Ковалевской такая задача в случае аналитичности всех рассматриваемых функций (разложимости в степенной ряд) имеет единственное локальное решение [5; 13; 15]. Мы видим, что в случае, если выполняется условие (15), краевые условия задавать не нужно на бесконечной незамкнутой поверхности. Заметим также, что если система уравнений (1), (2) линейна относительно неизвестных функций и их производных, то редуцированные системы уравнений, получающиеся из нее, также линейны относительно своих неизвестных и с ними гораздо легче оперировать.

Аналог утверждения 2 для подвижной поверхности приведен в работе авторов [3].

2. Теория метода. Обратная задача

Представляет интерес следующая обратная задача. Не ограничивая общности, рассмотрим в евклидовом пространстве (г, t) серию поверхностей {г = ^ —<х>< h <+<х>}, параметризованную параметром h. Рассмотрим переопределенную систему (1), (2) изр дифференциальных уравнений в частных производных первого порядка относительно неизвестных Sk (г, t), у = 1...р . Пусть нормальные производные дSk|дп на этих плоскостях выражаются в явном виде (5). Тогда, как было показано, уравнения (1), (2) можно преобразовать к системе уравнений (8), эволюционирующих исключительно на этих поверхностях. Обратно, если мы найдем решение уравнений (8) на каждой из плоскостей {г = Щ для любого h, то при каких условиях «сшитое» из них решение во всем пространстве (г, t) будет решением системы уравнений (1), (2).

Исследуем сначала переопределенные системы ОДУ [16]. Рассмотрим систему из р обыкновенных дифференциальных уравнений первого порядка относительно неизвестных ^), у = 1...р, эволюционирующих во «времени» ¿, и переопределенную одним независимым уравнением (на каком-нибудь решении, например):

¿у = Fk (8у,t), у = 1...р, к = 1...р, (18)

,0 = 0. (19)

Продифференцируем по времени t уравнение (19). Имеем или, подставляя выражения из (18),

р DG (Sk, t) , ч DG (Sk, t)

Z^sT^f (Sk, t у^г1- -G (1) (Sk' t )-°.

Аналогично

p dG(1) (S., t) . . dG(1) (S., t)

LSr1 F(Sk,th—^r1 -G (s,t)=

V-l v

p dG(p-2) (Sk, t) , , DG(p-2) (Sk, t) ( , ,

Z—d^Fv(Sk,th—d* -G (Sk,t

P DG (p-1) (Sk, t) . , DG<*-1) ( Sk, t) Z J k' j F (Sk, t) +-XüJ - G ® (Sk, t)-0.

Рассмотрим следующую переопределенную систему уравнений относительно Sk (t), к = 1...p:

G (Sk, t) = G(0) (Sk, t) = 0, (24)

G(1) (Sk, t ) = 0, (25)

G (p-1) (Sk, t ) = 0. (26)

G(p) (Sk, t ) = 0. (27)

Предположим, что переопределенная система уравнений (24)-(27) имеет решение Sk = S* (t), k = 1...p . Пусть определитель матрицы

Г -G(i) ^

— , k = 1...p, i = 0...(p -1), (28)

V -Sk )

не равен тождественно нулю на этом решении. Подставим его в (24)-(27) и продифференцируем по времени t первыеp - 1 уравнений. Тогда

^ -G(0) (S*, t) .. -G(0) (S*, t)

У-U^Ls* +-= 0 (29)

tí -Sv v -t ,

^ -G(1) (S*,t) .. -G(1) (S* t)

У-U^-S* +-^^ = 0 (30)

tf -Sv v -t ,

^ -G(p-1) (S*, t К -G(p-1) (S*, t)

У-L^i^* +-L^i = 0 (31)

t! -Sv v -t •

Подставим решение Sk = S* (t) , k = 1...p в выражения (20)-(23). Имеем:

^ -G(0) (S*, t) . . . -G(0) (S* t) . .

У -S k, ' F(S*, t) +-= G« (S*, t) = 0, (32)

^ -G(1) (S*,t) . . . -G(1) (S* t) . .

y --Sf-1 Fv ( s; , t)+=G (2) (S*, t)=0, (33)

^ -G(p-1) (S*, t) . . . -G(p-1) (S*, t) . .

У Sr1 F (S*, t)+-= G M (Sl t) = 0 . (34)

Почленно вычтем из системы уравнений (29)-(31) систему (32)-(34) соответственно. Получим

yGSS%'-Fv(s;,()H, (35)

У^^^^МИ, (36)

=1

p

ydG;P-SS*^(s^* - Fv (S*,t)) = 0. (37)

В силу нашего предположения относительно определителя матрицы (28) из системы линейных уравнений (35)-(37) относительно неизвестных (Б'* -FV (Б*,н)) , V = 1...р , следует, что [2; 16]:

Б* = Fk (Б*,Н), V = 1...р, к = 1...р . (38)

Таким образом, мы можем высказать утверждение.

Утверждение 3. Любую систему из р обыкновенных дифференциальных уравнений первого порядка (18), описывающих эволюционирующие во времени решения Sk (t), к = 1...р, и переопределенную любым дифференциальным уравнением (19) первого порядка, можно

решить в явном виде из переопределенной системы р +1 неявных уравнений (24)-(27), если определитель матрицы (28) тождественно не равен нулю на всех решениях системы неявных уравнений (24)-(27).

Очевидно, что в условиях этого утверждения все решения находятся из переопределенной системы уравнений (24)-(27) и множество решений не может зависеть от непрерывного параметра. Если система (24)-(27) не имеет решений, то не имеет решений и переопределенная система ОДУ (18), (19). Уравнение (27) может быть даже тождеством. Если система (18), (19) линейная, то система (24)-(27) также линейная. Если определитель матрицы (28) тождественно не равен нулю, то в этом случае линейная система (24)-(27), а значит и (18), (19) имеют не более одного решения [2; 16]. Заметим также, что если система ОДУ (18), (19) не зависит неявно от времени, то в условиях утверждения 3 она имеет постоянное решение (см.: (29)-(31)).

Пример 1. Рассмотрим следующую переопределенную систему ОДУ (осциллятор):

и = —у, (39)

У = и, (40)

и + у = ^^) + sin(t). (41)

Продифференцируем (41) по переменной t и вместо производных по t подставим их выражения из (39), (40). Имеем

и — у = ^^) — sin(t). (42)

Проделаем аналогичную процедуру с уравнением (42). Тогда

—и — у = — ^^) — sin(t). (43)

Переопределенная система неявных уравнений (41)-(43) имеет только одно общее решение

и = со8(0 , у = sin(t). (44)

Определитель матрицы (28) для этой системы равен -2. Следовательно, согласно утверждению 3, выражения (44) являются единственным решением (39)-(41). Интерес представляет некоторое обобщение этого результата. А именно, если к системе (39)-(41) добавить произвольную внешнюю силу, а также заменить в уравнении связи (41) правую часть на произвольную функцию от времени (взятую, например, из измерений в эксперименте), то аналогичными рассуждениями получим, что решение новой системы уравнений будет находиться из линейной системы уравнений (41)-(43), но с другой правой частью (в виде конкретной формулы). Причем, если это решение не существует, то не существует решение и исходной обобщенной системы. А если существует, то оно единственно [2].

Рассмотрим теперь п-мерный случай. Рассмотрим евклидово пространство (г, х, t) и серию поверхностей {г = И, —<х< И <+да}, параметризованную параметром И. Рассмотрим переопределенную систему из р дифференциальных уравнений в частных производных первого порядка относительно неизвестных Бк (г,х,t), у = 1...р,

Hk

(dS dS dS 0 —-,—-,—-, Sv, z, x, t dt dz dx

- v - 1...p, k - 1...p, (45)

а

(дБ,, дБ,, дБ

дН дг дх

, Б,, 2, х, Н

= 0

(46)

Пусть нормальные производные дБк /дг на этих плоскостях выражаются в явном виде

дг

(

= Ft

дБ дБ

дН дх

, Б,, 2, х, Н

, V = 1... р, к = 1... р.

(47)

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

а

(1)

(

дН дх

,Б„, г, х, Н

= 0, V = 1...р,

а

(р)

(д рБ дрБ дБ дБ

-,-г...—~ , Бу, г, х, Н

дНр дхдНр-1 дН дх г

= 0.

(48)

При этом, если соответствующий определитель не равен тождественно нулю (15) {на всяком решении уравнений (45), (46)}, то система (48) корректна [3]. Для наших целей получим этой процедурой еще два уравнения на этих поверхностях

а

(р+1)

(др+1Б д р+1Б дБ дБ

Л

V V

дНр+1 ' дхдНр "' дН ' дх '

а

( р+2)

(др+2 б др+2 Б дБ дБ

V V

Б,,,г, х, Н

, Б„, г, х, Н

= 0,

= 0 .

(49)

(50)

дНр+2 'дхдН^1'" дН ' дх

V у

Пусть система (48)-(50) имеет решение Б* = Б* (г,х,Н), к = 1...р . Подставим его в (48)-(50). Тогда

а1

(1)

(

* -л*

дБ*,^, Б*, г, х, Н дН дх

= 0, V = 1...р,

а

(р)

(дрБ* дрБ* дБ * дБ* Л

-,-г-.—- ,—, Б*, г, х, Н

дНр дхдНр-1 дН дх

а

(р+1)

( др+1Б* д р+1Б* дБ* дБ*

V

а

( р+2)

дНр+1 ' дхдНр "' дН ' дх

ч

(др+2б* др+2б* дБ* дБ*

-, Б*, г, х, Н

= 0

= 0,

V V

Л

Б*, г, х, Н

= 0.

дНр+2 'дхдН^1'" дН ' дх ' ' '

ч у

Продифференцируем первые р + 1 выражений (51)-(53) по г. Следовательно,

( Л

У

У

да(к) дк (дБ * Л да(к)

д V V

( дкБ* Л дНк

дНк

V дг у

д

к (дБ* Л

д

( дрБ* Л дхдНк

дхдНр-1

V дг у

+... + -

да(к) дБ*

дБ* дг

да(к)

дг

= 0,

(51)

(52)

(53)

(54)

где к = 1...(р + 1). С другой стороны, согласно самому методу получения уравнений (48)-(50) мы имеем также р + 1 выражений

+

+

=1

f

I

Л

dG(k) dkM„

8G(k)

д kM„

Jd ks*Л

д —г V (dt" У

dt

dpS* ^ dxdtk

д

- +... + -

dxdtp-1

dG(k)

dS,

M

+ dGl = Gk+ц..) = о, dz v '

(55)

где k = 1...(р +1) и М, = Fl (дБ*/дt,дБ*/дх,Б*,г,х, t), V = 1...р, I = 1...р . Вычтем ражений (55) выражения (54). Тогда

почленно из вы-

(

Г

I

Л

dG(k) д kQv dG(k) д kQv

i дkS* ^ dtk

д VV

д.tk

(

д ps

_v

дхд^-1

Л

дхдГ

+... +

дG(k)

дS,

* Qv

= 0, k = 1... (p +1),

(56)

где Ql = М, - дБ*, /дг, I = 1...р .

Рассмотрим (56) как переопределенную систему линейных уравнений в частных производных от неизвестных функций Qi = Qi (г, х, t), г = 1...р, на рассматриваемых поверхностях |г = h,—х> < h < +<х>}, то есть г - просто параметр. Предположим, что она нашим методом снижается в размерности вплоть до переопределенной линейной системы ОДУ, которая имеет решение, представимое в явном виде. Для этого необходимо, чтобы серия соответствующих определителей вида (15) и определителей матрицы вида (28) не равнялась нулю. Как было показано, в этом случае линейная система (56) имеет не более одного решения. Как можно легко заметить, она уже имеет нулевое решение. Следовательно, в этих условиях других решений она не имеет. Это означает, что

Q, = M, - дБ*/дz = F, (дS* /дt, дS*v/дх, S*, z, х, t)-дS"¡/дz = 0 , v = 1...p, l = 1...p

(57)

то есть Б* = Б* (г, х, t) удовлетворяет выражениям (47), следовательно, и (45).

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

Утверждение 4. Любая, описывающая эволюционирующие во времени решения в евклидовом пространстве Бк (г,t), k = 1...р, система из р дифференциальных уравнений в частных производных первого порядка (1), переопределенная любым дифференциальным уравнением (2) первого порядка, и редуцированная изложенным выше методом, переопределенная система из р + 2 дифференциальных уравнений от этих же неизвестных на некоторой параметризованной серии поверхностей, охватывающих все пространство, имеют одинаковое множество решений, если на этих поверхностях переопределенная система р + 1 линейных дифференциальных уравнений (56) от неизвестных Q , г = 1...р, имеет единственное нулевое решение для любых функций Б* = Б* (г,х,t), к = 1...р, стоящих в коэффициентах системы (56) и являющихся решениями этой редуцированной системы из р + 2 дифференциальных уравнений.

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

Пример 2. Рассмотрим следующую переопределенную систему из двух уравнений в частных производных:

да

— = а- x +1, дх

да ~öt

2 2 = а - х .

(58)

(59)

Редуцируем (58), (59) на серию поверхностей {x = const}. Запишем (59) в виде

+

=1

д

а® а] = ^-а2 + х2 = 0. (60)

1,дН у дН

Продифференцируем (60) по переменной х и вместо производных по х подставим их выражения из (58). Имеем

а(2) (да, а^^-2а2 + 2 ( х - 1)а + 2 х = 0. (61)

Проделаем аналогичную процедуру с уравнением (61). Тогда

сю (1'аУ=1-4аг+(6х-4)а-2х(х-2)=°. (62)

Мы получили 3 = 1 + 2 уравнения (60)-(62), эволюционирующих на этих плоскостях. Решая уравнения (60)-(62), находим, что они имеют только одно общее решение а = х . Соответствующая переопределенная система линейных уравнений выглядит следующим образом:

^ - 2а0 = 0, (63)

дН

+ (-4а + 2(х -1)) 0 = 0. (64)

Если подставить в (63), (64) это решение а = х, получим следующую систему:

дд - 2х0 = 0, (65)

дд+(-2х - 2) д = 0. (66)

Очевидно, переопределенная система линейных уравнений (65), (66) имеет единственное нулевое решение. Следовательно, согласно утверждению 4 переопределенная система уравнений (58), (59) имеет только одно решение а = х, которое, очевидно, удовлетворяет (58), (59).

Согласно утверждению 2 мы имеем некоторое достаточное условие независимости уравнения связи. Утверждение 4 также предлагает достаточное условие для независимости уравнения связи, но оно более сильное в том смысле, что не требует определенного преобразования к задаче Коши редуцированной системы уравнений.

Заметим, что изначально в условиях утверждения 4 требование о наличии общих решений у переопределенной системы уравнений (45), (46) заранее не содержится. Достаточно только выписать какое-нибудь дополнительное соотношение (46) к системе (45), потом найти решение переопределенной параметризованной системы дифференциальных уравнений на серии поверхностей (48)-(50) и проверить, что соответствующая переопределенная система линейных уравнений (56) для этого решения имеет единственное нулевое решение. Тогда это решение, распространенное через серию поверхностей на все пространство, будет решением и исходной системы уравнений (45).

3. Новый метод переопределения уравнений математической физики

Рассмотрим динамическую систему уравнений в каноническом виде [6]:

^ = Н (67)

Ж дq¡

^ = -дИ (68)

dt др. '

где i = 1...n , q = q0, p = p0 при t = 0, Н(q,p,t) - функция Гамильтона. Условие сохранения фазового объема означает, что [6]

^ =1 (69)

д( qo, Ро)

Мы имеем дополнительное уравнение связи (69) для уравнений (67), (68), если их рассматривать как систему уравнений в частных производных в евклидовом пространстве (q0, p0, t). Очевидно, на поверхностях вида {t = const} сократить размерность нашим методом нельзя, поскольку тогда из уравнения связи (69) получатся верные тождества (8). Произвольную систему ОДУ можно свести к динамической вида (67), (68). Действительно, рассмотрим произвольную систему ОДУ

f-= ^,0. (70)

где i, j = 1...n . Введем функцию

Н (q, p, t ) = -£f (q, t) р, (71)

i=1

и новые неизвестные p i = 1...n, которые определим из уравнений

dp дИ , ч

=--(72)

dt дq¡'

Очевидно, система ОДУ (70) ввиду определения (71) может быть записана в виде

dq дИ „ ч

=----(73)

dt др. '

Таким образом мы систему ОДУ (70) свели к динамической системе ОДУ (72), (73). Если функция Н (q, p, t) не зависит от времени t явно, то есть И = Н (q, p ) (а это всегда можно сделать, преобразовав систему (70) к автономной), то мы имеем еще одно дополнительное соотношение [6]

Н (q, p) = Н (q0, p0). (74)

Рассмотрим уравнения Эйлера несжимаемой жидкости в трехмерном потоке [7]

ди , „ч дР

+ (u • V) и + —= 0, (75)

д^ дг

divu = 0, (76)

где и| = U0 (r) - начальное распределение скорости u; Р - давление. Рассмотрим переменные Лагранжа r (r0, t) (разметку) [11]. Тогда для каждой движущейся частицы жидкости r0 выполняются следующие соотношения:

du дР

du = ~1РР, (77)

, (79)

dt дг

- _ и . (78)

dt

Очевидно, уравнения (77), (78) можно записать в виде

Ли __дН dt дг

^ _дН, (80)

dt ди

где Н _ Н (г, и, {)_ Р (г, ^ + и 2/2. Пусть г _ г (га, и а, ^ и и _ и (га, иа, ^ - решение системы (79), (80), рассматриваемой как система ОДУ с начальными данными г _ г(га,иа,0)_ га и и _ и (га,иа ,0) _ иа. Тогда для каждой частицы жидкости г0 выполняется: г _ га _ г0, и _ иа _ ио (г0) при t = 0, то есть движение всей жидкости можно представить как движение трехмерного многообразия в шестимерном фазовом пространстве (г иа). Уравнение (69) означает, что

д(г, и )

^ ' _ 1. (81)

д(га , иа )

Уравнения (77), (78) или (79), (80) можно записать в виде

ди (га , иа , t)_ дР (82)

(83)

дt дг

дг (га , иа , t) / А

—---_ и (г , и , t)

аа

д*

Все частные производные по г в выражениях (76) и (82) преобразовываются к производным по га и иа, например, согласно формулам [14]:

дР дг дР ди дг ди дР дг дР ди дг ди

-_--и-_--или-_--и-_----(84)

дг дг дг дг дг дг ди ди дг ди ди дг

а а а а а а а а

Таким образом, мы имеем переопределенную систему уравнений из 8 уравнений (76), (81)-(83) в частных производных и 7 неизвестных функций вида г _ г (га, иа, I), и _ и (г, tи (га,иа, t) и Р _ Р (г, I)_ Р (га, иа, I) в пространстве (га, иа, I). Конкретное решение фиксируется следующим образом: га _ г0, иа _ и0 (г0) [1; 3]. В этом случае полная производная скорости и находится как [11]

du ди (га, и а, t) ди (г (га, и а, t), t) ди дг (га, и а, t) ди ди ди

-_-_-_--I---_--Ъ и--( 85)

& дt дt дt дt дг дt дг

в полном соответствии с ее определением.

Рассмотрим уравнения Навье - Стокса несжимаемой жидкости в трехмерном потоке [7; 12]

ди / дР „

— + (и • V)и + — _ _уУх ш , (86)

д^ ' дг

ш _ V х и, (87)

divu _ 0, (88)

где и| _ и0 (г) ; V - вязкость. Введем «поправку» а к вектору скорости жидкости и с тем, чтобы обобщить систему уравнений (86), (87) следующим образом [1; 3]:

д(и + а) , (!„ . „Ч У7\Л. , „Ч , Ы О , 1 ««2

+ ((и + а)•V)(u + а) +VIР + 1и2 _ 1(и + а)2 |_ 0, (89)

— = [а х ш + и х (У х а) + а х (У х а )] + гУ х ш, (90)

д?

divu = 0. (91)

Начальные данные у функций а можно взять нулевыми [1].

Векторное поле и + а формально можно рассмотреть как невязкий гидродинамический поток

с давлением

P + u 72 -(u + а )2/2~L Пусть

1 1 2

V = и + а, Р = Р + 2и2 --(и + а)2. (92)

Тогда и = V - а . Подставим это в (89)-(91):

дУ

^ + (У-У)У + ур = 0, (93)

За = [а х (у х (У - а )) + (У - а )х(Ух а ) + а х(Ух а )] + уУх(Ух( У - а )), (94)

(У - а) = 0. (95)

Как было показано выше, уравнения (93) можно переопределить, и вместе с уравнениями (94), (95) мы получим переопределенную систему уравнений Навье - Стокса в пространстве переменных (га, Уа, ?) от неизвестных функций вида г = г (га, Уа, ?), У = У (г, ?) = У (га, Уа, ?) , а = а (г, ? ) = а (га, Уа, ?) и р = р(г, ? ) = р(га, Уа, ?) [1; 7]. Скорость и и давление Р можно тогда определить из (92). Конкретное решение (86)-(88) фиксируется следующим образом: г = га = г0, У = У = и0 (г) при ? = 0 или г„ = г0, У = и0 (г0).

Таким образом, если мы найдем нашим методом решения переопределенной системы уравнений, полученной из (93)-(95), в пространстве переменных (га,Уа,?) , и если этих решений, допустим, счетное количество, то среди них должны находиться решения уравнений Навье - Стокса в явном виде, зависящие от начальных данных (задача Коши), если таковые существуют.

Как известно, к решению системы уравнений (67), (68) сводятся уравнения в частных производных первого порядка от одной неизвестной функции [16]. Найдя решение в явном виде динамической системы (67), (68) в пространстве начальных данных, мы найдем решение и этого класса уравнений.

4. Задача Коши для редуцированной системы дифференциальных уравнений

Исследуем теперь задачу Коши при понижении размерности дифференциальных уравнений. Рассмотрим, не ограничивая общности, евклидово пространство (г, х, ?) и серию поверхностей |z = h,-<х>< h <+да}, параметризованную параметром И. Рассмотрим систему из р дифференциальных уравнений в частных производных первого порядка относительно неизвестных Sv (z, х,?),

V = 1...р,

Hk

f dS dS dS 0 } —,—-,—-, Sv, z, x, t dt dz dx

= 0, v = 1...p, k = 1...p . (96)

V у

Пусть к системе уравнений (96) заданы начальные данные Коши

^ = б: (^х), V = 1...р . (97)

Пусть нормальные производные д8к/дz на этих плоскостях выражаются в явном виде

дг

дБ дБ г,

_ Р _V__V. Б г х t

р дг ' дх , Б, ^ X, t

, V _ 1... р , к _ 1... р .

(98)

Выпишем пока произвольное дифференциальное соотношение

G

(дБ дБ дБ

V V V

с* дг дх

Б,,, г, х, г

_ 0.

(99)

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

G

(1)

(

с^, ^, Бу, г, х, * дг дх

_ 0 V _ 1...р

а

(р)

(дрБ дрБ дБ дБ

-,-г...—~ , Бг, г, х, г

дгр дхдгр_1 дг дх У

_0.

(100)

При этом, если соответствующий определитель не равен тождественно нулю (15), то система (100) корректна [3]. Для наших целей получим этой процедурой еще одно уравнение на этих поверхностях

а

(р+1)

(д ръ1Б д ръ1Б дБ дБ

\

V V

,Б„, г, х, г

_ 0.

дгръ1 ' дхдгр "' дг ' дх

ч у

Поставим для переопределенной системы уравнений (100), (101) следующую задачу Коши:

д 7 Б.

(101)

дг7

д ■'Б0

-^(г, х), V _ 1...р , ' _ 0...р,

д*' У ) и, J У,

(102)

где

'д'Б0^

дг

дг7

д7 д7

(

дБ дБ

V V

дг ' дх

, Б , г, х, г

V,к _ 1...р, ' _ 0...р _ 1.

(103)

Пусть система (100), (101) имеет решение Б* _ Б* (7, х, г), к _ 1...р . Тогда аналогично формуле (56) можно получить, что

где

г

I

са(к) скяу да(к) дkQv

да(к)

(дкБ* ^ дгк ' ( срб* \ дхд*к_х +... + дБ* ^

V V д*к у

дхдгр_'

_ 0, к _ 1...р,

(104)

Я, _ р

г дБ* дБ.

, Б*, г, х, г дг дх

_ дБ*/дг, V _ 1...р, I _ 1...р .

(105)

Рассмотрим выражения (104) как систему из р линейных дифференциальных уравнений в частных производных первого порядка относительно неизвестных (г, х, г), V _ 1...р . Из предположения (103) следует, что начальные данные у неизвестных (г, х, г), V _ 1...р нулевые. Сравнивая уравнения (104) и (100), мы также можем утверждать, что, если соответствующий определи-

г _0

д

г _0

тель не равен тождественно нулю (15), то система (104) корректна, и соответствующая задача Коши имеет единственное локальное решение [5; 13; 15]. Очевидно, это будет тогда нулевое решение, то есть Б* = Б* (7, х, ?) удовлетворяет выражениям (98), следовательно и (96). Таким образом, предлагается следующее утверждение.

Утверждение 5. Любая, описывающая эволюционирующие во времени решения в евклидовом пространстве (г,?), к = 1...р, система из р дифференциальных уравнений в частных производных первого порядка (96), переопределенная любым дифференциальным уравнением (99) первого порядка, для которой поставлена задача Коши (97), и редуцированная изложенным выше методом, переопределенная система (100), (101) из р + 1 дифференциальных уравнений от этих же неизвестных на некоторой параметризованной серии поверхностей, охватывающих все пространство, для которой поставлена задача Коши (102), (103), имеют локально одинаковые решения, если выполняется условие (15) для любых функций Б* = Б* (7,х,?), к = 1...р, являющихся решениями этой редуцированной переопределенной системы из р + 1 дифференциальных уравнений.

Это утверждение является некоторым уточнением утверждения 2. Отметим также, что изначально в условиях утверждения 5 (как и утверждения 4) требование о наличии общих решений у переопределенной системы уравнений (96), (99) заранее не содержится. Но в отличие от утверждения 4 утверждение 5 описывает только локальные свойства решений переопределенных систем дифференциальных уравнений. Заметим также, что вообще для переопределенной системы уравнений (96), (99) задача Коши (97) не может содержать произвольные начальные данные. Аналогично это можно утверждать для переопределенной системы (100), (101) и задачи Коши (102), (103).

5. Дальнейшая редукция уравнений вплоть до полного решения, представимого в явном виде

Требование на неравенство нулю некоторого определителя также возникает в следующей ситуации. Рассмотрим переопределенную систему из р + 1 дифференциальных уравнений в частных производных первого порядка относительно неизвестных (х) , V = 1...р, х = (х1,...хт):

Я„

öS

v

öx

, S, x

= 0, v = 1...p, k = 1... (p +1).

(106)

Будем сокращать размерность переопределенной системы уравнений (106). Продифференцируем выражения (106) p раз по переменной х1. Получим (p + 1)2 редуцированных уравнений от p(p + 2) неизвестных вида

S öSv 3p+'Sv

Sv , fcT ^ , V = P .

(107)

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

Pa(Qß, x )

Ö(i+

.)

x =

дхЧ.дх'т

относительно неизвестных вида

H

öS-, Sv, x öx

= 0, v = 1...P , k = 1... (P +1),

(108)

Öu+...+: )S

Qß =-:-~ , v = 1... p .

öx(1...öxJm~

(109)

Здесь а _ а(к, /1,../т)_ 1.. АН, Р _ Р (V,'1,...'т) _ 1...ИБ - функции от мультииндексов такие, что

Я^,0,...0) _ Б , V _ 1...р , и

/1 _ 0...р, /2 _ 0... (р +1)2 _ 1

/3 _ 0...

->2

(р + 1)2 _ 1

... / _ 0...

(т_1)

(р +1)2 _ 1

7 _ 0...(р + 1) , 72 _ 0...(р + 1)2, '3 _ 0...(р +1)2 , ... 7т _ 0... (р +1)

,(т_1

(110) (111)

Если решение системы (106) существует, то оно удовлетворяет системе неявных уравнений (108). Очевидно, количество уравнений (108), учитывая (110), равно

2(т_1)

Мн _(р + 1)-(р + 1>(р +1)2 •...•(р +1)2 _(р +1) а количество неизвестных (109), учитывая (111), равно

2+2+22 +...+2(т_1)) _

_( р + 1)2

(112)

N _

р •(( р +1)+1) • ((р +1)2 +1) •.. • ((р + 1)2(т1) +1) _

_ ((р +1) _ 1) • ((р +1)+1) • ((р +1)2 +1) •... • ((р+1)2(т1) +1) _ = ((р + 1)2 _ 1) • ((р + 1)2 + 1) • ... • ((р + 1)2(т1) + 1) _ (р + 1)2(т1)" _ 1 _ (р + 1)

(113)

то есть, сравнивая (112) и (113), ИБ _ Ин _ 1, что в полном соответствии с нашим предположением. Однако из (108)—(111) видно, что реально некоторые высшие производные (109) не входят в выражения (108), то есть на самом деле от некоторых неизвестных <2Р выражения Ра не зависят. Рассмотрим матрицу

АаР _

ГдР^

vдQpу

, а _ 1... Ин, р _ 1... N .

(114)

Рассмотрим расширенную переопределенную систему неявных уравнений вида

Ра(Яр , х)_

д(4 +... + 'т )

дх1..дхт

относительно неизвестных вида

н

дБ

V

дх

Б„, х

_ 0, V _ 1...р , к _ 1... (р +1),

(115)

д('1 + ...+7т Б

Яр_1.7 ' , V _ 1...р .

(116)

дх/1..^

Но функции от мультииндексов а _а(к, 11,...1т), Р _Р(V,'1,...7т), в отличие от (110), (111), определяются из условий:

/, _ 0... р+1, ,2 _ 0...( р+1)2, /3 _ 0...( р+1)22, ^ т _ 0... (р+1)2

'1 _ 0... (р + 1)+ 1, '2 _ 0... (р + 1)2 + 1, '3 _ 0...(р +1)22 +1, ^

'т _ 0... (р + 1) + 1.

(117)

(118)

Система уравнений (115) включает в себя систему уравнений (108). Неизвестных (116) несколько больше неизвестных (109). Выскажем и обоснуем следующее утверждение.

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

бым дифференциальным уравнением первого порядка (106), и редуцированная изложенным выше методом, переопределенная расширенная система неявных уравнений вида (115) от неизвестных вида (116) имеет одинаковое множество решений, если ранг матрицы (114) (на каждом решении этой расширенной системы неявных уравнений) равен количеству неизвестных Qр, реально присутствующих в уравнениях (108). Подсчет показывает, что это число не более

N"al < NH

( Р +!)

i+h

(p+1)

(119)

Поскольку мы допускаем, что ранг матрицы (114) равен N"а1, значит, существует N"а1 независимых уравнений из системы (108) от такого же числа неизвестных. Обозначим их как

pai(QPl, x ) = 0, l =1

= 1... n:

(120)

Здесь а, = а(к',,...¡'т), р. = р(V,^,...) . Продифференцируем уравнения (120) по переменной х , 1 < 5 < т . Имеем

f ^ + < Q, x) = 0

*=1 dQß, dx

dx„

, l = 1...NS

(121)

С другой стороны, из определения уравнений (115) следует, что

päl(Qß, x )=h^ß,

SP* п ,ÖP«i(öA' x)

f =1 dQß,

dx.

= o, i = 1... n:

(122)

где а, = а (к',/1',...(¡5 +1)..л'т), Д = р (V,^,...(^ +1).../т) . Вычтем почленно из выражений (122)

выражения (121). Тогда

h Q

Qs. -

dQß Л

_w

ex

i' v

= o, i = 1... n:

(123)

В силу нашего предположения относительно независимости уравнений (120) из выражений (123) следует

(

Qßi dx.

= 0 или Q ~ =

dQß

^ ßi 1 _ 1 Л Г real

Pi ~ dx

-, l = 1...N"ai.

(124)

Таким образом, мы пришли к выводу, что для каждой переменной х , 1 < 5 < т , выполняется соотношение (124). Отсюда, учитывая, что по определению QP(V,0,...0) = ^ (^), V = 1...р, следует

Q

dQ

ß(v, j -1,... jm )

dj1 Q

ß(v,0,...jm )

d( j1+j2 )Q

ß( v,0,0... jm )

ß( ^ Л-- jm )

dx

dx/1

dx/1 dxj

d( j1 +...+jm ),

Qß( v._

ß( v,0,...0)

d

i( j1 +...+jm )

S,,

dx/1 ...dxJm

1 m

dxj1...dxjm

1m

(125)

В итоге мы показали, что функции (x), V = 1...р, являются решением (106), а величины (109) являются соответствующими частными производными этого решения.

1

Для нашего примера переопределения уравнений Навье - Стокса из части 3 мы имеем р =10, т = 7. Следовательно, общее число неявных уравнений (108) равно Ин _ (р +1) « 10128. Это очень большое число. Оно возникает в основном из-за большой размерности евклидового пространства переменных т = 7. В работах [1; 3] приведены другие способы переопределения уравнений гидродинамики, где т = 4. Тем не менее утверждение 6 предлагает некоторое достаточное условие для существования сразу полного решения переопределенных систем дифференциальных уравнений.

Количество редуцированных уравнений и неизвестных можно сократить, если на каждом шаге сокращения размерности получать более одного дополнительного уравнения связи. Например, в случае, если мы имеем систему из р ОДУ, переопределенную I уравнениями связи, то достаточно продифференцировать уравнения порядка р/1 раз, и мы получим систему из порядка р2/1 неявных уравнений и такого же числа неизвестных вместо р2. Исследование данного вопроса дано в приложении А.

Заключение

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

В данной статье, в отличие от [1; 3], приводятся новые достаточные условия для независимости уравнения связи. А именно, требуется исследовать переопределенную систему линейных однородных уравнений в частных производных. Редуцированные системы уравнений, получающиеся из нее, также линейны относительно своих неизвестных, и все исследование сводится к существованию отличных от нуля определителей, составленных из их коэффициентов. Трудность заключается в их выписывании, понимании их структуры и преобразовании самих линейных уравнений при редукции. То есть необходимо произвести большое количество символьных операций. Для этого существуют специальные программы (МАТНЕМАТ1СА, MATHCAD и др.) и среды, например, Лисп. Если это достаточное условие будет выполнено, то исходную переопределенную систему дифференциальных уравнений в объеме можно будет преобразовать к системам дифференциальных уравнений на серии поверхностей, параметризованным некоторым параметром. «Сшитое» из них решение во всем пространстве будет решением этой исходной системы уравнений. Фактически одну из переменных рассматриваемой системы уравнений можно будет тогда преобразовать в параметр для редуцированной системы уравнений. Производные по этому параметру будут отсутствовать.

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

В утверждениях 3-6 требование о наличии заранее общих решений у переопределенных систем уравнений не содержится. В утверждениях 1, 2 говорится, что если переопределенная система имеет решение, то оно содержится среди решений редуцированной системы уравнений, и его можно выделить, если поставить соответствующую задачу Коши для редуцированных систем уравнений. Исследование на совместность необходимо в методе дифференциальных свя-

зей, где нужно подобрать так, чтобы уравнение связи содержало как можно больший произвол [13]. Там необходимо решить дополнительно систему уравнений, чтобы выписать уравнение связи. В нашем методе, наоборот, выгодно, чтобы переопределенная система имела как можно меньше общих решений (не более чем счетное множество, см. приложение Б). Тогда их возможно найти, например, с помощью утверждения 6. В наших работах [1; 3] предлагается переопределение уравнений гидродинамики таким образом, чтобы множество решений этих систем переопределенных уравнений содержало решение, удовлетворяющее заранее заданной задаче Коши для этих уравнений. Множество таких решений может состоять всего из одного элемента. Исследование на совместность там очевидным образом выполнено. Если мы все-таки выпишем уравнение связи, содержащее больший произвол, то с определенного момента редуцированная система уравнений просто не будет дальше сокращаться в размерности, и мы не сможем выписать общее решение в явном виде. Однако произвол может возникнуть в следующей ситуации. Если мы редуцируем переопределенную систему с помощью утверждения 6 и выпишем систему неявных уравнений, то может оказаться, что счетного множества решений оно не имеет (требование на ранг матрицы не выполняется). Но при этом может оказаться, что какие-то высшие производные от неизвестных функций Б^, (х), V _ 1...р, х _ (х1,...хт) (109) находятся в явном виде. Тогда, решая очевидные уравнения, можно учесть произвол в виде краевых и начальных условий.

Классическое исследование переопределенных систем на совместность, начатое Карта-ном, Спенсером, Яненко и др. [9; 13], в данной работе не проводится. Это и не нужно, потому что мы ссылаемся на полученные ранее в работах [1; 3] и в данной работе уже достаточно строго переопределенные системы уравнений математической физики, например, гидродинамики, и с ними работаем.

В 3-й части нашей работы приведены по-новому преобразованные, переопределенные системы уравнений Эйлера, Навье - Стокса и уравнений аналитической механики [8]. Главная особенность метода их переопределения заключается в том, что ищется решение сразу в пространстве начальных данных и времени, то есть при последовательной редукции этих уравнений может быть сразу найдено их решение (задачи Коши). В последней части статьи показывается, что это решение (или хотя бы частное) может быть чрезвычайно сложным, недоступным даже для решения современными ЭВМ. Можно предположить, что невозможность описания жидкости при больших числах Рейнольдса (развитой турбулентности) или при других условиях связана с чрезвычайно сложной структурой явного представления решения уравнений, описывающих движение этой жидкости.

В настоящее время создаются программные комплексы, которые уже эффективно решают сразу полную систему дифференциальных уравнений, где требуется как можно точнее найти решение задачи для конкретных систем дифференциальных уравнений сразу во всем объеме изменения пространства переменных или серьезное упрощение процесса численного получения этого решения. Интересуют частные, точные решения этой полной системы уравнений (они могут быть чрезвычайно сложными), в частности, для тестирования кодов [10]. В связи с этим в статью добавлено приложение А, которое систематически описывает способ получения явного представления решений любых систем дифференциальных уравнений, если эти системы каким-то образом переопределены.

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

Приложение А. Алгоритм нахождения решений переопределенных систем

дифференциальных уравнений

На основе утверждения 6 можно составить следующий общий алгоритм [4]. Пусть мы имеем переопределенную систему уравнений, получающуюся из (106), вида

Pa(QP, x ) = -

^U+...+m)

dxl..dxm

относительно неизвестных вида

H

^, Б,, x Эх

= 0, v = 1... p, k = 1... (p + n),

(А. 1)

Э(j + ...+jm

QP = j - / , V = 1... Р .

Эх/1... dxJm

( v,

(А. 2)

Здесь а = а(k,i1,..im) = 1...NH , Р = Р(v, j1,...jm) = 1...NS - функции от индексов (мультииндек-сов) такие, что

Qi

P(v,0,...0)

= S,, v = 1... p,

Pa( k ,0,...0) (QP , X ) = Hk (Qp( v,1,...0) ,...QP(v,0,...1), QP(v,0,...0), X) , k = 1... ( Р + П ) ,

¡1 = 0... (N1 - 1) , 1г = 0...( N2 - 1) , ... im = 0...( Nm - 1) , k = 1... (p + n) ,

j = 0...N1, j2 = 0...N2, ... jm = 0.. .N m , V = 1... p .

(А. 3) (А. 4) (А. 5) (А. 6)

Имеем также

NH =( p + n ) N1N 2 .. .Nm,

Ns = p•(N1 +1)(N2 +1)...(Nm +1) .

Пусть

(А. 7) (А. 8)

a=a(k, i1,...im )= k + ¡1 (p + n ) + ¡2 (p + n ) N1 + ... + ¡m (p + П) N1...Nm-1;

P = P( V, j1,...jm )= V + j1p + j2 p (N1 + 1) + ... + jmp (N1 + 1)...( Nm-1 + 1) . Рассмотрим матрицу

ЛхР =

v^pJ

a

= 1...Nh, P = 1...NS .

(А. 9)

Пусть ее ранг равен nrsea количеству неизвестных QP, реально присутствующих в уравнениях (А. 1). Это число не более (119)

p ( m 1 ^ Nr:aI < ы„. р .. 1+У —

I 1=1N J

S H

(p+n)

(А. 10)

Рассмотрим расширенную переопределенную систему неявных уравнений вида

P." (QP-, X ) =

+ ...+im )

H

^, Sv, X Эх

= 0, V = 1... p, k = 1... (p + n),

(А. 11)

относительно неизвестных вида

Л + ...+jm ) О

бр* •v=• (А 12)

Но функции от мультииндексов а* = а* (k,i1,...im) = 1...NH , р* = р* (v, j1,...jm) = L..NJ в отличие от (А. 3)-(А. 6) пусть определяются из условий:

бр-(v,0,...0)= О , V = 1...p , (А. 13)

P (к,«^)^ x)= (бр* (v,1,...0) ,..бр* (v,0,...1), бр*( v,0,...0), *) , k = (P + И) , (А- 14)

i1 = 0...N1, i2 = 0...N2, ... im = 0...Nm , к = 1...(p + n), (А. 15)

j = 0...(N1 +1), j2 = 0...(N2 +1), ... jm = 0...(Nm +1), v = 1...p . (А. 16)

Пусть

а* = а* (k,) = k+Л (р + п) + /2 (р + п)(N +1) +... + im (р + п)(У +1)...(Nm_l +1),

Р* = Р* (V, ) = V + Л р + 72 Р (N1 + 2) +... + ]тР (N1 + 2)... (Лт_1 + 2).

Имеем также

ЛН = ( Р + п)(N1 + 1)(Л +1)... (Nт +1), (А. 17)

У; = Р •(N1 + 2)(N2 + 2)...(Nт + 2). (А. 18)

Дифференцируя выражение (А. 11) (см. (122)) по переменной х ,, = 1...т, находим, что

р Л+1 N2+1 Лт+1 яр дР * (Q., х)

р. (о,. *)=хх х ...I =(а. 19)

V=1 71 = ° 72 = ° 7т =° дЦз* дх

где

а = а* (к./1....(/; +1)..Лт), Р = Р* (V.71....(Л +1)...7т), ^ = 1...т .

Здесь предполагается, что I, < Л, -1, < У,.

Любое из уравнений (А. 11) (см. (А. 19)) с индексом а =а (к, Л1,...Лт) содержит неизвестные вида (А. 12) с индексами Р* =Р*(V,71,...7), где 71 < Л +1, ... 7 < 1т +1 и ° < 7 +... + / < I +... + Л +1.

1 т 1 т

Очевидно, тогда

дР

а*(к=(-1)...<т ) = ° (А 2°)

Я0Р*(V,Л,...7 =У, +1,...]т )

Рекуррентное соотношение (А. 19), удобное для получения уравнений (А. 11), корректно.

Способ решения. Определяем из (А. 19) неявные уравнения, пользуясь (А. 13), (А. 14), и решаем их. Вычисляем ранг матрицы (А. 9). Вычисляем вручную количество переменных (А. 1°)! и сравниваем с рангом матрицы (А. 9). Если они совпадают, то мы имеем решение системы дифференциальных уравнений вида (1°6).

Очевидно, нас интересует только случай NН > у. Можно получить следующие оценки:

ЛН >(р + п)(тр/п)т,

где минимум реализуется при У1 и У2 и ... и Лт и тр/п.

Приложение Б. О множестве решений переопределенных систем дифференциальных уравнений

Рассмотрим следующий пример. Пример 3.

да да ^

дt дх

5а да

--1--= 2а.

д? дх

(Б. 1) (Б. 2)

Очевидно, переопределенная система (Б. 1), (Б. 2) имеет общее решение а = Const exp( x +1).

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

Рассмотрим переопределенную систему из p + n дифференциальных уравнений в частных производных первого порядка относительно неизвестных Sv (x) , v = 1...p, x = (xj,...xm)

H

dS

v

Эх

S.., x

= 0, v = 1...p, k = 1...(p + n) .

(Б. 3)

Пусть система (Б. 3) имеет общие решения вида = Sv (х,а), V = 1...р, где а - параметр. Подставим это решение в (Б. 3) и продифференцируем выражения (Б. 3) по параметру а. Тогда

г

I

днк dqv , дн,

d

dS

_v_

v dx у

+ -

dx dS,

-Qv

= 0, к = 1... (p + n),

(Б. 4)

где Qv = дSv/да, V = 1... р .

Рассмотрим (Б. 4) как переопределенную систему линейных однородных дифференциальных уравнений относительно неизвестных Qv, V = 1...р, где коэффициенты определяются конкретным решением системы (Б. 3). Очевидно, это система имеет общее нулевое решение. Согласно приложению А данную систему можно редуцировать до переопределенной системы линейных уравнений. И если требование на ранг матрицы (А. 9) будет выполнено, то это нулевое решение будет единственно. Следовательно, тогда Qv = дSv/да = 0, V = 1...р, то есть общие решения вида Sv = Б (х, а) существовать не могут.

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

Легко проверить, что для примеров 1 и 2 данное утверждение выполняется, а для примера 3 не выполняется.

СПИСОК ЛИТЕРА ТУРЫ

1. Аккерман, В. Б. Снижение размерности в уравнениях гидродинамики / В. Б. Аккерман, М. Л. Зайцев // Журнал вычислительной математики и математической физики. - 2011. - Т. 51, №2 8. - С. 1518-1530.

2. Беклемишев, Д. В. Курс аналитической геометрии и линейной алгебры / Д. В. Беклемишев. - М. : Физматлит, 2005. - 304 с.

3. Зайцев, М. Л. Гипотеза об упрощении переопределенных систем дифференциальных уравнений и ее применение к уравнениям гидродинамики / М. Л. Зайцев, В. Б. Аккерман // Вестник ВГУ Серия: Физика. Математика. - 2°15. - №№ 2. - С. 527.

4. Зайцев, М. Л. Еще один способ нахождения частных решений уравнений математической физики / М. Л. Зайцев, В. Б. Аккерман // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. - 2°16. - №> 6 (37). - С. 119-127.

5. Курант, Р. Уравнения с частными производными / Р. Курант. - М. : Мир, 1964. - 83° с.

6. Ландау, Л. Д. Теоретическая физика : учеб. пособие : в 1° т. Т. I. Механика / Л. Д. Ландау, Е. М. Лифшиц. - М. : Наука, 1988. - 216 с.

7. Ландау, Л. Д. Теоретическая физика: Гидродинамика / Л. Д. Ландау, Е. М. Лифшиц. - М. : Наука, 1986. -Т. VI. - 736 с.

8. Лурье, А. И. Аналитическая механика / А. И. Лурье. - М. : ГИФМЛ, 1961. - 824 с.

9. Полянин, А. Д. Методы решения нелинейных уравнений математической физики и механики / А. Д. Полянин, В. Ф. Зайцев, А. И. Журов. - М. : Физматлит, 2°°5. - 256 с.

1°. Полянин, А. Д. Справочник по нелинейным уравнениям математической физики: Точные решения / А. Д. Полянин, В. Ф. Зайцев. - М. : Физматлит, 2°°2. - 432 с.

11. Самарский, А. А. Разностные методы решения задач газовой динамики / А. А. Самарский, Ю. П. Попов. -М. : Наука, 198°. - 352 с.

12. Седов, Л. И. Механика сплошной среды : в 2 т. / Л. И. Седов. - М. : Наука, 1978. - Т. 1. - 492 с. ; Т. 2. -568 с.

13. Сидоров, А. Ф. Метод дифференциальных связей и его приложения к газовой динамике / А. Ф. Сидоров, В. П. Шапеев, Н. Н. Яненко. - Новосибирск : Наука, 1984. - 271 с.

14. Схоутен, Я. А. Тензорный анализ для физиков / Я. А. Схоутен. - М. : Наука, 1965. - 456 с.

15. Тихонов, А. Н. Уравнения математической физики / А. Н. Тихонов, А. А. Самарский. - М. : Наука, 1966. - 742 с.

16. Федорюк, М. В. Обыкновенные дифференциальные уравнения / М. В. Федорюк. - СПб. : Лань, 2°°3. - 448 с.

REFERENCES

1. Akkerman VB., Zaytsev M.L. Snizhenie razmernosti v uravneniyakh gidrodinamiki [Dimension Reduction in Fluid Dynamics Equations]. Zhurnal vychislitelnoy matematiki i matematicheskoy fiziki [Computational Mathematics and Mathematical Physics], 2011, vol. 8, no. 2, pp. 1418-1430.

2. Beklemishev D.V. Kurs analiticheskoy geometrii i lineynoy algebry [Course of Analytical Geometry and Linear Algebra]. Moscow, Fizmatlit Publ., 2005. 304 p.

3. Zaytsev M.L., Akkerman V.B. Gipoteza ob uproshchenii pereopredelennykh sistem differentsialnykh uravneniy i ee primenenie k uravneniyam gidrodinamiki [Hypothesis on Reduction of Overdetermined Systems of Differential Equations and Its Application to Equations of Hydrodynamics]. Vestnik VGU. Seriya: Fizika. Matematika, 2015, no. 2, pp. 5-27.

4. Zaytsev M.L., Akkerman V.B. Eshche odin sposob nakhozhdeniya chastnykh resheniy uravneniy matematicheskoy fiziki [Another Method for Finding Particular Solutions of Equations of Mathematical Physics]. Vestnik Volgogradskogo gosudarstvennogo universiteta. Seriya 1, Matematika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2016, no. 6 (37), pp. 119-127.

5. Kurant R. Uravneniya v chastnykh proizvodnykh [Partial Differential Equations]. Moscow, Mir Publ., 1964. 830 p.

6. Landau L.D., Lifshitz E.M. Teoreticheskaya fizika: ucheb. posobiye: v 101. T. I. Mekhanika [Course of Theoretical Physics, Vol. 1: Mechanics]. Moscow, Nauka Publ., 1988. 216 p.

7. Landau L.D., Lifshitz E.M. Teoreticheskaya fizika: Gidrodinamika. T. VI [Course of Theoretical Physics. Vol. 6: Fluid mechanics]. Moscow, Nauka Publ., 1986. 736 p.

8. Lurye A.I. Analiticheskayamekhanika [Analytical Mechanics]. Moscow, GIFML Publ., 1961. 824 p.

9. Polyanin A.D., Zaytsev V.F., Zhurov A.I. Metody resheniya nelineynyh uravneniy matematicheskoy fiziki Imekhaniki [Methods for Solution of Equations of Mathematical Physics and Mechanics]. Moscow, Fizmatlit Publ., 2005. 256 p.

10. Polyanin A.D., Zaytsev VF. Spravochnikpo nelineynym uravneniyam matematicheskoy fiziki: Tochnye resheniya [Handbook on Nonlinear Equations of Mathematical Physics: Exact Solutions]. Moscow, Fizmatlit Publ., 2002. 432 p.

11. Samarskiy A.A., Popov Yu.P. Raznostnye metodyresheniya zadachgazovoy dinamiki [Difference Methods for Solving Problems of Gas Dynamics]. Moscow, Nauka Publ., 1980. 352 p.

12. Sedov L.I. Mekhanika sploshnoy sredy: v 2 t. [Mechanics of Continuous Media. Vols. 1, 2]. Moscow, Nauka Publ., 1978, vol. 1. 492 p.; vol. 2. 568 p.

13. Sidorov A.F., Shapeev V.P., Yanenko N.N. Metod differentsialnykh svyazey i ego prilozheniya к gazovoy dinamike [Method of Differential Connections and Its Application to Gas Dynamics]. Novosibirsk, Nauka Publ., 1984. 271 p.

14. Schouten J.A. Tenzornyy analiz dlya fizikov [Tensor Analysis for Physicists]. Moscow, Nauka Publ., 1965. 456 p.

15. Tikhonov A.N., Samarskiy A.A. Uravneniyamatematicheskoy fiziki [Equations of Mathematical Physics]. Moscow, Nauka Publ., 1966. 742 p.

16. Fedoryuk M.V. Obyknovennye differentsialnye uravneniya [Ordinary Differential Equations]. Saint Petersburg, Lan Publ., 2003. 448 p.

REDUCTION OF OVERDETERMINED DIFFERENTIAL EQUATIONS OF MATHEMATICAL PHYSICS

Maxim Leonidovich Zaytsev

Postgraduate Student,

Nuclear Safety Institute of Russian Academy of Sciences mlzaytsev@gmail. com

Bolshaya Tulskaya St., 52, 115191 Moscow, Russian Federation

Vyacheslav Borisovich Akkerman

Candidate of Physical and Mathematical Sciences, Assistant Professor, West Virginia University Vyacheslav. Akkerman@mail.wvu. edu WV 26506-6106 Morgantown, USA

Abstract. A technical method of reducing the overdetermined systems of differential equations is further extended. Specifically, the fundamentals and validity limits of the method are identified, and the method is justified within its validity domain. Starting with an overview of the previous results, we subsequently employ them in deriving and justifying the new outcomes. In particular, overdetermined systems of ordinary differential equations (ODE) are studied as the simplest case. It is demonstrated that, if a determinant deviates from zero for an ODE system, then the solution to this system can be found. Based on this, we subsequently arrive to a more general statement for a system of partial differential equations (PDE). On a separate basis, a Cauchy problem for reduced overdetermined systems of differential equations is considered, and it is shown that such a problem cannot be with arbitrary initial conditions. It is also shown and substantiated how to employ a Cauchy problem to reduce the dimension of PDEs. A novel approach of how to transform ODE and PDE systems (such as Euler and Navies-Stokes equations as well as the analytical mechanics system of equations) into the overdetermined systems is presented. Finally, the results are generalized in such a manner that it is shown how to reduce an overdetermined system of deferential equations to that having a complete and explicit solution. The work also includes two appendices. The first appendix presents the algorithm of searching for a solution to an overdetermined system of differential equations, in particular, by means of the computational approaches. The second appendix is devoted to the study of the variety of the solutions to an overdetermined system of equations. In particular, it is shown that a certain condition for the determinant, associated with this system of equations, breaks the possibility, that such a variety of solutions can depend on

a continuous factor (for instance, from Cauchy conditions). For instance, it could be not more than a countable set. The paper is concluded with a brief summary, where the major results of the work are listed again and discussed, including their potential practical applications such as developing and testing of new computer codes to solving systems of differential equations.

Key words: overdetermined systems of differential equations, Euler and Navier-Stokes equations, differential equation on the surface, ODE, dimension of differential equations, Cauchy problem, partial differential equations.