Фундаментальное решение уравнения теплопроводности. Теплопроводность. математическое описание, частные задачи теплопроводности Неоднородное уравнение теплопроводности однородные краевые условия

с начальными условиями

и граничными условиями

Решение этой задачи будем искать в виде ряда Фурье по системе собственных функций (94)

т.е. в форме разложения

считая при этом t параметром.

Пусть функции f (x , t ) является непрерывной и имеет кусочно-непрерывную производную 1-го порядка по х и при всех t >0 выполняются условия

Предположим теперь, что функции f (x , t ) и
можно разложить в ряд Фурье по синусам

, (117)

(118)

, (119)

. (120)

Подставим (116) в уравнение (113) и с учетом (117), получим

Это равенство выполняется тогда, когда

, (121)

или, если
, то это уравнение (121) можно записать в виде

. (122)

Пользуясь начальным условием (114) с учетом (116), (117) и (119) получаем, что

. (123)

Таким образом, для нахождения искомой функции
приходим к задаче Коши (122), (123) для обыкновенного неоднородного дифференциального уравнения первого порядка. Пользуясь формулой Эйлера можно записать общее решение уравнения (122)

а с учетом (123) решение задачи Коши

Следовательно, когда мы подставим значение этой функции в выражение (116), в итоге получим решение исходной задачи

(124)

где функции f (x , t ) и
определены формулами (118) и (120).

Пример 14. Найти решение неоднородного уравнения параболического типа

при начальном условии

(14.2)

и граничных условиях

. (14.3)

▲ Подберем сначала такую функцию  , чтобы удовлетворяла граничным условиям (14.3). Пусть, например,  = xt 2 . Тогда

Следовательно, функция определяемая как

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

(14.5)

однородным граничным условиям

и нулевым начальным условиям

. (14.7)

Применяя метод Фурье для решения однородного уравнения

при условиях (14.6), (14.7), положим

Приходим к следующей задаче Штурма-Лиувилля:

,
.

Решая эту задачу, находим собственные значения

и соответствующие им собственные функции

. (14.8)

Решение задачи (14.5)-(14.7) ищем в виде ряда

, (14.9)

(14.10)

Подставив
из (14.9) в (14.5) получим

. (14.11)

Для нахождения функции T n (t ) разложим функцию (1-х ) в ряд Фурье по системе функций (14.8) на интервале (0,1):

. (14.12)

и из (14.11) и (14.12) получаем уравнение

, (14.13)

которое является обыкновенным неоднородным линейным дифференциальным уравнением первого порядка. Его общее решение найдем по формуле Эйлера

а с учетом условия (14.10), найдем решение задачи Коши

. (14.14)

Из (14.4), (14.9) и (14.14) находим решение исходной задачи (14.1)- (14.3)

Задания для самостоятельной работы

Решить начально-краевые задачи

3.4. Задача Коши для уравнения теплопроводности

В первую очередь рассмотрим задачу Коши для однородного уравнения теплопроводности.

удовлетворяющее

Начнем с того, что заменим переменные x и t на
и введем в рассмотрение функцию
. Тогда функции
будут удовлетворять уравнениям

где
- функция Грина, определяемая формулой

, (127)

и обладающая свойствами

; (130)

. (131)

Умножив первое уравнение на G * , а второе на и и затем сложив полученные результаты, получим равенство

. (132)

После интегрирования по частям равенства (132) по в пределах от -∞ до +∞ и пов пределах от 0 доt , получим

Если предполагать, что функция
и ее производнаяограничены при
, то в силу свойств (131) интеграл в правой части (133) равен нулю. Следовательно, можно записать

Заменив в этом равенстве на
, а
на
, получим соотношение

Отсюда, используя формулу (127) окончательно получим

. (135)

Формула (135) называется формулой Пуассона и определяет решение задачи Коши (125), (126) для однородного уравнения теплопроводности с неоднородным начальным условием.

Решение же задачи Коши для неоднородного уравнения теплопроводности

удовлетворяющее неоднородному начальному условию

представляет собой сумму решений:

где является решением задачи Коши для однородного уравнения теплопроводности. , удовлетворяющее неоднородному начальному условию, аявляется решением, удовлетворяющее однородному начальному условию. Таким образом, решение задачи Коши (136), (137) определяется формулой

