Помогите мне решить проблему теплопроводности/переноса излучения. Математика меня подвела

Моя проблема: Тонкостенная трубка (длина л , диаметр Д и толщина стенки т Д ) находится в вакууме. Держится за один конец (в Икс "=" 0 ) источником тепла при постоянной температуре Т ( 0 ) "=" Т 0 . Единственный способ, которым он может рассеивать тепло, — это излучение. Я предполагаю, что излучение происходит только с внешней поверхности трубки. Электропроводность трубки равна к в [ Вт / м К ] и коэффициент излучения ϵ . Что такое профиль равновесной температуры Т ( Икс ) в трубке? (численное приближение будет делать).

Моя попытка:

В устойчивом состоянии,

Вопрос я н "=" Вопрос о ты т

По закону теплопроводности Фурье тепло, поступающее через концевое сечение, равно

Вопрос я н "=" к г Т г Икс | Икс "=" 0 × π Д т

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

Вопрос о ты т "=" 0 л ϵ о Т 4 г Икс × π Д

Приравнивая эти два, проблема становится

к т ϵ о г Т г Икс | Икс "=" 0 "=" 0 л Т 4 г Икс ,       Т ( 0 ) "=" Т 0

Попытка решить это в Mathematica безнадежна. Я делаю что-то неправильно? Как найти локальную дифференциальную форму уравнения? Можно еще упростить?

Спасибо за вашу помощь.

Проблема в том, что вы вообще не построили тепловой баланс. См.: stealthskater.com/Personal/Thesis.pdf , стр. 3. Хотя для удилища его можно довольно легко приспособить для трубки.
Спасибо за ссылку, очень актуально. Но как же я не построил тепловой баланс? Как Вопрос я н "=" Вопрос о ты т не является условием устойчивого состояния?

Ответы (1)

Вам нужно сделать дифференциальный тепловой баланс на небольшом отрезке трубы между х и х + Δ Икс .

Нагрев при x = π Д т к ( Т Икс ) Икс

Нагрев при x + Δ Икс "=" + π Д т к ( Т Икс ) Икс + Δ Икс

Потери тепла за счет излучения = π Д Δ Икс ϵ о Т 4

Уравнение теплового баланса:

+ π Д т к ( Т Икс ) Икс + Δ Икс π Д т к ( Т Икс ) Икс "=" π Д Δ Икс ϵ о Т 4

Деление на Δ Икс и принимая предел как Δ Икс приближается к нулю, дает:

к т 2 Т Икс 2 "=" ϵ о Т 4

Что, по сути, просто дифференцирует исходное уравнение и более внимательно относится к знакам?
Точных форм решения этого уравнения, вероятно, не существует, но для физической интуиции это будет эквивалентно движению частицы под действием потенциала U ( Икс ) Икс 5 .
Спасибо @Farcher. Я отредактировал слова «тепло», чтобы читать «тепло в точке x +». Δ Икс . Я считаю, что это была единственная проблема со знаком. Я поддерживаю остальные знаки в уравнениях.
@ChesterMiller Я никоим образом не критиковал ваш вывод. Все, что я указывал, это то, что исходный вывод имеет знаковое отличие от вашего.
@MichaelSeifert моя интуиция кажется немного заржавевшей, не могли бы вы объяснить это немного подробнее (как траектория частицы в потенциале U (x) = -ax ^ 5 эквивалентна решению?) Спасибо
@uhoh: Если бы частица двигалась в таком потенциале, то она подчинялась бы уравнению движения м Икс ¨ "=" г U / г Икс "=" 5 а Икс 4 . Заменять Икс Т , т Икс , к т м , и 5 а ϵ о в этом уравнении, и вы получите одно и то же дифференциальное уравнение.
@Farcher Я не вижу никакого несоответствия между знаками в исходном выводе и в моем конечном результате.
@MichaelSeifert не является ли решением проблемы эволюция во времени распределения температуры: T (x, t), которое, как мы надеемся, будет иметь асимтотически постоянную равновесную форму, когда t стремится к бесконечности? Как здесь помогает одиночная частица, быстро ускоряющаяся влево x = f(t)?