Моделирование процесса передачи тепла в пористом материале МР Текст научной статьи по специальности «Физика»

УДК 536.071.7

МОДЕЛИРОВАНИЕ ПРОЦЕССА ПЕРЕДАЧИ ТЕПЛА В ПОРИСТОМ МАТЕРИАЛЕ МР

©2011 А. М. Жижкин

Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)

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

Модель, теплопередача, материал МР.

Настоящая статья является продолжением исследований, начатых в работе [1], целью которых является исследование и моделирование процесса передачи тепла в материале МР. Для проверки принятых гипотез и допущений, определения численных значений коэффициентов в аналитических моделях были использованы результаты экспериментальных исследований эффективной теплопроводности (в дальнейшем теплопроводности) материала МР, приведенные в работах [1,2].

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

При исследовании процесса теплопередачи в пористых материалах используется их тепловая характеристика -эффективная теплопроводность.

При определении эффективной теплопроводности материала МР были приняты следующие ограничения и допущения:

- значения параметров теплопрово-дностей материала проволоки X], вещества в порах Л.2 и пористости П у модели и структуры - одинаковые;

- структуру МР считать изотропной;

- агрегатное состояние компонент при передаче тепла в модели и структуре не изменяется;

- теплопроводность компонент модели не зависит от температуры;

- изменение температуры вдоль направления теплового потока линейно;

- массообмен внутри порового объема отсутствует.

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

Пн , (1~пн)

к

(1)

где Хэн и 1э - теплопроводности элементов с непроводящими и проводящими тепло контактами;

I/ I/, а,, - параметр, определяющийся отношением количества элементов с теплопроводностью Хэн к общему количеству элементов в цепи.

Теплопроводность /..,„ можно определить с помощью модели пористой структуры с взаимопроникающими элементами, полученной в работе [3] и приведенной в виде

ХЭН=Х}[С2 + Х^1-С)2 + +2Х0С(1-С)/(Х0С+1-С)~1], (2)

где Хо = /-2//4 - относительная теплопроводность (>ч и /~2 - теплопроводности материалов скелета и порового пространства соответственно); С = Ар /Ьр - относительный размер бруса (2Ар - ширина бруса

упорядоченной решетки; 2Ьр - длина ребра куба).

Зависимость для определения теплопроводности Хэ приведена автором в статье по модели структуры МР, публикуемой в настоящем сборнике в виде выражений:

=Х1[СгМ+Х0(1-Су

2к.,С(!-С) (/.„С 1-С) !/.

(3)

М = (Х1 и С- 4Кк + I)-1, где Кк - термическое сопротивление контактов проволок в материале МР.

Величина теплового сопротивления И, в МР может быть приближенно определена с помощью выражения, предложенного там же в виде

а-К)2

я.

(4)

тгб/||р/.2[/о-1 1пП//,,]'

где ёпр - диаметр проволоки.

Используя формулы для термического сопротивления теплопроводной модели структуры (1), теплопроводностей, составляющих структуру элементов (2), (3) и зависимость для определения термического сопротивления контакта в МР (4), можно получить выражение для определения теплопроводности материала МР в виде

зк.са-сг1

А — Aj

C'+kj1-СУ

AC+l-C

d-Ш)

СгМ + XJ\-C)2

2КС(\ С) АС+1-С

(5)

Используя формулы (3) и (4), можно записать выражение для параметра М в виде

М

4С~(1- Ао)2

■1

(6)

Из представления единичного объема МР можно определить параметр

пц = рГТз. Число контактов на виток может быть рассчитано по формуле

"К=4Д (7)

где ро = рмрфм = (1-П) - относительная плотность материала, рмр - плотность МР

Число контактов в единице объема МР равно

NK = nK N6 Щ = Рмр /пц , где рмр = (1-П) рм,- тв равна mr; (nd„ )2Дсрм единице объема.

(8)

— масса витка, которая 4. Nв - число витков в

Тогда число витков в единице объема

равно

Ne= 4(1-П):ШпГДс

выражения

С учетом последнего зависимость (8) примет вид

• (9)

материала МР

NK = 16(1-П)2 /ж2 dn3

