Например, применение двухциклического ускорения для i=5 осуществляется по алгоритму:
Надчеркнутые итерации
вычисляются по формуле (11), а ненадчеркнутые
,![]()
вычисляются по основному (т. е. ускоряемому) алгоритму.
Эксперименты показали, что выигрыш по времени имеет место для двухциклических ускорений с номерами
.
Надо иметь в виду, что невязки в граничных узлах при краевых условиях Дирихле заведомо равны нулю, при краевых условиях типа Неймана вычисляются по заданному алгоритму.
Для нелинейного оператора
, элементы которого зависят от номера итерации
,
, также возможна линейная комбинация (1) из предыдущих итераций. По аналогии с леммой 1 множители
вычисляются из системы
линейных уравнений
![]()
![]()
являющиеся необходимыми условиями существования экстре-мума функционала
,
где в связи с нелинейностью оператора
невязки вводятся следующим образом
![]()
Теорема 5. В двухциклических ускорениях (
)
коэффициенты
вычисляются по формулам
![]()
при этом имеет место неравенство
, т. е. невязки не возрастают.
Литература
1. . О циклической оптимизации итерационных про
цессов // Изв. АН КазССР, сер. физ.-мат. н. 1980, №5.
2. , . Итеративный процесс с
минимальными невязкими // Математический сборник, нов. серия
31(73), вып. 2, 1952, 315–334.


Доказательство переменности градиентов давления при смешении ламинарных струй в поперечном магнитном поле при числе Гартмана Ha=500. Здесь ![]()

Влияние 4-х точечных измерительных приборов на течение в канале
Глава 13. О ПРИМЕНЕНИИ СХЕМЫ ПРЕДИКТОР-КОРРЕКТОР В МЕХАНИКЕ СПЛОШНОЙ СРЕДЫ
В алгоритмах альтернативных направлений (типа Дугласа-Ракфорда, Писмана-Ракфорда, Марчука-Яненко и др.) в отли - чие от явных и неявных схем, возникает проблема постановки граничных условий для промежуточных сеточных функций. Эта проблема усугубляется для краевых условий типа Неймана, и становится неразрешимой в областях с криволинейной грани - цей. С этой точки зрения эффективным является применение ус - ловно устойчивых явных схем или абсолютно устойчивых неявных схем.
Для системы уравнений гидродинамики вязкой жидкости, в силу их нестандартности, идея применения неявных схем реали-
зуется в сочетании с полуявными схемами, технология констру - ирования которых изложена в главах 7, 8.
13.1. Предиктор-корректор в аэродинамике
Технология схемы предиктор-корректор демонстрируется для начально-краевой задачи:
(13.1.1)
![]()
![]()
Предиктором является монотонная схема:
![]()
(13.1.2)
![]()
![]()
(13.1.3)
![]()
(13.1.4)
(13.1.5)
![]()
![]()
,
где по главе 6
![]()
![]()
Реализация таких схем осуществляется глобальными итераци - ями, что приведено в главах 8 и 9. Итерационным методом из (13.1.3), (13.1.4), (13.1.5) вычисляется давление
, затем из (13.1.2) компоненты вектора скорости
на момент времени
.
Корректором является такая же монотонная, но неявная схе - ма на слое времени
. Корректор может состоять из двух видов.
Линейный корректор. Имеет вид:
(13.1.6)

![]()



![]()
![]()
![]()
В этой схеме известными сеточными величинами являются
,
,
и компоненты скорости
, они найдены по полуявной схеме (13.1.2). В (13.1.3) разностные производные взяты на момент времени
, а коэффициенты при них на слое
. В этом случае уравнения (13.1.3) являются системой линейных алгебраических уравнений относительно компонент скорости
и давления
на слое времени
. Причем система (13.1.3) имеет диагональное преобладание из-за членов вида
, т. е. выполняется условие сходимости итерационного метода Якоби.
|
Из за большого объема этот материал размещен на нескольких страницах:
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 |


