Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

emph_f

.pdf
Скачиваний:
77
Добавлен:
10.03.2016
Размер:
3.07 Mб
Скачать

рой поверхности S, описываемой уравнением

 

S = fx 2 Rn : (x) = 0; r (x) 6= 0g:

(4.5)

В окрестности поверхности S введем новые координаты 1; 2; ::; n, полагая

i =

i(x1; :::; xn); i = 1; 2; :::; n;

(4.6)

где 1 = , а функции

2; :::; n выберем так, чтобы якобиан преобразо-

вания (4.6) был отличен от нуля на S. Производные от u по xi выразятся через производные по i по следующим формулам:

 

@u

n

@u @ l

 

@2u

n

@2u @ l @ k

 

n

@u @2

l

 

 

 

=

X

 

 

 

;

 

=

X

 

 

 

 

 

+

Xl

 

 

 

 

:

 

 

 

@ l @xi

 

 

@ l@ k @xi @xj

 

@ l @xi@xj

@xi

l=1

 

@xi@xj

l;k=1

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя эти выражения в (4.1), приходим к следующему уравнению для u в новых переменных 1; :::; n:

a(x; r (x))

@2u(x)

= G x; u(x); ru(x);

@2u(x)

; i + j > 2: (4.7)

@ 12

 

@ i@ j

 

Здесь G – известная функция своих аргументов, зависящая от x, u, ru и всех вторых производных @2u=@ i@ j, кроме производной @2u=@ 12, а коэффициент a в (4.7) определяется формулой

 

n

@ (x) @ (x)

 

 

a(x; r (x)) =

X

 

(4.8)

aij(x)

 

 

 

:

i;j=1

@xi

 

@xj

 

 

 

 

 

 

 

Из (4.5) и (4.6) вытекает, что условия (4.2), (4.3) при l = n переходят при преобразовании (4.6) в условия Коши

 

@u

 

uj 1=0 = '~0( 2; :::; n);

 

j 1=0

= '~1( 2; :::; n)

@ 1

для уравнения (4.7), заданные на гиперплоскости 1 = 0 и имеющие рассмотренный выше специальный вид. Другими словами, замена (4.6) как бы “выпрямляет” поверхность S (в окрестности точки x0 2 S), переводя ее в участок гиперплоскости 1 = 0. В таком случае из предыдущего анализа следует, что для того, чтобы условия Коши на поверхности S приводили к несовместности или неопределенности при нахождении производных второго порядка на S, необходимо и достаточно, чтобы функция в (4.5) удовлетворяла условию

 

n

@ (x) @ (x)

 

 

a(x; r (x)) =

X

= 0 на S:

(4.9)

aij(x)

 

 

 

 

i;j=1

@xi

 

@xj

 

 

 

 

 

 

 

 

 

151

Определение 4.1. Поверхность S в Rn, описываемая уравнением (4.5), называется характеристической для уравнения (4.1), а соответствующее направление нормали r (x0) называется характеристическим, в точке x0 2 S, если

a(x0; r (x0)) = 0;

(4.10)

и нехарактеристической, если

a(x0; r (x0)) 6= 0:

(4.11)

Поверхность S называется характеристической поверхностью (или характеристикой) для уравнения (4.1), если

a(x; r (x)) = 0 8x 2 S:

(4.12)

Если при этом для некоторого j n (@=@xj) a(x; r (x)) 6= 0 8x 2 S, то характеристика называется простой. В противном случае она называется кратной. Поверхность S, не имеющая характеристического направления нормали ни в одной точке x 2 S, называется нехарактеристической или свободной поверхностью для уравнения (4.1).

Свободная поверхность S называется пространственно ориентирован-

ной (или пространственно подобной) поверхностью, если

 

a(x; r (x)) > 0 8x 2 S;

(4.13)

и временным образом ориентированной (или временно подобной) поверхностью, если

a(x; r (x)) < 0 8x 2 S:

(4.14)

В соответствии с определением 4.1 рассмотрим два случая.

1. Поверхность S является нехарактеристической в точке x0. В этом случае в окрестности точки x0 уравнение (4.7) можно записать в так называемой нормальной форме:

@2u(x)

= x; u(x); ru(x);

@2u(x)

; i + j > 2:

(4.15)

@ 12

 

@ i@ j

 

Здесь = G=a, а функция G определена в (4.7). Основываясь на этом, можно доказать существование и единственность аналитического решения задачи Коши в окрестности точки x0. Этот результат составляет содержание теоремы Ковалевской. В применении к задаче (4.1)–(4.3) она имеет вид (ср. с ее формулировкой для нормальных систем в §2.1): если S –