Пример 15. Найти решение уравнения

(15.1)

для следующего распределения температуры стержня:

▲ Стержень является бесконечным, поэтому решение можно записать, используя формулу (135)

Так как
в интервале
равна постоянной температуре, а вне этого интервала температура равна нулю, то решение принимает вид

. (15.3)

Полагая в (15.3)
, получим

Поскольку

представляет собой интеграл вероятностей, то окончательное решение исходной задачи (13.1), (13.2) можно выразить формулой

.▲

Формулы для расчета температурного поля и теплового потока в частных задачах стационарной и нестационарной теплопроводности получают исходя из математического описания (математической модели) процесса. Основу модели составляет дифференциальное уравнение теплопроводности, которое выводится с привлечением первого закона термодинамики для тел, не совершающих работы, и закона теплопроводности Фурье. Дифференциальное уравнение физического процесса обычно выводится при тех или иных допущениях, упрощающих процесс. Поэтому получаемое уравнение описывает класс процессов только в пределах принятых допущений. Каждая конкретная задача описывается соответствующими условиями однозначности. Таким образом, математическое описание процесса теплопроводности включает дифференциальное уравнение теплопроводности и условия однозначности.

Рассмотрим вывод дифференциального уравнения теплопроводности при следующих допущениях:

а) тело однородно и анизотропно;
б) коэффициент теплопроводности зависит от температуры;
в) деформация рассматриваемого объема, связанная с изменением температуры, очень мала по сравнению с самим объемом;
г) внутри тела имеются равномерно распределенные внутренние источники теплоты q v = f(x, у, z, т) = const;
д) перемещение макрочастиц тела относительно друг друга (конвекция) отсутствует.

В теле с принятыми характеристиками выделяем элементарный объем в форме параллелепипеда с ребрами dx, dy, dz, определенно ориентированный в ортогональной системе координат (рис. 14.1). В соответствии с первым законом термодинамики для тел, не совершающих работы, изменение внутренней энергии dU вещества в выделенном объеме за время dx равно сумме теплоты, поступающей

Рис. 14.1.

в объем вследствие теплопроводности dQ x , и теплоты, выделенной внутренними источниками dQ 2 ".

Из термодинамики известно, что изменение внутренней энергии вещества в объеме dV за время dx равно

где dG = рdV - масса вещества; р - плотность; с - удельная массовая теплоемкость (для сжимаемых жидкостей c = c v (изохорной теплоемкости)).

Количество энергии, выделенное внутренними источниками,

где q v - объемная плотность внутренних источников теплоты, Вт/м 3 .

Тепловой поток, поступающий в объем теплопроводностью, разделим на три составляющих соответственно направлению осей координат: Через противоположные грани теплота будет

отводиться в количестве соответственно Разница между количеством подведенной и отведенной теплоты эквивалентна изменению внутренней энергии вследствие теплопроводности dQ v Представим эту величину как сумму составляющих по осям координат:

Тогда в направлении оси х имеем

Поскольку -

плотности тепловых потоков на поотивоположных гоанях.

Функция q x+dx является непрерывной в рассматриваемом интервале dx и может быть разложена в ряд Тейлора:

Ограничиваясь двумя первыми членами ряда и подставляя в (14.6), получаем

Аналогичным образом получаем:

После подстановки (14.8)-(14.10) в (14.4) имеем

Подставляя (14.2), (14.3) и (14.11) в (14.1), получаем дифференциальное уравнение переноса теплоты теплопроводностью с учетом внутренних источников:

Согласно закону теплопроводности Фурье записываем выражения для проекций на оси координат плотности теплового потока:

где Х х, Х у, X z - коэффициенты теплопроводности в направлении координатных осей (тело анизотропное).

Подставляя эти выражения в (14.12), получаем

Уравнение (14.13) называют дифференциальным уравнением теплопроводности для анизотропных тел с независимыми от температуры физическими свойствами.

Если принять X = const, а тело изотропным, уравнение теплопроводности принимает вид

Здесь а = Х/(ср), м 2 /с, - коэффициент температуропроводности,

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

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

где г, z, Ф - соответственно радиальная, осевая и угловая координаты.

Уравнения (14.13), (14.14) и (14.15) описывают процесс теплопроводности в самом общем виде. Конкретные задачи отличаются условиями однозначности , т.е. описанием особенностей протекания рассматриваемого процесса.

