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

1

§13 Краевые задачи для уравнений Лапласа и Пуассона. Задача Дирихле в шаре.

Сформулируем теперь некоторые краевые задачи для уравнений Лапласа и Пуассона (заметим, что в x7 мы уже сформулировали так называемую задачу Дирихле в круге (n = 2)).

Определение.

Первой краевой задачей (или задачей Дирихле) для уравнения Пуассона (или Лапласа) называется задача о нахождении функции u(x) 2 C2(-) \ C(-), удовлетворяющей следующим условиям:

4xu = f(x)(= 0); x

-;

(1)

½ uj@- = '(x); x 2 @-;2

 

 

f, ' - заданные функции.

Третьей краевой задачей для уравнения Пуассона (или Лапласа) называется задача о нахождении функции u(x) 2 C2(-) \ C1(-), удовлетворяющей следующим условиям:

½

¡

xu = f(x)(= 0); x

2

-;

 

 

@

¢¯

 

(2)

4@u

 

 

 

 

 

@N

+ ¾(x)u @- = '(x); x 2 @-;

 

f, ', ¾ - заданные функции;

 

= (¯N; r), N - единичная внешняя нормаль к

@N

@-.

 

 

 

 

 

 

 

 

Если ¾(x) ´ 0, то задача (2) называется второй краевой задачей (или задачей Неймана) для уравнения Пуассона (или Лапласа).

Замечание.

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

б) При n = 2 задачи (1), (2) имеют четкий физический смысл: например, они описывают распределение прогиба упругой мембраны (при том или ином способе крепления края мембраны).

в) Давая определение краевых задач, мы полагали, что - ½ Rn - ограниченная область. Задачи (1), (2) при этом еще иногда называют внутренними (interior): внутренняя задача Дирихле (Di), внутренняя задача Неймана (Ni) и т.д. Если область - ½ Rn - неограниченная и является внешностью ограниченного множества, то задачи (1), (2) называют внешними: внешняя задача Дирихле (De), внешняя задача Неймана (Ne) и т.д.

Замечание о задаче Коши.

Рассмотрим сначала вопрос о единственности классического решения внутренних задач (1), (2).

Теорема 1.

Задача Дирихле (1) и задача (2) при ¾(x) ¸ 0, x 2 @- не могут иметь более одного классического решения. Классическое решение задачи Неймана определяется с точностью до постоянной.

Доказательство.

Рассмотрим задачу Di. Пусть u1;2(x) - два решения задачи (1), отвечающие одним и тем же функциям f, '. Тогда их разность u(x) = u1(x)¡u2(x) является гармонической функцией в -, непрерывной в - и uj@- = 0. Поэтому, в силу принципа максимума (вернее, в силу следствия из Теоремы 12 предыдущего x12)

jjujjC(-) · jjujjC(@-) = 0; т.е. u ´ 0 в -:

Значит в u1(x) ´ u2(x) области -.

Рассмотрим задачу (2i). Пусть u1;2(x) - два решения задачи (2i), отвечающие одним и тем же функциям f, ', ¾. Тогда их разность u(x) = u1(x) ¡u2(x)

 

 

 

 

2

удовлетворяет следующей задаче:

 

 

 

½

xu(x) = 0; x -;

i

¡@N + ¾(»)u¢¯@-

= 0; » 2 @-;

 

4@u

¯

2

(20)

 

 

 

 

Используем теперь первую формулу Грина (3) из x12 (убедитесь в том, что она применима):

0 = Z-

u ¢ 4xu ¢ dx =@Z-

u ¢ @N dS ¡ Z-

jruj2dx

 

 

 

 

 

 

@u

 

 

или

Z

jruj2dx + Z

 

 

 

 

 

 

 

¾(x)u2(x)dSx = 0:

(¤)

-@-

Из (¤) следует, что при ¾(x) ¸ 0 u(x) ´ 0 в -. Для задачи Ni из (¤) следует, что u(x) ´ const в -.

Замечание.

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

формулу (2) из x12):

Z

4udx = Z

 

@u

 

 

 

 

dS

 

 

 

@N

 

и

-

 

@-

 

 

 

 

 

Z

Z

 

 

 

 

 

 

 

 

 

(3)

 

 

fdx =

 

'dS:

-@-

Перейдем к вопросу о существовании классического решения задачи Di. Сформулируем сначала общий результат.

Теорема 2.

Если:

а) f 2 C1(-), б) ' 2 C(@-),

в) граница @- достаточно гладкая,

то задача Di имеет классическое решение.

Замечание.

На первый взгляд кажется, что условие а) завышено (не достаточно ли для разрешимости задачи Di выполнения условия f 2 C(-)?). Однако заменить а) на это условие нельзя (можно построить пример, показывающий, что задача Di не имеет классического решения, если f 2 C(-)).

Методы доказательства теоремы 2: метод функций Грина, метод потенциалов. Рассмотрим эти методы. Остановимся сначала на методе функций Грина.

Предварительные сведения.

Пусть u 2 C2(-), n ¸ 2. Тогда в силу теоремы 1 из x12 (см. формулы (7), (9)):

u(x) = u0(x) + u2(x) ¡ u1(x); x 2 -;

(+)

 

 

 

(++)

0 = u0(z) + u2(z) ¡ u1(z); z 2= -:

Умножим (++) на произвольную непрерывную функцию d(z), z 2= - и вычтем его из (+). В итоге мы получим еще одно представление произвольной функции u 2 C2(-):

u(x) = fU(x ¡ y) ¡ d(z)U(z ¡ y)g4yu(y)dy+

 

 

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

@u(»)

 

 

(4)

 

 

+ f[d(z)U(z ¡ ») ¡ U(x ¡ »)]

@N»

+

 

 

 

-

 

 

 

 

 

 

 

 

 

 

 

 

+

@

[U@(Rx

¡

»)

¡

d(z)U(z

¡

»)]u(»)

dS

; x

2

-:

 

@N»

 

 

 

g

 

»

 

 

 

3

Оказывается, что справедлив следующий факт: при весьма общих предположениях относительно области - существует такое отображение z = z(x), переводящее т.x 2 - в точку z 2= - и такая функция d(z(x)), что для всех x 2 - выполнено тождество

d(z(x))U(z(x) ¡ ») ¡ U(x ¡ ») ´ 0; » 2 @-:

(5)

Пусть мы нашли такие функции z(x), d(z(x)). Тогда из (4) следует еще одно представление для произвольной функции u 2 C2(-):

u(x) = fU(x ¡ y) ¡ d(z(x))U(z(x) ¡ y)g4yu(y)dy+

 

-

 

 

 

 

 

 

 

 

 

 

 

 

(6)

+

u(»)

R@

 

U(x

 

»)

 

d(z(x))U(z(x)

 

»)

dS

; x

 

-:

@N» f

¡

¡

¡

2

 

@R-

 

 

 

 

g

»

 

 

 

Дальнейшие рассуждения будем проводить для случая, когда область - - шар: SR;0. Простые рассуждения позволят нам найти функции z(x), d(z(x)) с нужными свойствами. В самом деле, пусть n > 2. Тогда условие (5) перепишется

так:

jz(x) ¡ »j2 ´ jx ¡ »j2d(z(x)); j»j = R

или

j

 

¡

 

2

´ j

 

¡

 

2b2

1

 

j

 

j

 

z(x)

»

x

»

(z(x)); b(z(x)) = [d(z(x))]

2

;

»

= R:

 

j

 

 

j

 

 

 

 

 

 

Функцию z(x) будем искать в следующем виде

z = a(x)x;

где a(x) - неизвестная пока функция. Тогда

jz(x) ¡ »j2 = (z ¡ »; z ¡ ») = jz(x)j2 ¡ 2(z(x); ») + R2 ´ ´ b2(z(x))jx ¡ »j2 = b2(z(x))fjxj2 ¡ 2(x; ») + R2g

или

 

 

 

 

 

[a2(x) ¡ b2(y(x))]jxj2 + (1 ¡ b2(z(x)))R2 ´

 

 

 

 

 

 

 

 

 

 

 

 

 

(

¤¤

)

Положим

 

 

 

 

 

 

´ 2[a(x) ¡ b2(z(x))](x; »);

j»j = R:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a(x) = b2(z(x)); b(z(x)) =

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jxj

 

 

 

 

 

При таком

выборе тождество ( ) выполняется. Далее, отображение

 

 

 

 

 