аналитическая поверхность в окрестности своей точки x0, коэффициенты aij, правая часть f уравнения (4.1) и данные Коши '0 и '1 в (4.2), (4.3) аналитичны в этой окрестности, и выполнено условие (4.11), то в

152

некоторой окрестности точки x0 существует аналитическое решение u задачи Коши (4.1)–(4.3), которое в классе аналитических функций единственно.

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

2. Поверхность S является характеристической в точке x0. В этом случае уравнение (4.7) в самой точке x0 сводится к соотношению

a(x0; r (x0))

@2u(x0)

= 0 = G x0; u(x0); ru(x0);

@2u(x0)

; i + j > 2:

@ 12

 

@ i@ j

 

 

 

 

 

(4.16)

Указанное соотношение представляет собой ограничение на данные Коши в (4.2), (4.3). При этом возможны два варианта:

а) данные Коши в (4.2), (4.3) таковы, что (при l = n) выполняется второе равенство в (4.16). Отсюда следует, что исходное уравнение (4.1) выполняется в точке x0, причем в качестве второй производной @2u=@ 12 в этой точке можно взять любое число. Последнее означает, что имеет место неединственность решения. Следовательно, условий Коши (4.2), (4.3) недостаточно для однозначного определения решения в точке x0;

б) данные Коши в (4.2), (4.3) (при l = n) не удовлетворяют второму равенству в (4.16). В этом случае уравнение (4.1) не может быть удовлетворено в точке x0, так что условия Коши (4.2), (4.3) переопределяют рассматриваемую задачу Коши (4.1)–(4.3).

Уравнение (4.9) называется уравнением характеристик или характеристическим уравнением для уравнения (4.1). Подчеркнем, что хотя (4.9) формально имеет вид уравнения в частных производных первого порядка относительно , оно по своему определению им не является. Действительно, уравнение (4.9) должно выполняться не тождественно относительно x1; :::; xn, а только для точек x = (x1; x2; :::; xn), удовлетворяющих условию

(x) = 0, т.е. лежащих на искомой характеристической поверхности. Потребуем теперь дополнительно, чтобы уравнение (4.9) выполнялось

не только на поверхности S (т. е. при (x) = 0), но и тождественно относительно x1; :::; xn. Тогда получим нелинейное уравнение в частных производных 1–го порядка

n

 

 

 

 

 

X

@ @

 

(4.17)

aij(x)

 

 

 

= 0;

i;j=1

@xi @xj

 

 

 

 

 

 

 

153

а из его вида вытекает, что всякое решение

уравнения (4.16), отличное

от константы, будет давать целое семейство характеристик

(x1; :::; xn) = C;

(4.18)

где C – произвольная константа.

Замечание 4.1. Справедливо и обратное, а именно: для того чтобы соотношение (4.18) определяло семейство характеристик, зависящих от произвольной постоянной C, необходимо и достаточно, чтобы функция удовлетворяла уравнению (4.17). Более того, можно показать [21, c. 111], что всякую характеристику уравнения (4.1) можно включить в семейство вида (4.18). Отсюда следует, что все возможные решения уравнения (4.17) определяют все характеристические поверхности.

2.4.2. Примеры нахождения характеристик. Займемся нахождением характеристических поверхностей для трех основных типов уравнений математической физики и некоторых их нелинейных аналогов.

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

 

 

@2u

a2 u

@2u

a2

 

@2u @2u

= f(x; y; t)

(4.19)

 

 

 

 

 

+

 

 

@t2

@t2

@x2

@y2

на плоскости R2 изменения переменных x; y и вещественной оси R1t изменения времени t. Это уравнение можно записать в виде (4.1), если положить

x1 = x; x2 = y; x3 = t; a11 = a22 = a2; a33 = 1; aij = 0; при i 6= j:

(4.20) Мы уже знаем, что волновое уравнение (4.19) имеет гиперболический тип в каждой точке (x; y; t) 2 R2 R1t . Найдем теперь его характеристики.

Подставляя (4.20) в (4.17), приходим к следующему характеристическому уравнению для уравнения (4.19):

a @x ; @y ; @t

@t

 

2

a2

" @x

 

2

+ @y

#

= 0 на S: (4.21)

@ @ @

 

@

 

 

 

@

 

 

@

 

 

 

Пусть (x0; y0; t0) 2 R2 R1t – произвольная точка. Рассмотрим поверхность