Условия однозначности. Исходя из физических представлений о теплопроводности можно выделить факторы, влияющие на процесс: физические свойства вещества; размеры и форма тела; начальное распределение температуры; условия теплообмена на поверхности (границе) тела. Таким образом, условия однозначности подразделяются на физические, геометрические, начальные и граничные (краевые).

Физическими условиями задаются физические параметры вещества X, с, р и распределение внутренних источников.

Геометрическими условиями задаются форма и линейные размеры тела, в котором протекает процесс.

Начальными условиями задается распределение температуры в теле в начальный момент времени t = /(х, у, z ) при т = 0. Начальные условия имеют значение при рассмотрении нестационарных процессов.

В зависимости от характера теплообмена на границе тела граничные (краевые) условия подразделяются на четыре рода.

Граничные условия первого рода. Задается распределение температуры на поверхности t n в течение процесса

В частном случае температура поверхности может оставаться постоянной (/ п = const).

Граничные условия первого рода имеют место, например, при контактном нагреве в процессах склеивания фанеры, прессования древесно-стружечных и древесно-волокнистых плит и т.п.

Граничные условия второго рода. Задается распределение значений плотности теплового потока на поверхности тела в течение процесса

В частном случае тепловой поток на поверхности может оставаться постоянным (

Граничные условия третьего рода соответствуют конвективному теплообмену на поверхности. При этих условиях должна задаваться температура жидкости, в которой находится тело, Г ж = /(т), и коэффициент теплоотдачи ос. В общем случае коэффициент теплоотдачи - переменная величина, поэтому должен задаваться закон его изменения а =/(т). Возможен частный случай: / ж = const; а = const.

Граничные условия четвертого рода характеризуют условия теплообмена тел с различными коэффициентами теплопроводности при их идеальном контакте, когда теплота передается теплопроводностью и тепловые потоки по разные стороны поверхности контакта равны:

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

Уравнения (14.14) и (14.15) решаются достаточно просто аналитически для одномерного стационарного теплового режима. Решения рассмотрены ниже. Для двумерных и трехмерных стационарных процессов применяются приближенные численные методы

Для решения уравнений (14.13)-(14.15) в условиях нестационарного теплового режима используется ряд методов, рассмотренных подробно в специальной литературе . Известны точные и приближенные аналитические методы, численные методы и др.

Численное решение уравнения теплопроводности осуществляется в основном методом конечных разностей . Выбор того или иного метода решения зависит от условий задачи. В результате решения аналитическими методами получают формулы, применимые для решения круга инженерных задач в соответствующих условиях. Численные методы дают возможность получить температурное поле t=f(x, у, z, т) в виде набора дискретных значений температуры в различных точках в фиксированные моменты времени для конкретной задачи. Поэтому использование аналитических методов предпочтительно, однако это не всегда возможно для многомерных задач и сложных граничных условий.

При построении математической модели распространения тепла в стержне сделаем следующие предположения:

1) стержень сделан из однородного проводящего материала с плотностью ρ ;

2) боковая поверхность стержня теплоизолирована, то есть тепло может распространяться только вдоль осиОХ ;

3) стержень тонкий - это значит, что температура во всех точках любого поперечного сечения стержня одна и та же.

Рассмотрим часть стержня на отрезке [х, х + ∆х ] (см. рис. 6) и воспользуемся законом сохранения количества тепла:

Общее количество тепла на отрезке [х, х + ∆х ] = полному количеству тепла, прошедшему через границы + полное количество тепла, образованного внутренними источниками.

Общее количество тепла, которое необходимо сообщить участку стержня, чтобы повысить его температуру на ∆U , вычисляется по формуле: ∆Q=CρS∆x∆U , где С -удельная теплоемкость материала (=количеству тепла, которое нужно сообщить 1 кг вещества, чтобы поднять его температуру на 1°), S - площадь поперечного сечения.

Количество тепла, прошедшее через левый конец участка стержня за время ∆t (тепловой поток) вычисляется по формуле: Q 1 = -kSU x (x, t)∆t , где k - коэффициент теплопроводности материала (= количеству тепла, протекающего в секунду через стержень единичной длины и единичной площади поперечного сечения при разности температур на противоположных концах, равной 1°). В этой формуле особого пояснения требует знак минус. Дело в том, что поток считается положительным, если он направлен в сторону увеличения х , а это, в свою очередь, означает, что слева от точки х температура больше, чем справа, то есть U x < 0 . Следовательно, чтобыQ 1 был положительным, в формуле стоит знак минус.

