Методы > Уравнения с частными производными, уравнения математической физики > А. Д. Полянин, К. Е. Казаков. Решение дифференциальных уравнений с помощью инвариантов

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

© 2008   А. Д. Полянин, К. Е. Казаков
Московский государственный технический университет имени Н.Э. Баумана

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

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

Оглавление

Введение

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

1. Симметрии. Общая схема использования инвариантов для решения математических уравнений

1.1. Симметрии. Преобразования, сохраняющие вид уравнений. Инварианты. В природе и технике часто встречаются геометрические симметрии: зеркальная симметрия, когда одна половина предмета зеркально-симметрична другой, и центральная симметрия, когда предмет переходит сам в себя при повороте относительно некоторого центра. Однако понятие симметрии гораздо шире, чем в указанных случаях, в общем случае под симметрией понимается неизменность при какой-либо процедуре не только предметов, но и физических явлений, математических формул, уравнений и т. д. В данном учебном пособии пойдет речь о симметрии алгебраических и дифференциальных уравнений и о том как использовать это свойство, чтобы находить решения или упрощать уравнения.

Под симметриями математических уравнений будут пониматься преобразования, сохраняющие вид уравнений. Ниже приведены примеры конкретных уравнений, которые сохраняют свой вид при некоторых простых преобразованиях.

Пример 1. Рассмотрим биквадратное уравнение

(1.1)

Замена

приводит к точно такому же уравнению (т. е. уравнение (1.1) сохраняет вид при преобразовании x = -).

Два других преобразования

также сохраняют вид уравнения (1.1), поскольку после умножения на ⁴ получим

Пример 2. Вид дифференциального уравнения

(1.2)

не изменится, если сделать любое преобразование

(1.3)

поскольку для всех трех указанных преобразований в результате получим точно такое же уравнение

Далее будет показано, что преобразования, сохраняющие вид уравнений, позволяют "размножать" решения уравнений.

Инвариант преобразования — функция (отличная от постоянной), которая сохраняет вид при действии данного преобразования. Инвариант преобразования может зависеть от независимых и зависимой переменной и их производных (если идет речь о дифференциальных уравнениях). Для пояснения понятия инвариантов, которые сохраняют вид при заданных преобразованиях, рассмотрим несколько простых конкретных примеров.

Пример 3. Преобразование одновременного сдвига по обоим осям

где a — любое число, имеет инвариант

Если x — независимая переменная, а y — зависимая переменная, то другими инвариантами этого преобразования являются производные

Пример 4. Преобразование одинакового изменения масштабов длины по обоим осям

где a ⁄= 0 — любое число, имеет инвариант

Если x — независимая переменная, а y — зависимая переменная, то имеются также более сложные инварианты, которые зависят от производных и сохраняют вид при действии данного преобразования, например,

1.2. Общая схема использования инвариантов для решения математических уравнений. Ниже изображена принципиальная схема исследования математических уравнений, которая основана на поиске преобразований, сохраняющих вид уравнений, с последующем переходом в уравнении от исходных переменных к новым переменным — инвариантам преобразований.

После выполнения указанных действий уравнение часто упрощается или приводится к разрешимому виду. Важно отметить, что данная схема с успехом может применяться для самых различных типов математических уравнений (см. далее разделы 2–4).

Для лучшего понимания и усвоения идей использования инвариантов для решения математических уравнений при изложении дальнейшего материала применяется схема "от простого к сложному": сначала кратко излагаются соответствующие результаты для алгебраических уравнений, потом — для обыкновенных дифференциальных уравнений, а затем — для нелинейных уравнений математической физики.

2. Алгебраические уравнения и системы уравнений

2.1. Алгебраические уравнения, содержащие четные степени. Рассмотрим алгебраическое уравнение

(2.1)

содержащие только четные степени. Биквадратное уравнение является частным случаем уравнения (2.1) при n = 2.

Замена

(2.2)

приводит точно к такому же уравнению для (говорят, что уравнение (2.1) инвариантно относительно преобразования (2.2)). Отсюда следует, что если уравнение (2.1) имеет решение x = x₁, то оно имеет также другое решение x = − x₁.

Возводя в квадрат обе части (2.2), получим простейшую алгебраическую функцию, которая сохраняет вид при преобразовании (2.2):

