Методы оценки температурного состояния

Рис.1.1 Пример прямоугольной равномерной сетки, построенной для прямоугольной области изменения переменных G (x,t).

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

Рис.1.2 Пример нерегулярной согласованной разностной сетки.

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

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

Разностная схема метода конечных элементов.

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

Рис.1.3 Равномерная конечноэлементная сетка, состоящая из треугольников, применяемая в методе конечных элементов.

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

Разностные схемы для нестационарных задач.

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

Для нестационарных задач в области вводится пространственная сетка, с которой связывается некоторое конечномерное пространство. Вводится сетка и по времени, для простоты, равномерная. Приближенное решение рассматривается как функция дискретного аргумента. Операторно-разностная схема связывает разностное решение на нескольких временных слоях. Такая разностная схема является многослойной.

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

Метод конечных разностей.

В качестве метода решения системы дифференциальных уравнений выбирается численный математический метод конечных разностей - широко известный и простейший метод интерполяции. Метод конечных разностей означает по сути обратный переход от дифференциальной модели к интегральной. При осуществлении данной методики строится конечно-разностная сетка и записываются конечно-разностные аналоги дифференциальных уравнений теплопроводности (разностная схема). Для аппроксимации дифференциальных уравнений теплопроводности применяется неявная консервативная итерационная разностная схема, реализуемая градиентным методом покоординатного спуска (Гаусса-Зейделя), являющимся классическим итерационным методом решения системы линейных уравнений. Неявной она является потому, что содержит несколько неизвестных значений функции на новом слое. Это означает, что значение функции нельзя явно выразить через значения функции на данном слое. Такая схема является безусловно устойчивой. Неявность разностной схемы достигается применением итерационной процедуры на каждом временном слое. Решение в узлах сетки получается приближенным (разностным). Поскольку одна из переменных имеет физический смысл времени t, то сетка строится так, чтобы среди ее линий были линии t=tm, где m - номер индекса дискретного момента времени. То есть переменная t не непрерывна, а увеличивается на дискретное значение. Решение численной задачи получается в виде таблицы.