Найти решение первой смешанной задачи методом Фурье для уравнения теплопроводности на отрезке: ut=16uxx; 0<x<8; t>0,

Найти решение первой смешанной задачи методом Фурье для уравнения теплопроводности на отрезке: ut=16uxx; 0<x<8; t>0,
(5.1)
ux,0=x24, 0≤x≤4,8-x, 4<x≤8.
(5.2)
u0,t=u8,t=0,
(5.3)
Решение:
Для решения начально-краевой задачи (5.1) — (5.3) применим метод Фурье разделения переменных. Будем искать решение задачи в виде
ux,t=Xx∙Tt.
Подставим предполагаемую форму решения в исходное уравнение (1)
Xx∙T’ (t)=16X»(x)∙T(t)
Разделим равенство на 16Xx∙T(t)
T’ (t)16T(t)=X»xXx=-λ=const , λ>0,
т.к. левая часть равенства зависит только от t, а правая – только от x.
В результате переменные разделяются, и получается два линейных обыкновенных дифференциальных уравнения
T'(t)+16λTt=0,
X»(x)+λXx=0.
Подставляя ux,t в виде Xx∙Tt в граничные условия (5.3), получим
X0⋅Tt=0, X8⋅Tt=0.
Поскольку равенства должны выполняться тождественно, то
X0=0, X8=0.
(5.4)
Таким образом, для функции X(x) получили задачу Штурма-Лиувилля
X»(x)+λXx=0X0=0, X8=0
Общее решение имеет вид
Xx=C1cosλx+C2 sinλx.
Неизвестные коэффициенты C1, C2 найдем из граничных условий (5.4)
X0=C1=0 X8=C2 sin8λ=0
Получили спектральное уравнение для нахождения собственных значений λ задачи Штурма-Лиувилля
sin8λ=0,
8λ=πn, n=1,2,…
Собственные значения задачи равны
λn=πn82, n=1,2,…
Им соответствуют собственные функции
Xnx=sinπnx8, n=1,2,…
Уравнение для функции Tt примет вид
Tn’t+16πn82Tnt=0.
Tn'(t)+πn22Tnt=0.
Общее решение этого уравнения имеет вид
Tnt=Cne-πn22t.
Решение ux,t исходной задачи записывается в виде
ux,t=n=1∞TntXnx=n=1∞Cne-πn22tsinπnx8.
Коэффициенты Cn этого ряда найдем из начального условия (5.2)
ux,0=n=1∞Cnsinπnx8=x24, 0≤x≤4,8-x, 4<x≤8.
Коэффициенты Cn представляют собой коэффициенты разложения функции в правой части в ряд Фурье по собственным функциям sinπnx8n=1∞
Cn =2804x24sinπnx8dx+488-xsinπnx8dx=
=14-8πn04x24dcosπnx8+488-xdcosπnx8=
=-2πnx24cosπnx804-04cosπnx82x4dx+8-xcosπnx848-48cosπnx8-1dx=
=-2πn4cosπn2-1204xcosπnx8dx-4cosπn2+48cosπnx8dx=
=-2πn8πn-1204x dsinπnx8+sinπnx848=
=-16π2n2-12xsinπnx804-04sinπnx8dx+sinπn-sinπn2=
=-16π2n2-124sinπn2+8πncosπnx804-sinπn2=
=-16π2n2-3sinπn2-4πncosπn2-1=
=-16π2n2-3sinπn2-4πncosπn2+4πn=
=16π3n33πnsinπn2+4cosπn2-4.
Таким образом, решение исходной начально-краевой задачи имеет вид
ux,t=n=1∞16π3n33πnsinπn2+4cosπn2-4e-πn22tsinπnx8
Ответ:
ux,t=n=1∞16π3n33πnsinπn2+4cosπn2-4e-πn22tsinπnx8

1 Звезда2 Звезды3 Звезды4 Звезды5 Звезд (Пока оценок нет)
Загрузка...