Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
otvety_Vovy_i_moi_nemnozhko.docx
Скачиваний:
15
Добавлен:
01.03.2025
Размер:
761 Кб
Скачать
☆

Вопрос 12) Метод разделения переем. (Фурье) для ур-я диффузии.

П окажем этот метод на конкретном при-мере. Пусть есть шар R, нагрет до Т0 . Пусть в момент t=0 этот шар в переме-шиваемую жидк-ть T1 < T0. Найти распред-ние темпер-ры внутри шара: T(M,t)=? Решение: В силу симметр.  что искомое темп/поле будет зависеть только от T(r,t). В силу интенсивного перемешив/я и выс. теплопроводно-сти шара можно считать, что t0-ра пов-ти шара T(R,t) = T1 . Уравнение для t0-ры : T = - (1/a)T/t = 0. Сферич система : 1/r2 (/r)(r2 T/r) - 1/a(T/t) = 0. Для реш-я задачи удобно ввести вспом ф-ю: (r,t)=T(r,t)-T1; { (1/r2)/r(r2 /r) - 1/a(/t)=0, { (RШ,t) = 0, { (r,0) = T0 – T1 , {  (r,t) <  . Если применить метод Ф. непоср-но к этому ур-ю, то реш. может быть выраж. через ф-и Бесселя. Можно упростить ур-е так, чтобы прийти к тригоном. ф-ям. Из 1/r (/r)(r2 /r) - 1/a(/t)(r) = 0, можно получить 2/r2(r) - 1/a(/t)(r) = 0. Введем новую ф-ю. (r,t)=r. Сформулир. краев. зад. для этой ф-и. { (2/r2) - 1/a(/t) = 0; { (R,t) = 0 : (r,0) = r(T0-T1); { (r,t) / r < .В соотв. с методом Фур. ищем реш в виде: (r,t)=R(r)(t). R``-mR = 0 : `-ma = 0 - m=const разделения. Решение имеет вид: (t)=C1emat. Поскольку реш не может неогр-но возраст  положить m = -2. В этом случае имеем: (t)=C1exp(-2at); R(r) = C2 sin r + C3 cos r. Заметим, что ф-ция R(r)/r <  откуда имеем: C3 = 0. Удовлетворяя краев. усл. : R | S = 0  sin R0=0  k = k/R0 , k = 1,2. Частное реш. для k можно записать в виде: k (r,t) = Dk exp(-2k at)sin k r. Общее реш. ищется наложением : (r,t)=k Dk exp(-k2at)sin kr. Удовлетворяя нач. усл. зад.: r (T0 - T1) = k Dk sin rk  находим коэф. Dk : Dk = (T0 - T1) 2R0 / [k(-1)k-1] Возвращаясь к темп/полю T(r,t) получ окончат решение.: { T(r,t) = T1 + 2R0(T0-T1)/(r) k[ (-1)k-1/ k]exp(-k2at) sin k r; { k = k / R0.

Вопрос 13) Примен. Преобр-я Лапласа для реш. Ур-я дифф

Используя преобр-е Лапл. по времени можно ур-е дифф-и привести к ур. содержащее только пространственые производные. При этом нач. краев. задача преобр в краев. задачу, методы реш. кот. р/ны раньше. Рассм ур-е дифф-и: { - 1/a(/t)=f. Предпол, что все ф-и в этой зад. { (M,0) = 0 (M) принадл. множ-ву оригиналов тогда{ (M,t)|Mt = 1(M,t). использ. преобр-е Л. по t к ур-ю и доп. усл, и use его св-ва получим краев. задачу для изобра-жений. { 0 - P/a (0) = F0 - 0(M)/a; {0(M,P)|MS = 01(M,P). Заметим, что нач. усл. вошло в прав. часть ур-я. Если реш. этой зад. построено, то оригин. можно восст. с пом-ю -ла:  (M,t) = 1/2i -i+i0(N,P)eptdP, Re p =. Этот -л можно вычислить с помощью вычетов.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]