Тема 2.1. Численный метод - метод конечных разностей (МКР)
Изложение численного метода
Из численных методов наибольшее распространение для решения задач теплопроводности получил метод конечных разностей (или, иначе, метод сеток). Его основная идея – это перевести дифференциальное уравнение в частных производных (в нашем случае это уравнение Лапласа) в алгебраические уравнения с заменой дифференциалов на конечные разности.
Алгоритм применения метода сеток для приближенного решения двумерных краевых задач стационарной теплопроводности заключается в следующем:
1. В плоской области L, в которой разыскивается решение, строится сеточная область L’, состоящая из одинаковых ячеек и приближающая заданную область (рис.1.). На пересечении линий, составляющих сетку, образуются так называемые узлы сетки i, j.
2. Заданное дифференциальное уравнение в частных производных – уравнение Лапласа
, (1)
которое описывает процесс стационарной теплопроводности, заменяется в каждом из внутренних узлов сетки, так называемым, конечно-разностным уравнением (см. ниже). Конечно-разностное уравнение представляет собой алгебраическое уравнение, связывающее значения температур в соседних узлах сетки.
Таким образом, дифференциальное уравнение (1) заменяется целой системой большого числа алгебраических уравнений.
3. На основании известных граничных условий на контуре M области L (например, при решении задач с граничными условиями 1-го рода) устанавливаются значения искомого решения в граничных узлах области L’, т.е. на контуре M’.
|
|
|
4. Значения во внутренних узлах области L’ находятся путем решения указанной системы алгебраических уравнений. На этом и заканчивается, в принципе, построение метода сеток.
Для решения систем алгебраических уравнений имеется большое число различных методов. Далее мы будем использовать здесь метод итераций.
Остановимся теперь подробнее на описании метода сеток. Для получения сеточной области L’ на заданную область L наносится сетка, состоящая из 2-х систем взаимно перпендикулярных прямых, отстоящих друг от друга на расстоянии h. Это расстояние h называется шагом сетки, выбор которого определяется требованиями расчета температурного поля двумерной области.
Контур M’ области L’ выбирается так, чтобы он был как можно ближе к контуру М области L.
В каждый из граничных узлов контура M’ задаются граничные значения, равные известным значениям точки контура М, ближайшей к данному узлу.
Р и с.1. Сеточная область двумерной задачи стационарной теплопроводности
Пусть температура вдоль оси х меняется так, как показано на рис.2.
|
|
|
Р и с.2. Схема определения разностных производных
Производную от t по x в узле i, j можно приближенно заменить на конечно-разностные соотношения следующими способами:
1) 
- разностная производная «вперед» для узла i, j;
2) 
- центральная разностная производная;
3) 
- разностная производная «назад» для узла i, j.
Вторую производную, очевидно, можно приближенно получить в виде
и
. (2)
Аналогичным образом запишем вторую производную по оси y
. (3)
Тогда уравнение (1) в конечных разностях будет иметь вид
(4)
или
, (5)
т.е. ti,j в узле i, j является среднеарифметическим значением температур в соседних узлах.
Такое уравнение может быть записано для каждой узловой точки i, j внутри области L’.
Записав уравнение вида (4) для всех точек внутренней сеточной области L’, заменим дифференциальное уравнение (1) целой системой линейных алгебраических уравнений типа (4).
Если внутренних точек области L’ было N, то получится система из N линейных алгебраических уравнений с N неизвестными. Решение этой системы уравнений может быть проведено методом итераций.
Дата добавления: 2018-02-28; просмотров: 417; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!