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

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

ВВЕДЕНИЕ

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

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

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

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

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

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

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

Уравнения параболического типа наиболее часто встречаются при изучении процессов теплопроводности и диффузии. К этим уравнениям приводятся также задачи о движении вязкой жидкости, например, нефти.

Обсудим процесс распространения тепла в неравномерно нагретом твердом теле. Если тело нагрето неравномерно, то в нем происходит передача тепла из мест с более высокой температурой в места с более низкой температурой. Процесс может быть описан функцией u = u (x, y, z, t) дающей температуру u в каждой точке M (x, y, z) тела и в любой момент времени t .

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

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

Рассмотрим однородный стержень длины l, теплоизолированный с боков (через поверхность не происходит теплообмена с окружающей средой) и достаточно тонкий, чтобы в любой момент времени температуру во всех точках поперечного сечения можно было считать одинаковой. Расположим ось Ox так, чтобы один конец стержня совпадал с точкой x = 0, а другой — с точкой x = l (рис. 1).

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

Чтобы найти функцию u=u(x, t), надо составить дифференциальное уравнение, которому она удовлетворяет.

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

1. Количество тепла DQ, которое необходимо сообщить однородному телу, чтобы повысить его температуру на Du , равно

где c — удельная теплоемкость, m — масса тела.

Для стержня имеем

где ρ — плотность материала стержня; S — площадь его поперечного сечения.

2. Перенос тепла в теле подчиняется эмпирическому закону Фурье количество тепла ΔQ, протекающее за время Δt через площадку ΔS в направлении нормали Одномерное нестационарное уравнение теплопроводности имеет видк этой площадке, равно

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

где k — коэффициент внутренней теплопроводности (зависит от точки и не зависит от направления, если тело изотропно).

Для стержня имеем

Одномерное нестационарное уравнение теплопроводности имеет вид, (2)

где коэффициент k будем считать постоянным в силу предположения о его однородности. Если стержень неоднороден, то k = k(x).

3. Если внутри тела есть источники тепла, то выделение тепла можно характеризовать плотностью тепловых источников, т.е. количеством выделяемого (или поглощаемого) тепла в единицу времени в единице объема.

Обозначим через F (x, t), плотность источников в точке x рассматриваемого стержня в момент t. Тогда в результате действия этих источников на участке (x, x +Δx) за промежуток Δt будет выделено количество тепла

Одномерное нестационарное уравнение теплопроводности имеет вид(3)

И, наконец, воспользуемся законом сохранения энергии.

Итак, приступим к выводу уравнения. Выделим элементарный участок стержня, заключенный между сечениями x = x1 и x = x2 (x2 — x1 = Δx), и составим уравнение теплового баланса на отрезке [x1, x2]. Так как боковая поверхность стержня теплоизолирована, то элемент стержня может получать тепло только через поперечные сечения. Согласно (2) количество тепла, прошедшее через сечение x = x1, равно

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

через сечение x = x2:

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

Найдем приток тепла Одномерное нестационарное уравнение теплопроводности имеет видв элемент стержня: Одномерное нестационарное уравнение теплопроводности имеет вид

(К разности частных производных применена теорема Лагранжа).

Кроме того, в результате действия внутренних источников тепла на этом участке в течение времени Δt выделится количество тепла согласно (3)

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

Все тепло за время Одномерное нестационарное уравнение теплопроводности имеет видΔt пойдет на изменение температуры выделенного элемента стержня на величину Δu .И поэтому сообщенное количество тепла ΔQ , с другой стороны, может быть найдено согласно формуле (1):

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

В силу закона сохранения энергии имеем равенство

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

Сокращая на общий множитель SΔxΔt , получим уравнение

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

Введя обозначения Одномерное нестационарное уравнение теплопроводности имеет вид, Одномерное нестационарное уравнение теплопроводности имеет видпридем к уравнению

Одномерное нестационарное уравнение теплопроводности имеет вид(4)

Это и есть искомое дифференциальное уравнение распространения тепла в однородном стержне. Уравнение (4) называют уравнением теплопроводности, в котором постоянную a² температуропроводности. Коэффициент a² называют коэффициентом имеет размерность м² /с.

