Математика к Математическое
моделирование
ХДК 536.2
Ссылка на статью:
// Математика и математическое моделирование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2016. №3. С. 18-28.
Б01:10.7463/шаШш.0316.0843802
Представлена в редакцию: 16.06.2016 © МГТУ им. Н.Э. Баумана
Температурное состояние пластины с зависящими от температуры теплопроводностью и энерговыделением
Зарубин В. С.1'*, Котович А. В.1, Кувыркин Г. Н.1
fn2@bmstu.ru
1МГТУ им. Н.Э. Баумана, Москва, Россия
На основе сформулированной нелинейной математической модели теплопроводности в пластине с зависящими от температуры коэффициентом теплопроводности материала пластины и объемной мощностью энерговыделения в аналитической форме представлено решение ряда конкретных задач, что позволяет выявить некоторые существенные эффекты, влияющие на температурное состояние пластины. Рассмотрены различные варианты зависимости объемной мощности энерговыделения от температуры. В частности, установлены условия, при которых в пластине невозможно существование установившегося распределения температуры, что соответствует явлению теплового взрыва. Полученные результаты могут быть использованы при тестировании компьютерных программ, используемых для численного решения нелинейных задач теплопроводности применительно к техническим устройствам, в объеме элементов которых происходит выделение или поглощение тепловой энергии
Ключевые слова: пластина; тепловой взрыв; температурное состояние; объемная мощность энерговыделения
Введение
Температурное состояние твердого тела помимо условий его теплообмена с окружающей средой может существенно зависеть от процессов тепловыделения (или теплопоглощения) в объеме тела. Среди возможных причин возникновения таких процессов следует отметить выделение энергии в тепловыделяющих элементах ядерных реакторов, протекание в материале твердого тела экзотермических или эндотермических химических реакций, сопровождаемых соответственно выделением или поглощением теплоты, переход в теплоту части электрической энергии при прохождении электрического тока через проводники (так называемая джоулева теплота) или энергии излучения, проникающего в тело из полупрозрачного материала, и т.п. Интенсивность этих процессов характеризуют объемной мощностью
Вт
энерговыделения ду, основной единицеи измерения которой является —. При поглощении
м3
теплоты в теле значение ду является отрицательным.
В обширной литературе по теории теплопроводности твердых тел [1, 2, 3] приведены решения задач по определению стационарного (установившегося во времени) и нестационарного температурного состояния твердых тел (как правило, канонической формы), в которых действуют источники объемного энерговыделения. При этом в общем случае учитывается возможность изменения величины ду по объему тела, а при решении нестационарных задач и возможная зависимость этой величины от времени.
Однако в реальных условиях величина ду часто зависит и от локального значения температуры, причем такая зависимость может быть нелинейной. Например, интенсивность выделения или поглощения теплоты при химических реакциях пропорционально скорости их протекания, которая, в свою очередь, чувствительна к значению температуры, причем зависимость от температуры имеет экспоненциальный характер [4]. Дополнительным фактором, усложняющим в таких случаях анализ температурного состояния твердого тела, является зависимость от температуры и коэффициента теплопроводности материала этого тела, особенно в случае существенной неоднородности распределения в нем температуры. Учет влияния перечисленных факторов приводит к необходимости использовать методы математического моделирования [5, 6, 7], позволяющие построить адекватную нелинейную математическую модель процесса теплопроводности в теле с объемным энерговыделением. Количественный анализ таких моделей требует, как правило, применения численных методов.
Вместе с тем для такой простой формы тела, какой является неограниченная пластина постоянной толщины, возможно при определенных допущениях получить в аналитическом виде решение нелинейной задачи теплопроводности, учитывающее зависимости от температуры коэффициента теплопроводности материала пластины и интенсивности энерговыделения. Такое решение позволяет выявить ряд существенных эффектов, влияющих на температурное состояние пластины, в том числе связанных с условиями существования установившегося распределения температуры, и может быть использовано для тестирования результатов, получаемых при помощи численных методов.
При действии в неограниченной пластине толщиной к внутренних источников (или стоков) тепловой энергии с объемной мощностью ду (г,Т), зависящей в общем случае как от координаты г, отсчитываемой по нормали к поверхности пластины, так и от температуры Т, одномерное распределение Т(г) температуры должно удовлетворять обыкновенному дифференциальному уравнению
где Л(г,Т) — коэффициент теплопроводности материала тела. В частном случае зависимости величины ду лишь от координаты можно найти первый интеграл этого уравнения в
1. Основные соотношения
(1)
виде
-(Т, г) ¿¿^ + 9о + / 9v(г) ¿г = 0, (2)
о
включающий постоянную 90, равную плотности теплового потока, поступающего в пластину через ее поверхность, от которой происходит отсчет координаты г.
Если для неоднородного тела коэффициент теплопроводности не зависит от температуры, но изменяется с изменением координаты г, то интегрирование равенства (2) приводит к соотношению
Т(г) = Т(0) - /Ч^о + / «V (С) (3)
о ( ) ^ о '
Отсюда при г = к следует
/ , г
т(0) - Т(к) -/ Ч^! 9V (г') ¿г' 9о =--. (4)
г ¿г
У -(г)
о
Для однородного тела, когда коэффициент теплопроводности зависит лишь от тепературы, после интегрирования равенства (2) получим
Т(г) г / \
| -(т) ¿Т = - ! ¿г'190 + / 9v(С) ¿() , (5)
Т(0) о ^ о '
Т(г) / г
9о к = I -(Т) ¿Т -У «V (г') ¿г'. (6)
Т( 0) о о
При граничных условиях первого рода по заданным значениям Т(0) = То и Т(к) = Т/ из формул (3) или (5) можно найти значение 90 и затем использовать это значение для расчета распределения температуры при помощи соотношений (4) или (6). Если граничные условия заданы в виде
¿Т
¿Т
--(Т(0))Тг + /о(Т(0)) = 0, -(Т(к)) — + /Н{Т(к)) = 0
¿г
х=/
ж=0
то значение 90 и неизвестные температуры Т(0) и Т(к) следует искать с учетом формул (4) или (6) из цепочки равенств
9о = -/о(Т(0)) = /Н(Т(к)) - ¡' qv(г) ¿г.
о
В частном случае идеальной теплоизоляции поверхности пластины при г = 0 имеем 90 = 0 и температуру Т(к) можно определить из равенства нулю правой части последнего соотношения, а затем по формуле
/ г
Т(0) = Т(к) + ! -^т/ 9v (г') ¿г'
п \ ' п
г
или
Т( 0) h z
f Л(T) dT = j dz f qv (z') dz'
T(h) 0 0
вычислить температуру T(0), что в итоге позволяет из соотношений (3) или (5) при q0 = 0 найти искомое распределение T(z) температуры. В случае постоянный величин Л и qv из последних формул следует T(0) — T(h) = qvh2/(2Л).
2. Зависимость энерговыделения от температуры
Если величина qv зависит от температуры и не зависит от координаты z, то подстановкой p = Л^)уравнение (1) можно привести к виду p+ qv(T)Л(T) = 0 и после
dz dT
интегрирования записать
T
p2 = Ci — 2 j qv (T') Л^) dT,
T*
где C1 — постоянная интегрирования, а значение T* может быть выбрано произвольно. Из этого равенства следует обыкновенное дифференциальное уравнение с разделяющимися переменными
dT / T \1/2
Л^) — = ± (ci — 2j qv (T' )Л^) dTj . (7)
т„
Выбор знака в правой части зависит от направления теплового потока в пластине. Интегрированием уравнения (7) получим
Т ( T Л-1/2
J Л^ )l Ci — 2 J qv (Г)Л^) dTj dT. (8)
TT
В случае задания на поверхностях стенки граничных условий первого рода целесообразно принять Т* = Тн, что приводит к равенству С2 = Тогда постоянную С\ можно найти из условия Т(г) = Т0 при г = 0, равносильного равенству
Tí ( T \-1/2 h ± Л(TД Ci — 2 qv(T')Л(Г) dT'j dT = 0. (9)
Tfc Th '
Ясно, что при произвольных зависимостях ду и Л от Т для поиска значения С\ необходимо использовать последовательные приближения.
Если при заданной температуре Т(к) = Тн поверхность пластины при г = 0 идеально теплоизолирована, то из формулы (7) с учетом выбора Т* = Тн следует
То
Ci = 2j qv (T') Л(T) dT,
Th
а вместо равенства (9) получим уравнение
Т( 0) . Т( 0)
Т0) ( Т0) .-1/2
к = ± Х(Т ) 2 / ду (Т')Х(Т') СТ') СТ
Т ^ Т '
для определения температуры Т(0), решение которого в общем случае также возможно лишь путем последовательный приближений. Уравнение (10) имеет смысл при положительном значении внутреннего интеграла, т.е. при ду (Т) > 0 должнобыть Т (0) > Т, и перед внешним интегралом следует выбрать знак плюс, а при ду (Т) < 0, наоборот, Т(0) < Тн и перед этим интегралом необходимо выбрать знак минус.
Положим ду (Т) > 0 и используем в уравнении (10) подстановку Сф = Х(Т) СТ, приняв
Т(0)
ф(0) = У Х(Т) СТ.
Тогда уравнение (10) примет вид
■ф(0) . ф(0)
-1/2
к = I 12 У ду(ф') сф'\ Сф
(11)
0 х -ф
где ду (ф) — зависимость обыемной мощности энерговыделения от ф, получаемая переходом в функции ду (Т) от температуры к функции ф в соответствии с подстановкой
ТТ
ф = I Х(Т) СТ.
Представим эту зависимость в виде ряда Маклорена [8]
Сд(ф)
д(ф) = д(0) +
Сф
ф +1 Яду (ф)
ф2 +
(12)
ф=0
ф=0 2 Сф
Если ограничиться двумя первыми членами ряда (12), то после подстановки в уравнение (11) и интегрирования получим
ф(0) = ду (0)к2
вес
V. А - 1
А
(13)
где А = к
2 ¿ду (ф)
йф
ф=0
. Очевидно, что при А ^ п2/4 правая часть формулы (13) не-
ограниченно возрастает, т.е. существование установившегося распределения температуры в пластине возможно лишь при условии
Сду (ф)
или
Сф Сду (Т)
ф=0
1( ^ 2 < 4Кк
СТ
Т=Тн
п
< Х<Гh)4кí•
(14)
При более высоком темпе возрастания объемной мощности энерговыделения с увеличением температуры существование установившегося распределения температуры в пластине невозможно. В этом случае малые приращения температуры вызывают рост выделения тепловой энергии, которое уже не может быть отведено к поверхности пластины с температурой Th путем теплопроводности без дальнейшего роста температуры. В итоге температура в пластине неограниченно возрастает, что соответствует явлению теплового взрыва. Для тел простой формы условия теплового взрыва исследованы в работе [4].
При выполнении условия (14) и использовании первых двух членов ряда (12) для рассматриваемой пластины с одной идеально теплоизолированной поверхностью из соотношения (8) следует
T(z)
ft, ITL )h2 / cos \
(15)
*>= /A(T)dT = ^ - 1)
Th
Отсюда как частный случай при A = 0 предельным переходом можно получить формулу
«V (Th)ft.2(i- (h )2), которая при A = const и qV = const переходит в известное равенство [9]
T(z) = T + «V(h2 - z2)
(z ) = Th +-2A-•
В случае A < 0 аргумент косинусов в соотношении (15) будет чисто мнимым и их следует заменить гиперболическими косинусами, записав
«и^ (i - ch #
A V ch Л Г\A
3. Экспоненциальная зависимость энерговыделения от температуры
Энерговыделение в пластине может быть вызвано протекающими в ее материале химическими реакциями, скорость которых описывает закон Аррениуса [4]. В этом случае зависимость объемной мощности энерговыделения от температуры можно представить в виде
9У (Т ) = ехр(- кЬ),
где 9V — предэкспоненциальный множитель; и — энергия активации химической реакции; Дж
к0 ~ 1,38 • 10-23 —--постоянная Больцмана. При 9у > 0 реакция является экзотермиче-
К
ской, а при 9у < 0 — эндотермической.
Если для пластины с постоянным коэффициентом теплопроводности -, идеально теплоизолированной поверхностью при г = 0 и заданной температурой Т/ на поверхности при г = к допустимо принять Т(г)/Т/ - 1 ^ 1, то дифференциальное уравнение (1), описывающее распределение температуры в такой пластине, можно привести к виду
¿2|Р + в ехр(0(С)) =0, (16)
где
< = * ) = т
коТ
7 =
о Тн
о * ( 1
в = ехр — 7>Т» V 7
Обозначив аЭ(г) = у, преобразуем уравнение (16) к виду У + в ехр(в) = 0 и после
ах аЭ
интегрирования получим У2 = Б. — 2в ехр(В), или
^ = — 2в ехр(в(С)),
(17)
где Б. — постоянная интегрирования. При цу > 0 (в > 0) направление теплового потока в данной пластине совпадает с направлением оси Oz. Поэтому в правой части равенства (17) следует выбрать знак минус. Из условия ОЭ = 0 при £ = 0 из этого равенства находим Б. = 2в ехр(в(0)^. В итоге интегрирование уравнения (17) с использованием граничного условия в(1) = 0 приводит к соотношению
4~в (1 — с),
(18)
где X = ехр(в(С)) иХо = ехр(6(0)
Для определения неизвестного значения Х0 из соотношения (18) при ( = 0 получаем трансцендентное уравнение
— Хо ^в
(19)
график зависимости Ш(Х0) левой части которого от Х0 представлен на рис. 1 (для сплошной кривой шкалы осей абсцисс и ординат расположены соответственно сверху и слева, а для штриховой — снизу и справа). Из этого графика и уравнения (19) следует, что существует
1,0
Ш
хк
0,6
0,2
0
Хп
X'
10
/ 1— .... --- --- » 1 ' 1
/
/
/
0,4
0,3
0,2
0,1
1 1,01 1,02 1,03 1,04 1,05 1,06 1,07 1,08 X
Рис. 1. График зависимости левой части уравнения (19) от Хо
наибольшее значение в * ~ 0,878 этого параметра, равное квадрату значения ординаты горизонтальной касательной к графику. Значения в > в * соответствуют тепловому взрыву в рассматриваемой пластине, а при 0 ^ в ^ в * в пластине может существовать установившееся распределение температуры, описываемое соотношением (18). При этом каждому значению в ° < в * формально отвечают два значения Х'0 и Х0, из которых может быть реализовано лишь меньшее. Предельно возможное значение X* ~ 3,274, равное абсциссе точки касания горизонтальной касательной к графику зависимости Ш(Х0), соответствует значению в* (0) « 1,186.
При < 0 (в < 0) направление теплового потока в пластине противоположно направлению координатной оси Oz. В этом случае направление теплового потока в пластине противоположно направлению координатной оси Oz. Поэтому в правой части равенства (17) следует выбрать знак плюс и после интегрирования полученного дифференциального уравнения с учетом граничного условия в(1) = 0 записать
Хо( агс^у Хо-1 - агсЧ Хо- ^ =
1 - С).
(20)
Для нахождения неизвестного значения Х0 из равенства (20) при ( = 0 получим транс-цендентое уравнение
'Х ^ Хо-1=
(21)
график зависимости ^1(Х0) левой части которого от Х0 приведен на рис. 2 (для сплошной кривой шкалы осей абсцисс и ординат расположены соответственно сверху и справа, а для штриховой — снизу и слева). Этот график по известному значению |в| позволяет определить единственное значение Х0 < 1, что дает возможность по формуле (18) построить распределение температуры в пластине. При этом следует отметить, что в силу отрицательности значения в(0) температура Т(0) на теплоизолированной поверхности пластины будет меньше температуры Т^.
А 100 80 60 40 20 0 ( »VI 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0Л
8 6 4 2 0 0
\
\
- - -
) 0,( )1 0,( )2 0,( )3 0,( )4 0,( 35 0,( 36 0,< 37 0,( 38 0, 39 X
Рис. 2. График зависимости левой части уравнения (21) от X
Заключение
На основе сформулированной нелинейной математической модели теплопроводности в пластине с зависящими (в общем случае) от температуры коэффициентом теплопроводности материала пластины и объемной мощностью энерговыделения в аналитической форме решен ряд конкретных задач, что позволяет выявить некоторые существенные эффекты, влияющие на температурное состояние пластины. В частности, установлены условия, при которых в пластине не может существовать установившееся распределение температуры. Полученные результаты могут быть использованы при тестировании компьютерных программ, используемых для численного решения нелинейных задач теплопроводности.
Работа выполнена в рамках проекта 1712 по государственному заданию 2014/104 Минобрнауки РФ и государственного задания по проекту 1.2640.2014.
Список литературы
1. Карслоу Г., Егер Д. Теплопроводность твердых тел: пер. с англ. М.: Наука, 1964. 488 с.
2. Лыков А.В. Теория теплпроводности. М.: Высшая школа, 1967. 600 с.
3. Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел. М.: Высшая школа, 2001. 550 с.
4. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. М.: Наука, 1987. 491 с.
5. Зарубин B.C. Моделирование. М.: Издательский центр «Академия», 2013. 336 с.
6. Зарубин B.C., Кувыркин Г.Н. Особенности математического моделирования технических устройств // Математическое моделирование и численные методы. 2014. Т. 1, № 1(1). C. 5-17.
7. Зарубин B.C., Кувыркин Г.Н. Математическое моделирование термомеханических процессов при интенсивном тепловом воздействии // Теплофизика высоких температур. 2003. Т. 41, №2. C. 300-309.
8. Власова Е.А. Ряды. М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. 616 с.
9. Теория тепломассообмена / Под ред. А.И. Леонтьева. М.: Изд-во МГТУ им. Н.Э. Баумана, 1997. 683 с.
Mathematics i Mathematical Modelling
Electronic journal of the Bauman MSTU
Mathematics and Mathematical Modelling of the Bauman MSTU, 2016, no. 3, pp. 18-28.
DOI: 10.7463/mathm.0316.0843802
Received: 16.06.2016
© Bauman Moscow State Technical University
http://mathmjournal.ru
The Temperature Condition of the Plate
with Temperature-Dependent Thermal Conductivity
and Energy Release
Zarubin V. S.1'*, Kotovich A. V.1, Kuvyrkin G. N.1 *fii2@bmstu.ru
1 Bauman Moscow State Technical University, Russia
Keywords: slice, thermal explosion, temperature state, the volume power of energy release
The temperature state of a solid body, in addition to the conditions of its heat exchange with the environment, can greatly depend on the heat release (or heat absorption) processes within the body volume. Among the possible causes of these processes should be noted such as a power release in the fuel elements of nuclear reactors, exothermic or endothermic chemical reactions in the solid body material, which respectively involve heat release or absorbtion, heat transfer of a part of the electric power in the current-carrying conductors (so-called Joule's heat) or the energy radiation penetrating into the body of a semitransparent material, etc. The volume power release characterizes an intensity of these processes.
The extensive list of references to the theory of heat conductivity of solids offers solutions to problems to determine a stationary (steady over time) and non-stationary temperature state of the solids (as a rule, of the canonical form), which act as the sources of volume power release. Thus, in general case, a possibility for changing power release according to the body volume and in solving the nonstationary problems also a possible dependence of this value on the time are taken into consideration.
However, in real conditions the volume power release often also depends on the local temperature, and such dependence can be nonlinear. For example, with chemical reactions the intensity of heat release or absorption is in proportion to their rate, which, in turn, is sensitive to the temperature value, and a dependence on the temperature is exponential. A further factor that in such cases makes the analysis of the solid temperature state complicated, is dependence on the temperature and the thermal conductivity of this body material, especially when temperature distribution therein is significantly non-uniform. Taking into account the influence of these factors requires the mathematical modeling methods, which allow us to build an adequate nonlinear mathematical model of the heat conductivity process in the volume power release body. Quantitative analysis of these models requires using the numerical methods, as a rule.
At the same time, such a simple body, which is an unlimited plate of the constant thickness allows us, under certain assumptions, to solve analatically a nonlinear heat conductivity problem taking into account the thermal conductivity of the plate material and the power release intensity versus temperature.
This solution enables us to reveal a number of significant effects that have impact on the thermal state of the plate, including those related to conditions of available steady temperature distribution, and it can be used to test the results obtained by numerical methods.
References
1. Carslaw H.S., Jaeger J.C. Conduction of Heat in Solids. Oxford University Press, 1959. (Russ. ed.: Karslou G., Eger D. Teploprovodnost' tverdykh tel. Moscow, Nauka Publ., 1964. 488 p.).
2. Lykov A.V. Teoriia teploprovodnosti [Heat conducton theory]. Moscow, Vysshaya shkola Publ., 1967. 600 p. (in Russian).
3. Kartashov E.M. Analiticheskie metody v teorii teploprovodnosti tverdykh tel [Analytical methods in heat conduction theory of solid bodies]. Moscow, Vysshaya shkola Publ., 2001. 550 p. (in Russian).
4. Frank-Kamenetskii D.A. Diffuziia i teploperedacha v khimicheskoi kinetike [Diffusion and heat transfer in chemical cinetics]. Moscow, Nauka Publ., 1987. 491 p. (in Russian).
5. Zarubin V.S. Modelirovanie [Modeling]. Moscow, Publishing Center "Academy", 2013. 336 p. (in Russian).
6. Zarubin V.S., Kuvyrkin G.N. Special features of mathematical modeling of technical instruments. Matematicheskoe modelirovanie i chislennye metody = Mathematical Modeling and Computational Methods, 2014, vol. 1, no. 1(1), pp. 5-17. (in Russian).
7. Zarubin V.S., Kuvyrkin G.N. Mathematical modeling of thermomechanical processes under intense thermal effect. Teplofizika vysokikh temperatur, 2003, vol.41, no. 2, pp. 300-309. (in Russian). (English version of journal: High Temperature, 2003, vol.41, no. 2, pp. 257-265. DOI: 10.1023/A:1023390021091)
8. Vlasova E.A. Riady [Series]. Moscow, Bauman MSTU Publ., 2006. 616 p. (in Russian).
9. Leont'ev A.I., ed. Teoriia teplomassoobmena [The theory of heat and mass transfer]. Moscow, Bauman MSTU Publ., 1997. 683 p. (in Russian).