Главная > Разностные схемы
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
След.
Макеты страниц

Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше

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

ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO

ДОПОЛНЕНИЕ. МЕТОД ВНУТРЕННИХ ГРАНИЧНЫХ УСЛОВИЙ

В теории краевых задач для аналитических функций, т. е. для решений системы уравнений Коши — Римана, а также для решений более общих систем уравнений с частными производными, применяется метод сингулярных интегральных уравнений. Он состоит в сведении краевых задач к некоторым интегральным уравнениям на границе рассматриваемой области. При этом в дополнение к заданным граничным условиям используются следствия самой системы дифференциальных уравнений — соотношения, которым Должны удовлетворять функции (и их нормальные производные) на границе области, чтобы их можно было доопределить внутри области до некоторого решения соответствующей системы. В случае аналитических функций — это классическое условие Сохоцкого — Племеля, которое возникает при переходе в интегральной формуле Коши

к пределу при стремлении к границе у. В случае дифференциальных уравнений второго порядка соответствующее условие возникает из формулы Грина, выражающей решение в каждой точке области через значения этого решения и его нормальной производной на границе. Чтобы получить это условие, также надо перейти к пределу при стремлении точки изнутри области к ее границе, воспользовавшись свойствами потенциалов простого и двойного слоев.

Метод внутренних граничных условий по идее аналогичен описанному методу редукции краевых задач для уравнений с частными производными к интегральным уравнениям на границе. Роль дополнительных граничных условий, аналогичных условию Сохоцкого — Племеля, играют внутренние граничные условия, возникающие из разностного аналога интегральной формулы Коши (или разностного аналога формулы Грина).

1. Класс систем разностных уравнений. Рассматриваются краевые задачи для общих систем разностных уравнений с постоянными коэффициентами, которые в векторной записи имеют вид

где квадратные матрицы, - заданная и - искомая вектор-функции, К — конечное множество (шаблон). Будем Предполагать, что система (1) удовлетворяет следующему алгебраическому условию: характеристическая матрица

где - комплексные параметры, не является тождественно по вырожденной:

Это ограничение естественно: можно показать, что в случае уравнение (1) имеет решение не при всякой финитной (по ) правой части

2. Фундаментальное решение.

Матричную функцию назовем фундаментальным решением системы (1), если она однорременно удовлетворяет следующим двум уравнениям:

Лемма. Пусть есть произвольный многочлен от произвольного числа t комплексных аргументов, не обращающийся тождественно в нуль. Тогда можно выбрать радиусы окружностей так, чтобы выполнялось неравенство если

Доказательство проведем индукцией по числу аргументов t. При число корней конечно и утверждение очевидно. Считая, что утверждение доказано для , установим его в случае . Многочлен расположим по степеням

где M — некоторое натуральное число и не обращается тождественно в нуль. Выберем так, чтобы при Это возможно по предположению индукции. Выбирая теперь достаточно большим, можно добиться, чтобы при выполнялось неравенство ,

Теорема 1. Матрица определяемая равенством

является фундаментальным решением.

Здесь в соответствии с леммой выбраны так, чтобы , если

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

3. Граница сеточной области.

Рассмотрим уравнение (1) на некотором ограниченном множестве

где — произвольная сеточная область определения правой части Тогда область определения решения есть множество D, которое пробегает точка

если n и к пробегают независимо и К соответственно. Сопоставим каждому подмножество множества К, состоящее из всех тех , для которых Границей Г назовем совокупность всех тех точек для которых непусто. Например, для простейшего разностного аналога уравнения Пуассона

множество состоит из тех точек которые попали внутрь квадрата

Рис. 57.

Множество К — из пяти векторов Множество D — совокупность всех целочисленных точек квадрата кроме четырех угловых Граница Г состоит из двух слоев точек, отмеченных на рис. 57 крестиками.

4. Разностные аналоги интегральных формул Коши и типа Коши.