R2

 

 

 

¤¤

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = a(x)x =

 

 

 

x переводит т.x

2

SR;0 (

x

j

< R) в т.z =

SR;0

( z

j

> R),

jxj

2

 

 

2

2

 

 

 

j

 

 

2

 

j

 

 

 

поскольку jzj =

R

jxj =

R

> R при jxj < R. Кроме этого:

 

 

 

 

 

jxj2

jxj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d(z(x)) = µjxj

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

2

 

 

 

 

 

С учетом найденных функций z(x), d(z(x)) преобразуем формулу (6), переписав ее сначала для шара SR;0.

Поскольку:

1) U(x ¡ y) ¡ d(z(x))U(z(x) ¡ y)n=2

 

 

1

1

 

 

 

R

 

¡

(7)

 

 

 

 

jxj

= ¡

 

½

 

¡ µ

 

 

 

¶ ¾ = ¡GR(x; y)

 

(2)¾n

jx¡yj2

j

R2

x¡yj

 

jxj2

 

2)

»

j

= R;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@ j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fU(x ¡ ») ¡ d(z(x))U(z(x) ¡ »)g = (N»; r»f: : :g) =

@N»

=

 

1

(»;

»

 

: : : ) =

 

 

 

 

1

 

 

 

1

(»;

 

 

 

 

 

 

 

1

)+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2

 

 

 

R

 

1

 

 

r1f

 

R gn 2

 

 

¡(2)¾n R

1

r»

©jx¡»j

¡

ª

+

 

 

 

 

 

 

 

 

 

¢

¡

(»; r»

 

 

 

 

 

 

 

 

 

 

 

ª

) =

 

 

(n 2)¾n

R

x

 

 

 

R2

x »

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

¡

 

 

 

 

 

¡j j

 

 

 

 

R n¡2©jjxj

¡ j

R2

 

 

 

 

 

 

1

 

 

(»;»¡x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jxj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(»; »

 

 

 

 

 

 

x) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

R2

 

 

 

 

 

 

2

 

 

 

 

 

n

½ jx¡»j

 

 

¡ j¡jxj2

x¢¡»jn

 

 

 

 

 

 

 

¡ jxj

 

 

 

 

¾

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2¡

R2

 

(»;x)

 

 

 

 

=

 

1

 

 

1

 

 

 

R

2

 

 

 

(»; x)

 

 

 

jxj2

 

=

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

 

 

 

n jx¡»j

 

©

 

 

 

¡

 

 

 

 

 

¡

 

 

 

 

 

 

jxj2

 

 

 

 

ª

 

 

4

(8)

 

R2¡jxj2

 

 

=

njx¡»jn

= PR(x; »); j»j = R:

Задача.

 

 

 

 

Покажите, что

 

 

@

 

 

PR(x; ») = ¡

GR(x; »); j»j = R:

 

 

 

@N»

Итак, для любой функции u 2 C2(-) где - - шар SR;0, имеет место представление:

u(x) =

Z

u(»)PR(x; »)dS» ¡ Z

GR(x; y)4yu(y)dy; x 2 -

(9)

 

j»j=R

jyj<R

 

 

Замечание.

1) При n = 2 представление (9) также имеет место, при этом PR(x; ») имеет тот же вид (8), а функция GR(x; y) принимает следующий вид:

 

1

ln½

y x

2

 

R2x

 

¾

 

GR(x; y) =

 

j j j ¡

 

y

j

j

(70)

2¼

R x

 

 

 

 

 

 

j

¡

 

 

 

 

2) Функция PR(x; »), jxj < R, j»j = R называется ядром Пуассона задачи Di

шаре SR;0;

Функция GR(x; y), jxj; jyj < R называется функцией Грина задачи Di в шаре

SR;0.

3) Напомним, что функция Грина GR(x; y) дается формулой (7) при n > 2 (и формулой (7’) при n = 2):

GR(x; y) = (n ¡ 2)¾n ½jx ¡ yj2

¡

µjjRxj22 x ¡ yj¶ ¾:

1

1

 

 

 

 

R

2

 

 

 

 

 

 

 

 

jxj

 

 

 

 

 

 

 

 

 

 

 

