19266

Организация итерационного процесса. Проблемы сходимости численных схем. Улучшенные итерационные методы. Внутренние и внешние итерации

Лекция

Энергетика

Лекция 14. Организация итерационного процесса. Проблемы сходимости численных схем. Улучшенные итерационные методы. Внутренние и внешние итерации. 14.1. Прямой метод решения уравнений в матричной форме. Систему конечноразностных уравнений записанную в матричной

Русский

2013-07-11

89.5 KB

11 чел.

Лекция 14.

«Организация итерационного процесса. Проблемы сходимости численных схем. Улучшенные итерационные методы. Внутренние и внешние итерации.»

14.1. Прямой метод решения уравнений в матричной форме.

Систему конечно-разностных уравнений, записанную в матричной форме, можно решить. Для этого необходимо ввести удобный порядок членов , так чтобы прямоугольную систему {k, т, g} можно было представить вектором . Очевидный выбор порядка членов состоит в том, что нумерация начинается с нижнего левого угла и производится по рядам. Все граничные точки исключаются с помощью граничных условий, например: (k, т, g) = 0, если k=0 или К. Для описания компонент вектора  используется единственный индекс:  j = 1,2, … , - l) М G.

Определив порядок компонент вектора , обозначив  =  +  –  – , рассмотрим уравнение в матричной форме в виде:

= .

Диагональные компоненты матрицы  положительны, в то время как недиагональные члены - отрицательны или равны нулю. Сумма недиагональных элементов в любом данном ряду меньше, чем диагональный элемент. Таким образом, матрица  является неприводимой диагонально преобладающей. Следовательно, для  существует обратная матрица , и решение уравнения можно записать в виде:

=

14.2. Организация итерационного процесса.

Прямые методы обращения матрицы  весьма громоздки, поэтому используют итерационные методы решения уравнения. Чтобы понять основные принципы, запишем матрицу  в виде суммы трех матриц:

=

где  – диагональная .матрица (отличные от нуля элементы находятся только на основной диагонали),  – верхняя треугольная матрица (отличные от нуля элементы находятся только выше основной диагонали) и  – нижняя треугольная матрица (отличные от нуля элементы находятся ниже основной диагонали).

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

= ( + ) + .

умножить обе части уравнения на  – матрицу, обратную , такую, что произведение этих матриц равно единичной матрице 1. Так как  представляет собой диагональную матрицу, то каждый элемент матрицы  равен обратной величине соответствующего элемента матрицы . Итерационный процесс можно определить следующим образом:

=  ( + ) + .

14.3. Проблемы сходимости численных схем.

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

14.4. Улучшенные итерационные методы.

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

( – ) =  + .

Так как матрица ( – ) треугольная, включая основную диагональ, то можно легко найти обратную ей или решить уравнение относительно .

14.5. Внутренние и внешние итерации.

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

PAGE  2


 

А также другие работы, которые могут Вас заинтересовать

35644. Творческий проект «Тигр». В технике холодный батик 246.84 KB
  Что означает слово батик. Влияние батика на новое поколение. На уроках технологии я научилась рисовать батиком.
35645. Бижутерия. Творческий проект 64.68 KB
  Бисер - (от араб. бусра – фальшивый жемчуг) - мелкие разноцветные стеклянные бусины с отверстиями, употребляемые главным образом для вышивания (на предметах быта, женской одежде и т. д.). Бисер плотно вошел в нашу жизнь, заявив о себе в разных сферах жизни человека. Бисер можно использовать как элемент декора интерьера, как эксклюзивное украшение, для некоторых людей бисер является неотъемлемой частью их жизни. Маленькие разноцветные шарики с отверстиями несут в себе магическую силу притяжения. И столкнувшись один раз с бисером, вы оказываетесь под его влиянием.
35648. Пошив фартука. Творческий проект 468.57 KB
  Техническая последовательность Сначала мы сняли мерки потом начертили части фартука размером 1к4 в тетради в настоящую величину специальными булавками прикрепили на ткань. Скоро мы сшили карманы и низ фартука.
35649. Вязаная сумочка под различные мелочи 800.13 KB
  Затраты на изделие=стоимость ниток стоимость Крючка Стоимость всей пряжи:60 рублей. Стоимость крючка:10 рублей.
35650. Краснодара Творческий проект Би. 22.55 KB
  170 руб. 2 200 руб. 3 220 руб. 4 250 руб.