Лемма. Пусть — произвольная матрица-функция, для которой имеет смысл умножение справа на квадратную матрицу порядка , определенная на всей целочисленной сетке. Справедливо следующее тождество:

Доказательство. Вектор-функцию можно записать в виде

Левая и правая части тождества (7) линейно зависят от и. Поэтому для доказательства достаточно проверить справедливость тождества (7) для вектор-функции

при каждом фиксированном

Теорема 2. Пусть , — произвольное решение уравнения произвольное фундаментальное решение. Тогда справедлива формула

Доказательство. Помножим обе части равенства (6) слева на матрицу и просуммируем по всем . Воспользовавшись тождеством (7), а затем равенством (4), получим формулу (8).

Следствие. Каждое решение уравнения (6) полностью определяется своими значениями на Г и восстанавливается по этим, значениям по формуле (8).

Теорема 3. Пусть произвольная вектор-функция размерности , определенная на Г, и пусть произвольное фундаментальное решение. Тогда формула

задает некоторое решение уравнения (6).

Доказательство. Применим оператор L к вектор-функции определенной формулой (9):

Вычислим правую часть. В силу (4) имеем

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

Равенство (8) аналогично интегральной формуле Коши для аналитических функций в ограниченной области d с границей у:

При этом роль аналитических функций, границы области и ядра Коши - играют соответственно решения задачи (6), граница точной области D и выражение учитывающее через множество по которому ведется суммирование, структуру границы вблизи точки

Формулу (9) в таком случае естественно сравнить с интегральной формулой типа Коши. Формула (8) аналогична также формуле Грина для уравнения Лапласа.

Подчеркнем, однако, следующее существенное различие между формулами (11) и (8): интегральная формула Коши справедлива только строго внутри области d, а разностная формула (8)- всюду на D, включая точки границы Г.

Аналогичное различие имеется также между формулой (9) и формулой Грина.

5. Внутренние граничные условия.

Теорема 4. Пусть - какое-нибудь фундаментальное решение уравнения (1). Для того чтобы заданную на Г вектор-функцию можно было доопределить всюду в ограниченной сеточной области D до некоторого решения уравнения (6), необходимо и достаточно, Чтобы при всех выполнялись равенства

Доказательство. Если можно доопределить всюду на D до некоторого решения , то, применив к этому решению формулу (8), а затем рассматривая полученное равенство только при , убедимся в выполнении (12). Обратно, если удовлетворяет (12), то примем и построим некоторое решение по формуле (9). В силу (12) граничные значения этого решения совпадут с заданными.

Доказанная теорема 4 дает основание назвать равенства (12) внутренними граничными условиями: эти условия не задаются извне, а являются следствиями самого разностного уравнения.

Если формулы (8) и (9) трактовать как аналоги интегральных формул Коши и типа Коши, то внутренние граничные условия аналогичны классическим условиям Сохоцкого — Племеля, при которых заданную на границе у области d на комплексной плоскости функцию можно доопределить всюду в области d до некоторой аналитической функции.

Формулу (8) можно понимать как адекватную системе (6) разностную формулу Грина, которая неявно учитывает «скачки потенциалов» на границе Г и приводит к внутренним граничным условиям (12).

6. Оператор граничного проектирования.

Возможна отличная от (12) запись внутренних граничных условий. Будем обозначать через линейное пространство всех сеточных вектор-функций , а через подпространство тех из них, которые можно доопределить всюду в D до решений однородного уравнения, соответствующего уравнению (6).

Определим линейное отображение пространства UF в себя следующей формулой:

Теорема 5. Оператор Р есть оператор проектирования на

Доказательство. Действительно, при любом в силу теоремы 3 элемент принадлежит Если то в силу теоремы 2 получим Теорема доказана.

Оператор Р, определенный формулой (13), будем называть оператором граничного проектирования. С его помощью внутренние граничные условия (12) в случае записываются в форме

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

7. Общая краевая задача.