Опишем более подробно эту функцию (имея в виду дальнейшие соображения). Эта функция определена в области f(x; y) 2 R2n; x =6 y; jRxj22 x 6= yg, непрерывна в ней и обладает следующими легко-проверяемыми свойствами:

а) GR(x; y) ´ 0 при jxj = R; б) GR(x; y) = GR(y; x);

в) 4x;yGR(x; y) = 0;

г) при jxj; jyj · R:

0 · GR(x; y) ·

 

 

1

 

; n > 2

 

 

 

 

 

¾njx ¡ yj2

или

1

 

2R

 

0 · GR(x; y) ·

ln

; n = 2:

 

 

 

2¼

jx ¡ yj

Свойство б) доказывается с помощью очевидного равенства (проверьте его справедливость!)

jR2y ¡ xjyj2j2 ¢ jxj2 = jR2x ¡ yjxj2j2 ¢ jyj2:

5

Из него вытекает, что условие jRxj22 x 6= y эквивалентно условию jRyj22 y 6= x и поэтому, если точка (x; y) принадлежит области определения функции GR(x; y) то точка (y; x) принадлежит той же области.

Свойство г) доказывается с помощью принципа максимума.

Задача.

Докажите свойство в).

Вернемся снова к представлению (9). Оказывается, такое же представление справедливо, если на функцию u(x), jxj < R наложены менее ограничительные требования. А именно, справедлива

Теорема 3.

Пусть u 2 C2(SR;0) \ C(SR;0), функция 4u(x) - непрерывна и ограничена в шаре SR;0. Тогда для любой точки jxj < R справедливо равенство (9).

Доказательство теоремы простое по идее, но громоздкое по технике. Поэтому я ограничусь изложением только идеи доказательства. Возьмем шар SR;0,

R < R

 

x

 

u

 

C2(S

 

)

 

шаре

S

 

 

, в котором лежит т.

 

(см. Рис.). Тогда

 

2

 

 

R;0

 

. Поэтому в

e

eR;0

справедливо представление (9). Записав это

 

 

e

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

e

R ! R, надо убедиться, что слагаемые в правой части формулы (9) имеют пределы для любой т.x: jxj < R. Закончим рассказ о предварительных сведе-

SR,0

x

S ~

R,0

0

ниях в этом параграфе ещё несколькими фактами. 1) Легко показать, что

4xPR(x; ») = 0; jxj < R; j»j = R:

(10)

2) Если уравнение 4xu(x) = 0 в шаре SR;0 имеет решение u 2 C2(SR;0), то это решение записывается в следующем виде:

Z

u(x) =

PR(x; »)u(»)dS»:

(11)

 

j»j=R

 

Формула (11) называется формулой Пуассона (при n = 2 формула Пуассона была получена нами в x7). В частности, задача Di следующего вида:

½ ujjxj=R = 1

 

4xu(x) = 0; jxj < R;

(Di)

имеет единственное решение (почему?) u(x) ´ 1. Тогда в силу (11) имеем:

Z

u(x) ´ 1 =

PR(x; »)dS»:

(12)

 

j»j=R

 

 

 

 

 

6

3) Пусть u 2 C2(

 

) финитная функция . Тогда из (9) получаем:

 

SR;0

 

 

 

Z

 

(13)

 

 

u(x) = ¡

GR(x; y)4yu(y)dy:

 

 

jyj<R

 

 

Из теоремы 3 следует, что если классическое решение задачи Di в шаре SR;0 существует (при непрерывной на сфере j»j = R функции '(») и ограниченной непрерывной в шаре jxj < R функции f(x)), то оно (это классическое решение)

представляется в виде

PR(x; »)'(»)dS» ¡ Z

 

 

u(x) =

Z

GR(x; y)f(y)dy:

(14)

 

j»j=R

jyj<R

 

 

Пусть теперь выполнены условия теоремы 2: а) f(x) 2 C1(SR;0),

б) '(») 2 C(@SR;0).

Покажем, что (14) и есть искомое решение задачи Di в шаре SR;0. Доказательство этого факта состоит из двух частей.

Часть 1.

Обозначим первое слагаемое в (14) через

Z

v(x) = PR(x; »)'(»)dS»; jxj < R:

j»j=R

Покажем, что

а) 4xv(x) = 0, jxj < R;

