u_practice
.pdfZ
(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