Автор Тема: Найти температуру пластины в космосе  (Прочитано 932 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Пластина металла толщиной в несколько миллиметров находится в космосе.
Её освещает солнце, лучи перпендикулярны поверхности.
Найти приблизительно равновесную температуру в среднем слое пластины.

\(\sigma~=~5.67036713 \cdot 10^{-8}~Вт/(м^2 \cdot K^4)~-\) постоянная Стефана-Больцмана.
\(I_0~=~1396~Вт/м^2\) - тепловой поток излучения солнца на орбите земли.

Рекомендуемое обозначение физических параметров.
\(A_s~-\) поглощательная способность;
\(I_0~-\) интенсивность излучения;
\(S~-\) площадь поверхности.
\(q~-\) тепловой поток.
\(\lambda~-\) теплопроводность.
\(T~-\) абсолютная температура.
\(h~-\) толщина стенки.



« Последнее редактирование: 03 Июнь 2021, 14:40:33 от Ost »

Большой Форум

Загрузка...

Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #1 : 03 Июнь 2021, 18:47:22 »
А какая цель задачи - найти точное значение температуры? Потому что если металл черный и хорошо проводящий, то будет около 330K, а для плохо проводящего - около 400K. Если же хочется найти точно, то нужно учитывать теплопроводность, а значение теплопроводности в условии не дано.

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #2 : 03 Июнь 2021, 19:53:53 »
А какая цель задачи - найти точное значение температуры? Потому что если металл черный и хорошо проводящий, то будет около 330K, а для плохо проводящего - около 400K. Если же хочется найти точно, то нужно учитывать теплопроводность, а значение теплопроводности в условии не дано.
Надо показать, что если материал не тепловой изолятор, т.е. по теплопроводности близок по порядку к металлу и даже несколько хуже,
то зависимость температуры в среднем слое от теплопроводности, отражающей способности, толщины стенки, весьма незначительна и
Вы правильно написали 330K. У меня незначительно больше. 

Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #3 : 03 Июнь 2021, 20:31:45 »
Надо показать, что если материал не тепловой изолятор, т.е. по теплопроводности близок по порядку к металлу и даже несколько хуже,
то зависимость температуры в среднем слое от теплопроводности, отражающей способности, толщины стенки, весьма незначительна и
Вы правильно написали 330K. У меня незначительно больше.

Ну да, 333.1K.
Если взять типичный металл с теплопроводностью 50-500 Вт/м/К, то, т.к. распределение температуры по толщине линейно, при толщине 1 см теплопередача в пластине в 1 кв.м. даже при разности температур в 1К будет 5-50кВт. Т.е. для железа разность температур в задаче будет около 0,2 градуса. При всей нелинейности формулы Стефана излучение с холодной стороны будет практически таким же, формулы можно и не выписывать.   

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #4 : 03 Июнь 2021, 21:25:13 »
Ну да, 333.1K.
Если взять типичный металл с теплопроводностью 50-500 Вт/м/К, то, т.к. распределение температуры по толщине линейно, при толщине 1 см теплопередача в пластине в 1 кв.м. даже при разности температур в 1К будет 5-50кВт. Т.е. для железа разность температур в задаче будет около 0,2 градуса. При всей нелинейности формулы Стефана излучение с холодной стороны будет практически таким же, формулы можно и не выписывать.   
Да температура 333.1K.
Однако без выписывания формул не многие смогут написать результирующую формулу температуры для этого случая.


Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #5 : 03 Июнь 2021, 23:05:03 »
Однако без выписывания формул не многие смогут написать результирующую формулу температуры для этого случая.

Это не сложно, если не в латехе:

(1)           \(I_0=\sigma(T_1^4+T_2^4)\)
(2)           \(A_s\sigma~T_2^4=\lambda(T_1-T_2)/h\)

если взять \(T_1=T_0+dT/2, T_2=T_0-dT/2 (T_0\)-температура в середине пластины), то можно всё найти:

dT~\(\sigma~T^4h/\lambda\)

(если выкинуть все малые члены порядка dT/T0 и выше и считать T1=T2=T0). При желании можно учесть поправки
« Последнее редактирование: 11 Июнь 2021, 17:23:52 от Andrey_R »

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #6 : 04 Июнь 2021, 15:43:56 »
Это не сложно, если не в латехе:

I0=sigma*(T1^4+T2^4)
sigma*T2^4=(lambda/h)*(T1-T2)

если взять T1=T0+dT/2, T2=T0-dT/2 (T0-температура в середине пластины), то можно всё найти:

dT~sigma*T^4*(h/lambda)

(если выкинуть все малые члены порядка dT/T и выше и считать T1=T2=T). При желании можно учесть поправки
Да это так, но для доказательства, что температура \(T_0\) практически не зависит от параметров \(h\), \(\lambda\) и \(A_s\)
мы должны найти такую предельную функцию, которая в диапазоне допустимой вариации этих параметров будет иметь
такой вид \(T_0=T(I_0,\sigma) = 333.1~K\).


Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #7 : 04 Июнь 2021, 18:00:35 »
Да это так, но для доказательства, что температура \(T_0\) практически не зависит от параметров \(h\), \(\lambda\) и \(A_s\)
мы должны найти такую предельную функцию, которая в диапазоне допустимой вариации этих параметров будет иметь
такой вид \(T_0=T(I_0,\sigma) = 333.1~K\).

Ну, если обозначить, например,

 \(T_c=A_sI_0h/2\lambda \), то это и есть приближённое значение dT и при значениях h от 1 до 10 мм и lambda от 50 до 500 Вт/(м*К) будет в диапазоне 0,0015-0,15 K.

 \(T_a=(I_0/2\sigma)^\frac{1}{4}=333,1 K \) - приближённое значение Т0.

если, кроме того, обозначить  \(\alpha=\frac {T_c}{T_a}= 0,000005..0,0005\) - малый коэффициент при условиях задачи,  то из уравнений, которые были выше, получим прилижённые значения и поправки

(3)          \(\frac{dT}{T_0}=\alpha (1-2\alpha+O(\alpha^2))=\alpha+O(\alpha^2) \), a

(4)            \(T_0=T_a(1+\frac{3}{2}(\frac{dT}{T_0})^2+\frac{1}{16}(\frac{dT}{T_0})^4)^{-\frac{1}{4}}=T_a(1-\frac{3}{8}\alpha^2+O(\alpha^3)) \)
« Последнее редактирование: 11 Июнь 2021, 22:35:32 от Andrey_R »

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #8 : 05 Июнь 2021, 17:14:01 »
Ну, если обозначить, например,

 \(T_c=I_0h/2\lambda \), то это и есть приближённое значение dT и при значениях h от 1 до 10 мм и lambda от 50 до 500 Вт/(м*К) будет в диапазоне 0,0015-0,15 K.

 \(T_a=(I_0/2\sigma)^\frac{1}{4}=333,1 K \) - приближённое значение Т0.

если, кроме того, обозначить  \(\alpha=\frac {T_c}{T_a}= 0,000005..0,0005\) - малый коэффициент при условиях задачи,  то из уравнений, которые были выше, получим прилижённые значения и поправки

 \(\frac{dT}{T_0}=\alpha (1-4\alpha+O(\alpha^2))=\alpha+O(\alpha^2) \), a

 \(T_0=T_a(1+\frac{3}{2}(\frac{dT}{T_0})^2+(\frac{dT}{T_0})^4)^{-\frac{1}{4}}=T_a(1-\frac{3}{8}\alpha^2+O(\alpha^3)) \)
Все верно.
Будет время распишу подробнее.
Вам +

\(\displaystyle a=\frac{A_s~\sigma~h}{\lambda}\);

\(\displaystyle T_2=\left(\frac{I_0}{2 \sigma}\right)^{1/4}-\frac{1}{2} \left(\frac{I_0}{2 \sigma}\right) a+\frac{5}{8} \left(\frac{I_0}{2 \sigma}\right)^{7/4} a^2 +~ ...~+\).

« Последнее редактирование: 05 Июнь 2021, 18:39:10 от Ost »

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #9 : 06 Июнь 2021, 22:45:08 »
Тепловой поток от солнца \(q=A_s~I_0~S\).
Тепловой поток через пластину \(\displaystyle q_3=\lambda~\frac{S}{h}~(T_1-T_2)\), где \(T_1~-\) температура нагреваемой поверхности;
\(T_2~-\) температура охлаждаемой поверхности.

Пластина излучает с двух сторон
\(q_2=A_s~\sigma~S~{T_2}^4\);

\(q_1=A_s~\sigma~S~{T_1}^4\).

Для тепловых потоков справедливы такие равенства
\(q=q_1+q_2\)\(q_3=q_2\).

Получили систему уравнений
\(\displaystyle T_1^4+T_2^4=\frac{I_0}{\sigma}=b\);            (1)

\(\displaystyle T_1-T_2=\frac{A_s~\sigma~h}{\lambda}~T_2^4=a~T_2^4\).   (2)

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

\(\displaystyle T_0=\frac{T_1+T_2}{2}\).

...

« Последнее редактирование: 08 Июнь 2021, 21:39:56 от Ost »

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #10 : 08 Июнь 2021, 19:03:34 »

рис. 1
« Последнее редактирование: 10 Июнь 2021, 14:22:31 от Ost »

Оффлайн Ost

  • Местный мудрец
  • *****
  • Сообщений: 2399
  • Страна: ru
  • Рейтинг: +287/-29
Re: Найти температуру пластины в космосе
« Ответ #11 : 10 Июнь 2021, 16:58:22 »
Дополнительная задача.
Построить ряд Тейлора для функций \(T_0(a);~~T_1(a);~~T_2(a)\), при максимальном \(n=4\), используя (1) и (2).
\(\displaystyle T(a)=\sum_{k=0}^n \frac{1}{n!} \frac{d^n T}{d a^n}~a^n\).
« Последнее редактирование: 10 Июнь 2021, 17:09:23 от Ost »

Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #12 : 11 Июнь 2021, 18:41:36 »
Дополнительная задача.
Построить ряд Тейлора для функций \(T_0(a);~~T_1(a);~~T_2(a)\), при максимальном \(n=4\), используя (1) и (2).
\(\displaystyle T(a)=\sum_{k=0}^n \frac{1}{n!} \frac{d^n T}{d a^n}~a^n\).

Буду решать потихоньку. Для начала:
1) сделал читаемым ответ №3 в ветке и добавил там \(A_s\), чтобы формулы не отличались от Ваших и чтобы не ограничиваться черным металлом;
2) ввел сквозную нумерацию в ответах №3 и №7, так что формулы (1) и (2) совпадают с Вашими;
3) исправил потерянный  коэффициент при четвертой степени в (2), что все равно не влияло на ответ в разложении в ряд до второй степени по а, а при разложении до четвертой степени будет.

