Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лекции УМФ.doc
Скачиваний:
35
Добавлен:
18.11.2019
Размер:
353.79 Кб
Скачать

Уравнения параболического типа.

Они описывают нестационарное распространение тепла и диффузии. U/t = k02U/x2, k0 – коэффициент теплопроводности, х(0,L), U(0,t) = T1, U(L,t) = T2, U|t=0=T0. k02U/x2 = 0. U(x) = ((T2 – T1)/L)*x + T1. При стационарном процессе потоки тепла, водящие в любое поперечное сечение и выходящие их него, равны между собой. Следовательно, поток должен быть постоянен в любой произвольной точке х. Это означает, что по закону Фурье это достижимо при линейном профиле температуры. W= –gradT. Для уравнения теплопроводности существует тип решений, который имеет вид бегущей волны в стационарных формах. U = f(x–wt) (1). CU/t = /x(kU/x) (2). Подставим (1) в (2), получим: –cwf=(kf) (3). Уравнение (3) это линейное дифференциальное уравнение с постоянными коэффициентами. f() = A + Bexp{–cw/k },  = x–wt.

Рис

U(x,t) = A+Bexp{(x–wt)*(–cw/k)}, А – это температура на бесконечности, т.е. там, куда тепло ещё не дошло. Прогрев вещества, по которому со скоростью w распространяется вправо тепловая волна. f(x–wt)|x–wt=0 = U0.

Общий подход к решению.

Рис.

S – граница , Г состоит из нижнего основания и части боковой поверхности.

Первая краевая задача.

Пусть  конечная область трёхмерного пространства XYZ. Q в пространстве XYZt цилиндр. Рассмотрим задачу: U/t = a2[2U/x2 + 2U/y2 + 2U/z2] (1), U|t=0 = f(x,y,z), (x,y,z) (2), U|s = (p,t), t[0,T] (3).

Теорема (о max (min)): Функция U(x,y,z,t), удовлетворяющая однородному уравнению (1) внутри цилиндра QT и непрерывна вплоть до её границы, принимает наибольшее (наименьшее) значение или при t = 0, или на боковой поверхности цилиндра.

Док–во: Пусть М – наибольшее значение функции U в цилиндре QT, а m – наибольшее значение функции U на границе Г. Предположим, что существует такое решение U(x,y,z,t), что M>m. Пусть эта функция принимает значение М в точке (x0, y0, z0, t0), где (x0, y0, z0); t0[0,T]. Рассмотрим вспомогательную функцию V(x,y,z,t) = U(x,y,z,t) – ((M–m)/6d2)*[(x–x0)2 + (y–y0)2 + (z–z0)2]; где d – диаметр области . На боковой поверхности цилиндра и на его нижнем основании V(x,y,z,t)  m + (M–m)/6 = M/6 + 5m/6 < M; V(x0, y0, z0, t0) = M. Следовательно функция V также как и U не принимает наибольшее значение ни на боковой поверхности цилиндра, ни на основании. В силу непрерывности V она должна принимать наибольшее значение в некоторой точке (x,y,z,t), (x,y,z), 0tT. Тогда в этой точке вторые производные 2V/x2, 2V/y2, 2V/z2 не положительны, а V/t  0 => V/t – a2[2V/x2 + 2V/y2 + 2V/z2]  0 (*). V/t = U/t, 2V/x2 + 2V/y2 + 2V/z2 = 2U/x2 + 2U/y2 + 2U/z2 + (M–m)/d2 => V/t – a2[2V/x2 + 2V/y2 + 2V/z2] = {U/t – a2[2U/x2 + 2U/y2 + 2U/z2]} – a2(M–m)/d2 0, т.к. величина, стоящая в фигурных скобках равна нулю. Получили противоречие (*).

Следствия:

  1. Решение первой краевой задачи (1) – (3) в цилиндре QT единственно.

Док–во: если бы имели два решения U1 и U2, то U1 – U2 = w, удовлетворяющая однородному уравнению, должна была обратиться в нуль при t = 0 или на поверхности S, тогда w|t=0 = 0 в силу доказанной теоремы. Значит, U1 = U2.

  1. Решение задачи (1) – (3) непрерывно зависит от правых частей краевых и начальных условий. U|t=0 = f, U|s = .

Док–во: Если разность f и  не превосходит по абсолютной величине некоторое , то и разность решений w = U1 – U2, w|t=0 = , w|s =  по абсолютной величине не превосходит .

Такие задачи относятся к классу задач, корректных по Адамару.