CC BY
40
ЭЛЕКТРОТЕХНИКА И ЭНЕРГЕТИКА
УДК 536.2.01
СИЛЬНАЯ ФИЗИЧЕСКАЯ НЕЛИНЕЙНОСТЬ И ПРОСТРАНСТВЕННАЯ НЕОДНОРОДНОСТЬ ТЕПЛОВЫХ ПОЛЕЙ В ОДНОКОМПОНЕНТНЫХ И ДВУХКОМПОНЕНТНЫХ СИСТЕМАХ
О. Н. ШАБЛОВСКИЙ, Д. Г. КРОЛЬ, И. А. КОНЦЕВОЙ
Учреждение образования «Гомельский государственный технический университет имени П. О. Сухого»,
Республика Беларусь
Введение
Фундаментальные уравнения теории теплопереноса (волновое уравнение, уравнение теплопроводности), содержащие нелинейные по температуре источники, позволяют моделировать разнообразные физические процессы, происходящие при энергетическом воздействии на материал. В настоящей работе рассматриваются закономерности формирования температурных полей в однокомпонентных и двухкомпонентных теплофизических системах, испытывающих конкурентное воздействие объемных источников энергии. Такие задачи являются важным элементом динамической теории неравновесных состояний вещества [1]. Данная работа продолжает исследования [2], [3] и имеет следующие цели: 1) изучить качественные и количественные свойства систем «среда - источник энергии», обладающих кубической, тригонометрической и экспоненциальной нелинейностями по температуре; 2) дать примеры аномального температурного отклика среды на тепловое воздействие («отрицательная теплоемкость»); 3) проанализировать тепловое состояние системы, в которой действует неоднородный по волновой координате источник энергии.
Волновое уравнение теплопереноса
Волновое уравнение теплопереноса с источником энергии имеет вид:
д2Т 2 д2Т ------ = ^
дг2 дх2 и
= + ^, (1)
где х - декартова координата; г - время; Т - температура; м>2 = Х/(су) - квадрат скорости распространения тепловых возмущений; X - коэффициент теплопроводности; с - объемная теплоемкость; у - время релаксации теплового потока; ^ /(су); qn - мощность внутренних источников тепла. Для двухкомпонентной
системы два зацепляющихся волновых уравнения с неоднородными по координате источниками имеют следующий вид:
д2Т д2Т
+е>(х,7;,Т2), о = 1,2, (2)
где Т1, Т2 - температуры взаимодействующих друг с другом компонентов сплошной среды; w1, w2 - две скорости распространения тепловых возмущений. В классе автомодельных решений типа распространяющейся волны £ = А*х + В*г; А*, В* - сом^ та = та(£) = Т(£) -Та°, уравнения (2) имеют вид:
(3)
(4)
где функции та характеризуют отклонения от равновесных температур Т° = сonst первой и второй компонент. Скорость перемещения £ -линии равна N = dx/dt = -В*/А*, причем N2 ^w2, т. е. £ - линия не является характеристикой; МI есть квадрат теплового числа Маха. С точки зрения аналитической механики функции Qa (т1з т2) - это нелинейные консервативные силы, которые действуют в автономной динамической системе с двумя степенями свободы. Если для обеих компонент процесс дозвуковой (М2 < 1) или сверхзвуковой (М2 > 1), то имеем инверсию: знаки функций ^а) в (4) меняются (при одних и тех же Qa) на противоположные при переходах «дозвук о- сверхзвук». Двухкомпонентную систему назовем контрастной, если М12 < 1; М2 > 1 или М12 > 1; М22 < 1.
Частным случаем (3) является динамическая система, соответствующая волновому уравнению (1):
Свойство инверсии проявляется в (5) следующим образом:
1) Q(T) > °, если в сверхзвуковом процессе (М2 > 1) имеем источник энергии & (Т) > 0), либо если в дозвуковом процессе (М2 < 1) имеем сток энергии
2) Q(T) < 0, если М2 > 1, ^ (Т) < 0, либо если М2 < 1; ^(Т) > 0.
Кроме того, имеем два важных частных случая:
1) одномерный стационарный процесс при А* = 1; В* = 0, т. е. £ = х;
2) нестационарный и однородный по координате х тепловой процесс при А* = 0;
Тригонометрическая нелинейность источника
Для уравнения (5) рассмотрим два примера источников, вид которых характерен для двойного уравнения синус-Г ордона.
Если Q(t) = - sin т + 2m sin(x /2), то
Точные решения (6) и (7) получены в работе [4]. Отличие между тепловыми режимами (6) и (7) в том, что решение (6) дает колебания, разрывные по отношению
d2т/d£ = Q(t); Q(t) = К(T)/[w2A*2(M2 -1)]; M2 = N2/w2; т = T - T0; т1 =т(£); T2 = °.
(5)
((К (T) < 0);
B* = 1, т. е. £ = t.
(6)
Если Q(t) = - sin т + 2(m - 1)sin(T /2), то
0 < m < 1.
(7)
к аргументу £ (рис. 1, а, б); линии, представляющие функцию Q(t), являются здесь фрагментами синусоид. Вместе с тем, решение (7) является непрерывным, а зависимость Q(t) имеет синусоидальный тип (рис. 1, в, г). Отметим здесь вариант, изображенный на рис. 1, д: монотонный по температуре источник с отрицательным наклоном dQ/di< 0 возбуждает тепловое поле, периодическое по аргументу £. Это
означает, что знакопеременный источник энергии kv (т) с отрицательным наклоном dkv(т)/dT< 0 наблюдается в сверхзвуковом (M2 > 1) процессе (5).
Экспоненциально-тригонометрическая нелинейность источника
Двухкомпонентная система (2) имеет следующее точное решение [4]:
С = exP(-T1 ) cos т2 [2m2 Aj (Aj + 2s cos(mx')) - A + 2Al (sm / w)2 sin2 (mx')] +
+ 2m2 exp(-Tj / 2) cos(T2 / 2)[(AjS / w2) cos(mx') - AJ; (8)
k(v2) = exp(-Tj)sin т2[-2m2 Aj(Aj + 2s cos(mx')) + A - 2A2(sm / w)2 sin2(mx')] +
+ 2m2 exp(-Tj /2)sin(T2 /2)[-(A2s/w2)cos(mx') + AJ; (9)
Tj(x', t) = ln[(Aj + (a + P)cos(mt) + s cos(mx'))2 + (a - P)2 sin2(mt)];
т2( x', t) = 2arctg D;
D = (a - P) sin(mt) /[ Aj + (a + P )cos(mt) + s cos(mx')].
Здесь три константы w, Aj, A2 связаны двумя равенствами:
w2 = wj2 - Aj = w2 - A2 > 0,
где Aj, A2 - параметры источников, w - параметр растяжения координаты, x'= x /w. Постоянную Aj ^ 0 задаем так, чтобы существовала конечная функция т1(х',t). Полученные источники k1^a)(Tj,т2,x); a = 1,2 являются нелинейными по температурам т1, т2 и обладают периодической неоднородностью по координате x. Величина (a-P) является параметром двухкомпонентности системы: если a = P, то т2 = 0; т1 =т1(x',t). Ясно, что s - параметр пространственной неоднородности источников. При s = 0 обе температуры зависят только от времени и представляют собой точное решение системы уравнений:
d2та(t)/dt2 = k(v])(т1, т2); а = 1,2.
Этот результат переформулируем для динамической системы (3):
Q1(t1, т2) = (2m2Aj2 - A)exp(-Tj)cos т2 - 2m2Aj exp(-Tj / 2)cos(t2 /2); (10)
Q2(t1, т2) = (A - 2m2Aj2)exp(-Tj)sin т2 + 2m2Aj exp(-Tj / 2)sin(T2 / 2); (11)
т1(£) = ln[(Aj + (a + P )cos(m£))2 + (a-P )2sin2(m£)]; (12)
т2(£) = 2arctg D; A = 8m2aP; (13)
D = (a-P ) sin(m£) /[ Aj + (a + P ) cos(m£)];
kua) (T, T2 ) = A*wa (Ma - 1)Qa (Tj, T2^ ° = ^ 2
т
5- '■w' ч
1 1 U? © 0 - 5 10 5 £
а) m = 0,1
б) m = 0,95
в) m = 0,05
г) m = 0,5
д) m = 0,9
Q
\ -1-2-
Рис. 1. Источник типа двойного синус-Гордона: а, б - разрывные колебания температуры; в, г, д - непрерывные колебания температуры
В дозвуковом и сверхзвуковом процессах, когда контрастности нет, полученные источники энергии обладают такими свойствами:
1) одинаковые знаки скоростей изменения источников по температуре «своей» компоненты:
2) противоположные знаки скоростей изменения источников по температуре «чужой» компоненты:
Характер конкуренции источников определяется следующими свойствами. Функция Q2 (т1, т2) обращается в ноль при т2 = 0; именно на этой изотерме первый источник имеет экстремум по «чужой» температуре:
В структуре решения (12), (13) важная роль принадлежит отношению Д/(а + Р): если числовые значения этих констант такие, что в отдельных точках выполняется равенство А1 + (а + Р)соБ(да£) = 0, то температура т2(£) совершает пилообразные колебания. Таким образом, решение (6) (однокомпонентная система, рис. 1, а, б) и решение (12), (13) (двухкомпонентная система, рис. 2) дают примеры нетривиального воздействия объемного источника энергии на среду. А именно: непрерывные источники ^ (Т)
и ^С)Т Т2); а = 1,2 возбуждают разрывные колебания температуры.
sgn[д^1) / дТ ] = sgn[д^2) / дТ2];
^2 = 0; Q2 = 0; дQl/ дх2 = 0.
X
X-
а)
б)
в)
г)
Рис. 2. Двухкомпонентная динамическая система (10), (11): разрывные колебания температуры одной из компонент: а, б (А1 = 1; т = 0,2; а = 1; р = 0,3); в, г (А1 = -1; т = 0,2; а = 1; р = 2)
Взяв в (12), (13) а = р, получаем, что поведение однокомпонентной системы (5) определяется решением:
Q1 = Q = (2т2Д2 - А) ехр(-т) - 2т2 А1 ехр(-т / 2); т1 = т(£) = 1п[А1 + 2аcos(m£)]2; х2 = 0; А = 8т2а2; А12 > 4а2. (14)
Для этого решения характерны источники с отрицательным наклоном dQ / dт < 0 и справедливы замечания, сделанные при обсуждении рис. 1, д. Точные решения (6), (7), (14) обладают свойством «отрицательной теплоемкости». Этот термин означает, что среда реагирует аномальным образом на воздействие источника энергии: подвод/отвод тепла дает снижение/рост температуры. Следовательно, при «отрицательной теплоемкости» необходимо иметь: 1) если qu > 0, то дТ/дt < 0; 2) если qu < 0, то
дТ/дt > 0. Обзор экспериментальных и теоретических работ по этой проблеме и
примеры «отрицательной теплоемкости» в задачах конвекции стратифицированных двухкомпонентных жидкостей в поле силы тяжести даны в [5]. Нелинейные аспекты явления «отрицательной теплоемкости» изучены в [6], [7] с позиции теории локально-неравновесного теплопереноса. Для представленных здесь автомодельных решений имеем дТ/д^ = В^Т/^, и можем сопоставить знаки dт/^ и Q(т). Например, на рис. 1, а, б имеем интервалы температур, для которых Q(т) > 0, но при этом производная dт / ^ знакопеременная: наблюдается чередование областей с нормальным и аномальным температурными откликами среды. Аналогичное явление демонстрируют графики на рис. 1, в, г, д.
Если а = Р, но в ^ 0, то из (8), (9) получаем однокомпонентную систему с источником ku (т, х), который нелинеен по температуре и неоднороден по координате х;
в - параметр пространственной неоднородности. Для этой однокомпонентной системы проведем расчеты на линии х' = МЛ; I > 0, где М - тепловое число Маха. Дозвуковой (М < 1); звуковой (М = 1) и сверхзвуковой (М > 1) варианты показаны на рис. 3; здесь А1 = 2,57; т = 1; а = 1; в = 0,5. Во всех трех случаях наблюдается аномальный температурный отклик: существуют конечные интервалы времени, для которых kudт/dt < 0.
а) M = 0,05
к
O
-20 -40 і -60
Ж
б) M = 1
dT
t5 d t I,
°\ \, 'o \ 20 t
-5- I 1 \
-10J
в) М = 1,5
,Рмс. 3. Аномальный температурный отклик среды на воздействие нелинейного и неоднородного по координате источника (8) [ [ а = р; є Ф 0]: а - дозвуковой режим; б - звуковой режим; в - сверхзвуковой режим колебаний температуры на линии X = Mt
Периодическая неоднородность источника энергии
Уравнение (1) имеет точное решение [8]:
T = [A(а)]ч[е(г)-A(а)]; г^-^-J
A (а)
dp
(15)
A (а) JB (p)
Є = -2/z; B(p) = B = const - 0; а = x'+t; p = x'-t; x' = x/w;
k. = ATb (a+AT )2.
где A (а) - произвольная функция; точка над символом функции означает дифференцирование по ее аргументу. В данном случае выделение/поглощение энергии зависит явным образом от волновой координаты а : линия а = const есть характеристика (волна), распространяющаяся с постоянной скоростью влево, в сторону отрицательных значений x. В таком тепловом процессе может наблюдаться нетривиальное поведение изотерм T (а, p)= T = const. Обсудим этот вопрос. Линия изотермы и ее безразмерная скорость такие:
AЫ* =(-2/z)-A(а); [£] = гШ; BD-1; D(а)= A
1 2 (А A t); а = аг (t).
2
Следует отметить, что величина ^J, = ~ ~x = ~ = M представляет собой тепловое число Маха, равное отношению скорости изотермы N = dx, / dt к скорости w распространения тепловых возмущений. Допустим, что неоднородность источника энергии по волновой координате имеет периодический характер: A (а) = A0 + A1 sin ша; A0, Д, ш - const. Остановки изотермы (N (а = а г ) = 0) существуют, если выполнены два условия: а = аr; A - AT{ = 0; A + B = 0.
Отсюда находим:
шаr = arctg(-T/ш)+лл0; B = -A0 + ATt(t2 +ш2)-1/2,
где n0 - положительное четное целое число; B > (- А) > A1 > 0. Координата точки остановки определяется в каждый фиксированный момент времени из формулы X / w = а r -1.
Проанализируем зависимость M = M (t); t > 0 при A0 + B = 0. Ясно, что
d а, (dX л л 2
’ ' 1 +1 = M +1 = -
dt I dt I 1 - ББ
V У ?
Представляет интерес процесс, в котором |л|>|А1, т*е* А (а) - функция, колеблющаяся около постоянного значения А0, и знак ее такой же, как у Д. Значит подходящим выбором констант А1, ш можем обеспечить условие 0 <М1 < ddа_ < М2 <го.
(Л I
Таким образом, на изотерме Т = Т функция а = аг. (t) монотонно возрастающая, ее можно считать аналогом времени t > 0. В итоге имеем:
dX^ = A1sinша(А - AT, J- 2B(A - AT)
,dt J, (A1 sin ша-2B )(A - AT, )2 + 2B (A - A T)'
(16)
Знаменатель данной дроби колеблется около постоянного значения 2A^T2 и сохраняет знак этой величины. Числитель дроби - суперпозиция двух периодических функций, таких, что скорость (16) колеблется около нуля с периодически по а меняющейся амплитудой. Это явление есть следствие периодической неоднородности источника энергии по волновой координате.
Изучим на изотерме T = T корреляцию «источник энергии - тепловое число Ма-
(-2)T
ха». Функция к (T = T, а) = , , , —гЧ-г Аш cos ша + T (A0 + A, sin ша)? < 0 оп-
и ' Д(А0 + A1 sin ша^ 1 01
ределяет периодическую неоднородность стока энергии. Здесь принимаем А0 > А1 > 0; берем входные параметры задачи, для которых выполнено условие
1 > BD = -ДD, т. е. M +1 > 0. Функция kD (T, а) обращается в ноль в отдельных точках а = а0, если Д = А1. Если 1 + sin ша 0 = 0, то cos ша 0 = 0; ша0 =тсп0 + (-1)”0—;
п0 = 0, ± 1, ±2,..., и тогда ^(Т = Т{,а = а0) = 0. Результаты расчета этого варианта даны на рис. 4, а.
г) Т, = 0,5
Рис. 4. Корреляция «источник энергии - тепловое число Маха».
Входные параметры: а - А0= 1; А1 = 1; ш =0,2; б - А0= 1; А1 = 0,7; ш =0,1; в - А0=1; А1 = 0,125; ш = 0,5; г - А0= 1; А1 = 0,5; ш = 0,2
Видим, что число Маха знакопеременное, т. е. периодически появляются точки остановки изотермы (й'х' / dt), = 0; тепловой процесс «дозвуковой» [М| < 1] на «холодных» [Т < 1] и на «горячих» [Т > 1] изотермах. Гладкому периодическому изменению стока энергии ku (Т, а) соответствуют гладкие периодические колебания скорости изотермы. А именно: график М = М(а) не имеет резких изломов. Ситуация меняется (рис. 4, б, в), если отсутствуют точки вырождения [^(Т = Т,а = а0) = 0] стока энергии. Для таких режимов всюду ^ (Т = Т, а) < 0, а на «холодных» изотермах [Т ~ 0,2] скорость перемещения, т. е. число Маха, имеет импульсно-
периодический характер. Следовательно, гладкой периодической неоднородности стока энергии соответствуют импульсы, профиль которых имеет острые пики. Причина этого явления - нелинейная (кубическая) зависимость от температуры стока ^ (Т, а):
д Ь 2
1УГ= АБ (А + АТ)(А + 3 АТ).
Формирование острых пиков импульсов скорости присуще сверхзвуковому [М > 1] режиму движения изотермы в «холодной» области. По мере перехода к изотермам с более высокой температурой (рис. 4, г) исчезают резкие пики импульсов скорости; в этих температурных интервалах движение дозвуковое; обе функции kv (Т, а) и М(а) изменяются по квазигармоническим законам.
Еще один важный параметр процесса - частота ш колебаний по волновой координате а. Результаты расчетов говорят о том, что во всех рассмотренных вариантах рост ш приводит к увеличению частоты колебаний М(а). Амплитудные параметры неоднородности стока энергии А0, А1 оказывают только количественное влияние на свойства теплопереноса; качественный характер корреляции «сток энергии - число Маха» сохраняется при варьировании А0, А1.
Заключение
1. Выполнен анализ новых точных решений волнового уравнения теплопереноса с источниками, содержащими тригонометрическую, экспоненциально-тригонометрическую и кубическую нелинейности.
2. Представлены примеры (рис. 1, а, б; рис. 2) нетривиального взаимодействия объемного источника энергии со средой, когда непрерывный источник kv (Т, х) генерирует разрывные колебания температуры в однокомпонентных и двухкомпонентных системах.
3. В рамках модели волнового теплопереноса обнаружен аномальный температурный отклик среды на тепловое воздействие: подвод/отвод тепла дает снижение/рост температуры, т. е. quдT / дt < 0 (рис. 1, 3).
4. Подробно изучена пространственно-периодическая неоднородность температурных полей с кубической нелинейностью источников.
Литература
1. Жоу, Д. Расширенная необратимая термодинамика / Д. Жоу, Х. Касас-Баскес, Дж. Лебон. - М. : Ижевск : НИЦ «Регулярная и хаотическая динамика», 2006. -528 с.
2. Шабловский, О. Н. Релаксационный теплоперенос в нелинейных средах / О. Н. Шаб-ловский. - Гомель : ГГТУ им. П. О. Сухого, 2003. - 382 с.
3. Шабловский, О. Н. Неравновесные тепловые структуры в средах с источниками энергии / О. Н. Шабловский, Д. Г. Кроль. - Гомель : ГГТУ им. П. О. Сухого, 2013. - 208 с.
4. Шабловский, О. Н. Точные решения волновых уравнений с нелинейными источниками / О. Н. Шабловский // Фундаментальные физико-математические проблемы и моделирование технико-технологических систем. - Вып. 14. - М. : Янус-К, 2011. - С. 382-391.
5. Ингель, Л. Х. «Отрицательная теплоемкость» стратифицированных жидкостей / Л. Х. Ингель // УФН. - 2002. - Т. 172, № 6. - С. 691-699.
6. Шабловский, О. Н. Тепловая неустойчивость и «отрицательная теплоемкость» в локально-неравновесной среде с источником энергии / О. Н. Шабловский // Фазовые переходы, критические и нелинейные явления в конденсированных средах : сб. тр. междунар. конф. - Махачкала : Ин-т физики Дагестан. науч. центра РАН, 2005. - С. 295-298.
7. Шабловский О. Н. «Отрицательная теплоемкость» в простой тепловой волне / О. Н. Шабловский // Труды IV Российской национальной конференции по теплообмену : в 8 т. - М. : Изд-во МЭИ, 2006. - Т. 7 : Теплопроводность, теплоизоляция. - С. 359-361.
8. Шабловский, О. Н. Нелинейные волновые уравнения и конкуренция источников энергии в двухкомпонентных системах / О. Н. Шабловский // Фундаментальные физико-математические проблемы и моделирование технико-технологических систем. - Вып. 13. - М. : Янус-К, 2010. - С. 78-89.
Получено 17.02.2014 г.
