CC BY
45
____________УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА
Том 157, кн. 2 Физико-математические науки
2015
УДК 519.63
АППРОКСИМАЦИЯ НАИМЕНЬШЕГО СОБСТВЕННОГО ЗНАЧЕНИЯ НЕЛИНЕЙНОЙ ЗАДАЧИ ШТУРМА-ЛИУВИЛЛЯ
В.С. Желтухин, С.И. Соловьёв, П.С. Соловьёв
Аннотация
Изучены свойства минимального собственного значения, отвечающего положительной собственной функции, нелинейной задачи на собственные значения для обыкновенного дифференциального уравнения. Эта задача аппроксимирована сеточной схемой метода конечных элементов на равномерной сетке с линейными конечными элементами. Исследована погрешность приближенных решений. Теоретические результаты подтверждены численными экспериментами для модельной задачи на собственные значения.
Ключевые слова: собственное значение, положительная собственная функция, нелинейная задача на собственные значения, обыкновенное дифференциальное уравнение, задача Штурма-Лиувилля, метод конечных элементов.
Введение
В статье исследуется задача определения наименьшего собственного значения А € Л, Л = [0, ж), и соответствующей положительной собственной функции u(x), x € О, О = (0, п), О = [0, п], удовлетворяющих обыкновенному дифференциальному уравнению и однородным граничным условиям Дирихле
— (p(Xs(x))u')' = r(As(x))u, x € О,
u(0) = u(n) = 0.
Коэффициенты дифференциального уравнения зависят от спектрального параметра А € Л. Пусть функции р(р), г(р), р € Л, s(x), x € О являются непрерывными положительными, функция р(р), р € Л - неубывающая ограниченная, а функция г(р), р € Л - неубывающая неограниченная. Предположим дополнительно, что производные р'(р), г'(р), р € Л, существуют и непрерывны.
Задачу (1) сформулируем в эквивалентном виде как задачу о нахождении корня линейного уравнения. Пусть р(р) для фиксированного р € Л обозначает наименьшее простое собственное значение, отвечающее положительной собственной функции u(x) = u^(x), x € О, параметрической линейной задачи на собственные значения для обыкновенного дифференциального уравнения
— ^р^у-и'У = y(р)r(рs(x))u, x € О, u(0) = u(n) = 0.
Тогда наименьшее собственное значение А € Л задачи (1) определяется как наименьший корень характеристического уравнения
Y(р) = l, р € Л.
(3)
40
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
41
В работе [1] предложено простое достаточное условие существования собственного значения Л € Л как корня характеристического уравнения (3) в терминах коэффициентов уравнения задачи (1). Это условие состоит в выполнении неравенства р(0) > г(0). В разд. 1 настоящей статьи доказывается непрерывная дифференцируемость функции y(p) , р € Л, и выводится формула для вычисления производной Y1 (р), р € Л. Доказательство основано на использовании вариационных формулировок задач (1) и (2) в гильбертовом пространстве. Установлено также необходимое и достаточное условие для существования собственного значения Л € Л задачи (1). В разд. 2 определяется сеточная схема метода конечных элементов для решения задачи (1) как конечномерная аппроксимация вариационной задачи в гильбертовом пространстве. Доказаны оценки погрешности аппроксимации собственного значения Л € Л задачи (1) и соответствующей собственной функции и(х), x € Q, при выполнении условия y '(Л) = 0 .В разд. 3 приведены результаты численных экспериментов для модельной задачи.
Задача вида (1) возникает при моделировании плазмы высокочастотного разряда пониженного давления [2-5]. Минимальное собственное значение, отвечающее положительной собственной функции, определяет условие, необходимое для поддержания стационарного высокочастотного индукционного разряда пониженного давления.
Приближенные методы решения задачи на собственные значения для обыкновенного дифференциального уравнения
— (р(Л,х)и')' + д(Л,х)и = Лг(Л,х)и, x € Q, и(0) = u(l) = 0,
с коэффициентами, зависящими от спектрального параметра, изучались в работах [6-8]. В этих работах предполагалось, что функции р(р,х) > 0, q(p,x) ^ 0, г(р,х) > 0, р € Л, x € Q, Q = (0,l), Q = [0, l], являются достаточно гладкими. Функции р(р,х), q(p,x), р € Л, х € Q, предполагались невозрастающими по первому аргументу р € Л, а функция г(р,х), р € Л, х € Q, - неубывающей по первому аргументу р € Л при фиксированном х € Q.
Чтобы сформулировать задачу (4) в эквивалентном виде обозначим через y(р) для фиксированного р € Л минимальное собственное значение, отвечающее положительной собственной функции и(х) = и^(х), х € Q, параметрической линейной задачи на собственные значения для обыкновенного дифференциального уравнения
— (р(р,х)и)' + q^x^ = д (р)г(р,х)и, х € Q, и(0) = и(1) = 0.
Тогда минимальное собственное значение Л задачи (4) определяется как наименьший корень характеристического уравнения
1(р) = P, р € Л.
(6)
В работах [6, 7] исследована погрешность метода конечных разностей для задачи (4) при г(р,х) = r(x), р € Л, х € Q. В [8] получены оценки погрешности метода конечных элементов для задачи (4). В этих случаях функция y(р), р € Л, невозрастающая и существует единственный корень уравнения (6). В отличие от работ [6-8], в настоящей статье функция y(р), р € Л, из характеристического уравнения (3) может быть немонотонной.
Задачи на собственные значения с нелинейной зависимостью от спектрального параметра имеют многочисленные приложения в науке и технике [9-18]. Вычислительные алгоритмы решения нелинейных спектральных задач исследовались
42
В.С. ЖЕЛТУХИН И ДР.
в работах [9, 19-22]. Погрешность приближенных методов решения задач на собственные значения с монотонным и немонотонным вхождением спектрального параметра изучалась в работах [23-32]. Исследование погрешности приближенных методов решения нелинейных спектральных задач опирается на применение результатов для линейных задач на собственные значения [33-43].
1. Вариационная постановка задачи
Пусть Л = [0, ж), 0 = (0, п), H = L2(0) - вещественное гильбертово пространство Лебега с нормой | • |о и со скалярным произведением (■, -)q при
т \ V2
2
М0=[ftvtx))2<fc) • м=/и(х)ф)d Vих"GH
Обозначим через V = {v : v,v' G H,v(0) = v(n) = 0} вещественное гильбертово пространство Соболева с нормой | ■ |i и со скалярным произведением (■, -)i при
т \ I/2
(v'(x))2
Wl = (/(v'(x))2 dx) • "vv)1 = /Axy>x) dx Vu,v G V
\o / 0
При фиксированном д G Л определим билинейные формы
= j pMWv' dx Ь{д, u.v) = j
для и, v G V и отношение Рэлея
R(^,v)
а(д, v, v) b(p,v,v)
Vv G V \ {0}.
Положим K = {v : v G V, v(x) > 0, x G 0}.
Дифференциальную задачу (1) представим в виде вариационной нелинейной задачи на собственные значения: найти минимальное число Л G Л и функцию и G K, Ь(Л, и, и) = 1, такие, что
а(Л, и, v) = Ь(Л, u,v) V v G V.
(7)
При фиксированном д G Л введем вариационную линейную параметрическую задачу на собственные значения: найти число д(д) и функцию и = им G K, Ь(д,и,и) = 1, такие, что
а(д,ир) = y(д)Ь(д,ир) Vv G V. (8)
Для собственного значения задачи (8) справедлив вариационный принцип минимума
Y (Д)= m\nn Жд^).
vev\{о}
Введем вспомогательную линейную задачу на собственные значения: найти минимальное число к и функцию и G K, (и,и)о = 1, такие, что
(u,v)i = k(u,v)q Vv G V.
(9)
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
43
Решения этой задачи определяются по формулам
к =1, u(x) = у 2 sin x, x € П,
и выполняется принцип минимума
к = min S(v)
«еП{0}
с отношением Рэлея
S(v) = j^Vi V v € V \{0}.
(v,v)0
Обозначим через Yi(p), ui(p) = ui(p,x), х € П, p € Л, i = 1, 2,..., собственные значения и собственные функции, удовлетворяющие уравнению (8) и такие, что
Yip < 72(м) < • • • < Yi(p) < ••• , a(p,ui(p),uj (p)) = Yi(p)sij, b(p,ui(p),uj (p)) = = hj, i,j = 1,2,..., ui(p) = uM € K, Yi(p) = Y(p), Yi(p) ^ ж при i ^ ж; функции ui(p), p € Л, i = 1, 2,..., образуют полную систему в пространстве V. Положим
s0 = min s(x), s1 = max s(x), p1 = sup p(p). xEQ xEQ mE Л
Для фиксированного v € Л обозначим Д = [0, v],
ai = p(0), a2 = max p(psi), аз = max p'(psi)si, meA meA
f3i = r(0), (32 = max r(psi), (33 = max r'(psi)si.
meA meA
Лемма 1. Для p,p € Д справедливы неравенства
aijvl2 ^ a(p,v,v) ^ a2|v|2 Vv € V,
AM2 < b(p,v,v) < e2|v|2 Vv € H,
la(p,v,v) - a(n,v,v) < a3p - рЦиЦ Vv € V,
\b(p,v,v) - b(p,v,v) < вз^ - П Hi Vv € H.
Доказательство. Используя соотношения ai ^ p(ps(x)) ^ a2, \p(ps(x)) - p(ps(x))| ^ а3Д - p^ p, p € Д, x € П, для любого v € V получаем
П
ai|v|i ^ a(p,v,v) = jp(ps(x))(v')2 dx ^ a2 v2,
0
Hp^^')
П
a(p,v,v)| <
|p(ps(x)) - p(ps(x))Kv')2 dx < аз|р - p| |v|21.
0
Аналогично устанавливаются свойства билинейной формы b(p, •, •).
□
Через а будем обозначать различные положительные постоянные, не зависящие от р,р € Д. Положим ||v||a(M) = a(p,v,v), ||v||2(m) = b(p,v,v) для v € V, p € Л.
44
В.С. ЖЕЛТУХИН И ДР.
Лемма 2. Если л,П € А и величина \л — п\ достаточно мала, то существует постоянная а, для которой справедливы оценки
\и(Ш) — ъШ p а \л — п\, * = 1, 2.
Доказательство. Утверждение вытекает из [8] с учетом леммы 1. □
Лемма 3. Если л,П € А и величина \л — п\ достаточно мала, то существует постоянная а , для которой выполняется оценка
\им P а \л — ц\.
Доказательство. Зафиксируем л € А. Представим элемент y = и^ в виде разложения
y = Pi(n)ui(n) + w1(n),
и\{р)
где
wi(n) = ^2Pi(n)ui(n), Pi (n) = Ъ(п,У,Щ(п)), * = 1 2,---, П € А■
i=2
Согласно лемме 2 выберем n € А так, чтобы 72(п) — Yi(^) > 0. Обозначим
S(n,v)
- sup
\{0}
\a(n,v,w) - ji(^)b(n,v,w)\
У wIL(n)
Докажем оценку
ll“’i(n)ll»«> p Sin-Vl
Нетрудно проверить, что
a(n,y,wi(n)) = a(n,wi(n),wi(n)),
b(n,y,wi(n)) = b(n,wi(n),wi (П)), a(n,Wi(n),Wi(n)) > Y2(n)b(n,wi(n),wi(n)). Поэтому приходим к соотношениям
a(n,y,wi(n)) — 7i(V)b(n, y,wi(n)) =
= a(n,wi(n),wi(n)) — Yi(^)b(n,wi(n),wi(n)) >
Y2(n) — И(л)
12(П)
a(n,wi(n),wi(n)),
которые приводят к требуемой оценке. Установленная оценка дает
Y2(n)
lly — Pi(n)ui(n)Ia(V) = Iwi(n)Ia(V) P ^(n) — Y1Ш) 5(n^ y) P а \л — n\■ Кроме того, имеем
Pi(V) = lei(n)ui(n)lb(n) P lui(^)yb(n) + Ь-Рл) — Pi(n)ui(n)yb(n) P
P IK^IIb^) + \ IK^Ikn) — lui(^)lb(^)\ + llui(Ш) — Pi(n)ui(n)yb(n) P
P 1+ а \л — n\,
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
45
Pi(v) = \\fji{n)ui{n)\Wri) > \ui(M)\b(n) - \ui(M) - fii(n)uiШ\ь(п) >
> \ui(M)\b(p) - | \ui(M)\b(n) - \ui(M)\b(M)| - \\ui(M) - ei(v)ui(v)\\b(ri) >
> 1 - a \m - n\,
поскольку
1 llui (M) \ ь(п) — llui(M)llb(p)
\b(n, щ(м), ui(p)) - Ь(м, ui(^), ui(p))\
IIui(m) llb(n) + llui (m) II b(p)
< a \m - n\.
Следовательно,
\1 - < a\m - n\-
Таким образом, собирая предыдущие оценки, заключаем
V^i \up - un\i ^ \ui(M) - ui(n)lla(n) <
< \\ui(M) - ei(n)ui(n)\\a(n) + \\ui(n) - ei(n)ui(n)\\a(n) =
= \\ui(M) - ei(n)ui(n)\\a(n) + VхY(n)\1 - Pi(n)\ a \m - n\-
Лемма доказана.
Лемма 4. Для м G Л имеет место равенство
l' (M)
d
R(M,v)
□
при v = up .
Доказательство. Для m, n G Л, m = n, справедливы соотношения
Y(M) - Y(n) =
M - n
= R(m, up) - R(n,uv) = R(m, u^) - R(n, up) + R(n, up) - R(n,uv) m - n m - n m - n
Докажем, что
G(n) = R(n,up) - R(n,un) . 0 M - n
при n ^ M. Предположим, что существует последовательность n' ^ M, для которой \G(n)\ > а при n = n' ^ M. Обозначим e = m - n, w = (up - un)/(m - n). Так как согласно лемме 3 выполняется неравенство \w\i ^ а при n = n' ^ M, то существует подпоследовательность n'' ^ M такая, что w ^ W в V при n = n'' ^ M, где W G V, символ ^ обозначает слабую сходимость в гильбертовом пространстве V. Тогда получим
R(n, up) - R(n, uv) = R(n, up) - R(n, up + ew) =
= a(n, up, up) a(n, up + ew,up + ew) = A(n)
b(n,up,up) b(n, up + ew, up + ew) B(n)
где
A(n) = 2eAi(n) + e2A2(n),
Ai(n) = a(n,up,up) b(n,up,w) - b(n,up,up) a(n,up,w),
46
В.С. ЖЕЛТУХИН И ДР.
A-2(n) = a(n,u^,u^) b(n,w,w) — b(n,Up,Up) a(n,w,w), B(n) = b(n, Up, Up) b(n, Up + ew, Up + ew).
Учитывая, что
Ai(n) - 0, 1A2 (n)\ <
получаем
A(n) 0, B(n) — h
при n = n" — И ■ В результате приходим к противоречию, так как
G(n)
2eAi(n) + £2A2(n) eB(n)
0
при n = n" — И■
Таким образом, заключаем
7Чр) = lim yW - Y(n) = lim R(l,‘’ UP) - R(n’UP) = Arm
n^p И — n n^p и — n Аи
при v = Up■ Лемма доказана.
Обозначим
П П
a{,„„,„) = ]P'(M*M*)UV V{P„U„v) = j
□
для и € Л, U, v € V■
Лемма 5. Функция y'(и) , И € Л является непрерывной и вычисляется по формуле
!'(И) = а'(И,v,v) — Y(И)ь'(И,v,v)
при и € Л, v = Up .
Доказательство. Для v = Up имеем
l' (И)
(a(V,v,v)\' =
\b(p,v,v))
а'(и, v, v) Ь(и, v, v) — а(и, v, v) b'(p, v, v)
ь2(иаа)
= а'(и, v,v) — 7(и)ь'(И, v> v)-
Отсюда с помощью соотношений
l(n) - 1(и), а (n, Un, Un) - а (и, up, Up), b Un, Un) - b (и, up, Up),
при n — И выводим, что y'(n) — l'(и) при n — И ■ Лемма доказана. □
Теорема 1. Для существования наименьшего простого собственного значения задачи (7), отвечающего положительной собственной функции, необходимо и достаточно, чтобы y(£) > 1 для некоторого £ € Л.
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
47
Доказательство. Согласно лемме 2 функция y(p) , p € Л, является непрерывной. Используя вариационные принципы минимума для наименьших собственных значений задач (8) и (9), получим
p(ps(x))(v')2 dx
y(p) = min R(p, v) = min
0
<
Pi
Pi
vEV\{0}
vEV\{0}
r(ps(x))v2 dx
"(ps o)
"(pso)
при p ^ ж, поскольку r(p) ^ ж при p ^ ж. Поэтому для существования наименьшего корня Л € Л уравнения (3) необходимо и достаточно, чтобы y(£) > 1 для некоторого £ € Л. Корень Л € Л определяет минимальное собственное значение задачи (7). Собственное значение Л € Л является простым и отвечает положительной собственной функции, поскольку y (p) - простое собственное значение параметрической задачи (8) при p = Л, соответствующее положительной собственной функции. Теорема доказана. □
2. Сеточная схема метода конечных элементов
Зададим разбиение отрезка [0, п] равноотстоящими точками Xi = ih, i = = 0,1,...,п, на элементы ei = (xi-i,xi), i = 1,2,...,n, h = п/n. Пусть Vh есть подпространство пространства V, состоящее из непрерывных функций vh, линейных на каждом элементе ei = (xi-i, Xj), i = 1, 2,... ,n. Размерность подпространства Vh есть Nh = n — 1. Обозначим Kh = {vh : vh € Vh, vh(x) > 0, x € 0}.
Метод конечных элементов бесконечномерной вариационной задаче (7) ставит в соответствие конечномерную задачу: найти минимальное число Л1'1 € Л и функ-
цию uh € Kh, Ь(Л1 ,uh ,uh) = 1, такие, что
(Л1, uh, vh) = Ь(Лh, uh, vh) У vh € Vh.
(10)
При фиксированном p € Л введем конечномерную линейную параметрическую задачу на собственные значения: найти минимальное число Yh(p) и функцию
= uh^ € Kh, b(p, uh, uh) = 1, такие, что
a(p, uh, vh) = Yh(p)b(p, uh, vh) У vh € Vh.
(11)
Наименьшее собственное значение задачи (11) удовлетворяет принципу минимума
Yh(p) = min R(p,vh).
vh evfc\{0}
Введем вспомогательную линейную задачу на собственные значения: найти минимальное число Kh и функцию uh € Kh, (uh,uh)o = 1, такие, что
(uh,vh)i = Kh(uh, vh)o У Vh € Vh. Для этой задачи также выполняется принцип минимума
(12)
Kh = min S(vh).
vhEVh\{0}
Известно, что 1 ^ Kh ^ 1 при h ^ 0.
Минимальное собственное значение Лh задачи (10) является наименьшим корнем характеристического уравнения
Yh(p) = 1, p € Л.
П
0
П
a
h
u
(13)
48
В.С. ЖЕЛТУХИН И ДР.
Лемма 6. Если л, п € А и величина \л — п\ достаточно мала, то существует постоянная а, для которой справедливы оценки
hi (л) — Yi (П)\ < а \л — nl i = 1, 2.
Лемма 7. Если л, П € А и величина \л — п\ достаточно мала, то существует постоянная а , для которой выполняется оценка
u— uhh\i < а \л — п\.
Лемма 8. Для л € Л имеет место равенство
Ьк(л)У
d
dл
R(p,v)
при v = u,y .
Лемма 9.
формуле
Функция (yh(p))', Л € Л является непрерывной и вычисляется по (Yh(л))' = а'(л, vh, vh) — Yh(v) Ъ'Д, vh, vh)
при p € Л, vh = .
Леммы 6-9 доказываются аналогично леммам 2-5 соответственно.
Теорема 2. Для существования наименьшего простого собственного значения задачи (10), отвечающего положительной собственной функции, необходимо и достаточно, чтобы yh(£) > 1 для некоторого £ € Л.
Доказательство. По лемме 6 функция Yh(h), л € Л, является непрерывной. Применяя вариационные принципы минимума для наименьших собственных значений задач (11) и (12), получаем
П
/ p(ys(x))((vh)')2 dx
1h(R)
min R(u,vh)
'еуд{0}
min -
^еУД{0}
/ r(ys(x))(vh)2 dx
<
Pi
r(yso)
h
0
V
0
при л ^ ж, поскольку г(л) ^ ж при л ^ ж. Поэтому для существования наименьшего корня Xh € Л уравнения (13) необходимо и достаточно, чтобы yh(£) > 1 для некоторого £ € Л. Корень Xh € Л определяет минимальное собственное значение задачи (10). Собственное значение Xh € Л является простым и отвечает положительной собственной функции, поскольку Yh(h) - простое собственное значение параметрической задачи (11) при л = Xh, соответствующее положительной собственной функции. Теорема доказана. □
Через c будем обозначать различные положительные постоянные, не зависящие от h. Для фиксированного л € Л введем оператор Ph(h) '■ V^Vh по правилу a(p,u — Ph(p)u,vh) = 0 для любого vh € Vh, где u € V, \u — Ph(y)u\i ^ ch. Положим
Ph = Ph(X). _
Обозначим через Yi^) и uh(y) = uh(y,x), x € £, л € Л, i = 1,2,...,Nh, собственные значения и собственные функции, удовлетворяющие уравнению (11) и такие что YihiH) < Y2M < ••• < Yn h Ы, (л)) = Yp^--
Jij ?
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
49
b(^,Uh(^),ullj(^)) = Sij , i,j = 1, 2,...,Nh, u’l(^) = u* e Kh, Y1 (ц) = 7*(ц); функции uh(^), ц e Л, i = 1, 2,..., Nh, образуют полную систему в пространстве Vh . Для ц e Л и достаточно малых h выполняются оценки 0 ^ 7*(ц) — 7г(ц) ^ ch2, i = 1, 2, U — uM|i < ch.
Теорема 3. Пусть 7'(X) = 0. Тогда для достаточно малых h выполняется оценка
0 < Xh — X < ch2.
Доказательство. Имеем, что (7 h(ц))' ^ 7 '(ц) при h ^ 0 .Действительно,
(Yh{^))' = а'(Ц, vh, vh) — Yh(Ц) Ь'(Ц, vh, vh) ^ а'(ц, v, v) — 7(ц) Ъ'(ц, v, v) = 7'(Ц)
при h ^ 0, ц e Л, vh = u^, v = иц, так как 7h(ц) ^ 7(ц), u^ ^ им в V при
h ^ 0.
Зафиксируем е > 0. Тогда Xh e [X — e,X + е] при достаточно малых h, и существует £h e [X — е,Х + е] такое, что
ci(Xh — X) < —(7h(Zh))'(Xh — X) = 7h(X) — 7h(Xh) = 7h(X) — 7(X) < C2h2
при достаточно малых h, поскольку — (7*(ц))' ^ ci для ц e [X — e,X + е] и достаточно малых h, 7h(Xh) = 7(X) = 1. Теорема доказана. □
Теорема 4. Пусть 7'(X) = 0, u - положительная собственная функция задачи (7), uh - приближение по схеме (10). Тогда для достаточно малых h выполняется оценка
luh — u|i ^ ch.
Доказательство. Положим pi1 = b(Xh, Phu,yh), i = 1,2,...,Nh, где yh = = uh(Xh), i = 1, 2,..., Nh. Поскольку элементы yh, i = 1, 2,..., Nh образуют ортонормированный базис в пространстве Vh, то элемент P*u e Vh можно представить в виде
Phu = ehyh + wh,
где
Nh
'Efty*.
Поскольку 72 (X) — 71 (X) > 0, то 7h(Xh) — 1 = 7h(Xh) — 7 1 (X) =
= (72(X) — 71 (X)) + (7h(Xh) — 7h(X)) + 7(X) — ^(X)) >
> (72(X) — 71 (X)) — ch2 > c
h
1
для достаточно малых h. Обозначим
Ch(u)
sup
гЕУн\{0}
la(Xh, Phu, vh) — Xb(Xh, Phu, vh
lvh 11
V
Тогда Zh(u) ^ ch. Докажем оценку
K|1 < cZh(u).
50
В.С. ЖЕЛТУХИН И ДР.
Нетрудно проверить, что
a(Xh, Phu, wh) = a(Xh, wh, wh),
b(Xh,Phu,w h) = b(Xh,w h,w h), a(Xh,wh,wh) > (Xh)b(Xh, wh,wh).
Поэтому
\wh\Zh(u) > a(Xh, Phu, wh) - b(Xh, Phu, wh) =
a(Xh, wh, wh) - b(Xh, wh, wh) > Yh(Xh) - 1
Y2(Xh)
a(Xh,wh,wh) > c 1 [w^j.
Отсюда получим требуемую оценку.
Теперь, используя доказанные выше оценки, имеем
\Phu - f^yh\i = \wh\i < cZh(u) < ch
при достаточно малых h. Кроме того,
eh = Wehyh\\b(\h) ^ ||uHb(A) + I llul|b(A) - ||uHb(Ah) I + llu - ehyh\\b(\h) ^ 1 + chi
eh = Wehyh\\b(Ah) ^ 11u \ b( A) - 1 ||u||b(A) - 11u \ b(Ah)1 - llu - ehyh\\b(Ah) ^ 1 - ch■,
поскольку
l\ullb(A) - llullb(Ah)1
<
I \u W b( A) llullb(Ah)
I \u W b( A) + ||u||b(Ah)
^ c (Xh - X) ^ ch2,
||u - Pvhyh\\b(Ah) < c^fp2 \u - ehyh\i < c^fp2(\u - Phu\ 1 + \Phu - ehyh\i) < ch.
Следовательно, выводим \1 - f3h\ ^ ch.
В результате заключаем
V&1 \u - u\i < Ilu - yh\\a(Ah) < Ilu
ehy1* II a(Ah) + \\yh - ehyh\\a(Ah) =
= \\u - ehy h\\a(Ah) + \1 - eh\ < ch.
Теорема доказана.
□
3. Численные эксперименты
Для дифференциальной задачи на собственные значения (1) определим коэффициенты по формулам s(x) = а/п - x + 1, x G Q,
P(P)
-p3 + p2 + )x + 2, p G [0,1],
3, p G (1, >x>),
r(P)
1, P G [0,1],
(P - 1)2/2+1, P G (1, to)-
В данном случае функции p(p), r(p), p G Л, являются непрерывно дифференцируемыми, p(0) = 2, r(0) = 1,
Y(0)
p(0)
r(0)
2 > 1.
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
51
Рис. 1. Минимальное собственное значение Л
Табл. 1 Скорость сходимости
n a
10 1.999309
20 1.999832
50 1.999973
100 1.999993
200 1.999998
300 2.000004
На рис. 1 представлен график функции д(д,) , М € Л, задачи (2) и отмечено наименьшее простое собственное значение Л = 1.3248 задачи (1). Этот рисунок служит иллюстрацией к теореме 1 о необходимом и достаточном условии существования собственного значения, отвечающего положительной собственной функции.
В табл. 1 приведены значения порядка сходимости приближенного собственного значения Лн, найденные по формуле
Л&__л^/2
“ = l0g2 ЛЬ/2 - ЛЧ4 '
для h = п/n, n = 10, 20, 50, 100, 200, 300, так, что Лн — Л и cha. Экспериментальные результаты табл. 1 согласуются с теоретическими результатами теоремы 3, а именно а и 2 при n = 10, 20, 50, 100, 200, 300.
Собственное значение Л = 1.3248 дифференциальной задачи получено как предельное значение приближений Лн при h ^ 0. Отметим, что точное собственное значение Л = 1.3248 совпадает с приближенным собственным значением Лн при n = 1000 с приведенной точностью.
Работа выполнена при финансовой поддержке фонда Гумбольдта (Alexander von Humboldt Foundation) и Российского фонда фундаментальных исследований (проекты № 13-01-00908, 14-01-00755, 15-41-02672).
Summary
V.S. Zheltukhin, S.I. Solov’ev, P.S. Solov’ev. Approximation of the Minimal Eigenvalue for a Nonlinear Sturm-Liouville Problem.
Properties of the minimal eigenvalue corresponding to the positive eigenfunction of a nonlinear eigenvalue problem for an ordinary differential equation are studied. This problem
52
В.С. ЖЕЛТУХИН И ДР.
is approximated by a mesh scheme of the finite element method. The error of approximate solutions is investigated. Theoretical results are illustrated by numerical experiments for a model eigenvalue problem.
Keywords: eigenvalue, positive eigenfunction, nonlinear eigenvalue problem, ordinary differential equation, Sturm-Liouville problem, finite element method.
Литература
1. Желтухин В.С., Соловьёв С.И., Соловьёв П.С., Чебакова В.Ю. Вычисление минимального собственного значения нелинейной задачи Штурма-Лиувилля // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2013. - Т. 155, кн. 3. - С. 91-104.
2. Абдуллин И.Ш., Желтухин В.С., Кашапов Н.Ф. Высокочастотная плазменноструйная обработка материалов при пониженных давлениях. Теория и практика применения. - Казань: Изд-во Казан. ун-та, 2000. - 348 с.
3. Желтухин В.С. О разрешимости одной нелинейной спектральной задачи теории высокочастотных разрядов пониженного давления // Изв. вузов. Матем. - 1999. -№ 5. - С. 26-31.
4. Желтухин В.С. Об условиях разрешимости системы краевых задач теории высокочастотной плазмы пониженного давления // Изв. вузов. Матем. - 2005. - № 1. -С. 52-57.
5. Желтухин В.С., Ольков Е.В. О разрешимости одной нелинейной задачи Штурма-Лиувилля // Исслед. по прикл. матем. и информатике. - Казань: Изд-во Казан. ун-та, 2004. - Вып. 25. - С. 59-65.
6. Гулин А.В., Крегжде А.В. Разностные схемы для некоторых нелинейных спектральных задач. Препринт № 153. - М.: ИПМ АН СССР, 1981. - 28 с.
7. Крегжде А.В. О разностных схемах для нелинейной задачи Штурма-Лиувилля // Дифференц. уравнения. - 1981. - Т. 17, № 7. - С. 1280-1284.
8. Соловьёв С.И. Аппроксимация дифференциальных задач на собственные значения с нелинейной зависимостью от параметра // Дифференц. уравнения. - 2014. - Т. 50, № 7. - С. 955-962.
9. Соловьёв С.И. Нелинейные задачи на собственные значения. Приближенные методы. - Saarbriicken: LAP Lambert Acad. Publ., 2011. - 256 с.
10. Goolin A.V., Kartyshov S.V. Numerical study of stability and nonlinear eigenvalue problems // Surv. Math. Ind. - 1993. - V. 3. - P. 29-48.
11. Apel Th., Sandig A.-M., Solov’ev S.I. Computation of 3D vertex singularities for linear elasticity: Error estimates for a finite element method on graded meshes // Math. Model. Numer. Anal. - 2002. - V. 36, No 6. - P. 1043-1070.
12. Lyashko A.D., Solov’yev S.I. Fourier method of solution of FE systems with Hermite elements for Poisson equation // Sov. J. Numer. Anal. Math. Model. - 1991. - V. 6, No 2. - P. 121-129.
13. Solov’ev S.I. A fast direct method of solving Hermitian fourth-order finite-element schemes for the Poisson equation // J. Math. Sci. - 1995. - V. 74, No 6. - P. 1371-1376.
14. Solov’ev S.I. Fast direct methods of solving finite-element grid schemes with bicubic elements for the Poisson equation // J. Math. Sci. - 1994. - V. 71, No 6. - P. 27992804.
15. Жигалко Ю.П., Ляшко А.Д., Соловьёв С.И. Колебания цилиндрической оболочки с присоединёнными жёсткими кольцевыми элементами // Моделирование в механике. - 1988. - Т. 2, № 2. - С. 68-85.
АППРОКСИМАЦИЯ СОБСТВЕННОГО ЗНАЧЕНИЯ
53
16. Жигалко Ю.П., Соловьёв С.И. Собственные колебания балки с гармоническим осциллятором // Изв. вузов. Матем. - 2001. - № 10. - С. 36-38.
17. Карчевский Е.М., Соловьёв С.И. Исследование спектральной задачи для оператора Гельмгольца на плоскости // Дифференц. уравнения. - 2000. - Т. 36, № 4. - С. 563565.
18. Карчевский Е.М., Соловьёв С.И. Существование собственных значений спектральной задачи теории диэлектрических волноводов // Изв. вузов. Матем. - 2003. - № 3. -С. 78-80.
19. Гулин А.В., Крегжде А.В. О применимости метода бисекции к решению нелинейных разностных задач на собственные значения. Препринт № 8. - М.: ИПМ АН СССР, 1982. - 22 с.
20. Гулин А.В., Яковлева С.А. О численном решении одной нелинейной задачи на собственные значения // Вычислительные процессы и системы. - М.: Наука, 1988. -Вып. 6. - С. 90-97.
21. Dautov R.Z., Lyashko A.D., Solov’ev S.I. The bisection method for symmetric eigenvalue problems with a parameter entering nonlinearly // Russ. J. Numer. Anal. Math. Model. -1994. - V. 9, No 5. - P. 417-427.
22. Solov’ev S.I. Preconditioned iterative methods for a class of nonlinear eigenvalue problems // Linear Algebra Appl. - 2006. - V. 415, No 1. - P. 210-229.
23. Вайникко Г.М., Карма О.О. О быстроте сходимости приближённых методов в проблеме собственных значений с нелинейным вхождением параметра // Журн. вычисл. матем. и матем. физики. - 1974. - Т. 14, № 6. - С. 1393-1408.
24. Карма О.О. Асимптотические оценки погрешности приближённых характеристических значений голоморфных фредгольмовых оператор-функций // Журн. вычисл. матем. и матем. физики. - 1971. - Т. 11, № 3. - С. 559-568.
25. Карма О.О. Об аппроксимации оператор-функций и сходимости приближённых собственных значений // Труды ВЦ Тарт. гос. ун-та. - Тарту: ВЦ ТГУ, 1971. - Вып. 24. -С. 3-143.
26. Карма О.О. О сходимости дискретизационных методов отыскания собственных значений интегральных и дифференциальных операторов, голоморфно зависящих от параметра // Труды ВЦ Тарт. гос. ун-та. - Тарту: ВЦ ТГУ, 1971. - Вып. 24. -С. 144-159.
27. Соловьёв С.И. Погрешность метода Бубнова-Галёркина с возмущениями для симметричных спектральных задач с нелинейным вхождением параметра // Журн. вы-числ. матем. и матем. физики. - 1992. - Т. 32, № 5. - С. 675-691.
28. Соловьёв С.И. Аппроксимация симметричных спектральных задач с нелинейным вхождением параметра // Изв. вузов. Матем. - 1993. - № 10. - С. 60-68.
29. Соловьёв С.И. Оценки погрешности метода конечных элементов для симметричных спектральных задач с нелинейным вхождением параметра // Изв. вузов. Матем. -1994. - № 9. - С. 70-77.
30. Соловьёв С.И. Метод конечных элементов для симметричных задач на собственные значения с нелинейным вхождением спектрального параметра // Журн. вычисл. матем. и матем. физики. - Т. 37, № 11. - С. 1311-1318.
31. Даутов Р.З., Ляшко А.Д., Соловьёв С.И. Сходимость метода Бубнова-Галёркина с возмущениями для симметричных спектральных задач с нелинейным вхождением параметра // Дифференц. уравнения. - 1991. - Т. 27, № 7. - С. 1144-1153.
54
В.С. ЖЕЛТУХИН И ДР.
32. Соловьёв С.И. Аппроксимация нелинейных спектральных задач в гильбертовом пространстве // Дифференц. уравнения. - 2015. - Т. 51, № 7. - С. 937-950.
33. Вайникко Г.М. Асимптотические оценки погрешности проекционных методов в проблеме собственных значений // Журн. вычисл. матем. и матем. физики. - 1964. -Т. 4, № 3. - С. 405-425.
34. Вайникко Г.М. Оценки погрешности метода Бубнова-Галёркина в проблеме собственных значений // Журн. вычисл. матем. и матем. физики. - 1965. - Т. 5, № 4. -С. 587-607.
35. Вайникко Г.М. О быстроте сходимости приближённых методов в проблеме собственных значений // Журн. вычисл. матем. и матем. физики. - 1967. - Т. 7, № 5. -С. 977-987.
36. Соловьёв С.И. Метод конечных элементов для несамосопряженных спектральных задач // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2006. - Т. 148, кн. 4. -С. 51-62.
37. Соловьёв С.И. Суперсходимость конечно-элементных аппроксимаций собственных функций // Дифференц. уравнения. - 1994. - Т. 30, № 7. - С. 1230-1238.
38. Соловьёв С.И. Суперсходимость конечно-элементных аппроксимаций собственных подпространств // Дифференц. уравнения. - 2002. - Т. 38, № 5. - С. 710-711.
39. Соловьёв С.И. Аппроксимация вариационных задач на собственные значения // Дифференц. уравнения. - 2010. - Т. 46, № 7. - С. 1022-1032.
40. Соловьёв С.И. Аппроксимация неотрицательно-определенных спектральных задач // Дифференц. уравнения. - 2011. - Т. 47, № 8. - С. 1075-1082.
41. Соловьёв С.И. Аппроксимация знаконеопределенных спектральных задач // Дифференц. уравнения. - 2012. - Т. 48, № 7. - С. 1042-1055.
42. Соловьёв С.И. Аппроксимация дифференциальных задач на собственные значения // Дифференц. уравнения. - 2013. - Т. 49, № 7. - С. 936-944.
43. Solov’ev S.I. Finite element approximation with numerical integration for differential eigenvalue problems // Appl. Numer. Math. - 2015. - V. 93. - P. 206-214.
Поступила в редакцию 06.04.15
Желтухин Виктор Семёнович - доктор физико-математических наук, главный научный сотрудник, Казанский национальный исследовательский технологический университет, г. Казань, Россия.
E-mail: [email protected]
Соловьёв Сергей Иванович - доктор физико-математических наук, доцент кафедры вычислительной математики, Казанский (Приволжский) федеральный университет, г. Казань, Россия.
E-mail: [email protected]
Соловьёв Павел Сергеевич - студент, Казанский (Приволжский) федеральный университет, г. Казань, Россия.