В силу следствия из теоремы 2 каждое решение уравнения (6) восстанавливается по его значениям на границе Г. Это

дает основание определить общую линейную краевую задачу для уравнения (6) как краевую задачу вида

где какой-нибудь линейный оператор, отображающий пространство на некоторое линейное пространство Ф.

Естественные разностные схемы, аппроксимирующие первую, вторую или третью краевые задачи для уравнения Пуассона, например, легко записать в виде (15).

Название «общая краевая задача» несколько условно: могут встретиться разностные краевые задачи, имеющие иной, чем (15), вид. Например, таковы естественные разностные схемы для дифференциальных краевых задач, в которых порядок дифференциального уравнения ниже порядка дифференциальных краевых условий.

8. Основная идея метода внутренних граничных условий.

Пусть для простоты Между разностной краевой задачей

и между задачей

существует тесная связь. Именно, граничные значения каждого решения , задачи (16) удовлетворяют в силу теоремы 2 уравнениям (17). Обратно, каждое решение задачи (17) в силу теоремы 4 и следствия из теоремы 2 можно единственным образом доопределить всюду в D до решения задачи (16). Основная идея метода внутренних граничных условий состоит в переходе от исходной разностной краевой задачи (16) к системе уравнений (17) на границе Г. Продвижения, которые при этом удается получить, основаны на двух обстоятельствах. Первое из них — малое по сравнению с задачей (16) число неизвестных, участвующих в задаче (17). Второе — специальный вид системы (17), в структуру которой органически включен оператор граничного проектирования.

9. Устойчивость внутренних граничных условий.

Можно опасаться, что внутренние граничные условия «почти вырождены», и поэтому задача (17) плохо обусловлена независимо от вида оператора так что переход от задачи (16) к задаче (17) связан с потерей вычислительной устойчивости.

Будем считать, что пространство Ф вложено в пространство UF, введем в пространстве значит, и в пространстве Ф с: ) норму IMI и докажем теорему, означающую, что при переходе от задачи (16) к задаче (17) не происходит потери вычислительной устойчивости.

Теорема 6. Пусть задача (17) имеет решение при любом причем выполнена оценка

где с от не зависит. Пусть, далее, произвольный элемент из Обозначим

Тогда справедлива оценка

Если рассматривать (19) как уравнения для определения , то оценка (20) означает, что чувствительность решения задачи (19) к возмущениям правой части внутренних граничных условий характеризуется постоянной с из оценки (18), т. е. чувствительностью решения к возмущениям правой части заданного граничного условия .

Доказательство. Обозначим

В силу теоремы 5 имеем . Поэтому

т. e. удовлетворяет системе вида силу (18) оценке

Отсюда с учетом тождества следует (20).

10. Дополнительная идея.

Изложим полезную при численном решении краевых задач для уравнений с частными производными идею применительно к следующей задаче.

Пусть функция определена в некоторой области d с достаточно гладкой границей у как решение задачи Дирихле

и требуется найти производную

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

Под s понимается длина дуги вдоль границы у. причем будем считать для определенности, что полная длина границы у есть . Функцию

будем искать приближено в виде частичной суммы

ее ряда Фурье. Для определения коэффициентов воспользуемся методом внутренних граничных условий.

Зададим построим сетку

и разностное уравнение

Отнесем к все те точки сетки, которые вместе со всеми четырьмя

соседними точками принадлежат Тогда определится сеточная область ее граница и внутренние граничные условия Идея состоит в том, чтобы по функции и фукнции

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

из условия минимизации невязки, возникающей при подстановке продолженной с границы Y в приграничную полоску функции во внутренние граничные условия.

11. Сопоставление метода внутренних граничных условий с методом сингулярных интегральных уравнений.

В начале Дополнения мы указывали на аналогию между методом внутренних граничных условий и методом сингулярных интегральных уравнений, которая не является полной. Здесь мы сопоставим эти методы, уточняя аналогию и выявляя существенные различия.

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