Уравнение (4) является линейным неоднородным уравнением параболического типа.

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

Процесс распределения температуры u = u (x, y, z, t) в изотропном теле описывается уравнением

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

которое кратко записывается так:

Одномерное нестационарное уравнение теплопроводности имеет вид(6)

где Одномерное нестационарное уравнение теплопроводности имеет вид— оператор Лаплас

Видео:Лекция №1.1 Явная и неявная схемы для уравнения теплопроводностиСкачать

Лекция №1.1 Явная и неявная схемы для уравнения теплопроводности

Теплопроводность при нестационарном режиме

Процессы передачи теплоты, в которых температурное поле и поле теплового потока изменяются во времени, называются неста­ционарными.

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

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

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

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

Все перечисленные задачи сводятся к нахождению температуры рассматриваемого тела как функции времени и координат t = f(t, x, у, z).

Эту зависимость можно найти, если проинтегрировать дифференциальное уравнение теплопроводности при заданных краевых условиях.

Для некоторых конкретных задач теплопроводности дифферен­циальное уравнение может быть упрощено: в случае передачи теп­лоты в одном направлении задача становится одномерной; при распростра­нении теплоты в двух направлениях задача является двухмерной. Для тел ци­линдрической формы удобно перейти к цилиндрическим координа­там, а для тел шаровой формы – к сферическим.

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

В качестве примера рассмотрим охлаждение неограниченной пластины.

Охлаждение неограниченной пластины

Будем рассматривать задачу теплопроводности при постоянных значениях теплофизических характеристик тела (l, с,r) с граничны­ми условиями третьего рода, так как они наиболее часто встреча­ются на практике. Задача формулируется следующим образом. Плоская неограни­ченная пластина толщиной d, имеющая во всех точках одинаковую начальную температуру tнч, в момент времени t = 0 помещается в среду, температура которой tж 0. Математически задачу можно сформули­ровать следующим образом. Дифференциальное уравнение теплопроводности для одномерной задачи без внутренних источников тепло­ты

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

Одномерное нестационарное уравнение теплопроводности имеет видгде х может изменяться в пределах 0 £ х £ d/2: так как охлаждение пластины происходит симметрич­но, целесообразно поместить на­чало координат в середину пласти­ны и рассматривать процесс толь­ко в одной ее половине (см. рисунок). Краевые условия:

1) начальное условие при t = 0 и 0 £ х £ d/2 t = tнч;

2) граничные условия: а) при х = 0и t > 0 (дt/дx)0 = 0, т. к. при симметричном охлаждении в середине пластины в любой момент времени температура будет максимальной; б) при х = l и t > 0 l(дt/дx)c = a(tctж).

Последнее выражение записано на основании равенства тепловых потоков на поверхности пластины: подходящего к поверхности из внутренних областей тела путем теплопроводности и отводимого от поверхности в процессе теплоотдачи.

Решение задачи в общем виде можно представить как функцию независимых переменных х и t и параметров процесса а,l, a, l, tж, tнч:

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

Для этого произведем сначала замену искомой величи­ны t так называемой избыточной температурой J = t – tж.

Так как dJ = dt,то запись дифференциального уравнения и гра­ничных условий от такой замены не изменится:

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

Приведем уравнение и граничные условия к безразмерному виду. Для этого еще раз произведем замену переменных: вместо избыточной температуры введем безразмерную избыточ­ную температуру Q = J/Jнч.Вместо координа­ты х введем безразмерную координату Х = х/l.Такая замена равно­сильна тому, что в качестве масштаба для измерения температуры используется величина Jнч, а в качестве масштаба длины – вели­чина l. Для сохранения равенств исходные уравнения в соответствую­щих местах необходимо умножить на масштабы температуры и длины. Тогда дифференциальное урав­нение будет иметь вид:

Одномерное нестационарное уравнение теплопроводности имеет вид, или после сокращения и преобразования Одномерное нестационарное уравнение теплопроводности имеет вид.

В такой форме дифференциальное уравнение безразмерно: величина l 2 имеет размерность времени и потому комплекс аt/l 2 безразмерен. Этот комплекс обознача­ется символом Fo и называется критерием Фурье:

Критерий Фурье можно трактовать как безразмерное время.

Окончательно дифференциальное уравнение теплопроводности в безразмерной записи получается в следующем виде:

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

Начальное условие: при Fo = 0, Qнч = 1;

где Qс = Jс/Jнч – безразмерная температура поверхности стенки; Bi = al/l – критерий Био.

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

Теперь искомая функция будет иметь вид Q = f(Fo,Bi, X).

Применяя метод разделения переменных решение дифференциального уравнения будет иметь вид

Одномерное нестационарное уравнение теплопроводности имеет вид, (1)

где Одномерное нестационарное уравнение теплопроводности имеет вид– коэффициенты уравнения;

mп – корни характеристического уравнения m/Bi = ctgm.

Значения mп и Ап приводятся в справочниках.

Результирующее выражение температурной функции, в форме произведения функции времени exp(-m 2 Fo) на некоторую функцию от координаты справедливо не только для пластины, но и для других тел, в которых распространение теплоты происходит в одном направлении, как, например, в бесконечно длинном цилиндре или шаре. Различаются результирующие выражения видом функции координаты: вместо cos – для пластины, для цилиндра появляется функция Бесселя, а для шара – гиперболическая. Для классических тел полу­чены аналитические решения задач нестационарной теплопроводнос­ти.

В соответствии с формой результирующих уравнений (1) порядок решения задачи нестационарной теплопроводности для тела классической формы следующий:

1. На основании исходных данных вычисляют безразмерную координату Х и критерии Bi и Fo. Здесь ха­рактерный размер тела: для пластины при симметричном охлажде­нии l = d/2,при одностороннем охлаждении l = d;для бесконечно длинного цилиндра и шара l = R,где R – радиус.

2. По величине критерия Bi в специальных таблицах нахо­дят значения mn и Ап для нескольких значений п.В обычных инже­нерных расчетах достаточно учитывать два-четыре члена суммы в формуле (1).

3. По формуле (1) или аналогичной ей для тел другой формы вычисляют значение безразмерной температуры Q в данной точке в заданный момент времени. Из Q определяют искомую температуру t = f(t, x).

Анализ решения (1) позволяет выявить влияние величины числа Bi на нестационарную теплопроводность. Рассмотрим два предельных случая: Bi ® ¥ и Bi ® 0.

Первый предельный случай:Bi ® ¥ (практически Bi >100). Для тела конечных размеров (l – конкретная конечная величина) этот случай соответствует условию a/l ® ¥, т. е. большим значе­ниям коэффициента теплоотдачи a и сравнительно малым значениям коэффициента теплопроводности l.В этом случае сразу после начала процесса температура поверхности тела при­нимает и в дальнейшем сохраняет постоянное значение tc = tж = const. Следовательно, интенсивность процесса охлаждения (нагрева) определяется внут­ренним процессом теплопроводности в теле и зависит только от фи­зических свойств и размеров тела.

При этом общее решение (1) упрощается: из числа определяющих критериев выпадает критерий Bi. Так, для точек, расположенных в средней плоскости пластины (при Х = 0), уравнение для безразмерной температуры при Fo > 0,3 приобретает вид

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

