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

u_practice

.pdf
Скачиваний:
63
Добавлен:
15.04.2015
Размер:
584.37 Кб
Скачать

Z

(f, g)I = [qfg + p(rf, rg)] dx

Ω

эквивалентны.

10.15. Пусть вещественные функции pij, pij(x) = pji(x), i, j = 1, 2, ..., m

 

q

 

C(Ω)

 

q(x)

 

0

 

x

 

Ω

 

вещественных

и

 

принадлежат

 

,

 

 

и для всех

 

 

 

и всехm

 

m

, ..., ξm)

Rm

 

 

 

 

 

P

векторов ξ = (ξ1

имеет место неравенство

pij(x)ξiξj

 

P

 

 

 

 

 

 

 

 

 

 

i,j=1

 

 

 

 

 

 

 

 

 

 

 

 

γ0

 

ξi2, где постоянная γ0 > 0. Доказать, что в H1 (Ω) можно определить

 

i=1

 

 

 

 

 

 

 

 

 

 

 

скалярное произведение

 

 

m pijfxi gxj + qfg! dx,

 

 

 

 

(f, g)I

=

 

 

 

 

 

 

 

Z

 

X

 

 

 

 

 

Ω i,j=1

эквивалентное скалярному произведению

Z

(f, g) = [fg + (rf, rg)] dx.

Ω

10.16. Пусть p, q C(Ω), p(x) ≥ p0 > 0, q(x) ≥ q0 > 0. Доказать, что скалярные произведения в H1(Ω)

Z

(f, g) = [fg + (rf, rg)] dx,

Ω

Z

(f, g)I = [qfg + p(rf, rg)] dx

Ω

эквивалентны.

10.17. Пусть p, q C(Ω), p(x) ≥ p0 > 0, q(x) ≥ 0 и q(x) 6≡0 в Ω. Доказать, что скалярные произведения в H1(Ω)

Z

(f, g) = [fg + (rf, rg)] dx,

Ω

Z

(f, g)I = [qfg + p(rf, rg)] dx

Ω

эквивалентны.

10.18. Доказать, что если σ(x) C(∂Ω) и σ(x) > 0, то выражение

Z Z

(f, g)I = (rf, rg) dx + σfg ds

Ω ∂Ω

61

задает скалярное произведение, причем оно будет эквивалентным скалярному произведению

Z

(f, g) = [fg + (rf, rg)] dx.

Ω

10.19. Пусть p, q C(Ω), σ(x) C(∂Ω), p(x) ≥ p0 > 0, q(x) ≥ 0 в Ω,

σ(x) ≥ 0 на ∂Ω или q(x) 6≡0, или σ(x) 6≡0 . Тогда скалярные произведения в H1(Ω)

Z

(f, g) = [fg + (rf, rg)] dx,

Ω

Z Z

(f, g)I = [qfg + p(rf, rg)] dx + σfg ds

Ω ∂Ω

эквивалентны.

11. Обобщенные решения уравнений эллиптического типа

Пусть Ω Em некоторая ограниченная область с гладкой границей

∂Ω.

В области Ω задано уравнение

−div(p(x)ru(x)) + q(x)u = f(x),

(11.1)

а на границе области граничные условия первого рода

 

 

u(x)|∂Ω = ϕ(x),

(11.2)

или третьего рода

∂n + σ(x)u ∂Ω = ϕ(x).

(11.3)

 

∂u(x)

 

 

 

 

 

Будем считать, что p(x) C1(Ω), p(x) ≥ p0 > 0, q(x) C(Ω), f(x) L2(Ω), σ(x) L2(∂Ω). В случае первой краевой задачи считаем, что ϕ(x) является следом на ∂Ω некоторой функции из H1(Ω).

Определение 11.1. Функция u(x) H1(Ω) называется обобщенным

решением первой краевой задачи (11.1),(11.2), если ее след на ∂Ω равен

ϕ(x) и она удовлетворяет при всех v H1 (Ω) интегральному тождеству

Z

(p(x)(r u, r v) + quv) dx = Z

fv dx.

Ω

Ω

 

62

Определение 11.2. Функция u(x) H1(Ω) называется обобщенным решением третьей краевой задачи (11.1),(11.3), если при всех v H1(Ω)

она удовлетворяет интегральному тождеству

 

 

Z (p(x)(r u, r v) + quv) dx + Z

p(x)σuv dS =

 

Ω

∂Ω

 

 

Z

Z

 

 

fv dx +

p(x)ϕv dS.

(11.4)

Ω∂Ω

Замечание. При σ ≡ 0 задача (11.1),(11.3) является второй краевой

задачей и формула (11.4) определяет ее обобщенное решение.

 

Пример 11.1. Пусть в области Ω задано уравнение

 

−div(p(x)ru(x)) + q(x)u = f(x),

(11.5)

а на границе области условие

 

u(x)|∂Ω = ϕ(x).

(11.6)

Доказать, что если p(x) C1(Ω), p(x) ≥ p0 > 0, q(x) C(Ω), f(x) L2(Ω), ϕ(x) C(∂Ω), то классическое решение u(x) C2(Ω) C1(Ω) задачи (11.5), (11.6) является е¸ обобщенным решением.

Решение. Так как классическое решение u(x) уравнения (11.5) – функция непрерывная в Ω, е¸ след на границе области Ω совпадает с е¸ значением. Поэтому краевое условие (11.6) будет верно и в смысле следа.

Умножим

уравнение

(11.5)

на

 

произвольную

функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v H1 (Ω) и проинтегрируем по области Ω

 

 

 

 

 

 

 

 

 

Z

div(pru)v dx + Z

quv dx = Z

fv dx.

(11.7)

 

Ω

 

 

 

 

Ω

 

 

 

Ω

 

 

 

 

Преобразуем первое слагаемое левой части тождества

 

Z

m

 

∂u

m

Z

 

 

∂u

 

 

 

 

(p

 

)v dx = − i=1

 

 

(p

 

)v dx.

 

i=1

∂xi

∂xi

∂xi

∂xi

 

Ω

X

 

 

X

Ω

 

 

 

 

 

 

 

В последнем выражении применим к каждому слагаемому формулу интегрирования по частям (см. задачу 10.10). Получим, что

m

Z

∂xi

(p∂xi )v dx =

i=1

X

 

 

∂u

Ω

 

 

 

 

63

Z

pv i=1

∂xi

cos(n, xi) ds + Z

p i=1

∂xi ∂xi

dx =

 

m

∂u

 

m

∂u ∂v

 

∂Ω

X

 

Ω

X

 

 

ZZ

pv∂n∂u ds + p(ru, rv) dx,

∂Ω Ω

1

так как v H (Ω), то след v|∂Ω = 0 и первый интеграл в этом выражении равен 0.

Здесь cos(n, xi), i = 1, 2, . . . , – направляющие косинусы единичной внеш-

ней нормали n к границе области Ω.

 

 

 

 

 

Согласно

последнему

равенству

первое слагаемое в ле-

вой части

тождества

(11.7) можно

заменить

на

выражение

R

 

 

 

 

 

 

 

Ω p(ru, rv) dx. После замены получим,

что для v

H1

(Ω) имеет

место интегральное тождество

quv dx = Z

 

 

 

 

Z

p(ru, rv) dx + Z

fv dx.

 

 

 

Ω

 

Ω

 

Ω

 

 

 

Пример 11.2. Рассмотрим в области Ω задачу (11.5), (11.6). Доказать, что если p(x) C1(Ω), p(x) ≥ p0 > 0, q(x) C(Ω), q(x) > 0,

f(x) L2(Ω), то задача (11.5), (11.6) имеет единственное обобщенное

решение u H1 (Ω).

Решение. Рассмотрим гильбертовое пространство H1 (Ω) со скалярным произведением

(u, v)I = Z

p(ru, rv) dx + Z

quv dx,

Ω

Ω

 

 

 

которое эквивалентно исходному скалярному произведению в H1 (Ω)

(u, v) = Z (ru, rv) dx + Z

uv dx.

Ω

Ω

 

Обозначим за k · kI норму, порожденную скалярным произведением (·, ·)I , за k · k норму, порожденную скалярным произведением (·, ·). Данные нормы эквивалентны (см. упражнение 10.14).

64

Для функции f(x) определим функционал

F (v) = Z

fv dx v H1 (Ω).

Ω

 

Заданный функционал линеен.

Покажем, что данный функционал ограничен.

R

|F (v)| ≤ |fg| dx ≤ (неравенство Коши – Буняковского)

 

 

 

 

 

Ω

 

 

 

v2 dx 1/2

 

 

 

 

 

 

Ω f2 dx 1/2 Ω

= kfkL2(Ω)kvkL2(Ω)

 

 

k

 

k

 

R

(Ω)

R

 

 

 

 

k

kH1

в силу эквивалентности норм в H1

 

 

f

 

L2(Ω) v

 

 

(Ω)

 

≤ CkfkL2(Ω)kvkI = CekvkI .

Из последнего неравенства следует, что F ограничен. Тогда согласно тео-

реме Рисса [3] о представлении линейного ограниченного функционала в

гильбертовом пространстве существует единственный элемент u H1 (Ω) такой, что линейный ограниченный функционал F (v) может быть представлен в виде F (v) = (u, v)I . Последнее означает, что функция u(x) удовлетворяет интегральному тождеству

Z

p(ru, rv) dx + Z

quv dx = Z

fv dx

Ω

Ω

Ω

 

для любой функции v H1 (Ω), т.е. функция u является обобщенным решением задачи (11.5), (11.6).

Пример 11.3. Доказать, что для всех функций v C2[0, 1] верно

 

 

 

 

1

2

+ 2xv) dx + v2(0) + v2(1)

 

 

 