Аналогично, тепловой поток через правый конец участка стержня вычисляется по формуле: Q 2 = -kSU x (x +∆x,t)∆t .

Если предположить, что внутренних источников тепла в стержне нет, и воспользоваться законом сохранения тепла, то получим:

∆Q = Q 1 - Q 2 => CpS∆x∆U = kSU x (x + ∆х, t) ∆t - kSU x (x, t)∆t .

Если это равенство поделить на S∆x∆t и устремить ∆х и ∆t к нулю, то будем иметь:

Отсюда уравнение теплопроводности имеет вид

U t =a 2 U xx ,

где - коэффициент температуропроводности.

В случае, когда внутри стержня имеются источники тепла, непрерывно распределенные с плотностью q(x,t) , получится неоднородное уравнение теплопроводности

U t = a 2 U xx + f(x,t) ,
где .

Начальные условия и граничные условия.

Для уравнения теплопроводности задается только одно начальное условие U| t=0 = φ(х) (или в другой записиU(x,0) = φ(х) ) и физически оно означает, что начальное распределение температуры стержня имеет вид φ(х) . Для уравнений теплопроводности на плоскости или в пространстве начальное условие имеет такой же вид, только функция φ будет зависеть, соответственно, от двух или трех переменных.

Граничные условия в случае уравнения теплопроводности имеют такой же вид, как и для волнового уравнения, но физический смысл их уже иной. Условия первого рода (5) означают, что на концах стержня задана температура. Если она не изменяется со временем, то g 1 (t) ≡ Т 1 и g 2 (t) ≡ Т 2 , где Т 1 и Т 2 - постоянные. Если концы поддерживаются все время при нулевой температуре, то Т 1 = Т 2 = 0 и условия будут однородными. Граничные условия второго рода (6) определяют тепловой поток на концах стержня. В частности, если g 1 (t) = g 2 (t) = 0 , то условия становятся однородными. Физически они означают, что через концы не происходит теплообмен с внешней средой (эти условия еще называют условиями теплоизоляции концов). Наконец, граничные условиятретьего рода (7) соответствуют случаю, когда через концы стержня происходит теплообмен с окружающей средой по закону Ньютона (напомним, что при выводе уравнения теплопроводности мы считали боковую поверхность теплоизолированной). Правда, в случае уравнения теплопроводности условия (7) записываются немного по-другому:

Физический закон теплообмена со средой (закон Ньютона) состоит в том, что поток тепла через единицу поверхности в единицу времени пропорционален разности температур тела и окружающей среды. Таким образом, для левого конца стержня он равен Здесь h 1 > 0 - коэффициент теплообмена с окружающей средой, g 1 (t) - температура окружающей среды на левом конце. Знак минус поставлен в формуле по той же причине, что и при выводе уравнения теплопроводности. С другой стороны, в силу теплопроводности материала поток тепла через этот же конец равен Применив закон сохранения количества тепла, получим:

Аналогично получается условие (14) на правом конце стержня, только постоянная λ 2 может быть другой, так как, вообще говоря, среды, окружающие левый и правый конец, бывают разные.

Граничные условия (14) являются более общими по сравнению с условиями первого и второго рода. Если предположить, что через какой-либо конец не происходит теплообмена со средой (то есть коэффициент теплообмена равен нулю), то получится условие второго рода. В другом случае предположим, что коэффициент теплообмена, например h 1 , очень большой.

Перепишем условие (14) при х = 0 в виде и устремим . В результате будем иметь условие первого рода:

Аналогично формулируются граничные условия и для большего числа переменных. Для задачи о распространении тепла в плоской пластине условие означает, что температура на ее краях поддерживается нулевой. Точно так же, условия и внешне очень похожи, но в первом случае оно означает, что рассматривается плоская пластина и края ее теплоизолированы, а во втором случае оно означает, что рассматривается задача о распространении тепла в теле и поверхность его теплоизолирована.

Решение первой начально-краевой задачи для уравнения теплопроводности.