Второй предельный случай:Bi ® 0 (практически при Bi 2 ; Biпл = al/l; Foпл = аt/l 2 .

Величины Qц иQпл могут быть найдены по графикам с учетом расположения рассматриваемой точки в безграничном теле. Так, для точки 1 (рис.)величина Qц находится по графику для цент­ральных точек неограниченного цилиндра, а величина Qпл – по гра­фику для средней плоскости пластины. Для точки 2величина Qц определяется по тому же графику, что и для точки 1, а Qпл – по графику для поверхностных точек пластины. Для точки 3обе вели­чины находятся по графикам для поверхностных точек цилиндра и пластины. Для точки 4величина Qц определяется по графику для поверхностных точек цилиндра, а величина Qпл – по графику для средней плоскости пластины. Перечисленные четыре точки являются характерными для огра­ниченного цилиндра. Температуры остальных точек ограниченного цилиндра по графикам не могут быть найдены, но для их определе­ния можно воспользоваться соответствующими формулами.

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

Регулярный режим охлаждения (нагревания) тел

При значении Fo > 0,3 в выражениях типа (1) достаточно ограничиться одним первым членом ряда. В этом случае для плас­тины

Одномерное нестационарное уравнение теплопроводности имеет вид. (2)

Режим охлаждения (или нагрева), определяемый формулой (2), называется регулярным. Этот результат обобщается и на более сложные задачи охлаждения (нагрева) тел любой геометрической формы при условии tж и a = const:

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

где m – темп охлаждения, [1/с].

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

Логарифмируя последнее уравнение, получаем:

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

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

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

Если экспериментально определить изменение избыточной темпера­туры Jво времени t и построить зависимость в полулогарифмических координатах, то темп охлаждения в стадии регу­лярного режима найдется как

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

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

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

где С – полная теплоемкость тела;

y = JF/JV – коэффициент неравномерности распределения температуры в теле;

JF, JV – средние по поверхности и по объему температуры тела.

Из уравнения следует, что темп охлаждения т, однородного тела при конечном значении коэффициента теплоотдачи a пропорционален коэффициенту теплоотдачи, поверхности тела и обратно пропорционален его теплоемкости (первая теорема Кондратьева ).

Коэффициент y зависит от числа Bi, учитывающего условия протекания процесса на поверхности тела. Рассмотрим два предельных случая:

а) Bi ® 0 (практически Bi 100). При этом условии задача становится внутренней, и процесс охлажде­ния определяется только размерами тела и его физическими свойствами. В силу большой интенсивности теплообмена температура на поверх­ности тела принимает постоянное значение, равное температуре окружаю­щей среды. Коэффициент неравномерности распределения тем­пературы y = 0.

При Bi ® ¥, или, что то же, a ® ¥, темп охлаждения т становится прямо пропорциональным коэффициенту температуропроводности тела а (вторая теорема Кондратьева):

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

Коэффициент пропорциональности К зависит от геометрической формы и размеров тела и определяется в зависимости от формы тела по выражениям:

для шара радиусом r Одномерное нестационарное уравнение теплопроводности имеет вид;

для параллелепипеда с длиной граней l1, l2, l3 Одномерное нестационарное уравнение теплопроводности имеет вид;

для цилиндра длиной l и радиусом r Одномерное нестационарное уравнение теплопроводности имеет вид.

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

Видео:Одномерное уравнение теплопроводности. Виды краевых задачСкачать

Одномерное уравнение теплопроводности. Виды краевых задач

Теплопроводность при нестационарном температурном поле

Решить задачу теплопроводности при нестационарном температурном поле – значить установить зависимость между температурой t, временем Одномерное нестационарное уравнение теплопроводности имеет види координатами тела x,y,z. Такая зависимость получается решением дифференциального уравнения теплопроводности при определенных условиях однозначности.

При отсутствии внутренних источников тепла дифференциальное уравнение теплопроводности имеет вид

Одномерное нестационарное уравнение теплопроводности имеет вид Одномерное нестационарное уравнение теплопроводности имеет вид. (54)

Уравнение (54) является линейным, однородным дифференциальным уравнением второго порядка в частных производных. Решения такого уравнения обладают свойством наложения аналогично решениям обыкновенного однородного дифференциального уравнения: если t1 и t2 — частные решения уравнения, то выражение Одномерное нестационарное уравнение теплопроводности имеет видявляется также его решением при произвольных значениях постоянных С1 и С2. Поскольку у постоянных С1, и С2 возможны различные значения, уравнение типа (54) может иметь бесконечно большое количество частных решений.

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

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

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

Одномерное нестационарное уравнение теплопроводности имеет вид, (55)

где С – произвольная постоянная.

Подставляя решение (55) в уравнение (54) получим

Одномерное нестационарное уравнение теплопроводности имеет вид. (56)

Уравнение (56) можно переписать так

Одномерное нестационарное уравнение теплопроводности имеет вид. (57)

Левая часть уравнения(57) может зависеть только от Одномерное нестационарное уравнение теплопроводности имеет видили быть постоянным числом; она не зависит от координат. Правая часть может зависеть от координат или быть постоянным числом; она не зависит от времени. Поскольку уравнение (57) справедливо при любых значениях времени и координат, то правая и левая части его равны постоянной величине, которую обозначим через D.