где — ограниченная область, ее граница. Краевое условие (22) связывает решение на границе области и его производную по направлению внутренней нормали Коэффициенты заданные операторы.

Выпишем классическую формулу Грина для уравнения (21):

где - фундаментальное решение уравнения (21), стремящееся к нулю на бесконечности. Устремим к границе у. Воспользовавшись свойствами потенциалов простого и двойного слоев, получим на границе Y соотношение вида

связывающее решение и его нормальную производную на границе области; и некоторые известные интегральные операторы. Переход от задачи (21), (22) к равносильной системе уравнений (22), (24) относительно функций определенных на границе и составляет сущность метода сингулярных интегральных уравнений.

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

Запишем внутренние граничные условия в удобной для дальнейшего форме.

Легко проверить, что формула (9) в этом случае может быть переписана в форме

где - совокупность точек Г, лежащих на сторонах квадрата т. е. на внешнем слое двухслойной границы Г сеточной квадратной области (рис. 57), а разностный аналог производной по направлению внутренней нормали.

Заметим, что формула (27) была бы полным аналогом классической формулы Грина (23), если бы в ее правой части отсутствовало «сингулярное слагаемое» Ьгпиг Однако в таком хлучае равенство (27) имело бы место не при всех , а лишь при Из него нельзя было бы получить тогда внутренние граничные условия Эти условия получаются из (27), если пробегает не всю область D, а только точки границы Г, и записываются двумя системами равенств

отвечающих, соответственно, точкам внешнего и точкам внутреннего слоев двойной границы Г. В качестве в равенствах (28) и (29) будем использовать ограниченное фундаментальное решение.

Можно показать, что внутренние граничные условия т. е. система уравнений (28), (29), алгебраически равносильна каждой из отдельно взятых подсистем (28) или (29).

Подсистема (28) аналогична интегральному соотношению (24), так что разностным аналогом задачи (22), (24) является задача (27), (28), но не задача

записываемая равенствами которая рассматривается в методе внутренних граничных условий.

Имеется очевидная разница между внутренними граничными условиями т. е. системой (28), (29), и одной только подсистемой (28). Внутренние граничные условия

содержат избыточные равенства (29). В этом смысле разностные внутренние граничные условия

больше похожи не на интегральное соотношение (24), а на условия Сохоцкого — Племеля для аналитических функций. Эти последние представляют собой два вещественных соотношения, связывающих две вещественных функции, но они не независимы, и многообразие удовлетворяющих условиям Сохоцкого — Племеля пар функций зависит от одной произвольной вещественной функции.

Отметим, что внутренние граничные условия

выгодно отличаются от равносильной им подсистемы (28) тем, что в их структуру входит оператор граничного проектирования. Благодаря этому обстоятельству задача

устойчива в смысле теоремы 6 относительно возмущения правой части . В общем случае при вычеркивании части уравнений из числа составляющих систему может получиться подсистема, алгебраически равносильная исходной, но уже не обладающая свойством устойчивости.

Можно показать, что в нашем примере (25), (26) вместо

удобнее использовать не подсистему (28), а подсистему (29), которая устойчива и в отличие от подсистемы (28) состоит из независимых уравнений — ее ранг равен числу составляющих ее уравнений.

Итак, в рассматриваемом примере аналогия между методом внутренних граничных условий и методом сингулярных интегральных уравнений, аналогичных условию Сохоцкого — Племеля, не является полной.

Тем более, нет полной аналогии с классическим методом интегральных уравнений, в котором искомой функцией является не само решение исходной задачи (21), (22) на границе, а некоторая вспомогательная плотность потенциала простого или двойного слоя.

В заключение заметим, что мы употребляем выражение «метод сингулярных интегральных уравнений», поскольку условие Сохоцкого — Племеля содержит сингулярный интеграл. В примере этого пункта условие (26) содержит сходящиеся несобственные интегралы.

1
Оглавление
email@scask.ru