S = (x; y; t) : (x; y; t) a2(t t0)2 (x x0)2 (y y0)2 = 0

;

(4.22)

 

0

0

 

0

 

являющуюся границей двух конусов с вершинами в точке (x

; y

; t ): верх-

него

 

 

 

 

 

 

 

 

 

K+(x0; y0; t0) = n(x; y; t) : a(t t0) >

 

 

 

o

(4.23)

 

(x x0)2 + (y y0)2

и нижнего

p

 

 

 

 

 

 

K (x0; y0; t0) = n(x; y; t) : a(t t0) > p

 

o;

(4.24)

(x x0)2 + (y y0)2

154

a

б

Рис. 4.2

называемых соответственно конусами будущего и прошлого (см. рис. 4.2,a). Легко проверить, что функция в (4.22) удовлетворяет уравнению (4.21)

в любой точке (x; y; t), принадлежащей поверхности S +(x0; y0; t0) [(x0; y0; t0), где (x0; y0; t0) = @K (x0; y0; t0). Действительно, так как

@

= 2a2(t t0);

@

= 2(x x0);

@

= 2(y y0);

(4.25)

 

 

 

 

@t

 

@x

@y

то

a

@

;

@

;

@

= 4a2 a2(t t0)2 (x x0)2 (y y0)2 : (4.26)

 

 

 

 

@x

@y

 

@t

Отсюда и (4.22) следует, что

a

@

;

@

;

@

0 на S +(x0; y0; t0) [ (x0; y0; t0):

(4.27)

 

 

 

 

 

@x

@y

@t

Это означает, что поверхность S, определяемая уравнением (4.22), для любой точки (x0; y0; t0) является характеристической. Ниже будем ссылаться на нее как на характеристический конус с вершиной в точке (x0; y0; t0), причем верхнюю ее часть +(x0; y0; t0) будем называть характеристическим конусом будущего, а нижнюю часть (x0; y0; t0) – характеристическим конусом прошлого.

Отметим, что поверхность S, описываемая уравнением

S = (x; y; z) : (x; y; t) a2(t t0)2 (x x0)2 (y y0)2 = const 6= 0 ;

155

не удовлетворяет (4.21). Следовательно, она не является характеристической.

Аналогичные выводы справедливы и для n-мерного волнового уравнения

@2u

@2u

n

@2u

 

 

 

 

 

Xi

 

= f(x; t); x = (x1; x2; :::; xn); (4.28)

@t2

a2 u @t2 a2

 

 

@x2

 

 

 

 

=1

i

 

 

 

 

 

 

 

рассматриваемого в Rnx R1t , для которого отвечающее ему характеристическое уравнение имеет вид:

 

 

a

r ;

@t

 

@t

2

n

@xi

2

= 0 на S:

(4.29)

 

 

a2

i=1

 

 

 

 

 

@

 

@

 

 

X

@

 

 

 

 

 

Рассмотрим поверхность вида

 

 

 

 

 

 

 

 

 

 

 

 

0 S

0= (x;0t) :

(x; t) a2(t t0)2 jx x0j2 = 0 ;

(4.30)

где x0 = (x1; x2; :::; xn), представляющую собой границу двух конусов (точ-

нее, гиперконусов) с вершинами в точке (x0; t0): верхнего

 

 

 

 

 

K+(x0; t0) = f(x; t) : a(t t0) > jx x0jg

(4.31)

и нижнего

K (x0; t0) = f(x; t) : a(t t0) > jx x0jg ;

(4.32)

 

 

 

 

называемых

соответственно

конусами

будущего

 

и прошлого.

Здесь

n

1

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 2

jx x0j = (x1 x10)2 + (x2 x20)2 + + (xn xn0 )2.

 

 

Рассуждая, как и выше, показываем, что для любой точки (x ; t )

R Rt

функция

в (4.30) удовлетворяет уравнению (4.29). Это означа-

ет, что для любой точки (x0; t0) поверхность S, определяемая уравнением (4.30), является характеристической поверхностью. Ее называют характеристическим конусом с вершиной в точке (x0; t0), причем верхнюю ее часть +(x0; t0) называют характеристическим конусом будущего, а нижнюю часть (x0; t0) – характеристическим конусом прошлого. Кроме построенного семейства характеристических поверхностей (конусов) уравнение (4.28) имеет и другое семейство характеристических поверхностей, а именно: семейство касательных плоскостей к характеристическим конусам (4.31) и (4.32). Уравнение второго семейства характеристик имеет вид

n

 

Xi

(4.33)

at + (b; x) at + bixi = C:

=1

 

Здесь b = (b1; b2; :::; bn), bi и C – любые вещественные числа, причем jbj = 1.

156

В частном случае n = 1 (4.28) переходит в одномерное волновое уравнение

 

 

 

 

@2u

a2

@2u

= f(x; t)

 

 

(4.34)

 

 

 

 

@t2

@x2

 

 

 

на плоскости R1 Rt1, а его характеристическое уравнение имеет вид

a @x

; @t

@t

2

 

 

 

 

2

(4.35)

 

a2 @x

= 0:

@

 

@

 

 

 

@

 

 

 

 

 

@

 

 

 

Записав указанное уравнение в виде (@

=@t+a@

=@x)(@ =@t a@ =@x) =

0, приходим к выводу, что функция

 

является решением уравнения (4.35)

тогда, когда удовлетворяет одному из уравнений:

 

 

@

 

@

 

 

 

 

и

@

@

 

 

(4.36)

 

 

 

+ a

 

 

= 0

 

 

 

 

 

a

 

= 0:

 

@t

@x

 

 

 

@t

 

@x

Из результатов п. 2.2.1 следует, что решение первого (либо второго) уравнения в (4.36) сохраняет постоянное значение на характеристиках

1(x; t) x at C1 = 0 либо 2(x; t) x + at C2 = 0; (4.37)

описывающих соответствующие прямые на плоскости x; t (см. рис. 2.4,б). Здесь C1 и C2 – произвольные константы. Отсюда и замечания 4.1 приходим к выводу, что указанные прямые (4.37) и представляют собой два искомых семейства характеристик уравнения (4.34). При этом через любую точку (x0; t0) 2 R1x R1t проходят две характеристики, определяемые уравнениями x at = x0 at0 и x + at = x0 + at0. Указанные прямые яв-

ляются границами двумерного конуса K(x0; t0) = K+(x0; t0) [ K (x0; t0), изображенного на рис. 4.2б.

Замечание 4.2. Все, что говорилось выше о характеристиках для n–мерного волнового уравнения (4.28), остается справедливым и для n– мерного квазилинейного уравнения вида

@2u

 

 

n

@2u

 

r

@u

 

 

@t2

 

i=1

@xi

@t

 

 

 

a2

X

 

 

 

 

(4.38)

 

 

 

2

= F

x; t; u;

u;

 

:

Здесь F – заданная функция своих аргументов. Другими словами, уравнение (4.38) так же, как и (4.28), имеет гиперболический тип в каждой точке (x; t) 2 Rn Rnt , а его характеристиками являются характеристические конусы +(x0; t0) и (x0; t0) так же, как и гиперплоскости (4.33).

Однако в случае более общего квазилинейного уравнения вида

@2u

n

 

@2u

@u

 

X

 

(4.39)

 

 

aij(x; t)

 

= F (x; t; u; ru;

 

);

@t2

i;j=1

@xi@xj

@t

 

 

 

 

 

 

 

 

157

гиперболического при выполнении (1.23), его характеристическое уравнение имеет вид нелинейного уравнения с переменными коэффициентами

@

@

 

2

n

@ @

 

 

 

 

X

 

 

 

 

; r )

 

 

 

 

 

 

 

(4.40)

a(x; t; ;

 

 

 

i;j=1 aij(x; t)

 

 

 

= 0:

@t

@t

 

@xi

@xj

Поэтому нахождение характеристик уравнения (4.39) в явном виде возможно лишь в исключительных случаях, хотя можно строго показать, что через каждую точку (x0; t0) проходят две характеристики уравнения (4.39).

То же самое справедливо и для одномерного аналога

@2u

a(x; t)

@2u

 

@u

 

@u

(4.41)

 

 

= F (x; t; u;

 

;

 

)

@t2

@x2

@x

@t

уравнения (4.39). Условием его гиперболичности в точке (x; t) (либо в обла-

сти D R Rt) является условие a(x; t) > 0 (либо a(x; t) > 0 8(x; t) 2 D), а характеристическое уравнение имеет вид

@t

2

 

@x

 

2

(4.42)

a2(x; t)

= 0:

@

 

 

 

@

 

 

 

Записав его в виде [@ =@t + a(x; t)@

=@x][@

=@t a(x; t)@

=@x] = 0, при-

ходим с учетом замечания 4.1 к выводу, что уравнение (4.41) имеет два

семейства характеристик

1(x; t) = C1 и

2(x; t) = C2. Здесь C1

и C2

произвольные константы, а функции

1 и

2 являются решениями следу-

ющих двух уравнений в частных производных 1-го порядка:

 

@ 1

 

@ 1

= 0 и

@ 2

 

@ 2

 

(4.43)

 

 

+ a(x; t)

 

 

 

a(x; t)

 

= 0:

 

@t

@x

@t

@x

Из результатов п. 2.2.4 и 2.2.5 вытекает, что нахождение решений каждого из уравнений 1-го порядка в (4.43) сводится к нахождению характеристик данного уравнения. Они определяются как интегральные кривые следующих задач Коши для обыкновенных дифференциальных уравнений:

dx

 

dx

(4.44)

 

 

= a(x; t);

xjt=0 = x0;

 

= a(x; t); xjt=t0 = x0:

 

dt

dt

При этом функция

i является решением i-го уравнения в (4.43) тогда и

только тогда, когда соотношение

i(x; t) = Ci определяет первый интеграл

соответствующего уравнения в (4.44), i = 1; 2.

 

Из теории обыкновенных дифференциальных уравнений известно [51], что в случае, когда функция a непрерывна по t и дифференцируема по x, через каждую точку (x0; t0) рассматриваемой области проходит единственная интегральная кривая первой задачи Коши в (4.44), определяющая

158

характеристику первого уравнения в (4.43), и единственная интегральная кривая второй задачи Коши в (4.44), определяющая характеристику второго уравнения в (4.43). Указанные характеристики по построению образуют два искомых семейства характеристик для исходного квазилинейного уравнения (4.41).

Пример 2. Рассмотрим n–мерное уравнение теплопроводности

@u

n

@2u

 

 

 

 

Xi

 

 

(4.45)

@t a2

 

= f

 

@x2

 

 

=1

i

 

 

 

 

 

 

 

в пространстве Rn R1t . Его характеристическое уравнение имеет вид

n

@

2

 

X

 

 

 

 

= 0:

(4.46)

i=1

@xi

Поскольку решениями уравнения (4.46) в Rn R1t являются функции (x; t) = t const, то единственными характеристиками уравнения теп-

лопроводности (4.45) являются гиперплоскости t = const.

Аналогичные утверждения справедливы и для n–мерного квазилинейного параболического уравнения вида

@u

 

 

n

@2u

 

r

@u

 

 

@t

 

i=1

@xi

@t

 

 

 

a2

X

 

 

 

 

(4.47)

 

 

 

2

= F

x; t; u;

u;

 

;

а также более общего квазилинейного уравнения вида

@u

n

 

@2u

@u

 

X

 

(4.48)

 

 

aij(x; t)

 

= F (x; t; u; ru;

 

);

@t

i;j=1

@xi@xj

@t

 

 

 

 

 

 

 

 

параболического при выполнении условия (1.23). Действительно, его характеристическое уравнение имеет вид

n

 

 

 

 

 

X

@ @

 

(4.49)

aij(x; t)

 

 

 

= 0:

i;j=1

@xi @xj

 

 

 

 

 

 

 

Ясно, что при выполнении (1.23) единственным решением уравнения (4.49) является функция (x; t) = t const. Отсюда следует, что уравнение (4.48), как и (4.47), при выполнении условия (1.23) имеют своими характеристиками гиперплоскости t = const.

Пример 3. Рассмотрим в пространстве Rn уравнение Пуассона

X

n @2u

u = @x2 = f: (4.50)

i=1 i

159

Его характеристическое уравнение по-прежнему имеет вид (4.46), однако, в отличие от уравнения теплопроводности (4.45), его следует рассматривать в пространстве Rn. С учетом этого приходим к выводу, что уравнение (4.50) не имеет вещественных характеристик. Подчеркнем, что речь идет о вещественных характеристиках уравнения (4.50), определяемых через вещественные решения уравнения (4.46). Можно показать, что в поле комплексных чисел уравнение (4.46) имеет два семейства “мнимых” решений, которым отвечают два семейства мнимых характеристик уравнения (4.50).

Аналогичный факт справедлив и для уравнения Гельмгольца

u + k2u = f;

рассматриваемого в пространстве Rn, а также для произвольного квазилинейного эллиптического уравнения u = f(x; u; ru), главная часть которого совпадает с оператором Лапласа, и, кроме того, для более общего эллиптического уравнения вида

n

@2u

 

X

= f(x; u; ru)

aij(x)@xi@xj

i;j=1

 

 

при условии, что коэффициенты aij(x) удовлетворяют условию эллиптичности (1.19).

160

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]