Информатика -продвинутый курс




МОДЕЛИРОВАНИЕ ПРОЦЕССА ТЕПЛОПРОВОДНОСТИ - часть 6


Рис. 7.34. Графическая иллюстрация решения задачи из примера 3

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

Покроем отрезок [а, b] одномерной сеткой (т.е. разобьем на n равных частей, рис. 7.35) с узлами в точках

Искомую функцию и(х)

будем аппроксимировать ее значениями в узлах сетки. Конечно, такое представление не дает полного описания, но в промежуточных точках, если сетка достаточно «мелкая», возможна интерполяция.

Рис. 7.35. Одномерная сетка

Остановимся на разностной аппроксимации производных. Производная дает информацию о локальном изменении функции в пространстве и, соответственно, связывает ее значения в соседних узлах сетки. Очевидная аппроксимация первой производной в точке х,

имеет вид

(7.54)

Для крайних точек, однако, такая аппроксимация невозможна, и простейший способ - ограничиться односторонними разностями:

(7.55)

Разумеется, (7.54) и (7.55) дают простейшие аппроксимации. Втягивая большое количество узлов, можно получить аппроксимации более высокого порядка, но часто бывает достаточно описанных выше. Аналогичная им аппроксимация вторых производных имеет вид

(7.56)

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

Далее будем индексы, соответствующие временной сетке, писать надстрочно (вверху), а пространственной - подстрочно (внизу).


Содержание  Назад  Вперед