Следовательно, решение смешанной задачи нужно искать в виде суммы

и функция w(x,t) , будет уже решением однородной смешанной задачи

 

которая решается в соответствии с простейшей схемой метода Фурье.

Как и при решении задачи 229 придем к ряду

где положительные нули цилиндрической функции J0(x). Что касается коэффициентов Аk, то из начального условия w/t=0=0 будет следовать

Второе начальное условие приведет к равенству

откуда в силу ортогональности собственных функций с весом r на[0,l] найдем

Чтобы вычислить J1 , нужно сделать замену переменной интегриро- вания и воспользоваться равенством (73), соответственно получим

Чтобы вычислить

заметим, что

,

удовлетворяют уравнениям

Умножив первое из этих уравнений на y2( r ), второе на y1( r), вычитая полученные результаты и интегрируя в пределах от 0 до l, будем иметь

Возвращаясь теперь к J2 , получим

Подставляя теперь найденные значения J1 и J2в выражение для Вk , найдем

В итоге придем к ответу

Заметим в заключение, что рассматривается так называемый нерезо- нансный случай, когда

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

Р е ш е н ие. После перевода текста на формулы придем к смешанной задаче

С помощью замены

u(r,t)=u0+v(r,t)

получим смешанную задачу для функции v(r,t):

решение которой уже ищем в виде произведения R(rT(t) так, что

v(r,t)=R(r)T(t).

Разделив переменные, придем к двум обыкновенным дифферен- циальным уравнениям:

С учетом граничных условий относительно радиальной функции R(r) придем к задаче Штурма Лиувилля

Такая задача уже встречалась ранее в задачах 229 и 230, ее собствен- ные значения и собственные функции будут

Собственные функции ортогональны между собой с весом р(r)=r. Дифференциальное решение для функции T(t) имеет вид

его можно решить как линейное с постоянным коэффициентом и мето- дом разделения переменных, тогда получим

Перемножая Tk(t) на собственные функции и суммируя, придем к ряду

Подставляя этот ряд в начальное условие, будем иметь

Итак решение исходной задачи будет выражаться рядом

234. Дан неограниченный цилиндр радиуса l , на поверхности которого поддерживается постоянная концентрация u0 вещества. Определить коли- чество вещества, продиффундировавшего внутрь цилиндра в момент вре- мени t, на единицу длины, если начальная концентрация u(r,t)/t=0=0.

Р е ш е н и е. Поскольку уравнение теплопроводность есть одновре -

менно и уравнение диффузии, то будем иметь смешанную задачу

Эта задача идентична встретившейся задаче 233 , поэтому ее решение сразу выписываем:

Теперь мы должны взять интеграл от концентрации, т. е.

235. Дана неограниченная цилиндрическая труба , и ее началь- ная температура u(r,0)=f(r). Внутренняя и наружная поверхности трубы поддерживаются при нулевой температуре. Найти распределение температуры в сечении u(r,t) при t>0.

Р е ш е н и е. Здесь возникает смешанная задача

Полагая u(r,t)=R(r)T(t), после разделения переменных придем к соот- ношениям

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

Собственные функции будут ортогональными на отрезке [r1 ,r2] с весом r.

Общее решение дифференциального уравнения будет иметь вид

Подставив его в граничные условия, придем к соотношениям

для определения С1, С2 и λ .Однородная система линейных уравнений имеет нетривиальные решения, если ее определитель

Корни уравнения

и будут собственными значениями, а собственные функции, как нетрудно проверить, будут иметь вид

как уже отмечалось, они образуют ортогональную систему с весом r на отрезке Найдем квадрат нормы

Через yn(r) обозначим функцию

Радиальная собственная функция Rn(r) и yn(r) суть соответственно решения уравнений

Умножим первое уравнение на yn(r),второе – на Rn(r), вычтем из пер- вого второе и проинтегрируем по отрезку [r1, r2.]. В итоге получим

откуда следует, что

Здесь использованы соотношения, вытекающие из определения соб-

ственных значений и формулы для определителя Вронского

Из дифференциального уравнения для временной функции

поэтому умножая Tn(t) на собственную функцию, найдем ряд

Подставляя этот ряд в начальное условие, будем иметь

Ответ запишется в форме

где

236. Найдите распределение температуры внутри бесконечного круглого цилиндра радиуса l , если начальная температура и на поверхности цилиндра поддерживается нуле- вая температура.

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

Здесь положительный корень уравнения

Решите следующие смешанные задачи

 



img src="images/image-1751-41.gif">