Теория нелинейной теплопроводности (стр. 1 из 5)

Министерство образования и науки Российской Федерации

Курсовая работа

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

Содержание

Аннотация

Введение

1. Теория нелинейной теплопроводности

2. Распространение тепловых возмущений в нелинейных средах

3. Пространственная локализация тепловых возмущений

4. Задача нелинейной теплопроводности с объемным поглощением

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

Заключение

Список используемой литературы

Аннотация

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

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

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

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

Введение

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

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

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

1. Теория нелинейной теплопроводности

Уравнение теплопроводности, учитывающее зависимость свойств среды от температуры и нелинейную зависимость от температуры мощности распределенных в объеме тепловых источников, является квазилинейным параболическим уравнением вида

(1.1)

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

Нелинейное граничное условие на поверхности тела может иметь вид

(1.2)

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

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

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

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

u1(P)=u2(P)=u*(1.3)

другое граничное условие:

(1.4)

где k1, k2 и и1, u2 - коэффициенты теплопроводности и температуры двух соприкасающихся фаз соответственно; q* - удельная массовая теплота фазового перехода; V - мгновенная скорость перемещения фронта фазового перехода в направлении нормали

поверхности∑ .

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

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

(1.5)

Здесь

= р(u) с(u) - теплоемкость единицы объема среды;

Q*=pq*;

импульсная функция Хевисайда, производная которой

есть дельта-функция

.

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

эф=

(u)+Q*

.

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

При таком описании фазового перехода уравнение теплопроводности в отсутствие объемных тепловых источников примет вид

[c(u)+q*

]p(u)

(1.6)

Здесь

Фронт фазового перехода в такой постановке задачи находится как изотермическая поверхность u = u* = const, положение которой в пространстве, а в общем случае и форма, изменяются с течением времени.

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

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

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

,-∞˂x˂+∞

Начальное условие: u|t=0=f(x)

Решение:

(интеграл Пуассона).

Распределение температуры в стержне, ограниченном с одной стороны

, 0˂x˂+∞

Начальное условие: u|t=0=f(x)

Краевое условие: u|t=0=φ(t)

Решение:

2. Распространение тепловых возмущений в нелинейных средах

В работах Г.И. Баренблатта, Я.Б. Зельдовича, С.П. Курдюмова, Л.К. Мартинсона, А.А. Самарского и других найдены точные аналитические решения некоторых задач нелинейной теплопроводности. Анализ свойств этих решений позволяет обнаружить ряд важных нелинейных эффектов при распространении тепловых возмущений в средах, коэффициент теплопроводности которых зависит от температуры.

Рассмотрим среду, коэффициент теплопроводности k которой изменяется в зависимости от температуры и по степенному закону

k=k0uб (2.1)

где

> 0 - параметр нелинейности среды. Плотность среды ρ и ее теплоемкость будем считать постоянными, не зависящими от температуры. Такую среду, в отличие от среды с постоянным коэффициентом теплопроводности (δ = 0), назовем нелинейной, так как процесс теплопроводности в такой среде в отсутствие объемных тепловых источников описывается нелинейным, точнее, квазилинейным параболическим уравнением

(2.2)

где

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

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

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

3. Пространственная локализация тепловых возмущений

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

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

(3.1)

Здесь

- коэффициент поглощения.

Поглощение энергии в объеме нелинейной среды приводит к уменьшению интегральной тепловой (внутренней) энергии среды. Поэтому при интегрировании (3.1) по пространственному переменному

в пределах от -∞до +∞ находим

(3.2)

где

Так как

, то, интегрируя уравнение (3.2), получаем

Для решения задачи (3.1) перейдем с помощью преобразования

(3.3)

к новой функции v(x,t) . Тогда уравнение для V принимает вид

Вводя новое независимое переменное (преобразованное время) по правилу

(3.4)

получаем для функции

задачу

(3.5)

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

(3.6)

(3.7)

Зависимости U(τ) и x0(τ) в (7.7) определены формулами в которых время t следует заменить на τ, понимая под τ = τ (t) преобразованное по закону

(3.8)

временное переменное. При этом существенно, что преобразование

отображает полубесконечный интервал [0, +∞) по переменному t в ограниченный отрезок [0, τm) по переменному τ .

Финитное решение (3.6) задачи (3.1) представляет собой фронтовое решение, описывающее распространение тепловой волны от мгновенного сосредоточенного источника с конечной скоростью перемещения фронтов x=±x0(

).

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

в любой момент времени t > 0 равна нулю вне области

, где