Рассмотрим однородную первую начально-краевую задачу для уравнения теплопроводности:

Найти решение уравнения

U t = U xx , 00,

удолетворяющее граничным условиям

U(0,t) = U(l,t)=0, t>0 ,

и начальному условию

Решим эту задачу методом Фурье.

Шаг 1 . Будем искать решения уравнения (15) в виде U(x,t) = X(x)T(t) .

Найдем частные производные:

Подставим эти производные в уравнение и разделим переменные:

По основной лемме получим

Отсюда следует

Теперь можно решить каждое из этих обыкновенных дифференциальных уравнений. Обратим внимание на то, что используя граничные условия (16), можно искать не общее решение уравнения б), а частные решения, удолетворяющие соответствующим граничным условиям:

Шаг 2. Решим задачу Штурма-Лиувилля

Эта задача совпадает с задачей Штурма-Лиувилля, рассмотренной в лекции 3. Напомним, что собственные значения и собственные функции этой задачи существуют только при λ>0.

Собственные значения равны

Собственные функции равны (См. решение задачи)

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

Мгновенный точечный источник

Для бесконечного тела, в начале координат которого действует мгновенный точечный источник, решение дифференциального уравнения теплопроводности следующее:

где T - температура точки с координатами x,y,z; Q - количество тепла, выделившееся в момент t = 0 в начале координат; t - время, прошедшее с момента введения тепла; R - расстояние от начала координат, где действует источник, до рассматриваемой точки (радиус - вектор). У равнение (4) является фундаментальным решением уравнения теплопроводности при действии мгновенного точечного источника в бесконечном теле.

В любой момент t ? 0 температура самого источника (R = 0) отлична от нуля и с течением времени уменьшается по закону t -3/2 , оставаясь выше температур других точек тела. Вместе с удалением от источника температура понижается по закону нормального распределения exp(-R 2 /4at). Изотермическими поверхностями являются сферы с центром в источнике, и температурное поле в данный момент времени зависит лишь от радиуса. В начальный момент времени (t = 0) температура не определена (T = ?), что связано со схемой сосредоточенного источника, в котором в бесконечно малом объеме в начальный момент времени содержится конечное количество тепла Q.

На основе решения для бесконечного тела (4) можно вывести уравнение температурного поля для схемы полубесконечного тела, которая применяется для описания тепловых процессов в массивных изделиях. Пусть в полубесконечном теле, ограниченном поверхность S - S действует мгновенный точечный источник Д (рис. 4). Для массивных тел тепловые потоки внутри значительно больше потока теплоотдачи с поверхности. Поэтому поверхность полубесконечного тела можно считать адиабатической границей, для которой (см.п. 1.4)

Дополним полубесконечную область z > 0 до бесконечной, дбавив область z < 0. В образовавшемся объеме введем дополнительный (фиктивный) источник нагрева Ф(-z), идентичный действительному источнику Д(z), но расположенный симметрично по другую сторону границы S. На рис. 4 приведено распределение температур в бесконечном теле отдельно для действительного (T Д) и фиктивного (T ф) источников. Суммарная температура от обоих источников T = T Д + T ф. При этом на границе, что соответствует определению адиабатической границы (5). Если действительный источник находится на поверхности полубесконечного тела, то фиктивный с ним совпадает, и T=2T Д. Тогда температурное поле мгновенного точечного источника на поверхности полубесконечного тела

По такой же схеме моделируется и изотермическая граница (граничное условие 1-го рода) T S =0, но в этом случае T = T Д - T Ф. Следует подчеркнуть, что источник нагрева не может действовать на изотермической поверхности.

Графическое изображение температурного поля (6) требует четкого понимания пространственного положения поверхности, на которой строится распределение температуры. В декартовой системе координат (x, y, z) контрольными сечениями полубесконечного тела при действии точечного источника являются плоскости xy, xz и yz (рис. 5, а). Для полубесконечного тела изотермические поверхности являются полусферами (температура зависит от радиуса - вектора R). В плоскости xy изотермы, как сечение поверхности плоскостью

z=const, являются окружностями, а в других плоскостях - полуокружностями (рис. 5, б). Температурное поле мгновенного точечного источника в разные моменты времени представлено на рис. (6) (см. П 1.1.). На рисунке температура графически ограничена значением T=1000K|.