Из модели

известно, что

Число определить выражением

п..

и.

структуры

= WJ3.

контактов в с учетом

цепи можно формулы (9)

_2_

ж

1/3

(10)

В работе [3] приведена зависимость между параметром С и объемной концентрацией, наполняющей компоненты (для МР это пористость IT):

П = 2С3 - ЗС2 + 1. (11)

Для расчета параметра С рекомендовано использовать первый корень уравнения (11):

С =0,5+A cos(p/3 , 270 <ср< 360 при 0 < П< 0,5 А -1. <р = arc cos(I-2II), при 0,5 < П < 1 А= 1, (р = arc Ш(2П-1).

Выражение (6) с учетом зависимостей (10) и (11) можно записать в виде

М

Ямр Вт/мК

\ (12С2)1,3(\- ч/

{[ф-2С)2]'3ф0- -Hhfl/Xj]

1 .(12)

Рис. 1. Коэффициент эффективной теплопроводности МР: 1,2,3,4 - расчетные зависимости при значениях параметра п„ соответственно 0; 0,5; 0,75; 0,8;

=14,6 Вт/мК А> = 2,19.1а3;

□ - экспериментальные данные для МР из проволоки ЭП-708А

2

-1

-1

Для материала МР, как было показано выше, параметр и = 0,5. Тогда зависимости (5) и (6) с учетом (12) можно представить в виде выражения

X — X.

м

С2+Хо(1-С)2

21оС(1-С) ХоС + 1-С

{1-п\с2М+Х (1-С)1+2Х°С(Х~С) V А и / ХоС + \-С

1,56 С2!ъ(\-Хо)2

,(13)

1

ГЗ-2С )21ъХо[Хо-\ + \п(\1Хо)] для С <0,5 (77 > 0,5).

Рассмотрим предельные случаи значений теплопроводности компонент, параметров контакта и значений пористости для теплопроводности модельной структуры МР по зависимостям (13).

Пусть пористость материала МР /7=1. В этом случае С = 0, М= 1, а теплопроводность Хмр = Аг-

В случае, когда в порах модели непроводящая компонента (вакуум), то имеем X], Х2 = 0, Х0 = 0. В этом случае М = 0 и теплопроводность Хмр = 0.

Если теплопроводность компонент одинаковая X] = Х2 = X, а Х0 = 7, то М = 7, при этом теплопроводность структуры модели равна Хмр = X.

Теплопроводности, полученные по зависимостям (13) при предельных значениях параметров модели, физически обоснованы, что свидетельствует о правдоподобности исходных предпосылок при построении модели структуры и выбранном способе описания процесса переноса тепла.

На рис. 1 представлены результаты расчета теплопроводности МР по зависимостям (13) при разных значениях параметра п. Точками показаны экспериментальные значения теплопроводности МР.

Расчетные данные, полученные при параметре п=0,5, хорошо сочетаются с экспериментальными значениями теплопроводности МР. Характер изменения расчетной кривой отличается в пределах погрешности измерений теплопроводности от экспериментальных данных. Такое различие может быть обусловлено неучетом в модели структуры МР закономерностей образования пятна контакта элементов скелета, а также возможным влиянием

пористости 77 на относительное количество непрерывных элементов п.

Зависимости (13) могут быть использованы для анализа и определения теплопроводности МР в диапазоне изменения пористости от 0,5 до 1.

С помощью разработанной модели теплопроводности можно проанализировать влияние параметров структуры материала МР на его теплопередающие свойства. Для обобщения влияния различных факторов на теплопередающие свойства МР используем её относительную теплопроводность.

Для этого, используя зависимость для вычисления коэффициента эффективной теплопроводности материала МР (13), получим выражение для вычисления параметра относительной теплопроводности в виде:

С2 +10(1-С)г +

(1-й)

С2м + Л0(1-С)2 +

2\оС(1-С) ЛоС + 1-С

2 ХоС(1-С)

.(14)

хос+1-с

Из анализа формул (13), (14) можно сделать вывод, что Х=/(С, п„, /.,,), где С, //„ -параметры структуры модели, а Х0 -относительная теплопроводность, характеристика соотношения теплопроводностей составляющих структуру МР компонент. Параметр модели С связан с параметром структуры МР 77 соотношением (11). Параметр структуры модели пн связан с особенностями строения структуры МР.

На рис. 2 представлены интерпретации зависимости теплопроводности от пористости 77 при разных значениях параметров пн и >ч>-

Анализ результатов расчета показывает, что теплопроводность МР существенно зависит от параметра пн при малых значениях параметра /,„ (Х0—>0). Зависимость теплопроводности от пористости при значениях параметра //„—> 0 становится малосущественной. Ее изменение составляет 50... 80% при изменении пористости 77 от 0,9 до 0,5. При значениях параметра п —► 1 зависимость теплопроводности от пористости становится существенной. Она изменяется более чем в 6 раз при изменении пористости 77от 0,9 до 0,5.

ЛоМРУСГ 50

а

05

0,1

п

0.7

Я

в

Рис. 2. Влияние параметров структуры на эффективную теплопроводность MP:

а - Хп

0; 2-п,

2,19 1(Г3; 4

1 -пя

2,19 Ж; 1

б -Ха в - Ха = 2,19 10 г - 1 — /?„ = 0; 2 - п„ = О,5; 3 - пн = 1 при Хо1

Влияние параметра Х0 на теплопроводность МР иллюстрировано на рис. 2,г при трех значениях параметра пн (1 - п„ = 0; 2 а,, 0,5; 3 п,1 1).

Анализ результатов расчетов позволяет сделать вывод о влиянии параметра л0 на относительную теплопроводность Хмр0 в зависимости от значений параметра пн. С ростом значений параметра 10 влияние параметра п„ на теплопроводность МР резко уменьшается.

Так при значениях Х0—>0,2 влиянием параметра п„ на теплопроводность МР можно пренебрегать, при этом погрешность вычисления эффективной теплопроводности не превысит 30% в рассматриваемом диапазоне изменения пористости П.

, = 0,5; 3 — п„ = 0,75; 4 п„ = 0,8; 5 - п„ = 1; 0; 2- п „ = 0,5; 3 -пн = 1;

п„

0,8;

2,19 Iff3; Хо2 = 0,2; Хо3 = 0,5; Хо4 = 0,8

основании постановления Правительства РФ №218 от 09.04.2010.

Работа выполнена при финансовой поддержке Правительства Российской Федерации (Минобрнауки), на

Библиографический список

1. Жижкин, A.M. Анализ физической картины процесса теплопередачи в материале МР [Текст] / A.M. Жижкин, А.И. Белоусов // Математическое моделирование информационных процессов и систем в науке, технике и обществе: межвуз. сб. науч. трудов. - Самара, 2004. - С. 43 - 52.

2. Дульнев, Г.Н. Теплопроводность смесей и композиционных материалов: справочная книга [Текст] / Г.Н. Дульнев, Ю.П. Заричняк - Л.: Энергия, 1974. - 264 с.

3. Жижкин, A.M. Экспериментальные исследования эффективной теплопроводности образцов из материала МР [Текст] / A.M. Жижкин // Вестн. Самар. гос. аэрокосм, ун-та. Вып. 2. -Самара: 2006. -4.2. - С. 15-18.

MODELING OF HEAT CONDUCTION PROCESS IN POROUS MR MATERIAL

©2011 A. M. Zhizhkin

Samara state aerospace university named after academician S. P. Korolyov (National Research University)

Principles of modeling of structure of MR material and heat conduction process for absence of convection of media in holes are presented in this article.

Model, heat conduction, material MR.

Информация об авторах

Жижкин Александр Михайлович, кандидат технических наук, доцент кафедры конструкции и проектирования двигателей летательных аппаратов, Самарский государственный аэрокосмический университет Тел.: (846) 267-46-87. Область научных интересов: виброизоляторы и уплотнения на основе материала MP.

Zhizhkin Alecsandr Michailovich, Candidate of Technical Sciences, Associate professor of Aircraft Engine Design Department, Samara state aerospace Phone: (846) 267-46-87. Area of research: sealing made of MR material.