Далее, мне не очень удобно манипулировать с введенными коэффициентами a и b, т.к. b  имеет размерность \(T^4\), а a - \(T^{-3}\).
Поэтому переопределим:
B = (\(\frac{b}{2})^\frac{1}{4}=T_a\)
A=  (\(\frac{ab}{2})=T_c\),
где  \(T_a\) и  \(T_с\) - те же, что и в ответе №7 в ветке. Здесь они заменены на B и А для упрощения написания формул и преемственности с b и а. Оба имеют размерность температуры,  A<<B, а
\(\alpha=A/B=a(\frac{b}{2})^\frac{3}{4}\) - малый безразмерный параметр, как в ответе №7 - по нему будет удобно раскладывать в ряд Тейлора.
\(\beta=(\frac{dT}{T_0})\) - еще один малый безразмерный параметр (обозначение dT=\(T_1-T_2\) введено в ответе №3).

Задача сводится к разложению \(T_0 ~по ~\alpha,для ~чего~ сначала ~разложим ~ T_0~ по ~ \beta , а~ потом  ~\beta ~по~ \alpha \)

Т.к. (1) можно записать в виде
\(2B^4=T_0^4((1+\frac{\beta}{2})^4+(1-\frac{\beta}{2})^4)\),
То (4) переписывается как:
\(T_0=B(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{1}{4}}=B(1-\frac{3}{8}\beta^2+\frac{43}{128}\beta^4-\frac{375}{1024}\beta^6+O(\beta^8)))\).

