CC BY
3Электронный журнал «Труды МАИ». Выпуск № 40
www.mai.ru/science/trudy/
УДК 536.22
Расчет эффективной теплопроводности элементарной ячейки с
помощью «действия».
Василевский Д.В., Симахин Е.А., Спирин Г.Г.
Аннотация
Предложено использовать «действие» (аналог интеграла энергии электростатического поля) для расчета эффективной теплопроводности элементарной ячейки. Предложенный подход применен для расчета теплопроводности дисперсных материалов, характеристики, необходимой при тепловом проектировании летательных аппаратов. Максимальное расхождение результатов расчета по сравнению с численным расчетом не превышает 10%.
Ключевые слова
Теплопроводность; «действие»; элементарная ячейка; дисперсная среда; температурное поле.
Метод обобщенной проводимости предполагает рассмотрение теплообмена в дисперсной среде в рамках одной элементарной ячейки, моделирующей основные черты переноса тепла в среде в целом. В 19 веке этот метод был предложен для расчета электрических и магнитных свойств дисперсных материалов [1], а затем распространен на тепловые свойства.
Расчет эффективной теплопроводности элементарной ячейки проводится численными методами [2], а в случае линеаризации температурного поля с помощью тепловых сопротивлений [3].
Вместе с тем, к расчету теплопроводности ячейки есть другой подход, основанный на аналогии стационарных электрических и температурных полей, описываемых уравнением Лапласа. Известно, что при заданном распределении зарядов энергия электростатического поля минимальна [4]. Отсюда можно говорить о принципе минимального действия для величины
где г-диэлектрическая проницаемость среды, ф - электрический потенциал. Действием применительно к теплопереносу является величина
5 = |ДУГ)2аV , (1)
V
X - коэффициент теплопроводности вещества.
Отметим, что принцип минимума «действия» может быть полезен при оценке точности численных расчетов стационарных температурных полей. В частности, с увеличением точности описания температурного поля, «действие» уменьшается, имея пределом минимальное значение, соответствующее истинному распределению температурного поля.
Остановимся на расчете эффективной теплопроводности элементарной ячейки при использовании «действия». Пусть элементарная ячейка, имеющая форму куба с ребром Ь, ограничена торцевыми изотермическими поверхностями с температурами Т1 и Т2 соответственно, и боковыми адиабатическими. Потребуем выполнения следующего равенства: 5 Ф = 5, (2)
где БЭф - «действие» для ячейки, заполненной однородным веществом с теплопроводностью ЯЭф, а 5 - «действие» реальной ячейки.
Рассмотрим элементарную ячейку, образованную чередующимися слоями, перпендикулярными тепловому потоку. В данном случае (2) имеет вид
2 f rm\2
Л
f T_т Л
12 1l
т ^ —
V L J i
dT
v dlr J
L2l, (3)
dTi _ 4o
где li - толщина i -ого слоя с теплопроводностью Л . Учитывая, что —= —, а
dh Л
T _ т q0
—-L = ——, (q0 = const), окончательно получим
L Лэф
f=2 Т w
Лэф i Л
_ li
Здесь да, — объемная концентрация / -ого слоя включения, I1'-объем l-ого слоя, I - объем ячейки.
В случае параллельно расположенных слоев относительно теплового потока имеем:
3 эф — Лэф
(Т Т 1 Ь3, 3 — 2ЛГ^ 1 '/Ь2. Отсюда
Лф - 2 Л
да, .
Соотношения (4) и (5) известны, они непосредственно следуют из закона Фурье, применимого к слоистым системам. Вместе с тем понятие «действия» имеет более общий характер и позволяет определить эффективную теплопроводность ячейки следующим интегральным соотношением:
И
дТ ]2 ( дт }2 (дт
дх
+
Лэф
V дУ ) Iдг.' У
2
Т2 -
(6)
V
Ь
Расчет теплопроводности существенно упрощается в случае линеаризациитемпера-турного поля, т.е. замене криволинейных линий тока прямыми. В качестве примера рассмотрим элементарную ячейку с включением сферической формы ( г - радиус включения, Л2 -
теплопроводность его материала, Л1 - теплопроводность среды, в которой находится включение) (рис. 1).
Сферу заключим внутри воображаемой адиабатической поверхности (цилиндра),
предполагая 3 эф — 3ад, имеем
'эф
л
( т1-т1
V Ь )
Ь —л
т - т
V Ь )
(( -пг2ь)+ Лэф(1,2)
т2 - т
V Ь )
(пг2 ь),
(7)
где Лэф(12) - эффективная теплопроводность цилиндра.
4 (г
Учитывая, что объемная концентрация включения т2 — -3- ж\ ^ I , перепишем (7) в ви
де
3
лад _ л
Лэф = Л1
1 -П — т2 I 4п
+ Лэф(и)П1^Пт2 )3
(8)
Далее рассмотрим «действие» для цилиндра 5^ 2) = 512, или
Л
эф(1,2)
( Т2_Т_У П2 Ь = Л
( ап ^2
^ил 1 у
п2 Ь - 4 п 3 | + Л
( ёТ. л 2
ёх
\ил2 у
—п3
,3
Отсюда имеем
1
1 — ( 3
1 —I — т2
3 I 4п
11
4 ( 3 ^3 -I — т.
Л,
3 I 4п
2
и'эф(1,2) Л
Объединяя (8) и (9) получим
2 '
Л
Дад л эф =Л
г
1 -п
4п
\
-т
2
3
+ ■
ЛЛ2
Л
1-
4
4п
-пи
■ П
+ Л
4
4п
■ т~
4п
/77-
(9)
(10)
Относительная величина эффективной теплопроводности равна
2
2
3
3
2
3
3
1
3
3
3
3
3
3
лад лэф
л
1 - п\ — т2 V 4п 2
п\ — т2 V 4п 2
1 4 ( 3
1 —\ — т2 3 V 4п
4 ( 3 ^ 3
+ V—\ — т0
3 V 4п
2
где V —
л
л
лад л
Рассмотрим отношение —эф. В случае--> да эта величина равна
л л2
ад лэф
л
1 -п\
4п
-пи
(12)
лад
Зависимость —эф (т2) в диапазоне 0 < т2 < —, представлена на рис. 3.
л 6
Рассмотрим случай расчета эффективной теплопроводности ячейки, для случая, когда в ячейке выделены две плоскости аа' и ЬЬ' с температурами соответственно т1' и т2' (рис. 2). По принятой терминологии это соответствует изотермическому дроблению ячейки.
Ь' Т2
а' Т1
!
ао
Рис. 2.
Имеем Б-иф — 3из, или
2
3
3
f Y
X
'эф
q0 L3 = X qo
л из \Лэф J VX J
f V
VX12 J
№),
(13)
Здесь X2 - эффективная теплопроводность слоя, заключенного между изотермическими поверхностями ad и ЬЬ'. Учитывая, что q0 = const, соотношение (13) перепишем в виде 1 1 "
лиз л
Хэф X
1 - 21
L
+—-— 21 Г
Х12эф V L,
(14)
Рассмотрим далее выделенный слой, для него имеем £12 эф = Би
X
2 эф
f т '_t 1Л f т —T
'2 T I-2rL2 = ' T T
\
2r
XI_nr22r)+X12эф
f T2' T \nr 22r ).
V
2r
X12 эф = X
1 _n( L j
+ M L J •
(15)
(16)
Здесь Л12эф - эффективная теплопроводность цилиндра, в котором заключен шар. Боковую поверхность цилиндра считаем адиабатической. В этом случае имеем 312эф = 312
2
q0
Х.2эф ~XX~ -nr 22r = х
V Х2эф J
1 = 1 2
X1' 2эф X1 3 X2 3
f qo Л2
VX J
n 2r _ з nr I + X2
f q£ Л2
VX2 J
— w I или
,3 J
(17)
Используя (14), (16), (17), получим
л из _
Хэф =■
X1 X1 2 1 f 3 Л 3 1 _nl — m2 I V 4П 2 J 2 3XX f 3 Л3 +-П — m2 I 2X + X2 V 4n J >
* X 2 1 f 3 1 3 1 _nl — m2 I V 4П 2 J 2 3\X2 f 3 Лз +-П — m2 I 2X + X2 V 4n J » - 1 1 _ 2|— m2 |3 V 4n 2 J 1 + 2x[— m2 13 V 4n 2 J
(18)
Относительная величина эффективной теплопроводности равна
X
'эф
1 _П — m2 V 4п
2
з f з Л з
+--П — m„
1 + 2v V 4п
2
1 _П — m2 V 4п
2
3 f 3 Л 3
+--П — m,
1 + 2v V 4п
2
1
f 3 Л3 1 _ 2l — m.
4П
2
1
+ 2| — m2 13
V 4п 2 J
Рассмотрим предельный случай. При V ^ да выражение (19) примет вид
2
2
2
4
> •
Л
1 13 I7
1 - п\ — т I
'эф
Л
3 13
1 - п\ — т I 4п
3 13
1 - 2\ — т I 4п
(20)
3 13
+ 2\ — т I 4п
Л
Зависимость —^ (т2), при V ^ ж представлена на рис. 3.
Л
Истинное значение теплопроводности Лэф лежит между Лф и Лф . Следовательно,
пиз лад
Лэф + Лэф
Лэф ^ Л Л Л ~ 2
Л
. Зависимость —эф (т2) также графически представлена на рис. 3. В том
Л
Л
числе, на этом же рисунке представлена зависимость —эф (т2), для случая сфера в кубе, поЛ
лученная численными методами [2]. Максимальное расхождение результата составляет не более 10%.
Представляется, что изложенный метод расчета теплопроводности дисперсных материалов может найти свое место в практике расчета теплофизических свойств материалов.
1.5
0.5
0.1
02 аз
т,
04 0:
Рис. 3.
4
2
2
2
2
Библиографический список
1.Maxwell C. Treatise on electricity and magnetism. Vol. 1.OxfordUniv. press, London
1892.
2.Дульнев Г.Н. Перенос тепла через твердые дисперсные системы. ИФЖ, Т.9, 1965, №3, с399-404.
3.Дульнев Г.Н. Заричняк Ю.П. Теплопроводность смесей и композиционных материалов, «Энергия», Ленинград , 1974, 265с.
4.Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике, т.6., «Мир», Москва, 1966, 342с.
Сведения об авторах
1)Василевский Дмитрий Валентинович, МАИ, в.н.с., к.т.н., 25080645.
2)Симахин Евгений Александрович, МАИ, студент, 89266254986.
3)Спирин Геннадий Георгиевич, МАИ, зав. кафедрой, д.т.н., профессор, 89857772487,spirinas@mail.ru.
