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


 

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

62397. Словарное богатство русского языка 29.16 KB
  Обратить внимание на информационное и словесное богатство русского языка. Поговорить о важности расширения словарного запаса. Развивать навыки работы со словарями.
62398. СОЦИАЛЬНЫЕ ОТНОШЕНИЯ И ВЗАИМОДЕЙСТВИЯ 15.87 KB
  В повседневной жизни мы совершаем множество элементарных актов социального взаимодействия например: мы приветствуем человека в автобусе мы уступаем место пожилым и др. Взаимодействия осуществляются на разных уровнях: 1 межличностные; 2 личностно-групповые; 3 межгрупповые.
62399. Социология как наука 24.46 KB
  Отраслевые теории формируются на стыках социологии с другими науками например экономическая социология политическая социология социология права и др. Объектом политологии то что исследуется является политическая реальность и ее различные стороны и отношения.
62403. OUR CLASSES 26.38 KB
  Becoming a doctor is a dream for millions of students around the world. Studying medicine can lead to a career as a doctor of course, but it doesnt have to be this cut and dry.
62405. Работа с разными видами глин: лепка по гипсовой форме 583.76 KB
  Из поколения в поколение переходят традиции ремесла и искусства игрушки передаются в народе представления о жизни труде красоте. Древнейшие глиняные игрушки найденные археологами на территории нашей страны относятся к эпохе бронзы ко II тысячелетию до Рождества Христова.