Теперь осталось выразить  \(\beta~через~\alpha\). Так как  \(\beta\) входит в выражение для \(T_0\) только в четной степени, то для для разложения \(T_0\) по \(\alpha\) до второй/третьей степени достаточно того, что   \(\beta~=~\alpha+O(\alpha^2)\) (что и было сделано ранее). Для разложения  \(T_0\) по \(\alpha\) до пятой степени нужно два члена разложения  \(\beta~по~\alpha\), седьмой - три и т.д.
« Последнее редактирование: 12 Июнь 2021, 20:14:11 от Andrey_R »

Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #13 : 11 Июнь 2021, 22:14:37 »
Продолжение.

(2) перепишем в виде:
\(dT=aT_2^4=a(T_0-\frac{dT}{2})^4=aT_0^4(1-\frac{\beta}{2})^4\)
или, разделив на \(T_0\),
\(\beta=aT_0^3(1-\frac{\beta}{2})^4\)
С учетом того, что
(5)               \(T_0=B(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{1}{4}}\), имеем
\(\beta=aB^3(1-\frac{\beta}{2})^4(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{3}{4}}=\alpha(1-\frac{\beta}{2})^4(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{3}{4}}\),
а
(6)        \(\alpha=\beta(1-\frac{\beta}{2})^{-4}(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{\frac{3}{4}}\), и придется сделать обращение этой формулы, предварительно разложив её в ряд Тейлора:
\(\alpha=\beta+2\beta^2+\frac{29}{8}\beta^3+O(\beta^4)\).
Делая обращение, например, напрямую:
\(\beta=\frac{\alpha}{1+2\beta+\frac{29}{8}\beta^2+O(\beta^3)}=\frac{\alpha}{1+2\beta(1+\frac{29}{16}\beta+O(\beta^2)}=\frac{\alpha}{1+\frac{2\alpha(1+\frac{29}{16}\beta+O(\beta^2)}{1+2\beta+O(\beta^2)}}=\frac{\alpha}{1+\frac{2\alpha(1+\frac{29}{16}\alpha+O(\alpha^2)}{1+2\alpha+O(\alpha^2)}}\),
получаем
(7)            \(\beta=\alpha-2\alpha^2+\frac{35}{8}\alpha^3+O(\alpha^4)\)

Теперь можно подставить \(\beta\) в выражение для \(T_0\):
\(T_0=B(1-\frac{3}{8}(\alpha-2\alpha^2+\frac{35}{8}\alpha^3)^2+\frac{43}{128}\alpha^4+O(\alpha^5)))\),
и, если нигде не вкралась ошибка,
то получается вот такое:
\(T_0=B(1-\frac{3}{8}\alpha^2 +\frac{3}{2}\alpha^3-\frac{569}{128}\alpha^4+O(\alpha^5)))\)
Осталось вспомнить, что
B = (\(\frac{b}{2})^\frac{1}{4}\),
\(\alpha=a(\frac{b}{2})^\frac{3}{4}\)







« Последнее редактирование: 13 Июнь 2021, 10:32:35 от Andrey_R »

Оффлайн Andrey_R

  • Пламенный трибун
  • ****
  • Сообщений: 379
  • Страна: su
  • Рейтинг: +65/-5
Re: Найти температуру пластины в космосе
« Ответ #14 : 13 Июнь 2021, 12:04:51 »
Построить ряд Тейлора для функций T0(a);  T1(a);  T2(a)

Чтобы найти с заданной точностью \(T_1\) и \(T_2\), с учетом (5) получим:

(8)   \(T_1=T_0(1+\frac{\beta}{2})= B(1+\frac{\beta}{2})(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{1}{4}}\)
(9)   \(T_2=T_0(1-\frac{\beta}{2})= B(1-\frac{\beta}{2})(1+\frac{3}{2}\beta^2+\frac{1}{16}\beta^4)^{-\frac{1}{4}}\),
Поэтому нужно найти разложение \(\beta\) до той же степени по a, а не на 1 меньше, как для \(T_0\) 
Вручную это делать лень, а с помощью Maple (6) и его обращение (7) разлагаются так:
\(\alpha=\beta+2\beta^2+\frac{29}{8}\beta^3+\frac{19}{4}\beta^4+\frac{619}{128}\beta^5+\frac{271}{64}\beta^6+O(\beta^7)\)
\(\beta=\alpha-2\alpha^2+\frac{35}{8}\alpha^3-\frac{17}{2}\alpha^4+\frac{1419}{128}\alpha^5+\frac{65}{8}\alpha^6+O(\alpha^7)\)
После посдстановки этих значений в (5), (8) и (9) получается:
\(T_0=B(1-\frac{3}{8}\alpha^2 +\frac{3}{2}\alpha^3-\frac{569}{128}\alpha^4+\frac{41}{4}\alpha^5-\frac{15019}{1024}\alpha^6+O(\alpha^7)))\)
\(T_1=B(1+\frac{1}{2}\alpha-\frac{11}{8}\alpha^2 +\frac{7}{2}\alpha^3-\frac{969}{128}\alpha^4+\frac{45}{4}\alpha^5+\frac{3933}{1024}\alpha^6+O(\alpha^7)))\)
\(T_2=B(1-\frac{1}{2}\alpha+\frac{5}{8}\alpha^2 -\frac{1}{2}\alpha^3-\frac{169}{128}\alpha^4+\frac{37}{4}\alpha^5-\frac{33971}{1024}\alpha^6+O(\alpha^7)))\)
« Последнее редактирование: 13 Июнь 2021, 13:23:22 от Andrey_R »

Большой Форум

Re: Найти температуру пластины в космосе
« Ответ #14 : 13 Июнь 2021, 12:04:51 »
Loading...