Температура в любой точке вне источника сначала возрастает, а затем убывает (рис.1.3). Момент достижения максимального значения температуры в данной точке найдется из условия

Дифференцируя выражение (6) по времени, получаем формулу для определения времени, когда температура максимальна

Максимальные темперы точек полубесконечного тела при действии точечного источника уменьшаются с расстоянием как R 3 .

АНАЛИТИЧЕСКИЕ МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

В настоящее время аналитическим путем решено очень большое количество одномерных задач теплопроводности.

А.В.Лыков, например, рассматривает четыре метода решения уравнения теплопроводности в условиях одномерной задачи: метод разделения переменных, метод источников, операционный метод, метод конечных интегральных преобразований.

В дальнейшем остановимся только на первом методе, получившем наибольшее распространение.

Метод разделения переменных при решении уравнения теплопроводности

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

T/?ф = a ? 2 t/?x 2 .(3.1)

Это уравнение является частным случаем однородного дифференциального уравнения с постоянными коэффициентами для некоторой функции t от двух переменных x и ф:

Легко проверить, что частным решением этого уравнения будет выражение

t = C exp (бx + вф).(3.3)

Действительно:

?t/?x = бС ехр (бx + вф);?t/?ф = вС ехр (бx + вф);
? 2 t/?x 2 = б 2 С ехр (бx + вф);
? 2 t/?ф 2 = в 2 С ехр (бx + вф);? 2 t/(?x ?ф) = бвС ехр (бx + вф).(3.4)

Совместное решение последних семи уравнении дает

a 1 б 2 + b 1 бв + c 1 в 2 + d 1 б + l 1 в + f 1 = 0.(3.5)

Последнее уравнение называется уравнением коэффициентов.

Переходя к уравнению (3.1) сопоставляя его с уравнением (3.2), заключаем, что

b 1 = c 1 = d 1 = f 1 = 0;a 1 = - a;l 1 = 1.(3.6)

Уравнение коэффициентов (3.5) для частного случая уравнения (3.1) приобретает вид

Б 2 a + в = 0(3.7)

в = б 2 a.(3.8)

Таким образом, частное решение (3.3) является интегралом дифференциального уравнения (3.1) и с учетом (3.8) приобретет вид

t = C exp (б 2 aф + бx).(3.9)

В этом уравнении можно задавать любые значения чисел для C, б, a.

Выражение (3.9) может быть представлено в виде произведения

t = C exp (б 2 aф) exp (бx),(3.10)

где сомножитель exp (б 2 aф) является функцией только времени ф, а сомножитель exp (бx) -- только расстояния x:

exp (б 2 aф) = f (ф);exp (бx) = ц (x).(3.11)

С увеличением времени ф температура во всех точках непрерывно растет и может стать выше наперед заданной, что в практических задачах не встречается. Поэтому обычно берут только такие значения б, при которых б 2 отрицательно, что возможно при б чисто мнимой величине. Примем

б = ± iq,(3.12)

где q -- произвольное действительное число (ранее значком q обозначали удельный тепловой поток),

В этом случае уравнение (3.10) приобретет следующий вид:

t = C exp (- q 2 aф) exp (± iqx).(3.13)

Обращаясь к известной формуле Эйлера

exp (± ix) = cos x ± i sin x(3.14)

и, пользуясь ею, преобразуем уравнение (3.13). Получим два решения в комплексном виде:

Суммируем левые и правые части уравнений (3.15), затем отделим действительные от мнимых частей в левой и правой частях суммы и приравняем их соответственно. Тогда получим два решения:

Введем обозначения:

(C 1 + C 2)/2 = D;(C 1 - C 2)/2 = C(3.17)

тогда получим два решения, удовлетворяющих дифференциальному уравнению теплопроводности (3.1):

t 1 = D exp (- q 2 aф) cos (qx);t 2 = C exp (- q 2 aф) sin (qx).(3.18)

Известно, что если искомая функция имеет два частных решения, то и сумма этих частных решений будет удовлетворять исходному дифференциальному уравнению (3.1), т. е. решением этого уравнения будет

t = C exp (- q 2 aф) sin (qx) + D exp (- q 2 aф) cos (qx),(3.19)

а общее решение, удовлетворяющее этому уравнению, можно записать в следующем виде:

Любые значения q m , q n , C i , D i в уравнении (3.20) будут удовлетворять уравнению (3.1). Конкретизация в выборе этих значений будет определяться начальными и граничными условиями каждой частной практической задачи, причем значения q m и q n определяются из граничных условий, а C i , и D i , -- из начальных.

Помимо общего решения уравнения теплопроводности (3.20) в котором имеет место произведение двух функций, одна из которых зависит от x, а другая - от ф, существуют еще решения, в которых такое разделение невозможно, например:

Оба решения удовлетворяют уравнению теплопроводности, в чем легко убедиться, продифференцировав их сначала по ф, а затем 2 раза по x и подставив результат в дифференциальное уравнение (3.1).

Частный пример нестационарного температурного поля в стенке

Рассмотрим пример применения полученного выше решения.

Исходные данные.

1. Дана бетонная стенка толщиной 2X = 0,80 м.
2. Температура окружающей стенку среды и = 0°С.
3. В начальный момент времени температура стенки во всех точках F(x)=1°C.
4. Коэффициент теплоотдачи стенки б=12,6Вт/(м 2 ·°С); коэффициент теплопроводности стенки л=0,7Вт/(м·°С); плотность материала стенки с=2000кг/м 3 ; удельная теплоемкость c=1,13·10 3 Дж/(кг·°С); коэффициент температуропроводности a=1,1·10 -3 м 2 /ч; относительный коэффициент теплоотдачи б/л = h=18,0 1/м. Требуется определить распределение температуры в стенке через 5 ч после начального момента времени.

Решение. Обращаясь к общему решению (3.20) и имея в виду, что начальное и последующие распределения температуры симметричны относительно оси стенки, заключаем, что ряд синусов в этом общем решении отпадает, и при x = Х оно будет иметь вид

Значения определены из граничных условий (без дополнительных здесь пояснений) и приведены в табл.3.1.

Располагая значениями из табл.3.1, находим искомый ряд значений по формуле

Таблица 3.1 Значения функций, входящих в формулу (3.24)


	0,982 0,189	--0,862 --0,507	0,713 0,701	10,03 --0,572 --0,820	13,08 0,488 0,874

т. е. Д1 = 1,250; Д2 = -- 0,373; Д3 = 0,188; Д4 = -- 0,109; Д5 = 0,072.

Начальное распределение температуры в рассматриваемой стенке приобретет следующий вид:

Чтобы получить расчетное распределение температуры через 5 ч после начального момента, необходимо определить ряд значений на время через 5 ч. Эти расчеты выполнены в табл.3.2.

Таблица 3.2 Значения функций, входящих в формулу (3.23)


A=(q ni X) 2 (aф/X 2)

Окончательное выражение для распределения температуры в толще стенки через 5 ч после начального момента

На рис.3.1 показано распределение температуры в толще стенки на начальный момент времени и через 5 ч. Наряду с общим решением здесь же изображены и частные, причем римскими цифрами указаны частные кривые, отвечающие последовательным слагаемым рядов (3.25) и (3.26).

Рис.3.1.

При решении практических задач обычно нет необходимости определять температуру во всех точках стенки. Можно ограничиться расчетом температуры лишь для какой-либо одной точки, например для точки в середине стенки. В этом случае объем вычислительных работ по формуле (3.23) значительно сократится.

Если начальная температура в рассмотренном выше случае равна не 1 °С, а Т с, то уравнение (3.20) примет вид

Решение уравнения теплопроводности при различных граничных условиях

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

Исходные данные. Стенка имеет толщину 2Х. В начальный момент во всех ее точках, кроме поверхности, температура Т с Температура на поверхности 0°С удерживается в течение всего расчетного периода.

Требуется найти t = f(x, ф).

Неподвижное водохранилище покрылось льдом при температуре наибольшей плотности воды (Т с = 4°С). Глубина водохранилища 5м (Х = 5 м). Рассчитать температуру воды в водохранилище через 3 месяца после ледостава. Температуропроводность неподвижной воды a = 4,8·10 -4 м 2 /ч. Тепловой поток у дна, т. е. при x = 0, отсутствует.