Таким образом, мы получим два дифференциальных уравнения для определения вида функций U(τ) и P(x,y,z):

Одномерное нестационарное уравнение теплопроводности имеет вид; Одномерное нестационарное уравнение теплопроводности имеет вид. (58)

Решением уравнения (58) является

Одномерное нестационарное уравнение теплопроводности имеет вид, (59)

где С – постоянная интегрирования.

Постоянная величина D выбирается из физических соображений. В большинстве случаев при нагревании или охлаждении тел по истечении длительного времени температура распределяется в теле определенным образом. Для тепловых процессов, стремящихся к тепловому равновесию, величина D не может быть положительной, потому что можно задать такой промежуток времени, при котором температура в теле будет стремиться к бесконечности, что физически невозможно. Величина D не может равняться нулю, так как при D=0 функция U(τ) в уравнении (59) имела бы постоянное значение, а температура тела не зависела бы от времени, как это следует из уравнения (55), что не реально.

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

Рассматривая случай, когда Одномерное нестационарное уравнение теплопроводности имеет видD

Одномерное нестационарное уравнение теплопроводности имеет вид(65)

Одномерное нестационарное уравнение теплопроводности имеет вид. (65а)

Таким образом, количество теплоты Q не зависит от времени. Оно равно произведению площади, ограниченной кривой G и осью абсцисс x, на объемную теплоемкость cp.

Функцию Одномерное нестационарное уравнение теплопроводности имеет видназывают фундаментальным решением уравнения теплопроводности, поскольку она удовлетворяет этому уравнению. В самом деле, для неограниченного тела при одномерном потоке теплоты уравнение (54) имеет вид

Одномерное нестационарное уравнение теплопроводности имеет вид. (66)

Если функция G является решением уравнения (66), его можно записать так

Одномерное нестационарное уравнение теплопроводности имеет вид. (67)

Пользуясь уравнением (64), найдем выражения для Одномерное нестационарное уравнение теплопроводности имеет види Одномерное нестационарное уравнение теплопроводности имеет вид:

Одномерное нестационарное уравнение теплопроводности имеет вид(68)

Одномерное нестационарное уравнение теплопроводности имеет вид. (69)

Сопоставление последних двух выражений показывает, что действительно справедливо уравнение (67).

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

Преобразование Лапласа функции Одномерное нестационарное уравнение теплопроводности имеет вид, обозначаемое символом Одномерное нестационарное уравнение теплопроводности имеет вид, называется операцией умножения Одномерное нестационарное уравнение теплопроводности имеет видна Одномерное нестационарное уравнение теплопроводности имеет видс последующим интегрированием по в интервале от 0 до Одномерное нестационарное уравнение теплопроводности имеет вид

Одномерное нестационарное уравнение теплопроводности имеет вид. (70)

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

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

Одномерное нестационарное уравнение теплопроводности имеет вид. (71)

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

Преобразования Лапласа первой и второй производных функций Одномерное нестационарное уравнение теплопроводности имеет видопределяются соотношениями:

Одномерное нестационарное уравнение теплопроводности имеет вид(72)

Одномерное нестационарное уравнение теплопроводности имеет вид(73)

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

Интеграл Лапласа (71) и соотношения (72) и (73) можно использовать для интегрирования дифференциальных уравнений.

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

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

Одномерное нестационарное уравнение теплопроводности имеет вид. (74)

Решение уравнения (74) может быть выполнено аналитический и графически.

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

Рассмотрим применение численного метода к расчету распределения температуры в плоской стенке. Разбивая стенку на элементарные объемы V=δ·δ·1=δ 2 (рис. 4а,б), где δ – сторона элементарного объема.

Количество теплоты, подводимое к узловой точке в соответствии с законом Фурье, равно Одномерное нестационарное уравнение теплопроводности имеет вид. При малой величине δ тепловой поток q можно выразить через конечные разности

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

Рис. 4. Разбиение и числовая сетка определения нестационарного температурного поля а – одномерное температурное поле; б – двухмерное температурное поле

Одномерное нестационарное уравнение теплопроводности имеет вид(75)