(2.3)

Эта функция и является инвариантом преобразования (2.2). Если выбирать инвариант (2.3) за новую переменную, z = x², то уравнение (2.1) порядка 2n преобразуется к уравнению порядка n:

Таким образом в данном случае переход от исходной переменной x к инварианту преобразования (2.3) z = x² позволяет упростить рассматриваемое уравнение (понизить его порядок в два раза).

2.2. Возвратные уравнения. Возвратное алгебраическое уравнение четной степени имеет вид

(2.4)

Левая часть этого уравнения называется возвратным многочленом.

Замена

(2.5)

преобразует уравнение (2.4) в точно такое же уравнение Отсюда следует, что если уравнение (2.4) имеет корень x = b, то оно имеет также другое решение x = 1∕b.

Простейшее возвратное уравнение — квадратное уравнение

Делим его на x²: Результат удобно представить в форме уравнения первого порядка где

(2.6)

— простейший инвариант преобразования (2.5).

Теорема 1. В общем случае возвратное уравнение порядка 2n допускает понижение порядка с помощью подстановки (2.6). В результате получается алгебраическое уравнение порядка n.

Пример 1. Рассмотрим возвратное уравнение четвертого порядка

Делим его на x²:

(2.7)

Учитываем, что

В результате уравнение  (2.7) путем замены (2.6) (основанной на инварианте преобразования (2.5) сводится к уравнению второго порядка

Теорема 2. В общем случае возвратное алгебраическое уравнение нечетного порядка

имеет корень x = -1, а его левая часть допускает представление где Q_2n(x) — возвратный многочлен степени 2n.

Пример 2. Возвратное уравнение третьего порядка

может быть представлено в виде

Из теоремы 2 следует, что возвратное уравнение порядка 2n + 1 после деления на (x + 1) и введения новой переменной (2.6) сводится к уравнению порядка n.

2.3. Системы алгебраических уравнений, симметричные относительно перестановки аргументов. Многочлен P(x,y) от двух переменных называетс симметрическим, если он не меняется при перестановке аргументов: P(x,y) = P(y,x).

Замечание. В терминах преобразований симметрический многочлен определяется как многочлен, сохраняющий вид при преобразовании x = , y = .

Простейшие симметричные многочлены

(2.8)

называются элементарными. Эти многочлены являются простейшими алгебраическими инвариантами при перестановке аргументов.

Любой симметрический многочлен от двух переменных может быть единственным образом выражен через элементарные многочлены.

Для решения систем двух алгебраических уравнений

где P и Q — симметричные многочлены, полезно в качестве новых переменных использовать элементарные симметрические многочлены (2.8). Подобные системы обладают следующим свойством: если система имеет решение x = x₀, y = y₀, то она имеет также решение x = y₀, y = x₀.

Пример 3. Рассмотрим нелинейную систему алгебраических уравнений

(2.9)

Система (2.9) не меняется при перестановке аргументов.

Переходя в (2.9) от x, y к переменным (2.8) и учитывая формулы

получим

(2.10)

Исключая из уравнений u, приходим к квадратному уравнению

Дальнейшая процедура определения решений элементарна и мы ее опускаем.

Пример 4. Рассмотрим нелинейную систему алгебраических уравнений

(2.11)

Переходя к переменным (2.8) и учитывая равенство x³ + y³ = (x + y)³ - 3xy(x + y), получим

Исключив w, приходим к кубическому уравнению

(2.12)

Отметим, что прямое исключение из системы (2.11) переменной y приводит к существенно более сложному чем (2.12) уравнению шестого порядка:

Замечание. Более подробное изложение теоретического материала по алгебраическим уравнениям и системам уравнений, обладающими свойствами симметрии (т. е. допускающими преобразования, сохраняющими вид уравнений), и многочисленные примеры можно найти в [1, 2].

2.4. Выводы. Описанный метод исследования, основанный на использовании инвариантов, позволяет решать и упрощать некоторые типы алгебраических уравнений (понижать их порядки) и систем уравнений.

Упражнения к разделу 2.

1. Свести алгебраическое уравнение пятой степени x⁵ + ax⁴ + bx³ + bx² + ax + x = 0 к квадратному уравнению. Найти корни исходного уравнения при a = b = -1/2.

2. Решить уравнение четвертой степени ax⁴ + bx³ + cx² - bx + a = 0. Указание: учесть, что преобразование x = -1∕ сохраняет вид уравнения; найти и использовать инвариант этого преобразования.

3. Найти общий вид алгебраического уравнения четной степени, которое сохраняет вид при преобразовании x = -1∕. Указать замену, понижающую порядок полученного уравнения. Подробно рассмотреть случай уравнения шестого порядка и свести его к уравнению третьего порядка.

4. Решить систему уравнений

5. Решить систему уравнений

6. Решить систему уравнений

7. Решить систему уравнений

3. Обыкновенные дифференциальные уравнения

3.1. Преобразования, сохраняющие вид уравнения, и их инварианты. Говорят, что обыкновенное дифференциальное уравнение

(3.1)

инвариантно относительно (обратимого) преобразования

(3.2)

если в результате подстановки (3.2) в (3.1) получим точно такое же уравнение

(3.3)

Функция F в уравнениях (3.1) и (3.3) одинакова.

Преобразования, сохраняющие вид уравнения, позволяют "размножать" его решения. Действительно, пусть известно частное решение

(3.4)

уравнения (3.1). Поскольку уравнение (3.1) после перехода к новым переменным (3.2) сохраняет такой же вид, то преобразованное уравнение (3.3) имеет решение

(3.5)

Возвращаясь в (3.5) к старым переменным по формулам (3.2) (их предварительно надо разрешить относительно и ) получим решение уравнения (3.1), которое, в случае общего положения будет отличаться от исходного решения (3.4).

Пример 1. Уравнение третьего порядка

(3.6)

имеет частное решение

Преобразование

(3.7)

сохраняет вид этого уравнения, поэтому преобразованное уравнение

имеет решение = e. Подставляя сюда старые переменные, полученные обращением формул (3.7), имеем новое решение уравнения (3.6): содержащее две произвольные постоянные A и b.

Функция I(x,y,y′) (отличная от постоянной) называется инвариантом преобразования (3.2), если она сохраняется при этом преобразовании

Замечание. Если I = I(x,y,y′) инвариант преобразования (3.2), то Ψ(I), где Ψ — произвольная функция, тоже является инвариантом данного преобразования.

3.2. Процедура понижения порядка уравнений при n ≥ 2 (приведение уравнений к разрешимому виду при n = 1). Конкретизируем описанную в разделе 1 общую схему использования инвариантов для анализа математических уравнений применительно к обыкновенным дифференциальным уравнениям.

На первом этапе ищется преобразование

(3.8)

которое сохраняет вид уравнения (3.1). Преобразование (3.8) обязательно должно зависеть от одного свободного параметра a [a₁,a₂] (само исходное уравнение (3.1) от этого параметра не зависит).

На втором этапе для уравнений второго и более высоких порядков (при n ≥ 2) строятся два функционально-независимых инварианта преобразования (3.8):

(3.9)

На третьем этапе инварианты (3.9) выбираются в качестве новых переменных для уравнения (3.1), т. е. делается преобразование

Указанная процедура приводит к понижению порядка уравнения на единицу.

Для уравнений первого порядка (при n = 1) на третьем этапе в (3.1) надо сделать замену

В результате уравнение приводится к разрешимому виду (к уравнению с разделяющимися переменными).

3.3. Используемые преобразования. Процедура определения инвариантов. Далее будем использовать только простейшие преобразования вида

и их комбинации

(3.10)

В этом случае между производными имеют место линейные соотношения 

(3.11)

Коэффициенты преобразования A₁, A₂, B₁, B₂ определяются из условия инвариантности рассматриваемого уравнения (коэффициенты A₁, A₂, B₁, B₂ должны зависеть от одного свободного параметра a).

Справедливо следующее утверждение. Пусть преобразование (3.10) сохраняет вид некоторого уравнения, имеющего частное решение (3.4). Тогда это уравнение имеет также решение

Первый инвариант I₁ получаем путем исключения параметра a из равенств (3.10), а второй инвариант I₂ — путем исключения параметра a из одного из соотношений (3.10) и первого соотношения (3.11).

3.4. Анализ конкретных обыкновенных дифференциальных уравнений.

Пример 2. Уравнение второго порядка, не зависящее явно от y:

(3.12)

Это уравнение не меняется при произвольном сдвиге по зависимой переменной: y y + a (соответствует y = + a), где a — свободный параметр. При этом из трех величин x, y, y′ две остаются неизменными:

Это инварианты уравнения (3.12), т. е. Их выбираем за новые переменные: В результате получим уравнение первого порядка

Пример 3. Автономное уравнение второго порядка, не зависящее явно от x:

(3.13)

Это уравнение не меняется при произвольном сдвиге по независимой переменной: x x + a, где a — свободный параметр. При этом из трех величин x, y, y′ две остаются неизменными:

Это инварианты уравнения (3.13), т. е. I₁ = y, I₂ = y′. Их выбираем за новые переменные: В результате получим уравнение первого порядка: uu′ = F(z,u).

Пример 4. Нелинейное уравнение второго порядка

(3.14)

не меняется при растяжении зависимой переменной: y ay. При этом две комбинации из трех величин x,y,y′ остаются остаются неизменными:

Это инварианты уравнения (3.14). Их выбираем за новые переменные: Дифференцируя u по x, имеем Исключая с помощью этого соотношения вторую производную y″ в (3.14), приходим к уравнению первого порядка

Замечание. При F(x,u) = g(x) + f(x)u исходное уравнение (3.14) является общим линейным однородным уравнением второго порядка. Указанным преобразованием оно сводится к уравнению первого порядка с квадратичной нелинейностью.

Пример 5. Рассмотрим нелинейное уравнение второго порядка

(3.15)

Ищем инвариантное преобразование в виде

(3.16)

Подставив (3.16) в (3.15), а затем сократив на a, получим

Требуя совпадения с уравнением (3.15), получим соотношение для определения параметра b:

(3.17)

При этом параметр a остается свободным.

Подставим (3.17) в (3.16), а затем из второго соотношения исключим параметр a с помощью первого соотношения. Имеем

Таким образом получен первый инвариант преобразования (3.16)–(3.17) в виде

(3.18)

Чтобы найти второй инвариант, вычислим производную

Исключая здесь параметр a с помощью второго соотношения (3.16), получим второй инвариант

(3.19)

Для понижения порядка исходного уравнения надо принять за новые переменные инварианты (3.18)–(3.19):

(3.20)

С одной стороны

(3.21)

с другой стороны

(3.22)

Приравнивая (3.21)–(3.22), имеем

Подставив в (3.15), после элементарных преобразований получим уравнение первого порядка с разделяющимися переменными:

Замечание. Аналогичными свойствами обладает более общее нелинейное уравнение второго порядка

Преобразование (3.20) приводит его к уравнению первого порядка:

Пример 6. Рассмотрим теперь нелинейное уравнение первого порядка

(3.23)

Левая часть этого уравнения не меняется при преобразованиях вида

(3.24)

где a, b, c — произвольные постоянные. Подставим (3.24) в аргумент правой части уравнения (3.23):

(3.25)

Чтобы уравнение (3.23) было инвариантным относительно преобразования (3.24) в (3.25) надо положить

(3.26)

Эти соотношения можно рассматривать как систему двух линейных алгебраических уравнений первого порядка для определения коэффициентов b и c. При этом коэффициент a остается произвольным. Таким образом уравнение (3.23) инвариантно относительно преобразования (3.24), (3.26), а первым инвариантом этого преобразования является аргумент правой части уравнения (3.23). Поэтому в уравнении (3.23) надо сделать замену

(3.27)

Разрешим (3.27) относительно y, а затем продифференцируем по x. Заменив затем производную y′ на f(z) (это следует из (3.23) и (3.27)), после элементарных преобразований приходим к уравнению с разделяющимися переменными

В таблице 1 представлены некоторые обыкновенные дифференциальные уравнения второго порядка, допускающие понижение порядка путем использования простейших инвариантных преобразований. Для уравнений первого порядка (когда функция F(u,v,w) не зависит от третьего аргумента) уравнения, указанные в таблице 1, решаются путем перехода от y к новой зависимой переменной z = I₁(x,y), где I₁ — первый инвариант.

Таблица 1. Некоторые обыкновенные дифференциальные уравнения, допускающие понижение порядка, и их инварианты

Уравнение	Инвариантное преобразование	Первый инвариант	Второй инвариант

Приведенные в таблице 1 результаты легко обобщаются на нелинейные уравнения произвольного порядка.

3.5. Выводы. Описанный метод исследования обыкновенных дифференциальных уравнений использует идеи метода группового анализа, но значительно проще последнего. Содержание разд. 3.2 (основанного на результатах [3–5]) на начальном этапе может излагаться как общий эвристический подход, применимость которого проверяется на конкретных примерах. Для применения данного метода надо уметь решать лишь самые простые алгебраические уравнения (и системы) и уметь дифференцировать, в то время как при использовании метода группового анализа на промежуточных этапах приходится рассматривать уравнения с частными производными (т. е. приходится выходить за рамки стандартного курса по обыкновенным дифференциальных уравнениям). В качестве еще одного достоинства этого простого метода следует отметить, что здесь практически не приходится вводить новых понятий, которыми изобилует метод группового анализа. Данный метод позволяет единообразно описать разрешимых (или допускающих понижение порядка) обыкновенных дифференциальных уравнений значительно больше, чем разобрано в подавляющем большинстве существующих учебников.

Упражнения к разделу 3.

1. Решить уравнение y′ = y(kx + ln y).

2. Решить уравнение xy′ = f(ky + ln x), где f = f(z) — произвольная функция.

3. Решить обобщенно-однородное уравнение y′ = ax^k−1exp(yx^{− k}).

4. Решить обобщенно-однородное уравнение xy′ = yf(xⁿy^m).

5. Понизить порядок уравнения (1.2) тремя различными способами, используя преобразования (1.3).

6. Понизить порядок уравнени xy″ − y′ = 0 тремя различными способами.

7. Понизить порядок уравнени y″ − (y′)² = 0 тремя различными способами.

8. Понизить порядок уравнения Эмдена — Фаулера y″ = Axⁿy^m.

9. Понизить порядок уравнени y″ = Axⁿe^λy.

10. Понизить порядок уравнени y″ = Ayⁿe^λx.

11. Найти преобразование, понижающее порядок обобщенного уравнения Эмдена — Фаулера y″ = Axⁿy^m(y′)^k.

12. Найти преобразование, понижающее порядок уравнения  y″ = Ayⁿe^λx(y′)^m.

13. Найти преобразование, понижающее порядок уравнения  y″ = Axⁿe^λy(y′)^m.

14. Написать обобщения уравнений, приведенных в табл. 1, на случай уравнений произвольного порядка.

4. Дифференциальные уравнения с частными производными

4.1. Преобразования, сохраняющие вид уравнения, и их инварианты. Будем рассматривать точечные (обратимые) преобразования вида

(4.1)

Говорят, что уравнение с частными производными

(4.2)

инвариантно относительно преобразования (4.1), если в результате подстановки (4.1) в (4.2) получим точно такое же уравнение

(4.3)

Функция F в уравнениях (4.1) и (4.2) одинакова.

Величина I(x,t,w) называется инвариантом преобразования (4.1), если она сохраняется при этом преобразовании

Замечание. Если I = I(x,t,w) — инвариант преобразования (4.1), то Ψ(I), где Ψ — произвольная функция, тоже является инвариантом данного преобразования.

4.2. Процедура построения точных решений. Конкретизируем описанную в разделе 1 общую схему использования инвариантов для анализа математических уравнений применительно к дифференциальным уравнениям с частными производными.

На первом этапе ищется преобразование

(4.4)

которое сохраняет вид уравнения (4.2). Преобразование (4.4) обязательно должно зависеть от одного свободного параметра C [C₁,C₂] (причем само уравнение (4.2) от этого параметра не зависит).

На втором этапе строятся два функционально-независимых инварианта преобразования (4.4):

(4.5)

На третьем этапе инварианты (4.5) выбираются в качестве новых переменных и точное решение уравнения (4.2) ищется в виде

(4.6)

где Φ(z) — функция, удовлетворяющая обыкновенному дифференциальному уравнению, которое получается путем подстановки (4.6) в (4.2) (соответствующее решение называется инвариантным).

Как и ранее, в качестве преобразований (4.4) будем рассматривать только простейшие преобразования 

(4.7)

(4.8)

и их комбинации.

4.3. Примеры построения точных решений нелинейных уравнений математической физики.

Пример 1. Рассмотрим нелинейное уравнение нестационарной теплопроводности

(4.9)

Это уравнение не меняется при произвольном сдвиге по зависимым переменным

(4.10)

где C — свободный параметр, а k и λ — некоторые заданные числа (которые могут выбираться произвольно). Исключая из первых двух соотношений (4.10) параметр C находим один из инвариантов

Вторым инвариантом здесь является I₂ = w = . В силу (4.6) уравнение (4.9) допускает решение вида

(4.11)

Подставив (4.11) в (4.9), приходим к обыкновенному дифференциальному уравнению

Интегрируя дважды, получим его решение в неявном виде где C₁ и C₂ — произвольные постоянные.

Замечание 1. Решения вида (4.11) называются решениями типа бегущей волны [6, 7].

Замечание 2. Решения типа бегущей волны допускают уравнения общего вида, которые не зависят явно от независимых переменных:

(4.12)

Подставляя (4.11) в (4.12), получим автономное обыкновенное дифференциальное уравнение относительно функции Φ(z):

Пример 2. Рассмотрим уравнение нестационарной теплопроводности с нелинейным источником степенного вида

(4.13)

Ищем инвариантное преобразование масштабирования. Подставив (4.8) в (4.13) и умножив все члены на b∕c, имеем

Требуя совпадения с уравнением (4.13), получим два соотношения Выразим параметры a и c через b: Подставляя эти величины в (4.8), находим инвариантное преобразование

(4.14)

содержащее свободный параметр C = b.

Перейдем теперь к определению инвариантов. Из первых двух соотношений (4.14) исключаем b:

(4.15)

Из второго и последнего соотношений (4.14) исключаем b:

(4.16)

Соотношение (4.6) с учетом (4.15), (4.16) определяет вид автомодельное решения

(4.17)

Подставив (4.17) в (4.13), после несложных преобразований приходим к обыкновенному дифференциальному уравнению

Замечание. Преобразование (4.8) часто удобнее записывать в виде

а затем определять значения постоянных k и m, при которых сохраняется вид исходного уравнения (C — свободный параметр). Соответствующие таким преобразованиям точные решения называются автомодельными [6–8].

Примеры некоторых автомодельных решений, указанные в [9–11], приведены в таблице 2 (используется сокращение ОДУ — обыкновенное дифференциальное уравнение).

Таблица 2. Некоторые нелинейные уравнения математической физики, которые допускают автомодельные решения

Уравнение	Название уравнения	Вид решения	Итоговое ОДУ
	Нестационарное уравнение теплопроводности
	Уравнение Бюргерса
	Уравнение фильтрации
	Волновое уравнение
	Уравнение теплопроводности с источником
	Уравнение околозвукового газового потока
	Уравнение Кортевега – де Фриза
	Уравнение пограничного слоя

Аналогичный подход с успехом может использоваться для построения точных решений систем уравнений с частными производными.

Пример 3. Рассмотрим систему уравнений стационарного ламинарного гидродинамического пограничного слоя на плоской пластине:

(4.18)

Сделаем в (4.18) растяжение независимых и зависимых переменных по правилу

(4.19)

Умножив полученные уравнения на подходящие постоянные множители, имеем

(4.20)

Потребуем, чтобы вид уравнений преобразованной системы (4.20) совпал с видом уравнений исходной системы (4.18). Это условие дает два линейных алгебраических уравнения: n - m - k + 1 = 0, -2k - m + 1 = 0. Разрешив их относительно m и n, получим

(4.21)

где показатель k может быть выбран произвольно. Подставим (4.21) в (4.19):

Находим отсюда инварианты путем исключения параметра C:

Решение ищется в виде I₂ = U(I₁), I₃ = V (I₁), т. е.

или

(4.22)

где k — произвольная постоянная. Подставив (4.22) в исходную систему (4.18), получим систему обыкновенных дифференциальных уравнений для функций U = U(ζ), V = V (ζ):

Пример 4. Рассмотрим уравнение Клейна–Гордона с нелинейностью экспоненциального вида

(4.23)

Ищем инвариантное преобразование в виде комбинации преобразований сдвига и масштабирования:

(4.24)

Подставив (4.24) в (4.23) и умножив на b², получим

Чтобы это уравнение совпало с (4.23), надо положить Отсюда имеем и преобразование (4.24) принимает вид

(4.25)

где b > 0 — произвольная постоянная.

Исключая из соотношений (4.24) параметр b, находим инварианты

(4.26)

Соотношение (4.6) с учетом (4.26) определяет вид автомодельное решения

(4.27)

Подставляя (4.27) в (4.23), получим обыкновенное дифференциальное уравнение

Замечание. Точно таким же образом ищется точное решение стационарного уравнения теории горения с экспоненциальным тепловыделением

Пример 5. Рассмотрим нелинейное уравнение нестационарной теплопроводности

(4.28)

Ищем инвариантное преобразование в виде комбинации преобразований сдвига и масштабирования:

(4.29)

Подставив (4.29) в (4.28) и поделив на c, получим

Требование совпадения с уравнением (4.28) дает одно условие где c — произвольная постоянная. При этом другой параметр b также остается произвольным, т. е. полученное инвариантное преобразование

(4.30)

является двухпараметрическим.

Поскольку для применимости используемого метода требуется однопараметрическое инвариантное преобразование положим

(4.31)

где вид функции f будет определяться далее. Исключая параметр c из двух последних соотношений (4.30) с учетом (4.31), имеем

(4.32)

Это равенство будет определять инвариант I(t,w) = I( ,), если выполняется соотношение

(4.33)

которое можно рассматривать как функциональное уравнение. Нетрудно проверить, что уравнение (4.33) имеет решение

(4.34)

где A — произвольная постоянная. Подставив (4.31) в (4.30) с учетом (4.34), получим

(4.35)

При A ⁄= 0 выразим из второго соотношения (4.35) параметр c через t и подставим в оставшиеся соотношения. В результате находим инварианты

Используя теперь формулу (4.6), определяем вид точного решения уравнения (4.28):

(4.36)

Подставив (4.36) в (4.28), получим обыкновенное дифференциальное уравнение для функции Φ(z):

Замечание. Преобразование (4.30) при b = c сохраняет вид уравнения (4.28), но не имеет инвариантов. Отсюда следует, что не каждое преобразование, сохраняющее вид уравнения, имеет инварианты.

Таблица 3. Инвариантные решения, которые могут быть получены путем использования комбинаций преобразований сдвига и масштабирования, сохраняющих вид уравнений (C > 0 — произвольная постоянная)

N	Инвариантные преобразования	Вид инвариантных решений	Примеры уравнений
1
2			Уравнения см. в таблице 2
3
4
5
6
7
8

4.4. Выводы. Описанный метод исследования нелинейных уравнений с частными производными является естественным обобщением метода подобия [6, 8] и использует идеи метода группового анализа, но значительно проще последнего. Содержание разд. 4.2 (основанного на результатах [7, 10–14]) на начальном этапе может излагаться как общий эвристический подход, применимость которого проверяется на конкретных примерах. Для успешного использования данного метода надо уметь решать лишь простейшие алгебраические уравнения и уметь дифференцировать (и не требуется понимать и использовать сложную терминологию группового анализа дифференциальных уравнений).

Основные ограничения описанного метода:

1) не всегда можно найти преобразование (отличное от тождественного), которое сохраняет вид уравнения;