41

 

 

 

 

 

 

 

(v0

 

 

 

 

 

 

неравенство

 

 

≥ −270 . Имеет ли место знак

0

 

 

 

 

 

 

 

 

равенства

для какой-либо функции?

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выражение

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(v02 + 2xv) dx + v2(0) + v2(1) =

1

Решение. Рассмотрим

 

1

R0

 

 

 

 

 

 

 

 

 

R0

v02 dx + v2(0) + v2(1) + 2

xv dx = ( v

k

 

1

 

)2 + 2(x, v)

 

= Φ(v).

 

 

 

 

 

 

1

 

R0

k

 

 

1/2

 

 

 

L2(0,1)

 

 

 

 

 

 

 

 

R

 

 

 

 

 

H (0,1)

 

 

 

 

 

 

v

 

C

,

 

 

 

 

 

 

 

 

 

 

 

v

Здесь kvkH1(0,1) =

 

1 0

v02 dx + v2(0) + v2(1)

 

. Известно, что существу-

ет функция

 

0

 

 

(0 1), реализующая минимум функционала Φ( ) [20].

65

Рассмотрим функцию F (t) = Φ(v0 + tw), где v0 - функция, на которой функционал Φ(v) принимает минимальное значение, w(x) - произвольная функция из H1(0, 1), t - произвольное вещественное число, тогда очевидно, что F (0) = Φ(v0), а так как Φ(v) ≥ Φ(v0), то F 0(0) = 0. То есть

F 0(t) =

 

1

(v00 + tw0)2 dx + (v(0) + tw(0))2 + (v(1) + tw(1))2 + 2 Z0

1

Z0

x(v + tw) dx =

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

Z0

v002 dx + v2(0) + v2(1) + 2 Z0

xv dx+

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

2t(Z0

v00 w0 dx + v0(0)w(0) + v0(1)w(1) + Z0

xw dx)+

 

 

 

 

 

t2(Z0

1

 

 

 

 

 

 

 

 

 

 

w02 dx + w2(0) + w2(1))

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

и1

F 0(0) = 2(R0

v00 w0 dx + v0(0)w(0) +1v0(1)w(1) + R0

xw dx). Отсюда имеем

R0

v00 w0 dx + v0(0)w(0) + v0(1)w(1) + R0

xw dx = 0 для любой функции w

H1(0, 1).

Считая, что v0 C2([0, 1]), проинтегрируем по частям первый интеграл в этой сумме, получим

 

1

1

 

v00 w|01 Z0

v000w dx + v0(0)w(0) + v0(1)w(1) + Z0

xw dx = 0.

После очевидных преобразований имеем

 

 

 

1

 

 

(v00 (1) + v0(1))w(1) + (−v00 (0) + v0(0))w(0) + Z0

(−v000 + x)w dx = 0.

Рассмотрим w H1 (0, 1). След функции w в граничных точках интервала (0, 1) равен 0 и

66

1

Z

(−v000 + x)w dx = 0,

0

то есть

(−v000 + x, w)L2(0,1) = 0.

Так как H1 (0, 1) всюду плотно в L2(0, 1), то по теореме об ортогональности −v000 + x = 0, и мы имеем следующее равенство:

(v00 (1) + v0(1))w(1) + (−v00 (0) + v0(0))w(0) = 0,

где w(0) и w(1) могут принимать произвольные значения, но тогда v00 (1) + v0(1) = 0, −v00 (0) + v0(0) = 0.

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

v000 − x = 0,

v00 (1) + v0(1) = 0, −v00 (0) + v0(0) = 0.

Решая ее, получим, что v(x) = 16 x3 + c1x + c2. В силу краевых условий

найдем

v(x) = 16x3 29(x + 1).

1

Вычислим значение функционала на функции v(x): R [(12 x2 29 )2 + 13 x4

0

49 (x2 + x)] dx + (−29 )2 + (16 49 )2 = −27041 . Таким образом, мы доказали требуемое неравенство и нашли функцию, для которой верно равенство.

Задачи и упражнения

11.1.Пусть u(x) классическое решение задачи

u = f,

u|∂Ω = ϕ.

Показать, что если u C1(Ω) и f L2(Ω), то u(x) является обобщенным решением этой задачи.

67

11.2.Пусть u(x) классическое решение задачи

u = f,

 

∂u(x)

+ σ(x)u ∂Ω

= ϕ(x).

∂n

 

 

 

 

 

 

 

 

 

 

 

 

Показать, что u(x) является обобщенным решением этой задачи.

11.3.Если u(x) обобщенное решение задачи

u = f,

u|∂Ω = ϕ

и u(x) C2(Ω) ∩ C(Ω), то u(x) является классическим решением этой задачи.

11.4.Если u(x) обобщенное решение задачи

u = f,

 

∂u(x)

+ σ(x)u ∂Ω

= ϕ(x)

∂n

 

 

 

 

 

 

 

 

 

 

 

 

и u(x) C2(Ω) ∩ C1(Ω), то u(x) является классическим решением этой задачи.

11.5. Пусть u(x) классическое решение задачи

−div(p(x)ru(x)) + q(x)u = f(x),

u(x)|∂Ω = ϕ(x),

принадлежащее H1(Ω). Показать, что u(x) является обобщенным решением этой задачи.

11.6. Пусть задана функция f(x) L2(Ω). Доказать, что функционал

Z

F (v) = f(x)v(x) dx

Ω

является линейным непрерывным функционалом над H1(Ω) (∂Ω C1). 11.7. Пусть задана функция u(x) H1(Ω). Доказать, что функционал

Z

F (v) = u(s)v(s) ds

∂Ω

является линейным непрерывным функционалом над H1(Ω) (∂Ω C1).

68

11.8. Пусть задана функция u(x) H1(Ω). Доказать, что функционал

F (v) = Z u(s)v(s) ds + Z

v(x) dx

∂Ω

Ω

 

является линейным непрерывным функционалом над H1(Ω) (∂Ω C1). 11.9. Доказать существование и единственность обобщенного решения задачи

− u + u = f,

∂u

∂n|∂Ω = 0

в классе функций H1(Ω).

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

− u = f,

u|∂Ω = 0.

11.11. Рассмотрим при f L2(Ω) функционал

ZZ

E(v) = |rv|2 dx − 2 fv dx

Ω Ω

на множестве функций v H1(Ω), для которых v|∂Ω = ϕ(x), где функция ϕ(x) является следом на ∂Ω некоторой функции из H1(Ω). Показать, что функция u(x), на которой функционал E(v) достигает минимального значения, есть обобщенное решение задачи

− u = f,

u|∂Ω = ϕ.

11.12. Рассмотрим при f L2(Ω), p C(Ω), q C(Ω), p(x) ≥ p0 > 0, q(x) ≥ 0 функционал

Z Z Z

E1(v) = p|rv|2 dx + qv2 dx − 2 fv dx

Ω Ω Ω

на множестве функций v H1(Ω), для которых v|∂Ω = ϕ(x), где функция ϕ(x) является следом на ∂Ω некоторой функции из H1(Ω). Показать,

69

что функция u(x), на которой функционал E1(v) достигает минимального значения, есть обобщенное решение задачи

−div(p(x)ru(x)) + q(x)u = f(x),

u|∂Ω = ϕ.

11.13. Рассмотрим при f L2(Ω), ϕ L2(∂Ω), σ C(∂Ω), σ(x) ≥ 0 на ∂Ω, σ(x) 6≡0, функционал

Z Z Z Z

˜ |r |2 2 − −

E1(v) = v dx + σv dS 2 fv dx 2 ϕv dS

Ω ∂Ω Ω ∂Ω

на множестве функций v H1(Ω). Показать, что функция u(x), на которой

˜

функционал E1(v) достигает минимума, есть обобщенное решение задачи

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− u = f,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂n

+ σ(x)u ∂Ω = ϕ(x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂u(x)

 

 

 

 

 

 

 

 

 

 

 

11.14.

1

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

минимум функционала

R0

(v02 +

 

Найти функцию v , реализующую

 

v2) dx + 2 R0

v dx в классе H1 (0, 1).

 

v

 

C2

,

v

(1) = 0 имеет место

11.15. Доказать, что для всех функций

 

 

[0 1],

 

 

 

 

1

 

 

 

 

5

 

 

 

v

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

неравенство

 

v dx ≤

+

(0)

+ 41

 

v02 dx. Найти функцию из этого класса,

0

24

 

4

0

для которой Rдостигается

 

равенство.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

H1(Ω)

1

 

|r |

 

 

 

 

 

 

 

{

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

11.16. Найти

 

 

 

R

(

 

 

v 2 + 2 sin x1 sin x2v) dx , где Ω = 0

 

x1

 

 

inf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π, 0 ≤ x1 ≤ π}.

70

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