книги из ГПНТБ / Годунов С.К. Разностные схемы. Введение в теорию учеб. пособие
.pdf360 |
ГЛ. 12. КОНСТРУКЦИЯ ОПЕРАТОРА ПЕРЕХОДА |
достаточно большую квадратную сеточную область и рассмотрим оператор Л ^ + Л т о над функциями из U'h, определенными на сеточном квадрате и обращающимися в нуль на его границе.
Доопределим каждую функцию и' из U'n до некоторой функ ции и" из U'h, положив ее тождественно равной нулю во всех тех точках сеточного квадрата, которые не принадлежат перво начальной области. Легко видеть, что для каждой такой функ ции благодаря ее обращению в нуль на границе первоначальной области выполнено равенство
(Rhu', |
и') |
_ |
{Rhu", |
и") |
п |
_ |
д |
, д |
(и', |
и') |
~ |
(и", |
и") |
' * H |
~ L |
Y X X |
1УУУ |
Поэтому при переходе от формул (18) к формулам
Пт", |
и") |
min ^*,/ \ , , - } — X m i n , |
|
(и", |
и") |
"& vh
ш а х |
t,,rr |
,,n\ |
Лщах |
|
|
„ |
,, |
\U , |
U ) |
|
|
получим числа Хтт и Ятах, |
которые удовлетворяют оценкам |
||||
. ^min ^ ^min> |
Ятах ^ Ятах> |
(19) |
|||
Но в случае квадратной области собственные значения |
известны, |
||||
так что k'min и /Сах известны, |
и мы получили оценки |
(19) гра |
|||
ниц спектра оператора Л ж ж + Л т а |
над функциями из U'h, опреде |
||||
ленными в первоначальной |
области. |
|
|||
Во многих случаях можно применять вариационные методы оценки собственных значений, аналогичные вариационным ме тодам для дифференциальных уравнений. Например, первое
собственное значение задачи |
|
|
||
|
|
Ли№>, |
«(«Lл= 0 , |
|
где Г„ — граница сеточной |
области Dh и во внутренних |
точках |
||
задано выражение |
|
|
|
|
Rhumn — "д2"{ а |
"Ь "2 ' У") (Um+1. « ~*" итп) |
|
||
— a(xm--J> |
Уп) (Umn — " т - 1 . п) } + |
|
||
+ |
- j - j Ь (хт, |
уп + -j) («т. л+1 — «ли) — |
|
|
|
|
|
b [хт, уп — -j) (итп — Um, |
| , |
|
а (х, у)>0, |
b (х, у) > 0, |
|
|
362 |
ГЛ. 12. КОНСТРУКЦИЯ ОПЕРАТОРА |
ПЕРЕХОДА |
|
и что |
расположение всех |
собственных |
чисел Xh на отрезке |
— 1 ^ |
X ^ 1 необходимо и |
достаточно |
для выполнения нера |
венств |
|
|
|
|
1 |
№ й < 1 . |
(23) |
где норма оператора задается с помощью скалярного умноже
ния |
(21). |
|
|
|
6. |
Критерии |
устойчивости |
Самарского. |
В предложенной |
А. А. Самарским |
[17], [18] теории |
устойчивости |
широкого класса |
|
разностных схем в гильбертовом пространстве указаны необхо димые и достаточные условия устойчивости в терминах линей
ных неравенств между |
операторными коэффициентами этих |
схем, а также получены |
другие результаты. |
Приведем здесь лишь два результата из этой теории.
Пусть U'h — евклидово пространство с некоторым скалярным
умножением |
(и, v) == [и, v], и пусть оператор RH, V = Rhu, |
u , t i s |
s U'h, задан |
равенством |
|
|
B n ^ - + Anu=0, |
(24) |
где Ah и Bh — самосопряженные операторы, причем Bh > 0. Определим энергетическую норм$ || и \\в^ в пространстве U'h,
положив
|
|
|
II « t |
= |
[f i A"> |
и] |
= (". |
")в |
• |
|
(26) |
|
|
|
|
|
h |
|
|
|
|
h |
|
|
|
Тогда справедлива |
следующая |
|
|
|
|
|
|
|||||
Т е о р е м а |
2. |
Условия |
|
|
|
|
|
|
|
|
||
|
|
|
|
0 < Л й < - | я л |
|
|
|
(26) |
||||
необходимы и достаточны |
для |
выполнения |
|
неравенств |
|
|
||||||
|
|
|
\Rl\<\, |
|
Р>0. |
|
|
(27) |
||||
Д о к а з а т е л ь с т в о . |
Определим |
самосопряженный |
опера |
|||||||||
тор Ah равенством |
Ah = Bh — xAh. |
Тогда |
(24) равносильно |
ра |
||||||||
венству |
(20), |
а условия |
(26) |
равносильны |
условиям |
—Bh |
^ |
|||||
Ah ^ |
Вн, т. е. условиям |
|
|
|
|
|
|
|
|
|||
|
|
— [Bhu, |
ы ] < [ Л л ы , |
« ] < [ В л и , |
и]. |
(28) |
||||||
Как показано в п. 5, оператор Rh является самосопряжен ным в смысле скалярного умножения (21) и утверждение тео ремы равносильно утверждению, что все собственные числа Xh оператора Rh лежат на отрезке —1 Я ^ 1 в том и только том случае, если выполнены условия (28). Докажем это по следнее утверждение.
ЗАДАЧИ |
363 |
Пусть выполнены условия (28). Умножим равенство (22) скалярно на собственную функцию грМ оператора Rh. Получим
**UWW. *W J-[^+W , Л откуда
[ВН^\ *<*>] < |
1. |
|
Обратно, пусть m a x | A * | ^ 1. Покажем, |
что выполнены |
усло |
вия (28). Пусть и = 2 ск$к)—разложение |
произвольного |
эле |
мента u e t / J по ортонормальному в смысле скалярного умно
жения |
(21) |
базису |
{i|>< f t ) }. |
Тогда |
|
|
|
|
|
|
|
||
| [Апи, |
и] | = |
| [~Ah2 |
|
2 с*ЛI = I [2 ckAh^k\ |
2 <*Н I = |
||||||||
-1 [2 cft |
|
2 с»ЛI = I [Bh 2 c |
* M w . |
2 <*i>w] I = |
|||||||||
|
|
= |(2*A+«. 2^>)J==24N< |
|
|
|
||||||||
|
|
|
|
|
|
|
<2с| |
= (ы, w)B |
= [ B f |
t u , и]. |
|||
Отсюда |
[BhU,u]^\[AhU,u]\, |
что равносильно |
условиям |
(28). |
|||||||||
Теорема |
доказана. |
|
|
|
|
|
|
|
|
|
|||
Заметим, |
что проверка |
условий |
(28) равносильна |
проверке |
|||||||||
того, будут |
ли |
неотрицательны все собственные числа |
самосо |
||||||||||
пряженных |
в смысле |
скалярного |
умножения |
[и, v] |
операторов |
||||||||
Bh — Ah |
и Bh |
+ Ah. |
|
|
|
|
|
|
|
|
|
|
|
Приведем без доказательства еще один критерий устойчивости, приме |
|||||||||||||
нимый к разностным схемам (24), для которых B f t |
> 0 , Ап |
= Ап |
>0. |
|
Введем |
||||||||
в пространстве U'h |
энергетическую норму |
|| и \\, |
, |
положив || и [\2 |
|
=^[Апи,и]. |
|||||||
Т е о р е м а |
3. |
Выполнение |
условия Вп~^ -g- Л„ |
необходимо |
и достаточно |
||||||||
для того, чтобы имело |
место неравенство |
|| Rh II д |
<^ 1. |
|
|
|
|
||||||
|
|
|
|
|
|
|
"h |
|
|
|
|
|
|
Теорема 3 содержится |
в п. 4 § 1 гл. VI книги [17] и доказывается без |
||||||||||||
помощи спектрального подхода, который здесь не удается применить из-за несамосопряженности оператора Вн.
ЗА Д А Ч И
1.Пусть оператор b = Rna, задан формулами
bm^V — г) am + ram+1, |
m = 0, 1 |
М — 1 , "> |
|
& ж = о . |
|
; |
|
Показать, что в пространстве lfh сеточных функций {a m } , m = 0, 1, |
М, |
||
нельзя задать скалярное произведение |
так, чтобы |
оператор Rh стал |
само |
сопряженным. |
|
|
|
Г Л А В А 13
СПЕКТРАЛЬНЫЙ ПРИЗНАК УСТОЙЧИВОСТИ НЕСАМОСОПРЯЖЕННЫХ ЭВОЛЮЦИОННЫХ КРАЕВЫХ ЗАДАЧ
Здесь мы покажем, что по спектру несамосопряженного опе ратора Rh нельзя судить об устойчивости разностной краевой задачи в ограниченной области, введем понятие спектра семей ства операторов {Rn} и рассмотрим спектральную постановку вопроса об устойчивости, остающуюся разумной и в случае не самосопряженных разностных краевых задач в ограниченных об ластях. Будет указан необходимый и близкий к достаточному спектральный признак устойчивости.
§42. Спектр семейства операторов {Rn}
1.Необходимость усовершенствования спектрального при знака устойчивости. В гл. 12 было показано, что обычно эволю
ционные разностные краевые задачи можно привести к виду ыр+i = RhUp + трр,
ц° задано
так, чтобы устойчивость на интервале времени 0 ^ t ^ Т была равносильна равномерной по h ограниченности норм степеней оператора перехода Rh, т. е. оценке
|
\\Rph\\<K, |
р = 1 , |
2, . . . , [ 7 / т ] , |
(2) |
||
где т — шаг сетки по |
времени, т = |
т(л). |
|
|
||
|
Было установлено, что расположение собственных значений |
|||||
оператора Rh внутри круга |
|
|
|
|
||
|
|
U |
К 1 + сх |
|
(3) |
|
на |
комплексной плоскости |
необходимо для |
выполнения усло |
|||
вия |
(1), т. е. для устойчивости. В § 41 было показано, что в |
слу |
||||
чае |
самосопряженного |
оператора |
Rh условие |
(3) является |
не |
|
S 42, СПЕКТР СЕМЕЙСТВА ОПЕРАТОРОВ |
365 |
только необходимым, но и достаточным условием |
равномерной |
ограниченности (2) норм степеней оператора Rh. |
Этот же факт |
установлен в § 25 и для разностной задачи Коши с постоян ными коэффициентами в случае двуслойных разностных схем относительно одной неизвестной функции независимо от того, имеет ли место самосопряженность. Однако в общем случае не самосопряженных разностных краевых задач в ограниченных областях необходимый признак (3) очень далек от достаточного и совершенно неадекватен вопросу о равномерной ограниченно
сти |
(2) |
норм |
I Rh\ |
степеней |
оператора |
Rh- |
|
|
||||
|
Это показывает |
следующий |
пример. |
|
|
|||||||
|
П р и м е р . |
Для |
разностной |
краевой |
задачи |
|
||||||
|
|
„р+1 |
_ |
»Р |
*т+1 |
|
= 4>(хт, |
tp), |
|
|||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
0,1, ...... |
|
|
(4) |
|
|
|
и°т = У(хт)> |
т = |
М; |
Mh=\, |
|||||||
аппроксимирующей |
задачу |
|
|
|
|
|
||||||
|
|
|
|
|
— их = |
у(х, |
t), |
|
|
|
||
|
|
|
|
и(1, |
t) = |
0, |
0 < * < 1, |
0 < / < |
Г, |
|||
|
|
|
и (х, |
0) = |
(х) |
|
|
|
|
|||
при естественном приведении к каноническому |
виду (1) опера |
|||||||||||
тор |
Rh, |
v = |
Rnu, |
задается |
формулами |
|
|
|
||||
|
|
vm = |
(l—r)um |
+ rum+u |
т = |
0, |
1, .. |
|
||||
vM = Q, r = x//i. Его матрица имеет вид
1 - г |
г |
0 . . . |
0 |
0) |
|
0 |
l—r |
г . . . |
0 |
о |
|
|
|
|
|
|
(5) |
о |
о |
о |
|
1-г |
г |
о |
о |
о |
|
0 |
0J |
Спектр матрицы состоит из ее |
собственных |
значений, т. е. из |
|||
корней уравнения |
|
|
|
|
|
det (Rh — IE) = 0 или |
- |
к (1 - |
г — Х)м = 0. |
||
Корни этого уравнения |
Х = |
0 |
и Я = |
1 — г образуют спектр |
|
оператора Rh при любом h. Этот спектр |
лежит в единичном |
|
круге |
при 0 _ < / - < 2 . Между тем |
для схемы (4) при |
366 ГЛ. 13. УСТОЙЧИВОСТЬ НЕСАМОСОПРЯЖЕННЫХ ЗАДАЧ
1 < |
г ^ |
2 не выполнено условие Куранта, |
Фридрихса |
и |
Леви, |
||||||||||||||
так |
что |
устойчивости |
||/?А||</ С |
ни |
в |
какой |
разумной |
норме |
|||||||||||
быть не |
может. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
В самом |
деле, |
покажем, |
что |
в случае |
Г > 1 |
и |
нормы II и II = |
max | ит | |
||||||||||
справедливо |
неравенство |
|
|
|
|
|
|
|
|
|
|
|
т |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
Р=1. |
max |
|
\Rl\>\ |
|
1 - 2 г | 1 / л |
= р , / \ |
|
|
(6) |
|||||
|
|
|
|
|
2 |
|
[Г/т]" |
я " |
|
|
|
|
|
|
|
|
|
||
При |
г > 1 |
также р > 1 , |
так |
что при |
h->Q |
и т = |
гА - >0 |
величина |
max | / ? ^ | |
||||||||||
экспоненциально |
возрастает |
и |
условие |
|
|
I |
К |
грубо |
|
Р |
|
||||||||
|
| R ^ || < |
нарушается. |
|||||||||||||||||
Для |
доказательства |
неравенства |
(6) |
заметим, |
что |
в |
случае |
ит |
— |
(—\)т, |
|||||||||
т = |
0, 1, |
|
М, значения и?т функции |
|
|
|
|
|
|
|
|
|
|||||||
|
|
и" = Rphu°, |
р = |
1, |
2, . . . . |
М, |
при |
nt = |
О, I |
|
Af - |
р |
|
|
|||||
задаются |
формулой |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
< = ( - 1 ) т ( 1 - 2 г ) " , ' т = 0, 1 |
|
М - р . |
|
|
|
|||||||||||
Поэтому |
|
|
|
Jj ?gB0 [>|l-2r|"||«°IL |
р<М, |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
так |
что при |
этих |
значениях |
р р = 1, 2, |
|
|
М, |
|
|
|
|
|
|
||||||
|
|
|
| | / ? Р | | > | 1 - 2 л | р , |
р = 1 , 2 , |
|
A f - 1 / А , |
|
|
|
||||||||||
и неравенство (6) |
доказано. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Итак, |
установлено, |
что |
необходимый |
спектральный |
признак |
|||||||||||||
(3) |
равномерной |
ограниченности || Rhl |
< |
К, |
использующий |
соб |
|||||||||||||
ственные значения операторов Rn, слишком груб в случае неса мосопряженных операторов Rh: в нашем примере он не улавли вает неустойчивость, имеющую место при 1 < г 2.
2. Определение спектра семейства операторов. Пусть линей' ный оператор Rh определен на линейном нормированном про странстве U'h-
Будем обозначать через {Rh} совокупность операторов Rh при всех тех значениях h, которые принимает параметр h, характери зующий густоту сетки. По самой природе разностных схем шаг сетки h может принимать сколь угодно малые положительные значения.
Комплексное число К будем называть точкой спектра семей ства операторов {Rh}, если для любых положительных h0 и в можно указать такое h, h •< ho, что неравенство
| | Д А и - Я и | 1 < в | | и Ц
{шеет некоторое решение и, к е U'h.
|
|
§ 42. СПЕКТР СЕМЕЙСТВА ОПЕРАТОРОВ |
|
|
367 |
||||
Совокупность всех |
таких |
чисел К будем |
называть |
спектром |
|||||
семейства |
операторов {Rh}- |
|
|
|
|
|
|
||
3. Необходимое условие устойчивости. |
|
|
|
|
|||||
Т е о р е м а |
1. Пусть хотя бы одна точка Ко спектра |
семейства |
|||||||
операторов |
{Rn} лежит |
вне единичного |
круга |
комплексной |
пло |
||||
скости, так что | Я о | > 1. В таком случае |
нельзя указать |
общую |
|||||||
для всех h постоянную |
К такую, чтобы выполнялось |
неравенство |
|||||||
|
|
|
\Rl\<K, |
|
|
|
|
(7) |
|
в котором |
р |
пробегает |
целые |
значения |
от 0 до p = |
p0(h), |
где |
||
po{h) —* 0 0 |
П Р И п ~* О- |
|
|
|
|
|
|
|
|
Д о к а з а т е л ь с т в о . Допустим сначала, |
что не существует |
||||||||
чисел Но >_ 0 и с > 0 таких, |
что при всех h < Н0 |
справедлива |
|||||||
оценка |
|
|
Ц Я * 1 1 < с . |
|
|
|
|
(8 ) |
|
|
|
|
|
|
|
|
|||
При этом допущении доказываемое утверждение очевидно. По
этому |
остается рассмотреть |
случай, когда существуют Н0 > |
0 и |
|||||
с > |
0 |
такие, |
что |
при Н < |
/г0 неравенство |
(8) справедливо. |
||
|
Положим |
\\о\ |
= 1 + 6 , где Ко— та точка |
спектра, для кото |
||||
рой |
|Аю|> 1. Задавшись |
произвольно числом К, выберем |
рае |
|||||
так, чтобы выполнялись |
неравенства |
|
|
|||||
( l + 6 ) p > 2 t f ,
1 - ( 1 + с + с 2 + . . . + С - 1 ) 8 > | .
По определению точки спектра семейства операторов {Rh}, можно указать сколь угодно малые положительные Н, при кото рых существует вектор и е £//,, являющийся решением нера венства
Положим |
\\Rhu-k0u\\<e\\u\\. |
(9) |
||
|
Rhu = %oU + z. |
(10) |
||
|
|
|||
Ясно, что | | 2 | | < е | | « | | . Далее, |
из (10) можно вывести, что |
|||
^ ы - Я о Р |
« + (ЯоР - 1 2 + Я о р |
- 2 ^ + . . . |
+RVz). |
|
Поскольку | Я01 > |
1, то |
|
|
|
W^-'z+K-'RhZ-r- |
. . . |
+/?Г'г||< |
|
|
а следовательно, |
< U 0 | p |
( l + | | t f A | | + IUy+ |
... +1яГ , !)в||и||, |
|
|
|
|
|
|
| | ^ " | > U 0 | p [ l - e ( l + c + C 2 + . . . + с р - 1 ) ] | | « | | > |
||||
|
|
> ( 1 + |
б)" J I I и || > 2К111"11 e К 1 1 и И- |
|
368 |
ГЛ. 13. УСТОЙЧИВОСТЬ НЕСАМОСОПРЯЖЕННЫХ ЗАДАЧ |
Число h во всем этом построении можно считать настолько малым, чтобы р было меньше, чем Po(h).
Ввиду произвольности К доказано наше утверждение о том,
что |
расположение |
всех точек спектра |
семейства |
операторов |
|
{Rh} |
внутри или на |
границе |
единичного |
круга \Х\ ^ |
1 необхо |
димо |
для выполнения |
оценки |
|| Rn || < К. |
|
|
4. Обсуждение понятия спектра семейства операторов {Rh).
Начнем с того, что обратим внимание на аналогию между опре делением точки спектра семейства операторов {Rh} и следующим определением точки спектра какого-либо оператора R, которое приводится в курсах функционального анализа. Мы будем в ка честве R брать оператор Rh при некотором фиксированном h.
Точка X на комплексной плоскости называется точкой спек тра оператора Rh, если при любом положительном е неравенство
|
|
||/? А и - А,и||<в||и|| |
|
|
|
имеет |
решение и, |
принадлежащее пространству Un, на котором |
|||
определен оператор Rh. |
|
|
|
||
При сравнении |
определений точки спектра семейства |
опера |
|||
торов |
{Rh} и точки спектра |
оператора Rh может |
возникнуть |
||
мысль, что спектр |
семейства |
{Rh} состоит из тех точек |
комплекс |
||
ной плоскости, которые получаются путем предельного перехода при /г-* О из точек спектра оператора Rh, когда И—+0 по все возможным подпоследовательностям. Но, вообще говоря, это
предположение |
ошибочно. |
|
|
|
|
|
|
|
|
|
|
||||||
Рассмотрим |
оператор |
Rh, v = |
RhU, задаваемый |
равенствами |
|||||||||||||
|
fm |
= |
(1 — г) um + |
|
rum+l, |
|
m = 0, |
1,..., M - l , |
J |
|
|
||||||
|
vM |
= |
0, |
|
|
|
|
|
|
|
|
Mh = |
\. |
J |
( |
' |
|
Оператор (11) действует в (M -f- 1) -мерном линейном про |
|||||||||||||||||
странстве |
и |
описывается |
матрицей |
(5). Известно, что спектр |
|||||||||||||
матрицы состоит из ее собственных значений, т. е. из корней X |
|||||||||||||||||
уравнения |
det (/?л — ХЕ) = |
0. |
Мы |
вычислили эти |
собственные |
||||||||||||
значения |
в п. |
1. Это л. = |
0 |
и А = |
1 — г . |
Таким |
образом, |
спектр |
|||||||||
оператора |
Rh |
при любом |
h |
состоит |
из двух |
точек |
X = |
0 и |
|||||||||
X = 1—г, |
не зависящих от h. Однако спектр семейства |
опера |
|||||||||||||||
торов {Rh}, как будет показано в |
§ 43, состоит не только из |
||||||||||||||||
этих двух точек, чего, казалось бы, можно было ожидать, |
а |
еще |
|||||||||||||||
и из |
всех |
точек |
круга |
\Х— |
\-\-г\<.г |
|
радиуса г с центром в |
||||||||||
точке |
X = |
1 — |
г |
(рис. 27, |
стр. 248). |
При г sg; 1 спектр |
семейства |
||||||||||
операторов |
{Rh} лежит |
в единичном |
круге |Я| ^ |
1, а |
при |
r~> 1 |
|||||||||||
этот необходимый спектральный признак устойчивости не вы полнен: неравенство \ Rph\ < ^ н е может выполняться равномерно ро h.