2) не всякое преобразование, сохраняющее вид уравнения, можно эффективно использовать (преобразование должно зависеть от свободного параметра);

3) не всегда существует инвариант заданного преобразования.

Упражнения к разделу 4.

1. Найти автомодельное решение нелинейного уравнения теплопроводности (4.28).

2. Найти решение нелинейного уравнения теплопроводности (4.28), основанное на инвариантном преобразовании вида

3. Найти автомодельное решение нелинейного уравнения теплопроводности с источником

4. Найти точные решения нелинейного уравнения теплопроводности

5. Найти автомодельное решение потенциального уравнение Бюргерса

6. Найти автомодельное и другие точные решения уравнения нелинейной фильтрации

7. Найти автомодельное и другие точные решения нелинейного волнового уравнения 

8. Найти точные решения нелинейного волнового уравнения 

9. Найти точные решения уравнения околозвукового газового потока

основанные на поиске инвариантных комбинаций преобразований сдвига и масштабирования.

10. Найти точное решение уравнения 

11. Найти точные решения уравнения Монжа — Ампера

12. Найти точные решения уравнения Буссинеска:

5. Общие выводы и замечания

Описанный метод с успехом может использоваться для решения самых различных типов математических уравнений. В случае обыкновенных дифференциальных уравнений и уравнений математической физики данный метод намного проще дл понимания, чем классический метод группового анализа [3–5, 7, 10–14], и может трактоваться как его частный случай. При использовании описанного простого метода надо уметь решать лишь самые простые алгебраические уравнения (и системы) и уметь дифференцировать. Важно отметить, что эффективность данного метода для обыкновенных дифференциальных уравнений практически не уступает классическому методу группового анализа [см. 4, 5], а для нелинейных уравнений математической физики позволяет находить все наиболее распространенные решения (которые по самым грубым оценкам составляют более половины решений, которые можно найти с помощью метода группового анализа; эту оценку можно получить путем исследования представленных в [9–11] результатов). Указанные обстоятельства говорят о целесообразности введения данного метода в стандартные курсы лекций по обыкновенным дифференциальным уравнениям и уравнениям математической физики, которые читаются студентам, специализирующимся в области прикладной математики, физики и механики (этот метод с успехом можно включать также в соответствующие специальные курсы, читаемые в некоторых технических и педагогических вузах). Кроме того, этот простой метод полезно излагать на предварительном этапе (в качестве введения) до чтения основного материала в существенно более сложных для восприятия спецкурсах по методам группового анализа дифференциальных уравнений.