В течение расчетного периода (ф=3·30·24=2160ч) температура на поверхности удерживается постоянной и равной нулю, т. е. при x = Х Т п = 0°С. Весь расчет сводим в табл. 3 и 4. Эти таблицы позволяют вычислить значения температуры через 3 месяца после начального момента для глубин у дна, а затем выше через 1 м, т. е. t 0(дно) = 4°С; t 1 = 4°С; t 2 = 3,85°С; t 3 = 3,30°С; t 4 = 2,96°С; t 5(пов) = 0°С.

Таблица 3.3

Таблица 3.4

Как видим, в абсолютно неподвижной воде температурные возмущения весьма медленно проникают вглубь. В природных условиях в водоемах под ледяным покровом всегда наблюдаются течения либо гравитационные (проточные), либо конвективные (разноплотностные), либо, наконец, вызванные поступлением грунтовых вод. Все многообразие указанных природных особенностей следует учитывать при практических расчетах, а рекомендации к этим расчетам можно найти в пособиях и в работах К.И.Россинского .

Тело ограничено с одной стороны (полуплоскость). В момент времени ф = 0 во всех точках температура тела равна Т с. Для всех моментов времени ф > 0 на поверхности тела поддерживается температура Т п = 0°С.

Требуется найти распределение температуры в толще тела и потерю теплоты через свободную поверхность как функцию времени: t = f (x, ф),

Решение. Температура в любой точке тела и в любой момент времени

где есть интеграл Гаусса. Его значения в зависимости от функции даны в табл.3.5.

Таблица 3.5

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

Количество теплоты, теряемой единицей поверхности тела в окружающую среду, определяется по закону Фурье. За весь расчетный период с начального момента до расчетного

В начальный момент времени температура почвы от поверхности до значительной глубины была постоянной и равной 6°С. В этот момент температура на поверхности почвы упала до 0°С.

Требуется определить температуру почвы на глубине 0,5 м через 48 ч при значении коэффициента температуропроводности почвы a = 0,001 м 2 /ч, а также оценить количество теплоты, теряемое поверхностью за это время.

По формуле (3.29) температура почвы на глубине 0,5 м через 48 ч t=6·0,87=5,2°С.

Общее же количество теплоты, потерянной единицей поверхности почвы, при коэффициенте теплопроводности л = 0,35 Вт/(м·°С), удельной теплоемкости c = 0,83·10 3 Дж/(кг·°С) и плотности с = 1500 кг/м 3 определим по формуле (3.30) Q=l,86·10 6 Дж/м 2 .

интегральный теплопроводность теплота тело

Рис.3.2

Вследствие некоторого внешнего воздействия температура поверхности тела, ограниченного с одной стороны (полуплоскость), претерпевает периодические колебания около нуля. Будем считать, что эти колебания гармонические, т. е. температура поверхности меняется по косинусоиде:

где -- продолжительность колебания (период), T 0 -- температура поверхности,

T 0 макс -- ее максимальное отклонение,.

Требуется определить температурное поле как функцию времени.

Амплитуда колебаний температуры меняется с x по следующему закону (рис.3.2):

Пример к задаче № 3. Изменение температуры на поверхности сухой песчаной почвы в течение года характеризуется косинусоидальным ходом. Средняя годовая температура при этом равна 6°С при максимальных отклонениях от средней летом и зимой, достигающих 24 °С.

Требуется определить температуру грунта на глубине 1 м в момент, когда температура на поверхности равна 30°С (условно 1/VII).

Выражение косинусоиды (3.31) применительно к данному случаю (температуре поверхности) при T 0 макс = 24 0 С примет вид

Т 0 = 24 cos (2рф/8760) + 6.

Ввиду того, что поверхность грунта имеет среднюю годовую температуру 6°С, а не нуль, как в уравнении (3.32), расчетное уравнение примет следующий вид:

Приняв для грунта коэффициент температуропроводности a = 0,001 м 2 /ч и имея в виду, что по условию задачи необходимо определить температуру на конец расчетного периода (через 8760 ч от начального момента), найдем

Расчетное выражение (3.34) приобретет следующий вид: t = 24e -0,6 ·0,825 + 6 = 16,9 °С.

На той же глубине 1м максимальная амплитуда годового колебания температуры, согласно выражению (3.33), составит

T 1 макс = 24e -0,6 = 13,2 °С,

а максимальная температура на глубине 1 м

t 1 макс = T x макс + 6 = 13,2 + 6 =19, 2 °С.

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