б) v(x) ! '(³), если x ! ³, j³j = R.

Свойство а) выполнимо в силу формулы (10). Докажем свойство б). Очевидно, что функция v(x) 2 C2(SR;0) по крайней мере. Составим разность (с учетом формулы (12)):

v(x) ¡ '(³) =

PR(x; »)f'(») ¡ '(³)gdS» =

 

 

=R

 

 

j»jR

 

= j»j=R;jR»¡³j<± PR(x; »)f'(») ¡ '(³)gdS»+

(¤)

+j»j=R;jR»¡³j¸± PR(x; »)f'(») ¡ '(³)gdS»; j³j = R;

 

± > 0 - некоторое число.

 

 

 

 

 

jv(x) ¡ '(³)j · j1 слагj + j2 слагj:

 

 

 

 

 

 

 

 

'

»)

C(@S

 

)

, то 8

" > 0

± > 0

 

»

¡

³

j

< ±

 

Поскольку

(

 

"

2

 

R;0

 

, 9

, что как только j

 

 

 

, то

j'(») ¡ '(³)j ·

2

и поэтому

Z

 

 

 

 

 

 

 

 

 

 

j1 слагj = j

 

 

PR(x; »)f'(») ¡ '(³)gdS»j ·

"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

j»j=R;j»¡³j<±

(в силу (12)).

Рассмотрим теперь второе слагаемое. Поскольку x ! ³, то можно считать,

что jx¡³j · 2± . Поэтому j»¡xj ¸ j»¡³j¡j³ ¡xj ¸ ±¡ 2± = 2± > 0 и следовательно второе слагаемое в (¤) стремится к нулю (?), когда x ! ³ и станет меньше 2"

(по абсолютной величине).

Итак, v(x) ! '(³), когда x ! ³, j³j = R, jxj = R.

Часть 2.

Обозначим второе слагаемое в (14) через

Z

w(x) = ¡ GR(x; y)f(y)dy; jxj < R:

jyj<R

7

x

Для функции w(x) справедливо следующее утверждение:

Если f(x) 2 C1(SR;0), то w(x) 2 C2(SR;0) \ C1(SR;0) и для всех x: jxj < R

4xw(x) = f(x);

(¤¤)

при этом w(x) ! 0, когда x ! ³, j³j = R.

Последнее свойство очевидно в силу того, что GR(x; y) = 0 при jxj = R. Первая часть утверждения доказывается путем громоздких выкладок (и мы их опустим). Докажем свойство (¤¤). Возьмем произвольную финитную в SR;0 функцию ª(x) 2 C2(SR;0) В силу (13):

Z

ª(x) = ¡ GR(x; y)4yª(y)dy; jxj < R:

jyj<R

К функциям ª(x), w(x) применим вторую формулу Грина (см. x12):

Z

fª(x)4xw(x) ¡ w(x)4xª(x)gdx =

jxj<R

=

Z ½ª(»)

@w(»)

¡ w(»)

@ª(»)

 

 

 

@N»

@N»

 

j»j=R

 

 

 

 

¾

dS» = 0;

т.е.

Z

 

 

Z

 

 

ª(x)4xw(x)dx =

w(x)4xª(x)dx =

 

jxj<R

4xª(x)½¡ Z

jxj<R

 

 

=

Z

GR(x; y)f(y)dy¾dx =

 

jxj<R

jyj<R

¾

 

 

Z

½ Z

 

 

=

 

f(y) ¡ GR(x; y)4xª(x)dx dy =

jyj<R jxj<R

ZZ

=

f(y)ª(y)dy =

ª(x)f(x)dx

jyj<R

jxj<R

8

или

Z

 

f4xw(x) ¡ f(x)gª(x)dx:

(+)

jxj<R

При получении (+) использовались: теорема Фубини, свойство симметрии функции Грина GR(x; y). Поскольку ª(x) - произвольная финитная функция из C2(SR;0), то

4xw(x) = f(x);

что и требовалось доказать.

Итак, решение задачи Di в шаре SR;0

½4xu(x) = f(x); f(x) 2 C1(SR;0); x 2 SR;0; ujjxj=R = '(x); '(x) 2 C(@SR;0); jxj = R

дается формулой (14).

Соседние файлы в папке 2011-02-12