Лекция №4. Тема 1.3 Методы решения краевых задач

Тема 1.3 Методы решения краевых задач. Метод разделения переменных (Метод Фурье).

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

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

Метод применим для конечных областей.

Рассмотрим метод Фурье применительно к следующей задаче

(23)

( ).

В случае декартовых координат ;

;

- конечная пространственная область.

Предположим, что ГУ приведены к однородным. Тогда, частное решение уравнения (23) ищем в виде произведения двух функций

, (24)

одна из которых зависит только от времени, а другая - только от пространственных координат; А- произвольная постоянная.

Если (24) подставить в (23), то получим два дифференциальных уравнения относительно

; (25)

относительно

; (26)

где - постоянная разделения.

Решение (25) элементарно

Решение (26) получено лишь для некоторых частных случаев.

Задача нахождения тех значений постоянной , для которых существуют нетривиальные решения уравнения (26) (называется собственными функциями), удовлетворяющие граничным условиям, называется задачей Штурма-Лиувилля.

При постоянном коэффициенте задача Штурма-Лиувилля решена для тел, образованных пересечением координатных поверхностей в различных системах координат. Например, для одномерной задачи решение уравнения (26) имеет вид

в прямоугольных координатах

; (27)

в сферических координатах

; (28)

в цилиндрических координатах

, (29)

где С, D – произвольные постоянные; а числа определяются из граничных условий задачи; - функция Бесселя первого ряда нулевого порядка; - функция Бесселя второго рода нулевого порядка.

Определив выражения для функций и , решение уравнения (23) с соответствующими ГУ представится в виде

где - собственные функции, отвечающие собственным числам .

Определим коэффициенты из начального условия [при ]

где - рассматриваемая конечная область; N – норма собственной функции , равная

При этом используется ортогональность собственных функций , является свойством собственных функций задачи Штурма-Лиувилля.

Окончательно решение краевой задачи имеет вид

. (30)

При условии, что ряд (30) допускает почленное дифференцирование дважды по пространственным координатам и один раз по времени.

Нагрев неограниченной пластины.

Дана неограниченная пластина, толщина которой равна 2R. В начальный момент времени пластина помещается в среду с постоянной температурой . Между ограничивающими поверхностями пластины и окружающей средой происходит теплообмен по закону Ньютона. Найти распределение температуры по толщине пластины в любой момент времени.

Дифференциальное уравнение теплопроводности и его краевые условия имеют вид

, , -R<x<R, (31)

T(x,0)=f(x); (32)

; (33)

. (34)

Решение проведем методом разделения переменных. Предположим, что функция четная, то есть , поэтому . Тогда вместо граничного условия (34) можно записать

. (35)

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