Литература

1. Болтянский В. Г., Виленкин Н. Я. Симметрия в алгебре, 2-е изд. М.: Наука, 2002.
2. Кудряшов Н. А. Симметрия алгебраических и дифференциальных уравнений. Соросовский образовательный журнал, 1998, N 9, с. 104–110.
3. Ибрагимов Н. Х. Азбука группового анализа. М.: Знание, 1989, N 8.
4. Ибрагимов Н. Х. Опыт группового анализа обыкновенных дифференциальных уравнений. М.: Знание, 1991, N 7.
5. Polyanin A. D., Zaitsev V. F. Handbook of Exact Solutions for Ordinary Differential Equations, 2nd Edition. Boca Raton: Chapman & Hall/CRC Press, 2003.
6. Баренблатт Г. И. Подобие, автомодельность, промежуточная асимптотика. М.: Гидрометеоиздат, 1978.
7. Полянин А. Д., Зайцев В. Ф., Журов А. И. Методы решения нелинейных уравнений математической физики и механики. М.: Физматлит, 2005.
8. Седов Л. И. Методы подобия и размерности в механике. М.: Наука, 1972.
9. Полянин А. Д., Зайцев В. Ф. Справочник по нелинейным уравнениям математической физики: Точные решения. М.: Физматлит, 2002.
10. Ibragimov N. H. (ed.). CRC Handbook of Lie Group to Differential Equations, Vol. 1. Symmetries, Exact Solutions and Conservation Laws. Boca Raton: CRC Press, 1994.
11. Polyanin A. D., Zaitsev V. F. Handbook of Nonlinear Partial Differential Equations. Boca Raton: Chapman & Hall/CRC Press, 2004.
12. Овсянников Л. В. Групповой анализ дифференциальных уравнений. М.: Наука, 1978.
13. Дородницын В. А., Еленин Г. Г. Симметрия в решениях уравнений математической физики. М.: Знание, 1984, N 4.
14. Polyanin A. D., Manzhirov, A. V. Handbook of Mathematics for Engineers and Scientists. Boca Raton: Chapman & Hall/CRC Press, 2007.

Адрес статьи в Интернете: http://eqworld.ipmnet.ru/ru/methods/lpde/Invarianty/
Статья выложена на сайт: 19.03.2010 г.