где Δt – разность температур между смежными узловыми точками

Общее количество теплоты за время Δτ равно

Одномерное нестационарное уравнение теплопроводности имеет вид Одномерное нестационарное уравнение теплопроводности имеет вид(76)

Изменение внутренней энергии в данной узловой точке за время Δτ согласно первому началу термодинамики определяется следующим образом

Одномерное нестационарное уравнение теплопроводности имеет вид(77)

где t – температура в рассматриваемой узловой точке в момент времени τ;

Одномерное нестационарное уравнение теплопроводности имеет вид– температура в той же точке в момент времени Одномерное нестационарное уравнение теплопроводности имеет вид; V – объем элементарного участка.

Уравнение теплового баланса в конечных разностях для узловой точки 1 (см. рис. 4а) можно записать в виде

Одномерное нестационарное уравнение теплопроводности имеет вид. (78)

С учетом (76) уравнение (78) принимает вид

Одномерное нестационарное уравнение теплопроводности имеет вид(79)

Разделим уравнение (79) на Одномерное нестационарное уравнение теплопроводности имеет види с учетом того, что Одномерное нестационарное уравнение теплопроводности имеет види Одномерное нестационарное уравнение теплопроводности имеет вид— критерий Фурье (безразмерное время) искомая температура в рассматриваемой точке 1 в последующий интервал времени Одномерное нестационарное уравнение теплопроводности имеет видбудет равна

Одномерное нестационарное уравнение теплопроводности имеет вид. (80)

В случае двухмерного температурного поля тело разбивается на элементарные объемы с размерами ячеек Одномерное нестационарное уравнение теплопроводности имеет вид; расчетная схема узловых точек показана на рис. 4б.

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

Одномерное нестационарное уравнение теплопроводности имеет вид Одномерное нестационарное уравнение теплопроводности имеет вид. (81)

Уравнения (80 и 81) являются основой численного метода расчета нестационарной теплопроводности одномерного и двухмерного тела.

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

🎦 Видео

Метод Фурье для неоднородного уравнения теплопроводностиСкачать

Метод Фурье для неоднородного уравнения теплопроводности

Решение нестационарного уравнения теплопроводности в двухмерной постановке в ExcelСкачать

Решение нестационарного уравнения теплопроводности в двухмерной постановке в Excel

Решение уравнения теплопроводности в одномерной постановке в ExcelСкачать

Решение уравнения теплопроводности в одномерной постановке в Excel

Уравнение в частных производных Уравнение теплопроводностиСкачать

Уравнение в частных производных  Уравнение теплопроводности

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

8.1 Решение уравнения теплопроводности на отрезке

Стационарное решение одномерного уравнения теплопроводности.Скачать

Стационарное решение одномерного уравнения теплопроводности.

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

Решение Пуассона одномерного уравнения теплопроводности

6-1. Уравнение теплопроводностиСкачать

6-1. Уравнение теплопроводности

Двумерное нестационарное уравнение теплопроводности в MatLab l 2D Heat transfer equation in MatLabСкачать

Двумерное нестационарное уравнение теплопроводности в MatLab l 2D Heat transfer equation in MatLab

Общее уравнение динамики. Задача 1Скачать

Общее уравнение динамики. Задача 1

Теплопроводность плоской стенкиСкачать

Теплопроводность плоской стенки

Нестационарное уравнение теплопроводности в матлабеl Time dependent heat transfer equation in MatLabСкачать

Нестационарное уравнение теплопроводности в матлабеl Time dependent heat transfer equation in MatLab

Интуитивное понимание формулы теплопроводности (часть 11) | Термодинамика | ФизикаСкачать

Интуитивное понимание формулы теплопроводности (часть 11) | Термодинамика | Физика

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

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

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

Линейное однородное дифференциальное уравнение 2-го порядка с постоянными коэффициентами.

Решение уравнения теплопроводности в одномерной постановке в Excel с применением неявной схемыСкачать

Решение уравнения теплопроводности в одномерной постановке в Excel с применением неявной схемы

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

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

12. Как остывает шар (решение уравнения теплопроводности)Скачать

12. Как остывает шар (решение уравнения теплопроводности)
Поделиться или сохранить к себе: